Updated annotations

c/aguilar-2017.mmt

@@ -392,6 +392,7 @@ dot(xr) = (inf - xr) / tau
 gKr = 0.029411765 [nS/pF] # Extra digits 11765 from CellML file
     in [nS/pF]
     desc: Maximal IKr conductance
+    label g_Kr
 IKr = gKr * xr * (V - nernst.EK) / (1 + exp((V + 15 [mV]) / 22.4 [mV]))
     in [A/F]
     label I_Kr

c/bai-2018.mmt

@@ -91,6 +91,7 @@ i_ion = (+ ina.INa
          + inab.INab
         )
     in [A/F]
+    label cellular_current
 
 #
 # Stimulus current
@@ -196,10 +197,11 @@ inf = alpha / (alpha + beta)
 #
 [ikr]
 use membrane.V, nernst.EK
-IKr = gKr * xr1 * xr2 * (V - EK)
+IKr = gKr * sqrt(extra.Ko / 5.4 [mM]) * xr1 * xr2 * (V - EK)
     in [A/F]
-gKr = 0.1989 [mS/uF] * sqrt(extra.Ko / 5.4 [mM])    # Rescaled
+gKr = 0.1989 [mS/uF]    # Rescaled
     in [mS/uF]
+    label g_Kr
 dot(xr1) = (inf - xr1) / tau
     inf = 1 / (1 + exp((-26 [mV] - V) / 7 [mV]))
     tau = 1 [ms] * alpha * beta

c/bartolucci-2020.mmt

@@ -612,6 +612,7 @@ dot(Cbound) = -(r3 * Cbound - Kt * Obound) - (r3 * Cbound - Kt * IObound)
 fKr = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.1, 0.8) * 1.2
 gKr = 0.046 [mS/uF]
     in [mS/uF]
+    label g_Kr
 IKr = fKr * gKr * sqrt(extra.K_o / 5.4 [mM]) * O * (V - rev.EK)
     in [A/F]
 

c/carro-2011.mmt

@@ -306,10 +306,11 @@ KmNaip = 11 [mM]
 #
 [ikr]
 use membrane.V
-IKr = gKr * x * r * (V - rev.EK)
+IKr = gKr * sqrt(ion.K_o / 5.4 [mM]) * x * r * (V - rev.EK)
     in [A/F]
-gKr = 0.035 [mS/uF] * sqrt(ion.K_o / 5.4 [mM])
+gKr = 0.035 [mS/uF]
     in [mS/uF]
+    label g_Kr
 dot(x) = (inf - x) / tau
     inf = 1 / (1 + exp(-(V + 10 [mV]) / 5 [mV]))
     tau = (+ 550 [ms] / (1 + exp((-22 [mV] - V) / 9 [mV])) * 6 / (1 + exp((V + 11 [mV]) / 9 [mV]))

c/courtemanche-1998.mmt

@@ -401,6 +401,7 @@ dot(xr) = (inf - xr) / tau
 gKr = 0.029411765 [nS/pF] # Extra digits 11765 from CellML file
     in [nS/pF]
     desc: Maximal IKr conductance
+    label g_Kr
 IKr = gKr * xr * (V - nernst.EK) / (1 + exp((V + 15 [mV]) / 22.4 [mV]))
     in [A/F]
     label I_Kr

c/ellinwood-2017.mmt

@@ -380,11 +380,12 @@ KmNaip = 11 [mM] * (1 - 0.25 * mode.ISO)
 #
 [ikr]
 use membrane.V
-IKr = gKr * x * r * (V - rev.EK)
+IKr = gKr * sqrt(ion.K_o / 5.4 [mM]) * x * r * (V - rev.EK)
     in [A/F]
     label I_Kr
-gKr = 0.035 [mS/uF] * sqrt(ion.K_o / 5.4 [mM])
+gKr = 0.035 [mS/uF]
     in [mS/uF]
+    label g_Kr
 dot(x) = (inf - x) / tau
     inf = 1 / (1 + exp(-(V + 10 [mV]) / 5 [mV]))
     tau = (+ 550 [ms] / (1 + exp((-22 [mV] - V) / 9 [mV])) * 6 / (1 + exp((V + 11 [mV]) / 9 [mV]))

c/fabbri-2017.mmt

@@ -1,6 +1,6 @@
 [[model]]
 name: fabbri-2017
-version: 20240903
+version: 20240904
 mmt_authors: Michael Clerx, Aditi Agrawal
 display_name: Fabbri et al., 2017
 desc: """
@@ -97,6 +97,7 @@ i_tot = (
         + ikur.IKur
     )
     in [nA]
+    label cellular_current
 
 #
 # Acetylcholine
@@ -302,6 +303,7 @@ dot(r) = (inf - r) / tau
 use membrane.V
 gKr = 0.00424 [uS]
     in [uS]
+    label g_Kr
 IKr = gKr * (0.9 * paf + 0.1 * pas) * piy * (V - rev.E_K)
     in [nA]
 a_inf = 1 / (1 + exp(-(V + 10.0144 [mV]) / 7.6607 [mV]))

c/fink-2008.mmt

@@ -1,6 +1,6 @@
 [[model]]
 name: fink-2008
-version: 20240903
+version: 20240904
 mmt_authors: Michael Clerx, Chon Lok Lei
 display_name: Fink et al., 2008
 desc: """
@@ -9,8 +9,8 @@ desc: """
 
     (Note that there is a 2006 model accompanying a publication with a similar
     name.)
-    
-    In this Myokit implementation, the stimulus was set to 0.5ms and 
+
+    In this Myokit implementation, the stimulus was set to 0.5ms and
     approximately twice the threshold for depolarisation.
 
     Adapted from a CellML implementation by Catherine Lloyd [3].
@@ -29,7 +29,7 @@ desc: """
         https://doi.org/10.1152/ajpheart.00109.2006
 
     [3] https://models.physiomeproject.org/exposure/eeb81adc372c2f172399ec7160b0331e
-    
+
     """
 # Initial values
 membrane.V    = -86.45
@@ -77,17 +77,17 @@ pace = 0
 dot(V) = -(i_ion + stimulus.i_stim)
     in [mV]
     label membrane_potential
-i_ion = (+ ik1.i_K1 
-         + ito.i_to 
-         + ikr.i_Kr 
-         + iks.i_Ks 
-         + ical.i_CaL 
-         + inak.i_NaK 
-         + ina.i_Na 
-         + inab.i_b_Na 
-         + inaca.i_NaCa 
-         + icab.i_b_Ca 
-         + ipk.i_p_K 
+i_ion = (+ ik1.i_K1
+         + ito.i_to
+         + ikr.i_Kr
+         + iks.i_Ks
+         + ical.i_CaL
+         + inak.i_NaK
+         + ina.i_Na
+         + inab.i_b_Na
+         + inaca.i_NaCa
+         + icab.i_b_Ca
+         + ipk.i_p_K
          + ipca.i_p_Ca)
     in [A/F]
     label cellular_current
@@ -246,6 +246,7 @@ xK1_inf = phi * rec1 + (1 - phi) * rec2
 use membrane.V
 g_Kr_0 = 0.024 [mS/uF]
     in [mS/uF]
+    label g_Kr
 g_Kr = g_Kr_0 * (phys.T / 35 [K] - 55 / 7)
     in [mS/uF]
 i_Kr = g_Kr * sqrt(ion.K_o / 5.4 [mM]) * Or4 * (V - rev.E_K)

c/grandi-2010.mmt

@@ -302,10 +302,11 @@ KmNaip = 11 [mM]
 #
 [ikr]
 use membrane.V
-IKr = gKr * x * r * (V - rev.EK)
+IKr = gKr * sqrt(ion.K_o / 5.4 [mM]) * x * r * (V - rev.EK)
     in [A/F]
-gKr = 0.035 [mS/uF] * sqrt(ion.K_o / 5.4 [mM])
+gKr = 0.035 [mS/uF]
     in [mS/uF]
+    label g_Kr
 dot(x) = (inf - x) / tau
     inf = 1 / (1 + exp(-(V + 10 [mV]) / 5 [mV]))
     tau = (+ 550 [ms] / (1 + exp((-22 [mV] - V) / 9 [mV])) * 6 / (1 + exp((V + 11 [mV]) / 9 [mV]))

c/grandi-2011.mmt

@@ -337,11 +337,12 @@ KmNaip = 11 [mM] * (1 - 0.25 * mode.ISO)
 #
 [ikr]
 use membrane.V
-IKr = gKr * x * r * (V - rev.EK)
+IKr = gKr * sqrt(ion.K_o / 5.4 [mM]) * x * r * (V - rev.EK)
     in [A/F]
     label I_Kr
-gKr = 0.035 [mS/uF] * sqrt(ion.K_o / 5.4 [mM])
+gKr = 0.035 [mS/uF]
     in [mS/uF]
+    label g_Kr
 dot(x) = (inf - x) / tau
     inf = 1 / (1 + exp(-(V + 10 [mV]) / 5 [mV]))
     tau = (+ 550 [ms] / (1 + exp((-22 [mV] - V) / 9 [mV])) * 6 / (1 + exp((V + 11 [mV]) / 9 [mV]))

c/iyer-2004.mmt

@@ -420,6 +420,7 @@ IKr = GKr * fKo * O * (V - rev.EK)
     in [A/F]
 GKr = 0.0186 [mS/uF]
     in [mS/uF]
+    label g_Kr
 TC =  5.32000000100000037
 fKo = sqrt(extra.Ko / 4 [mM])
 # Rate parameters

c/koivumaki-2011.mmt

@@ -73,7 +73,6 @@ calcium.CaSR2 =  6.10104205027038882e-01
 calcium.CaSR3 =  5.90926659120631581e-01
 calcium.CaSR4 =  5.72078107967216365e-01
 
-
 #
 # External variables
 #
@@ -113,7 +112,6 @@ i_diff = 0 [pA]
     in [pA]
     bind diffusion_current
 
-
 #
 # Stimulus current
 #
@@ -128,7 +126,6 @@ amplitude = 2 * -2943 [pA]
         approximately the smallest value that triggers three successive APs.
         """
 
-
 #
 # Fast sodium current
 #
@@ -152,7 +149,6 @@ INa = PNa * m^3 * (0.9 * h1 + 0.1 * h2) * extra.Nao * V * F * FRT * (exp((V - ne
     in [pA]
     label I_Na
 
-
 #
 # L-type calcium current
 #
@@ -184,7 +180,6 @@ ICaL = gCaL * d * fca * f1 * f2 * (V - ECa_app)
     in [pA]
     label I_CaL
 
-
 #
 # Transient outward potassium current
 #
@@ -206,7 +201,6 @@ It = gt * r * s * (V - nernst.EK)
     in [pA]
     label I_to
 
-
 #
 # Ultra-rapid potassium current (IKur or Isus)
 #
@@ -243,7 +237,6 @@ gKs = 1 [nS]
 IKs = gKs * n * (V - nernst.EK)
     in [pA]
 
-
 #
 # Rapid delayed-rectifier potassium current
 #
@@ -256,11 +249,11 @@ dot(pa) = (inf - pa) / tau
 pi = 1 / (1 + exp((V + 55 [mV]) / 24 [mV]))
 gKr = 0.5 [nS]
     in [nS]
+    label g_Kr
 IKr = gKr * pa * pi * (V - nernst.EK)
     in [pA]
     label I_Kr
 
-
 #
 # Hyperpolarisation-activated potassium current
 #
@@ -280,7 +273,6 @@ IfK  = gIf * y * ((1 - 0.2677) * (V - nernst.EK))
 If = IfK + IfNa
     in [pA]
 
-
 #
 # Inward-rectifier potassium current
 # Note: The factor 0.90 is due to Koivumaki et al.
@@ -294,7 +286,6 @@ IK1 = gK1 * (extra.Ko * 1[1/mM])^0.4457 * (V - EK) / (1 + exp(1.5 * (V - EK + 3.
     in [pA]
     label I_K1
 
-
 #
 # Background sodium current
 #
@@ -304,7 +295,6 @@ gNab = 0.060599 [nS]
 INab = gNab * (membrane.V - nernst.ENa)
     in [pA]
 
-
 #
 # Background calcium current
 #
@@ -314,7 +304,6 @@ gCab = 0.0952 [nS]
 ICab = gCab * (membrane.V - nernst.ECa)
     in [pA]
 
-
 #
 # Sodium-potassium pump
 #
@@ -331,7 +320,6 @@ INaK = INaKmax * Ko / (Ko + kNaKK) * Nass^1.5 / (Nass^1.5 + kNaKNa^1.5) * (V + 1
     in [pA]
     label I_NaK
 
-
 #
 # Sodium-calcium exchanger
 #
@@ -349,7 +337,6 @@ INaCa = kNaCa * (exp(gam * V * FRT) * Nass^3 * Cao - exp((gam - 1) * V * FRT) *
     in [pA]
     label I_NaCa
 
-
 #
 # Calcium pump
 #
@@ -361,7 +348,6 @@ ICaPmax = 2 [pA]
 ICaP = ICaPmax * calcium.Cass / (kCaP + calcium.Cass)
     in [pA]
 
-
 #
 # External concentrations
 #
@@ -374,7 +360,6 @@ Ko = 5.4 [mM]
     in [mM]
     label K_o
 
-
 #
 # Reversal potentials
 #
@@ -387,7 +372,6 @@ EK  = RTF * log(extra.Ko / potassium.Ki)
 ECa = RTF * log(extra.Cao / calcium.Cass) / 2
     in [mV]
 
-
 #
 # Physical constants
 #
@@ -403,7 +387,6 @@ RTF = R * T / F
 FRT = 1 / RTF
     in [1/mV]
 
-
 #
 # Cell geometry
 #
@@ -494,7 +477,6 @@ xj_nj_Nai = 0.02 [um] / 2 + 2 * dx
     in [um]
     desc: diffusion distance from center of junct to center of njunct (between 2nd and 3rd njunct)
 
-
 #
 # Sodium concentrations
 #
@@ -529,7 +511,6 @@ betaNass = 1 / (1 + BNa * KdBNa / (Nass + KdBNa)^2)
 dot(Nass) = betaNass * (-JNa / cell.Vss - i_ss / (cell.Vss * phys.F))
     in [mM]
 
-
 #
 # Potassium concentration
 #
@@ -549,7 +530,6 @@ dot(Ki) = -i_tot / (cell.Vcytosol * phys.F)
     in [mM]
     label K_i
 
-
 #
 # Calcium concentrations
 #
@@ -718,7 +698,6 @@ dot(a3) = 0.5 * (-J_SERCASR3 + J_bulkSERCA3) / cell.Vnonjunct3
 dot(ass) = 0.5 * (-J_SERCASRss + J_bulkSERCAss) / cell.Vss
     in [mM]
 
-
 #
 # Ryanodine receptors
 #