#penkler2aa



#reactions
vPFvALD:
 {1.0}f16bpPF = {1.0}dhapPF + {1.0}gapPF
  (VPFvALD*f16bpPF*(1 - (dhapPF*gapPF)/(KeqPFvALD*Vpf*f16bpPF)))/(Kf16bpPFvALD*(1 + dhapPF/(KdhapPFvALD*Vpf) + f16bpPF/(Kf16bpPFvALD*Vpf) + gapPF/(KgapPFvALD*Vpf) + (dhapPF*gapPF)/(KdhapPFvALD*KgapPFvALD*Vpf**2)))
vPFvATPASE:
 {1.0}atpPF = {1.0}adpPF
  (VPFvATPASE*atpPF**5)/(KmPFvATPASE**5*Vpf**4*(1 + atpPF**5/(KmPFvATPASE**5*Vpf**5)))
vPFvENO:
 {1.0}p2gPF = {1.0}pepPF
  (VPFvENO*p2gPF*(1 - pepPF/(KeqPFvENO*p2gPF)))/(Kp2gPFvENO*(1 + p2gPF/(Kp2gPFvENO*Vpf) + pepPF/(KpepPFvENO*Vpf)))
vPFvG3PDH:
 {1.0}nadhPF + {1.0}dhapPF = {1.0}g3pPF + {1.0}nadPF
  (VfPFvG3PDH*dhapPF*nadhPF*(1 - (g3pPF*nadPF)/(KeqPFvG3PDH*dhapPF*nadhPF)))/(KdhapPFvG3PDH*KnadhPFvG3PDH*Vpf*(1 + dhapPF/(KdhapPFvG3PDH*Vpf) + g3pPF/(Kg3pPFvG3PDH*Vpf))*(1 + nadhPF/(KnadhPFvG3PDH*Vpf) + nadPF/(KnadPFvG3PDH*Vpf)))
vPFvGAPDH:
 {1.0}gapPF + {1.0}nadPF = {1.0}nadhPF + {1.0}b13pgPF
  (-((VrPFvGAPDH*b13pgPF*nadhPF)/(Kb13pgPFvGAPDH*KnadhPFvGAPDH*Vpf)) + (VPFvGAPDH*gapPF*nadPF)/(KgapPFvGAPDH*KnadPFvGAPDH*Vpf))/((1 + b13pgPF/(Kb13pgPFvGAPDH*Vpf) + gapPF/(KgapPFvGAPDH*Vpf))*(1 + nadhPF/(KnadhPFvGAPDH*Vpf) + nadPF/(KnadPFvGAPDH*Vpf)))
vPFvGLCtr:
 {1.0}glcEX = {1.0}glcPF
  (Vpf*VPFvGLCtr*(glcEX/Vex - glcPF/Vpf))/(KmPFvGLCtr*(1 + glcEX/(KmPFvGLCtr*Vex) + glcPF/(KmPFvGLCtr*Vpf) + (alphaPFvGLCtr*glcEX*glcPF)/(KmPFvGLCtr**2*Vex*Vpf)))
vPFvGLYtr:
 {1.0}g3pPF = {1.0}glyEX
  kPFvGLYtr*Vpf*(g3pPF/Vpf - glyEX/Vex)
vPFvHK:
 {1.0}glcPF + {1.0}atpPF = {1.0}g6pPF + {1.0}adpPF
  (VPFvHK*atpPF*(1 - (adpPF*g6pPF)/(KeqPFvHK*atpPF*glcPF))*glcPF)/(KatpPFvHK*KglcPFvHK*Vpf*(1 + adpPF/(KadpPFvHK*Vpf) + atpPF/(KatpPFvHK*Vpf))*(1 + atpPF/(KiatpPFvHK*Vpf))*(1 + g6pPF/(Kg6pPFvHK*Vpf) + glcPF/(KglcPFvHK*Vpf)))
vPFvLACtr:
 {1.0}lacPF = {1.0}lacEX
  KdPFvLACtr*(1 - (hEX*Vpf*lacEX)/(hPF*Vex*lacPF))*lacPF + (VPFvLACtr*(1 - (hEX*Vpf*lacEX)/(hPF*Vex*lacPF))*lacPF)/(KlacPFvLACtr*(1 + lacEX/(KlacPFvLACtr*Vex) + lacPF/(KlacPFvLACtr*Vpf) + pyrEX/(KipyrPFvLACtr*Vex) + pyrPF/(KipyrPFvLACtr*Vpf)))
vPFvLDH:
 {1.0}nadhPF + {1.0}pyrPF = {1.0}lacPF + {1.0}nadPF
  (-((VrPFvLDH*lacPF*nadPF)/(KlacPFvLDH*KnadPFvLDH*Vpf)) + (VfPFvLDH*nadhPF*pyrPF)/(KnadhPFvLDH*KpyrPFvLDH*Vpf))/((1 + nadhPF/(KnadhPFvLDH*Vpf) + nadPF/(KnadPFvLDH*Vpf))*(1 + lacPF/(KlacPFvLDH*Vpf) + pyrPF/(KpyrPFvLDH*Vpf)))
vPFvPFK:
 {1.0}atpPF + {1.0}f6pPF = {1.0}adpPF + {1.0}f16bpPF
  (VPFvPFK*atpPF*f6pPF)/(KatpPFvPFK*Kf6pPFvPFK*Vpf*(1 + atpPF/(KatpPFvPFK*Vpf))*(1 + adpPF/(KadpPFvPFK*Vpf) + atpPF/(KatpPFvPFK*Vpf))*(1 + f16bpPF/(Kf16bpPFvPFK*Vpf) + f6pPF/(Kf6pPFvPFK*Vpf)))
vPFvPGI:
 {1.0}g6pPF = {1.0}f6pPF
  (VPFvPGI*(1 - f6pPF/(KeqPFvPGI*g6pPF))*g6pPF)/(Kg6pPFvPGI*(1 + f6pPF/(Kf6pPFvPGI*Vpf) + g6pPF/(Kg6pPFvPGI*Vpf)))
vPFvPGK:
 {1.0}adpPF + {1.0}b13pgPF = {1.0}atpPF + {1.0}p3gPF
  (KadpPFvPGK*Kb13pgPFvPGK*KeqPFvPGK*VrPFvPGK*atpPF*(-1 + (KeqPFvPGK*adpPF*b13pgPF)/(atpPF*p3gPF))*p3gPF)/(KatpPFvPGK**2*Kp3gPFvPGK**2*Vpf*(1 + adpPF/(KadpPFvPGK*Vpf) + atpPF/(KatpPFvPGK*Vpf))*(1 + b13pgPF/(Kb13pgPFvPGK*Vpf) + p3gPF/(Kp3gPFvPGK*Vpf)))
vPFvPGM:
 {1.0}p3gPF = {1.0}p2gPF
  (VPFvPGM*(1 - p2gPF/(KeqPFvPGM*p3gPF))*p3gPF)/(Kp3gPFvPGM*(1 + p2gPF/(Kp2gPFvPGM*Vpf) + p3gPF/(Kp3gPFvPGM*Vpf)))
vPFvPK:
 {1.0}adpPF + {1.0}pepPF = {1.0}atpPF + {1.0}pyrPF
  (VPFvPK*adpPF*pepPF)/(KadpPFvPK*KpepPFvPK*Vpf*(1 + adpPF/(KadpPFvPK*Vpf))*(1 + pepPF/(KpepPFvPK*Vpf))*(1 + LPFvPK*(1 + (adpPF/(KiadpPFvPK*Vpf) + atpPF/(KiatpPFvPK*Vpf))**h)*(1 + (pepPF/(KipepPFvPK*Vpf) + pyrPF/(KipyrPFvPK*Vpf))**h)))
vPFvPYRtr:
 {1.0}pyrPF = {1.0}pyrEX
  KdPFvPYRtr*(1 - (hEX*Vpf*pyrEX)/(hPF*Vex*pyrPF))*pyrPF + (VPFvPYRtr*(1 - (hEX*Vpf*pyrEX)/(hPF*Vex*pyrPF))*pyrPF)/(KpyrPFvPYRtr*(1 + lacEX/(KilacPFvPYRtr*Vex) + lacPF/(KilacPFvPYRtr*Vpf) + pyrEX/(KpyrPFvPYRtr*Vex) + pyrPF/(KpyrPFvPYRtr*Vpf)))
vPFvTPI:
 {1.0}dhapPF = {1.0}gapPF
  (VPFvTPI*dhapPF*(1 - gapPF/(KeqPFvTPI*dhapPF)))/(KdhapPFvTPI*(1 + dhapPF/(KdhapPFvTPI*Vpf) + gapPF/(KgapPFvTPI*Vpf) + pepPF/(KipepPFvTPI*Vpf)))


#parameters
ConcAdpPF = 0.0015
ConcAtpPF = 0.0015
ConcB13pgPF = 0.0001
ConcDhapPF = 0.0001
ConcF16bpPF = 0.0001
ConcF6pPF = 0.0001
ConcG6pPF = 0.0001
ConcGapPF = 0.0001
ConcGlcEX = 0.005
ConcGlcPF = 0.0001
ConcGly3pPF = 0.0001
ConcGlyEX = 3.79121e-11
ConcLacEX = 3.79121e-11
ConcLacPF = 0.0001
ConcNadPF = 0.0015
ConcNadhPF = 0.0015
ConcP2gPF = 0.0001
ConcP3gPF = 0.0001
ConcPepPF = 1e-05
ConcPyrEX = 3.79121e-11
ConcPyrPF = 0.0001
KadpPFvHK = 0.000846735
KadpPFvPFK = 0.00074176
KadpPFvPGK = 0.00015
KadpPFvPK = 0.000317
KatpPFvHK = 0.00069647
KatpPFvPFK = 0.0007862
KatpPFvPGK = 0.00077
Kb13pgPFvGAPDH = 2.359e-05
Kb13pgPFvPGK = 1.34e-05
KdPFvLACtr = 0.005
KdPFvPYRtr = 0.0007
KdhapPFvALD = 0.00011
KdhapPFvG3PDH = 0.00034
KdhapPFvTPI = 0.001954
KeqPFvALD = 9e-05
KeqPFvENO = 4.6
KeqPFvG3PDH = 32600.0
KeqPFvHK = 1310.0
KeqPFvPGI = 0.33
KeqPFvPGK = 3200.0
KeqPFvPGM = 0.19
KeqPFvTPI = 0.04545
Kf16bpPFvALD = 6.84e-05
Kf16bpPFvPFK = 0.003626
Kf6pPFvPFK = 0.00109454
Kf6pPFvPGI = 9.6651e-05
Kg3pPFvG3PDH = 0.00398
Kg6pPFvHK = 4.3e-05
Kg6pPFvPGI = 0.00100774
KgapPFvALD = 4.6e-05
KgapPFvGAPDH = 0.000917
KgapPFvTPI = 0.000337
KglcPFvHK = 0.000168613
KiadpPFvPK = 0.002
KiatpPFvHK = 0.026
KiatpPFvPK = 0.00182
KilacPFvPYRtr = 0.000358
KipepPFvPK = 0.00292
KipepPFvTPI = 1.59e-05
KipyrPFvLACtr = 0.00163
KipyrPFvPK = 0.105
KlacPFvLACtr = 0.0038
KlacPFvLDH = 0.003611
KmPFvATPASE = 0.0045
KmPFvGLCtr = 0.000213
KnadPFvG3PDH = 0.00051
KnadPFvGAPDH = 0.0005662
KnadPFvLDH = 0.000234
KnadhPFvG3PDH = 9e-05
KnadhPFvGAPDH = 2.87e-05
KnadhPFvLDH = 4.6e-05
Kp2gPFvENO = 0.000521
Kp2gPFvPGM = 0.000318
Kp3gPFvPGK = 0.000267
Kp3gPFvPGM = 0.00173
KpepPFvENO = 0.00129
KpepPFvPK = 0.000406
KpyrPFvLDH = 0.000133
KpyrPFvPYRtr = 0.0157
LPFvPK = 0.255
VPFvALD = 0.0569593147718
VPFvATPASE = 0.345
VPFvENO = 0.286295503195
VPFvGAPDH = 0.656263383259
VPFvGLCtr = 0.0468
VPFvHK = 0.083
VPFvLACtr = 0.598
VPFvPFK = 0.41
VPFvPGI = 0.800428265477
VPFvPGM = 0.186937901488
VPFvPK = 0.762
VPFvPYRtr = 0.216
VPFvTPI = 1.50599571726
Vex = 2637.68
VfPFvG3PDH = 0.012
VfPFvLDH = 0.542612419668
Vpf = 28.0
VrPFvGAPDH = 0.243199143455
VrPFvLDH = 0.290792291203
VrPFvPGK = 0.353533190557
alphaPFvGLCtr = 0.91
h = 4.0
kPFvGLYtr = 2.0
pHConversionFactor = 1.0
phEX = 7.1
phPF = 7.2
default = 1.0


#initial values
adpPF = 0.042
atpPF = 0.042
b13pgPF = 0.0028
dhapPF = 0.0028
f16bpPF = 0.0028
f6pPF = 0.0028
g3pPF = 0.0028
g6pPF = 0.0028
gapPF = 0.0028
glcEX = 13.1884
glcPF = 0.0028
glyEX = 0.0001
lacEX = 0.0001
lacPF = 0.0028
nadPF = 0.042
nadhPF = 0.042
p2gPF = 0.0028
p3gPF = 0.0028
pepPF = 0.00028
pyrEX = 0.0001
pyrPF = 0.0028


#assignment rules
!F hPF = pHConversionFactor/10.0^phPF
!F hEX = pHConversionFactor/10.0^phEX