Component: environment

Component: membrane

$d Vm d time = − (Ik + Ica + Ikca + Ikatp) / cm$

Component: Ik

$Ik = gK n (Vm − VK)$

Component: n

$d n d time = (n_infinity − n) / tau_n n_infinity = 1.0 / (1.0 + e (vn − Vm) / sn)$

Component: Ica

$Ica = gCa m_infinity (Vm − VCa)$

Component: m

$m_infinity = 1.0 / (1.0 + e (v − Vm) / sm)$

Component: Ikca

$Ikca = gkCa / (1.0 + (kd / c) 2.0) (Vm − VK)$

Component: Ikatp

$Ikatp = gkATP_katpo (Vm − VK) katpo = topo / bottomo topo = 0.08 (1.0 + 2.0 mgadp / 17.0) + 0.89 (mgadp / 17.0) 2.0 bottomo = (1.0 + mgadp / 17.0) 2.0 (1.0 + adp3m / 26.0 + atp4m / 1.0) mgadp = 0.165 adp adp3m = 0.135 adp atp4m = 0.05 atp$

Component: JGPDH

$JGPDH = kGPDH | FBP |$

Component: G6P

$d G6P d time = JGK_ms − JPFK_ms$

Component: FBP

$d FBP d time = JPFK_ms − 0.5 JGPDH$

Component: F6P

$F6P = 0.3 G6P$

Component: JPFK

$JPFK_ms = kappa JPFK JPFK = (lambda VmaxPFK topa16 + VmaxPFK topb) / bottom16$

$weight2 = atp 2 / k4 topa2 = topa1 bottom2 = bottom1 + weight2$

$weight3 = F6P 2 / k3 topa3 = topa2 + weight3 bottom3 = bottom2 + weight3$

$weight4 = (F6P atp) 2 / f43 k3 k4 topa4 = topa3 + weight4 bottom4 = bottom3 + weight4$

$weight5 = FBP / k2 topa5 = topa4 bottom5 = bottom4 + weight5$

$weight6 = FBP atp 2 / k2 k4 f42 topa6 = topa5 bottom6 = bottom5 + weight6$

$weight7 = FBP F6P 2 / k2 k3 f23 topa7 = topa6 + weight7 bottom7 = bottom6 + weight7$

$weight8 = FBP F6P 2 atp 2 / k2 k3 k4 f23 f42 f43 topa8 = topa7 + weight8 bottom8 = bottom7 + weight8$

$weight9 = AMP / k1 topa9 = topa8 bottom9 = bottom8 + weight9$

$weight10 = AMP atp 2 / k1 k4 f41 topa10 = topa9 bottom10 = bottom9 + weight10$

$weight11 = AMP F6P 2 / k1 k3 f13 topa11 = topa10 + weight11 bottom11 = bottom10 + weight11$

$weight12 = AMP F6P 2 atp 2 / k1 k3 k4 f13 f41 f43 topa12 = topa11 + weight12 bottom12 = bottom11 + weight12$

$weight13 = AMP FBP / k1 k2 topa13 = topa12 bottom13 = bottom12 + weight13$

$weight14 = AMP FBP atp 2 / k1 k2 k4 f42 f41 topa14 = topa13 bottom14 = bottom13 + weight14$

$weight15 = AMP FBP F6P 2 / k1 k2 k3 f23 f13 topa15 = topa14 bottom15 = bottom14 + weight15 topb = weight15$

$weight16 = AMP FBP F6P 2 atp 2 / k1 k2 k3 k4 f23 f13 f42 f41 f43 topa16 = topa15 + weight16 bottom16 = bottom15 + weight16$

Component: JPDH

$JPDH = p1 / (p2 + NADHm / NADm) Cam / (p3 + Cam) (JGPDH + JGPDHbas)$

Component: NADHm

$d NADHm d time = gamma (JPDH − JO)$

Component: JO

$JO = p4 NADHm / (p5 + NADHm) / (1.0 + e (delta_psi − p6) / p7)$

Component: NADm

$NADm = NADmtot − NADHm$

Component: delta_psi

$d delta_psi d time = (JHres − JHatp + JANT + JHleak + JNaCa + 2.0 Juni) / Cmito$

Component: JHres

$JHres = p8 NADHm / (p9 + NADHm) / (1.0 + e (delta_psi − p10) / p11)$

Component: JF1F0

$JF1F0 = p16 p13 / (p13 + ATPm) / (1.0 + e (p14 − delta_psi) / p15)$

Component: JHatp

$JHatp = 3.0 JF1F0$

Component: JGK

$JGK_ms = kappa JGK$

Component: JHleak

$JHleak = p17 delta_psi + p18$

Component: JANT

$JANT = p19 RATm / (RATm + p20) / e -0.5 FRT delta_psi$

Component: RATm

$RATm = ATPm / ADPm$

Component: Juni

$Juni = (p21 delta_psi − p22) c 2.0$

Component: JNaCa

$JNaCa = p23 / c Cam e p24 delta_psi$

Component: Cam

$d Cam d time = − fmito Jmito$

Component: Jmito

$Jmito = JNaCa − Juni$

Component: ADPm

$d ADPm d time = gamma (JANT − JF1F0)$

Component: ATPm

$ATPm = Amtot − ADPm$

Component: adp

$d adp d time = delta JANT + Jhyd$

Component: Jhyd

$Jhyd = (khyd c + khydbas) atp$

Component: atp

$atp = atot − adp$

Component: c

$d c d time = fcyt (Jmem + Jer + delta Jmito)$

Component: Jmem

$Jmem = − (alpha Ica + kPMCA (c − Cbas))$

Component: Jleak

$Jleak = pleak (Caer − c)$

Component: JSERCA

$JSERCA = kSERCA c$

Component: Jer

$Jer = Jleak − JSERCA$

Component: Caer

$d Caer d time = − fer Vc_Ver Jer$

Component: model_parameters

$delta = 3.9 / 53.2$