Model Mathematics

Component: environment

Component: membrane

ddtimeV=-1Cm⁢i_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_L⁢d⁢f⁢V-E_Ca_L E_Ca_L=R⁢T2⁢F⁢ln⁡Ca_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.74⁢ff+0.26

Component: calcium_activated_potassium_current

i_K_Ca=g_K_Ca⁢P_K_Ca⁢V-E_K E_K=R⁢TF⁢ln⁡K_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=-45⁢log10⁡Ca_i-198.55 Pf_infinity=Ps_infinity-Po_infinity ddtimePf=Pf_infinity-Pftau_Pf ddtimePs=Ps_infinity-Pstau_Ps P_K_Ca=0.65⁢Pf+0.35⁢Ps

Component: delayed_rectifier_current

i_K=g_K⁢Pk2⁢V-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.58⁢P1+0.42⁢P2

Component: stretch_sensitive_current

i_M=i_M_K+i_M_Na+i_M_Ca i_M_K=Am⁢Pm⁢F2⁢VR⁢T⁢P_K⁢K_o-K_i⁢ⅇF⁢VR⁢T1-ⅇF⁢VR⁢T i_M_Na=Am⁢Pm⁢F2⁢VR⁢T⁢P_Na⁢Na_o-Na_i⁢ⅇF⁢VR⁢T1-ⅇF⁢VR⁢T i_M_Ca=Am⁢Pm⁢F2⁢VR⁢T⁢22⁢P_Ca⁢Ca_o-Ca_i⁢ⅇF⁢VR⁢T1-ⅇF⁢VR⁢T Pm=11+ⅇ-sigma-sigma_1_2k_sigma

Component: inward_rectifier_current

i_Ki=g_max_Ki⁢V-E_K1+ⅇ-V-V_1_2_Ki28.89 g_max_Ki=G_Ki⁢K_o⁢n_Ki V_1_2_Ki=A⁢log10⁡K_i+B

Component: calcium_pump

i_CaP=I_Ca_P⁢CaCMCaCM+K_CaCM

Component: sodium_potassium_pump

i_NaK=I_NaK⁢K_oK_o+Km_Ko⁢Na_i1.5Na_i1.5+Km_Nai1.5⁢V+150V+200

Component: sodium_calcium_exchanger

i_NaCa=g_NaCa⁢Na_i3⁢Ca_o⁢phi_F-Na_o3⁢Ca_i1+Na_o3⁢Ca_i⁢phi_R+Na_i3⁢Ca_o phi_F=ⅇgamma⁢V⁢FR⁢T phi_R=ⅇ1-gamma⁢V⁢FR⁢T

Component: background_currents

i_B=i_B_Na+i_B_Ca+i_B_K i_B_Na=g_Nab⁢V-E_Na i_B_Ca=g_Cab⁢V-E_Ca i_B_K=g_Kb⁢V-E_K E_Na=R⁢TF⁢ln⁡Na_oNa_i E_Ca=R⁢T2⁢F⁢ln⁡Ca_oCa_i

Component: extracellular_ionic_concentrations

Component: intracellular_ionic_concentrations

ddtimeNa_i=-3⁢i_NaK+3⁢i_NaCa+i_B_Na+i_M_NaF⁢Vol_i ddtimeK_i=-i_K+i_K_Ca+i_Ki+i_B_K+i_M_K-2⁢i_NaKF⁢Vol_i ddtimeCa_i=-i_Ca_L+i_CaP+i_B_Ca+i_M_Ca+i_up-2⁢i_NaCa+i_rel2⁢F⁢Vol_Ca+delta_S_CM_cal+delta_B_F_cal ddtimeS_CM_cal=delta_S_CM_cal delta_S_CM_cal=k_1⁢CaS_CM-k1⁢Ca_i⁢S_CM_cal CaS_CM=CaCM CaCM=S_CM_Ca-S_CM_cal ddtimeB_F_cal=delta_B_F_cal delta_B_F_cal=k_d⁢CaB_F-kd⁢Ca_i⁢B_F_cal CaB_F=B_F_Ca-B_F_cal

Component: calcium_fluxes_in_the_SR

i_up=I_up⁢Ca_iCa_i+Km_up i_tr=Ca_u-Ca_r⁢2⁢F⁢Vol_utau_tr i_rel=R102⁢Ca_u-Ca_r⁢2⁢F⁢Vol_rtau_rel ddtimeCa_up=i_up-i_tr2⁢F⁢Vol_u ddtimeCa_r=i_tr-i_rel2⁢F⁢Vol_r

Component: CICR_mechanism_of_the_SR

R00=Kr1_⁢R10+Kr2_⁢R01-Kr1⁢Ca_i2+Kr2⁢Ca_i⁢R00 R10=Kr1⁢Ca_i2⁢R00+Kr2_⁢R11-Kr2⁢Ca_i+Kr1_⁢R10 R11=Kr1⁢Ca_i2⁢R01+Kr2⁢Ca_i⁢R10-Kr2_+Kr1_⁢R11 R01=Kr1_⁢R11+Kr2⁢Ca_i⁢R00-Kr1⁢Ca_i2+Kr2_⁢R01

Component: smooth_muscle_contraction

M=K2⁢Mp+K7⁢AM-K1⁢M Mp=K1⁢M+K4⁢AMp-K2+K3⁢Mp AMp=K6⁢AM+K3⁢Mp-K5+K4⁢AMp AM=K5⁢AMp-K7+K6⁢AM K1=K6 K6=CaCM2CaCM2+K_CaCM2

Component: passive_force

Fp=kp⁢ⅇalpha_p⁢lc-l0l0-1

Component: cross_bridge_elasticity

Fx=kx1⁢AMp+kx2⁢AM⁢lx⁢ⅇ-beta⁢la-loptlopt2

Component: active_force_generation

Fa=fAMp⁢AMp⁢vx+ia+fAM⁢AM⁢ia⁢ⅇ-beta⁢la-loptlopt2

Component: series_viscoelastic_force

Fs=mu_s⁢delta_ls+ks⁢ⅇalpha_s⁢ls-ls0ls0-1 ddtimels=delta_ls