Skip to content

Commit

Permalink
Changed ohara 2011 state order
Browse files Browse the repository at this point in the history
  • Loading branch information
MichaelClerx committed Aug 15, 2024
1 parent aa9565b commit 62a06b8
Show file tree
Hide file tree
Showing 3 changed files with 9 additions and 9 deletions.
14 changes: 7 additions & 7 deletions c/bartolucci-2020.mmt
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@ version: 20240815
mmt_authors: Michael Clerx
desc: """
The 2020 "BPS" model of the human ventricular AP by Bartolucci et al. [1].

The model is based on the cipa adaptation [2] of the O'Hara model [3].
Changes to [2] are indicated throughout the code.

Expand All @@ -13,11 +13,11 @@ desc: """
by comparing the calculates derivatives to those calculated from both the
original matlab code and the author-provided CellML, both of which matched
to within machine precision.

After verification, the INaK formulation was corrected as described in [6].
Note that the CellML version uses an elevated external calcium of 2.7,
while the paper and the matlab code use the default value of 1.8.

References:

[1] Bartolucci C, Passini E, Hyttinen J, Paci M & Severi S (2020).
Expand Down Expand Up @@ -411,7 +411,7 @@ dot(jnca) = (inf - jnca) / 1 [ms]
dot(nca) = 1000 [1/ms] * anca - 150 [1/ms] * jnca * nca
anca = (1 - nca) / (1 + 0.05 [mM] / Ca_ss)^4
r_down = 0.1 [1/ms]
in [1/ms] # or ms ?
in [1/ms]
r_up = r_down * nca / (1 - nca)
in [1/ms]
# Omega and psi rate
Expand Down Expand Up @@ -494,7 +494,7 @@ ICaL = (icalp.ICaLp * camk.f + ICaLnp * (1 - camk.f))
in [A/F]
ICaNa = (icalp.ICaNap * camk.f + ICaNanp * (1 - camk.f))
in [A/F]

#
# Phosphorylated ICaL
#
Expand Down Expand Up @@ -1006,7 +1006,7 @@ g_irel_max = 0.02 [1/ms]
in [1/ms]
Jrelnp = frel * g_irel_max * SRCass * o * c * (Ca_sr - Ca_ss)
in [mM/ms]
Jrelp = 1.25 * frel * g_irel_max * SRCass * o * cp * (Ca_sr - Ca_ss)
Jrelp = 1.25 * frel * g_irel_max * SRCass * o * cp * (Ca_sr - Ca_ss)
in [mM/ms]
Jrel = (1 - camk.f) * Jrelnp + camk.f * Jrelp
desc: SR Calcium release flux via Ryanodine receptor
Expand All @@ -1022,7 +1022,7 @@ f = if(cell.mode == 1, 1.3, 1)
j = 0.004375 [mM/ms]
in [mM/ms]
k = 0.00092 [mM]
in [mM]
in [mM]
Jupnp = j * f * Ca_i / (Ca_i + k)
in [mM/ms]
Jupp = 2.75 * j * f * Ca_i / (Ca_i + k - 0.00017 [mM])
Expand Down
2 changes: 1 addition & 1 deletion c/ohara-2011.mmt
Original file line number Diff line number Diff line change
Expand Up @@ -83,9 +83,9 @@ ical.fs = 1
ical.fcaf = 1
ical.fcas = 1
ical.jca = 1
ical.nca = 0.2
ical.ffp = 1
ical.fcafp = 1
ical.nca = 0.2
ikr.xf = 0
ikr.xs = 0
iks.x1 = 0
Expand Down
2 changes: 1 addition & 1 deletion c/trovato-2020.mmt
Original file line number Diff line number Diff line change
Expand Up @@ -95,7 +95,7 @@ ryr.Jrel2 = 1.25045800437316998e-69
IP3.u = 4.66236137183557997e-1

#
# Simulation variables
# Simulator variables
#
[engine]
time = 0 [ms]
Expand Down

0 comments on commit 62a06b8

Please sign in to comment.