#teusink

FIX: CO2 ETOH GLCo GLY Glyc SUCC Trh X

#reactions
v_1:
 {1.0}Prb + {1.0}GLCi = {1.0}G6P
  VmGLK * (-(G6P * (SUMAXP - (SUMAXP**2 - 2 * SUMAXP * Prb + 8 * KeqAK * SUMAXP * Prb + Prb**2 - 4 * KeqAK * Prb**2)**0.5) / ((1 - 4 * KeqAK) * KeqGLK)) + GLCi * (-SUMAXP + Prb - 4 * KeqAK * Prb + (SUMAXP**2 - 2 * SUMAXP * Prb + 8 * KeqAK * SUMAXP * Prb + Prb**2 - 4 * KeqAK * Prb**2)**0.5) / (2 - 8 * KeqAK)) / (KmGLKATP * KmGLKGLCi * (1 + G6P / KmGLKG6P + GLCi / KmGLKGLCi) * (1 + (SUMAXP - (SUMAXP**2 - 2 * SUMAXP * Prb + 8 * KeqAK * SUMAXP * Prb + Prb**2 - 4 * KeqAK * Prb**2)**0.5) / ((1 - 4 * KeqAK) * KmGLKADP) + (-SUMAXP + Prb - 4 * KeqAK * Prb + (SUMAXP**2 - 2 * SUMAXP * Prb + 8 * KeqAK * SUMAXP * Prb + Prb**2 - 4 * KeqAK * Prb**2)**0.5) / ((2 - 8 * KeqAK) * KmGLKATP)))
v_10:
 {1.0}P2G = {1.0}PEP
  VmENO/KmENOP2G*(P2G - PEP/KeqENO)/(1 + P2G/KmENOP2G + PEP/KmENOPEP)
v_11:
 {1.0}PEP = {1.0}PYR + {1.0}Prb
  VmPYK/(KmPYKPEP*KmPYKADP)*(PEP*(SUMAXP - (Prb**2 - 4*KeqAK*Prb**2 - 2*Prb*SUMAXP + 8*KeqAK*Prb*SUMAXP + SUMAXP**2)**0.5)/(1 - 4*KeqAK) - PYR*((Prb - 4*KeqAK*Prb - SUMAXP + (Prb**2 - 4*KeqAK*Prb**2 - 2*Prb*SUMAXP + 8*KeqAK*Prb*SUMAXP + SUMAXP**2)**0.5)/(2 - 8*KeqAK))/KeqPYK)/((1 + PEP/KmPYKPEP + PYR/KmPYKPYR)*(1 + ((Prb - 4*KeqAK*Prb - SUMAXP + (Prb**2 - 4*KeqAK*Prb**2 - 2*Prb*SUMAXP + 8*KeqAK*Prb*SUMAXP + SUMAXP**2)**0.5)/(2 - 8*KeqAK))/KmPYKATP + (SUMAXP - (Prb**2 - 4*KeqAK*Prb**2 - 2*Prb*SUMAXP + 8*KeqAK*Prb*SUMAXP + SUMAXP**2)**0.5)/(1 - 4*KeqAK)/KmPYKADP))
v_12:
 {1.0}PYR = {1.0}CO2 + {1.0}ACE
  VmPDC*(PYR**(nPDC)/KmPDCPYR**(nPDC))/(1 + PYR**(nPDC)/KmPDCPYR**(nPDC))
v_13:
 {2.0}ACE + {3.0}NAD = {1.0}SUCC + {3.0}NADH
  KSUCC*ACE
v_14:
 {1.0}GLCo = {1.0}GLCi
  VmGLT * (GLCo - GLCi / KeqGLT) / (KmGLTGLCo * (1 + GLCo / KmGLTGLCo + GLCi / KmGLTGLCi + 0.91 * GLCo * GLCi / (KmGLTGLCi * KmGLTGLCo)))
v_15:
 {1.0}ACE + {1.0}NADH = {1.0}NAD + {1.0}ETOH
  -(VmADH/(KiADHNAD*KmADHETOH)*(NAD*ETOH - NADH*ACE/KeqADH)/(1 + NAD/KiADHNAD + KmADHNAD*ETOH/(KiADHNAD*KmADHETOH) + KmADHNADH*ACE/(KiADHNADH*KmADHACE) + NADH/KiADHNADH + NAD*ETOH/(KiADHNAD*KmADHETOH) + KmADHNADH*NAD*ACE/(KiADHNAD*KiADHNADH*KmADHACE) + KmADHNAD*ETOH*NADH/(KiADHNAD*KmADHETOH*KiADHNADH) + NADH*ACE/(KiADHNADH*KmADHACE)+ NAD*ETOH*ACE/(KiADHNAD*KmADHETOH*KiADHACE) + ETOH*NADH*ACE/(KiADHETOH*KiADHNADH*KmADHACE)))
v_16:
 {1.0}NADH + {1.0}TRIO = {1.0}GLY + {1.0}NAD
  VmG3PDH * (-(GLY * NAD / KeqG3PDH) + NADH * TRIO / (1 + KeqTPI)) / (KmG3PDHDHAP * KmG3PDHNADH * (1 + NAD / KmG3PDHNAD + NADH / KmG3PDHNADH) * (1 + GLY / KmG3PDHGLY + TRIO / ((1 + KeqTPI) * KmG3PDHDHAP)))
v_17:
 {1.0}Prb = {1.0}X
  KATPASE*((Prb - 4*KeqAK*Prb - SUMAXP + (Prb**2 - 4*KeqAK*Prb**2 - 2*Prb*SUMAXP + 8*KeqAK*Prb*SUMAXP + SUMAXP**2)**0.5)/(2 - 8*KeqAK))
v_2:
 {1.0}G6P = {1.0}F6P
  VmPGI/KmPGIG6P*(G6P - F6P/KeqPGI)/(1 + G6P/KmPGIG6P + F6P/KmPGIF6P)
v_3:
 {1.0}G6P + {1.0}Prb = {1.0}Glyc
  KGLYCOGEN
v_4:
 {2.0}G6P + {1.0}Prb = {1.0}Trh
  KTREHALOSE
v_5:
 {1.0}F6P + {1.0}Prb = {1.0}F16P
  (gR*VmPFK*F6P*(-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)*(1 + F6P/KmPFKF6P + (-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)/((2 - 8*KeqAK)*KmPFKATP) + (gR*F6P*(-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P)))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P*((L0*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P)/KPFKF16BP)**2*(1 + (2*CPFKAMP*KeqAK*(SUMAXP - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)**2)/((-1 + 4*KeqAK)*KPFKAMP*(SUMAXP - Prb + 4*KeqAK*Prb - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)))**2*(1 + (CiPFKATP*(-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5))/((2 - 8*KeqAK)*KiPFKATP))**2*(1 + (CPFKATP*(-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5))/((2 - 8*KeqAK)*KmPFKATP))**2)/((1 + F26BP/KPFKF26BP + F16P/KPFKF16BP)**2*(1 + (2*KeqAK*(SUMAXP - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)**2)/((-1 + 4*KeqAK)*KPFKAMP*(SUMAXP - Prb + 4*KeqAK*Prb - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)))**2*(1 + (-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)/((2 - 8*KeqAK)*KiPFKATP))**2) + (1 + F6P/KmPFKF6P + (-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)/((2 - 8*KeqAK)*KmPFKATP) + (gR*F6P*(-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5))/((2 - 8*KeqAK)*KmPFKATP*KmPFKF6P))**2))
v_6:
 {1.0}F16P = {2.0}TRIO
  (VmALD*(F16P - (KeqTPI*TRIO**2)/(KeqALD*(1 + KeqTPI)**2)))/(KmALDF16P*(1 + F16P/KmALDF16P + TRIO/((1 + KeqTPI)*KmALDDHAP) + (KeqTPI*TRIO)/((1 + KeqTPI)*KmALDGAP) + (KeqTPI*F16P*TRIO)/((1 + KeqTPI)*KmALDF16P*KmALDGAPi) + (KeqTPI*TRIO**2)/((1 + KeqTPI)**2*KmALDDHAP*KmALDGAP)))
v_7:
 {1.0}TRIO + {1.0}NAD = {1.0}BPG + {1.0}NADH
  (-((VmGAPDHr*BPG*NADH)/(KmGAPDHBPG*KmGAPDHNADH)) + (KeqTPI*VmGAPDHf*NAD*TRIO)/((1 + KeqTPI)*KmGAPDHGAP*KmGAPDHNAD))/((1 + NAD/KmGAPDHNAD + NADH/KmGAPDHNADH)*(1 + BPG/KmGAPDHBPG + (KeqTPI*TRIO)/((1 + KeqTPI)*KmGAPDHGAP)))
v_8:
 {1.0}BPG = {1.0}P3G + {1.0}Prb
  (VmPGK*((KeqPGK*BPG*(SUMAXP - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5))/(1 - 4*KeqAK) - ((-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)*P3G)/(2 - 8*KeqAK)))/(KmPGKATP*KmPGKP3G*(1 + (SUMAXP - (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)/((1 - 4*KeqAK)*KmPGKADP) + (-SUMAXP + Prb - 4*KeqAK*Prb + (SUMAXP**2 - 2*SUMAXP*Prb + 8*KeqAK*SUMAXP*Prb + Prb**2 - 4*KeqAK*Prb**2)**0.5)/((2 - 8*KeqAK)*KmPGKATP))*(1 + BPG/KmPGKBPG + P3G/KmPGKP3G))
v_9:
 {1.0}P3G = {1.0}P2G
  VmPGM/KmPGMP3G*(P3G - P2G/KeqPGM)/(1 + P3G/KmPGMP3G + P2G/KmPGMP2G)


#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
default_compartment = 1.0


#initial values
ACE = 0.04
BPG = 0.0
F16P = 0.1
F6P = 0.28
G6P = 1.39
GLCi = 0.087
NAD = 1.2
NADH = 0.39
P2G = 0.1
P3G = 0.1
PEP = 0.1
PYR = 3.36
Prb = 5.0
TRIO = 5.17


#assignment rules