begin name
vanEunen2
end name

begin reactions
v1 {1.0}$EXTERNAL = {1.0}GLCi
v10 {1.0}BPG = {1.0}P3G
v11 {1.0}P3G = {1.0}P2G
v12 {1.0}P2G = {1.0}PEP
v13 {1.0}PEP = {1.0}PYR
v14 {1.0}PYR = {1.0}ACALD
v15 {2.0}PYR = {3.0}NADH
v16 {1.0}ACALD + {1.0}NADH = {1.0}$EXTERNAL
v17 {1.0}ACALD = {1.0}NADH
v2 {1.0}GLCi = {1.0}G6P
v3 {1.0}$EXTERNAL = {2.0}GLCi
v4 {1.0}G6P = {1.0}F6P
v5 {2.0}G6P = {1.0}$EXTERNAL
v6 {1.0}F6P = {1.0}F16P
v7 {1.0}F16P = {2.0}TRIO
v8 {1.0}NADH + {1.0}TRIO = {1.0}$EXTERNAL
v9 {1.0}TRIO = {1.0}BPG + {1.0}NADH
end reactions

begin rate equations
v1 = (VMAXGLT*(GLCo/KMGLTGLCo - GLCi[t]/(KEQGLT*KMGLTGLCo)))/(1 + GLCo/KMGLTGLCo + GLCi[t]/KMGLTGLCi + (0.91*GLCo*GLCi[t])/(KMGLTGLCi*KMGLTGLCo))
v10 = (VMAXPGK*((ADP*KEQPGK*BPG[t])/(KMPGKATP*KMPGKP3G) - (ATP*P3G[t])/(KMPGKATP*KMPGKP3G)))/((1 + ADP/KMPGKADP + ATP/KMPGKATP)*(1 + BPG[t]/KMPGKBPG + P3G[t]/KMPGKP3G))
v11 = (VMAXPGM*(-(P2G[t]/(KEQPGM*KMPGMP3G)) + P3G[t]/KMPGMP3G))/(1 + P2G[t]/KMPGMP2G + P3G[t]/KMPGMP3G)
v12 = (VMAXENO*(P2G[t]/KMENOP2G - PEP[t]/(KEQENO*KMENOP2G)))/(1 + P2G[t]/KMENOP2G + PEP[t]/KMENOPEP)
v13 = (ADP*VMAXPYK*PEP[t]*(1 + PEP[t]/KMPYKPEP)^(-1 + n10))/((ADP + KMPYKADP)*KMPYKPEP*(L10*((1 + ATP/KMPYKATP)/(1 + F16P[t]/KMPYKF16P))^n10 + (1 + PEP[t]/KMPYKPEP)^n10))
v14 = (VMAXPDC*PYR[t]^NHPDC)/(KMPDCPYR^NHPDC*(1 + PYR[t]^NHPDC/KMPDCPYR^NHPDC))
v15 = KSUC
v16 = -((VMAXADH*((ETOH*(NADt - NADH[t]))/(KIADHNAD*KMADHETOH) - (ACALD[t]*NADH[t])/(KEQADH*KIADHNAD*KMADHETOH)))/(1 + (ETOH*KMADHNAD)/(KIADHNAD*KMADHETOH) + (KMADHNADH*ACALD[t])/(KIADHNADH*KMADHACALD) + (NADt - NADH[t])/KIADHNAD + (ETOH*(NADt - NADH[t]))/(KIADHNAD*KMADHETOH) + (ETOH*ACALD[t]*(NADt - NADH[t]))/(KIADHACALD*KIADHNAD*KMADHETOH) + (KMADHNADH*ACALD[t]*(NADt - NADH[t]))/(KIADHNAD*KIADHNADH*KMADHACALD) + NADH[t]/KIADHNADH + (ETOH*KMADHNAD*NADH[t])/(KIADHNAD*KIADHNADH*KMADHETOH) + (ACALD[t]*NADH[t])/(KIADHNADH*KMADHACALD) + (ETOH*ACALD[t]*NADH[t])/(KIADHETOH*KIADHNADH*KMADHACALD)))
v17 = KACE*ACALD[t]
v2 = (VMAXHK*(-((ADP*G6P[t])/(KEQHK*KMHKATP*KMHKGLCi)) + (ATP*GLCi[t])/(KMHKATP*KMHKGLCi)))/((1 + ADP/KMHKADP + ATP/KMHKATP)*(1 + T6P/KIHKT6P + G6P[t]/KMHKG6P + GLCi[t]/KMHKGLCi))
v3 = KTRE1
v4 = (VMAXPGI*(-(F6P[t]/(KEQPGI*KMPGIG6P)) + G6P[t]/KMPGIG6P))/(1 + F6P[t]/KMPGIF6P + G6P[t]/KMPGIG6P)
v5 = KTRE2
v6 = (ATP*gR*VMAXPFK*F6P[t]*(1 + (ATP*F6P[t])/(KMPFKATP*KMPFKF6P) + (ATP*gR*F6P[t])/(KMPFKATP*KMPFKF6P)))/(KMPFKATP*KMPFKF6P*(((1 + (ATP*CiPFKATP)/KiPFKATP)^2*(1 + (ATP*CPFKATP)/KMPFKATP)^2*(1 + (AMP*CPFKAMP)/KPFKAMP)^2*L0*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P[t])/KPFKF16BP)^2)/((1 + ATP/KiPFKATP)^2*(1 + AMP/KPFKAMP)^2*(1 + F26BP/KPFKF26BP + F16P[t]/KPFKF16BP)^2) + (1 + (ATP*F6P[t])/(KMPFKATP*KMPFKF6P) + (ATP*gR*F6P[t])/(KMPFKATP*KMPFKF6P))^2))
v7 = (VMAXALD*(F16P[t]/KMALDF16P - (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))
v8 = KGLY
v9 = (CGAPDH*(-((VMAXGAPDHr*BPG[t]*NADH[t])/(KMGAPDHBPG*KMGAPDHNADH)) + (KEQTPI*VMAXGAPDHf*(NADt - NADH[t])*TRIO[t])/((1 + KEQTPI)*KMGAPDHGAP*KMGAPDHNAD)))/((1 + (NADt - NADH[t])/KMGAPDHNAD + NADH[t]/KMGAPDHNADH)*(1 + BPG[t]/KMGAPDHBPG + (KEQTPI*TRIO[t])/((1 + KEQTPI)*KMGAPDHGAP)))
end rate equations

begin parameters
ACE = 10.0
ADP = 0.81
AMP = 0.25
ATP = 3.92
CGAPDH = 1.0
CO2 = 1.0
CPFKAMP = 0.0845
CPFKATP = 3.0
CPFKF16BP = 0.397
CPFKF26BP = 0.0174
CiPFKATP = 100.0
ETOH = 25.0
EXTERNAL = 1.0
F26BP = 0.009
GLCo = 50.0
GLY = 10.0
KACE = 0.5
KEQADH = "6.9e-05"
KEQALD = 0.069
KEQENO = 6.7
KEQGLT = 1.0
KEQHK = 3800.0
KEQPGI = 0.314
KEQPGK = 3200.0
KEQPGM = 0.19
KEQTPI = 0.045
KGLY = 24.9
KIADHACALD = 1.1
KIADHETOH = 90.0
KIADHNAD = 0.92
KIADHNADH = 0.031
KIHKT6P = 0.04
KMADHACALD = 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
KMGAPDHBPG = 0.51
KMGAPDHGAP = 0.39
KMGAPDHNAD = 2.85
KMGAPDHNADH = 0.007
KMGLTGLCi = 11.0
KMGLTGLCo = 11.0
KMHKADP = 0.23
KMHKATP = 0.15
KMHKG6P = 30.0
KMHKGLCi = 0.08
KMPDCPYR = 6.36
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.3
KMPYKATP = 9.3
KMPYKF16P = 0.2
KMPYKPEP = 0.19
KPFKAMP = 0.0995
KPFKF16BP = 0.111
KPFKF26BP = 0.000682
KSUC = 0.0
KTRE1 = 2.2
KTRE2 = 0.0
KiPFKATP = 0.65
L0 = 0.66
L10 = 60000.0
NADt = 1.59
NHPDC = 1.9
T6P = 0.36
TREH = 2.0
VMAXADH = 744.0
VMAXALD = 153.0
VMAXENO = 285.0
VMAXGAPDHf = 1514.0
VMAXGAPDHr = 877.0
VMAXGLT = 121.0
VMAXHK = 223.0
VMAXPDC = 297.0
VMAXPFK = 165.0
VMAXPGI = 852.0
VMAXPGK = 3030.0
VMAXPGM = 748.0
VMAXPYK = 636.0
gR = 5.12
n10 = 4.0
defaultcompartment = 1.0
end parameters

begin initial conditions
ACALD[0] = ACALDi
BPG[0] = BPGi
F16P[0] = F16Pi
F6P[0] = F6Pi
G6P[0] = G6Pi
GLCi[0] = GLCii
NADH[0] = NADHi
P2G[0] = P2Gi
P3G[0] = P3Gi
PEP[0] = PEPi
PYR[0] = PYRi
TRIO[0] = TRIOi
end initial conditions

begin initial values
ACALDi = 0.04
BPGi = 1e-05
F16Pi = 14.63
F6Pi = 0.8
G6Pi = 4.22
GLCii = 0.1
NADHi = 0.29
P2Gi = 0.13
P3Gi = 0.96
PEPi = 0.12
PYRi = 3.48
TRIOi = 1.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units