From 91fb907970b9df82fa02279a2e9deb1ddc630479 Mon Sep 17 00:00:00 2001 From: Michael Clerx Date: Wed, 4 Sep 2024 19:07:04 +0100 Subject: [PATCH] Updated annotations --- c/aguilar-2017.mmt | 1 + c/bai-2018.mmt | 6 ++++-- c/bartolucci-2020.mmt | 1 + c/carro-2011.mmt | 5 +++-- c/courtemanche-1998.mmt | 1 + c/ellinwood-2017.mmt | 5 +++-- c/fabbri-2017.mmt | 4 +++- c/fink-2008.mmt | 31 ++++++++++++++++--------------- c/grandi-2010.mmt | 5 +++-- c/grandi-2011.mmt | 5 +++-- c/iyer-2004.mmt | 1 + c/koivumaki-2011.mmt | 23 +---------------------- c/loewe-2019.mmt | 14 ++++++++------ c/ni-2017.mmt | 2 ++ c/nygren-1998.mmt | 1 + c/ohara-2011.mmt | 1 + c/ohara-cipa-v1-2017.mmt | 1 + c/priebe-1998.mmt | 1 + c/sampson-2010.mmt | 2 ++ c/stewart-2009.mmt | 4 +++- c/tentusscher-2004.mmt | 5 +++-- c/tentusscher-2006.mmt | 6 ++++-- c/tomek-2020.mmt | 1 + c/trovato-2020.mmt | 2 ++ c/voigt-2013.mmt | 5 +++-- 25 files changed, 72 insertions(+), 61 deletions(-) diff --git a/c/aguilar-2017.mmt b/c/aguilar-2017.mmt index c45ce71..630b050 100644 --- a/c/aguilar-2017.mmt +++ b/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 diff --git a/c/bai-2018.mmt b/c/bai-2018.mmt index 213e7b7..b14464f 100644 --- a/c/bai-2018.mmt +++ b/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 diff --git a/c/bartolucci-2020.mmt b/c/bartolucci-2020.mmt index 414e7d7..ae0f098 100644 --- a/c/bartolucci-2020.mmt +++ b/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] diff --git a/c/carro-2011.mmt b/c/carro-2011.mmt index cbb7b21..65a8144 100644 --- a/c/carro-2011.mmt +++ b/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])) diff --git a/c/courtemanche-1998.mmt b/c/courtemanche-1998.mmt index 1515aaa..5dea739 100644 --- a/c/courtemanche-1998.mmt +++ b/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 diff --git a/c/ellinwood-2017.mmt b/c/ellinwood-2017.mmt index 8b4d30d..3039082 100644 --- a/c/ellinwood-2017.mmt +++ b/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])) diff --git a/c/fabbri-2017.mmt b/c/fabbri-2017.mmt index cd9f06a..6419bef 100644 --- a/c/fabbri-2017.mmt +++ b/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])) diff --git a/c/fink-2008.mmt b/c/fink-2008.mmt index 0872c35..444bcac 100644 --- a/c/fink-2008.mmt +++ b/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) diff --git a/c/grandi-2010.mmt b/c/grandi-2010.mmt index 52c5ae9..76d2bb0 100644 --- a/c/grandi-2010.mmt +++ b/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])) diff --git a/c/grandi-2011.mmt b/c/grandi-2011.mmt index 3202c65..28830df 100644 --- a/c/grandi-2011.mmt +++ b/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])) diff --git a/c/iyer-2004.mmt b/c/iyer-2004.mmt index 809e4b9..31773c2 100644 --- a/c/iyer-2004.mmt +++ b/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 diff --git a/c/koivumaki-2011.mmt b/c/koivumaki-2011.mmt index 4c13feb..e81e49d 100644 --- a/c/koivumaki-2011.mmt +++ b/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 # diff --git a/c/loewe-2019.mmt b/c/loewe-2019.mmt index 2a9cec3..f2e619f 100644 --- a/c/loewe-2019.mmt +++ b/c/loewe-2019.mmt @@ -1,6 +1,6 @@ [[model]] name: loewe-2019 -version: 20240903 +version: 20240904 mmt_authors: Michael Clerx display_name: Loewe et al., 2019 desc: """ @@ -9,8 +9,8 @@ desc: """ This Myokit implementation is based on the CellML code by Axel Loewe [4]. Two errors in the CellML file were corrected before import. The Myokit - implementation has extensive reformatting and unit corrections, and - includes a small offset in IKACh to avoid divide-by-zero errors, as + implementation has extensive reformatting and unit corrections, and + includes a small offset in IKACh to avoid divide-by-zero errors, as suggested by Alan Fabbri. It was tested against the original CellML by comparing the calculated derivatives, which matched to within machine precision. @@ -102,6 +102,7 @@ i_tot = (+ if.If + isk.ISK ) in [nA] + label cellular_current # # Acetylcholine @@ -314,10 +315,11 @@ dot(r) = (inf - r) / tau # [ikr] use membrane.V -gKr = 4.989099e-3 [uS] * sqrt(extra.K_o / 5.4 [mM]) - in [uS] -IKr = gKr * (0.9 * paf + 0.1 * pas) * piy * (V - rev.E_K) +IKr = gKr * sqrt(extra.K_o / 5.4 [mM]) * (0.9 * paf + 0.1 * pas) * piy * (V - rev.E_K) in [nA] +gKr = 4.989099e-3 [uS] + in [uS] + label g_Kr a_inf = 1 / (1 + exp(-(V + 10.0144 [mV]) / 7.6607 [mV])) dot(paf) = (a_inf - paf) / tau tau = 1 [s] / (30 * exp(V / 10 [mV]) + exp(-V / 12 [mV])) diff --git a/c/ni-2017.mmt b/c/ni-2017.mmt index 13d0014..b3837de 100644 --- a/c/ni-2017.mmt +++ b/c/ni-2017.mmt @@ -138,6 +138,7 @@ i_ion = ( + icap.ICap ) in [pA] + label cellular_current dot(V) = -(i_ion + i_stim + i_diff) / cell.Cm in [mV] label membrane_potential @@ -315,6 +316,7 @@ dot(xr) = (inf - xr) / tau inf = 1 / (1 + exp((V + 14.1 [mV]) / -6.5 [mV])) gKr = 0.8 * 0.029411765 [nS/pF] in [nS/pF] + label g_Kr IKr = Cm * fKr * gKr * xr * (V - nernst.EK) / (1 + exp((V + 15 [mV]) / 22.4 [mV])) in [pA] diff --git a/c/nygren-1998.mmt b/c/nygren-1998.mmt index bb22dfb..6cd945c 100644 --- a/c/nygren-1998.mmt +++ b/c/nygren-1998.mmt @@ -301,6 +301,7 @@ 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 - rev.EK) in [pA] diff --git a/c/ohara-2011.mmt b/c/ohara-2011.mmt index e855231..f9eafb5 100644 --- a/c/ohara-2011.mmt +++ b/c/ohara-2011.mmt @@ -515,6 +515,7 @@ fKr = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.3, 0.8) gKr = 0.046 [mS/uF] in [mS/uF] desc: Maximum conductance of IKr + label g_Kr IKr = fKr * gKr * sqrt(extra.K_o / 5.4 [mM]) * x * r * (V - rev.EK) desc: Rapid delayed potassium current in [A/F] diff --git a/c/ohara-cipa-v1-2017.mmt b/c/ohara-cipa-v1-2017.mmt index 4d23f91..391c4bc 100644 --- a/c/ohara-cipa-v1-2017.mmt +++ b/c/ohara-cipa-v1-2017.mmt @@ -611,6 +611,7 @@ dot(Cbound) = -(r3 * Cbound - Kt * Obound) - (r3 * Cbound - Kt * IObound) fKr = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.3, 0.8) gKr = 4.65854545454545618e-2 [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] diff --git a/c/priebe-1998.mmt b/c/priebe-1998.mmt index a629c0c..da6efe9 100644 --- a/c/priebe-1998.mmt +++ b/c/priebe-1998.mmt @@ -241,6 +241,7 @@ i_Kr = g_Kr_max * xr * rik * (V - rev.E_K) in [A/F] g_Kr_max = 0.015 [mS/uF] in [mS/uF] + label g_Kr dot(xr) = alpha * (1 - xr) - beta * xr alpha = 0.005 [1/ms] * exp(0.0005266 [1/mV] * (V + 4.067 [mV])) / (1 + exp(-0.1262 [1/mV] * (V + 4.067 [mV]))) in [1/ms] diff --git a/c/sampson-2010.mmt b/c/sampson-2010.mmt index 6168c5c..c7eb231 100644 --- a/c/sampson-2010.mmt +++ b/c/sampson-2010.mmt @@ -142,6 +142,7 @@ dot(V) = -(i_ion + i_diff + i_stim) label membrane_potential i_ion = a1 + a2 + a3 in [A/F] + label cellular_current a1 = nav15.INa + nav11.INa1 + ical.ICa + icat.I + ical.ICaK + ikr.IKr + iks.IKs + ihcn.IHCN in [A/F] a2 = ik1.IK1 + inaca.INaCa + inak.INaK + ito.Ito1 + isus.Isus @@ -453,6 +454,7 @@ dot(I) = a1 - a2 desc: HERG channel state I1(Inactivated) GKr = 0.0383724 [mS/uF] : peak IKr conductance in [mS/uF] + label g_Kr fKo = sqrt(extra.Ko / 4 [mM]) IKr = GKr * fKo * O * (V - nernst.EK) desc: rapid activating delayed rectifier K+ current diff --git a/c/stewart-2009.mmt b/c/stewart-2009.mmt index 2285014..bf19dd8 100644 --- a/c/stewart-2009.mmt +++ b/c/stewart-2009.mmt @@ -16,7 +16,7 @@ desc: """ with links to the relevant document section. In addition, incorrect units in F, R, Cm, gCaL, Vc, Vss, and Vsr were corrected. - + The stimulus was set to 0.5ms and approximately twice the threshold for (immediate) depolarisation. @@ -97,6 +97,7 @@ dot(V) = -(i_ion + i_stim + i_diff) label membrane_potential i_ion = sodium.INa_tot + potassium.IK_tot + calcium.ICa_tot in [A/F] + label cellular_current i_diff = 0 [A/F] in [A/F] bind diffusion_current @@ -297,6 +298,7 @@ dot(Xr2) = (xr2_inf - Xr2) / tau_xr2 xr2_inf = 1 / (1 + exp((V + 88 [mV]) / 24 [mV])) g_Kr = 0.0918 [mS/uF] in [mS/uF] + label g_Kr i_Kr = g_Kr * sqrt(extra.K_o / 5.4 [mM]) * Xr1 * Xr2 * (V - nernst.E_K) in [A/F] diff --git a/c/tentusscher-2004.mmt b/c/tentusscher-2004.mmt index b45e72a..300a99d 100644 --- a/c/tentusscher-2004.mmt +++ b/c/tentusscher-2004.mmt @@ -186,10 +186,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.096 [mS/uF] * sqrt(extra.Ko / 5.4 [mM]) +gKr = 0.096 [mS/uF] 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 diff --git a/c/tentusscher-2006.mmt b/c/tentusscher-2006.mmt index c0955f4..a12f628 100644 --- a/c/tentusscher-2006.mmt +++ b/c/tentusscher-2006.mmt @@ -90,6 +90,7 @@ i_ion = (+ ina.INa + inab.INab ) in [A/F] + label cellular_current i_diff = 0 [A/F] in [A/F] bind diffusion_current @@ -209,10 +210,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.153 [mS/uF] * sqrt(extra.Ko / 5.4 [mM]) +gKr = 0.153 [mS/uF] 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 diff --git a/c/tomek-2020.mmt b/c/tomek-2020.mmt index 11b36e9..931b14f 100644 --- a/c/tomek-2020.mmt +++ b/c/tomek-2020.mmt @@ -550,6 +550,7 @@ dot(I) = alpha_C2ToI * C1 + alpha_i * O - (beta_ItoC2 + beta_i) * I fKr = piecewise(cell.mode == 0, 1, cell.mode == 1, 1.3, 0.8) gKr = 0.0321 [mS/uF] in [mS/uF] + label g_Kr IKr = fKr * gKr * sqrt(extra.K_o / 5 [mM]) * O * (V - rev.EK) in [A/F] diff --git a/c/trovato-2020.mmt b/c/trovato-2020.mmt index 821b641..0915cec 100644 --- a/c/trovato-2020.mmt +++ b/c/trovato-2020.mmt @@ -130,6 +130,7 @@ i_ion = (ina.INa + inal.INaL + inaca.INaCa + inacass.INaCa_ss + inak.INaK + inab.INab + ipca.IpCa + icab.ICab) in [A/F] + label cellular_current # # Stimulus current @@ -423,6 +424,7 @@ r = 1 / (1 + exp((V + 55 [mV]) / (0.32 * 75 [mV]))) * 1 / (1 + exp((V - 10 [mV]) desc: Instantaneous IKr inactivation gKr = 0.0342 [mS/uF] in [mS/uF] + label g_Kr IKr = gKr * sqrt(extra.K_o / 5.4 [mM]) * x * r * (V - rev.E_K) in [A/F] diff --git a/c/voigt-2013.mmt b/c/voigt-2013.mmt index 4682f60..1d3d348 100644 --- a/c/voigt-2013.mmt +++ b/c/voigt-2013.mmt @@ -350,9 +350,10 @@ dot(x) = (inf - x) / tau + 230 [ms] / (1 + exp((V + 40 [mV]) / 20 [mV]))) in [ms] r = 1 / (1 + exp((V + 74 [mV]) / 24 [mV])) -gKr = 0.035 [mS/uF] * sqrt(ion.K_o / 5.4 [mM]) +gKr = 0.035 [mS/uF] in [mS/uF] -IKr = gKr * x * r * (V - nernst.EK) + label g_Kr +IKr = gKr * sqrt(ion.K_o / 5.4 [mM]) * x * r * (V - nernst.EK) in [A/F] label I_Kr