#smallbone16

FIX: ACE EtOH GLCx GLY SUC TRH

#reactions
ADHADH1:
 {1.0}AcAld + {1.0}NADH = {1.0}NAD + {1.0}EtOH
  (ADH1*ADHADH1kcat*(-((EtOH*NAD)/(ADHADH1Kacald*ADHADH1Kinadh*KeqADH)) + (AcAld*NADH)/(ADHADH1Kacald*ADHADH1Kinadh)))/(1 + (ADHADH1Knad*EtOH)/(ADHADH1Ketoh*ADHADH1Kinad) + (ADHADH1Knadh*AcAld)/(ADHADH1Kacald*ADHADH1Kinadh) + NAD/ADHADH1Kinad + (EtOH*NAD)/(ADHADH1Ketoh*ADHADH1Kinad) + (ADHADH1Knadh*AcAld*NAD)/(ADHADH1Kacald*ADHADH1Kinad*ADHADH1Kinadh) + (EtOH*AcAld*NAD)/(ADHADH1Ketoh*ADHADH1Kiacald*ADHADH1Kinad) + NADH/ADHADH1Kinadh + (ADHADH1Knad*EtOH*NADH)/(ADHADH1Ketoh*ADHADH1Kinad*ADHADH1Kinadh) + (AcAld*NADH)/(ADHADH1Kacald*ADHADH1Kinadh) + (EtOH*AcAld*NADH)/(ADHADH1Kacald*ADHADH1Kietoh*ADHADH1Kinadh))
ADHADH5:
 {1.0}AcAld + {1.0}NADH = {1.0}NAD + {1.0}EtOH
  (ADH5*ADHADH5kcat*(-((EtOH*NAD)/(ADHADH5Kacald*ADHADH5Kinadh*KeqADH)) + (AcAld*NADH)/(ADHADH5Kacald*ADHADH5Kinadh)))/(1 + (ADHADH5Knad*EtOH)/(ADHADH5Ketoh*ADHADH5Kinad) + (ADHADH5Knadh*AcAld)/(ADHADH5Kacald*ADHADH5Kinadh) + NAD/ADHADH5Kinad + (EtOH*NAD)/(ADHADH5Ketoh*ADHADH5Kinad) + (ADHADH5Knadh*AcAld*NAD)/(ADHADH5Kacald*ADHADH5Kinad*ADHADH5Kinadh) + (EtOH*AcAld*NAD)/(ADHADH5Ketoh*ADHADH5Kiacald*ADHADH5Kinad) + NADH/ADHADH5Kinadh + (ADHADH5Knad*EtOH*NADH)/(ADHADH5Ketoh*ADHADH5Kinad*ADHADH5Kinadh) + (AcAld*NADH)/(ADHADH5Kacald*ADHADH5Kinadh) + (EtOH*AcAld*NADH)/(ADHADH5Kacald*ADHADH5Kietoh*ADHADH5Kinadh))
AK:
 {2.0}ADP = {1.0}AMP + {1.0}ATP
  AKk*(ADP**2 - (AMP*ATP)/AKKeq)
ATPase:
 {1.0}ATP = {1.0}ADP
  (ATPaseVmax*ATP)/(ATPaseKatp*(1 + ATP/ATPaseKatp))
ENOENO1:
 {1.0}P2G = {1.0}PEP
  (ENO1*ENOENO1kcat*(P2G/ENOENO1Kp2g - PEP/(ENOENO1Kp2g*KeqENO)))/(1 + P2G/ENOENO1Kp2g + PEP/ENOENO1Kpep)
ENOENO2:
 {1.0}P2G = {1.0}PEP
  (ENO2*ENOENO2kcat*(P2G/ENOENO2Kp2g - PEP/(ENOENO2Kp2g*KeqENO)))/(1 + P2G/ENOENO2Kp2g + PEP/ENOENO2Kpep)
FBA:
 {1.0}F16bP = {1.0}DHAP + {1.0}GAP
  (FBA1*FBAkcat*(F16bP/FBAKf16bp - (DHAP*GAP)/(FBAKeq*FBAKf16bp)))/(1 + DHAP/FBAKdhap + F16bP/FBAKf16bp + GAP/FBAKgap + (DHAP*GAP)/(FBAKdhap*FBAKgap) + (F16bP*GAP)/(FBAKf16bp*FBAKigap))
GPD:
 {1.0}DHAP + {1.0}NADH = {1.0}G3P + {1.0}NAD
  GPDVmax * (-(G3P * NAD / GPDKeq) + DHAP * NADH) / (GPDKdhap * GPDKnadh * (1 + ADP / GPDKadp + ATP / GPDKatp + F16bP / GPDKfbp) * (1 + DHAP / GPDKdhap + G3P / GPDKg3p) * (1 + NAD / GPDKnad + NADH / GPDKnadh))
GPM:
 {1.0}P3G = {1.0}P2G
  (GPM1*GPMkcat*(-(P2G/(GPMKeq*GPMKp3g)) + P3G/GPMKp3g))/(1 + P2G/GPMKp2g + P3G/GPMKp3g)
GPP:
 {1.0}G3P = {1.0}GLY
  (GPPVmax*G3P)/(GPPKg3p*(1 + G3P/GPPKg3p))
HXKGLK1:
 {1.0}ATP + {1.0}GLC = {1.0}ADP + {1.0}G6P
  (GLK1*HXKGLK1kcat*(-((ADP*G6P)/(HXKGLK1Katp*HXKGLK1Kglc*KeqHXK)) + (ATP*GLC)/(HXKGLK1Katp*HXKGLK1Kglc)))/((1 + ADP/HXKGLK1Kadp + ATP/HXKGLK1Katp)*(1 + G6P/HXKGLK1Kg6p + GLC/HXKGLK1Kglc))
HXKHXK1:
 {1.0}ATP + {1.0}GLC = {1.0}ADP + {1.0}G6P
  HXK1 * HXKHXK1kcat * (-(ADP * G6P / (HXKHXK1Katp * HXKHXK1Kglc * KeqHXK)) + ATP * GLC / (HXKHXK1Katp * HXKHXK1Kglc)) / ((1 + ADP / HXKHXK1Kadp + ATP / HXKHXK1Katp) * (1 + G6P / HXKHXK1Kg6p + GLC / HXKHXK1Kglc + T6P / HXKHXK1Kit6p))
HXKHXK2:
 {1.0}ATP + {1.0}GLC = {1.0}ADP + {1.0}G6P
  HXK2 * HXKHXK2kcat * (-(ADP * G6P / (HXKHXK2Katp * HXKHXK2Kglc * KeqHXK)) + ATP * GLC / (HXKHXK2Katp * HXKHXK2Kglc)) / ((1 + ADP / HXKHXK2Kadp + ATP / HXKHXK2Katp) * (1 + G6P / HXKHXK2Kg6p + GLC / HXKHXK2Kglc + T6P / HXKHXK2Kit6p))
HXT:
 {1.0}GLCx = {1.0}GLC
  (HXTVmax*(GLCx - GLC))/(HXTKglc*(1 + GLCx/HXTKglc + GLC/HXTKglc + (GLCx*HXTKi*GLC)/HXTKglc**2))
PDCPDC1:
 {1.0}PYR = {1.0}AcAld
  (PDC1*PDCPDC1kcat*PYR)/(PDCPDC1Kpyr*(1 + PYR/PDCPDC1Kpyr))
PDCPDC5:
 {1.0}PYR = {1.0}AcAld
  (PDC5*PDCPDC5kcat*PYR)/(PDCPDC5Kpyr*(1 + PYR/PDCPDC5Kpyr))
PDCPDC6:
 {1.0}PYR = {1.0}AcAld
  (PDC6*PDCPDC6kcat*PYR)/(PDCPDC6Kpyr*(1 + PYR/PDCPDC6Kpyr))
PFK:
 {1.0}ATP + {1.0}F6P = {1.0}ADP + {1.0}F16bP
  PFKgR * PFKkcat * (PFK1 + PFK2 - Abs(PFK1 - PFK2)) * ATP * (1 - ADP * F16bP / (PFKKeq * ATP * F6P)) * F6P * (1 + ADP / PFKKadp + ATP / PFKKatp + F16bP / PFKKf16 + PFKgR * ADP * F16bP / (PFKKadp * PFKKf16) + F6P / PFKKf6p + PFKgR * ATP * F6P / (PFKKatp * PFKKf6p)) / (2 * PFKKatp * PFKKf6p * (PFKL0 * (1 + PFKCamp * AMP / PFKKamp)**2 * (1 + PFKCatp * ATP / PFKKatp)**2 * (1 + PFKCiatp * ATP / PFKKiatp)**2 * (1 + F26bP * PFKCf26 / PFKKf26 + PFKCf16 * F16bP / PFKKf16)**2 / ((1 + AMP / PFKKamp)**2 * (1 + ATP / PFKKiatp)**2 * (1 + F26bP / PFKKf26 + F16bP / PFKKf16)**2) + (1 + ADP / PFKKadp + ATP / PFKKatp + F16bP / PFKKf16 + PFKgR * ADP * F16bP / (PFKKadp * PFKKf16) + F6P / PFKKf6p + PFKgR * ATP * F6P / (PFKKatp * PFKKf6p))**2))
PGI:
 {1.0}G6P = {1.0}F6P
  (PGI1*PGIkcat*(-(F6P/(PGIKeq*PGIKg6p)) + G6P/PGIKg6p))/(1 + F6P/PGIKf6p + G6P/PGIKg6p)
PGK:
 {1.0}ADP + {1.0}BPG = {1.0}ATP + {1.0}P3G
  (PGK1*PGKkcat*(ADP/PGKKadp)**(-1 + PGKnHadp)*((ADP*BPG)/(PGKKadp*PGKKbpg) - (ATP*P3G)/(PGKKadp*PGKKbpg*PGKKeq)))/((1 + (ADP/PGKKadp)**PGKnHadp + ATP/PGKKatp)*(1 + BPG/PGKKbpg + P3G/PGKKp3g))
PGM:
 {1.0}G6P = {1.0}G1P
  (PGMVmax*rescaleTRH*(-(G1P/(PGMKeq*PGMKg6p)) + G6P/PGMKg6p))/(1 + G1P/PGMKg1p + G6P/PGMKg6p)
PYKCDC19:
 {1.0}ADP + {1.0}PEP = {1.0}ATP + {1.0}PYR
  CDC19 * PYKCDC19kcat * (ADP * PEP - ATP * PYR / KeqPYK) / (PYKCDC19Kadp * PYKCDC19Kpep * (1 + ADP / PYKCDC19Kadp + ATP / PYKCDC19Katp) * (1 + PYKCDC19L0 * (1 + ATP / PYKCDC19Kiatp) / (1 + F16bP / PYKCDC19Kf16p) + PEP / PYKCDC19Kpep + PYR / PYKCDC19Kpyr))
PYKPYK2:
 {1.0}ADP + {1.0}PEP = {1.0}ATP + {1.0}PYR
  PYK2 * PYKPYK2kcat * (ADP * PEP - ATP * PYR / KeqPYK) / (PYKPYK2Kadp * PYKPYK2Kpep * (1 + ADP / PYKPYK2Kadp + ATP / PYKPYK2Katp) * (1 + PYKPYK2L0 * (1 + ATP / PYKPYK2Kiatp) / (1 + F16bP / PYKPYK2Kf16p) + PEP / PYKPYK2Kpep + PYR / PYKPYK2Kpyr))
TDHTDH1:
 {1.0}GAP + {1.0}NAD = {1.0}BPG + {1.0}NADH
  (TDH1*TDHTDH1kcat*((GAP*NAD)/(TDHTDH1Kgap*TDHTDH1Knad) - (BPG*NADH)/(KeqTDH*TDHTDH1Kgap*TDHTDH1Knad)))/((1 + BPG/TDHTDH1Kbpg + GAP/TDHTDH1Kgap)*(1 + NAD/TDHTDH1Knad + NADH/TDHTDH1Knadh))
TDHTDH2:
 {1.0}GAP + {1.0}NAD = {1.0}BPG + {1.0}NADH
  (TDH2*TDHTDH2kcat*((GAP*NAD)/(TDHTDH2Kgap*TDHTDH2Knad) - (BPG*NADH)/(KeqTDH*TDHTDH2Kgap*TDHTDH2Knad)))/((1 + BPG/TDHTDH2Kbpg + GAP/TDHTDH2Kgap)*(1 + NAD/TDHTDH2Knad + NADH/TDHTDH2Knadh))
TDHTDH3:
 {1.0}GAP + {1.0}NAD = {1.0}BPG + {1.0}NADH
  (TDH3*TDHTDH3kcat*((GAP*NAD)/(TDHTDH3Kgap*TDHTDH3Knad) - (BPG*NADH)/(KeqTDH*TDHTDH3Kgap*TDHTDH3Knad)))/((1 + BPG/TDHTDH3Kbpg + GAP/TDHTDH3Kgap)*(1 + NAD/TDHTDH3Knad + NADH/TDHTDH3Knadh))
TPI:
 {1.0}DHAP = {1.0}GAP
  TPIk*(DHAP - GAP/TPIKeq)
TPP:
 {1.0}T6P = {1.0}TRH
  (rescaleTRH*TPPVmax*T6P)/(TPPKt6p*(1 + T6P/TPPKt6p))
TPS:
 {1.0}G6P + {1.0}UDG = {1.0}T6P + {1.0}UDP
  (rescaleTRH*TPSVmax*G6P*UDG)/(TPSKg6p*TPSKudg*(1 + G6P/TPSKg6p)*(1 + UDG/TPSKudg))
UGP:
 {1.0}G1P + {1.0}UTP = {1.0}UDG
  (rescaleTRH*UGPVmax*G1P*UTP)/(UGPKg1p*UGPKutp*(UGPKiutp/UGPKutp + G1P/UGPKg1p + (UGPKiutp*UDG)/(UGPKiudg*UGPKutp) + (G1P*UDG)/(UGPKg1p*UGPKiudg) + UTP/UGPKutp + (G1P*UTP)/(UGPKg1p*UGPKutp)))
acetatebranch:
 {1.0}AcAld + {1.0}NAD = {1.0}NADH + {1.0}ACE
  acetatebranchk*AcAld*NAD
succinatebranch:
 {1.0}PYR + {3.0}NAD = {0.75}SUC + {3.0}NADH
  succinatebranchk*NAD*PYR
udptoutp:
 {1.0}ATP + {1.0}UDP = {1.0}ADP + {1.0}UTP
  rescaleTRH*udptoutpk*ATP*UDP


#parameters
ADH1 = 0.0409771277320022
ADH5 = 0.00106249605621922
ADHADH1Kacald = 0.4622
ADHADH1Ketoh = 17.0
ADHADH1Kiacald = 1.1
ADHADH1Kietoh = 90.0
ADHADH1Kinad = 0.92
ADHADH1Kinadh = 0.031
ADHADH1Knad = 0.17
ADHADH1Knadh = 0.11
ADHADH1kcat = 176.0
ADHADH5Kacald = 1.11
ADHADH5Ketoh = 17.0
ADHADH5Kiacald = 1.1
ADHADH5Kietoh = 90.0
ADHADH5Kinad = 0.92
ADHADH5Kinadh = 0.031
ADHADH5Knad = 0.17
ADHADH5Knadh = 0.11
ADHADH5kcat = 0.0
AKKeq = 0.45
AKk = 0.75
ATPaseKatp = 3.0
ATPaseVmax = 1.13333333333333
CDC19 = 0.512097526792801
ENO1 = 0.171592988538958
ENO2 = 0.493611573294543
ENOENO1Kp2g = 0.043
ENOENO1Kpep = 0.5
ENOENO1kcat = 7.6
ENOENO2Kp2g = 0.104
ENOENO2Kpep = 0.5
ENOENO2kcat = 19.87
EXTERNAL = 0.0
F26bP = 0.003
FBA1 = 0.334598664262206
FBAKdhap = 2.0
FBAKeq = 0.069
FBAKf16bp = 0.4507
FBAKgap = 2.4
FBAKigap = 10.0
FBAkcat = 4.139
GLK1 = 0.0112717406104807
GPD1 = 0.00170877794272468
GPD2 = 0.000198351416606057
GPDKadp = 2.0
GPDKatp = 0.73
GPDKdhap = 0.54
GPDKeq = 10000.0
GPDKfbp = 4.8
GPDKg3p = 1.2
GPDKnad = 0.93
GPDKnadh = 0.023
GPDVmax = 0.783333333333333
GPM1 = 0.182500074724267
GPMKeq = 0.19
GPMKp2g = 1.41
GPMKp3g = 1.2
GPMkcat = 400.0
GPPKg3p = 3.5
GPPVmax = 0.883333333333333
HOR2 = 0.0013683673909939
HXK1 = 0.0041951864287446
HXK2 = 0.0153328384926289
HXKGLK1Kadp = 0.23
HXKGLK1Katp = 0.865
HXKGLK1Kg6p = 30.0
HXKGLK1Kglc = 0.0106
HXKGLK1kcat = 0.0721
HXKHXK1Kadp = 0.23
HXKHXK1Katp = 0.293
HXKHXK1Kg6p = 30.0
HXKHXK1Kglc = 0.15
HXKHXK1Kit6p = 0.2
HXKHXK1kcat = 10.2
HXKHXK2Kadp = 0.23
HXKHXK2Katp = 0.195
HXKHXK2Kg6p = 30.0
HXKHXK2Kglc = 0.2
HXKHXK2Kit6p = 0.04
HXKHXK2kcat = 63.1
HXTKglc = 0.9
HXTKi = 0.91
HXTVmax = 3.35
KeqADH = 14492.7536231884
KeqENO = 6.7
KeqHXK = 2000.0
KeqPYK = 6500.0
KeqTDH = 0.00533412710224736
NA = 6.02214e+20
PDC1 = 0.266952694557084
PDC5 = 0.00308868608169189
PDC6 = 0.0016352160527653
PDCPDC1Kpyr = 8.5
PDCPDC1kcat = 12.14
PDCPDC5Kpyr = 7.08
PDCPDC5kcat = 10.32
PDCPDC6Kpyr = 2.92
PDCPDC6kcat = 9.21
PFK1 = 0.0116962574765781
PFK2 = 0.00975915538330228
PFKCamp = 0.0845
PFKCatp = 3.0
PFKCf16 = 0.397
PFKCf26 = 0.0174
PFKCiatp = 100.0
PFKKadp = 1.0
PFKKamp = 0.0995
PFKKatp = 0.71
PFKKeq = 800.0
PFKKf16 = 0.111
PFKKf26 = 0.000682
PFKKf6p = 0.1
PFKKiatp = 0.65
PFKL0 = 0.66
PFKgR = 5.12
PFKkcat = 209.6
PGI1 = 0.0345726768225249
PGIKeq = 0.29
PGIKf6p = 0.307
PGIKg6p = 1.0257
PGIkcat = 487.36
PGK1 = 0.0644142281647388
PGKKadp = 0.2
PGKKatp = 1.99
PGKKbpg = 0.003
PGKKeq = 3200.0
PGKKp3g = 4.58
PGKkcat = 58.6
PGKnHadp = 2.0
PGM1 = 0.000815573865768647
PGM2 = 0.00031467219294138
PGMKeq = 0.1667
PGMKg1p = 0.023
PGMKg6p = 0.05
PGMVmax = 0.00590833333333333
PYK2 = 0.00151748381804475
PYKCDC19Kadp = 0.243
PYKCDC19Katp = 1.5
PYKCDC19Kf16p = 0.2
PYKCDC19Kiatp = 9.3
PYKCDC19Kpep = 0.281
PYKCDC19Kpyr = 21.0
PYKCDC19L0 = 100.0
PYKCDC19kcat = 20.146
PYKPYK2Kadp = 0.3
PYKPYK2Katp = 1.5
PYKPYK2Kf16p = 0.2
PYKPYK2Kiatp = 9.3
PYKPYK2Kpep = 0.19
PYKPYK2Kpyr = 21.0
PYKPYK2L0 = 100.0
PYKPYK2kcat = 0.0
RHR2 = 0.0127951193429578
TDH1 = 0.087716160700349
TDH2 = 0.0
TDH3 = 1.05110118662137
TDHTDH1Kbpg = 0.0098
TDHTDH1Kgap = 0.495
TDHTDH1Knad = 0.09
TDHTDH1Knadh = 0.06
TDHTDH1kcat = 19.12
TDHTDH2Kbpg = 0.0098
TDHTDH2Kgap = 0.77
TDHTDH2Knad = 0.09
TDHTDH2Knadh = 0.06
TDHTDH2kcat = 8.633
TDHTDH3Kbpg = 0.909
TDHTDH3Kgap = 0.423
TDHTDH3Knad = 0.09
TDHTDH3Knadh = 0.06
TDHTDH3kcat = 18.162
TPI1 = 0.073589454911377
TPIKeq = 0.045
TPIk = 7500.0
TPPKt6p = 0.5
TPPVmax = 0.108333333333333
TPS1 = 0.000848120435592663
TPS2 = 0.000664962953368736
TPSKg6p = 3.8
TPSKudg = 0.886
TPSVmax = 0.02285
UGP1 = 0.00155052854965178
UGPKg1p = 0.32
UGPKiudg = 0.0035
UGPKiutp = 0.11
UGPKutp = 0.11
UGPVmax = 0.613666666666667
acetatebranchk = 0.00554339592436782
cell = 1.0
extracellular = 1.0
rescaleTRH = 22.0
succinatebranchk = 0.00549663630823988
sumAXP = 6.02
sumNAD = 1.59
sumUXP = 1.39784619487425
udptoutpk = 0.00345026988011002
volume = 2e-14
ACE = 55.7500633495734
EtOH = 55.472577854384
GLCx = 74.0
GLY = 0.15
SUC = 0.0
TRH = 0.00384697134241316
default_compartment = 1.0


#initial values
ADP = 1.29
AMP = 0.44
ATP = 4.29
AcAld = 0.178140579850657
BPG = 0.000736873499865602
DHAP = 0.290344213186674
F16bP = 1.14580464751733
F6P = 0.0588603054728053
G1P = 0.13481212658623
G3P = 0.068500732297821
G6P = 0.193120800911304
GAP = 0.0789727571926259
GLC = 1.5700004483456
NAD = 1.50329030201531
NADH = 0.0867096979846952
P2G = 0.0169344770274836
P3G = 0.117456252783611
PEP = 0.152501353339511
PYR = 0.527117851793548
T6P = 0.02
UDG = 0.467246194874247
UDP = 0.2815
UTP = 0.6491


#assignment rules
!F energycharge = (ADP/2 + ATP)/sumAXP
!F sumPXG = P2G + P3G
!F fitconc = Sqrt((1 - (NA*sumPXG*volume)/1618640)^2 + (1 - (NA*volume*DHAP)/3496987)^2 + (1 - (NA*volume*F16bP)/13800392)^2 + (1 - (NA*volume*F6P)/708930)^2 + (1 - (NA*volume*G6P)/2326001)^2 + (1 - (NA*volume*GAP)/951170)^2 + (1 - (NA*volume*GLC)/18909525)^2 + (1 - (NA*volume*PEP)/1836769)^2 + (1 - (NA*volume*PYR)/6348755)^2)/3