XPP model

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

# Cardiac ventricular action potential (Luo Rudy dynamic model) (Luo Rudy 1994) used 
# in (Wang et al 2006)
# LRdt01.ode
# 

# Initial Conditions
init V=-85.2119207874627, Cai=0.000117482029668194
init m=0.00149183115674257, h=0.985596581239651, j=0.990898461370389
init d=5.82597094505446e-6
init f=0.997765362821995
init b=0.00136737866785149
init g=0.98881442877378
init xr=0.000204700363126417
init xs1=0.00660746743356887
init xs2=0.0303768241233812
init zdv=0.0144622472219576
init ydv=0.999945568566232
init Ca_JSR=1.12791401197882, Ca_NSR=1.76731003671612
init APtrack=9.65910542308504e-196, APtrack2=5.33944967562997e-195, APtrack3=0.000129515197402902
init CaFluxtr=2.69380318286645e-196
init O_track=0, O_track2=0, O_track3=0
init Na_i=13.3649235394859, Ki=141.056872392446

# Stimulation Protocol
Par period=200
Par pulse=25.5
Par tf=0, tp=2, tstart=50
ts=t-tstart
Iapp = -pulse*(heav(mod(ts,period)-tf)-heav(mod(ts,period)-(tf+tp)))

# Constant Variables
Number Rgas=8314, Temp=310, Fara=96485
Number delta_m=1e-5
Number gg_Nai=0.75, gg_Nao=0.75, gg_Ki=0.75, gg_Ko=0.75, gg_Cai=1, gg_Cao=0.341
Number Km_Ca=0.0006
Number g_CaT=0.05
Number K_mpCa=0.0005
Number I_pCa=1.15
Number II_NaK=2.25
Number K_mNai=10
Number K_mKo=1.5
Number K_m_ns_Ca=0.0012
Number A_cap=0.0001534
Number G_rel_max=60
Number Grelover=4
Number tau_tr=180
Number K_mrel=0.0008
Number delCaith=0.00018
Number CSQN_max=10
Number KmCSQN=0.8
Number K_mup=0.00092
Number Iup=0.00875
Number CaNSRmax=15
Number K_mTn=0.0005
Number KmCMDN=0.00238
Number Tn_max=0.07
Number CMDNmax=0.05
Number CSQN_Th=0.7, Log_Th=0.98

# Parameters
Par tauT1=0.5, tauT2=0.5
Par GG_Kr=0.02614, GG_Ks=0.433
Par Nao=140, Ko=5.4, Cao=1.8
Par g_Na=16
Par P_Ca=0.00054, P_Na=6.75e-7, P_K=1.93e-7
Par PNaK=0.01833
Par GG_K1=0.75
Par g_Kp=0.00552
Par g_K_Na=0.12848
Par nKNa=2.8
Par kdKNa=66
Par iKATP_on=1
Par nATP=0.24
Par N_area=5e-5
Par ATPi=3
Par hATP=2
Par kATP=0.00025
Par g_Nab=0.004
Par g_Cab=0.003016
Par P_ns_Ca=1.75e-7
Par c1=0.00025, c2=0.0001, gamma=0.15

# Functions
E_Na = Rgas*Temp/Fara*ln(Nao/Na_i)
alpha_m=heav(abs(V+47.13)-delta_m)*(0.32*(V+47.13)/(1-exp(-0.1*(V+47.13))))+heav(delta_m-abs(V+47.13))*3.2
beta_m = 0.08*exp(-V/11)
alpha_h=heav(-40-V)*(0.135*exp((80+V)/(-6.8)))+heav(V-(-40))*0
beta_h=heav(-40-V)*(3.56*exp(0.079*V)+310000*exp(0.35*V))+heav(V-(-40))*(1/(0.13*(1+exp((V+10.66)/(-11.1)))))
alpha_j=heav(-40-V)*(-(127140*exp(0.2444*V)+3.474e-5*exp(-0.04391*V))*(V+37.78)/(1+exp(0.311*(V+79.23))))+heav(V-(-40))*0
beta_j= heav(-V-40)*((0.1212  * exp( - (0.01052 * V)) / (1.0 + exp( - (0.1378 * (40.14 + V)))))) + heav(V+40)*((0.3 * exp( - (2.535E-7 * V)) / (1.0 + exp( - (0.1 * (32.0 + V))))))
d_inf = 1/(1+exp(-(V+10)/6.24))
tau_d=heav(1e-5-abs(V+10))*1/(0.035*6.24*2)+heav(abs(V+10)-1e-5)*1*d_inf*(1-exp(-(V+10)/6.24))/(0.035*(V+10))
alpha_d = d_inf/tau_d
beta_d = (1-d_inf)/tau_d
f_inf = 1/(1+exp((V+35.06)/8.6))+0.6/(1+exp((50-V)/20))
tau_f = 1/(0.0197*exp(-(0.0337*(V+10))^2)+0.02)
alpha_f = f_inf/tau_f
beta_f = (1-f_inf)/tau_f
f' = alpha_f*(1-f)-beta_f*f
f_Ca = 1/(1+Cai/Km_Ca)
b_inf = 1/(1+exp(-(V+14)/10.8))
tau_b = 3.7+6.1/(1+exp((V+25)/4.5))
b' = (b_inf-b)/tau_b
g_inf = 1/(1+exp((V+60)/5.6))
tau_g=heav(0-V)*(-0.875*V+12)+heav(V-0)*12
g_Kr = GG_Kr*(Ko/5.4)^(0.5)
Rect = 1/(1+exp((V+9)/22.4))
xr_inf = 1/(1+exp(-(V+21.5)/7.5))
tau_xr = 1/(0.00138*(V+14.2)/(1-exp(-0.123*(V+14.2)))+0.00061*(V+38.9)/(exp(0.145*(V+38.9))-1))
E_Ks = Rgas*Temp/Fara*ln((Ko+PNaK*Nao)/(Ki+PNaK*Na_i))
g_Ks = GG_Ks*(1+0.6/(1+((3.8e-5/Cai)^1.4)))
xs1_inf = 1/(1+exp(-(V-1.5)/16.7))
tau_xs1 = 1/(7.19e-5*(V+30)/(1-exp(-0.148*(V+30)))+0.000131*(V+30)/(exp(0.0687*(V+30))-1))
xs2_inf = 1/(1+exp(-(V-1.5)/16.7))
tau_xs2 = 4/(7.19e-5*(V+30)/(1-exp(-0.148*(V+30)))+0.000131*(V+30)/(exp(0.0687*(V+30))-1))
g_K1 = GG_K1*(Ko/5.4)^(0.5)
E_K = Rgas*Temp/Fara*ln(Ko/Ki)
alpha_K1 = 1.02/(1+exp(0.2385*(V-E_K-59.215)))
beta_K1 = 1*(0.49124*exp(0.08032*(V-E_K+5.476))+exp(0.06175*(V-E_K-594.31)))/(1+exp(-0.5143*(V-E_K+4.753)))
K1_inf = alpha_K1/(alpha_K1+beta_K1)
Kp = 1/(1+exp((7.488-V)/5.98))
i_Kp = g_Kp*Kp*(V-E_K)   
pona = 0.85/(1+((kdKNa/Na_i)^nKNa))
pov = 0.8-0.65/(1+exp((V+125)/15))
g_K_ATP = iKATP_on*0.000193/N_area
pATP = 1/(1+((ATPi/kATP)^hATP))
GKbarT = g_K_ATP*pATP*((Ko/4)^nATP)
g_to = 0*0.5
rvdv = exp(V/100)
alpha_zdv = 10*exp((V-40)/25)/(1+exp((V-40)/25))
beta_zdv = 10*exp(-(V+90)/25)/(1+exp(-(V+90)/25))
tau_zdv = 1/(alpha_zdv+beta_zdv)
zdv_ss = alpha_zdv/(alpha_zdv+beta_zdv)
alpha_ydv = 0.015/(1+exp((V+60)/5))
beta_ydv = 0.1*exp((V+25)/5)/(1+exp((V+25)/5))
tau_ydv = 1/(alpha_ydv+beta_ydv)
ydv_ss = alpha_ydv/(alpha_ydv+beta_ydv)
E_Ca = Rgas*Temp/(2*Fara)*ln(Cao/Cai)
sigma = 1/7*(exp(Nao/67.3)-1)
f_NaK = 1/(1+0.1245*exp(-0.1*V*Fara/(Rgas*Temp))+0.0365*sigma*exp(-V*Fara/(Rgas*Temp)))
Ins_Na = P_ns_Ca*(1^2)*V*(Fara^2)/(Rgas*Temp)*(gg_Nai*Na_i*exp(1*V*Fara/(Rgas*Temp))-gg_Nao*Nao)/(exp(1*V*Fara/(Rgas*Temp))-1)
Ins_K = P_ns_Ca*(1^2)*V*(Fara^2)/(Rgas*Temp)*(gg_Ki*Ki*exp(1*V*Fara/(Rgas*Temp))-gg_Ko*Ko)/(exp(1*V*Fara/(Rgas*Temp))-1)
Vmyo = 0.68*pi*0.1*(0.011^2)
VJSR = 0.0048*pi*0.1*(0.011^2)
VNSR = 0.0552*pi*0.1*(0.011^2)
Kleak = Iup/CaNSRmax

# Currents
i_Na = g_Na*(m^3)*h*j*(V-E_Na)
ICaCa = (P_Ca*(2^2)*V*(Fara^2)/(Rgas*Temp))*(gg_Cai*Cai*exp(2*V*Fara/(Rgas*Temp))-gg_Cao*Cao)/(exp(2*V*Fara/(Rgas*Temp))-1)
ICaNa = P_Na*(1^2)*V*(Fara^2)/(Rgas*Temp)*(gg_Nai*Na_i*exp(1*V*Fara/(Rgas*Temp))-gg_Nao*Nao)/(exp(1*V*Fara/(Rgas*Temp))-1)
ICaK = P_K*(1^2)*V*(Fara^2)/(Rgas*Temp)*(gg_Ki*Ki*exp(1*V*Fara/(Rgas*Temp))-gg_Ko*Ko)/(exp(1*V*Fara/(Rgas*Temp))-1)
i_CaCa = d*f*f_Ca*ICaCa
i_CaNa = d*f*f_Ca*ICaNa
i_CaK = d*f*f_Ca*ICaK
iCaL = i_CaCa+i_CaK+i_CaNa
iCaT = g_CaT*b*b*g*(V-E_Ca)
i_Kr = g_Kr*xr*Rect*(V-E_K)
i_Ks = g_Ks*xs1*xs2*(V-E_Ks)
i_K1 = g_K1*K1_inf*(V-E_K)
iKNa = g_K_Na*pona*pov*(V-E_K)
iKATP = GKbarT*(V-E_K)
i_to = g_to*(zdv^3)*ydv*rvdv*(V-E_K)
ipCa = I_pCa*Cai/(K_mpCa+Cai)
i_Ca_b = g_Cab*(V-E_Ca)
iNab = g_Nab*(V-E_Na)
iNaK = II_NaK*f_NaK*1/(1+(K_mNai/Na_i)^2)*Ko/(Ko+K_mKo)
i_NaCa = c1*exp((gamma-1)*V*Fara/(Rgas*Temp))*(exp(V*Fara/(Rgas*Temp))*(Na_i^3)*Cao-(Nao^3)*Cai)/(1+c2*exp((gamma-1)*V*Fara/(Rgas*Temp))*(exp(V*Fara/(Rgas*Temp))*(Na_i^3)*Cao+(Nao^3)*Cai))
insNa = Ins_Na*1/(1+((K_m_ns_Ca/Cai)^3))
insK = Ins_K*1/(1+((K_m_ns_Ca/Cai)^3))
insCa = insNa+insK
i_rel = G_relVis*(Ca_JSR-Cai)
i_up = Iup*Cai/(Cai+K_mup)
i_leak = Kleak*Ca_NSR
i_tr = (Ca_NSR-Ca_JSR)/tau_tr

# Control of Calcium Fluxes During an Action Potential
APtrack' = heav(dV_dt-150)*(100*(1-APtrack)-tauT1*APtrack)+  \
heav(150-dV_dt)*(-tauT2*APtrack)
APtrack2' = heav(0.2-APtrack)*heav(APtrack-0.18)*(100*(1-APtrack2)-0.5*APtrack2)+heav(APtrack-0.2)*(-0.5*APtrack2)+heav(0.18-APtrack)*(-0.5*APtrack2)  
APtrack3' = heav(0.2-APtrack)*heav(APtrack-0.18)*(100*(1-APtrack3)-0.5*APtrack3)+heav(APtrack-0.2)*(-0.01*APtrack3)+heav(0.18-APtrack)*(-0.01*APtrack3)
CaFluxtr'=heav(APtrack-0.2)*(-1*A_cap*(i_CaCa+iCaT-2*i_NaCa+ipCa+i_Ca_b)/(2*Vmyo*Fara))+heav(APtrack2-0.01)*heav(0.2-APtrack)*0+ \
heav(0.01-APtrack2)*(-0.5*CaFluxtr)
O_track'=heav(1/(1+KmCSQN/Ca_JSR)-CSQN_Th)*heav(0.37-O_track3)*heav(0.37-APtrack3)*(50*(1-O_track))+ \
heav(CSQN_Th-1/(1+KmCSQN/Ca_JSR))*(-0.5*O_track)+heav(O_track3-0.37)*(-0.5*O_track)+heav(APtrack3-0.37)*(-0.5*O_track)
O_track2'=heav(O_track - Log_Th)*heav(Log_Th- O_track2)*(50*(1-O_track2))+heav(Log_Th-O_track)*(-0.5*O_track2)+heav(O_track2-Log_Th)*(-0.5*O_track2)
O_track3'=heav(O_track - Log_Th)*heav(Log_Th-O_track3)*(50*(1-O_track3))+heav(Log_Th-O_track)*(-0.01*O_track3)+heav(O_track3-Log_Th)*(-0.01*O_track3)
G_relVis=heav(CaFluxtr-delCaith)*((G_rel_max*(CaFluxtr-delCaith)/(K_mrel+CaFluxtr-delCaith))*(1-APtrack2)*APtrack2)+ \
 heav(delCaith-CaFluxtr)*heav(O_track2-0)*(Grelover*(1-O_track2)*O_track2)+heav(delCaith-CaFluxtr)*0+heav(CaFluxtr-delCaith)*0+heav(0-O_track2)*0

# Differential Wquations
dV_dt = -(i_Na+iCaL+iCaT+i_Kr+i_Ks+iKNa+iKATP+i_to+i_K1+i_Kp+i_NaCa+ipCa+iNab+i_Ca_b+iNaK+insCa+Iapp)
V' = dV_dt
m' = alpha_m*(1-m)-beta_m*m
h' = alpha_h*(1-h)-beta_h*h
j' = alpha_j*(1-j)-beta_j*j
d' = alpha_d*(1-d)-beta_d*d
g' = (g_inf-g)/tau_g
xr' = (xr_inf-xr)/tau_xr
xs1' = (xs1_inf-xs1)/tau_xs1
xs2' = (xs2_inf-xs2)/tau_xs2
zdv' = (zdv_ss-zdv)/tau_zdv
ydv' = (ydv_ss-ydv)/tau_ydv
Cai' = 1/(1+CMDNmax*KmCMDN/(KmCMDN+Cai)^2+ \
Tn_max*K_mTn/(K_mTn+Cai)^2)*(-1*A_cap*(i_CaCa+iCaT-2*i_NaCa+ipCa+i_Ca_b)/(2*Vmyo*Fara)+i_rel*VJSR/Vmyo+(i_leak-i_up)*VNSR/Vmyo)
Ca_JSR' = (1/(1+CSQN_max*KmCSQN/(KmCSQN+Ca_JSR)^2))*(i_tr-i_rel)
Ca_NSR' = -i_tr*VJSR/VNSR-i_leak+i_up
Na_i' = -1*(i_Na+i_CaNa+iNab+insNa+i_NaCa*3+iNaK*3)*A_cap/(Vmyo*Fara)
Ki' = -1*(-Iapp+i_CaK+i_Kr+i_Ks+i_K1+i_Kp+iKNa+iKATP+i_to+insK+(-iNaK)*2)*A_cap/(Vmyo*Fara)

aux ikr=i_Kr
aux iks=i_Ks
# only t, Ca_JSR, V

# Numerical and plotting parameters for xpp
@ method=euler, bounds=1000000, dt=0.01, maxstores=1000000, total=3000
@ ylo=-95, yhi=55, yp=V, xlo=0, xhi=3000, nout=1

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 [189] name constant value unit derived unit sbo cvterm
v v = -85.2119207874627 -85.2119207874627 None
cai cai = 0.000117482029668194 0.000117482029668194 None
m m = 0.00149183115674257 0.00149183115674257 None
h h = 0.985596581239651 0.985596581239651 None
j j = 0.990898461370389 0.990898461370389 None
d d = 5.82597094505446e-6 5.82597094505446e-06 None
f f = 0.997765362821995 0.997765362821995 None
b b = 0.00136737866785149 0.00136737866785149 None
g g = 0.98881442877378 0.98881442877378 None
xr xr = 0.000204700363126417 0.000204700363126417 None
xs1 xs1 = 0.00660746743356887 0.00660746743356887 None
xs2 xs2 = 0.0303768241233812 0.0303768241233812 None
zdv zdv = 0.0144622472219576 0.0144622472219576 None
ydv ydv = 0.999945568566232 0.999945568566232 None
ca_jsr ca_jsr = 1.12791401197882 1.12791401197882 None
ca_nsr ca_nsr = 1.76731003671612 1.76731003671612 None
aptrack aptrack = 9.65910542308504e-196 9.65910542308504e-196 None
aptrack2 aptrack2 = 5.33944967562997e-195 5.33944967562997e-195 None
aptrack3 aptrack3 = 0.000129515197402902 0.000129515197402902 None
cafluxtr cafluxtr = 2.69380318286645e-196 2.69380318286645e-196 None
o_track o_track = 0 0.0 None
o_track2 o_track2 = 0 0.0 None
o_track3 o_track3 = 0 0.0 None
na_i na_i = 13.3649235394859 13.3649235394859 None
ki ki = 141.056872392446 141.056872392446 None
period period = 200 200.0 None
pulse pulse = 25.5 25.5 None
tf tf = 0 0.0 None
tp tp = 2 2.0 None
tstart tstart = 50 50.0 None
rgas rgas = 8314 8314.0 None
temp temp = 310 310.0 None
fara fara = 96485 96485.0 None
delta_m delta_m = 1e-5 1e-05 None
gg_nai gg_nai = 0.75 0.75 None
gg_nao gg_nao = 0.75 0.75 None
gg_ki gg_ki = 0.75 0.75 None
gg_ko gg_ko = 0.75 0.75 None
gg_cai gg_cai = 1 1.0 None
gg_cao gg_cao = 0.341 0.341 None
km_ca km_ca = 0.0006 0.0006 None
g_cat g_cat = 0.05 0.05 None
k_mpca k_mpca = 0.0005 0.0005 None
i_pca i_pca = 1.15 1.15 None
ii_nak ii_nak = 2.25 2.25 None
k_mnai k_mnai = 10 10.0 None
k_mko k_mko = 1.5 1.5 None
k_m_ns_ca k_m_ns_ca = 0.0012 0.0012 None
a_cap a_cap = 0.0001534 0.0001534 None
g_rel_max g_rel_max = 60 60.0 None
grelover grelover = 4 4.0 None
tau_tr tau_tr = 180 180.0 None
k_mrel k_mrel = 0.0008 0.0008 None
delcaith delcaith = 0.00018 0.00018 None
csqn_max csqn_max = 10 10.0 None
kmcsqn kmcsqn = 0.8 0.8 None
k_mup k_mup = 0.00092 0.00092 None
iup iup = 0.00875 0.00875 None
cansrmax cansrmax = 15 15.0 None
k_mtn k_mtn = 0.0005 0.0005 None
kmcmdn kmcmdn = 0.00238 0.00238 None
tn_max tn_max = 0.07 0.07 None
cmdnmax cmdnmax = 0.05 0.05 None
csqn_th csqn_th = 0.7 0.7 None
log_th log_th = 0.98 0.98 None
taut1 taut1 = 0.5 0.5 None
taut2 taut2 = 0.5 0.5 None
gg_kr gg_kr = 0.02614 0.02614 None
gg_ks gg_ks = 0.433 0.433 None
nao nao = 140 140.0 None
ko ko = 5.4 5.4 None
cao cao = 1.8 1.8 None
g_na g_na = 16 16.0 None
p_ca p_ca = 0.00054 0.00054 None
p_na p_na = 6.75e-7 6.75e-07 None
p_k p_k = 1.93e-7 1.93e-07 None
pnak pnak = 0.01833 0.01833 None
gg_k1 gg_k1 = 0.75 0.75 None
g_kp g_kp = 0.00552 0.00552 None
g_k_na g_k_na = 0.12848 0.12848 None
nkna nkna = 2.8 2.8 None
kdkna kdkna = 66 66.0 None
ikatp_on ikatp_on = 1 1.0 None
natp natp = 0.24 0.24 None
n_area n_area = 5e-5 5e-05 None
atpi atpi = 3 3.0 None
hatp hatp = 2 2.0 None
katp katp = 0.00025 0.00025 None
g_nab g_nab = 0.004 0.004 None
g_cab g_cab = 0.003016 0.003016 None
p_ns_ca p_ns_ca = 1.75e-7 1.75e-07 None
c1 c1 = 0.00025 0.00025 None
c2 c2 = 0.0001 0.0001 None
gamma gamma = 0.15 0.15 None
ts 0.0 dimensionless None
iapp 0.0 dimensionless None
e_na 0.0 dimensionless None
alpha_m 0.0 dimensionless None
beta_m 0.0 dimensionless None
alpha_h 0.0 dimensionless None
beta_h 0.0 dimensionless None
alpha_j 0.0 dimensionless None
beta_j 0.0 dimensionless None
d_inf 0.0 dimensionless None
tau_d 0.0 dimensionless None
alpha_d 0.0 dimensionless None
beta_d 0.0 dimensionless None
f_inf 0.0 dimensionless None
tau_f 0.0 dimensionless None
alpha_f 0.0 dimensionless None
beta_f 0.0 dimensionless None
f_ca 0.0 dimensionless None
b_inf 0.0 dimensionless None
tau_b 0.0 dimensionless None
g_inf 0.0 dimensionless None
tau_g 0.0 dimensionless None
g_kr 0.0 dimensionless None
rect 0.0 dimensionless None
xr_inf 0.0 dimensionless None
tau_xr 0.0 dimensionless None
e_ks 0.0 dimensionless None
g_ks 0.0 dimensionless None
xs1_inf 0.0 dimensionless None
tau_xs1 0.0 dimensionless None
xs2_inf 0.0 dimensionless None
tau_xs2 0.0 dimensionless None
g_k1 0.0 dimensionless None
e_k 0.0 dimensionless None
alpha_k1 0.0 dimensionless None
beta_k1 0.0 dimensionless None
k1_inf 0.0 dimensionless None
kp 0.0 dimensionless None
i_kp 0.0 dimensionless None
pona 0.0 dimensionless None
pov 0.0 dimensionless None
g_k_atp 0.0 dimensionless None
patp 0.0 dimensionless None
gkbart 0.0 dimensionless None
g_to 0.0 dimensionless None
rvdv 0.0 dimensionless None
alpha_zdv 0.0 dimensionless None
beta_zdv 0.0 dimensionless None
tau_zdv 0.0 dimensionless None
zdv_ss 0.0 dimensionless None
alpha_ydv 0.0 dimensionless None
beta_ydv 0.0 dimensionless None
tau_ydv 0.0 dimensionless None
ydv_ss 0.0 dimensionless None
e_ca 0.0 dimensionless None
sigma 0.0 dimensionless None
f_nak 0.0 dimensionless None
ins_na 0.0 dimensionless None
ins_k 0.0 dimensionless None
vmyo 0.0 dimensionless None
vjsr 0.0 dimensionless None
vnsr 0.0 dimensionless None
kleak 0.0 dimensionless None
i_na 0.0 dimensionless None
icaca 0.0 dimensionless None
icana 0.0 dimensionless None
icak 0.0 dimensionless None
i_caca 0.0 dimensionless None
i_cana 0.0 dimensionless None
i_cak 0.0 dimensionless None
ical 0.0 dimensionless None
icat 0.0 dimensionless None
i_kr 0.0 dimensionless None
i_ks 0.0 dimensionless None
i_k1 0.0 dimensionless None
ikna 0.0 dimensionless None
ikatp 0.0 dimensionless None
i_to 0.0 dimensionless None
ipca 0.0 dimensionless None
i_ca_b 0.0 dimensionless None
inab 0.0 dimensionless None
inak 0.0 dimensionless None
i_naca 0.0 dimensionless None
insna 0.0 dimensionless None
insk 0.0 dimensionless None
insca 0.0 dimensionless None
i_rel 0.0 dimensionless None
i_up 0.0 dimensionless None
i_leak 0.0 dimensionless None
i_tr 0.0 dimensionless None
g_relvis 0.0 dimensionless None
dv_dt 0.0 dimensionless None
ikr 0.0 dimensionless None
iks 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [120]   assignment name derived units sbo cvterm
d f/dt = alpha_f 1 f beta_f f None
d b/dt = b_inf b tau_b None
d aptrack/dt = heav dv_dt 150 100 1 aptrack taut1 aptrack heav 150 dv_dt taut2 aptrack None
d aptrack2/dt = heav 0.2 aptrack heav aptrack 0.18 100 1 aptrack2 0.5 aptrack2 heav aptrack 0.2 0.5 aptrack2 heav 0.18 aptrack 0.5 aptrack2 None
d aptrack3/dt = heav 0.2 aptrack heav aptrack 0.18 100 1 aptrack3 0.5 aptrack3 heav aptrack 0.2 0.01 aptrack3 heav 0.18 aptrack 0.01 aptrack3 None
d cafluxtr/dt = heav aptrack 0.2 1 a_cap i_caca icat 2 i_naca ipca i_ca_b 2 vmyo fara heav aptrack2 0.01 heav 0.2 aptrack 0 heav 0.01 aptrack2 0.5 cafluxtr None
d o_track/dt = heav 1 1 kmcsqn ca_jsr csqn_th heav 0.37 o_track3 heav 0.37 aptrack3 50 1 o_track heav csqn_th 1 1 kmcsqn ca_jsr 0.5 o_track heav o_track3 0.37 0.5 o_track heav aptrack3 0.37 0.5 o_track None
d o_track2/dt = heav o_track log_th heav log_th o_track2 50 1 o_track2 heav log_th o_track 0.5 o_track2 heav o_track2 log_th 0.5 o_track2 None
d o_track3/dt = heav o_track log_th heav log_th o_track3 50 1 o_track3 heav log_th o_track 0.01 o_track3 heav o_track3 log_th 0.01 o_track3 None
d v/dt = dv_dt None
d m/dt = alpha_m 1 m beta_m m None
d h/dt = alpha_h 1 h beta_h h None
d j/dt = alpha_j 1 j beta_j j None
d d/dt = alpha_d 1 d beta_d d None
d g/dt = g_inf g tau_g None
d xr/dt = xr_inf xr tau_xr None
d xs1/dt = xs1_inf xs1 tau_xs1 None
d xs2/dt = xs2_inf xs2 tau_xs2 None
d zdv/dt = zdv_ss zdv tau_zdv None
d ydv/dt = ydv_ss ydv tau_ydv None
d cai/dt = 1 1 cmdnmax kmcmdn kmcmdn cai 2 tn_max k_mtn k_mtn cai 2 1 a_cap i_caca icat 2 i_naca ipca i_ca_b 2 vmyo fara i_rel vjsr vmyo i_leak i_up vnsr vmyo None
d ca_jsr/dt = 1 1 csqn_max kmcsqn kmcsqn ca_jsr 2 i_tr i_rel None
d ca_nsr/dt = i_tr vjsr vnsr i_leak i_up None
d na_i/dt = 1 i_na i_cana inab insna i_naca 3 inak 3 a_cap vmyo fara None
d ki/dt = 1 iapp i_cak i_kr i_ks i_k1 i_kp ikna ikatp i_to insk inak 2 a_cap vmyo fara None
ts = t tstart None
iapp = pulse heav mod ts period tf heav mod ts period tf tp None
e_na = rgas temp fara nao na_i None
alpha_m = heav v 47.13 delta_m 0.32 v 47.13 1 0.1 v 47.13 heav delta_m v 47.13 3.2 None
beta_m = 0.08 v 11 None
alpha_h = heav 40 v 0.135 80 v 6.8 heav v 40 0 None
beta_h = heav 40 v 3.56 0.079 v 310000 0.35 v heav v 40 1 0.13 1 v 10.66 11.1 None
alpha_j = heav 40 v 127140 0.2444 v 3.474 -5 0.04391 v v 37.78 1 0.311 v 79.23 heav v 40 0 None
beta_j = heav v 40 0.1212 0.01052 v 1 0.1378 40.14 v heav v 40 0.3 2.535 -7 v 1 0.1 32 v None
d_inf = 1 1 v 10 6.24 None
tau_d = heav 1 -5 v 10 1 0.035 6.24 2 heav v 10 1 -5 1 d_inf 1 v 10 6.24 0.035 v 10 None
alpha_d = d_inf tau_d None
beta_d = 1 d_inf tau_d None
f_inf = 1 1 v 35.06 8.6 0.6 1 50 v 20 None
tau_f = 1 0.0197 0.0337 v 10 2 0.02 None
alpha_f = f_inf tau_f None
beta_f = 1 f_inf tau_f None
f_ca = 1 1 cai km_ca None
b_inf = 1 1 v 14 10.8 None
tau_b = 3.7 6.1 1 v 25 4.5 None
g_inf = 1 1 v 60 5.6 None
tau_g = heav 0 v 0.875 v 12 heav v 0 12 None
g_kr = gg_kr ko 5.4 0.5 None
rect = 1 1 v 9 22.4 None
xr_inf = 1 1 v 21.5 7.5 None
tau_xr = 1 0.00138 v 14.2 1 0.123 v 14.2 0.00061 v 38.9 0.145 v 38.9 1 None
e_ks = rgas temp fara ko pnak nao ki pnak na_i None
g_ks = gg_ks 1 0.6 1 3.8 -5 cai 1.4 None
xs1_inf = 1 1 v 1.5 16.7 None
tau_xs1 = 1 7.19 -5 v 30 1 0.148 v 30 0.000131 v 30 0.0687 v 30 1 None
xs2_inf = 1 1 v 1.5 16.7 None
tau_xs2 = 4 7.19 -5 v 30 1 0.148 v 30 0.000131 v 30 0.0687 v 30 1 None
g_k1 = gg_k1 ko 5.4 0.5 None
e_k = rgas temp fara ko ki None
alpha_k1 = 1.02 1 0.2385 v e_k 59.215 None
beta_k1 = 1 0.49124 0.08032 v e_k 5.476 0.06175 v e_k 594.31 1 0.5143 v e_k 4.753 None
k1_inf = alpha_k1 alpha_k1 beta_k1 None
kp = 1 1 7.488 v 5.98 None
i_kp = g_kp kp v e_k None
pona = 0.85 1 kdkna na_i nkna None
pov = 0.8 0.65 1 v 125 15 None
g_k_atp = ikatp_on 0.000193 n_area None
patp = 1 1 atpi katp hatp None
gkbart = g_k_atp patp ko 4 natp None
g_to = 0 0.5 None
rvdv = v 100 None
alpha_zdv = 10 v 40 25 1 v 40 25 None
beta_zdv = 10 v 90 25 1 v 90 25 None
tau_zdv = 1 alpha_zdv beta_zdv None
zdv_ss = alpha_zdv alpha_zdv beta_zdv None
alpha_ydv = 0.015 1 v 60 5 None
beta_ydv = 0.1 v 25 5 1 v 25 5 None
tau_ydv = 1 alpha_ydv beta_ydv None
ydv_ss = alpha_ydv alpha_ydv beta_ydv None
e_ca = rgas temp 2 fara cao cai None
sigma = 1 7 nao 67.3 1 None
f_nak = 1 1 0.1245 0.1 v fara rgas temp 0.0365 sigma v fara rgas temp None
ins_na = p_ns_ca 1 2 v fara 2 rgas temp gg_nai na_i 1 v fara rgas temp gg_nao nao 1 v fara rgas temp 1 None
ins_k = p_ns_ca 1 2 v fara 2 rgas temp gg_ki ki 1 v fara rgas temp gg_ko ko 1 v fara rgas temp 1 None
vmyo = 0.68 0.1 0.011 2 None
vjsr = 0.0048 0.1 0.011 2 None
vnsr = 0.0552 0.1 0.011 2 None
kleak = iup cansrmax None
i_na = g_na m 3 h j v e_na None
icaca = p_ca 2 2 v fara 2 rgas temp gg_cai cai 2 v fara rgas temp gg_cao cao 2 v fara rgas temp 1 None
icana = p_na 1 2 v fara 2 rgas temp gg_nai na_i 1 v fara rgas temp gg_nao nao 1 v fara rgas temp 1 None
icak = p_k 1 2 v fara 2 rgas temp gg_ki ki 1 v fara rgas temp gg_ko ko 1 v fara rgas temp 1 None
i_caca = d f f_ca icaca None
i_cana = d f f_ca icana None
i_cak = d f f_ca icak None
ical = i_caca i_cak i_cana None
icat = g_cat b b g v e_ca None
i_kr = g_kr xr rect v e_k None
i_ks = g_ks xs1 xs2 v e_ks None
i_k1 = g_k1 k1_inf v e_k None
ikna = g_k_na pona pov v e_k None
ikatp = gkbart v e_k None
i_to = g_to zdv 3 ydv rvdv v e_k None
ipca = i_pca cai k_mpca cai None
i_ca_b = g_cab v e_ca None
inab = g_nab v e_na None
inak = ii_nak f_nak 1 1 k_mnai na_i 2 ko ko k_mko None
i_naca = c1 gamma 1 v fara rgas temp v fara rgas temp na_i 3 cao nao 3 cai 1 c2 gamma 1 v fara rgas temp v fara rgas temp na_i 3 cao nao 3 cai None
insna = ins_na 1 1 k_m_ns_ca cai 3 None
insk = ins_k 1 1 k_m_ns_ca cai 3 None
insca = insna insk None
i_rel = g_relvis ca_jsr cai None
i_up = iup cai cai k_mup None
i_leak = kleak ca_nsr None
i_tr = ca_nsr ca_jsr tau_tr None
g_relvis = heav cafluxtr delcaith g_rel_max cafluxtr delcaith k_mrel cafluxtr delcaith 1 aptrack2 aptrack2 heav delcaith cafluxtr heav o_track2 0 grelover 1 o_track2 o_track2 heav delcaith cafluxtr 0 heav cafluxtr delcaith 0 heav 0 o_track2 0 None
dv_dt = i_na ical icat i_kr i_ks ikna ikatp i_to i_k1 i_kp i_naca ipca inab i_ca_b inak insca iapp None
ikr = i_kr None
iks = i_ks None
t = time None