Model Mathematics

Component: main

k_fca=0ifx_fca_inf_inact>X_fca_inact∧v>-601otherwise v_clamp=v E_K=R⁢TF⁢ln⁡KoKi g_k1=var_1 i_k1=g_k1⁢Ko5.4⁢x_inf_act_Ik1⁢v-E_K ddtX_kr_act=x_inf_act-X_kr_acttau_kr_act ddtX_kr_inact=x_inf_inact-X_kr_inacttau_kr_inact g_Kr=Xr1_0 i_Kr=g_Kr⁢Ko5.4⁢X_kr_act⁢X_kr_inact⁢v-E_K ddtX_ks_act=x_ks_inf_act-X_ks_acttau_ks_act g_Ks=ks_0 i_Ks=g_Ks⁢X_ks_act2⁢v-E_K ddtX_to_act=x_to_inf_act-X_to_acttau_to_act ddtX_to_inact=x_to_inf_inact-X_to_inacttau_to_inact g_to=r_0 i_to=g_to⁢X_to_act⁢X_to_inact⁢v-E_K ddtX_ca_act=x_ca_inf_act-X_ca_acttau_ca_act ddtX_ca_inact=x_ca_inf_inact-X_ca_inacttau_ca_inact ddtX_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⁢ⅇ2.0⁢v⁢FR⁢T-0.75⁢Naoⅇ2.0⁢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⁢ⅇ2.0⁢v⁢FR⁢T-0.75⁢Koⅇ2.0⁢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 ddtX_cat_act=x_cat_inf_act-X_cat_acttau_cat_act ddtX_cat_inact=x_cat_inf_inact-X_cat_inacttau_cat_inact i_cat=g_cat⁢X_cat_act⁢X_cat_inact⁢v-E_Ca E_Ca=0.5⁢R⁢TF⁢ln⁡Caoca_i E_Na=R⁢TF⁢ln⁡NaoNa_i ddtX_na_h_inact=x_na_h_inf_inact-X_na_h_inacttau_na_h_inact ddtX_na_m_act=x_na_m_inf_act-X_na_m_acttau_na_m_act ddtX_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 ddtX_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 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 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 i_up=VmaxUp1.0+Kup2ca_i2 j_leak=ca_SR-ca_i⁢V_leak kCaSR=MaxSR-MaxSR-MinSR1+ec50SRca_SR2.5 koSRCa=koCakCaSR kiSRCa=kiCa⁢kCaSR CI=1-C-O-I ddtC=kim⁢CI-kiSRCa⁢ca_i⁢C-koSRCa⁢ca_i2⁢C-kom⁢O ddtO=koSRCa⁢ca_i2⁢C-kom⁢O-kiSRCa⁢ca_i⁢O-kim⁢I ddtI=kiSRCa⁢ca_i⁢O-kim⁢I-kom⁢I-koSRCa⁢ca_i2⁢CI j_rel=ks⁢O⁢ca_SR-ca_i⁢V_SRVc i_b_Na=g_b_Na⁢v-E_Na i_b_Ca=g_b_Ca⁢v-E_Ca i_PCa=g_PCa⁢ca_ica_i+KPCa Ca_SR_bufSR=11+Buf_SR⁢Kbuf_SRca_SR+Kbuf_SR2 ddtca_SR=Ca_SR_bufSR⁢VcV_SR⁢i_up-j_rel+j_leak Cai_buf=11.0+Buf_C⁢Kbuf_Cca_i+Kbuf_C2 ddtca_ligand=0 ddtca_i=Cai_buf⁢j_leak-i_up+j_rel-ca_ligand-i_CaL_Ca+i_cat+i_b_Ca+i_PCa-2⁢i_naca⁢1⁢Cm2⁢Vc⁢F ddtNa_i=-Cm⁢i_na+i_b_Na+i_fNa+3⁢i_nak+3⁢i_naca+i_CaL_NaF⁢Vc ddtKi=-Cm⁢i_k1+i_to+i_Kr+i_Ks+i_fK-2⁢i_nak+i_CaL_KF⁢Vc i_voltageclamp=v_clamp-vR_clamp ddtv=-i_k1+i_to+i_Kr+i_Ks+i_CaL+i_cat+i_nak+i_na+i_naca+i_PCa+i_f+i_b_Na+i_b_Ca-i_stim-i_voltageclamp

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_Ik1

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: parameter_Ikr

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_kr_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_kr_inact=1alpha_inact+beta_inact+Xr2_7

Component: parameter_Iks

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: parameter_Ito

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: parameter_Ica

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: parameter_Ina

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: parameter_If

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