Model Mathematics

Component: Environment

Component: Parameters

D_K=D_Ca1.960.79D_Na=D_Ca1.330.79

Component: Potassium_Conc_cyt

ddtime conc_K_cyt =-I_K_cyt_b +I_K_cyt_ir +I_K_cyt_ATP +I_K_cyt_v +I_K_cyt_Ca -2 I_NaK_cyt_alpha1 +F J_cyt_KCl +J_cyt_NKCC +J_Kdiff F vol_cyt

Component: Potassium_Conc_md

ddtime conc_K_md =-I_K_md_b +I_K_md_ir +I_K_md_ATP +I_K_md_v +I_K_md_Ca -2 I_NaK_md_alpha2 +F J_md_KCl +J_md_NKCC -J_Kdiff F vol_md

Component: Sodium_Conc_cyt

ddtime conc_Na_cyt =-I_Na_cyt_b +3 I_NaK_cyt_alpha1 +I_VONa_cyt +I_Na_cyt_SOC +F J_cyt_NaCl +J_cyt_NKCC +J_Nadiff F vol_cyt

Component: Sodium_Conc_md

ddtime conc_Na_md =-I_Na_md_b +3 I_NaK_md_alpha2 +I_VONa_md +I_Na_md_SOC +3 I_NaCa_md +F J_md_NaCl +J_md_NKCC -J_Nadiff F vol_md

Component: Chloride_Conc_cyt

ddtime conc_Cl_cyt =I_Cl_cyt_b +I_Cl_cyt_Ca +F J_cyt_NaCl +F J_cyt_KCl +2 J_cyt_NKCC +J_Cldiff F vol_cyt

Component: Chloride_Conc_md

ddtime conc_Cl_md =I_Cl_md_b +I_Cl_md_Ca +F J_md_NaCl +F J_md_KCl +2 J_md_NKCC -J_Cldiff F vol_md

Component: Calcium_Conc_cyt

ddtimeconc_Ca_cyt=-I_Ca_cyt_b+I_Ca_cyt_P+I_Ca_cyt_SOC+I_Ca_cyt_L+I_cyt_SERCa_IP3R+I_cyt_SERCa_RyR-I_cyt_RyR-I_cyt_IP3R+2FJ_Cadiff2Fvol_cyt_Ca-2ddtimeconc_CaMN__+ddtimeconc_CaM_C_+ddtimeconc_CaM_CM-ddtimeconc_BfCa_cyt

Component: Bf_Kinetics

ddtime conc_BfCa_cyt =k_Bf_on conc_Ca_cyt conc_Bf_tot -conc_BfCa_cyt -k_Bf_off conc_BfCa_cyt ddtime conc_BfCa_md =k_Bf_on conc_Ca_md conc_Bf_tot -conc_BfCa_md -k_Bf_off conc_BfCa_md

Component: Calcium_Conc_md

ddtimeconc_Ca_md=-I_Ca_md_b+I_Ca_md_P+I_Ca_md_SOC+I_Ca_md_L-2I_NaCa_md+I_md_SERCa_RyR+I_md_SERCa_IP3R-I_md_RyR-I_md_IP3R-2FJ_Cadiff2Fvol_md-2ddtimeconc_CaMN__+ddtimeconc_CaM_C_+ddtimeconc_CaM_CM-ddtimeconc_BfCa_md

Component: Calcium_Conc_SR

ddtimeconc_Ca_IP3R=I_cyt_SERCa_IP3R-I_cyt_IP3R+I_md_SERCa_IP3R-I_md_IP3RFvol_SR-ddtimeconc_CalseqCa_IP3Rddtimeconc_Ca_RyR=I_cyt_SERCa_RyR-I_cyt_RyR+I_md_SERCa_RyR-I_md_RyRFvol_SR-ddtimeconc_CalseqCa_RyR

Component: Background_Currents

I_K_cyt_b =G_K_b f_cyt V_m_cyt -E_K_cyt I_K_md_b =G_K_b f_md V_m_md -E_K_md I_Na_cyt_b =G_Na_b f_cyt V_m_cyt -E_Na_cyt I_Na_md_b =G_Na_b f_md V_m_md -E_Na_md I_Cl_cyt_b =G_Cl_b f_cyt V_m_cyt -E_Cl_cyt I_Cl_md_b =G_Cl_b f_md V_m_md -E_Cl_md I_Ca_cyt_b =G_Ca_b f_cyt V_m_cyt -E_Ca_cyt I_Ca_md_b =G_Ca_b f_md V_m_md -E_Ca_md

Component: Nernst_Potentials

E_K_cyt =RTF z_K lnconc_K_out conc_K_cyt E_K_md =RTF z_K lnconc_K_out conc_K_md E_Na_cyt =RTF z_Na lnconc_Na_out conc_Na_cyt E_Na_md =RTF z_Na lnconc_Na_out conc_Na_md E_Cl_cyt =RTF z_Cl lnconc_Cl_out conc_Cl_cyt E_Cl_md =RTF z_Cl lnconc_Cl_out conc_Cl_md E_Ca_cyt =RTF z_Ca lnconc_Ca_out conc_Ca_cyt E_Ca_md =RTF z_Ca lnconc_Ca_out conc_Ca_md RTF =R T F

Component: SERCA_Currents_IP3R

I_cyt_SERCa_IP3R=I_SERCA_maxf_cytconc_Ca_cytK_mfn_H-conc_Ca_IP3RK_mrn_H1+conc_Ca_cytK_mfn_H+conc_Ca_IP3RK_mrn_HI_md_SERCa_IP3R=I_SERCA_maxf_mdconc_Ca_mdK_mfn_H-conc_Ca_IP3RK_mrn_H1+conc_Ca_mdK_mfn_H+conc_Ca_IP3RK_mrn_H

Component: SERCA_Currents_RyR

I_cyt_SERCa_RyR=I_SERCA_maxf_cytconc_Ca_cytK_mfn_H-conc_Ca_RyRK_mrn_H1+conc_Ca_cytK_mfn_H+conc_Ca_RyRK_mrn_H I_md_SERCa_RyR=I_SERCA_maxf_mdconc_Ca_mdK_mfn_H-conc_Ca_RyRK_mrn_H1+conc_Ca_mdK_mfn_H+conc_Ca_RyRK_mrn_H

Component: RyR_Currents

I_cyt_RyR=v_RyRP_cyt_RyRconc_Ca_RyR-conc_Ca_cyt2Fvol_SRI_md_RyR=v_RyRP_md_RyRconc_Ca_RyR-conc_Ca_md2Fvol_SRP_cyt_RyR=w_cyt1+conc_Ca_cytK_b31+K_aconc_Ca_cyt4+conc_Ca_cytK_b3P_md_RyR=w_md1+conc_Ca_mdK_b31+K_aconc_Ca_md4+conc_Ca_mdK_b3ddtimew_cyt=K_c_ratew_inf_cyt-w_cytw_inf_cytddtimew_md=K_c_ratew_inf_md-w_mdw_inf_mdw_inf_cyt=1+K_aconc_Ca_cyt4+conc_Ca_cytK_b31+1K_c+K_aconc_Ca_cyt4+conc_Ca_cytK_b3w_inf_md=1+K_aconc_Ca_md4+conc_Ca_mdK_b31+1K_c+K_aconc_Ca_md4+conc_Ca_mdK_b3

Component: IP3R_Currents

ddtime conc_IP3_cyt =v_IP3 conc_IP3_ref conc_Ca_cyt +1 -alpha_4 k4 conc_Ca_cyt +k4 -I_r conc_IP3_cyt ddtime conc_IP3_md =v_IP3 conc_IP3_ref conc_Ca_md +1 -alpha_4 k4 conc_Ca_md +k4 -I_r conc_IP3_md I_cyt_IP3R =v_IP3R x110_cyt 3 conc_Ca_IP3R -conc_Ca_cyt 2 F vol_SR I_md_IP3R =v_IP3R x110_md 3 conc_Ca_IP3R -conc_Ca_md 2 F vol_SR

Component: IP3R_Binding_Sites

ddtime x000_cyt =-a4 conc_Ca_cyt x000_cyt -d4 x001_cyt -a5 conc_Ca_cyt x000_cyt -d5 x010_cyt ddtime x000_md =-a4 conc_Ca_md x000_md -d4 x001_md -a5 conc_Ca_md x000_md -d5 x010_md ddtime x001_cyt =a4 conc_Ca_cyt x000_cyt -d4 x001_cyt -a5 conc_Ca_cyt x001_cyt -d5 x011_cyt ddtime x001_md =a4 conc_Ca_md x000_md -d4 x001_md -a5 conc_Ca_md x001_md -d5 x011_md ddtime x010_cyt =a5 conc_Ca_cyt x000_cyt -d5 x010_cyt -a4 conc_Ca_cyt x010_cyt -d4 x011_cyt ddtime x010_md =a5 conc_Ca_md x000_md -d5 x010_md -a4 conc_Ca_md x010_md -d4 x011_md x100_cyt =conc_IP3_cyt d1 x000_cyt x100_md =conc_IP3_md d1 x000_md x101_cyt =conc_IP3_cyt d3 x001_cyt x101_md =conc_IP3_md d3 x001_md x110_cyt =conc_IP3_cyt d1 x010_cyt x110_md =conc_IP3_md d1 x010_md x111_cyt =conc_IP3_cyt d3 x011_cyt x111_md =conc_IP3_md d3 x011_md x011_cyt =1 -x000_cyt +x001_cyt +x010_cyt +x100_cyt +x101_cyt +x110_cyt +x111_cyt x011_md =1 -x000_md +x001_md +x010_md +x100_md +x101_md +x110_md +x111_md

Component: SOC_Current_Ca

I_Ca_cyt_SOC=G_CaSOC_maxK_SOCK_SOC+conc_Ca_SRV_m_cyt-E_Ca_cytI_Ca_md_SOC=G_CaSOC_maxK_SOCK_SOC+conc_Ca_SRV_m_md-E_Ca_mdconc_Ca_SR=conc_Ca_IP3R+conc_Ca_RyR2

Component: SOC_Current_Na

I_Na_cyt_SOC =I_Ca_cyt_SOC z_Na 2 z_Ca 2 P_SOC_ratio conc_Na_cyt -conc_Na_out -z_Na F V_m_cyt R T conc_Ca_cyt -conc_Ca_out -z_Ca F V_m_cyt R T 1 --z_Ca F V_m_cyt R T 1 --z_Na F V_m_cyt R T I_Na_md_SOC =I_Ca_md_SOC z_Na 2 z_Ca 2 P_SOC_ratio conc_Na_md -conc_Na_out -z_Na F V_m_md R T conc_Ca_md -conc_Ca_out -z_Ca F V_m_md R T 1 --z_Ca F V_m_md R T 1 --z_Na F V_m_md R T

Component: Cl_flux_NaCl

J_cyt_NaCl =J_NaCl_max f_cyt E_Na_cyt -E_Cl_cyt 4 E_Na_cyt -E_Cl_cyt 4 +delta_CT 4 J_md_NaCl =J_NaCl_max f_md E_Na_md -E_Cl_md 4 E_Na_md -E_Cl_md 4 +delta_CT 4

Component: Cl_flux_KCl

J_cyt_KCl =J_KCl_max f_cyt E_K_cyt -E_Cl_cyt E_K_cyt -E_Cl_cyt +delta_CT J_md_KCl =J_KCl_max f_md E_K_md -E_Cl_md E_K_md -E_Cl_md +delta_CT

Component: Cl_flux_NKCC

J_cyt_NKCC =J_NKCC_max f_cyt conc_Na_out conc_K_out conc_Cl_out 2 -conc_Na_cyt conc_K_cyt conc_Cl_cyt 2 K_NKCC_Na K_NKCC_K K_NKCC_Cl1 K_NKCC_Cl2 conc_Na_out +K_NKCC_Na conc_K_out +K_NKCC_K conc_Cl_out +K_NKCC_Cl1 conc_Cl_out +K_NKCC_Cl2 conc_Na_cyt +K_NKCC_Na conc_K_cyt +K_NKCC_K conc_Cl_cyt +K_NKCC_Cl1 conc_Cl_cyt +K_NKCC_Cl2 J_md_NKCC =J_NKCC_max f_md conc_Na_out conc_K_out conc_Cl_out 2 -conc_Na_md conc_K_md conc_Cl_md 2 K_NKCC_Na K_NKCC_K K_NKCC_Cl1 K_NKCC_Cl2 conc_Na_out +K_NKCC_Na conc_K_out +K_NKCC_K conc_Cl_out +K_NKCC_Cl1 conc_Cl_out +K_NKCC_Cl2 conc_Na_md +K_NKCC_Na conc_K_md +K_NKCC_K conc_Cl_md +K_NKCC_Cl1 conc_Cl_md +K_NKCC_Cl2

Component: Cl_Channels_cyt

ddtime y_C1_cyt =-k_on conc_Ca_cyt y_C1_cyt -y_C0_cyt +k_off y_C2_cyt -y_C1_cyt -alpha1 y_C1_cyt +Beta_cyt y_O1_cyt ddtime y_C2_cyt =-k_on conc_Ca_cyt y_C2_cyt -y_C1_cyt +k_off y_C3_cyt -y_C2_cyt -alpha2 y_C2_cyt +Beta_cyt y_O2_cyt ddtime y_C3_cyt =k_on conc_Ca_cyt y_C2_cyt -k_off y_C3_cyt -alpha3 y_C3_cyt +Beta_cyt y_O3_cyt ddtime y_O1_cyt =alpha1 y_C1_cyt -Beta_cyt y_O1_cyt ddtime y_O2_cyt =alpha2 y_C2_cyt -Beta_cyt y_O2_cyt ddtime y_O3_cyt =alpha3 y_C3_cyt -Beta_cyt y_O3_cyt 1 =y_C0_cyt +y_C1_cyt +y_C2_cyt +y_C3_cyt +y_O1_cyt +y_O2_cyt +y_O3_cyt

Component: Cl_Channels_md

ddtime y_C1_md =-k_on conc_Ca_md y_C1_md -y_C0_md +k_off y_C2_md -y_C1_md -alpha1 y_C1_md +Beta_md y_O1_md ddtime y_C2_md =-k_on conc_Ca_md y_C2_md -y_C1_md +k_off y_C3_md -y_C2_md -alpha2 y_C2_md +Beta_md y_O2_md ddtime y_C3_md =k_on conc_Ca_md y_C2_md -k_off y_C3_md -alpha3 y_C3_md +Beta_md y_O3_md ddtime y_O1_md =alpha1 y_C1_md -Beta_md y_O1_md ddtime y_O2_md =alpha2 y_C2_md -Beta_md y_O2_md ddtime y_O3_md =alpha3 y_C3_md -Beta_md y_O3_md 1 =y_C0_md +y_C1_md +y_C2_md +y_C3_md +y_O1_md +y_O2_md +y_O3_md

Component: Cl_Channel_Rates

Beta_cyt =Lambda_Beta V1 +V2 V_m_cyt Beta_md =Lambda_Beta V1 +V2 V_m_md

Component: ClCa_Current

I_Cl_cyt_Ca =G_ClCa_max f_cyt y_O1_cyt +y_O2_cyt +y_O3_cyt V_m_cyt -E_Cl_cyt I_Cl_md_Ca =G_ClCa_max f_md y_O1_md +y_O2_md +y_O3_md V_m_md -E_Cl_md

Component: K_Inward_Rectifier_Current

I_K_cyt_ir =G_Kir_max f_cyt x_Kir_cyt conc_K_out conc_K_ref V_m_cyt -E_K_cyt I_K_md_ir =G_Kir_max f_md x_Kir_md conc_K_out conc_K_ref V_m_md -E_K_md x_Kir_cyt =alpha_Kir_cyt alpha_Kir_cyt +Beta_Kir_cyt x_Kir_md =alpha_Kir_md alpha_Kir_md +Beta_Kir_md

Component: K_ir_Rates

alpha_Kir_cyt =3.241 1 +V_m_cyt -E_K_cyt -185.72 1.596 alpha_Kir_md =3.241 1 +V_m_md -E_K_md -185.72 1.596 Beta_Kir_cyt =1 0.013 V_m_cyt -E_K_cyt +5.7922 0.294 +V_m_cyt -E_K_cyt -653.73299999999995 0.244 1 +-V_m_cyt -E_K_cyt -7.7775 0.0928 Beta_Kir_md =1 0.013 V_m_md -E_K_md +5.7922 0.294 +V_m_md -E_K_md -653.73299999999995 0.244 1 +-V_m_md -E_K_md -7.7775 0.0928

Component: K_Delayed_Rectifier_Current

I_K_cyt_v =G_Kv_max f_cyt P_Kv_cyt 2 V_m_cyt -E_K_cyt I_K_md_v =G_Kv_max f_md P_Kv_md 2 V_m_md -E_K_md

Component: Kv_Activations

P_Kv_cyt =0.58 P1_cyt +0.42 P2_cyt P_Kv_md =0.58 P1_md +0.42 P2_md ddtime P1_cyt =P_bar_Kv_cyt -P1_cyt Tau_P1_cyt ddtime P1_md =P_bar_Kv_md -P1_md Tau_P1_md ddtime P2_cyt =P_bar_Kv_cyt -P2_cyt Tau_P2_cyt ddtime P2_md =P_bar_Kv_md -P2_md Tau_P2_md P_bar_Kv_cyt =1 1 +-V_m_cyt +0.00177 0.01452 P_bar_Kv_md =1 1 +-V_m_md +0.00177 0.01452

Component: Kv_Time_Constants

Tau_P1_cyt =0.2109873 -V_m_cyt +0.21434 0.19535 2 -0.02059 Tau_P1_md =0.2109873 -V_m_md +0.21434 0.19535 2 -0.02059 Tau_P2_cyt =0.82139 -V_m_cyt +0.03159 0.02746 2 +0.00019 Tau_P2_md =0.82139 -V_m_md +0.03159 0.02746 2 +0.00019

Component: K_ATP_Current

I_K_cyt_ATP =G_KATP_max f_cyt conc_K_out conc_K_ref n_ATP V_m_cyt -E_K_cyt I_K_md_ATP =G_KATP_max f_md conc_K_out conc_K_ref n_ATP V_m_md -E_K_md

Component: K_Ca_Current

I_K_cyt_Ca =G_KCa_max f_cyt P_KCa_cyt V_m_cyt -E_K_cyt I_K_md_Ca =G_KCa_max f_md P_KCa_md V_m_md -E_K_md

Component: KCa_Activations

P_KCa_cyt =0.65 PF_cyt +0.35 PS_cyt P_KCa_md =0.65 PF_md +0.35 PS_md ddtime PF_cyt =P_bar_KCa_cyt -PF_cyt Tau_PF ddtime PF_md =P_bar_KCa_md -PF_md Tau_PF ddtime PS_cyt =P_bar_KCa_cyt -PS_cyt Tau_PS ddtime PS_md =P_bar_KCa_md -PS_md Tau_PS P_bar_KCa_cyt =1 1 +-V_m_cyt -V_KCa_cyt 0.0308 P_bar_KCa_md =1 1 +-V_m_md -V_KCa_md 0.0308

Component: NaK_alpha_Current

I_NaK_cyt_alpha1 =I_NaK_alpha1_max Psi_NaK_cyt conc_K_out conc_K_out +K_m_K conc_Na_cyt 1.5 conc_Na_cyt 1.5 +K_m_Na_alpha1 1.5 I_NaK_md_alpha2=I_NaK_alpha1_maxPsi_NaK_mdconc_K_outconc_K_out+K_m_Kconc_Na_md1.5conc_Na_md1.5+K_m_Na_alpha21.5 I_NaK_alpha2_max=I_NaK_alpha1_max

Component: NaK_alpha_Waveform

Psi_NaK_cyt=11+0.1245-0.1V_m_cytFRT+0.0365Sigma-V_m_cytFRTPsi_NaK_md=11+0.1245-0.1V_m_mdFRT+0.0365Sigma-V_m_mdFRTSigma=17conc_Na_out67.299999999999997-1

Component: NCX_Current

I_NaCa_md =I_NaCa_max K_a_NaCa Psi_F conc_Na_md 3 conc_Ca_out -Psi_R conc_Na_out 3 conc_Ca_md G 1 +k_sat Psi_R

Component: NCX_Waveform

K_a_NaCa =1 1 +K_mCa_act conc_Ca_md 2 Psi_F =gamma V_m_md F R T Psi_R =gamma -1 V_m_md F R T G =conc_Na_out 3 conc_Ca_md +conc_Na_md 3 conc_Ca_out +K_m_Nao 3 conc_Ca_md +K_m_Cao conc_Na_md 3 +K_m_Nai 3 conc_Ca_out 1 +conc_Ca_md K_m_Cai +K_m_Cai conc_Na_out 3 1 +conc_Na_md K_m_Nai 3

Component: VONa_Current

I_VONa_cyt =G_VONa_max f_cyt m_VONa_cyt h_VONa_cyt V_m_cyt -E_Na_cyt I_VONa_md =G_VONa_max f_md m_VONa_md h_VONa_md V_m_md -E_Na_md

Component: VONa_Channels

ddtime m_VONa_cyt =m_bar_cyt -m_VONa_cyt Tau_m ddtime m_VONa_md =m_bar_md -m_VONa_md Tau_m ddtime h_VONa_cyt =h_bar_cyt -h_VONa_cyt Tau_h ddtime h_VONa_md =h_bar_md -h_VONa_md Tau_h m_bar_cyt =1 1 +-V_m_cyt +0.0239 0.0048 m_bar_md =1 1 +-V_m_md +0.0239 0.0048 h_bar_cyt =1 1 +V_m_cyt +0.0661 0.0065 h_bar_md =1 1 +V_m_md +0.0661 0.0065

Component: Ca_Pump_Current

I_Ca_cyt_P =I_CaP_max f_cyt conc_Ca_cyt K_m_CaP +conc_Ca_cyt I_Ca_md_P =I_CaP_max f_md conc_Ca_md K_m_CaP +conc_Ca_md

Component: CaL_Current

I_Ca_cyt_L =G_CaL_max f_cyt d_L_cyt f_L_cyt V_m_cyt -E_Ca_cyt I_Ca_md_L =G_CaL_max f_md d_L_md f_L_md V_m_md -E_Ca_md

Component: CaL_Activations

f_L_cyt =0.74 f_F_cyt +0.26 f_L_md =0.74 f_F_md +0.26 ddtime d_L_cyt =d_bar_L_cyt -d_L_cyt Tau_d_cyt ddtime d_L_md =d_bar_L_md -d_L_md Tau_d_md ddtime f_F_cyt =f_bar_F_cyt -f_F_cyt Tau_f_cyt ddtime f_F_md =f_bar_F_md -f_F_md Tau_f_md d_bar_L_cyt =1 1 +-V_m_cyt +0.00188 0.00757 d_bar_L_md =1 1 +-V_m_md +0.00188 0.00757 f_bar_F_md =1 1 +V_m_md +0.02932 0.00154 f_bar_F_cyt =1 1 +V_m_cyt +0.02932 0.00154

Component: CaL_Time_Constants

Tau_d_cyt =0.00289 -V_m_cyt +0.00863 0.01239 2 +0.00243 Tau_d_md =0.0289 -V_m_md +0.00863 0.01239 2 +0.00243 Tau_f_cyt =0.29559 -V_m_cyt -0.00472 0.11255 2 +0.02319 Tau_f_md =0.29559 -V_m_md -0.00472 0.11255 2 +0.02319

Component: Calseq_Kinetics

ddtime conc_CalseqCa_IP3R =k_Calseq_on conc_Ca_IP3R conc_Calseq_SR_tot -conc_CalseqCa_IP3R -k_Calseq_off conc_CalseqCa_IP3R ddtime conc_CalseqCa_RyR =k_Calseq_on conc_Ca_RyR conc_Calseq_SR_tot -conc_CalseqCa_RyR -k_Calseq_off conc_CalseqCa_RyR

Component: Md_Cyt_Flux_K

J_Kdiff=Volsecconc_K_md-conc_K_cytifV_m_md-V_m_cyt=0VolsecXi_Kconc_K_md-conc_K_cyt-Xi_K1--Xi_KotherwiseXi_K=z_KFRTV_m_md-V_m_cytP_K_diff=hconD_KLVolsec=DiffAreaP_K_diff

Component: Md_Cyt_Flux_Na

J_Nadiff=Volsecconc_Na_md-conc_Na_cytifV_m_md-V_m_cyt=0VolsecXi_Naconc_Na_md-conc_Na_cyt-Xi_Na1--Xi_NaotherwiseXi_Na=z_NaFRTV_m_md-V_m_cytP_Na_diff=hconD_NaLVolsec=DiffAreaP_Na_diff

Component: Md_Cyt_Flux_Ca

J_Cadiff=Volsecconc_Ca_md-conc_Ca_cytifV_m_md-V_m_cyt=0VolsecXi_Caconc_Ca_md-conc_Ca_cyt-Xi_Ca1--Xi_CaotherwiseXi_Ca=z_CaFRTV_m_md-V_m_cytP_Ca_diff=hconD_CaLVolsec=DiffAreaP_Ca_diff

Component: Md_Cyt_Flux_Cl

J_Cldiff=Volsecconc_Cl_md-conc_Cl_cytifV_m_md-V_m_cyt=0VolsecXi_Clconc_Cl_md-conc_Cl_cyt-Xi_Cl1--Xi_ClotherwiseXi_Cl=z_ClFRTV_m_md-V_m_cytP_Cl_diff=hconD_ClLVolsec=DiffAreaP_Cl_diff

Component: Membrane_Potential_Cyt

I_cyt_tot =I_K_cyt_b +I_K_cyt_ir +I_K_cyt_ATP +I_K_cyt_v +I_K_cyt_Ca -F J_Kdiff +I_NaK_cyt_alpha1 +I_Na_cyt_b +I_VONa_cyt +I_Na_cyt_SOC -F J_Nadiff +I_Cl_cyt_b +I_Cl_cyt_Ca +F J_Cldiff +I_Ca_cyt_b +I_Ca_cyt_P +I_Ca_cyt_L +I_Ca_cyt_SOC -2 F J_Cadiff +I_cyt_up -I_cyt_RyR -I_cyt_IP3R I_cyt_up =I_cyt_SERCa_IP3R +I_cyt_SERCa_RyR

Component: Membrane_Voltage

ddtime V_m_cyt =-I_cyt_tot f_cyt C_m ddtime V_m_md =-I_md_tot f_md C_m

Component: Membrane_Potential_md

I_md_tot =I_K_md_b +I_K_md_ir +I_K_md_ATP +I_K_md_v +I_K_md_Ca +F J_Kdiff +I_NaK_md_alpha2 +I_Na_md_b +I_VONa_md +I_Na_md_SOC +I_NaCa_md +F J_Nadiff +I_Cl_md_b +I_Cl_md_Ca -F J_Cldiff +I_Ca_md_b +I_Ca_md_P +I_Ca_md_L +I_Ca_md_SOC +2 F J_Cadiff +I_md_up -I_md_RyR -I_md_IP3R I_md_up =I_md_SERCa_IP3R +I_md_SERCa_RyR

Component: Parameters_VR

Component: MLCK_Activation_CaM

ddtimeconc_CaM___=-k1_CAM_onconc_Ca_cyt2conc_CaM___+k1_CAM_offconc_CaM_C_-k2_CAM_onconc_Ca_cyt2conc_CaM___+k2_CAM_offconc_CaMN__ddtimeconc_CaMN__=k2_CAM_onconc_Ca_cyt2conc_CaM___-k2_CAM_offconc_CaMN__-k3_CAM_onconc_Ca_cyt2conc_CaMN__+k3_CAM_offconc_CaMNC_ddtimeconc_CaM_C_=k1_CAM_onconc_Ca_cyt2conc_CaM___-k1_CAM_offconc_CaM_C_-k4_CAM_onconc_Ca_cyt2conc_CaM_C_+k4_CAM_offconc_CaMNC_-k5_CAM_onconc_MLCK_freeconc_CaM_C_+k5_CAM_offconc_CaM_CMddtimeconc_CaM_CM=k5_CAM_onconc_MLCK_freeconc_CaM_C_-k5_CAM_offconc_CaM_CM-k6_CAM_onconc_Ca_cyt2conc_CaM_CM+k6_CAM_offconc_CaMNCMddtimeconc_CaMNCM=k6_CAM_onconc_Ca_cyt2conc_CaM_CM-k6_CAM_offconc_CaMNCM+k7_CAM_onconc_MLCK_freeconc_CaMNC_-k7_CAM_offconc_CaMNCMconc_CaMNC_=conc_CaM_tot-conc_CaM___+conc_CaMN__+conc_CaM_C_+conc_CaM_CM+conc_CaMNCMconc_MLCK_free=conc_MLCK_tot-conc_CaM_CM+conc_CaMNCM

Component: MLCK_Phospho_Myosin

ddtimeconc_Myo=-k1_6_Myoconc_Myo+k2_Myoconc_MyoP+k7_Myoconc_AMyoddtimeconc_MyoP=k1_6_Myoconc_Myo-k2_Myoconc_MyoP-k3_Myoconc_MyoP+k4_Myoconc_AMyoPddtimeconc_AMyo=-k7_Myoconc_AMyo+k5_Myoconc_AMyoP-k1_6_Myoconc_AMyoconc_AMyoP=conc_Myo_tot-conc_Myo+conc_MyoP+conc_AMyo

Component: MLCP_Rate_Constants

k1_6_Myo=k_Myo_MLCKconc_CaMNCMk2_Myo=k_Myo_MLCP1+3.65conc_cGMP2K_M_cGMP2+conc_cGMP2R_MLCPk5_Myo=k2_MyoR_MLCP=1+alpha_MLCPconc_H2O2_cyt2K_MLCP_H2O22+conc_H2O2_cyt2

Component: NO_Kinetics

ddtimeconc_NO_cyt=Area2vol_cytJ_NO-P_NOconc_NO_cyt-k_NO_O2conc_NO_cytconc_O2_cyt-k_NO_consconc_NO_cyt

Component: O2_Kinetics

ddtimeconc_O2_cyt=G_O2-Area2vol_cytP_O2conc_O2_cytzeta1--zeta-k_NO_O2conc_NO_cytconc_O2_cyt-k_SOD_O2conc_SODconc_O2_cyt

Component: Zeta

zeta=-FV_m_cytRT

Component: H2O2_Kinetics

ddtimeconc_H2O2_cyt=Area2vol_cytJ_H2O2-P_H2O2conc_H2O2_cyt+k_SOD_O2conc_SODconc_O2_cyt-k_catconc_H2O2_cytconc_CAT

Component: sGC_Activation

ddtimeconc_Eb=-k_sGC_forconc_Ebconc_NO_cyt+k_sGC_backconc_E6c+k4_sGC_forconc_E5cddtimeconc_E5c=k3_sGC_forconc_E6cconc_NO_cyt+k2_sGC_forconc_E6c-k4_sGC_forconc_E5cconc_E6c=1-conc_Eb+conc_E5ck4_sGC_for=K4_sGCconc_cGMP

Component: cGMP_Kinetics

ddtimeconc_cGMP=V_sGC_maxconc_E5c-k_PDEconc_cGMP2K_M_PDE+conc_cGMP

Component: cGMP_Ca_Interaction

V_KCa_cyt=V_Calog10conc_Ca_cyt1-V_B-V_KCa_NOconc_NO_cytK_KCa_NO+conc_NO_cyt-V_KCa_cGMPconc_cGMP2K_KCa_cGMP2+conc_cGMP2convertV_KCa_md=V_Calog10conc_Ca_md1-V_B-V_KCa_NOconc_NO_cytK_KCa_NO+conc_NO_cyt-V_KCa_cGMPconc_cGMP2K_KCa_cGMP2+conc_cGMP2convert

Component: SERCA_Activation

K_mf=K_mf_restK_mod1-conc_cGMP2K_SERCA+conc_cGMP

Component: CaP_Current_Max

I_CaP_max=I_PMCA_rest_max1+3.6conc_cGMPK_PMCA+conc_cGMP

Component: ROS_SERCA_Interaction

theta_SERCA_O2=conc_O2_cytK_SERCA_O2+conc_O2_cyttheta_SERCA_H2O2=conc_H2O2_cytK_SERCA_H2O2+conc_H2O2_cytI_SERCA_max=I_SERCA_rest_max1-theta_SERCA_O21-theta_SERCA_H2O2

Component: ROS_NaK_Interaction

I_NaK_alpha1_max=I_NaK_alpha1_rest_max1-conc_O2_cytK_NaK_O2+conc_O2_cyt1-conc_H2O2_cytK_NaK_H2O2+conc_H2O2_cyt

Component: MLCK_Activation_CaM_md

ddtimeconc_CaM___=-k1_CAM_onconc_Ca_md2conc_CaM___+k1_CAM_offconc_CaM_C_-k2_CAM_onconc_Ca_md2conc_CaM___+k2_CAM_offconc_CaMN__ddtimeconc_CaMN__=k2_CAM_onconc_Ca_md2conc_CaM___-k2_CAM_offconc_CaMN__-k3_CAM_onconc_Ca_md2conc_CaMN__+k3_CAM_offconc_CaMNC_ddtimeconc_CaM_C_=k1_CAM_onconc_Ca_md2conc_CaM___-k1_CAM_offconc_CaM_C_-k4_CAM_onconc_Ca_md2conc_CaM_C_+k4_CAM_offconc_CaMNC_-k5_CAM_onconc_MLCK_freeconc_CaM_C_+k5_CAM_offconc_CaM_CMddtimeconc_CaM_CM=k5_CAM_onconc_MLCK_freeconc_CaM_C_-k5_CAM_offconc_CaM_CM-k6_CAM_onconc_Ca_md2conc_CaM_CM+k6_CAM_offconc_CaMNCMddtimeconc_CaMNCM=k6_CAM_onconc_Ca_md2conc_CaM_CM-k6_CAM_offconc_CaMNCM+k7_CAM_onconc_MLCK_freeconc_CaMNC_-k7_CAM_offconc_CaMNCMconc_CaMNC_=conc_CaM_tot-conc_CaM___+conc_CaMN__+conc_CaM_C_+conc_CaM_CM+conc_CaMNCMconc_MLCK_free=conc_MLCK_tot-conc_CaM_CM+conc_CaMNCM

Component: Contractile_Force

f_p=conc_MyoP+conc_AMyoPconc_Myo_totF_contract=a-bf_p+cf_p2-df_p3