begin name
teusink
end name

begin reactions
v1 {1.0}Prb + {1.0}GLCi = {1.0}G6P
v10 {1.0}P2G = {1.0}PEP
v11 {1.0}PEP = {1.0}PYR + {1.0}Prb
v12 {1.0}PYR = {1.0}$CO2 + {1.0}ACE
v13 {2.0}ACE + {3.0}NAD = {1.0}$SUCC + {3.0}NADH
v14 {1.0}$GLCo = {1.0}GLCi
v15 {1.0}ACE + {1.0}NADH = {1.0}NAD + {1.0}$ETOH
v16 {1.0}NADH + {1.0}TRIO = {1.0}$GLY + {1.0}NAD
v17 {1.0}Prb = {1.0}$X
v2 {1.0}G6P = {1.0}F6P
v3 {1.0}G6P + {1.0}Prb = {1.0}$Glyc
v4 {2.0}G6P + {1.0}Prb = {1.0}$Trh
v5 {1.0}F6P + {1.0}Prb = {1.0}F16P
v6 {1.0}F16P = {2.0}TRIO
v7 {1.0}TRIO + {1.0}NAD = {1.0}BPG + {1.0}NADH
v8 {1.0}BPG = {1.0}P3G + {1.0}Prb
v9 {1.0}P3G = {1.0}P2G
end reactions

begin rate equations
v1 = (VmGLK*(-((G6P[t]*(SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/((1 - 4*KeqAK)*KeqGLK)) + (GLCi[t]*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/(2 - 8*KeqAK)))/(KmGLKATP*KmGLKGLCi*(1 + G6P[t]/KmGLKG6P + GLCi[t]/KmGLKGLCi)*(1 + (SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((1 - 4*KeqAK)*KmGLKADP) + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KmGLKATP)))
v10 = (VmENO*(P2G[t] - PEP[t]/KeqENO))/(KmENOP2G*(1 + P2G[t]/KmENOP2G + PEP[t]/KmENOPEP))
v11 = (VmPYK*((PEP[t]*(SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/(1 - 4*KeqAK) - ((-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)*PYR[t])/((2 - 8*KeqAK)*KeqPYK)))/(KmPYKADP*KmPYKPEP*(1 + (SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((1 - 4*KeqAK)*KmPYKADP) + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KmPYKATP))*(1 + PEP[t]/KmPYKPEP + PYR[t]/KmPYKPYR))
v12 = (VmPDC*PYR[t]^nPDC)/(KmPDCPYR^nPDC*(1 + PYR[t]^nPDC/KmPDCPYR^nPDC))
v13 = KSUCC*ACE[t]
v14 = (VmGLT*(GLCo - GLCi[t]/KeqGLT))/(KmGLTGLCo*(1 + GLCo/KmGLTGLCo + GLCi[t]/KmGLTGLCi + (0.91*GLCo*GLCi[t])/(KmGLTGLCi*KmGLTGLCo)))
v15 = -((VmADH*(ETOH*NAD[t] - (ACE[t]*NADH[t])/KeqADH))/(KiADHNAD*KmADHETOH*(1 + (ETOH*KmADHNAD)/(KiADHNAD*KmADHETOH) + (KmADHNADH*ACE[t])/(KiADHNADH*KmADHACE) + NAD[t]/KiADHNAD + (ETOH*NAD[t])/(KiADHNAD*KmADHETOH) + (ETOH*ACE[t]*NAD[t])/(KiADHACE*KiADHNAD*KmADHETOH) + (KmADHNADH*ACE[t]*NAD[t])/(KiADHNAD*KiADHNADH*KmADHACE) + NADH[t]/KiADHNADH + (ETOH*KmADHNAD*NADH[t])/(KiADHNAD*KiADHNADH*KmADHETOH) + (ACE[t]*NADH[t])/(KiADHNADH*KmADHACE) + (ETOH*ACE[t]*NADH[t])/(KiADHETOH*KiADHNADH*KmADHACE))))
v16 = (VmG3PDH*(-((GLY*NAD[t])/KeqG3PDH) + (NADH[t]*TRIO[t])/(1 + KeqTPI)))/(KmG3PDHDHAP*KmG3PDHNADH*(1 + NAD[t]/KmG3PDHNAD + NADH[t]/KmG3PDHNADH)*(1 + GLY/KmG3PDHGLY + TRIO[t]/((1 + KeqTPI)*KmG3PDHDHAP)))
v17 = (KATPASE*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/(2 - 8*KeqAK)
v2 = (VmPGI*(-(F6P[t]/KeqPGI) + G6P[t]))/(KmPGIG6P*(1 + F6P[t]/KmPGIF6P + G6P[t]/KmPGIG6P))
v3 = KGLYCOGEN
v4 = KTREHALOSE
v5 = (gR*VmPFK*F6P[t]*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)*(1 + F6P[t]/KmPFKF6P + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KmPFKATP) + (gR*F6P[t]*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P)))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P*((L0*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P[t])/KPFKF16BP)^2*(1 + (2*CPFKAMP*KeqAK*(SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)^2)/((-1 + 4*KeqAK)*KPFKAMP*(SUMAXP - Prb[t] + 4*KeqAK*Prb[t] - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)))^2*(1 + (CiPFKATP*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/((2 - 8*KeqAK)*KiPFKATP))^2*(1 + (CPFKATP*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/((2 - 8*KeqAK)*KmPFKATP))^2)/((1 + F26BP/KPFKF26BP + F16P[t]/KPFKF16BP)^2*(1 + (2*KeqAK*(SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)^2)/((-1 + 4*KeqAK)*KPFKAMP*(SUMAXP - Prb[t] + 4*KeqAK*Prb[t] - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)))^2*(1 + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KiPFKATP))^2) + (1 + F6P[t]/KmPFKF6P + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KmPFKATP) + (gR*F6P[t]*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P))^2))
v6 = (VmALD*(F16P[t] - (KeqTPI*TRIO[t]^2)/(KeqALD*(1 + KeqTPI)^2)))/(KmALDF16P*(1 + F16P[t]/KmALDF16P + TRIO[t]/((1 + KeqTPI)*KmALDDHAP) + (KeqTPI*TRIO[t])/((1 + KeqTPI)*KmALDGAP) + (KeqTPI*F16P[t]*TRIO[t])/((1 + KeqTPI)*KmALDF16P*KmALDGAPi) + (KeqTPI*TRIO[t]^2)/((1 + KeqTPI)^2*KmALDDHAP*KmALDGAP)))
v7 = (-((VmGAPDHr*BPG[t]*NADH[t])/(KmGAPDHBPG*KmGAPDHNADH)) + (KeqTPI*VmGAPDHf*NAD[t]*TRIO[t])/((1 + KeqTPI)*KmGAPDHGAP*KmGAPDHNAD))/((1 + NAD[t]/KmGAPDHNAD + NADH[t]/KmGAPDHNADH)*(1 + BPG[t]/KmGAPDHBPG + (KeqTPI*TRIO[t])/((1 + KeqTPI)*KmGAPDHGAP)))
v8 = (VmPGK*((KeqPGK*BPG[t]*(SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/(1 - 4*KeqAK) - (P3G[t]*(-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5))/(2 - 8*KeqAK)))/(KmPGKATP*KmPGKP3G*(1 + BPG[t]/KmPGKBPG + P3G[t]/KmPGKP3G)*(1 + (SUMAXP - (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((1 - 4*KeqAK)*KmPGKADP) + (-SUMAXP + Prb[t] - 4*KeqAK*Prb[t] + (SUMAXP^2 - 2*SUMAXP*Prb[t] + 8*KeqAK*SUMAXP*Prb[t] + Prb[t]^2 - 4*KeqAK*Prb[t]^2)^0.5)/((2 - 8*KeqAK)*KmPGKATP)))
v9 = (VmPGM*(-(P2G[t]/KeqPGM) + P3G[t]))/(KmPGMP3G*(1 + P2G[t]/KmPGMP2G + P3G[t]/KmPGMP3G))
end rate equations

begin parameters
CPFKAMP = 0.0845
CPFKATP = 3.0
CPFKF16BP = 0.397
CPFKF26BP = 0.0174
CPFKF6P = 0.0
CiPFKATP = 100.0
F26BP = 0.02
KATPASE = 39.5
KGLYCOGEN = 6.0
KPFKAMP = 0.0995
KPFKF16BP = 0.111
KPFKF26BP = 0.000682
KSUCC = 21.4
KTREHALOSE = 2.4
KeqADH = "6.9e-05"
KeqAK = 0.45
KeqALD = 0.069
KeqENO = 6.7
KeqG3PDH = 4300.0
KeqGLK = 3800.0
KeqGLT = 1.0
KeqPGI = 0.314
KeqPGK = 3200.0
KeqPGM = 0.19
KeqPYK = 6500.0
KeqTPI = 0.045
KiADHACE = 1.1
KiADHETOH = 90.0
KiADHNAD = 0.92
KiADHNADH = 0.031
KiPFKATP = 0.65
KmADHACE = 1.11
KmADHETOH = 17.0
KmADHNAD = 0.17
KmADHNADH = 0.11
KmALDDHAP = 2.4
KmALDF16P = 0.3
KmALDGAP = 2.0
KmALDGAPi = 10.0
KmENOP2G = 0.04
KmENOPEP = 0.5
KmG3PDHDHAP = 0.4
KmG3PDHGLY = 1.0
KmG3PDHNAD = 0.93
KmG3PDHNADH = 0.023
KmGAPDHBPG = 0.0098
KmGAPDHGAP = 0.21
KmGAPDHNAD = 0.09
KmGAPDHNADH = 0.06
KmGLKADP = 0.23
KmGLKATP = 0.15
KmGLKG6P = 30.0
KmGLKGLCi = 0.08
KmGLTGLCi = 1.1918
KmGLTGLCo = 1.1918
KmPDCPYR = 4.33
KmPFKATP = 0.71
KmPFKF6P = 0.1
KmPGIF6P = 0.3
KmPGIG6P = 1.4
KmPGKADP = 0.2
KmPGKATP = 0.3
KmPGKBPG = 0.003
KmPGKP3G = 0.53
KmPGMP2G = 0.08
KmPGMP3G = 1.2
KmPYKADP = 0.53
KmPYKATP = 1.5
KmPYKPEP = 0.14
KmPYKPYR = 21.0
L0 = 0.66
SUMAXP = 4.1
VmADH = 810.0
VmALD = 322.258
VmENO = 365.806
VmG3PDH = 70.15
VmGAPDHf = 1184.52
VmGAPDHr = 6549.68
VmGLK = 226.452
VmGLT = 97.264
VmPDC = 174.194
VmPFK = 182.903
VmPGI = 339.677
VmPGK = 1306.45
VmPGM = 2525.81
VmPYK = 1088.71
gR = 5.12
nPDC = 1.9
CO2 = 1.0
ETOH = 50.0
GLCo = 50.0
GLY = 0.15
Glyc = 0.0
SUCC = 0.0
Trh = 0.0
X = 0.0
defaultcompartment = 1.0
end parameters

begin initial conditions
ACE[0] = ACEi
BPG[0] = BPGi
F16P[0] = F16Pi
F6P[0] = F6Pi
G6P[0] = G6Pi
GLCi[0] = GLCii
NAD[0] = NADi
NADH[0] = NADHi
P2G[0] = P2Gi
P3G[0] = P3Gi
PEP[0] = PEPi
PYR[0] = PYRi
Prb[0] = Prbi
TRIO[0] = TRIOi
end initial conditions

begin initial values
ACEi = 0.04
BPGi = 0.0
F16Pi = 0.1
F6Pi = 0.28
G6Pi = 1.39
GLCii = 0.087
NADi = 1.2
NADHi = 0.39
P2Gi = 0.1
P3Gi = 0.1
PEPi = 0.1
PYRi = 3.36
Prbi = 5.0
TRIOi = 5.17
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
v_1 http://identifiers.org/interpro/IPR001312
v_1 http://identifiers.org/ec-code/2.7.1.1
v_10 urn:miriam:kegg.reaction:R01518
v_11 urn:miriam:kegg.reaction:R00658
v_12 urn:miriam:kegg.reaction:R00200
v_13 urn:miriam:kegg.reaction:R00224
v_15 urn:miriam:kegg.reaction:R00754
v_16 http://identifiers.org/ec-code/1.1.1.94
v_17 urn:miriam:obo.go:0006200
v_2 urn:miriam:kegg.reaction:R00299
v_5 urn:miriam:kegg.reaction:R02740
v_6 urn:miriam:kegg.reaction:R04779
v_6 http://identifiers.org/ec-code/2.7.1.11
v_7 urn:miriam:kegg.reaction:R01068
v_8 urn:miriam:kegg.reaction:R01061
v_9 urn:miriam:kegg.reaction:R01512
end process annotations

begin species annotations
ACE http://identifiers.org/obo.chebi/CHEBI%3A15343
BPG urn:miriam:obo.chebi:28907
BPG urn:miriam:kegg.compound:C00236
CO2 http://identifiers.org/obo.chebi/CHEBI%3A16526
ETOH http://identifiers.org/obo.chebi/CHEBI%3A16236
F16P urn:miriam:obo.chebi:16905
F16P urn:miriam:kegg.compound:C00354
F6P urn:miriam:kegg.compound:C00085
F6P urn:miriam:obo.chebi:15946
G6P urn:miriam:kegg.compound:C00092
G6P urn:miriam:obo.chebi:17665
GLCi http://identifiers.org/obo.chebi/CHEBI%3A17234
GLCo http://identifiers.org/obo.chebi/CHEBI%3A17234
GLY http://identifiers.org/obo.chebi/CHEBI%3A17754
Glyc http://identifiers.org/obo.chebi/CHEBI%3A28087
NAD urn:miriam:obo.chebi:15846
NADH urn:miriam:obo.chebi:16908
NADH urn:miriam:kegg.compound:C00004
P2G urn:miriam:kegg.compound:C00631
P2G urn:miriam:obo.chebi:17835
P3G urn:miriam:kegg.compound:C00197
P3G urn:miriam:obo.chebi:17794
PEP urn:miriam:obo.chebi:18021
PEP urn:miriam:kegg.compound:C00074
PYR urn:miriam:obo.chebi:15361
PYR urn:miriam:kegg.compound:C00022
SUCC http://identifiers.org/obo.chebi/CHEBI%3A15741
SUCC http://identifiers.org/kegg.compound/C00042
TRIO http://identifiers.org/obo.chebi/CHEBI%3A29052
TRIO http://identifiers.org/kegg.compound/C00111
Trh http://identifiers.org/obo.chebi/CHEBI%3A27082
end species annotations

begin units
end units