Model Mathematics

\frac{d}{d time} (NADHm) = gamma (Jpdh - Jo) NADm = NADtot - NADHm ATPm = Amtot - ADPm \frac{d}{d time} (ADPm) = gamma (Jant - Jf1f0) RATm = \frac{ATPm}{ADPm} \frac{d}{d time} (PSIm) = \frac{Jh_res - Jh_atp - Jant - Jh_leak - Jnaca - (Juni 2)}{Cmito} \frac{d}{d time} (Cam) = (- fm) Jmito

delta = \frac{3.9}{53.2} ATPc = Ac_tot - ADPc Jh_yd = \frac{(khyd c + Jhydbas) ATPc}{1} \frac{d}{d time} (ADPc) = (- delta) Jant + Jh_yd 1

\frac{d}{d time} (FBP) = \frac{Fproto - FBP}{0.1} \frac{d}{d time} (c) = \frac{Cproto - c}{0.1}

pulse1 = (Ctest - Chold) (heav_Cton1 - heav_Ctoff1) pulse2 = (Ctest - Chold) (heav_Cton2 - heav_Ctoff2) pulse3 = (Ctest - Chold) (heav_Cton3 - heav_Ctoff3) Cproto = Chold + pulse1 + pulse2 + pulse3 Fproto = Fhold + (Ftest - Fhold) (heav_on - heav_off) heav_on = \{\begin{matrix} 1 if time - Fton \geq 0 \\ 0 otherwise \end{matrix} heav_off = \{\begin{matrix} 1 if time - Ftoff \geq 0 \\ 0 otherwise \end{matrix} heav_Cton1 = \{\begin{matrix} 1 if time - Cton1 \geq 0 \\ 0 otherwise \end{matrix} heav_Cton2 = \{\begin{matrix} 1 if time - Cton2 \geq 0 \\ 0 otherwise \end{matrix} heav_Cton3 = \{\begin{matrix} 1 if time - Cton3 \geq 0 \\ 0 otherwise \end{matrix} heav_Ctoff1 = \{\begin{matrix} 1 if time - Ctoff1 \geq 0 \\ 0 otherwise \end{matrix} heav_Ctoff2 = \{\begin{matrix} 1 if time - Ctoff2 \geq 0 \\ 0 otherwise \end{matrix} heav_Ctoff3 = \{\begin{matrix} 1 if time - Ctoff3 \geq 0 \\ 0 otherwise \end{matrix}

Juni = (p21 PSIm - p22) c^{2} Jpdh = \frac{p1}{p2 + \frac{NADHm}{NADm}} \frac{Cam}{p3 + Cam} Jgpdh Jnaca = \frac{p23 Cam}{c} ⅇ^{p24 PSIm} Jmito = Jnaca - Juni Jh_leak = p17 PSIm + p18 MM1 = \frac{p4 NADHm}{p5 + NADHm} MM2 = \frac{p8 NADHm}{p9 + NADHm} Jo = \frac{MM1}{1 + ⅇ^{\frac{PSIm - p6}{p7}}} Jh_res = \frac{MM2}{1 + ⅇ^{\frac{PSIm - p10}{p11}}} b13 = \frac{p12 p13}{p13 + ATPm} Jh_atp = \frac{b13}{1 + ⅇ^{\frac{p14 - PSIm}{p15}}} b2 = \frac{p16 p13}{p13 + ATPm} Jf1f0 = \frac{b2}{1 + ⅇ^{\frac{p14 - PSIm}{p15}}} Jant = \frac{\frac{p19 RATm}{RATm + p20}}{ⅇ^{(- 0.5) FRT PSIm}} Jgpdh = kgpdh \sqrt{\frac{FBP}{1}}