Model Mathematics

Component based math viewer is available

Component: environment

Component: membrane

past = ⌊ \frac{time}{stim_period} ⌋ stim_period i_Stim = \{\begin{matrix} stim_amplitude if time - past \geq stim_offset \land time - past ≦ stim_offset + stim_duration \\ 0 otherwise \end{matrix} \frac{d}{d time} (V) = - (i_CaL + i_pCa + i_NaCa + i_Cab + i_Na + i_Nab + i_NaK + i_Kto_f + i_Kto_s + i_K1 + i_Ks + i_Kur + i_Kss + i_Kr + i_ClCa + i_Stim)

Component: calcium_concentration

\frac{d}{d time} (Cai) = Bi (J_leak + J_xfer - (\frac{(i_Cab + i_pCa - (2 i_NaCa)) Acap Cm}{2 Vmyo F}) - J_serca) \frac{d}{d time} (Cass) = Bss (\frac{J_rel VJSR}{Vss} - (\frac{J_xfer Vmyo}{Vss} + \frac{i_CaL Acap Cm}{2 Vss F})) \frac{d}{d time} (CaJSR) = BJSR (J_tr - J_rel) \frac{d}{d time} (CaNSR) = \frac{(J_serca - J_leak) Vmyo}{VNSR} - (\frac{J_tr VJSR}{VNSR}) Bi = {(1 + \frac{Bmax Kd}{{(Kd + Cai)}^{2}})}^{- 1} Bss = {(1 + \frac{Bmax Kd}{{(Kd + Cass)}^{2}})}^{- 1} BJSR = {(1 + \frac{CSQN_tot Km_CSQN}{{(Km_CSQN + CaJSR)}^{2}})}^{- 1}

Component: calcium_fluxes

J_rel = v1 (P_O1 + P_O2) (CaJSR - Cass) P_RyR J_tr = \frac{CaNSR - CaJSR}{tau_tr} J_xfer = \frac{Cass - Cai}{tau_xfer} J_leak = v2 (CaNSR - Cai) CaMKb = \frac{0.05 (1 - CaMKt) 1}{1 + \frac{0.7}{Cass}} \frac{d}{d time} (CaMKt) = on_rate CaMKb (CaMKb + CaMKt) - (off_rate CaMKt) CaMKa = CaMKb + CaMKt vmup = (\frac{3.1512 {CaMKa}^{2.2062}}{{0.1588}^{2.2062} + {CaMKa}^{2.2062}} + 1) vmup_init J_serca = \frac{vmup {Cai}^{2}}{{Km_up}^{2} + {Cai}^{2}} \frac{d}{d time} (P_RyR) = P_ryr_const1 P_RyR + \frac{P_ryr_const2 i_CaL}{i_CaL_max} ⅇ^{\frac{- ({(V - 5)}^{2})}{648}}

Component: ryanodine_receptors

\frac{d}{d time} (P_O1) = k_plus_a {Cass}^{n} P_C1 + k_minus_b P_O2 + k_minus_c P_C2 - (k_minus_a P_O1 + k_plus_b {Cass}^{m} P_O1 + k_plus_c P_O1) P_C1 = 1 - (P_C2 + P_O1 + P_O2) \frac{d}{d time} (P_O2) = k_plus_b {Cass}^{m} P_O1 - (k_minus_b P_O2) \frac{d}{d time} (P_C2) = k_plus_c P_O1 - (k_minus_c P_C2)

Component: L_type_calcium_current

FVRT = \frac{F V}{R T} FVRT_Ca = 2 FVRT expVL = ⅇ^{\frac{V - V_L}{delta_V_L}} alpha_p = \frac{expVL}{t_L (expVL + 1)} alpha_m = \frac{phi_L}{t_L} epsilon_p = \frac{expVL + a}{tau_L K_L (expVL + 1)} epsilon_m = \frac{b (expVL + a)}{tau_L (b expVL + a)} y_gate_inf = \frac{1}{1 + ⅇ^{\frac{V + 16.6577}{const5}}} + \frac{0.1}{1 + ⅇ^{\frac{(- V) + 40}{6}}} y_gate_tau = 20 + \frac{600}{1 + ⅇ^{\frac{V + 30}{9.6}}} C = 1 - O - I \frac{d}{d time} (O) = alpha_p C - (alpha_m O) \frac{d}{d time} (I) = epsilon_p Cass C - (epsilon_m I) \frac{d}{d time} (y_gate) = \frac{y_gate_inf - y_gate}{y_gate_tau} i_CaL = \{\begin{matrix} \frac{\frac{(- P_CaL) 2 Vss F}{Acap Cm} O y_gate FVRT_Ca}{1 - (ⅇ^{- FVRT_Ca})} (Cao ⅇ^{- FVRT_Ca} - Cass) if | FVRT_Ca | > 0.00001 \\ \frac{\frac{(- P_CaL) 2 Vss F}{Acap Cm} O y_gate 0.00001}{1 - (ⅇ^{- 0.00001})} (Cao ⅇ^{- 0.00001} - Cass) otherwise \end{matrix}

Component: calcium_pump_current

i_pCa = \frac{i_pCa_max {Cai}^{2}}{{Km_pCa}^{2} + {Cai}^{2}} J_pCa = \frac{(- i_pCa) Acap Cm}{2 Vmyo F}

Component: sodium_calcium_exchange_current

i_NaCa = \frac{\frac{\frac{k_NaCa 1}{{K_mNa}^{3} + {Nao}^{3}} 1}{K_mCa + Cao} 1}{1 + k_sat ⅇ^{\frac{(eta - 1) V F}{R T}}} (ⅇ^{\frac{eta V F}{R T}} {Nai}^{3} Cao - (ⅇ^{\frac{(eta - 1) V F}{R T}} {Nao}^{3} Cai)) J_ncx = \frac{i_NaCa Acap Cm}{Vmyo F}

Component: calcium_background_current

i_Cab = g_Cab (V - E_CaN) E_CaN = \frac{R T}{2 F} \ln (\frac{Cao}{Cai}) J_Cab = \frac{(- i_Cab) Acap Cm}{2 Vmyo F}

Component: sodium_concentration

\frac{d}{d time} (Nai) = \frac{(- (i_Na + i_Nab + 3 i_NaK + 3 i_NaCa)) Acap Cm}{Vmyo F}

Component: fast_sodium_current

i_Na = g_Na O_Na (V - E_Na) E_Na = \frac{R T}{F} \ln (\frac{0.9 Nao + 0.1 Ko}{0.9 Nai + 0.1 Ki}) C_Na3 = 1 - (O_Na + C_Na1 + C_Na2 + IF_Na + I1_Na + I2_Na + IC_Na2 + IC_Na3) \frac{d}{d time} (C_Na2) = alpha_Na11 C_Na3 + beta_Na12 C_Na1 + alpha_Na3 IC_Na2 - (beta_Na11 C_Na2 + alpha_Na12 C_Na2 + beta_Na3 C_Na2) \frac{d}{d time} (C_Na1) = alpha_Na12 C_Na2 + beta_Na13 O_Na + alpha_Na3 IF_Na - (beta_Na12 C_Na1 + alpha_Na13 C_Na1 + beta_Na3 C_Na1) \frac{d}{d time} (O_Na) = alpha_Na13 C_Na1 + beta_Na2 IF_Na - (beta_Na13 O_Na + alpha_Na2 O_Na) \frac{d}{d time} (IF_Na) = alpha_Na2 O_Na + beta_Na3 C_Na1 + beta_Na4 I1_Na + alpha_Na12 IC_Na2 - (beta_Na2 IF_Na + alpha_Na3 IF_Na + alpha_Na4 IF_Na + beta_Na12 IF_Na) \frac{d}{d time} (I1_Na) = alpha_Na4 IF_Na + beta_Na5 I2_Na - (beta_Na4 I1_Na + alpha_Na5 I1_Na) \frac{d}{d time} (I2_Na) = alpha_Na5 I1_Na - (beta_Na5 I2_Na) \frac{d}{d time} (IC_Na2) = alpha_Na11 IC_Na3 + beta_Na12 IF_Na + beta_Na3 IC_Na2 - (beta_Na11 IC_Na2 + alpha_Na12 IC_Na2 + alpha_Na3 IC_Na2) \frac{d}{d time} (IC_Na3) = beta_Na11 IC_Na2 + beta_Na3 C_Na3 - (alpha_Na11 IC_Na3 + alpha_Na3 IC_Na3) alpha_Na11 = \frac{3.802}{0.1027 ⅇ^{\frac{- (V + 2.5)}{17}} + 0.2 ⅇ^{\frac{- (V + 2.5)}{150}}} alpha_Na12 = \frac{3.802}{0.1027 ⅇ^{\frac{- (V + 2.5)}{15}} + 0.23 ⅇ^{\frac{- (V + 2.5)}{150}}} alpha_Na13 = \frac{3.802}{0.1027 ⅇ^{\frac{- (V + 2.5)}{12}} + 0.25 ⅇ^{\frac{- (V + 2.5)}{150}}} beta_Na11 = 0.1917 ⅇ^{\frac{- (V + 2.5)}{20.3}} beta_Na12 = 0.2 ⅇ^{\frac{- (V - 2.5)}{20.3}} beta_Na13 = 0.22 ⅇ^{\frac{- (V - 7.5)}{20.3}} alpha_Na3 = 7 e -7 ⅇ^{\frac{- (V + 7)}{7.7}} beta_Na3 = 0.0084 + 0.00002 (V + 7) alpha_Na2 = \frac{1}{0.188495 ⅇ^{\frac{- (V + 7)}{16.6}} + 0.393956} beta_Na2 = \frac{alpha_Na13 alpha_Na2 alpha_Na3}{beta_Na13 beta_Na3} alpha_Na4 = \frac{alpha_Na2}{1000} beta_Na4 = alpha_Na3 alpha_Na5 = \frac{alpha_Na2}{95000} beta_Na5 = \frac{alpha_Na3}{50}

Component: sodium_background_current

i_Nab = g_Nab (V - E_Na)

Component: potassium_concentration

\frac{d}{d time} (Ki) = \frac{(- (i_Stim + i_Kto_f + i_Kto_s + i_K1 + i_Ks + i_Kss + i_Kur + i_Kr - (2 i_NaK))) Acap Cm}{Vmyo F}

Component: fast_transient_outward_K_I

i_Kto_f = g_Kto_f {ato_f}^{3} ito_f (V - E_K) E_K = \frac{R T}{F} \ln (\frac{Ko}{Ki}) \frac{d}{d time} (ato_f) = alpha_a (1 - ato_f) - (beta_a ato_f) \frac{d}{d time} (ito_f) = alpha_i (1 - ito_f) - (beta_i ito_f) alpha_a = 0.18064 ⅇ^{0.03577 (V + 30)} beta_a = 0.3956 ⅇ^{(- 0.06237) (V + 30)} alpha_i = \frac{0.000152 ⅇ^{\frac{- (V + 13.5)}{7}}}{0.0067083 ⅇ^{\frac{- (V + 33.5)}{7}} + 1} beta_i = \frac{0.00095 ⅇ^{\frac{V + 33.5}{7}}}{0.051335 ⅇ^{\frac{V + 33.5}{7}} + 1}

Component: slow_transient_outward_K_I

i_Kto_s = g_Kto_s ato_s ito_s (V - E_K) \frac{d}{d time} (ato_s) = \frac{ass - ato_s}{tau_ta_s} \frac{d}{d time} (ito_s) = \frac{iss - ito_s}{tau_ti_s} ass = \frac{1}{1 + ⅇ^{\frac{- (V + 22.5)}{7.7}}} iss = \frac{1}{1 + ⅇ^{\frac{V + 45.2}{5.7}}} tau_ta_s = 0.493 ⅇ^{(- 0.0629) V} + 2.058 tau_ti_s = 270 + \frac{1050}{1 + ⅇ^{\frac{V + 45.2}{5.7}}}

Component: time_independent_K_I

i_K1 = \frac{\frac{g_K1 Ko}{Ko + 210} (V - E_K)}{1 + ⅇ^{0.0896 (V - E_K)}}

Component: slow_delayed_rectifier_K_I

i_Ks = g_Ks {nKs}^{2} (V - E_K) \frac{d}{d time} (nKs) = alpha_n (1 - nKs) - (beta_n nKs) alpha_n = \frac{0.00000481333 (V + 26.5)}{1 - (ⅇ^{(- 0.128) (V + 26.5)})} beta_n = 0.0000953333 ⅇ^{(- 0.038) (V + 26.5)}

Component: ultra_rapidly_activating_delayed_rectifier_K_I

i_Kur = g_Kur aur iur (V - E_K) \frac{d}{d time} (aur) = \frac{ass - aur}{tau_aur} \frac{d}{d time} (iur) = \frac{iss - iur}{tau_iur} tau_aur = 0.493 ⅇ^{(- 0.0629) V} + 2.058 tau_iur = 1200 - (\frac{170}{1 + ⅇ^{\frac{V + 45.2}{5.7}}})

Component: non_inactivating_steady_state_K_I

i_Kss = g_Kss aKss iKss (V - E_K) \frac{d}{d time} (aKss) = \frac{ass - aKss}{tau_Kss} \frac{d}{d time} (iKss) = 0 tau_Kss = 39.3 ⅇ^{(- 0.0862) V} + 13.17

Component: rapid_delayed_rectifier_K_I

i_Kr = g_Kr O_K (V - (\frac{R T}{F} \ln (\frac{0.98 Ko + 0.02 Nao}{0.98 Ki + 0.02 Nai}))) C_K0 = 1 - (C_K1 + C_K2 + O_K + I_K) \frac{d}{d time} (C_K2) = kf C_K1 + beta_a1 O_K - (kb C_K2 + alpha_a1 C_K2) \frac{d}{d time} (C_K1) = alpha_a0 C_K0 + kb C_K2 - (beta_a0 C_K1 + kf C_K1) \frac{d}{d time} (O_K) = alpha_a1 C_K2 + beta_i I_K - (beta_a1 O_K + alpha_i O_K) \frac{d}{d time} (I_K) = alpha_i O_K - (beta_i I_K) alpha_a0 = 0.022348 ⅇ^{0.01176 V} beta_a0 = 0.047002 ⅇ^{(- 0.0631) V} alpha_a1 = 0.013733 ⅇ^{0.038198 V} beta_a1 = 0.0000689 ⅇ^{(- 0.04178) V} alpha_i = 0.090821 ⅇ^{0.023391 (V + 5)} beta_i = 0.006497 ⅇ^{(- 0.03268) (V + 5)}

Component: sodium_potassium_pump_current

i_NaK = \frac{\frac{i_NaK_max f_NaK 1}{1 + {(\frac{Km_Nai}{Nai})}^{1.5}} Ko}{Ko + Km_Ko} f_NaK = \frac{1}{1 + 0.1245 ⅇ^{\frac{(- 0.1) V F}{R T}} + 0.0365 sigma ⅇ^{\frac{(- V) F}{R T}}} sigma = \frac{1}{7} (ⅇ^{\frac{Nao}{67300}} - 1)

Component: calcium_activated_chloride_current

i_ClCa = \frac{g_ClCa O_ClCa Cai}{Cai + Km_Cl} (V - E_Cl) O_ClCa = \frac{0.2}{1 + ⅇ^{\frac{- (V - 46.7)}{7.8}}}

Model Mathematics

Component: environment

Component: membrane

Component: calcium_concentration

Component: calcium_fluxes

Component: ryanodine_receptors

Component: L_type_calcium_current

Component: calcium_pump_current

Component: sodium_calcium_exchange_current

Component: calcium_background_current

Component: sodium_concentration

Component: fast_sodium_current

Component: sodium_background_current

Component: potassium_concentration

Component: fast_transient_outward_K_I

Component: slow_transient_outward_K_I

Component: time_independent_K_I

Component: slow_delayed_rectifier_K_I

Component: ultra_rapidly_activating_delayed_rectifier_K_I

Component: non_inactivating_steady_state_K_I

Component: rapid_delayed_rectifier_K_I

Component: sodium_potassium_pump_current

Component: calcium_activated_chloride_current

Component: transmembrane_currents

Component: intracellular_currents

Component: Ions_n_reversal_potentials