Model Mathematics

Component: environment

Component: membrane

ddtimeV=-1Cmi_Ca_L+i_K+i_K_Ca+i_Ki+i_M+i_NaCa+i_NaK+i_CaP+i_B

Component: L_type_calcium_current

i_Ca_L=g_Ca_LdfV-E_Ca_L E_Ca_L=RT2FlnCa_oCa_i

Component: L_type_calcium_current_d_gate

d_infinity=11+-V+1.8787.5704 tau_d=2.8928-V+8.634412.38842+2.4343 ddtimed=d_infinity-dtau_d

Component: L_type_calcium_current_f_gate

f_infinity=11+V+29.31881.5389 tau_ff=295.5937-V-4.7187112.5452+23.1907 ddtimeff=f_infinity-fftau_ff f=0.74ff+0.26

Component: calcium_activated_potassium_current

i_K_Ca=g_K_CaP_K_CaV-E_K E_K=RTFlnK_oK_i

Component: calcium_activated_potassium_current_P_K_Ca_gate

Po_infinity=11+-V-V_1_2_K_Ca21.7 V_1_2_K_Ca=-45log10Ca_i-198.55 Pf_infinity=Ps_infinity-Po_infinity ddtimePf=Pf_infinity-Pftau_Pf ddtimePs=Ps_infinity-Pstau_Ps P_K_Ca=0.65Pf+0.35Ps

Component: delayed_rectifier_current

i_K=g_KPk2V-E_K

Component: delayed_rectifier_current_Pk_gate

Pk_infinity=11+-V-V_1_2k tau_P1=210.9873-V+214.3355195.35022-20.5866 tau_P2=821.3949-V+31.589127.45682-0.1892 ddtimeP2=Pk_infinity-P2tau_P2 ddtimeP1=Pk_infinity-P1tau_P1 Pk=0.58P1+0.42P2

Component: stretch_sensitive_current

i_M=i_M_K+i_M_Na+i_M_Ca i_M_K=AmPmF2VRTP_KK_o-K_iFVRT1-FVRT i_M_Na=AmPmF2VRTP_NaNa_o-Na_iFVRT1-FVRT i_M_Ca=AmPmF2VRT22P_CaCa_o-Ca_iFVRT1-FVRT Pm=11+-sigma-sigma_1_2k_sigma

Component: inward_rectifier_current

i_Ki=g_max_KiV-E_K1+-V-V_1_2_Ki28.89 g_max_Ki=G_KiK_on_Ki V_1_2_Ki=Alog10K_i+B

Component: calcium_pump

i_CaP=I_Ca_PCaCMCaCM+K_CaCM

Component: sodium_potassium_pump

i_NaK=I_NaKK_oK_o+Km_KoNa_i1.5Na_i1.5+Km_Nai1.5V+150V+200

Component: sodium_calcium_exchanger

i_NaCa=g_NaCaNa_i3Ca_ophi_F-Na_o3Ca_i1+Na_o3Ca_iphi_R+Na_i3Ca_o phi_F=gammaVFRT phi_R=1-gammaVFRT

Component: background_currents

i_B=i_B_Na+i_B_Ca+i_B_K i_B_Na=g_NabV-E_Na i_B_Ca=g_CabV-E_Ca i_B_K=g_KbV-E_K E_Na=RTFlnNa_oNa_i E_Ca=RT2FlnCa_oCa_i

Component: extracellular_ionic_concentrations

Component: intracellular_ionic_concentrations

ddtimeNa_i=-3i_NaK+3i_NaCa+i_B_Na+i_M_NaFVol_i ddtimeK_i=-i_K+i_K_Ca+i_Ki+i_B_K+i_M_K-2i_NaKFVol_i ddtimeCa_i=-i_Ca_L+i_CaP+i_B_Ca+i_M_Ca+i_up-2i_NaCa+i_rel2FVol_Ca+delta_S_CM_cal+delta_B_F_cal ddtimeS_CM_cal=delta_S_CM_cal delta_S_CM_cal=k_1CaS_CM-k1Ca_iS_CM_cal CaS_CM=CaCM CaCM=S_CM_Ca-S_CM_cal ddtimeB_F_cal=delta_B_F_cal delta_B_F_cal=k_dCaB_F-kdCa_iB_F_cal CaB_F=B_F_Ca-B_F_cal

Component: calcium_fluxes_in_the_SR

i_up=I_upCa_iCa_i+Km_up i_tr=Ca_u-Ca_r2FVol_utau_tr i_rel=R102Ca_u-Ca_r2FVol_rtau_rel ddtimeCa_up=i_up-i_tr2FVol_u ddtimeCa_r=i_tr-i_rel2FVol_r

Component: CICR_mechanism_of_the_SR

R00=Kr1_R10+Kr2_R01-Kr1Ca_i2+Kr2Ca_iR00 R10=Kr1Ca_i2R00+Kr2_R11-Kr2Ca_i+Kr1_R10 R11=Kr1Ca_i2R01+Kr2Ca_iR10-Kr2_+Kr1_R11 R01=Kr1_R11+Kr2Ca_iR00-Kr1Ca_i2+Kr2_R01

Component: smooth_muscle_contraction

M=K2Mp+K7AM-K1M Mp=K1M+K4AMp-K2+K3Mp AMp=K6AM+K3Mp-K5+K4AMp AM=K5AMp-K7+K6AM K1=K6 K6=CaCM2CaCM2+K_CaCM2

Component: passive_force

Fp=kpalpha_plc-l0l0-1

Component: cross_bridge_elasticity

Fx=kx1AMp+kx2AMlx-betala-loptlopt2

Component: active_force_generation

Fa=fAMpAMpvx+ia+fAMAMia-betala-loptlopt2

Component: series_viscoelastic_force

Fs=mu_sdelta_ls+ksalpha_sls-ls0ls0-1 ddtimels=delta_ls