Model Mathematics

Component: main

ddta=1

Component: parameter

V_tot_tenT=Vc_tenT+VSR_tenT Vc=V_tot⁢Vc_tenTV_tot_tenT V_SR=V_tot⁢VSR_tenTV_tot_tenT

Component: parameter_Ikr

Component: parameter_Ik1

Component: parameter_Iks

Component: parameter_Ito

Component: parameter_Ica

Component: parameter_Ina

Component: parameter_If

Component: gating_Ikr

Xr1_3=Xr1_5⁢Xr1_1 Xr2_3=Xr2_5⁢Xr2_1 Xr1_4=11Xr1_2+1Xr1_6 Xr2_4=11Xr2_2+1Xr2_6 alpha_act=Xr1_1⁢ⅇvXr1_2 beta_act=Xr1_3⁢ⅇvXr1_4 x_inf_act=alpha_actalpha_act+beta_act tau_act=1alpha_act+beta_act+Xr1_7 alpha_inact=Xr2_1⁢ⅇvXr2_2 beta_inact=Xr2_3⁢ⅇvXr2_4 x_inf_inact=alpha_inactalpha_inact+beta_inact tau_inact=1alpha_inact+beta_inact+Xr2_7

Component: gating_Ik1

alpha_act=var_2⁢ⅇv+var_4var_3 beta_act=ⅇv+var_6var_5 x_inf_act_Ik1=alpha_actalpha_act+beta_act

Component: gating_Iks

ks_3=ks_5⁢ks_1 ks_4=11ks_2+1ks_6 alpha_act=ks_1⁢ⅇvks_2 beta_act=ks_3⁢ⅇvks_4 x_ks_inf_act=alpha_actalpha_act+beta_act tau_ks_act=1alpha_act+beta_act+tau_ks_const

Component: gating_Ito

r_3=r_5⁢r_1 s_3=s_5⁢s_1 r_4=11r_2+1r_6 s_4=11s_2+1s_6 alpha_act=r_1⁢ⅇvr_2 beta_act=r_3⁢ⅇvr_4 x_to_inf_act=alpha_actalpha_act+beta_act tau_to_act=1alpha_act+beta_act+tau_r_const alpha_inact=s_1⁢ⅇvs_2 beta_inact=s_3⁢ⅇvs_4 x_to_inf_inact=alpha_inactalpha_inact+beta_inact tau_to_inact=1alpha_inact+beta_inact+tau_s_const

Component: gating_Ica

d_3=d_5⁢d_1 f_3=f_5⁢f_1 d_4=11d_2+1d_6 f_4=11f_2+1f_6 alpha_act=d_1⁢ⅇvd_2 beta_act=d_3⁢ⅇvd_4 x_ca_inf_act=alpha_actalpha_act+beta_act tau_ca_act=1alpha_act+beta_act+tau_d_const alpha_inact=f_1⁢ⅇvf_2 beta_inact=f_3⁢ⅇvf_4 x_ca_inf_inact=alpha_inactalpha_inact+beta_inact tau_ca_inact=1alpha_inact+beta_inact+tau_f_const scale_Ical_Fca_Cadep=1.2 alpha_fCa=1.01.0+scale_Ical_Fca_Cadep⁢ca_i.0003258 beta_fCa=0.11.0+ⅇscale_Ical_Fca_Cadep⁢ca_i-.00050.0001 gamma_fCa=0.21.0+ⅇscale_Ical_Fca_Cadep⁢ca_i-0.000750.0008 x_fca_inf_inact=alpha_fCa+beta_fCa+gamma_fCa+0.231.46

Component: gating_Icat

x_cat_inf_act=11+ⅇ-v+26.36 tau_cat_act=11.068⁢ⅇv+26.330+1.068⁢ⅇ-v+26.330 x_cat_inf_inact=11+ⅇv+61.75.6 tau_cat_inact=10.0153⁢ⅇ-v+61.783.3+0.015⁢ⅇv+61.715.38

Component: gating_Ina

m_3=m_5⁢m_1 h_3=h_5⁢h_1 j_3=j_5⁢j_1 m_4=11m_2+1m_6 h_4=11h_2+1h_6 j_4=11j_2+1j_6 j_5=h_5 j_6=h_6 alpha_m_act=m_1⁢ⅇvm_2 beta_m_act=m_3⁢ⅇvm_4 x_na_m_inf_act=alpha_m_actalpha_m_act+beta_m_act tau_na_m_act=1alpha_m_act+beta_m_act+tau_m_const alpha_h_inact=h_1⁢ⅇvh_2 beta_h_inact=h_3⁢ⅇvh_4 x_na_h_inf_inact=alpha_h_inactalpha_h_inact+beta_h_inact tau_na_h_inact=1alpha_h_inact+beta_h_inact+tau_h_const alpha_j_inact=j_1⁢ⅇvj_2 beta_j_inact=j_3⁢ⅇvj_4 x_na_j_inf_inact=alpha_j_inactalpha_j_inact+beta_j_inact tau_na_j_inact=1alpha_j_inact+beta_j_inact+tau_j_const

Component: gating_If

xf_3=xf_5⁢xf_1 xf_4=11xf_2+1xf_6 alpha_act=xf_1⁢ⅇvxf_2 beta_act=xf_3⁢ⅇvxf_4 x_f_inf_act=alpha_actalpha_act+beta_act tau_f_act=1alpha_act+beta_act+xf_const

Component: Protocol

Ki=150ifprotocol=1148ifprotocol=2k_newifprotocol=0125otherwise

Component: Nernst

E_K=R⁢TF⁢ln⁡KoKi E_Ca=0.5⁢R⁢TF⁢ln⁡Caoca_i E_Na=R⁢TF⁢ln⁡NaoNa_i

Component: Current_Ik1

g_k1=var_1 i_k1=g_k1⁢Ko5.4⁢x_inf_act_Ik1⁢v-E_K

Component: Current_Ikr

dX_kr_act=x_inf_act-X_kr_acttau_act dX_kr_inact=x_inf_inact-X_kr_inacttau_inact g_Kr=Xr1_0 i_Kr=g_Kr⁢Ko5.4⁢X_kr_act⁢X_kr_inact⁢v-E_K

Component: Current_Iks

dX_ks_act=x_ks_inf_act-X_ks_acttau_ks_act g_Ks=ks_0 i_Ks=g_Ks⁢Ko5.4⁢X_ks_act2⁢v-E_K

Component: Current_Ito

dX_to_act=x_to_inf_act-X_to_acttau_to_act dX_to_inact=x_to_inf_inact-X_to_inacttau_to_inact g_to=r_0 i_to=g_to⁢Ko5.4⁢X_to_act⁢X_to_inact⁢v-E_K

Component: Current_IcaL

k_fca=0ifx_fca_inf_inact>X_fca_inact∧v>-601otherwise dX_ca_act=x_ca_inf_act-X_ca_acttau_ca_act dX_ca_inact=x_ca_inf_inact-X_ca_inacttau_ca_inact dX_fca_inact=k_fca⁢x_fca_inf_inact-X_fca_inacttau_fca p_CaL=d_0 p_CaL_shannonTot=p_CaL_shannonCa+p_CaL_shannonNa+p_CaL_shannonK p_CaL_shannonCap=p_CaL_shannonCap_CaL_shannonTot p_CaL_shannonNap=p_CaL_shannonNap_CaL_shannonTot p_CaL_shannonKp=p_CaL_shannonKp_CaL_shannonTot p_CaL_Ca=p_CaL_shannonCap⁢p_CaL p_CaL_Na=p_CaL_shannonNap⁢p_CaL p_CaL_K=p_CaL_shannonKp⁢p_CaL ibarca=p_CaL_Ca⁢4.0⁢v⁢F2R⁢T⁢0.341⁢ca_i⁢ⅇ2.0⁢v⁢FR⁢T-0.341⁢Caoⅇ2.0⁢v⁢FR⁢T-1.0 i_CaL_Ca=ibarca⁢X_ca_act⁢X_ca_inact⁢X_fca_inact ibarna=p_CaL_Na⁢v⁢F2R⁢T⁢0.75⁢Na_i⁢ⅇv⁢FR⁢T-0.75⁢Naoⅇv⁢FR⁢T-1.0 i_CaL_Na=ibarna⁢X_ca_act⁢X_ca_inact⁢X_fca_inact ibark=p_CaL_K⁢v⁢F2R⁢T⁢0.75⁢Ki⁢ⅇv⁢FR⁢T-0.75⁢Koⅇv⁢FR⁢T-1.0 i_CaL_K=ibark⁢X_ca_act⁢X_ca_inact⁢X_fca_inact i_CaL=i_CaL_Ca+i_CaL_Na+i_CaL_K

Component: Current_Icat

dX_cat_act=x_cat_inf_act-X_cat_acttau_cat_act dX_cat_inact=x_cat_inf_inact-X_cat_inacttau_cat_inact i_cat=g_cat⁢X_cat_act⁢X_cat_inact⁢v-E_Ca

Component: Current_Ina

dX_na_m_act=x_na_m_inf_act-X_na_m_acttau_na_m_act dX_na_h_inact=x_na_h_inf_inact-X_na_h_inacttau_na_h_inact dX_na_j_inact=x_na_j_inf_inact-X_na_j_inacttau_na_j_inact i_na=g_Na⁢X_na_m_act3⁢X_na_h_inact⁢X_na_j_inact⁢v-E_Na

Component: Current_If

dX_f_act=x_f_inf_act-X_f_acttau_f_act Na_frac=NatoK_ratioNatoK_ratio+1 i_fNa=Na_frac⁢g_f⁢X_f_act⁢v-E_Na i_fK=1-Na_frac⁢g_f⁢X_f_act⁢v-E_K i_f=i_fNa+i_fK

Component: Current_Inaca

i_naca=kNaCa⁢ⅇgamma⁢v⁢FR⁢T⁢Na_i3.0⁢Cao-ⅇgamma-1.0⁢v⁢FR⁢T⁢Nao3⁢ca_i⁢alphaKmNai3.0+Nao3⁢KmCa+Cao⁢1.0+Ksat⁢ⅇgamma-1.0⁢v⁢FR⁢T

Component: Current_Inak

i_nak=PNaK⁢Ko⁢Na_iKo+Km_K⁢Na_i+Km_Na⁢1.0+0.1245⁢ⅇ-0.1⁢v⁢FR⁢T+0.0353⁢ⅇ-v⁢FR⁢T

Component: Current_Iup

i_up=VmaxUp1.0+Kup2ca_i2

Component: Current_leak

j_leak=ca_SR-ca_i⁢V_leak

Component: Current_release

kCaSR=MaxSR-MaxSR-MinSR1+ec50SRca_SR2.5 koSRCa=koCakCaSR kiSRCa=kiCa⁢kCaSR CI=1-C-O-I dC=kim⁢CI-kiSRCa⁢ca_i⁢C-koSRCa⁢ca_i2⁢C-kom⁢O dO=koSRCa⁢ca_i2⁢C-kom⁢O-kiSRCa⁢ca_i⁢O-kim⁢I dI=kiSRCa⁢ca_i⁢O-kim⁢I-kom⁢I-koSRCa⁢ca_i2⁢CI j_rel=ks⁢O⁢ca_SR-ca_i⁢V_SRVc

Component: Current_bNa

i_b_Na=g_b_Na⁢v-E_Na

Component: Current_bCa

i_b_Ca=g_b_Ca⁢v-E_Ca

Component: Current_pCa

i_PCa=g_PCa⁢ca_ica_i+KPCa

Component: ca_SR

Ca_SR_bufSR=11+Buf_SR⁢Kbuf_SRca_SR+Kbuf_SR2 dca_SR=Ca_SR_bufSR⁢VcV_SR⁢i_up-j_rel+j_leak

Component: ca_buf

Cai_buf=11.0+Buf_C⁢Kbuf_Cca_i+Kbuf_C2