Model Mathematics

Component based math viewer is available

Component: main

Component: TT04

phi_Cai = (T_o - 1.2) K_Cai pC_50 = pC_50_ref (1.0 + beta_2 (lambda - 1.0)) C_50 = 1.0 e -3 10^{6 - pC_50} n = n_ref (1.0 + beta_1 (lambda - 1.0)) T_o = T_ref (1.0 + beta_0 (lambda - 1.0)) z \frac{d}{d time} (Ca_b) = rho_0 Ca_i (Ca_b_max - Ca_b) - (rho_1 Ca_b) \frac{d}{d time} (z) = alpha_0 ({(\frac{Ca_b}{C_50})}^{n} (1.0 - z) - z)

Component: cvs

u_lv_max = phi E_lv_max (q_lv - q_lv_u) + (1 - phi) u_lv_0 (ⅇ^{k_lv_E q_lv} - 1) u_rv_max = phi E_rv_max (q_rv - q_rv_u) + (1 - phi) u_rv_0 (ⅇ^{k_rv_E q_rv} - 1) \frac{d}{d t} (q_ra) = v_sv - v_TV \frac{d}{d t} (q_rv) = v_TV - v_PV \frac{d}{d t} (q_pa) = v_PV - v_pa \frac{d}{d t} (q_pc) = v_pa - v_pv \frac{d}{d t} (q_pv) = v_pv - v_la \frac{d}{d t} (q_la) = v_la - v_MV \frac{d}{d t} (q_lv) = v_MV - v_AV \frac{d}{d t} (q_sa) = v_AV - v_sa \frac{d}{d t} (q_sc) = v_sa - v_sc q_sv = q_tot - q_tot_u - (q_ra + q_rv + q_pa + q_pc + q_pv + q_la + q_lv + q_sa + q_sc) u_ra = \frac{q_ra}{C_ra} u_rv = u_rv_max - (R_PV v_PV) u_pa = \frac{q_pa}{C_pa} u_pc = \frac{q_pc}{C_pc} u_pv = \frac{q_pv}{C_pv} u_la = \frac{q_la}{C_la} u_lv = u_lv_max - (R_AV v_AV) u_sa = \frac{q_sa}{C_sa} u_sc = \frac{q_sc}{C_sc} u_sv = \frac{q_sv}{C_sv} du_sa = \frac{v_AV - v_sa}{C_sa} v_TV = \{\begin{matrix} 0 if u_ra < u_rv \\ \frac{u_ra - u_rv}{R_TV} otherwise \end{matrix} v_PV = \{\begin{matrix} 0 if u_rv_max < u_pa \\ \frac{u_rv_max - u_pa}{R_PV} otherwise \end{matrix} R_PV = k_PV_R u_rv_max \frac{d}{d t} (v_pa) = \frac{u_pa - u_pc - (R_pa v_pa)}{L_pa} v_pv = \frac{u_pc - u_pv}{R_pv} v_la = \frac{u_pv - u_la}{R_la} v_MV = \{\begin{matrix} 0 if u_la < u_lv \\ \frac{u_la - u_lv}{R_MV} otherwise \end{matrix} v_AV = \{\begin{matrix} 0 if u_lv_max < u_sa \\ \frac{u_lv_max - u_sa}{R_AV} otherwise \end{matrix} R_AV = k_AV_R u_lv_max \frac{d}{d t} (v_sa) = \frac{u_sa - u_sc - (R_sa v_sa)}{L_sa} v_sc = \frac{u_sc - u_sv}{R_sc} v_sv = \frac{u_sv - u_ra}{R_sv}

Component: environment

Component: membrane

stim_period = \{\begin{matrix} rest_period if time < (⌊ \frac{rest_duration}{rest_period} ⌋ + 1) rest_period \\ act_period if time \geq (⌊ \frac{rest_duration}{rest_period} ⌋ + 1) rest_period \land time < ⌊ \frac{rest_duration + act_duration}{act_period} ⌋ act_period \\ rest_period otherwise \end{matrix} i_Stim = \{\begin{matrix} stim_amplitude if time - (⌊ \frac{time}{stim_period} ⌋ stim_period) \geq stim_start \land time - (⌊ \frac{time}{stim_period} ⌋ stim_period) ≦ stim_start + stim_duration \\ 0 otherwise \end{matrix} \frac{d}{d time} (V) = - (i_K1 + i_to + i_Kr + i_Ks + i_CaL + i_NaK + i_Na + i_b_Na + i_NaCa + i_b_Ca + i_p_K + i_p_Ca + i_Stim)

Component: activation

u = \frac{t}{T_0} - (⌊ \frac{t}{T_0} ⌋) T = t - (T_0 ⌊ \frac{t}{T_0} ⌋) T_sys = \{\begin{matrix} T_sys_0 - (\frac{k_sys}{T}) if T > 0 \\ T_sys_0 - (\frac{k_sys}{T_0}) otherwise \end{matrix} phi = \{\begin{matrix} {(\sin (\frac{π u T}{T_sys_0}))}^{2} if u \geq 0 \land u ≦ \frac{T_sys_0}{T} \\ 0 otherwise \end{matrix}

Component: fast_sodium_current

i_Na = g_Na m^{3} h j (V - E_Na)

Component: L_type_Ca_current

i_CaL = \frac{\frac{g_CaL d f fCa 4 V F^{2}}{R T} (Ca_i ⅇ^{\frac{2 V F}{R T}} - (0.341 Ca_o))}{ⅇ^{\frac{2 V F}{R T}} - 1}

Component: rapid_time_dependent_potassium_current

i_Kr = g_Kr \sqrt{\frac{K_o}{5.4}} Xr1 Xr2 (V - E_K)

Component: slow_time_dependent_potassium_current

i_Ks = g_Ks {Xs}^{2} (V - E_Ks)

Component: sodium_background_current

i_b_Na = g_bna (V - E_Na)

Component: calcium_background_current

i_b_Ca = g_bca (V - E_Ca)

Component: transient_outward_current

i_to = g_to r s (V - E_K)

Component: sodium_potassium_pump_current

i_NaK = \frac{\frac{\frac{P_NaK K_o}{K_o + K_mk} Na_i}{Na_i + K_mNa}}{1 + 0.1245 ⅇ^{\frac{(- 0.1) V F}{R T}} + 0.0353 ⅇ^{\frac{(- V) F}{R T}}}

Component: sodium_calcium_exchanger_current

i_NaCa = \frac{K_NaCa (ⅇ^{\frac{gamma V F}{R T}} {Na_i}^{3} Ca_o - (ⅇ^{\frac{(gamma - 1) V F}{R T}} {Na_o}^{3} Ca_i alpha))}{({Km_Nai}^{3} + {Na_o}^{3}) (Km_Ca + Ca_o) (1 + K_sat ⅇ^{\frac{(gamma - 1) V F}{R T}})}

Component: calcium_pump_current

i_p_Ca = \frac{g_pCa Ca_i}{Ca_i + K_pCa}

Component: potassium_pump_current

i_p_K = \frac{g_pK (V - E_K)}{1 + ⅇ^{\frac{25 - V}{5.98}}}

Component: calcium_dynamics

i_rel = (\frac{a_rel {Ca_SR}^{2}}{{b_rel}^{2} + {Ca_SR}^{2}} + c_rel) d g i_up = \frac{Vmax_up}{1 + \frac{{K_up}^{2}}{{Ca_i}^{2}}} i_leak = V_leak (Ca_SR - Ca_i) g_inf = \{\begin{matrix} \frac{1}{1 + {(\frac{Ca_i}{0.00035})}^{6}} if Ca_i < 0.00035 \\ \frac{1}{1 + {(\frac{Ca_i}{0.00035})}^{16}} otherwise \end{matrix} d_g = \frac{g_inf - g}{tau_g} \frac{d}{d time} (g) = \{\begin{matrix} 0 if g_inf > g \land V > - 60 \\ d_g otherwise \end{matrix} ddt_Ca_i_total = (\frac{- (i_CaL + i_b_Ca + i_p_Ca - (2 i_NaCa))}{2 V_c F} Cm + i_leak - i_up) + i_rel ddt_Ca_sr_total = \frac{V_c}{V_sr} (i_up - (i_rel + i_leak)) f_JCa_i_free = \frac{1}{1 + \frac{Buf_c K_buf_c}{{(Ca_i + K_buf_c)}^{2}}} f_JCa_sr_free = \frac{1}{1 + \frac{Buf_sr K_buf_sr}{{(Ca_SR + K_buf_sr)}^{2}}} \frac{d}{d time} (Ca_i) = ddt_Ca_i_total f_JCa_i_free \frac{d}{d time} (Ca_SR) = ddt_Ca_sr_total f_JCa_sr_free

Component: sodium_dynamics

\frac{d}{d time} (Na_i) = \frac{(- (i_Na + i_b_Na + 3 i_NaK + 3 i_NaCa)) Cm}{V_c F}

Component: potassium_dynamics

\frac{d}{d time} (K_i) = \frac{(- (i_K1 + i_to + i_Kr + i_Ks + i_p_K + i_Stim - (2 i_NaK))) Cm}{V_c F}

Component: reversal_potentials

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_Ks = \frac{R T}{F} \ln (\frac{K_o + P_kna Na_o}{K_i + P_kna Na_i}) E_Ca = \frac{0.5 R T}{F} \ln (\frac{Ca_o}{Ca_i})

Component: inward_rectifier_potassium_current

alpha_K1 = \frac{0.1}{1 + ⅇ^{0.06 (V - E_K - 200)}} beta_K1 = \frac{3 ⅇ^{0.0002 ((V - E_K) + 100)} + 1 ⅇ^{0.1 (V - E_K - 10)}}{1 + ⅇ^{(- 0.5) (V - E_K)}} xK1_inf = \frac{alpha_K1}{alpha_K1 + beta_K1} i_K1 = g_K1 xK1_inf \sqrt{\frac{K_o}{5.4}} (V - E_K)

Component: fast_sodium_current_m_gate

m_inf = \frac{1}{{(1 + ⅇ^{\frac{- 56.86 - V}{9.03}})}^{2}} alpha_m = \frac{1}{1 + ⅇ^{\frac{- 60 - V}{5}}} beta_m = \frac{0.1}{1 + ⅇ^{\frac{V + 35}{5}}} + \frac{0.1}{1 + ⅇ^{\frac{V - 50}{200}}} tau_m = 1 alpha_m beta_m \frac{d}{d time} (m) = \frac{m_inf - m}{tau_m}

Component: fast_sodium_current_h_gate

h_inf = \frac{1}{{(1 + ⅇ^{\frac{V + 71.55}{7.43}})}^{2}} alpha_h = \{\begin{matrix} 0.057 ⅇ^{\frac{- (V + 80)}{6.8}} if V < - 40 \\ 0 otherwise \end{matrix} beta_h = \{\begin{matrix} 2.7 ⅇ^{0.079 V} + 310000 ⅇ^{0.3485 V} if V < - 40 \\ \frac{0.77}{0.13 (1 + ⅇ^{\frac{V + 10.66}{- 11.1}})} otherwise \end{matrix} tau_h = \frac{1}{alpha_h + beta_h} \frac{d}{d time} (h) = \frac{h_inf - h}{tau_h}

Component: fast_sodium_current_j_gate

j_inf = \frac{1}{{(1 + ⅇ^{\frac{V + 71.55}{7.43}})}^{2}} alpha_j = \{\begin{matrix} \frac{\frac{((- 25428) ⅇ^{0.2444 V} - (6.948 e -6 ⅇ^{(- 0.04391) V})) (V + 37.78)}{1}}{1 + ⅇ^{0.311 (V + 79.23)}} if V < - 40 \\ 0 otherwise \end{matrix} beta_j = \{\begin{matrix} \frac{0.02424 ⅇ^{(- 0.01052) V}}{1 + ⅇ^{(- 0.1378) (V + 40.14)}} if V < - 40 \\ \frac{0.6 ⅇ^{0.057 V}}{1 + ⅇ^{(- 0.1) (V + 32)}} otherwise \end{matrix} tau_j = \frac{1}{alpha_j + beta_j} \frac{d}{d time} (j) = \frac{j_inf - j}{tau_j}

Component: L_type_Ca_current_d_gate

d_inf = \frac{1}{1 + ⅇ^{\frac{- 5 - V}{7.5}}} alpha_d = \frac{1.4}{1 + ⅇ^{\frac{- 35 - V}{13}}} + 0.25 beta_d = \frac{1.4}{1 + ⅇ^{\frac{V + 5}{5}}} gamma_d = \frac{1}{1 + ⅇ^{\frac{50 - V}{20}}} tau_d = 1 alpha_d beta_d + gamma_d \frac{d}{d time} (d) = \frac{d_inf - d}{tau_d}

Component: L_type_Ca_current_f_gate

f_inf = \frac{1}{1 + ⅇ^{\frac{V + 20}{7}}} tau_f = 1125 ⅇ^{\frac{- ({(V + 27)}^{2})}{240}} + 80 + \frac{165}{1 + ⅇ^{\frac{25 - V}{10}}} \frac{d}{d time} (f) = \frac{f_inf - f}{tau_f}

Component: L_type_Ca_current_fCa_gate

alpha_fCa = \frac{1}{1 + {(\frac{Ca_i}{0.000325})}^{8}} beta_fCa = \frac{0.1}{1 + ⅇ^{\frac{Ca_i - 0.0005}{0.0001}}} gama_fCa = \frac{0.2}{1 + ⅇ^{\frac{Ca_i - 0.00075}{0.0008}}} fCa_inf = \frac{alpha_fCa + beta_fCa + gama_fCa + 0.23}{1.46} tau_fCa = 2 d_fCa = \frac{fCa_inf - fCa}{tau_fCa} \frac{d}{d time} (fCa) = \{\begin{matrix} 0 if fCa_inf > fCa \land V > - 60 \\ d_fCa otherwise \end{matrix}

Component: rapid_time_dependent_potassium_current_Xr1_gate

xr1_inf = \frac{1}{1 + ⅇ^{\frac{- 26 - V}{7}}} alpha_xr1 = \frac{450}{1 + ⅇ^{\frac{- 45 - V}{10}}} beta_xr1 = \frac{6}{1 + ⅇ^{\frac{V + 30}{11.5}}} tau_xr1 = 1 alpha_xr1 beta_xr1 \frac{d}{d time} (Xr1) = \frac{xr1_inf - Xr1}{tau_xr1}

Component: rapid_time_dependent_potassium_current_Xr2_gate

xr2_inf = \frac{1}{1 + ⅇ^{\frac{V + 88}{24}}} alpha_xr2 = \frac{3}{1 + ⅇ^{\frac{- 60 - V}{20}}} beta_xr2 = \frac{1.12}{1 + ⅇ^{\frac{V - 60}{20}}} tau_xr2 = 1 alpha_xr2 beta_xr2 \frac{d}{d time} (Xr2) = \frac{xr2_inf - Xr2}{tau_xr2}

Component: slow_time_dependent_potassium_current_Xs_gate

xs_inf = \frac{1}{1 + ⅇ^{\frac{- 5 - V}{14}}} alpha_xs = \frac{1100}{\sqrt{1 + ⅇ^{\frac{- 10 - V}{6}}}} beta_xs = \frac{1}{1 + ⅇ^{\frac{V - 60}{20}}} tau_xs = 1 alpha_xs beta_xs \frac{d}{d time} (Xs) = \frac{xs_inf - Xs}{tau_xs}

Component: transient_outward_current_s_gate

s_inf = \frac{1}{1 + ⅇ^{\frac{V + 20}{5}}} tau_s = 85 ⅇ^{\frac{- ({(V + 45)}^{2})}{320}} + \frac{5}{1 + ⅇ^{\frac{V - 20}{5}}} + 3 \frac{d}{d time} (s) = \frac{s_inf - s}{tau_s}

Component: transient_outward_current_r_gate

r_inf = \frac{1}{1 + ⅇ^{\frac{20 - V}{6}}} tau_r = 9.5 ⅇ^{\frac{- ({(V + 40)}^{2})}{1800}} + 0.8 \frac{d}{d time} (r) = \frac{r_inf - r}{tau_r}