Model Mathematics

Component: environment

Component: cell

FonRT=FRT past=timestim_periodstim_period i_Stim=stim_amplitudeiftime-paststim_offsettime-paststim_offset+stim_duration0otherwise I_tot=i_Na+i_Ca+i_to+i_Kr+i_Ks+i_K1+i_NaCa+i_NaK+i_b_Na+i_b_Ca+i_Stim ddtimeV=dVdt dVdt=V_hold-V1ifV_clamp=1time500time>800V_step-V1ifV_clamp=1time>500time800-1I_tototherwise

Component: INa

i_Na=g_Nam3hjV-E_Na E_Na=lnNaoNaiFonRT

Component: INa_m_gate

alpha_m=0.32V+47.131--0.1V+47.13if|V+47.13|>0.0013.2otherwise beta_m=0.08-V11 ddtimem=alpha_m1-m-beta_mm

Component: INa_h_gate

alpha_h=0.13580+V-6.8ifV<-400otherwise beta_h=3.560.079V+3.1e50.35VifV<-4010.131+-V+10.6611.1otherwise ddtimeh=alpha_h1-h-beta_hh

Component: INa_j_gate

alpha_j=-1271400.244V-3.474e-5-0.04391VV+37.781+0.311V+79.23ifV<-400otherwise beta_j=0.1212-0.01052V1+-0.1378V+40.14ifV<-400.3-2.535e-7V1+-0.1V+32otherwise ddtimej=alpha_j1-j-beta_jj

Component: ICa

i_Ca=g_Ca_maxdff_CaV-E_Ca E_Ca=lnCaoCai2FonRT

Component: ICa_d_gate

alpha_d=14.985916.68132π-V-22.3616.681322 beta_d=0.1471-5.314.932π-V-6.274414.9322 ddtimed=alpha_d1-d-beta_dd

Component: ICa_f_gate

alpha_f=6.872e-31+V-6.15466.1223 beta_f=0.0687-0.1081V+9.8255+0.01121+-0.2779V+9.8255+5.474e-4 ddtimef=alpha_f1-f-beta_ff

Component: ICa_f_Ca_gate

f_Ca=Km_CaKm_Ca+Cai

Component: Ito

g_to_max=0.3iffailing=00.191otherwise i_to=g_to_maxrtV-E_to E_to=ln0.043Nao+Ko0.043Nai+KiFonRT

Component: Ito_r_gate

ddtimer=alpha_r1-r-beta_rr alpha_r=0.5266-0.0166V-42.29121+-0.0943V-42.2912 beta_r=5.186e-5V+0.5149-0.1344V-5.00271+-0.1348V-5.186e-5

Component: Ito_t_gate

ddtimet=alpha_t1-t-beta_tt alpha_t=5.612e-5V+0.0721-0.173V+34.25311+-0.1732V+34.2531 beta_t=1.215e-4V+0.0767-1.66e-9V+34.02351+-0.1604V+34.0235

Component: IKs

i_Ks=g_Ks_maxXs2V-E_Ks E_Ks=ln0.01833Nao+Ko0.01833Nai+KiFonRT

Component: IKs_Xs_gate

ddtimeXs=alpha_Xs1-Xs-beta_XsXs alpha_Xs=3.013e-31+7.4454-V+1014.3171 beta_Xs=5.87e-31+5.95+V+1015.82

Component: IKr

i_Kr=g_Kr_maxXrrikV-E_K E_K=lnKoKiFonRT rik=11+V+2623

Component: IKr_Xr_gate

ddtimeXr=alpha_Xr1-Xr-beta_XrXr alpha_Xr=0.0055.266e-4V+4.0671+-0.1262V+4.067 beta_Xr=0.0161.6e-3V+65.661+0.0783V+65.66

Component: IK1

g_K1_max=2.5iffailing=02otherwise i_K1=g_K1_maxK1_infinityV-E_K1 E_K1=lnKoKiFonRT

Component: IK1_K1_gate

K1_infinity=alpha_K1alpha_K1+beta_K1 alpha_K1=0.11+0.06V-E_K1+200 beta_K1=32e-4V+100+-E_K1+10.1V-10+E_K11+-0.5V-E_K1

Component: ICab

g_b_Ca_max=0.00085iffailing=00.0013otherwise i_b_Ca=g_b_Ca_maxV-E_Ca

Component: INab

g_b_Na_max=0.001iffailing=00otherwise i_b_Na=g_b_Na_maxV-E_Na

Component: INaK

I_NaK_max=1.3iffailing=00.75otherwise f_NaK=11+0.1245-0.1VFonRT+0.0365sigma-VFonRT sigma=17Nao67.3-1 i_NaK=I_NaK_maxf_NaK11+K_mNaiNai1.5KoKo+K_mKo

Component: INaCa

K_NaCa=1000iffailing=01650otherwise i_NaCa=K_NaCa1K_mNa3+Nao31K_mCa+Cao11+K_sateta-1VFonRTetaVFonRTNai3Cao-eta-1VFonRTNao3Cai

Component: Irel

ddtimeAPtrack=1001-APtrack-0.5APtrackifdVdt>150-0.5APtrackotherwise ddtimeAPtrack2=1001-APtrack2-0.5APtrack2ifAPtrack<0.2APtrack>0.18-0.5APtrack2otherwise ddtimeAPtrack3=1001-APtrack3-0.5APtrack3ifAPtrack<0.2APtrack>0.18-0.01APtrack3otherwise ddtimeCainfluxtrack=-Cmi_Ca-i_NaCa+i_b_Ca2V_myoFifAPtrack>0.20ifAPtrack2>0.01APtrack0.2-0.5Cainfluxtrackotherwise ddtimeOVRLDtrack=0501-OVRLDtrackif11+K_mCSQNCa_JSR>CSQNthreshOVRLDtrack3<0.37APtrack3<0.37-00.5OVRLDtrackotherwise ddtimeOVRLDtrack2=0501-OVRLDtrack2ifOVRLDtrack>LogicthreshOVRLDtrack2<Logicthresh-00.5OVRLDtrack2otherwise ddtimeOVRLDtrack3=0501-OVRLDtrack3ifOVRLDtrack>LogicthreshOVRLDtrack3<Logicthresh-00.01OVRLDtrack3otherwise G_rel=1G_rel_maxCainfluxtrack-delta_Ca_ithK_mrel+Cainfluxtrack-delta_Ca_ith1-APtrack2APtrack2ifCainfluxtrack>delta_Ca_ith0G_rel_overload1-OVRLDtrack2OVRLDtrack2ifCainfluxtrackdelta_Ca_ithOVRLDtrack2>00otherwise i_rel=G_relCa_JSR-Cai

Component: Iup

I_up_max=0.0045iffailing=00.0015otherwise i_up=I_up_maxCaiCai+K_mup

Component: Ileak

K_leak=0.00026iffailing=00.00017otherwise i_leak=K_leakCa_NSR

Component: Itr

i_tr=1Ca_NSR-Ca_JSRtau_tr

Component: calcium_buffers_in_the_myoplasm

Cai_bufc=11+buffonCMDN_maxK_mCMDNK_mCMDN+Cai2+Tn_maxK_mTnK_mTn+Cai2

Component: calcium_buffers_in_the_JSR

Ca_JSR_bufc=11+CSQN_maxK_mCSQNK_mCSQN+Ca_JSR2

Component: Ionic_concentrations

volume=πpreplengthradius2 V_myo=0.68volume V_NSR=0.0552volume V_JSR=0.0048volume ddtimeNai=-1Cmi_Na+i_b_Na+i_NaCa3+i_NaK3V_myoF ddtimeKi=-1Cmi_to+i_Kr+i_K1+i_Stim+i_Ks-2i_NaKV_myoF ddtimeCai=Cai_bufc-Cmi_Ca-2i_NaCa+i_b_Ca2V_myoF+i_relV_JSRV_myo+i_leak-i_upV_NSRV_myo ddtimeCa_JSR=Ca_JSR_bufci_tr-i_rel ddtimeCa_NSR=-1i_leak+V_JSRV_NSRi_tr-i_up

Component: transmembrane_currents

Component: intracellular_currents

Component: Ions_n_reversal_potentials