XPP model

This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6.

# RatAPtest01.ode
# Rat action potential is simulated (bondarenko_model_2004)
# (cf. Bondarenko et al. Am J Physiol Heart Circ Physiol 2004;287: H1378-H1403)
# (cf. Wang et al. Toxicol Sci 2008;106:454-463)
# implemented by Dr. Sheng-Nan Wu

# Constant values
Number  Temp=298.0, R=8.314, Fara=96.5
Number  E_CaL=63.0, E_Cl=-40.0
Number Cm=1.0, Vmyo=2.584E-5, VNSR=2.098E-6, VJSR=1.2E-7

# Initial values
Initial V=-82.4202, Cai=0.115001
Initial Nai=14237.1, Ki=143720.0
Initial C_Na1=2.79132E-4, C_Na2=0.020752, O_Na=7.13483E-7
Initial IF_Na=1.53176E-4, I1_Na=6.73345E-7, I2_Na=1.55787E-9
Initial IC_Na2=0.0113879, IC_Na3=0.34278
Initial ato_f=0.00265563, ito_f=0.999977
Initial ato_s=4.17069E-4, ito_s=0.998543
Initial nKs=2.62753E-4, aur=4.17069E-4, iur=0.998543
Initial aKss=4.17069E-4, iKss=1.0
Initial C_K1=9.92513E-4, C_K2=6.41229E-4
Initial O_K=1.75298E-4, I_K=3.19129E-5
Initial Cass=0.115001, CaJSR=1299.5, CaNSR=1299.5
Initial P_RyR=0.0
Initial LTRPN_Ca=11.2684, HTRPN_Ca=125.29
Initial O=9.30308E-19, O1=1.49102E-5, O2=9.51726E-11
Initial P_C2=1.6774E-4
Initial C2=1.24216E-4, C3=5.78679E-9, C4=1.19816E-13
Initial I1=4.97023E-19, I2=3.45847E-14, I3=1.85106E-14

# Stimulus protocol
Par r0=1, period=300, pulse=10
Par tf=0, tp=5, tstart=10
ts=t-tstart
rstar = r0 + pulse*(heav(mod(ts,period)-tf)-heav(mod(ts,period)-(tf+tp)))

Par  Nao=140000.0
Par  LTRPN_tot=70.0
Par  kf=0.023761
Par  kb=0.036778
Par  Km_CSQN=800.0
Par  i_pCa_max=1.0
Par  CMDN_tot=50.0
Par  Kpc_max=0.23324
Par  Km_CMDN=0.238
Par  K_mCa=1380.0
Par  Km_Nai=21.0
Par  v1=4.5, v2=1.7E-5, v3=0.45
Par  k_plus_c=0.0090
Par  k_minus_ltrpn=0.196
Par  g_Cab=3.67E-4
Par  k_plus_b=0.00405
Par  k_plus_a=0.006075
Par  Km_up=0.5
Par  K_mNa=87500.0
Par  Km_Cl=10.0
Par  CSQN_tot=15000.0
Par  k_minus_htrpn	=3.2E-5
Par  Km_pCa=0.5
Par  k_plus_ltrpn=0.0327
Par  Acap=1.534E-4
Par  g_Na=13.0, g_CaL=1.51729, g_Ks=0.000575, g_Kr=0.0078, g_Kur=0.0016
Par  g_ClCa=10.0
Par  Vss=1.485E-9
Par  k_NaCa=992.8
Par  k_plus_htrpn=0.00237
Par  Ko=5400.0, Cao=1800.0
Par  tau_xfer=8.0
Par  k_sat	=0.1
Par  n=4.0
Par  g_Kto_f=0.4067, g_Kto_s=0.01, g_Nab=0.0026
Par  m=3.0
Par  Km_Ko=1.5
Par  g_Kss=0.05
Par  HTRPN_tot=140.0
Par  i_NaK_max=0.88
Par  i_CaL_max=7.0
Par  tau_tr=20.0
Par  k_minus_c=8.0E-4
Par  k_minus_b=0.965
Par  k_minus_a=0.07125
Par  Kpc_half=20.0
Par  eta=0.35
Par  Kpcb=5.0E-4

ass = (1.0 / (1.0 + exp( - (0.12987012987012986 * (22.5 + V)))))
E_K = (ln((Ko / Ki)) * R * Temp / Fara)
i_Ks = (g_Ks * (nKs^2.0) * (V - E_K))
E_Na = (ln((((0.9 * Nao) + (0.1 * Ko)) / ((0.9 * Nai) + (0.1 * Ki)))) * R * Temp / Fara)
i_Kr = (g_Kr * O_K * (V - (ln((((0.98 * Ko) + (0.02 * Nao)) / ((0.98 * Ki) + (0.02 * Nai)))) * R * Temp / Fara)))
i_Kur = (g_Kur * aur * iur * (V - E_K))
sigma = (0.14285714285714285 * (-1.0 + exp((1.4858841010401188E-5 * Nao))))
f_NaK = (1.0 / (1.0 + (0.1245 * exp( - (0.1 * V * Fara / (R * Temp)))) + (0.0365 * sigma * exp( - (Fara * V / (R * Temp))))))
J_leak = (v2 * (CaNSR - Cai))
i_NaCa =	 (k_NaCa * ((exp((eta * V * Fara / (R * Temp))) * (Nai^3.0) * Cao) - (exp(((-1.0 + eta) * V * Fara / (R * Temp))) * (Nao^3.0) * Cai)) / ((K_mNa^3.0) + (Nao^3.0)) / (K_mCa + Cao) / (1.0 + (k_sat * exp(((-1.0 + eta) * V * Fara / (R * Temp))))))
O_ClCa = (0.2 / (1.0 + exp( - (0.12820512820512822 * (-46.7 + V)))))
i_ClCa = (g_ClCa * O_ClCa * (V - E_Cl) * Cai / (Cai + Km_Cl))
i_Nab = (g_Nab * (V - E_Na))
tau_Kss = (13.17 + (39.3 * exp( - (0.0862 * V))))
beta_i1 = (9.5E-4 * exp((0.14285714285714285 * (33.5 + V))) / (1.0 + (0.051335 * exp((0.14285714285714285 * (33.5 + V))))))
beta = (0.05 * exp( - (0.07692307692307693 * (12.0 + V))))
i_Kss = (g_Kss * aKss * iKss * (V - E_K))
i_K1 = (0.2938 * Ko * (V - E_K) / (210.0 + Ko) / (1.0 + exp((0.0896 * (V - E_K)))))
i_NaK = (i_NaK_max * f_NaK * Ko / (1.0 + ((Km_Nai / Nai)^0.66667)) / (Ko + Km_Ko))
J_rel = (v1 * (O1 + O2) * (CaJSR - Cass) * P_RyR)
alpha_n = (4.81333E-6 * (26.5 + V) * (1.0 - exp( - (0.128 * (26.5 + V)))))
alpha_i=	 (0.090821 * exp((0.023391 * (5.0 + V))))
alpha_a =	 (0.18064 * exp((0.03577 * (30.0 + V))))
i_pCa = (i_pCa_max * (Cai^2.0) / ((Km_pCa^2.0) + (Cai^2.0)))
alpha_i1 = (1.52E-4 * exp( - (0.14285714285714285 * (13.5 + V))) / (1.0 + (0.0067083 * exp( - (0.14285714285714285 * (33.5 + V))))))
# i_stim =  - (10.0 * (t > 10.0) * (t < 15.0))
beta_Na13 = (0.22 * exp( - (0.04926108374384236 * (-7.5 + V))))
beta_Na12 =  (0.2 * exp( - (0.04926108374384236 * (-2.5 + V))))
beta_Na11 = (0.1917 * exp( - (0.04926108374384236 * (2.5 + V))))
alpha = (0.4 * exp((0.1 * (12.0 + V))) * (1.0 - (0.75 * exp( - (0.0025 * ((20.0 + V)^2.0)))) + (0.7 * exp( - (0.1 * ((40.0 + V)^2.0))))) / (1.0 + (0.12 * exp((0.1 * (12.0 + V))))))
Bi = ((1.0 + (CMDN_tot * Km_CMDN / ((Km_CMDN + Cai)^2.0)))^(-1.0))
gamma =	 (Kpc_max * Cass / (Kpc_half + Cass))
P_C1= (1.0 - (P_C2 + O1 + O2))
C1= (1.0 - (O + C2 + C2 + C3 + C4 + I1 + I2 + I3))
beta_n= (9.53333E-5 * exp( - (0.038 * (26.5 + V))))
beta_i= (0.006497 * exp( - (0.03268 * (5.0 + V))))
beta_a= (0.3956 * exp( - (0.06237 * (30.0 + V))))
i_Kto_s = (g_Kto_s * ato_s * ito_s * (V - E_K))
i_Kto_f = (g_Kto_f * (ato_f^3.0) * ito_f * (V - E_K))
alpha_Na3 = (7.0E-7 * exp( - (0.12987012987012986 * (7.0 + V))))
beta_Na5 = (0.02 * alpha_Na3)
beta_Na4 = alpha_Na3
beta_Na3 = (0.0084 + (2.0E-5 * (7.0 + V)))
alpha_Na13 = (3.802 / ((0.1027 * exp( - (0.08333333333333333 * (2.5 + V)))) + (0.25 * exp( - (0.006666666666666667 * (2.5 + V))))))
alpha_Na2 = (1.0 / (0.393956 + (0.188495 * exp( - (0.06024096385542168 * (7.0 + V))))))
beta_Na2 = (alpha_Na13 * alpha_Na2 * alpha_Na3 / (beta_Na13 * beta_Na3))
BJSR =	((1.0 + (CSQN_tot * Km_CSQN / ((Km_CSQN + CaJSR)^2.0)))^(-1.0))
tau_iur=	 (1200.0 - (170.0 / (1.0 + exp((0.17543859649122806 * (45.2 + V))))))
E_CaN=	 (ln((Cao / Cai)) * R * Temp / (2.0 * Fara))
i_Cab=	 (g_Cab * (V - E_CaN))
tau_ti_s=	 (270.0 + (1050.0 / (1.0 + exp((0.17543859649122806 * (45.2 + V))))))
Bss=	 ((1.0 + (CMDN_tot * Km_CMDN / ((Km_CMDN + Cass)^2.0)))^(-1.0))
alpha_Na12=	 (3.802 / ((0.1027 * exp( - (0.06666666666666667 * (2.5 + V)))) + (0.23 * exp( - (0.006666666666666667 * (2.5 + V))))))
alpha_Na11=	 (3.802 / ((0.1027 * exp( - (0.058823529411764705 * (2.5 + V)))) + (0.2 * exp( - (0.006666666666666667 * (2.5 + V))))))
tau_aur=	 (2.058 + (0.493 * exp( - (0.0629 * V))))
beta_a1=	 (6.89E-5 * exp( - (0.04178 * V)))

beta_a0=	 (0.047002 * exp( - (0.0631 * V)))
J_xfer=	 ((Cass - Cai) / tau_xfer)
tau_ta_s=	 (2.058 + (0.493 * exp( - (0.0629 * V))))
J_up= (v3 * (Cai^2.0) / ((Km_up^2.0) + (Cai^2.0)))
alpha_Na5 = (1.0526315789473684E-5 * alpha_Na2)
alpha_Na4 = (0.0010 * alpha_Na2)
iss = (1.0 / (1.0 + exp((0.17543859649122806 * (45.2 + V)))))
Kpcf = (13.0 * (1.0 - exp( - (0.01 * ((14.5 + V)^2.0)))))
C_K0 = (1.0 - (C_K1 + C_K2 + O_K + I_K))
J_tr = ((CaNSR - CaJSR) / tau_tr)
alpha_a1= (0.013733 * exp((0.038198 * V)))
alpha_a0 = (0.022348 * exp((0.01176 * V)))
J_trpn = ( - ((k_minus_htrpn * HTRPN_Ca) + (k_minus_ltrpn * LTRPN_Ca)) + (k_plus_htrpn * Cai * (HTRPN_tot - HTRPN_Ca)) + (k_plus_ltrpn * Cai * (LTRPN_tot - LTRPN_Ca)))
C_Na3 = (1.0 - (O_Na + C_Na1 + C_Na2 + IF_Na + I1_Na + I2_Na + IC_Na2 + IC_Na3))

i_Na=	 (g_Na * O_Na * (V - E_Na))
i_CaL=	 (g_CaL * O * (V - E_CaL))

V'= - ((i_CaL + i_pCa + i_NaCa + i_Cab + i_Na + i_Nab + i_NaK + i_Kto_f + i_Kto_s + i_K1 + i_Ks + i_Kur + i_Kss + i_Kr + i_ClCa - rstar) / Cm)
Cai'=(Bi * ( - (J_up + J_trpn + (( - (2.0 * i_NaCa) + i_Cab + i_pCa) * Acap * Cm / (2.0 * Vmyo * Fara))) + J_leak + J_xfer))
Cass'=(Bss * ((J_rel * VJSR / Vss) - ((J_xfer * Vmyo / Vss) + (i_CaL * Acap * Cm / (2.0 * Vss * Fara)))))
CaJSR'=(BJSR * (J_tr - J_rel))
CaNSR'=(((J_up - J_leak) * Vmyo / VNSR) - (J_tr * VJSR / VNSR))
P_RyR'=( - (0.04 * P_RyR) - (0.1 * exp( - (0.0015432098765432098 * ((-5.0 + V)^2.0))) * i_CaL / i_CaL_max))
LTRPN_Ca'=((k_plus_ltrpn * Cai * (LTRPN_tot - LTRPN_Ca)) - (k_minus_ltrpn * LTRPN_Ca))
HTRPN_Ca'=((k_plus_htrpn * Cai * (HTRPN_tot - HTRPN_Ca)) - (k_minus_htrpn * HTRPN_Ca))
O1'=( - ((k_minus_a * O1) + (k_plus_b * (Cass^m) * O1) + (k_plus_c * O1)) + (k_plus_a * (Cass^n) * P_C1) + (k_minus_b * O2) + (k_minus_c * P_C2))
O2'=((k_plus_b * (Cass^m) * O1) - (k_minus_b * O2))
P_C2'=((k_plus_c * O1) - (k_minus_c * P_C2))
O'=( - ((4.0 * beta * O) + (gamma * O)) + (alpha * C4) + (Kpcb * I1) + (0.0010 * ((alpha * I2) - (Kpcf * O))))
C2'=( - ((beta * C2) + (3.0 * alpha * C2)) + (4.0 * alpha * C1) + (2.0 * beta * C3))
C3'=( - ((2.0 * beta * C3) + (2.0 * alpha * C3)) + (3.0 * alpha * C2) + (3.0 * beta * C4))
C4'=( - ((3.0 * beta * C4) + (alpha * C4) + (gamma * Kpcf * C4)) + (2.0 * alpha * C3) + (4.0 * beta * O) + (0.01 * ((4.0 * Kpcb * beta * I1) - (alpha * gamma * C4))) + (0.0020 * ((4.0 * beta * I2) - (Kpcf * C4))) + (4.0 * beta * Kpcb * I3))
I1'=( - (Kpcb * I1) + (gamma * O) + (0.0010 * ((alpha * I3) - (Kpcf * I1))) + (0.01 * ((alpha * gamma * C4) - (4.0 * beta * Kpcf * I1))))
I2'=( - (gamma * I2) + (0.0010 * ((Kpcf * O) - (alpha * I2))) + (Kpcb * I3) + (0.0020 * ((Kpcf * C4) - (4.0 * beta * I2))))
I3'=( - ((4.0 * beta * Kpcb * I3) + (Kpcb * I3)) + (0.0010 * ((Kpcf * I1) - (alpha * I3))) + (gamma * I2) + (gamma * Kpcf * C4))
Nai'= - (Acap * Cm * (i_Na + i_Nab + (3.0 * i_NaK) + (3.0 * i_NaCa)) / (Vmyo * Fara))
C_Na2'=( - ((beta_Na11 * C_Na2) + (alpha_Na12 * C_Na2) + (beta_Na3 * C_Na2)) + (alpha_Na11 * C_Na3) + (beta_Na12 * C_Na1) + (alpha_Na3 * IC_Na2))
C_Na1'=( - ((beta_Na12 * C_Na1) + (alpha_Na13 * C_Na1) + (beta_Na3 * C_Na1)) + (alpha_Na12 * C_Na2) + (beta_Na13 * O_Na) + (alpha_Na3 * IF_Na))
O_Na'=( - ((beta_Na13 * O_Na) + (alpha_Na2 * O_Na)) + (alpha_Na13 * C_Na1) + (beta_Na2 * IF_Na))
IF_Na'=( - ((beta_Na2 * IF_Na) + (alpha_Na3 * IF_Na) + (alpha_Na4 * IF_Na) + (beta_Na12 * IF_Na)) + (alpha_Na2 * O_Na) + (beta_Na3 * C_Na1) + (beta_Na4 * I1_Na) + (alpha_Na12 * IC_Na2))
I1_Na'=( - ((beta_Na4 * I1_Na) + (alpha_Na5 * I1_Na)) + (alpha_Na4 * IF_Na) + (beta_Na5 * I2_Na))
I2_Na'=((alpha_Na5 * I1_Na) - (beta_Na5 * I2_Na))
IC_Na2'=( - ((beta_Na11 * IC_Na2) + (alpha_Na12 * IC_Na2) + (alpha_Na3 * IC_Na2)) + (alpha_Na11 * IC_Na3) + (beta_Na12 * IF_Na) + (beta_Na3 * IC_Na2))
IC_Na3'=( - ((alpha_Na11 * IC_Na3) + (alpha_Na3 * IC_Na3)) + (beta_Na11 * IC_Na2) + (beta_Na3 * C_Na3))
Ki'= - (Acap * Cm * (i_Kto_f + i_Kto_s + i_K1 + i_Ks + i_Kss + i_Kur + i_Kr + (2.0 * i_NaK)) / (Vmyo * Fara))
ato_f'=((alpha_a * (1.0 - ato_f)) - (beta_a * ato_f))
ito_f'=((alpha_i1 * (1.0 - ito_f)) - (beta_i1 * ito_f))
ato_s' = ((ass - ato_s) / tau_ta_s)
ito_s' = ((iss - ito_s) / tau_ti_s)
nKs'=((alpha_n * (1.0 - nKs)) - (beta_n * nKs))
aur'=((ass - aur) / tau_aur)
iur'=((iss - iur) / tau_iur)
aKss'=((ass - aKss) / tau_Kss)
iKss'=0.0
C_K2'=( - ((kb * C_K2) + (alpha_a1 * C_K2)) + (kf * C_K1) + (beta_a1 * O_K))
C_K1'=( - ((beta_a0 * C_K1) + (kf * C_K1)) + (alpha_a0 * C_K0) + (kb * C_K2))
O_K'=( - ((beta_a1 * O_K) + (alpha_i * O_K)) + (alpha_a1 * C_K2) + (beta_i * I_K))
I_K'=((alpha_i * O_K) - (beta_i * I_K))

aux ina=i_Na
aux ical=i_CaL
aux ikur=i_Kur

# Numerical and plotting parameters for xpp
@ maxstor=5000000, bounds=10000000, total=1100, xp=t, yp=V, trans=860
@ meth=Gear, dt=0.05, toler=0.01, xlo=860, xhi=1100, ylo=-100, yhi=60     

done
This file has been produced by sbmlutils.

Terms of use

Copyright © 2017 Matthias Koenig

Redistribution and use of any part of this model, with or without modification, are permitted provided that the following conditions are met:

  1. Redistributions of this SBML file must retain the above copyright notice, this list of conditions and the following disclaimer.
  2. Redistributions in a different form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
This model is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.


Model :

id
name
time
substance
extent
volume
area
length
Access SBML model  L3V1

FunctionDefinitions [4] name math sbo cvterm
max minimum x y x x y y
min maximum x y x x y y
heav heavyside x 0 x 0 0.5 x 0 1 x 0 0
mod modulo x y x y x y x 0 y 0 x y x y

Parameters [190] name constant value unit derived unit sbo cvterm
temp temp = 298.0 298.0 None
r r = 8.314 8.314 None
fara fara = 96.5 96.5 None
e_cal e_cal = 63.0 63.0 None
e_cl e_cl = -40.0 -40.0 None
cm cm = 1.0 1.0 None
vmyo vmyo = 2.584e-5 2.584e-05 None
vnsr vnsr = 2.098e-6 2.098e-06 None
vjsr vjsr = 1.2e-7 1.2e-07 None
v v = -82.4202 -82.4202 None
cai cai = 0.115001 0.115001 None
nai nai = 14237.1 14237.1 None
ki ki = 143720.0 143720.0 None
c_na1 c_na1 = 2.79132e-4 0.000279132 None
c_na2 c_na2 = 0.020752 0.020752 None
o_na o_na = 7.13483e-7 7.13483e-07 None
if_na if_na = 1.53176e-4 0.000153176 None
i1_na i1_na = 6.73345e-7 6.73345e-07 None
i2_na i2_na = 1.55787e-9 1.55787e-09 None
ic_na2 ic_na2 = 0.0113879 0.0113879 None
ic_na3 ic_na3 = 0.34278 0.34278 None
ato_f ato_f = 0.00265563 0.00265563 None
ito_f ito_f = 0.999977 0.999977 None
ato_s ato_s = 4.17069e-4 0.000417069 None
ito_s ito_s = 0.998543 0.998543 None
nks nks = 2.62753e-4 0.000262753 None
aur aur = 4.17069e-4 0.000417069 None
iur iur = 0.998543 0.998543 None
akss akss = 4.17069e-4 0.000417069 None
ikss ikss = 1.0 1.0 None
c_k1 c_k1 = 9.92513e-4 0.000992513 None
c_k2 c_k2 = 6.41229e-4 0.000641229 None
o_k o_k = 1.75298e-4 0.000175298 None
i_k i_k = 3.19129e-5 3.19129e-05 None
cass cass = 0.115001 0.115001 None
cajsr cajsr = 1299.5 1299.5 None
cansr cansr = 1299.5 1299.5 None
p_ryr p_ryr = 0.0 0.0 None
ltrpn_ca ltrpn_ca = 11.2684 11.2684 None
htrpn_ca htrpn_ca = 125.29 125.29 None
o o = 9.30308e-19 9.30308e-19 None
o1 o1 = 1.49102e-5 1.49102e-05 None
o2 o2 = 9.51726e-11 9.51726e-11 None
p_c2 p_c2 = 1.6774e-4 0.00016774 None
c2 c2 = 1.24216e-4 0.000124216 None
c3 c3 = 5.78679e-9 5.78679e-09 None
c4 c4 = 1.19816e-13 1.19816e-13 None
i1 i1 = 4.97023e-19 4.97023e-19 None
i2 i2 = 3.45847e-14 3.45847e-14 None
i3 i3 = 1.85106e-14 1.85106e-14 None
r0 r0 = 1 1.0 None
period period = 300 300.0 None
pulse pulse = 10 10.0 None
tf tf = 0 0.0 None
tp tp = 5 5.0 None
tstart tstart = 10 10.0 None
nao nao = 140000.0 140000.0 None
ltrpn_tot ltrpn_tot = 70.0 70.0 None
kf kf = 0.023761 0.023761 None
kb kb = 0.036778 0.036778 None
km_csqn km_csqn = 800.0 800.0 None
i_pca_max i_pca_max = 1.0 1.0 None
cmdn_tot cmdn_tot = 50.0 50.0 None
kpc_max kpc_max = 0.23324 0.23324 None
km_cmdn km_cmdn = 0.238 0.238 None
k_mca k_mca = 1380.0 1380.0 None
km_nai km_nai = 21.0 21.0 None
v1 v1 = 4.5 4.5 None
v2 v2 = 1.7e-5 1.7e-05 None
v3 v3 = 0.45 0.45 None
k_plus_c k_plus_c = 0.0090 0.009 None
k_minus_ltrpn k_minus_ltrpn = 0.196 0.196 None
g_cab g_cab = 3.67e-4 0.000367 None
k_plus_b k_plus_b = 0.00405 0.00405 None
k_plus_a k_plus_a = 0.006075 0.006075 None
km_up km_up = 0.5 0.5 None
k_mna k_mna = 87500.0 87500.0 None
km_cl km_cl = 10.0 10.0 None
csqn_tot csqn_tot = 15000.0 15000.0 None
k_minus_htrpn k_minus_htrpn = 3.2e-5 3.2e-05 None
km_pca km_pca = 0.5 0.5 None
k_plus_ltrpn k_plus_ltrpn = 0.0327 0.0327 None
acap acap = 1.534e-4 0.0001534 None
g_na g_na = 13.0 13.0 None
g_cal g_cal = 1.51729 1.51729 None
g_ks g_ks = 0.000575 0.000575 None
g_kr g_kr = 0.0078 0.0078 None
g_kur g_kur = 0.0016 0.0016 None
g_clca g_clca = 10.0 10.0 None
vss vss = 1.485e-9 1.485e-09 None
k_naca k_naca = 992.8 992.8 None
k_plus_htrpn k_plus_htrpn = 0.00237 0.00237 None
ko ko = 5400.0 5400.0 None
cao cao = 1800.0 1800.0 None
tau_xfer tau_xfer = 8.0 8.0 None
k_sat k_sat = 0.1 0.1 None
n n = 4.0 4.0 None
g_kto_f g_kto_f = 0.4067 0.4067 None
g_kto_s g_kto_s = 0.01 0.01 None
g_nab g_nab = 0.0026 0.0026 None
m m = 3.0 3.0 None
km_ko km_ko = 1.5 1.5 None
g_kss g_kss = 0.05 0.05 None
htrpn_tot htrpn_tot = 140.0 140.0 None
i_nak_max i_nak_max = 0.88 0.88 None
i_cal_max i_cal_max = 7.0 7.0 None
tau_tr tau_tr = 20.0 20.0 None
k_minus_c k_minus_c = 8.0e-4 0.0008 None
k_minus_b k_minus_b = 0.965 0.965 None
k_minus_a k_minus_a = 0.07125 0.07125 None
kpc_half kpc_half = 20.0 20.0 None
eta eta = 0.35 0.35 None
kpcb kpcb = 5.0e-4 0.0005 None
ts 0.0 dimensionless None
rstar 0.0 dimensionless None
ass 0.0 dimensionless None
e_k 0.0 dimensionless None
i_ks 0.0 dimensionless None
e_na 0.0 dimensionless None
i_kr 0.0 dimensionless None
i_kur 0.0 dimensionless None
sigma 0.0 dimensionless None
f_nak 0.0 dimensionless None
j_leak 0.0 dimensionless None
i_naca 0.0 dimensionless None
o_clca 0.0 dimensionless None
i_clca 0.0 dimensionless None
i_nab 0.0 dimensionless None
tau_kss 0.0 dimensionless None
beta_i1 0.0 dimensionless None
beta 0.0 dimensionless None
i_kss 0.0 dimensionless None
i_k1 0.0 dimensionless None
i_nak 0.0 dimensionless None
j_rel 0.0 dimensionless None
alpha_n 0.0 dimensionless None
alpha_i 0.0 dimensionless None
alpha_a 0.0 dimensionless None
i_pca 0.0 dimensionless None
alpha_i1 0.0 dimensionless None
beta_na13 0.0 dimensionless None
beta_na12 0.0 dimensionless None
beta_na11 0.0 dimensionless None
alpha 0.0 dimensionless None
bi 0.0 dimensionless None
gamma 0.0 dimensionless None
p_c1 0.0 dimensionless None
c1 0.0 dimensionless None
beta_n 0.0 dimensionless None
beta_i 0.0 dimensionless None
beta_a 0.0 dimensionless None
i_kto_s 0.0 dimensionless None
i_kto_f 0.0 dimensionless None
alpha_na3 0.0 dimensionless None
beta_na5 0.0 dimensionless None
beta_na4 0.0 dimensionless None
beta_na3 0.0 dimensionless None
alpha_na13 0.0 dimensionless None
alpha_na2 0.0 dimensionless None
beta_na2 0.0 dimensionless None
bjsr 0.0 dimensionless None
tau_iur 0.0 dimensionless None
e_can 0.0 dimensionless None
i_cab 0.0 dimensionless None
tau_ti_s 0.0 dimensionless None
bss 0.0 dimensionless None
alpha_na12 0.0 dimensionless None
alpha_na11 0.0 dimensionless None
tau_aur 0.0 dimensionless None
beta_a1 0.0 dimensionless None
beta_a0 0.0 dimensionless None
j_xfer 0.0 dimensionless None
tau_ta_s 0.0 dimensionless None
j_up 0.0 dimensionless None
alpha_na5 0.0 dimensionless None
alpha_na4 0.0 dimensionless None
iss 0.0 dimensionless None
kpcf 0.0 dimensionless None
c_k0 0.0 dimensionless None
j_tr 0.0 dimensionless None
alpha_a1 0.0 dimensionless None
alpha_a0 0.0 dimensionless None
j_trpn 0.0 dimensionless None
c_na3 0.0 dimensionless None
i_na 0.0 dimensionless None
i_cal 0.0 dimensionless None
ina 0.0 dimensionless None
ical 0.0 dimensionless None
ikur 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [118]   assignment name derived units sbo cvterm
d v/dt = i_cal i_pca i_naca i_cab i_na i_nab i_nak i_kto_f i_kto_s i_k1 i_ks i_kur i_kss i_kr i_clca rstar cm None
d cai/dt = bi j_up j_trpn 2 i_naca i_cab i_pca acap cm 2 vmyo fara j_leak j_xfer None
d cass/dt = bss j_rel vjsr vss j_xfer vmyo vss i_cal acap cm 2 vss fara None
d cajsr/dt = bjsr j_tr j_rel None
d cansr/dt = j_up j_leak vmyo vnsr j_tr vjsr vnsr None
d p_ryr/dt = 0.04 p_ryr 0.1 0.00154320987654321 5 v 2 i_cal i_cal_max None
d ltrpn_ca/dt = k_plus_ltrpn cai ltrpn_tot ltrpn_ca k_minus_ltrpn ltrpn_ca None
d htrpn_ca/dt = k_plus_htrpn cai htrpn_tot htrpn_ca k_minus_htrpn htrpn_ca None
d o1/dt = k_minus_a o1 k_plus_b cass m o1 k_plus_c o1 k_plus_a cass n p_c1 k_minus_b o2 k_minus_c p_c2 None
d o2/dt = k_plus_b cass m o1 k_minus_b o2 None
d p_c2/dt = k_plus_c o1 k_minus_c p_c2 None
d o/dt = 4 beta o gamma o alpha c4 kpcb i1 0.001 alpha i2 kpcf o None
d c2/dt = beta c2 3 alpha c2 4 alpha c1 2 beta c3 None
d c3/dt = 2 beta c3 2 alpha c3 3 alpha c2 3 beta c4 None
d c4/dt = 3 beta c4 alpha c4 gamma kpcf c4 2 alpha c3 4 beta o 0.01 4 kpcb beta i1 alpha gamma c4 0.002 4 beta i2 kpcf c4 4 beta kpcb i3 None
d i1/dt = kpcb i1 gamma o 0.001 alpha i3 kpcf i1 0.01 alpha gamma c4 4 beta kpcf i1 None
d i2/dt = gamma i2 0.001 kpcf o alpha i2 kpcb i3 0.002 kpcf c4 4 beta i2 None
d i3/dt = 4 beta kpcb i3 kpcb i3 0.001 kpcf i1 alpha i3 gamma i2 gamma kpcf c4 None
d nai/dt = acap cm i_na i_nab 3 i_nak 3 i_naca vmyo fara None
d c_na2/dt = beta_na11 c_na2 alpha_na12 c_na2 beta_na3 c_na2 alpha_na11 c_na3 beta_na12 c_na1 alpha_na3 ic_na2 None
d c_na1/dt = beta_na12 c_na1 alpha_na13 c_na1 beta_na3 c_na1 alpha_na12 c_na2 beta_na13 o_na alpha_na3 if_na None
d o_na/dt = beta_na13 o_na alpha_na2 o_na alpha_na13 c_na1 beta_na2 if_na None
d if_na/dt = beta_na2 if_na alpha_na3 if_na alpha_na4 if_na beta_na12 if_na alpha_na2 o_na beta_na3 c_na1 beta_na4 i1_na alpha_na12 ic_na2 None
d i1_na/dt = beta_na4 i1_na alpha_na5 i1_na alpha_na4 if_na beta_na5 i2_na None
d i2_na/dt = alpha_na5 i1_na beta_na5 i2_na None
d ic_na2/dt = beta_na11 ic_na2 alpha_na12 ic_na2 alpha_na3 ic_na2 alpha_na11 ic_na3 beta_na12 if_na beta_na3 ic_na2 None
d ic_na3/dt = alpha_na11 ic_na3 alpha_na3 ic_na3 beta_na11 ic_na2 beta_na3 c_na3 None
d ki/dt = acap cm i_kto_f i_kto_s i_k1 i_ks i_kss i_kur i_kr 2 i_nak vmyo fara None
d ato_f/dt = alpha_a 1 ato_f beta_a ato_f None
d ito_f/dt = alpha_i1 1 ito_f beta_i1 ito_f None
d ato_s/dt = ass ato_s tau_ta_s None
d ito_s/dt = iss ito_s tau_ti_s None
d nks/dt = alpha_n 1 nks beta_n nks None
d aur/dt = ass aur tau_aur None
d iur/dt = iss iur tau_iur None
d akss/dt = ass akss tau_kss None
d ikss/dt = 0 None
d c_k2/dt = kb c_k2 alpha_a1 c_k2 kf c_k1 beta_a1 o_k None
d c_k1/dt = beta_a0 c_k1 kf c_k1 alpha_a0 c_k0 kb c_k2 None
d o_k/dt = beta_a1 o_k alpha_i o_k alpha_a1 c_k2 beta_i i_k None
d i_k/dt = alpha_i o_k beta_i i_k None
ts = t tstart None
rstar = r0 pulse heav mod ts period tf heav mod ts period tf tp None
ass = 1 1 0.12987012987013 22.5 v None
e_k = ko ki r temp fara None
i_ks = g_ks nks 2 v e_k None
e_na = 0.9 nao 0.1 ko 0.9 nai 0.1 ki r temp fara None
i_kr = g_kr o_k v 0.98 ko 0.02 nao 0.98 ki 0.02 nai r temp fara None
i_kur = g_kur aur iur v e_k None
sigma = 0.142857142857143 1 1.48588410104012 -5 nao None
f_nak = 1 1 0.1245 0.1 v fara r temp 0.0365 sigma fara v r temp None
j_leak = v2 cansr cai None
i_naca = k_naca eta v fara r temp nai 3 cao 1 eta v fara r temp nao 3 cai k_mna 3 nao 3 k_mca cao 1 k_sat 1 eta v fara r temp None
o_clca = 0.2 1 0.128205128205128 46.7 v None
i_clca = g_clca o_clca v e_cl cai cai km_cl None
i_nab = g_nab v e_na None
tau_kss = 13.17 39.3 0.0862 v None
beta_i1 = 9.5 -4 0.142857142857143 33.5 v 1 0.051335 0.142857142857143 33.5 v None
beta = 0.05 0.0769230769230769 12 v None
i_kss = g_kss akss ikss v e_k None
i_k1 = 0.2938 ko v e_k 210 ko 1 0.0896 v e_k None
i_nak = i_nak_max f_nak ko 1 km_nai nai 0.66667 ko km_ko None
j_rel = v1 o1 o2 cajsr cass p_ryr None
alpha_n = 4.81333 -6 26.5 v 1 0.128 26.5 v None
alpha_i = 0.090821 0.023391 5 v None
alpha_a = 0.18064 0.03577 30 v None
i_pca = i_pca_max cai 2 km_pca 2 cai 2 None
alpha_i1 = 1.52 -4 0.142857142857143 13.5 v 1 0.0067083 0.142857142857143 33.5 v None
beta_na13 = 0.22 0.0492610837438424 7.5 v None
beta_na12 = 0.2 0.0492610837438424 2.5 v None
beta_na11 = 0.1917 0.0492610837438424 2.5 v None
alpha = 0.4 0.1 12 v 1 0.75 0.0025 20 v 2 0.7 0.1 40 v 2 1 0.12 0.1 12 v None
bi = 1 cmdn_tot km_cmdn km_cmdn cai 2 1 None
gamma = kpc_max cass kpc_half cass None
p_c1 = 1 p_c2 o1 o2 None
c1 = 1 o c2 c2 c3 c4 i1 i2 i3 None
beta_n = 9.53333 -5 0.038 26.5 v None
beta_i = 0.006497 0.03268 5 v None
beta_a = 0.3956 0.06237 30 v None
i_kto_s = g_kto_s ato_s ito_s v e_k None
i_kto_f = g_kto_f ato_f 3 ito_f v e_k None
alpha_na3 = 7 -7 0.12987012987013 7 v None
beta_na5 = 0.02 alpha_na3 None
beta_na4 = alpha_na3 None
beta_na3 = 0.0084 2 -5 7 v None
alpha_na13 = 3.802 0.1027 0.0833333333333333 2.5 v 0.25 0.00666666666666667 2.5 v None
alpha_na2 = 1 0.393956 0.188495 0.0602409638554217 7 v None
beta_na2 = alpha_na13 alpha_na2 alpha_na3 beta_na13 beta_na3 None
bjsr = 1 csqn_tot km_csqn km_csqn cajsr 2 1 None
tau_iur = 1200 170 1 0.175438596491228 45.2 v None
e_can = cao cai r temp 2 fara None
i_cab = g_cab v e_can None
tau_ti_s = 270 1050 1 0.175438596491228 45.2 v None
bss = 1 cmdn_tot km_cmdn km_cmdn cass 2 1 None
alpha_na12 = 3.802 0.1027 0.0666666666666667 2.5 v 0.23 0.00666666666666667 2.5 v None
alpha_na11 = 3.802 0.1027 0.0588235294117647 2.5 v 0.2 0.00666666666666667 2.5 v None
tau_aur = 2.058 0.493 0.0629 v None
beta_a1 = 6.89 -5 0.04178 v None
beta_a0 = 0.047002 0.0631 v None
j_xfer = cass cai tau_xfer None
tau_ta_s = 2.058 0.493 0.0629 v None
j_up = v3 cai 2 km_up 2 cai 2 None
alpha_na5 = 1.05263157894737 -5 alpha_na2 None
alpha_na4 = 0.001 alpha_na2 None
iss = 1 1 0.175438596491228 45.2 v None
kpcf = 13 1 0.01 14.5 v 2 None
c_k0 = 1 c_k1 c_k2 o_k i_k None
j_tr = cansr cajsr tau_tr None
alpha_a1 = 0.013733 0.038198 v None
alpha_a0 = 0.022348 0.01176 v None
j_trpn = k_minus_htrpn htrpn_ca k_minus_ltrpn ltrpn_ca k_plus_htrpn cai htrpn_tot htrpn_ca k_plus_ltrpn cai ltrpn_tot ltrpn_ca None
c_na3 = 1 o_na c_na1 c_na2 if_na i1_na i2_na ic_na2 ic_na3 None
i_na = g_na o_na v e_na None
i_cal = g_cal o v e_cal None
ina = i_na None
ical = i_cal None
ikur = i_kur None
t = time None