Model Mathematics

FCell = \frac{1.07 (3 dCell - 0.1)}{3 (1 + 0.7745 ⅇ^{\frac{- (3 dCell - 2.05)}{0.295}})} Cm = CmCentre + FCell (CmPeriphery - CmCentre) \frac{d}{d time} (V) = \frac{- 1}{Cm} (i_Na + i_Ca_L + i_Ca_T + i_to + i_sus + i_K_r + i_K_s + i_f_Na + i_f_K + i_b_Na + i_b_Ca + i_b_K + i_NaCa + i_p + i_Ca_p)

g_Na = g_Na_Centre + FCell (g_Na_Periphery - g_Na_Centre) i_Na = \frac{\frac{g_Na m^{3} h Na_o F^{2}}{R T} (ⅇ^{\frac{(V - E_Na) F}{R T}} - 1)}{ⅇ^{\frac{V F}{R T}} - 1} V

\frac{d}{d time} (m) = \frac{m_infinity - m}{tau_m} m_infinity = {(\frac{1}{1 + ⅇ^{\frac{- V}{5.46}}})}^{\frac{1}{3}} tau_m = \frac{0.0006247}{0.832 ⅇ^{(- 0.335) (V + 56.7)} + 0.627 ⅇ^{0.082 (V + 65.01)}} + 4 e -5

F_Na = \frac{0.0952 ⅇ^{(- 0.063) (V + 34.4)}}{1 + 1.66 ⅇ^{(- 0.225) (V + 63.7)}} + 0.0869 h = (1 - F_Na) h1 + F_Na h2 \frac{d}{d time} (h1) = \frac{h1_infinity - h1}{tau_h1} \frac{d}{d time} (h2) = \frac{h2_infinity - h2}{tau_h2} h1_infinity = \frac{1}{1 + ⅇ^{\frac{V + 66.1}{6.4}}} h2_infinity = h1_infinity tau_h1 = \frac{3.717 e -6 ⅇ^{(- 0.2815) (V + 17.11)}}{1 + 0.003732 ⅇ^{(- 0.3426) (V + 37.76)}} + 0.0005977 tau_h2 = \frac{3.186 e -8 ⅇ^{(- 0.6219) (V + 18.8)}}{1 + 7.189 e -5 ⅇ^{(- 0.6683) (V + 34.07)}} + 0.003556

g_Ca_L = g_Ca_L_Centre + FCell (g_Ca_L_Periphery - g_Ca_L_Centre) i_Ca_L = g_Ca_L (f_L d_L + \frac{0.006}{1 + ⅇ^{\frac{- (V + 14.1)}{6}}}) (V - E_Ca_L)

\frac{d}{d time} (d_L) = \frac{d_L_infinity - d_L}{tau_d_L} alpha_d_L = \frac{(- 28.38) (V + 35)}{ⅇ^{\frac{- (V + 35)}{2.5}} - 1} - (\frac{84.9 V}{ⅇ^{(- 0.208) V} - 1}) beta_d_L = \frac{11.42 (V - 5)}{ⅇ^{0.4 (V - 5)} - 1} tau_d_L = \frac{2}{alpha_d_L + beta_d_L} d_L_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 23.1)}{6}}}

\frac{d}{d time} (f_L) = \frac{f_L_infinity - f_L}{tau_f_L} alpha_f_L = \frac{3.12 (V + 28)}{ⅇ^{\frac{V + 28}{4}} - 1} beta_f_L = \frac{25}{1 + ⅇ^{\frac{- (V + 28)}{4}}} tau_f_L = \frac{1}{alpha_f_L + beta_f_L} f_L_infinity = \frac{1}{1 + ⅇ^{\frac{V + 45}{5}}}

g_Ca_T = g_Ca_T_Centre + FCell (g_Ca_T_Periphery - g_Ca_T_Centre) i_Ca_T = g_Ca_T d_T f_T (V - E_Ca_T)

\frac{d}{d time} (d_T) = \frac{d_T_infinity - d_T}{tau_d_T} alpha_d_T = 1068 ⅇ^{\frac{V + 26.3}{30}} beta_d_T = 1068 ⅇ^{\frac{- (V + 26.3)}{30}} tau_d_T = \frac{1}{alpha_d_T + beta_d_T} d_T_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 37)}{6.8}}}

\frac{d}{d time} (f_T) = \frac{f_T_infinity - f_T}{tau_f_T} alpha_f_T = 15.3 ⅇ^{\frac{- (V + 71.7)}{83.3}} beta_f_T = 15 ⅇ^{\frac{V + 71.7}{15.38}} tau_f_T = \frac{1}{alpha_f_T + beta_f_T} f_T_infinity = \frac{1}{1 + ⅇ^{\frac{V + 71}{9}}}

g_to = g_to_Centre + FCell (g_to_Periphery - g_to_Centre) g_sus = g_sus_Centre + FCell (g_sus_Periphery - g_sus_Centre) i_to = g_to q r (V - E_K) i_sus = g_sus r (V - E_K)

\frac{d}{d time} (q) = \frac{q_infinity - q}{tau_q} q_infinity = \frac{1}{1 + ⅇ^{\frac{V + 59.37}{13.1}}} tau_q = 0.0101 + \frac{0.06517}{0.57 ⅇ^{(- 0.08) (V + 49)}} + 2.4 e -5 ⅇ^{0.1 (V + 50.93)}

\frac{d}{d time} (r) = \frac{r_infinity - r}{tau_r} r_infinity = \frac{1}{1 + ⅇ^{\frac{- (V - 10.93)}{19.7}}} tau_r = 0.001 (2.98 + \frac{15.59}{1.037 ⅇ^{0.09 (V + 30.61)} + 0.369 ⅇ^{(- 0.12) (V + 23.84)}})

g_K_r = g_K_r_Centre + FCell (g_K_r_Periphery - g_K_r_Centre) i_K_r = g_K_r P_a P_i (V - E_K) P_a = 0.6 P_af + 0.4 P_as

\frac{d}{d time} (P_af) = \frac{P_af_infinity - P_af}{tau_P_af} P_af_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 14.2)}{10.6}}} tau_P_af = \frac{1}{37.2 ⅇ^{\frac{V - 9}{15.9}} + 0.96 ⅇ^{\frac{- (V - 9)}{22.5}}}

\frac{d}{d time} (P_as) = \frac{P_as_infinity - P_as}{tau_P_as} P_as_infinity = P_af_infinity tau_P_as = \frac{1}{4.2 ⅇ^{\frac{V - 9}{17}} + 0.15 ⅇ^{\frac{- (V - 9)}{21.6}}}

\frac{d}{d time} (P_i) = \frac{P_i_infinity - P_i}{tau_P_i} P_i_infinity = \frac{1}{1 + ⅇ^{\frac{V + 18.6}{10.1}}}

g_K_s = g_K_s_Centre + FCell (g_K_s_Periphery - g_K_s_Centre) i_K_s = g_K_s {xs}^{2} (V - E_K_s)

\frac{d}{d time} (xs) = alpha_xs (1 - xs) - (beta_xs xs) alpha_xs = \frac{14}{1 + ⅇ^{\frac{- (V - 40)}{9}}} beta_xs = 1 ⅇ^{\frac{- V}{45}}

g_f_Na = g_f_Na_Centre + FCell (g_f_Na_Periphery - g_f_Na_Centre) i_f_Na = g_f_Na y (V - E_Na) g_f_K = g_f_K_Centre + FCell (g_f_K_Periphery - g_f_K_Centre) i_f_K = g_f_K y (V - E_K)

\frac{d}{d time} (y) = alpha_y (1 - y) - (beta_y y) alpha_y = 1 ⅇ^{\frac{- (V + 78.91)}{26.62}} beta_y = 1 ⅇ^{\frac{V + 75.13}{21.25}}

g_b_Na = g_b_Na_Centre + FCell (g_b_Na_Periphery - g_b_Na_Centre) i_b_Na = g_b_Na (V - E_Na)

g_b_K = g_b_K_Centre + FCell (g_b_K_Periphery - g_b_K_Centre) i_b_K = g_b_K (V - E_K)

g_b_Ca = g_b_Ca_Centre + FCell (g_b_Ca_Periphery - g_b_Ca_Centre) i_b_Ca = g_b_Ca (V - E_Ca)

k_NaCa = k_NaCa_Centre + FCell (k_NaCa_Periphery - k_NaCa_Centre) i_NaCa = \frac{k_NaCa ({Na_i}^{3} Ca_o ⅇ^{0.03743 V gamma_NaCa} - ({Na_o}^{3} Ca_i ⅇ^{0.0374 V (gamma_NaCa - 1)}))}{1 + d_NaCa (Ca_i {Na_o}^{3} + Ca_o {Na_i}^{3})}

i_p_max = i_p_max_Centre + FCell (i_p_max_Periphery - i_p_max_Centre) i_p = \frac{i_p_max {(\frac{Na_i}{K_m_Na + Na_i})}^{3} {(\frac{K_o}{K_m_K + K_o})}^{2} 1.6}{1.5 + ⅇ^{\frac{- (V + 60)}{40}}}

i_Ca_p_max = i_Ca_p_max_Centre + FCell (i_Ca_p_max_Periphery - i_Ca_p_max_Centre) i_Ca_p = \frac{i_Ca_p_max Ca_i}{Ca_i + 0.0004}

E_Na = \frac{R T}{F} \ln (\frac{Na_o}{Na_i}) E_K = \frac{R T}{F} \ln (\frac{K_o}{K_i}) E_Ca = \frac{R T}{2 F} \ln (\frac{Ca_o}{Ca_i}) E_K_s = \frac{R T}{F} \ln (\frac{K_o + 0.12 Na_o}{K_i + 0.12 Na_i})