v1
∅ > Pyr
v2
Pyr + NAD > AcCoA + NADH
v3
OAA + AcCoA > Cit
v4
Cit + NAD > KG + NADH
v5
KG + ADP + {2.0}NAD > OAA + ATP + {2.0}NADH
v6
OAA = KG
v7
Pyr + ATP > OAA + ADP
v8
OAA > ∅
vANT
ATP > ADP
vATP
ADP + iP + {3.0}He = ATP + H2O + {3.0}H
vleak
He > H
vresp
NADH + {0.5}O2 + {11.0}H > NAD + H2O + {10.0}He
NADH = Nt - NAD
ADP = At - ATP
DeltaGtransport = 1.2 * F * DeltaPsi
ATPcrit = At / (1.0 + exp(-3.0 * DeltaGtransport / (R * T)) / (Kapp * iP))
JATP = kATP * (2.0 / (1.0 + exp(b * (ATP - ATPcrit))) - 1.0)
Jresp = kresp * ((Nt - NAD) / (K + Nt - NAD)) / (1.0 + exp(a * (DeltaPsi - DeltaPsim)))
Jleak = kleak * DeltaPsi
JANT = kANT * ATP
Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.