#dupreez1



#reactions
v_ADH:
 {1.0}ACE + {1.0}NADH = {1.0}NAD
  -((VmADH*(ETOH*NAD - (ACE*NADH)/KeqADH))/(KiADHNAD*KmADHETOH*(1 + (ETOH*KmADHNAD)/(KiADHNAD*KmADHETOH) + (KmADHNADH*ACE)/(KiADHNADH*KmADHACE) + NAD/KiADHNAD + (ETOH*NAD)/(KiADHNAD*KmADHETOH) + (ETOH*ACE*NAD)/(KiADHACE*KiADHNAD*KmADHETOH) + (KmADHNADH*ACE*NAD)/(KiADHNAD*KiADHNADH*KmADHACE) + NADH/KiADHNADH + (ETOH*KmADHNAD*NADH)/(KiADHNAD*KiADHNADH*KmADHETOH) + (ACE*NADH)/(KiADHNADH*KmADHACE) + (ETOH*ACE*NADH)/(KiADHETOH*KiADHNADH*KmADHACE))))
v_AK:
 {2.0}ADP = {1.0}AMP + {1.0}ATP
  133.333*(ADP**2 - (AMP*ATP)/KeqAK)
v_ALD:
 {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_ATP:
 {1.0}ATP = {1.0}ADP
  (KATPASE*ATP**nATP)/(KmATP**nATP + ATP**nATP)
v_ENO:
 {1.0}P2G = {1.0}PEP
  (VmENO*(P2G - PEP/KeqENO))/(KmENOP2G*(1 + P2G/KmENOP2G + PEP/KmENOPEP))
v_G3PA:
 {1.0}G3P = $pool
  (VmG3PA*G3P)/(KmG3PAG3P*(1 + Phi/KmG3PAPhi)*(1 + G3P/KmG3PAG3P))
v_G3PDH:
 {1.0}NADH + {1.0}TRIO = {1.0}G3P + {1.0}NAD
  (VmG3PDH*(-((G3P*NAD)/KeqG3PDH) + (NADH*TRIO)/(1 + KeqTPI)))/(KmG3PDHDHAP*KmG3PDHNADH*(1 + ADP/KmG3PDHADP + ATP/KmG3PDHATP + F16P/KmG3PDHF16P)*(1 + NAD/KmG3PDHNAD + NADH/KmG3PDHNADH)*(1 + G3P/KmG3PDHG3P + TRIO/((1 + KeqTPI)*KmG3PDHDHAP)))
v_GAPDH:
 {1.0}NAD + {1.0}TRIO = {1.0}BPG + {1.0}NADH
  (-(VmGAPDHf * BPG * NADH / (KeqGAPDH * KmGAPDHGAP * KmGAPDHNAD)) + KeqTPI * VmGAPDHf * NAD * TRIO / ((1 + KeqTPI) * KmGAPDHGAP * KmGAPDHNAD)) / ((1 + NAD / KmGAPDHNAD + NADH / KmGAPDHNADH) * (1 + BPG / KmGAPDHBPG + KeqTPI * TRIO / ((1 + KeqTPI) * KmGAPDHGAP)))
v_GLK:
 {1.0}ATP + {1.0}GLCi = {1.0}ADP + {1.0}G6P
  (VmGLK*(-((ADP*G6P)/KeqGLK) + ATP*GLCi))/(KmGLKATP*KmGLKGLCi*(1 + ADP/KmGLKADP + ATP/KmGLKATP)*(1 + G6P/KmGLKG6P + GLCi/KmGLKGLCi))
v_GLT:
 $pool = {1.0}GLCi
  ((VmGLT*(GLCo - GLCi/KeqGLT))/(KmGLTGLCo*(1 + GLCo/KmGLTGLCo + GLCi/KmGLTGLCi + (alpha*GLCo*GLCi)/(KmGLTGLCi*KmGLTGLCo))))
v_GLYCO:
 {1.0}ATP + {1.0}G6P = {1.0}ADP
  KGLYCOGEN*ATP*G6P
v_PDC:
 {1.0}PYR = {1.0}ACE
  (VmPDC*PYR**nPDC)/(KmPDCPYR**nPDC*(1 + PYR**nPDC/KmPDCPYR**nPDC))
v_PFK:
 {1.0}ATP + {1.0}F6P = {1.0}ADP + {1.0}F16P
  (gR*VmPFK*ATP*F6P*(1 + ATP/KmPFKATP + F6P/KmPFKF6P + (gR*ATP*F6P)/(KmPFKATP*KmPFKF6P)))/(KmPFKATP*KmPFKF6P*((L0*(1 + (CPFKAMP*AMP)/KPFKAMP)**2*(1 + (CiPFKATP*ATP)/KiPFKATP)**2*(1 + (CPFKATP*ATP)/KmPFKATP)**2*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P)/KPFKF16BP)**2)/((1 + AMP/KPFKAMP)**2*(1 + ATP/KiPFKATP)**2*(1 + F26BP/KPFKF26BP + F16P/KPFKF16BP)**2) + (1 + ATP/KmPFKATP + F6P/KmPFKF6P + (gR*ATP*F6P)/(KmPFKATP*KmPFKF6P))**2))
v_PGI:
 {1.0}G6P = {1.0}F6P
  (VmPGI*(-(F6P/KeqPGI) + G6P))/(KmPGIG6P*(1 + F6P/KmPGIF6P + G6P/KmPGIG6P))
v_PGK:
 {1.0}ADP + {1.0}BPG = {1.0}ATP + {1.0}P3G
  (VmPGK*(KeqPGK*ADP*BPG - ATP*P3G))/(KmPGKATP*KmPGKP3G*(1 + ADP/KmPGKADP + ATP/KmPGKATP)*(1 + BPG/KmPGKBPG + P3G/KmPGKP3G))
v_PGM:
 {1.0}P3G = {1.0}P2G
  (VmPGM*(-(P2G/KeqPGM) + P3G))/(KmPGMP3G*(1 + P2G/KmPGMP2G + P3G/KmPGMP3G))
v_PYK:
 {1.0}ADP + {1.0}PEP = {1.0}ATP + {1.0}PYR
  (VmPYK*(ADP*PEP - (ATP*PYR)/KeqPYK))/(KmPYKADP*KmPYKPEP*(1 + ADP/KmPYKADP + ATP/KmPYKATP)*(1 + PEP/KmPYKPEP + PYR/KmPYKPYR))
v_SUC:
 {2.0}ACE + {3.0}NAD + {4.0}ATP = {3.0}NADH + {4.0}ADP
  KSUCC*ACE
v_Treha:
 {1.0}ATP + {2.0}G6P = {1.0}ADP
  KTREHALOSE*ATP*G6P


#parameters
CPFKAMP = 0.0845
CPFKATP = 3.0
CPFKF16BP = 0.397
CPFKF26BP = 0.0174
CiPFKATP = 100.0
ETOH = 50.0
EXTERNAL = 0.0
F26BP = 0.02
GLCo = 50.0
KATPASE = 93.1853
KGLYCOGEN = 2.376965061
KPFKAMP = 0.0995
KPFKF16BP = 0.111
KPFKF26BP = 0.000682
KSUCC = 21.4
KTREHALOSE = 0.9507860244
KeqADH = 6.9e-05
KeqAK = 0.45
KeqALD = 0.069
KeqENO = 6.7
KeqG3PDH = 4300.0
KeqGAPDH = 0.00562639062770364
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
KmATP = 0.263159506089159
KmENOP2G = 0.04
KmENOPEP = 0.5
KmG3PAG3P = 4.21421727097407
KmG3PAPhi = 0.798017147070638
KmG3PDHADP = 1.60537927993171
KmG3PDHATP = 0.568725492710765
KmG3PDHDHAP = 0.4
KmG3PDHF16P = 4.77110176421347
KmG3PDHG3P = 1.088994646912
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
Phi = 1.20470265921072
VmADH = 810.0
VmALD = 322.258
VmENO = 365.806
VmG3PA = 538.371
VmG3PDH = 477.424
VmGAPDHf = 1184.52
VmGLK = 226.452
VmGLT = 97.264
VmPDC = 174.194
VmPFK = 182.903
VmPGI = 339.677
VmPGK = 1306.45
VmPGM = 2525.81
VmPYK = 1088.71
alpha = 0.91
gR = 5.12
nATP = 1.0
nPDC = 1.9
default_compartment = 1.0


#initial values
ACE = 0.0320190578411128
ADP = 1.0419077376654
AMP = 0.169126793522056
ATP = 2.88896546881255
BPG = 0.000272378466415018
F16P = 5.12303060629701
F6P = 1.35316959717949
G3P = 1.51464838664408
G6P = 5.23300353550081
GLCi = 0.0264527982731908
NAD = 1.01975659448225
NADH = 0.570243405517748
P2G = 0.0219411145027765
P3G = 0.182294273856502
PEP = 0.0180777173496713
PYR = 2.26719823411809
TRIO = 2.65068578040418


#assignment rules