Model Mathematics

Component: environment

vol_cyt=fr_cyt⁢vol_pmu vol_pmu=π⁢dia36⁢1 e 91⁢1⁢1 V_tau=R⁢TF

Component: Membrane

SVR_pmu=svrdia C_m=C_sp⁢A_pmu VD=VV_tau V=F⁢vol_cytC_sp⁢A_pmu⁢K_i-K_o+2⁢Ca_i-Ca_o+Na_i-Na_o+anoff anoff=-52.31⁢C_sp⁢A_pmuF⁢vol_cyt-2⁢1.5 e -4-Ca_o-6-Na_o-140-K_o A_pmu=SVR_pmu⁢vol_pmu⁢1 e -3⁢1 e -3⁢1 e -31 V_Ca=0.5⁢ln⁡Ca_oCa_i V_Na=ln⁡Na_oNa_i V_K=ln⁡K_oK_i

Component: Cytosol

ddtCa_i=J_Ca ddtNa_i=J_Na ddtK_i=J_K

Component: calcium_dynamics

J_Ca=J_ca-J_calb+4⁢J_cam J_ca=-12⁢F⁢vol_cyt⁢I_CaL+I_CaT+I_CaHVA+2⁢I_pmca-2⁢I_xm

Component: sodium_dynamics

J_Na=-1F⁢vol_cyt⁢3⁢I_nk+3⁢I_xm+I_Na+I_Nalk+I_NaHCN

Component: potassium_dynamics

I_K=I_Ksk+I_KHCN+I_Ka+I_Kdr+I_Kir J_K=-1F⁢vol_cyt⁢I_K-2⁢I_nk

Component: calcium_buffer_dynamics

J_calb=kcal_1⁢Calb⁢Ca_i-kcal_2⁢CaCalb J_cam=alpha_cam⁢Cam-beta_cam⁢CaCam kcam_cb=12000⁢Ca_i2 kcam_nb=3.7 e 6⁢Ca_i2 alpha_cam=kcam_cb⁢kcam_nb⁢1kcam_cb+kcam_nd+1kcam_cd+kcam_nd beta_cam=kcam_cd⁢kcam_nd⁢1kcam_cb+kcam_nd+1kcam_cd+kcam_nd CaCalb=Calbtot-Calb CaCam=Camtot-Cam ddtCalb=-J_calb ddtCam=-J_cam

Component: L_type_Ca_channel

I_CaL=g_cal⁢m_cal⁢h_cal⁢Ca_i⁢Ca_o⁢sinh⁡VD-V_Casinh⁡VDVD ddtm_cal=11+ⅇ-V+157-m_cal10⁢ⅇ-V+86.423.22+0.943 h_cal=0.000450.00045+Ca_i

Component: T_type_Ca_channel

I_CaT=g_cat⁢m_cat⁢h_cat⁢Ca_i⁢Ca_o⁢sinh⁡VD-V_Casinh⁡VDVD ddtm_cat=11+ⅇ-V+631.5-m_cat65⁢ⅇ-V+6862+12 ddth_cat=11+ⅇV+76.23-h_cat50⁢ⅇ-V+72102+10

Component: HVA_Ca_channel

I_CaHVA=g_cah⁢m_cah⁢h_cah⁢Ca_i⁢Ca_o⁢sinh⁡VD-V_Casinh⁡VDVD ddtm_cah=11+ⅇ-V+1010-m_cah0.1⁢ⅇ-V+62132+0.05 ddth_cah=11+ⅇV+485-h_cah0.5⁢ⅇ-V+55.6182+0.5

Component: transient_Na_channel

O_na=m_na3⁢h_na ddtm_na=A_mna⁢ⅇza_mna⁢VD⁢1-m_na-B_mna⁢ⅇ-zb_mna⁢VD⁢m_na ddth_na=A_hna⁢ⅇza_hna⁢VD⁢1-h_na-B_hna⁢ⅇ-zb_hna⁢VD⁢h_na I_Na=g_na⁢O_na⁢Na_i⁢Na_o⁢sinh⁡0.5⁢VD-V_Nasinh⁡0.5⁢VD0.5⁢VD

Component: Leak_Na_channel

I_Nalk=g_nalk⁢Na_i⁢Na_o⁢sinh⁡0.5⁢VD-V_Nasinh⁡0.5⁢VD0.5⁢VD

Component: HCN_channel

ddtO_hcn=kf_hcn⁢1-O_hcn-kr_hcn⁢O_hcn kf_hcn=kf_free⁢P_c+kf_bnd⁢1-P_c kr_hcn=kr_free⁢P_o+kr_bnd⁢1-P_o kf_free=0.0061+ⅇV+87.76.45 kr_free=0.081+ⅇ-V+51.77 kf_bnd=0.02681+ⅇV+94.213.3 kr_bnd=0.081+ⅇ-V+35.57 P_c=11+cAMP1.163 e -3 P_o=11+cAMP1.45 e -5 I_NaHCN=g_nahcn⁢O_hcn⁢Na_i⁢Na_o⁢sinh⁡0.5⁢VD-V_Nasinh⁡0.5⁢VD0.5⁢VD I_KHCN=g_khcn⁢O_hcn⁢K_i⁢K_o⁢sinh⁡0.5⁢VD-V_Ksinh⁡0.5⁢VD0.5⁢VD

Component: SK_K_channel

O_sk=Ca_i4.20.35 e -34.2+Ca_i4.2 I_Ksk=g_ksk⁢O_sk⁢K_i⁢K_o⁢sinh⁡0.5⁢VD-V_Ksinh⁡0.5⁢VD0.5⁢VD

Component: DR_K_channel

I_Kdr=g_kdr⁢O_kdr⁢V-V_K⁢V_tau O_kdr=m_kdr3 ddtm_kdr=11+ⅇ-V+2512-m_kdr181+ⅇV+398+1

Component: Atype_K_channel

I_Ka=g_ka⁢O_ka⁢V-V_K⁢V_tau O_ka=p_ka3⁢q_ka ddtp_ka=11+ⅇ-V+4324-p_ka2⁢ⅇ-V+5023.452+1.1 ddtq_ka=11+ⅇV+568-q_ka20

Component: IR_K_channel

I_Kir=g_kir⁢O_kir⁢V-V_K⁢V_tau O_kir=11+ⅇV+9012.1

Component: sodium_pump

I_nk=k_nk⁢k_1nk⁢P_E1Snk⁢y_nk-k_2nk⁢P_E2Snk⁢1-y_nk⁢1 k_1nk=0.371+0.094atp Na_eff=Na_o⁢ⅇ-0.82⁢VD ddty_nk=beta_nk⁢1-y_nk-alpha_nk⁢y_nk beta_nk=k_2nk⁢P_E2Snk+k_4nk⁢P_E2Dnk alpha_nk=k_1nk⁢P_E1Snk+k_3nk⁢P_E1Dnk P_E1Snk=11+K_nknaiNa_i⁢1+K_iK_nkki P_E1Dnk=11+K_nkkiK_i⁢1+Na_iK_nknai P_E2Snk=11+K_nknaoNa_eff⁢1+K_oK_nkko P_E2Dnk=11+K_nkkoK_o⁢1+Na_effK_nknao

Component: PMCA

I_pmca=K_pmca⁢k_1pc⁢P_E1Spc⁢y_pc-k_2pc⁢P_E2Spc⁢1-y_pc⁢1 K_pmca=k_pmca⁢10.56⁢CaCamCaCam+0.00005+1.2 ddty_pc=beta_pc⁢1-y_pc-alpha_pc⁢y_pc k_1pc=11+0.1atp P_E1Spc=11+K_pciCa_i P_E1pc=1-P_E1Spc P_E2Spc=11+K_pcoCa_o P_E2pc=1-P_E2Spc beta_pc=k_2pc⁢P_E2Spc+k_4pc⁢P_E2pc alpha_pc=k_1pc⁢P_E1Spc+k_3pc⁢P_E1pc K_pci=173.61+CaCam0.00005+6.4⁢1 e -5

Component: Membrane_NaCa

I_xm=k_xm⁢k_3xm⁢P_E1Dxm⁢y_xm-k_4xm⁢P_E2Dxm⁢1-y_xm ddty_xm=beta_xm⁢1-y_xm-alpha_xm⁢y_xm beta_xm=k_2xm⁢P_E2Sxm+k_4xm⁢P_E2Dxm alpha_xm=k_1xm⁢P_E1Sxm+k_3xm⁢P_E1Dxm P_E1Sxm=11+K_xmcaiCa_i⁢1+Na_iK_xmnai3 P_E1Dxm=11+K_xmnaiNa_i3⁢1+Ca_iK_xmcai P_E2Sxm=11+K_xmcaoCa_o⁢1+Na_oK_xmnao3 P_E2Dxm=11+K_xmnaoNa_o3⁢1+Ca_oK_xmcao k_1xm=1⁢ⅇ-delta_xmca⁢VD k_2xm=1⁢ⅇ1-delta_xmca⁢VD k_3xm=1⁢ⅇ1-delta_xmna⁢VD k_4xm=1⁢ⅇ-delta_xmna⁢VD