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