Model Mathematics

Component: environment

Component: membrane

ddtimeV=-1Cmi_Na+i_Ca_L+i_Ca_T+i_Kr+i_Ks+i_K1+i_Kp+i_NaCa+i_p_Ca+i_Na_b+i_Ca_b+i_NaK+i_ns_Ca+I_st

Component: fast_sodium_current

i_Na=g_Nam3hjV-E_Na E_Na=RTFlnNaoNai

Component: fast_sodium_current_m_gate

alpha_m=0.32V+47.131--0.1V+47.13 beta_m=0.08-V11 ddtimem=alpha_m1-m-beta_mm

Component: fast_sodium_current_h_gate

alpha_h=0.13580+V-6.8ifV<-400otherwise beta_h=3.560.079V+3100000.35VifV<-4010.131+V+10.66-11.1otherwise ddtimeh=alpha_h1-h-beta_hh

Component: fast_sodium_current_j_gate

alpha_j=-1271400.2444V-0.00003474-0.04391VV+37.781+0.311V+79.23ifV<-400otherwise beta_j=0.1212-0.01052V1+-0.1378V+40.14ifV<-400.3-0.0000002535V1+-0.1V+32otherwise ddtimej=alpha_j1-j-beta_jj

Component: L_type_Ca_channel

i_CaCa=dff_CaI_CaCa i_CaNa=dff_CaI_CaNa i_CaK=dff_CaI_CaK I_CaCa=P_Ca22VF2RTgamma_CaiCai2VFRT-gamma_CaoCao2VFRT-1 I_CaNa=P_Na12VF2RTgamma_NaiNai1VFRT-gamma_NaoNao1VFRT-1 I_CaK=P_K12VF2RTgamma_KiKi1VFRT-gamma_KoKo1VFRT-1 i_Ca_L=i_CaCa+i_CaK+i_CaNa

Component: L_type_Ca_channel_d_gate

alpha_d=d_infinitytau_d d_infinity=11+-V+106.24 tau_d=d_infinity1--V+106.240.035V+10 beta_d=1-d_infinitytau_d ddtimed=alpha_d1-d-beta_dd

Component: L_type_Ca_channel_f_gate

alpha_f=f_infinitytau_f f_infinity=11+V+35.068.6+0.61+50-V20 tau_f=10.0197-0.0337V+102+0.02 beta_f=1-f_infinitytau_f ddtimef=alpha_f1-f-beta_ff

Component: L_type_Ca_channel_f_Ca_gate

f_Ca=11+CaiKm_Ca2

Component: T_type_Ca_channel

i_Ca_T=g_Ca_Tb2gV-E_Ca E_Ca=RT2FlnCaoCai

Component: T_type_Ca_channel_b_gate

ddtimeb=b_infinity-btau_b b_infinity=11+-V+1410.8 tau_b=3.7+6.11+25+V4.5

Component: T_type_Ca_channel_g_gate

ddtimeg=g_infinity-gtau_g g_infinity=11+V+605.6 tau_g=12ifV>0-0.875V+12otherwise

Component: rapid_time_dependent_potassium_current

g_Kr=0.02614Ko5.4 E_Kr=RTFlnKoKi i_Kr=g_KrXrRrV-E_Kr

Component: rapid_time_dependent_potassium_current_Xr_gate

ddtimeXr=Xr_infinity-Xrtau_Xr Xr_infinity=11+-V+21.57.5 tau_Xr=10.00138V+14.21--0.123V+14.2+0.00061V+38.90.145V+38.9-1

Component: rapid_time_dependent_potassium_current_Rr_gate

Rr=11+V+922.4

Component: slow_time_dependent_potassium_current

g_Ks=0.057+0.191+-7.2+P_Ca0.6 E_Ks=RTFlnKo+P_NaKNaoKi+P_NaKNai P_Ca=-log10Cai+3 i_Ks=g_KsXs2V-E_Ks

Component: slow_time_dependent_potassium_current_Xs_gate

ddtimeXs=Xs_infinity-Xstau_Xs Xs_infinity=11+-V-1.516.7 tau_Xs=10.0000719V+301--0.148V+30+0.000131V+300.0687V+30-1

Component: time_independent_potassium_current

g_K1=0.75Ko5.4 E_K1=RTFlnKoKi i_K1=g_K1K1_infinityV-E_K1

Component: time_independent_potassium_current_K1_gate

alpha_K1=1.021+0.2385V-E_K1-59.215 beta_K1=0.491240.08032V+5.476-E_K1+0.06175V-E_K1+594.311+-0.5143V-E_K1+4.753 K1_infinity=alpha_K1alpha_K1+beta_K1

Component: plateau_potassium_current

E_Kp=E_K1 Kp=11+7.488-V5.98 i_Kp=g_KpKpV-E_Kp

Component: sarcolemmal_calcium_pump

i_p_Ca=I_pCaCaiK_mpCa+Cai

Component: sodium_background_current

E_NaN=E_Na i_Na_b=g_NabV-E_NaN

Component: calcium_background_current

E_CaN=RT2FlnCaoCai i_Ca_b=g_CabV-E_CaN

Component: sodium_potassium_pump

f_NaK=11+0.1245-0.1VFRT+0.0365sigma-VFRT sigma=17Nao67.3-1 i_NaK=I_NaKf_NaK11+K_mNaiNai3KoKo+K_mKo

Component: non_specific_calcium_activated_current

i_ns_Na=I_ns_Na11+K_m_ns_CaCai3 i_ns_K=I_ns_K11+K_m_ns_CaCai3 i_ns_Ca=i_ns_Na+i_ns_K I_ns_Na=P_ns_Ca12VF2RTgamma_NaiNai1VFRT-gamma_NaoNao1VFRT-1 I_ns_K=P_ns_Ca12VF2RTgamma_KiKi1VFRT-gamma_KoKo1VFRT-1

Component: Na_Ca_exchanger

i_NaCa=K_NaCa1K_mNa3+Nao31K_mCa+Cao11+K_sateta-1VFRTetaVFRTNai3Cao-eta-1VFRTNao3Cai

Component: calcium_buffers_in_the_myoplasm

TRPN_buff=TRPN_maxCaiCai+K_mTRPN CMDN_buff=CMDN_maxCaiCai+K_mCMDN Ca_JSR_new=b12+4c1-b12 delta_Ca_JSR=Ca_JSR_old-Ca_JSR_new delta_Cai=Cai_old-Cai b1=K_mCSQN-CSQN_max+CSQN_buff+delta_Ca_JSR+Ca_JSR_old c1=K_mCSQNCSQN_buff+delta_Ca_JSR+Ca_JSR_old Cai=23b2-3ccosarccos9bc-2b3+27d2b2-3c323-b3 b=CMDN_max+TRPN_max+K_mTRPN+K_mCMDN-Ca_total c=K_mCMDNK_mTRPN+TRPN_maxK_mCMDN+CMDN_maxK_mTRPN-Ca_totalK_mTRPN+K_mCMDN d=-K_mTRPNK_mCMDNCa_total Ca_total=TRPN_buff+CMDN_buff+delta_Cai+Cai_old

Component: calcium_fluxes_in_the_SR

i_rel=G_relCa_JSR_new-Cai G_rel=G_rel_maxdelta_Ca_i2-delta_Ca_ithK_mrel+delta_Ca_i2-delta_Ca_ith1--ttau_on-ttau_offifcalcium_overload=0G_rel_max1--ttau_on-ttau_offotherwise G_rel_max=0ifcalcium_overload=0delta_Ca_i2<delta_Ca_ith60ifcalcium_overload=0delta_Ca_i2delta_Ca_ith0ifCSQN_buff<CSQN_th4otherwise CSQN_buff=CSQN_maxCa_JSR_newCa_JSR_new+K_mCSQN i_up=I_upCaiCai+K_mup i_leak=K_leakCa_NSR K_leak=I_upCa_NSR_max i_tr=Ca_NSR-Ca_JSR_newtau_tr

Component: ionic_concentrations

ddtimeNai=-i_Na+i_CaNa+i_Na_b+i_ns_Na+i_NaCa3+i_NaK3A_capV_myoF ddtimeKi=-i_CaK+i_Kr+i_Ks+i_K1+i_Kp+i_ns_K+-i_NaK2A_capV_myoF ddtimeKo=i_CaK+i_Kr+i_Ks+i_K1+i_Kp+i_ns_K+-i_NaK2A_capV_cleftF ddtimeCa_NSR=-i_leak+i_tr-i_up ddtimeCa_foot=-i_CaCaA_cap2V_myoFR_A_V