Model Mathematics

Component based math viewer is available

Component: environment

Component: Membrane

i_tot = i_f + i_Kr + i_Ks + i_to + i_NaK + i_NaCa + i_Na + i_CaL + i_CaT + i_KACh \frac{d}{d time} (V_ode) = \frac{- i_tot}{C} RTONF = \frac{R T}{F} V = \{\begin{matrix} V_clamp if clamp_mode \geq 1 \\ V_ode otherwise \end{matrix}

Component: Voltage_clamp

V_clamp = \{\begin{matrix} V_test if time > t_holding \land time < t_holding + t_test \\ V_holding otherwise \end{matrix}

Component: Rate_modulation_experiments

Component: Ionic_values

E_Na = RTONF \ln (\frac{Nao}{Nai}) E_K = RTONF \ln (\frac{Ko}{Ki}) E_Ca = 0.5 RTONF \ln (\frac{Cao}{Ca_sub})

Component: Nai_concentration

\frac{d}{d time} (Nai_) = \frac{(- 1) (i_Na + i_fNa + i_siNa + 3 i_NaK + 3 i_NaCa)}{1 (V_i + V_sub) F} Nai = \{\begin{matrix} 7.5 if BAPTA_10_mM > 0 \\ Nai_otherwise \end{matrix}

Component: i_f

g_f_Na = \{\begin{matrix} 0.03 (1 - 0.66) if Iva_3_uM \geq 1 \\ 0.03 otherwise \end{matrix} g_f_K = \{\begin{matrix} 0.03 (1 - 0.66) if Iva_3_uM \geq 1 \\ 0.03 otherwise \end{matrix} ICs_on_Icontrol = \{\begin{matrix} \frac{\frac{10.6015}{5}}{\frac{10.6015}{5} + ⅇ^{\frac{(- 0.71) V}{25}}} if Cs_5_mM \geq 1 \\ 1 otherwise \end{matrix} i_fNa = \frac{y^{2} Ko}{Ko + Km_f} g_f_Na (V - E_Na) ICs_on_Icontrol i_fK = \frac{y^{2} Ko}{Ko + Km_f} g_f_K (V - E_K) ICs_on_Icontrol i_f = i_fNa + i_fK

Component: i_f_y_gate

ACh_shift = \{\begin{matrix} - 1 - (\frac{9.898 {(1 e -6 ACh)}^{0.618}}{{(1 e -6 ACh)}^{0.618} + 0.00122423}) if ACh > 0 \\ 0 otherwise \end{matrix} Iso_shift = \{\begin{matrix} 7.5 if Iso_1_uM > 0 \\ 0 otherwise \end{matrix} ANS_shift_i_f_y_gate = \{\begin{matrix} \frac{Ki_f {cAMP}^{ni_f}}{{K_05i_f}^{ni_f} + {cAMP}^{ni_f}} - 18.76 if ANS > 0 \\ 0 otherwise \end{matrix} tau_y = \frac{0.7166529}{0.0708 ⅇ^{\frac{- (V + 5 - ACh_shift - Iso_shift - ANS_shift_i_f_y_gate)}{20.2791}} + 10.6 ⅇ^{\frac{V - ACh_shift - Iso_shift - ANS_shift_i_f_y_gate}{18}}} y_infinity = \frac{1}{1 + ⅇ^{\frac{V + 52.5 - ACh_shift - Iso_shift - ANS_shift_i_f_y_gate}{9}}} \frac{d}{d time} (y) = \frac{y_infinity - y}{tau_y}

Component: i_NaK

Iso_increase = \{\begin{matrix} 1.2 if Iso_1_uM > 0 \\ 1 otherwise \end{matrix} ANS_cond_i_NaK = \{\begin{matrix} (- 0.2152) + \frac{K_i_NaK {PKA}^{10.0808}}{{K_05i_NaK}^{10.0808} + {PKA}^{10.0808}} if ANS > 0 \\ 0 otherwise \end{matrix} i_NaK = Iso_increase (1 + ANS_cond_i_NaK) i_NaK_max {(1 + {(\frac{Km_Kp}{Ko})}^{1.2})}^{- 1} {(1 + {(\frac{Km_Nap}{Nai})}^{1.3})}^{- 1} {(1 + ⅇ^{\frac{- ((V - E_Na) + 110)}{20}})}^{- 1}

Component: i_NaCa

i_NaCa = \frac{K_NaCa (x2 k21 - (x1 k12))}{x1 + x2 + x3 + x4} x1 = k41 k34 (k23 + k21) + k21 k32 (k43 + k41) x2 = k32 k43 (k14 + k12) + k41 k12 (k34 + k32) x3 = k14 k43 (k23 + k21) + k12 k23 (k43 + k41) x4 = k23 k34 (k14 + k12) + k14 k21 (k34 + k32) k43 = \frac{Nai}{K3ni + Nai} k12 = \frac{\frac{Ca_sub}{Kci} ⅇ^{\frac{(- Qci) V}{RTONF}}}{di} k14 = \frac{\frac{\frac{Nai}{K1ni} Nai}{K2ni} (1 + \frac{Nai}{K3ni}) ⅇ^{\frac{Qn V}{2 RTONF}}}{di} k41 = ⅇ^{\frac{(- Qn) V}{2 RTONF}} di = 1 + \frac{Ca_sub}{Kci} (1 + ⅇ^{\frac{(- Qci) V}{RTONF}} + \frac{Nai}{Kcni}) + \frac{Nai}{K1ni} (1 + \frac{Nai}{K2ni} (1 + \frac{Nai}{K3ni})) k34 = \frac{Nao}{K3no + Nao} k21 = \frac{\frac{Cao}{Kco} ⅇ^{\frac{Qco V}{RTONF}}}{do} k23 = \frac{\frac{\frac{Nao}{K1no} Nao}{K2no} (1 + \frac{Nao}{K3no}) ⅇ^{\frac{(- Qn) V}{2 RTONF}}}{do} k32 = ⅇ^{\frac{Qn V}{2 RTONF}} do = 1 + \frac{Cao}{Kco} (1 + ⅇ^{\frac{Qco V}{RTONF}}) + \frac{Nao}{K1no} (1 + \frac{Nao}{K2no} (1 + \frac{Nao}{K3no}))

Component: i_Na

E_mh = RTONF \ln (\frac{Nao + 0.12 Ko}{Nai + 0.12 Ki}) i_Na = g_Na m^{3} h (V - E_mh)

Component: i_Na_m_gate

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)

Component: i_Na_h_gate

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

Component: i_CaL

Iso_increase = \{\begin{matrix} 1.23 if Iso_1_uM > 0 \\ 1 otherwise \end{matrix} ANS_cond_i_CaL = \{\begin{matrix} (- 0.2152) + \frac{K_i_CaL1 {PKA}^{10.0808}}{{K_05i_CaL1}^{10.0808} + {PKA}^{10.0808}} if ANS > 0 \\ 0 otherwise \end{matrix} i_siCa = \frac{2 P_CaL (1 - ACh_block) Iso_increase (1 + ANS_cond_i_CaL) (V - 0)}{RTONF (1 - (ⅇ^{\frac{(- 1) (V - 0) 2}{RTONF}}))} (Ca_sub - (Cao ⅇ^{\frac{(- 2) (V - 0)}{RTONF}})) dL fL fCa i_siK = \frac{0.000365 P_CaL (1 - ACh_block) Iso_increase (1 + ANS_cond_i_CaL) (V - 0)}{RTONF (1 - (ⅇ^{\frac{(- 1) (V - 0)}{RTONF}}))} (Ki - (Ko ⅇ^{\frac{(- 1) (V - 0)}{RTONF}})) dL fL fCa i_siNa = \frac{0.0000185 P_CaL (1 - ACh_block) Iso_increase (1 + ANS_cond_i_CaL) (V - 0)}{RTONF (1 - (ⅇ^{\frac{(- 1) (V - 0)}{RTONF}}))} (Nai - (Nao ⅇ^{\frac{(- 1) (V - 0)}{RTONF}})) dL fL fCa ACh_block = \frac{0.31 1 e -6 ACh}{1 e -6 ACh + 0.00009} i_CaL = i_siCa + i_siK + i_siNa

Component: i_CaL_dL_gate

Iso_shift = \{\begin{matrix} - 8 if Iso_1_uM > 0 \\ 0 otherwise \end{matrix} ANS_shift_i_CaL_dL_gate = \{\begin{matrix} (- K_ANS_shift_red_ICaL_dL_gate) (\frac{K_i_CaL2 {PKA}^{9.281}}{{K_05i_CaL2}^{9.281} + {PKA}^{9.281}} - 18.76) if ANS > 0 \\ 0 otherwise \end{matrix} K_ANS_shift_red_ICaL_dL_gate = \{\begin{matrix} 1 - ((0.741 - PKA) 32.258) if PKA ≦ 0.741 \land PKA \geq 0.71 \\ 0 if PKA < 0.71 \\ 1 if PKA > 0.741 \end{matrix} Iso_slope = \{\begin{matrix} 0.69 if Iso_1_uM > 0 \\ 1 otherwise \end{matrix} ANS_slope_i_CaL = \{\begin{matrix} (- 1.199) PKA + 1.8526 if ANS > 0 \\ 1 otherwise \end{matrix} dL_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 20.3 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{Iso_slope ANS_slope_i_CaL 4.2}}} tau_dL = \frac{0.001}{alpha_dL + beta_dL} alpha_dL = \frac{(- 0.02839) (adVm + 41.8 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{ⅇ^{\frac{- (adVm + 41.8 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{2.5}} - 1} - (\frac{0.0849 (adVm + 6.8 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{ⅇ^{\frac{- (adVm + 6.8 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{4.8}} - 1}) adVm = \{\begin{matrix} - 41.80001 if V = - 41.8 \\ 0 if V = 0 \\ - 6.80001 if V = - 6.8 \\ V otherwise \end{matrix} beta_dL = \frac{0.01143 (bdVm + 1.8 - Iso_shift - ANS_shift_i_CaL_dL_gate)}{ⅇ^{\frac{bdVm + 1.8 - Iso_shift - ANS_shift_i_CaL_dL_gate}{2.5}} - 1} bdVm = \{\begin{matrix} - 1.80001 if V = - 1.8 \\ V otherwise \end{matrix} \frac{d}{d time} (dL) = \frac{dL_infinity - dL}{tau_dL}

Component: i_CaL_fL_gate

fL_infinity = \frac{1}{1 + ⅇ^{\frac{V + 37.4}{5.3}}} tau_fL = 0.001 (44.3 + 230 ⅇ^{- ({(\frac{V + 36}{10})}^{2})}) \frac{d}{d time} (fL) = \frac{fL_infinity - fL}{tau_fL}

Component: i_CaL_fCa_gate

fCa_infinity = \frac{Km_fCa}{Km_fCa + Ca_sub} tau_fCa = \frac{0.001 fCa_infinity}{alpha_fCa} \frac{d}{d time} (fCa) = \frac{fCa_infinity - fCa}{tau_fCa}

Component: i_CaT

i_CaT = \frac{2 P_CaT V}{RTONF (1 - (ⅇ^{\frac{(- 1) V 2}{RTONF}}))} (Ca_sub - (Cao ⅇ^{\frac{(- 2) V}{RTONF}})) dT fT

Component: i_CaT_dT_gate

dT_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 38.3)}{5.5}}} tau_dT = \frac{0.001}{1.068 ⅇ^{\frac{V + 38.3}{30}} + 1.068 ⅇ^{\frac{- (V + 38.3)}{30}}} \frac{d}{d time} (dT) = \frac{dT_infinity - dT}{tau_dT}

Component: i_CaT_fT_gate

fT_infinity = \frac{1}{1 + ⅇ^{\frac{V + 58.7}{3.8}}} tau_fT = \frac{1}{16.67 ⅇ^{\frac{- (V + 75)}{83.3}} + 16.67 ⅇ^{\frac{V + 75}{15.38}}} \frac{d}{d time} (fT) = \frac{fT_infinity - fT}{tau_fT}

Component: Ca_SR_release

koCa = \{\begin{matrix} koCa_max ((RyR_min - (\frac{RyR_max {PKA}^{nRyR}}{{K_05RyR}^{nRyR} + {PKA}^{nRyR}})) + 1) if ANS > 0 \\ 10000 otherwise \end{matrix} j_SRCarel = ks O (Ca_jsr - Ca_sub) kCaSR = MaxSR - (\frac{MaxSR - MinSR}{1 + {(\frac{EC50_SR}{Ca_jsr})}^{HSR}}) koSRCa = \frac{koCa}{kCaSR} kiSRCa = kiCa kCaSR \frac{d}{d time} (R) = kim RI - (kiSRCa Ca_sub R) - (koSRCa {Ca_sub}^{2} R - (kom O)) \frac{d}{d time} (O) = koSRCa {Ca_sub}^{2} R - (kom O) - (kiSRCa Ca_sub O - (kim I)) \frac{d}{d time} (I) = kiSRCa Ca_sub O - (kim I) - (kom I - (koSRCa {Ca_sub}^{2} RI)) \frac{d}{d time} (RI) = kom I - (koSRCa {Ca_sub}^{2} RI) - (kim RI - (kiSRCa Ca_sub R))

Component: Ca_intracellular_fluxes

b_up = \{\begin{matrix} - 0.25 if Iso_1_uM > 0 \\ \frac{0.7 1 e -6 ACh}{1 e -6 ACh + 0.00009} if ACh > 0 \\ 0 otherwise \end{matrix} P_up = P_up_basal (1 - b_up) j_Ca_dif = \frac{Ca_sub - Cai}{tau_dif_Ca} j_tr = \frac{Ca_nsr - Ca_jsr}{tau_tr} F_PLBp = \{\begin{matrix} \frac{3.3931 {PLBp}^{4.0695}}{{0.2805}^{4.0695} + {PLBp}^{4.0695}} if PLBp > 0.23 \\ \frac{1.698 {PLBp}^{13.5842}}{{0.2240}^{13.5842} + {PLBp}^{13.5842}} otherwise \end{matrix} j_up = \{\begin{matrix} \frac{0.9 P_up F_PLBp}{1 + \frac{K_up}{Cai}} if ANS > 0 \\ \frac{P_up}{1 + \frac{K_up}{Cai}} otherwise \end{matrix}

Component: Ca_buffering

\frac{d}{d time} (fTC) = delta_fTC delta_fTC = kf_TC Cai (1 - fTC) - (kb_TC fTC) \frac{d}{d time} (fTMC) = delta_fTMC delta_fTMC = kf_TMC Cai (1 - (fTMC + fTMM)) - (kb_TMC fTMC) \frac{d}{d time} (fTMM) = delta_fTMM delta_fTMM = kf_TMM Mgi (1 - (fTMC + fTMM)) - (kb_TMM fTMM) \frac{d}{d time} (fCMi) = delta_fCMi delta_fCMi = kf_CM Cai (1 - fCMi) - (kb_CM fCMi) \frac{d}{d time} (fCMs) = delta_fCMs delta_fCMs = kf_CM Ca_sub (1 - fCMs) - (kb_CM fCMs) \frac{d}{d time} (fCQ) = delta_fCQ delta_fCQ = kf_CQ Ca_jsr (1 - fCQ) - (kb_CQ fCQ)

Component: Ca_dynamics

BAPTA = \{\begin{matrix} 10 if BAPTA_10_mM > 0 \land time > T \\ 0 otherwise \end{matrix} \frac{d}{d time} (Cai) = \frac{1 (j_Ca_dif V_sub - (j_up V_nsr))}{V_i} - (CM_tot delta_fCMi + TC_tot delta_fTC + TMC_tot delta_fTMC) - (kfBAPTA Cai (BAPTA - fBAPTA) - (kbBAPTA fBAPTA)) \frac{d}{d time} (fBAPTA) = kfBAPTA Cai (BAPTA - fBAPTA) - (kbBAPTA fBAPTA) \frac{d}{d time} (Ca_sub) = \frac{j_SRCarel V_jsr}{V_sub} - (\frac{i_siCa + i_CaT - (2 i_NaCa)}{2 F V_sub} + j_Ca_dif + CM_tot delta_fCMs) - (kfBAPTA Ca_sub (BAPTA - fBAPTA_sub) - (kbBAPTA fBAPTA_sub)) \frac{d}{d time} (fBAPTA_sub) = kfBAPTA Ca_sub (BAPTA - fBAPTA_sub) - (kbBAPTA fBAPTA_sub) \frac{d}{d time} (Ca_nsr) = j_up - (\frac{j_tr V_jsr}{V_nsr}) \frac{d}{d time} (Ca_jsr) = j_tr - (j_SRCarel + CQ_tot delta_fCQ)

Component: Cell_parameters

V_cell = 0.000000001 π {R_cell}^{2} L_cell V_sub = 0.000000001 2 π L_sub (R_cell - (\frac{L_sub}{2})) L_cell V_jsr = V_jsr_part V_cell V_i = V_i_part V_cell - V_sub V_nsr = V_nsr_part V_cell

Component: i_to

i_to = g_to (V - E_K) q r

Component: i_to_q_gate

q_infinity = \frac{1}{1 + ⅇ^{\frac{V + 49}{13}}} tau_q = 0.001 0.6 (\frac{65.17}{0.57 ⅇ^{(- 0.08) (V + 44)} + 0.065 ⅇ^{0.1 (V + 45.93)}} + 10.1) \frac{d}{d time} (q) = \frac{q_infinity - q}{tau_q}

Component: i_to_r_gate

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

Component: i_Kr

i_Kr = g_Kr (V - E_K) (0.9 paF + 0.1 paS) piy

Component: i_Kr_pa_gate

alfapaF = \frac{\frac{1}{1 + ⅇ^{\frac{- (V + 23.2)}{6.6}}}}{\frac{0.84655354}{37.2 ⅇ^{\frac{V}{11.9}} + 0.96 ⅇ^{\frac{- V}{18.5}}}} betapaF = 4 (\frac{37.2 ⅇ^{\frac{V}{15.9}} + 0.96 ⅇ^{\frac{- V}{22.5}}}{0.84655354} - (\frac{\frac{1}{1 + ⅇ^{\frac{- (V + 23.2)}{10.6}}}}{\frac{0.84655354}{37.2 ⅇ^{\frac{V}{15.9}} + 0.96 ⅇ^{\frac{- V}{22.5}}}})) pa_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 14.8)}{8.5}}} tau_paS = \frac{0.84655354}{4.2 ⅇ^{\frac{V}{17}} + 0.15 ⅇ^{\frac{- V}{21.6}}} tau_paF = \frac{1}{30 ⅇ^{\frac{V}{10}} + ⅇ^{\frac{- V}{12}}} \frac{d}{d time} (paS) = \frac{pa_infinity - paS}{tau_paS} \frac{d}{d time} (paF) = \frac{pa_infinity - paF}{tau_paF}

Component: i_Kr_pi_gate

tau_pi = \frac{1}{100 ⅇ^{\frac{- V}{54.645}} + 656 ⅇ^{\frac{V}{106.157}}} pi_infinity = \frac{1}{1 + ⅇ^{\frac{V + 28.6}{17.1}}} \frac{d}{d time} (piy) = \frac{pi_infinity - piy}{tau_pi}

Component: i_Ks

g_Ks = \{\begin{matrix} 1.2 0.0016576 if Iso_1_uM > 0 \\ 0.0016576 otherwise \end{matrix} ANS_cond_i_Ks = \{\begin{matrix} (- 0.2152) + \frac{K_i_Ks1 {PKA}^{10.0808}}{{K_05i_Ks1}^{10.0808} + {PKA}^{10.0808}} if ANS > 0 \\ 0 otherwise \end{matrix} E_Ks = RTONF \ln (\frac{Ko + 0 Nao}{Ki + 0 Nai}) i_Ks = g_Ks (1 + ANS_cond_i_Ks) (V - E_Ks) n^{2}

Component: i_Ks_n_gate

Iso_shift = \{\begin{matrix} - 14 if Iso_1_uM > 0 \\ 0 otherwise \end{matrix} ANS_shift_i_Ks_n_gate = \{\begin{matrix} - (\frac{K_i_Ks2 {PKA}^{9.281}}{{K_05i_Ks2}^{9.281} + {PKA}^{9.281}} - 18.76) if ANS > 0 \\ 0 otherwise \end{matrix} n_infinity = \frac{\frac{14}{1 + ⅇ^{\frac{- (V - 40 - Iso_shift - ANS_shift_i_Ks_n_gate)}{12}}}}{\frac{14}{1 + ⅇ^{\frac{- (V - 40 - Iso_shift - ANS_shift_i_Ks_n_gate)}{12}}} + 1 ⅇ^{\frac{- (V - Iso_shift - ANS_shift_i_Ks_n_gate)}{45}}} tau_n = \frac{1}{alpha_n + beta_n} alpha_n = \frac{28}{1 + ⅇ^{\frac{- (V - 40 - Iso_shift - ANS_shift_i_Ks_n_gate)}{3}}} beta_n = 1 ⅇ^{\frac{- (V - Iso_shift - ANS_shift_i_Ks_n_gate - shift - 5)}{25}} \frac{d}{d time} (n) = \frac{n_infinity - n}{tau_n}

Component: i_KACh

i_KACh = \{\begin{matrix} K_g_ACh_max g_KACh (V - E_K) a if ANS > 0 \land ACh_cas > 0 \\ g_KACh (V - E_K) (1 + ⅇ^{\frac{V + 20}{20}}) a if ACh > 0 \\ 0 otherwise \end{matrix}

Component: i_KACh_a_gate

alpha_a = \{\begin{matrix} \frac{K_alpha 0.001 12.32}{K_alpha_add + \frac{K_alpha_div}{ACh_cas 1 e -6}} if ANS > 0 \\ \frac{3.5988 - 0.025641}{1 + \frac{0.0000012155}{{(1 e -6 ACh)}^{1.6951}}} + 0.025641 otherwise \end{matrix} beta_a = \{\begin{matrix} 0.001 K_beta ⅇ^{0.0133 (V + 40)} if ANS > 0 \\ 10 ⅇ^{0.0133 (V + 40)} otherwise \end{matrix} a_infinity = \frac{alpha_a}{alpha_a + beta_a} tau_a = \frac{1}{alpha_a + beta_a} \frac{d}{d time} (a) = \frac{a_infinity - a}{tau_a}

Component: cAMP

\frac{d}{d time} (cAMP) = diff_cAMP diff_cAMP = \frac{(kiso - kach) ATPi + k1 ATPi - (k2 cAMP) - (k3 cAMP)}{60} kiso = K_iso + \frac{0.1181 {Iso_cas}^{niso}}{{K_05iso}^{niso} + {Iso_cas}^{niso}} kach = \frac{K_ach 0.0146 {ACh_cas}^{nach}}{{K_05ach}^{nach} + {ACh_cas}^{nach}} k1 = K_k1 (K_ACI + \frac{K_AC}{1 + ⅇ^{\frac{K_Ca - (\frac{kb_CM fCMi}{kf_CM (1 - fCMi)})}{K_ACCa}}}) k2 = \frac{1.1 237.9851 {(cAMP 600)}^{5.101}}{{20.1077}^{6.101} + {(cAMP 600)}^{6.101}} k3 = \frac{kPKA {cAMP}^{nPKA - 1}}{{kPKA_cAMP}^{nPKA} + {cAMP}^{nPKA}}

Component: PLBp

\frac{d}{d time} (PLBp) = diff_PLBp diff_PLBp = \frac{k4 - k5}{60} k4 = \frac{kPLBp {PKA}^{nPLB}}{{kPKA_PLB}^{nPLB} + {PKA}^{nPLB}} k5 = \frac{kPP1 PP1 PLBp}{kPP1_PLB + PLBp}

Component: PKA

PKA = \{\begin{matrix} (- 0.9483029) ⅇ^{(- cAMP) 600 0.06561479} + 0.97781646 if cAMP 600 < 25.87 \\ (- 0.45260509) ⅇ^{(- cAMP) 600 0.03395094} + 0.99221714 otherwise \end{matrix}

Component: ATPi

ATPi = \frac{ATPi_max (\frac{kATP {(\frac{cAMP 100}{cAMPb})}^{nATP}}{kATP05 + {(\frac{cAMP 100}{cAMPb})}^{nATP}} - KATP_min)}{100}