Model Mathematics

Component: environment

Component: membrane

dd time V =- 1.0 Cm ⁢ 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 ⁢ T 2.0 ⁢ F ⁢ln⁡ Ca_o Ca_i

Component: L_type_calcium_current_d_gate

d_infinity = 1.0 1.0 +ⅇ- V + 1.878 7.5704 tau_d = 2.8928 ⁢ⅇ- V + 8.6344 12.3884 2.0 + 2.4343 dd time d = d_infinity - d tau_d

Component: L_type_calcium_current_f_gate

f_infinity = 1.0 1.0 +ⅇ V + 29.3188 1.5389 tau_ff = 295.5937 ⁢ⅇ- V - 4.7187 112.545 2.0 + 23.1907 dd time ff = f_infinity - ff tau_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 ⁢ T F ⁢ln⁡ K_o K_i

Component: calcium_activated_potassium_current_P_K_Ca_gate

Po_infinity = 1.0 1.0 +ⅇ- V - V_1_2_K_Ca 21.70 V_1_2_K_Ca =- 45.0 ⁢log10⁡ Ca_i - 198.55 Pf_infinity = Ps_infinity - Po_infinity dd time Pf = Pf_infinity - Pf tau_Pf dd time Ps = Ps_infinity - Ps tau_Ps P_K_Ca = 0.65 ⁢ Pf + 0.35 ⁢ Ps

Component: delayed_rectifier_current

i_K = g_K ⁢ Pk 2.0 ⁢ V - E_K

Component: delayed_rectifier_current_Pk_gate

Pk_infinity = 1.0 1.0 +ⅇ- V - V_1_2 k tau_P1 = 210.9873 ⁢ⅇ- V + 214.3355 195.3502 2.0 - 20.5866 tau_P2 = 821.3949 ⁢ⅇ- V + 31.5891 27.4568 2.0 - 0.1892 dd time P2 = Pk_infinity - P2 tau_P2 dd time P1 = Pk_infinity - P1 tau_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 ⁢ F 2.0 ⁢ V R ⁢ T ⁢ P_K⁢ K_o - K_i ⁢ⅇ F ⁢ V R ⁢ T 1.0 -ⅇ F ⁢ V R ⁢ T i_M_Na = Am ⁢ Pm ⁢ F 2.0 ⁢ V R ⁢ T ⁢ P_Na ⁢ Na_o - Na_i ⁢ⅇ F ⁢ V R ⁢ T 1.0 -ⅇ F ⁢ V R ⁢ T i_M_Ca = Am ⁢ Pm ⁢ F 2.0 ⁢ V R ⁢ T ⁢ 2.0 2.0 ⁢ P_Ca ⁢ Ca_o - Ca_i ⁢ⅇ F ⁢ V R ⁢ T 1.0 -ⅇ F ⁢ V R ⁢ T Pm = 1.0 1.0 +ⅇ- sigma - sigma_1_2 k_sigma

Component: inward_rectifier_current

i_Ki = g_max_Ki ⁢ V - E_K 1.0 +ⅇ- V - V_1_2_Ki 28.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 ⁢ CaCM CaCM + K_CaCM

Component: sodium_potassium_pump

i_NaK = I_NaK ⁢ K_o K_o + Km_Ko ⁢ Na_i 1.5 Na_i 1.5 + Km_Nai 1.5 ⁢ V + 150.0 V + 200.0

Component: sodium_calcium_exchanger

i_NaCa = g_NaCa ⁢ Na_i 3.0 ⁢ Ca_o ⁢ phi_F - Na_o 3.0 ⁢ Ca_i 1.0 + Na_o 3.0 ⁢ Ca_i ⁢ phi_R + Na_i 3.0 ⁢ Ca_o phi_F =ⅇ gamma ⁢ V ⁢ F R ⁢ T phi_R =ⅇ 1.0 - gamma ⁢ V ⁢ F R ⁢ 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 ⁢ T F ⁢ln⁡ Na_o Na_i E_Ca = R ⁢ T 2.0 ⁢ F ⁢ln⁡ Ca_o Ca_i

Component: extracellular_ionic_concentrations

Component: intracellular_ionic_concentrations

dd time Na_i =- 3.0 ⁢ i_NaK + 3.0 ⁢ i_NaCa + i_B_Na + i_M_Na F ⁢ Vol_i dd time K_i =- i_K + i_K_Ca + i_Ki + i_B_K + i_M_K - 2.0 ⁢ i_NaK F ⁢ Vol_i dd time Ca_i =- i_Ca_L + i_CaP + i_B_Ca + i_M_Ca + i_up - 2.0 ⁢ i_NaCa + i_rel 2.0 ⁢ F ⁢ Vol_Ca + delta_S_CM_cal + delta_B_F_cal dd time S_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 dd time B_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_i Ca_i + Km_up i_tr = Ca_u - Ca_r ⁢ 2.0 ⁢ F ⁢ Vol_u tau_tr i_rel = R10 2.0 ⁢ Ca_u - Ca_r ⁢ 2.0 ⁢ F ⁢ Vol_r tau_rel dd time Ca_up = i_up - i_tr 2.0 ⁢ F ⁢ Vol_u dd time Ca_r = i_tr - i_rel 2.0 ⁢ F ⁢ Vol_r

Component: CICR_mechanism_of_the_SR

R00 = Kr1_ ⁢ R10 + Kr2_ ⁢ R01 - Kr1 ⁢ Ca_i 2.0 + Kr2 ⁢ Ca_i ⁢ R00 R10 = Kr1 ⁢ Ca_i 2.0 ⁢ R00 + Kr2_ ⁢ R11 - Kr2 ⁢ Ca_i + Kr1_ ⁢ R10 R11 = Kr1 ⁢ Ca_i 2.0 ⁢ R01 + Kr2 ⁢ Ca_i ⁢ R10 - Kr2_ + Kr1_ ⁢ R11 R01 = Kr1_ ⁢ R11 + Kr2 ⁢ Ca_i ⁢ R00 - Kr1 ⁢ Ca_i 2.0 + 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 = CaCM 2.0 CaCM 2.0 + K_CaCM 2.0

Component: passive_force

Fp = kp ⁢ⅇ alpha_p ⁢ lc - l0 l0 - 1.0

Component: cross_bridge_elasticity

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

Component: active_force_generation

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

Component: series_viscoelastic_force

Fs = mu_s ⁢ delta_ls + ks ⁢ⅇ alpha_s ⁢ ls - ls0 ls0 - 1.0 dd time ls = delta_ls