Updated IK1 comments in O'Hara-derived models.

c/bartolucci-2020.mmt

@@ -649,14 +649,12 @@ use extra.K_o
 V = membrane.V + rev.EKshift
     in [mV]
 dot(x) = (inf - x) / tau
-    desc: Activation of IK1
+    desc: Time-dependent inactivation at hyperpolarised (low) potentials
     inf = 1 / (1 + exp(-(V + 2.5538 [mV/mM] * K_o + 144.59 [mV]) / (1.5692 [mV/mM] * K_o + 3.8115 [mV])))
-        desc: Steady-state value for activation of IK1
     tau = 122.2 [ms] / (exp((V + 127.2 [mV]) / -20.36 [mV]) + exp((V + 236.8 [mV]) / 69.33 [mV]))
-        desc: Time constant for activation of IK1
         in [ms]
 r = 1 / (1 + exp((V + 105.8 [mV] - 2.6 [mV/mM] * K_o) / (1.09 * 9.493 [mV])))
-    desc: Inactivation of IK1
+    desc: Instantaneous rectification
 fK1 = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.2, 1.3) * 0.71
 gK1 = 0.1908 [mS/uF]
     in [mS/uF]

c/ohara-cipa-v1-2017.mmt

@@ -628,14 +628,14 @@ r2 = Ku * k53 / k63
 r3 = Kt / (1 + exp(-(V - Vhalf) / 6.789 [mV]))
     in [1/ms]
 # States
-dot(C2) =  (k1 * C1 - k2 * C2) - (k31 * C2 - k41 * O) - (k52 * C2 - k62 * IC2)
-dot(C1) = -(k1 * C1 - k2 * C2) - (k51 * C1 - k61 * IC1)
-dot(O) =  (k31 * C2 - k41 * O) - (k53 * O - k63 * IO) - (r1 * O - Ku * Obound)
-dot(IC2) =  (k11 * IC1 - k21 * IC2) - (k3 * IC2 - k4 * IO) + (k52 * C2 - k62 * IC2)
 dot(IC1) = -(k11 * IC1 - k21 * IC2) + (k51 * C1 - k61 * IC1)
+dot(IC2) =  (k11 * IC1 - k21 * IC2) - (k3 * IC2 - k4 * IO) + (k52 * C2 - k62 * IC2)
 dot(IO) = (k3 * IC2 - k4 * IO) + (k53 * O - k63 * IO) - (r1 * IO - r2 * IObound)
-dot(Obound) = (r1 * O - Ku * Obound) + (r3 * Cbound - Kt * Obound)
+dot(C1) = -(k1 * C1 - k2 * C2) - (k51 * C1 - k61 * IC1)
+dot(C2) =  (k1 * C1 - k2 * C2) - (k31 * C2 - k41 * O) - (k52 * C2 - k62 * IC2)
+dot(O) =  (k31 * C2 - k41 * O) - (k53 * O - k63 * IO) - (r1 * O - Ku * Obound)
 dot(IObound) = (r1 * IO - r2 * IObound) + (r3 * Cbound - Kt * IObound)
+dot(Obound) = (r1 * O - Ku * Obound) + (r3 * Cbound - Kt * Obound)
 dot(Cbound) = -(r3 * Cbound - Kt * Obound) - (r3 * Cbound - Kt * IObound)
 # Current
 fKr = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.3, 0.8)
@@ -682,14 +682,12 @@ IKs = fKs * gKs * KsCa * x1 * x2 * (V - rev.EKs)
 use membrane.V
 use extra.K_o
 dot(x) = (inf - x) / tau
-    desc: Activation of IK1
+    desc: Time-dependent inactivation at hyperpolarised (low) potentials
     inf = 1 / (1 + exp(-(V + 2.5538 [mV/mM] * K_o + 144.59 [mV]) / (1.5692 [mV/mM] * K_o + 3.8115 [mV])))
-        desc: Steady-state value for activation of IK1
     tau = 122.2 [ms] / (exp((V + 127.2 [mV]) / -20.36 [mV]) + exp((V + 236.8 [mV]) / 69.33 [mV]))
-        desc: Time constant for activation of IK1
         in [ms]
 r = 1 / (1 + exp((V + 105.8 [mV] - 2.6 [mV/mM] * K_o) / 9.493 [mV]))
-    desc: Inactivation of IK1
+    desc: Instantaneous rectification
 fK1 = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.2, 1.3) * 1.698
 gK1 = 0.1908 [mS/uF]
     in [mS/uF]