Model Mathematics

RTONF = \frac{R T}{F} i_Stim = \{\begin{matrix} stim_amplitude if time \geq stim_start \land time ≦ stim_end \land time - stim_start - (⌊ \frac{time - stim_start}{stim_period} ⌋ stim_period) ≦ stim_duration \\ 0 otherwise \end{matrix} \frac{d}{d time} (V) = \frac{- 1}{C_m} (i_Stim + i_b_K + i_K1 + i_to + i_b_Na + i_b_Ca + i_NaK + i_NaCa + i_Na + i_Ca_L)

E_mh = RTONF \ln (\frac{Na_o + 0.12 K_c}{Na_i + 0.12 K_i}) i_Na = g_Na m^{3} h (V - E_mh)

E0_m = V + 41 alpha_m = \{\begin{matrix} 2000 if | E0_m | < delta_m \\ \frac{200 E0_m}{1 - (ⅇ^{(- 0.1) E0_m})} otherwise \end{matrix} beta_m = 8000 ⅇ^{(- 0.056) (V + 66)} \frac{d}{d time} (m) = alpha_m (1 - m) - (beta_m m)

alpha_h = 20 ⅇ^{(- 0.125) (V + 75)} beta_h = \frac{2000}{1 + 320 ⅇ^{(- 0.1) (V + 75)}} \frac{d}{d time} (h) = alpha_h (1 - h) - (beta_h h)

E_K = RTONF \ln (\frac{K_c}{K_i}) i_to = g_to (g_to_s + s (1 - g_to_s)) r (V - E_K)

\frac{d}{d time} (r) = 333 (\frac{1}{1 + ⅇ^{\frac{- (V + 4)}{5}}} - r)

alpha_s = 0.033 ⅇ^{\frac{- V}{17}} beta_s = \frac{33}{1 + ⅇ^{(- 0.125) (V + 10)}} \frac{d}{d time} (s) = alpha_s (1 - s) - (beta_s s)

i_NaK = \frac{\frac{i_NaK_max K_c}{K_mK + K_c} Na_i}{K_mNa + Na_i}

E_Na = RTONF \ln (\frac{Na_o}{Na_i}) i_b_Na = g_b_Na (V - E_Na)

E_Ca = 0.5 RTONF \ln (\frac{Ca_o}{Ca_i}) i_b_Ca = g_b_Ca (V - E_Ca)

i_NaCa = \frac{k_NaCa (ⅇ^{\frac{gamma (n_NaCa - 2) V}{RTONF}} {Na_i}^{n_NaCa} Ca_o - (ⅇ^{\frac{(gamma - 1) (n_NaCa - 2) V}{RTONF}} {Na_o}^{n_NaCa} Ca_i))}{(1 + d_NaCa (Ca_i {Na_o}^{n_NaCa} + Ca_o {Na_i}^{n_NaCa})) (1 + \frac{Ca_i}{0.0069})}

i_b_K = g_b_K (V - E_K)

i_K1 = \frac{\frac{g_K1 K_c}{K_c + K_m_K1} (V - E_K)}{1 + ⅇ^{\frac{(V - E_K - 10) 2}{RTONF}}}

i_Ca_L = i_Ca_L_Ca + i_Ca_L_K + i_Ca_L_Na i_Ca_L_Ca = \frac{\frac{4 P_Ca_L d CaChon (V - 50)}{RTONF}}{1 - (ⅇ^{\frac{(- (V - 50)) 2}{RTONF}})} (Ca_i ⅇ^{\frac{100}{RTONF}} - (Ca_o ⅇ^{\frac{(- (V - 50)) 2}{RTONF}})) i_Ca_L_K = \frac{\frac{0.002 P_Ca_L d CaChon (V - 50)}{RTONF}}{1 - (ⅇ^{\frac{- (V - 50)}{RTONF}})} (K_i ⅇ^{\frac{50}{RTONF}} - (K_c ⅇ^{\frac{- (V - 50)}{RTONF}})) i_Ca_L_Na = \frac{\frac{0.01 P_Ca_L d CaChon (V - 50)}{RTONF}}{1 - (ⅇ^{\frac{- (V - 50)}{RTONF}})} (Na_i ⅇ^{\frac{50}{RTONF}} - (Na_o ⅇ^{\frac{- (V - 50)}{RTONF}}))

E0_d = V + 19 alpha_d = \{\begin{matrix} 120 if | E0_d | < 0.0001 \\ \frac{30 E0_d}{1 - (ⅇ^{\frac{- E0_d}{4}})} otherwise \end{matrix} beta_d = \{\begin{matrix} 120 if | E0_d | < 0.0001 \\ \frac{12 E0_d}{ⅇ^{\frac{E0_d}{10}} - 1} otherwise \end{matrix} \frac{d}{d time} (d) = alpha_d (1 - d) - (beta_d d)

E0_f = V + 34 alpha_f_Ca = \{\begin{matrix} 25 if | E0_f | < 0.0001 \\ \frac{6.25 E0_f}{ⅇ^{\frac{E0_f}{4}} - 1} otherwise \end{matrix} beta_f_Ca = \frac{12}{1 + ⅇ^{\frac{- E0_f}{4}}} \frac{d}{d time} (f_Ca) = (120 (1 - f_Ca) CaChoff + (1 - f_Ca) (1 - CaChoff)) beta_f_Ca - (alpha_f_Ca f_Ca) CaChoff = \frac{Ca_i}{0.001 + Ca_i} CaChon = (1 - f_Ca) (1 - CaChoff)

K_1 = \frac{K_cyca K_xcs}{K_srca} K_2 = Ca_i + Ca_up K_1 + K_cyca K_xcs + K_cyca i_up = \frac{Ca_i}{K_2} alpha_up - (\frac{Ca_up K_1}{K_2} beta_up)

PrecFrac = 1 - ActFrac - ProdFrac VoltDep = ⅇ^{0.08 (V - 40)} RegBindSite = {(\frac{Ca_i}{Ca_i + 0.0005})}^{2} ActRate = 600 VoltDep + 500 RegBindSite InactRate = 60 + 500 RegBindSite \frac{d}{d time} (ActFrac) = PrecFrac ActRate - (ActFrac InactRate) \frac{d}{d time} (ProdFrac) = ActFrac InactRate - (ProdFracRate ProdFrac) i_rel = ({(\frac{ActFrac}{ActFrac + 0.25})}^{2} K_m_rel + K_leak_rate) Ca_rel

i_trans = (Ca_up - Ca_rel) alpha_tr

\frac{d}{d time} (Na_i) = \frac{- 1}{1 V_i F} (i_Na + i_b_Na + 3 i_NaK + 3 i_NaCa + i_Ca_L_Na)

\frac{d}{d time} (K_i) = \frac{- 1}{1 V_i F} (i_K1 + i_Ca_L_K + i_to + i_b_K - (2 i_NaK))

V_Cell = 3.141592654 {radius}^{2} length V_i_ratio = 1 - V_e_ratio - V_up_ratio - V_rel_ratio V_i = V_Cell V_i_ratio \frac{d}{d time} (Ca_i) = \frac{- 1}{2 1 V_i F} (i_Ca_L_Ca + i_b_Ca - (2 i_NaCa)) + \frac{i_rel V_rel_ratio}{V_i_ratio} - (\frac{d}{d time} (Ca_Calmod)) - (\frac{d}{d time} (Ca_Trop)) - i_up \frac{d}{d time} (Ca_up) = \frac{V_i_ratio}{V_up_ratio} i_up - i_trans \frac{d}{d time} (Ca_rel) = \frac{V_up_ratio}{V_rel_ratio} i_trans - i_rel \frac{d}{d time} (Ca_Calmod) = alpha_Calmod Ca_i (Calmod - Ca_Calmod) - (beta_Calmod Ca_Calmod) \frac{d}{d time} (Ca_Trop) = alpha_Trop Ca_i (Trop - Ca_Trop) - (beta_Trop Ca_Trop)