Model Mathematics

Component: environment

Component: extracellular

Component: physical_constants

Component: cell_geometry

vcell=10003.14159265radradL Ageo=2greekpiradrad+2greekpiradL Acap=2Ageo vmyo=0.6vcell vnsr=0.04vcell vjsr=0.002vcell vcsr=0.008vcell vss=0.02vcell vsl=0.15vcell

Component: membrane

Istim=ampiftimeduration0otherwise vffrt=vFFRT vfrt=vFRT ddtimev=-INa+INaL+Ito+Isus+ICaL+ICaT+ICaNa+ICaK+IKr+IKs+If+IK1+INaCa_i+INaCa_ss+INaK+INab+IpCa+ICab+Istim

Component: CaMK

CaMKb=CaMKo1-CaMKt1+KmCaMcass CaMKa=CaMKb+CaMKt ddtimeCaMKt=aCaMKCaMKbCaMKb+CaMKt-bCaMKCaMKt

Component: intracellular_ions

ddtimenai=JgapNavslvmyo ddtimenasl=-INa+INaL+3INaCa_i+3INaK+IfNa+INabAcapcmFvsl+JdiffNavssvsl+-JgapNa ddtimenass=-ICaNa+3INaCa_sscmAcapFvss-JdiffNa ddtimeki=JgapKvslvmyo ddtimeksl=-Ito+Isus+IKr+IKs+IfK+IK1+Istim-2INaKcmAcapFvsl+JdiffKvssvsl+-JgapK ddtimekss=-ICaKcmAcapFvss-JdiffK Bcai=11+cmdnmaxkmcmdnkmcmdn+cai2+trpnmaxkmtrpnkmtrpn+cai2 ddtimecai=BcaiJgapvslvmyo+Jrel2vcsrvmyo-Jup2vnsrvmyo Bcass=11+BSRmaxKmBSRKmBSR+cass2+BSLmaxKmBSLKmBSL+cass2 ddtimecass=Bcass-ICaL-2INaCa_sscmAcap2Fvss+Jrel1+Jip3vjsrvss-Jdiff Bcasl=11+cmdnmaxslkmcmdnkmcmdn+casl2+trpnmaxslkmtrpnkmtrpn+casl2 ddtimecasl=Bcasl-IpCa+ICab+ICaT-2INaCa_icmAcap2Fvsl-Jup1vnsrvsl+-Jgap+Jdiffvssvsl ddtimecansr=Jup1+Jup2-Jtr1vjsrvnsr+Jtr2vcsrvnsr Bcajsr=11+csqnmaxslkmcsqnkmcsqn+cajsr2 ddtimecajsr=BcajsrJtr1-Jrel1+Jip3 Bcacsr=11+csqnmaxkmcsqnkmcsqn+cacsr2 ddtimecacsr=BcacsrJtr2-Jrel2

Component: reversal_potentials

ENa=RTFlnnaonasl EK=RTFlnkoksl ECa=RT2Flncaocasl EKs=RTFlnko+PKNanaoksl+PKNanasl

Component: INa

mss=11+-v+mssV1mssV2 tm=1mtD1v+mtV1mtV2+mtD2-v+mtV3mtV4 ddtimem=mss-mtm hss=11+v+hssV1hssV2 thf=13.686e-6-v+3.88757.8579+16v-0.49639.1843 ths=10.009794-v+17.9528.05+0.3343v+5.7356.66 Ahs=1-Ahf ddtimehf=hss-hfthf ddtimehs=hss-hsths h=Ahfhf+Ahshs jss=hss tj=4.859+10.8628-v+116.72587.6005+1.1096v+6.27199.0358 ddtimej=jss-jtj hssp=11+v+84.76.22 thsp=3ths ddtimehsp=hssp-hspthsp hp=Ahfhf+Ahshsp tjp=1.46tj ddtimejp=jss-jptjp fINap=11+KmCaMKCaMKa INa=GNav-ENam31-fINaphj+fINaphpjp

Component: INaL

mLss=11+-v+42.855.264 tmL=tm ddtimemL=mLss-mLtmL hLss=11+v+87.617.488 ddtimehL=hLss-hLthL hLssp=11+v+93.817.488 thLp=3thL ddtimehLp=hLssp-hLpthLp fINaLp=11+KmCaMKCaMKa INaL=GNaLv-ENamL1-fINaLphL+fINaLphLp

Component: Ito

ass=11+20-v13 iss=11+27+v13 taua=1.051511.20891+-v-18.409929.3814+3.51+v+10029.3814 tauis=43+10.001416-v+96.5259.05+0.0000000178v+114.18.079 tauif=6.162+10.3933-v+100100+0.08004v-88.59 ddtimea=ass-ataua ddtimei1=iss-i1tauis ddtimei2=iss-i2tauif Ito=Gtoai1i2v-EK

Component: Isus

asus=11+-v-1216 Isus=Gsusasusv-EK

Component: ICaL

dss=11+-v+3.94+3.34.23 td=0.6+1-0.05v+6+0.09v+14 ddtimed=dss-dtd fss=11+v+19.58+3.33.696 tff=7+10.0045-v+20+15.1910+0.0045v+20+15.1910 tfs=1000+10.000035-v+5+15.194+0.000035v+5+15.196 Aff=0.6 Afs=1-Aff ddtimeff=fss-fftff ddtimefs=fss-fstfs f=Affff+Afsfs fcass=fss tfcaf=0.727+10.04-v+15.19-47+0.04v+15.19-47 tfcas=0.49100+10.00012-v+15.193+0.00012v+15.197 Afcaf=0.3+0.61+v-1010 Afcas=1-Afcaf ddtimefcaf=fcass-fcaftfcaf ddtimefcas=fcass-fcastfcas fca=Afcaffcaf+Afcasfcas tjca=75 ddtimejca=fcass-jcatjca tffp=2.5tff ddtimeffp=fss-ffptffp fp=Affffp+Afsfs tfcafp=2.5tfcaf ddtimefcafp=fcass-fcafptfcafp fcap=Afcaffcafp+Afcasfcas km2n=jca1 anca=1k2nkm2n+1+Kmncass4 ddtimenca=ancak2n-ncakm2n PhiCaL=4vffrtcass2vfrt-0.341cao2vfrt-1 PhiCaNa=1vffrt0.75nass1vfrt-0.75nao1vfrt-1 PhiCaK=1vffrt0.75kss1vfrt-0.75ko1vfrt-1 PCap=1.1PCa PCaNa=0.00125PCa PCaK=3.574e-4PCa PCaNap=0.00125PCap PCaKp=3.574e-4PCap fICaLp=11+KmCaMKCaMKa ICaL=1-fICaLpPCaPhiCaLdf1-nca+jcafcanca+fICaLpPCapPhiCaLdfp1-nca+jcafcapnca ICaNa=1-fICaLpPCaNaPhiCaNadf1-nca+jcafcanca+fICaLpPCaNapPhiCaNadfp1-nca+jcafcapnca ICaK=1-fICaLpPCaKPhiCaKdf1-nca+jcafcanca+fICaLpPCaKpPhiCaKdfp1-nca+jcafcapnca

Component: ICaT

bss=11+-v+307 gss=11+v+615 taub=11.068v+16.330+1.068-v+16.330 taug=10.015v+71.715.4+0.015-v+71.783.3 ddtimeb=bss-btaub ddtimeg=gss-gtaug ICaT=GCaTbgv-ECa

Component: IKr

xrss=11+-v+8.3376.789 txrf=12.98+10.3652v+17.6-31.663.869+4.123e-5-v+17.6-47.7820.38 txrs=1.865+10.06629v+17.2-34.77.355+1.128e-5-v+17.2-29.7425.94 Axrf=11+v+54.8138.21 Axrs=1-Axrf ddtimexrf=xrss-xrftxrf ddtimexrs=xrss-xrstxrs xr=Axrfxrf+Axrsxrs rkr=11+v+550.327511+v-100.3230 IKr=GKrko5.4xrrkrv-EK

Component: IKs

xs1ss=11+-v+11.68.932 txs1=817.3+12.326e-4v+48.2817.8+0.001292-v+210230 ddtimexs1=xs1ss-xs1txs1 xs2ss=xs1ss txs2=10.01v-5020+0.0193-v+66.5431 ddtimexs2=xs2ss-xs2txs2 KsCa=1+0.61+3.8e-5casl1.4 IKs=GKsKsCaxs1xs2v-EKs

Component: If

yss=11+v+879.5 tauy=2000v+5760+-v+13210 ddtimey=yss-ytauy IfNa=GfNayyv-ENa IfK=GfKyyv-EK If=IfNa+IfK

Component: IK1

xk1ss=11+-v+2.5538ko+144.591.5692ko+3.8115 txk1=122.2-v+127.220.36+v+236.869.33 ddtimexk1=xk1ss-xk1txk1 rk1=11+v+116-5.5ko11 IK1=GK12.3238ko5.4rk1xk1v-EK

Component: INaCa_i

hca=qcavFRT hna=qnavFRT h1_i=1+naslkna31+hna h2_i=naslhnakna3h1_i h3_i=1h1_i h4_i=1+naslkna11+naslkna2 h5_i=naslnaslh4_ikna1kna2 h6_i=1h4_i h7_i=1+naokna31+1hna h8_i=naokna3hnah7_i h9_i=1h7_i h10_i=kasymm+1+naokna11+naokna2 h11_i=naonaoh10_ikna1kna2 h12_i=1h10_i k1_i=h12_icaokcaon k2_i=kcaoff k3p_i=h9_iwca k3pp_i=h8_iwnaca k3_i=k3p_i+k3pp_i k4p_i=h3_iwcahca k4pp_i=h2_iwnaca k4_i=k4p_i+k4pp_i k5_i=kcaoff k6_i=h6_icaslkcaon k7_i=h5_ih2_iwna k8_i=h8_ih11_iwna x1_i=k2_ik4_ik7_i+k6_i+k5_ik7_ik2_i+k3_i x2_i=k1_ik7_ik4_i+k5_i+k4_ik6_ik1_i+k8_i x3_i=k1_ik3_ik7_i+k6_i+k8_ik6_ik2_i+k3_i x4_i=k2_ik8_ik4_i+k5_i+k3_ik5_ik1_i+k8_i E1_i=x1_ix1_i+x2_i+x3_i+x4_i E2_i=x2_ix1_i+x2_i+x3_i+x4_i E3_i=x3_ix1_i+x2_i+x3_i+x4_i E4_i=x4_ix1_i+x2_i+x3_i+x4_i allo_i=11+KmCaActcasl2 JncxNa_i=3E4_ik7_i-E1_ik8_i+E3_ik4pp_i-E2_ik3pp_i JncxCa_i=E2_ik2_i-E1_ik1_i INaCa_i=0.8Gncxallo_iznaJncxNa_i+zcaJncxCa_i h1_ss=1+nasskna31+hna h2_ss=nasshnakna3h1_ss h3_ss=1h1_ss h4_ss=1+nasskna11+nasskna2 h5_ss=nassnassh4_sskna1kna2 h6_ss=1h4_ss h7_ss=1+naokna31+1hna h8_ss=naokna3hnah7_ss h9_ss=1h7_ss h10_ss=kasymm+1+naokna11+naokna2 h11_ss=naonaoh10_sskna1kna2 h12_ss=1h10_ss k1_ss=h12_sscaokcaon k2_ss=kcaoff k3p_ss=h9_sswca k3pp_ss=h8_sswnaca k3_ss=k3p_ss+k3pp_ss k4p_ss=h3_sswcahca k4pp_ss=h2_sswnaca k4_ss=k4p_ss+k4pp_ss k5_ss=kcaoff k6_ss=h6_sscasskcaon k7_ss=h5_ssh2_sswna k8_ss=h8_ssh11_sswna x1_ss=k2_ssk4_ssk7_ss+k6_ss+k5_ssk7_ssk2_ss+k3_ss x2_ss=k1_ssk7_ssk4_ss+k5_ss+k4_ssk6_ssk1_ss+k8_ss x3_ss=k1_ssk3_ssk7_ss+k6_ss+k8_ssk6_ssk2_ss+k3_ss x4_ss=k2_ssk8_ssk4_ss+k5_ss+k3_ssk5_ssk1_ss+k8_ss E1_ss=x1_ssx1_ss+x2_ss+x3_ss+x4_ss E2_ss=x2_ssx1_ss+x2_ss+x3_ss+x4_ss E3_ss=x3_ssx1_ss+x2_ss+x3_ss+x4_ss E4_ss=x4_ssx1_ss+x2_ss+x3_ss+x4_ss allo_ss=11+KmCaActcass2 JncxNa_ss=3E4_ssk7_ss-E1_ssk8_ss+E3_ssk4pp_ss-E2_ssk3pp_ss JncxCa_ss=E2_ssk2_ss-E1_ssk1_ss INaCa_ss=0.2Gncxallo_ssznaJncxNa_ss+zcaJncxCa_ss

Component: INaK

Knai=Knai0deltavF3RT Knao=Knao01-deltavF3RT P=eP1+HKhp+naslKnap+kslKxkur a1=k1pnaslKnai31+naslKnai3+1+kslKki2-1 b1=k1mMgADP a2=k2p b2=k2mnaoKnao31+naoKnao3+1+koKko2-1 a3=k3pkoKko21+naoKnao3+1+koKko2-1 b3=k3mPH1+MgATPKmgatp a4=k4pMgATPKmgatp1+MgATPKmgatp b4=k4mkslKki21+naslKnai3+1+kslKki2-1 x1=a4a1a2+b2b4b3+a2b4b3+b3a1a2 x2=b2b1b4+a1a2a3+a3b1b4+a2a3b4 x3=a2a3a4+b3b2b1+b2b1a4+a3a4b1 x4=b4b3b2+a3a4a1+b2a4a1+b3b2a1 E1=x1x1+x2+x3+x4 E2=x2x1+x2+x3+x4 E3=x3x1+x2+x3+x4 E4=x4x1+x2+x3+x4 JnakNa=3E1a3-E2b3 JnakK=2E4b1-E3a1 INaK=PnakznaJnakNa+zkJnakK

Component: INab

INab=PNabvffrtnaslvfrt-naovfrt-1

Component: ICab

ICab=PCab4vffrtcasl2vfrt-0.341cao2vfrt-1

Component: IpCa

IpCa=GpCacaslKmCap+casl

Component: diff

JdiffNa=nass-naslsstau JgapNa=nasl-naigaptau JdiffK=kss-kslsstau JgapK=ksl-kigaptau Jdiff=cass-caslsstau Jgap=casl-caigaptau

Component: ryr

REL=-ICaLAcapvssF2+-Jrel1+Jip3vjsrvss+Jdiff ireltau=21+111+0.28CaMKa81+0.0123cajsr irelss=151+111+0.28CaMKa8REL1+1cajsr8ifREL>00otherwise ddtimeJrel1=irelss-Jrel1ireltau REL2=Jgapvslvmyo+-Jup2vnsrvmyo+Jrel2vcsrvmyo ireltau2=61+111+0.28CaMKa81+0.0123cacsr irelss2=911+111+0.28CaMKa8REL21+1cacsr8ifREL2>00otherwise ddtimeJrel2=irelss2-Jrel2ireltau2

Component: SERCA

dkmplb=dkmplbbarCaMKaCaMKa+KmCaMK dqupcamk=dqupcamkbarCaMKaCaMKa+KmCaMK Jup1=0.0002dqupcamk+11+kmup-dkmplbcasl-0.00105cansrnsrbar Jup2=0.0026dqupcamk+11+kmup-dkmplbcai-0.0042cansrnsrbar

Component: IP3

ddtimeu=cassk21-u-k2au POip3=tauip3IP3cass1-u1+IP3k0k0a1+cassk1k1a Jip3=10.92POip3cajsr-cass

Component: trans_flux

Jtr1=cansr-cajsr120 Jtr2=cansr-cacsr120