#kerkhovenAA



#reactions
AKc:
 {0.366643}ADPc = {0.183321}AMPc + {0.183321}ATPc
  AKcAKck1*ADPc**2 - AKcAKck2*AMPc*ATPc
AKg:
 {8.15993}ADPg = {4.07997}AMPg + {4.07997}ATPg
  AKgAKgk1*ADPg**2 - AKgAKgk2*AMPg*ATPg
ALDg:
 {4.07997}Fru16BPg = {4.07997}DHAPg + {4.07997}GA3Pg
  (ALDgALDgVmax*Fru16BPg*(1 - (DHAPg*GA3Pg)/(ALDgALDgKeq*Fru16BPg)))/(ALDgALDgKmFru16BP*(1 + ADPg/ALDgALDgKiADP + AMPg/ALDgALDgKiAMP + ATPg/ALDgALDgKiATP)*(1 + DHAPg/ALDgALDgKmDHAP + Fru16BPg/(ALDgALDgKmFru16BP*(1 + ADPg/ALDgALDgKiADP + AMPg/ALDgALDgKiAMP + ATPg/ALDgALDgKiATP)) + GA3Pg/ALDgALDgKmGA3P + (DHAPg*GA3Pg)/(ALDgALDgKmDHAP*ALDgALDgKmGA3P) + (Fru16BPg*GA3Pg)/(ALDgALDgKiGA3P*ALDgALDgKmFru16BP*(1 + ADPg/ALDgALDgKiADP + AMPg/ALDgALDgKiAMP + ATPg/ALDgALDgKiATP))))
ATPuc:
 {0.183321}ATPc = {0.183321}ADPc
  (ATPucATPuck*ATPc)/ADPc
ENOc:
 {0.183321}P2GAc = {0.183321}PEPc
  (ENOcENOcVmax*(1 - PEPc/(ENOcENOcKeq*P2GAc))*P2GAc)/(ENOcENOcKm2PGA*(1 + PEPc/ENOcENOcKmPEP + P2GAc/ENOcENOcKm2PGA))
G3PDHg:
 {4.07997}DHAPg + {4.07997}NADHg = {4.07997}Gly3Pg + {4.07997}NADg
  (G3PDHgG3PDHgVmax*DHAPg*NADHg*(1 - (Gly3Pg*NADg)/(G3PDHgG3PDHgKeq*DHAPg*NADHg)))/(G3PDHgG3PDHgKmDHAP*G3PDHgG3PDHgKmNADH*(1 + DHAPg/G3PDHgG3PDHgKmDHAP + Gly3Pg/G3PDHgG3PDHgKmGly3P)*(1 + NADHg/G3PDHgG3PDHgKmNADH + NADg/G3PDHgG3PDHgKmNAD))
GAPDHg:
 {4.07997}GA3Pg + {4.07997}NADg = {4.07997}NADHg + {4.07997}B13PGAg
  (GAPDHgGAPDHgVmax*GA3Pg*NADg*(1 - (NADHg*B13PGAg)/(GAPDHgGAPDHgKeq*GA3Pg*NADg)))/(GAPDHgGAPDHgKmGA3P*GAPDHgGAPDHgKmNAD*(1 + NADHg/GAPDHgGAPDHgKmNADH + NADg/GAPDHgGAPDHgKmNAD)*(1 + GA3Pg/GAPDHgGAPDHgKmGA3P + B13PGAg/GAPDHgGAPDHgKm13BPGA))
GDAg:
 {0.183321}DHAPc + {4.07997}Gly3Pg = {4.07997}DHAPg + {0.183321}Gly3Pc
  -(GDAgGDAgk*DHAPg*Gly3Pc) + GDAgGDAgk*DHAPc*Gly3Pg
GKg:
 {4.07997}ADPg + {4.07997}Gly3Pg = {4.07997}ATPg
  (GKgGKgVmax*ADPg*(1 - (Glye*ATPg)/(GKgGKgKeq*ADPg*Gly3Pg))*Gly3Pg)/(GKgGKgKmADP*GKgGKgKmGly3P*(1 + ADPg/GKgGKgKmADP + ATPg/GKgGKgKmATP)*(1 + Glye/GKgGKgKmGly + Gly3Pg/GKgGKgKmGly3P))
GPOc:
 {0.183321}Gly3Pc = {0.183321}DHAPc
  (GPOcGPOcVmax*Gly3Pc)/(GPOcGPOcKmGly3P*(1 + Gly3Pc/GPOcGPOcKmGly3P))
GlcTc:
 $pool = {0.183321}Glcc
  (GlcTcGlcTcVmax*(Glce - Glcc))/(GlcTcGlcTcKmGlc + Glce + Glcc + (GlcTcGlcTcalpha*Glce*Glcc)/GlcTcGlcTcKmGlc)
GlcTg:
 {0.183321}Glcc = {4.07997}Glcg
  GlcTgGlcTgk1*Glcc - GlcTgGlcTgk2*Glcg
HXKg:
 {4.07997}ATPg + {4.07997}Glcg = {4.07997}ADPg + {4.07997}Glc6Pg
  (HXKgHXKgVmax*ATPg*(1 - (ADPg*Glc6Pg)/(HXKgHXKgKeq*ATPg*Glcg))*Glcg)/(HXKgHXKgKmATP*HXKgHXKgKmGlc*(1 + ADPg/HXKgHXKgKmADP + ATPg/HXKgHXKgKmATP)*(1 + Glc6Pg/HXKgHXKgKmGlc6P + Glcg/HXKgHXKgKmGlc))
P3GATg:
 {4.07997}P3GAg = {0.183321}P3GAc
  -(P3GATgP3GATgk*P3GAc) + P3GATgP3GATgk*P3GAg
PFKg:
 {4.07997}ATPg + {4.07997}Fru6Pg = {4.07997}ADPg + {4.07997}Fru16BPg
  (PFKgPFKgKi1*PFKgPFKgVmax*ATPg*(1 - (ADPg*Fru16BPg)/(PFKgPFKgKeq*ATPg*Fru6Pg))*Fru6Pg)/(PFKgPFKgKmATP*PFKgPFKgKmFru6P*(PFKgPFKgKi1 + Fru16BPg)*(PFKgPFKgKsATP/PFKgPFKgKmATP + ADPg/PFKgPFKgKmADP + ATPg/PFKgPFKgKmATP + (ADPg*Fru16BPg)/(PFKgPFKgKi2*PFKgPFKgKmADP) + Fru6Pg/PFKgPFKgKmFru6P + (ATPg*Fru6Pg)/(PFKgPFKgKmATP*PFKgPFKgKmFru6P)))
PGAMc:
 {0.183321}P3GAc = {0.183321}P2GAc
  (PGAMcPGAMcVmax*(1 - P2GAc/(PGAMcPGAMcKeq*P3GAc))*P3GAc)/(PGAMcPGAMcKm3PGA*(1 + P2GAc/PGAMcPGAMcKm2PGA + P3GAc/PGAMcPGAMcKm3PGA))
PGIg:
 {4.07997}Glc6Pg = {4.07997}Fru6Pg
  (PGIgPGIgVmax*(1 - Fru6Pg/(PGIgPGIgKeq*Glc6Pg))*Glc6Pg)/(PGIgPGIgKmGlc6P*(1 + PGIg6PGg/PGIgPGIgKi6PG + Fru6Pg/PGIgPGIgKmFru6P + Glc6Pg/PGIgPGIgKmGlc6P))
PGKg:
 {4.07997}ADPg + {4.07997}B13PGAg = {4.07997}ATPg + {4.07997}P3GAg
  (PGKgPGKgVmax*ADPg*B13PGAg*(1 - (ATPg*P3GAg)/(PGKgPGKgKeq*ADPg*B13PGAg)))/(PGKgPGKgKm13BPGA*PGKgPGKgKmADP*(1 + ADPg/PGKgPGKgKmADP + ATPg/PGKgPGKgKmATP)*(1 + B13PGAg/PGKgPGKgKm13BPGA + P3GAg/PGKgPGKgKm3PGA))
PYKc:
 {0.183321}ADPc + {0.183321}PEPc = {0.183321}ATPc + {0.183321}Pyrc
  (PYKcPYKcVmax*ADPc*(PEPc/(PYKcPYKcKmPEP*(1 + ADPc/PYKcPYKcKiADP + ATPc/PYKcPYKcKiATP)))**PYKcPYKcn*(1 - (ATPc*Pyrc)/(PYKcPYKcKeq*ADPc*PEPc)))/(PYKcPYKcKmADP*(1 + ADPc/PYKcPYKcKmADP + ATPc/PYKcPYKcKmATP)*(1 + (PEPc/(PYKcPYKcKmPEP*(1 + ADPc/PYKcPYKcKiADP + ATPc/PYKcPYKcKiATP)))**PYKcPYKcn + Pyrc/PYKcPYKcKmPyr))
PyrTc:
 {0.183321}Pyrc = $pool
  (PyrTcPyrTcVmax*Pyrc)/(PyrTcPyrTcKmPyr*(1 + Pyrc/PyrTcPyrTcKmPyr))
TPIg:
 {4.07997}DHAPg = {4.07997}GA3Pg
  (TPIgTPIgVmax*DHAPg*(1 - GA3Pg/(TPIgTPIgKeq*DHAPg)))/(TPIgTPIgKmDHAP*(1 + DHAPg/TPIgTPIgKmDHAP + GA3Pg/TPIgTPIgKmGA3P))


#parameters
AKcAKck1 = 480.0
AKcAKck2 = 1000.0
AKgAKgk1 = 480.0
AKgAKgk2 = 1000.0
ALDgALDgKeq = 0.084
ALDgALDgKiADP = 1.51
ALDgALDgKiAMP = 3.65
ALDgALDgKiATP = 0.68
ALDgALDgKiGA3P = 0.098
ALDgALDgKmDHAP = 0.015
ALDgALDgKmFru16BP = 0.009
ALDgALDgKmGA3P = 0.067
ALDgALDgVmax = 560.0
ATPucATPuck = 50.0
ENOcENOcKeq = 4.17
ENOcENOcKm2PGA = 0.054
ENOcENOcKmPEP = 0.24
ENOcENOcVmax = 598.0
G3PDHgG3PDHgKeq = 17085.0
G3PDHgG3PDHgKmDHAP = 0.1
G3PDHgG3PDHgKmGly3P = 2.0
G3PDHgG3PDHgKmNAD = 0.4
G3PDHgG3PDHgKmNADH = 0.01
G3PDHgG3PDHgVmax = 465.0
GAPDHgGAPDHgKeq = 0.066
GAPDHgGAPDHgKm13BPGA = 0.1
GAPDHgGAPDHgKmGA3P = 0.15
GAPDHgGAPDHgKmNAD = 0.45
GAPDHgGAPDHgKmNADH = 0.02
GAPDHgGAPDHgVmax = 720.9
GDAgGDAgk = 600.0
GKgGKgKeq = 0.000837
GKgGKgKmADP = 0.56
GKgGKgKmATP = 0.24
GKgGKgKmGly = 0.44
GKgGKgKmGly3P = 3.83
GKgGKgVmax = 200.0
GPOcGPOcKmGly3P = 1.7
GPOcGPOcVmax = 368.0
GlcTcGlcTcKmGlc = 1.0
GlcTcGlcTcVmax = 111.7
GlcTcGlcTcalpha = 0.75
GlcTgGlcTgk1 = 250000.0
GlcTgGlcTgk2 = 250000.0
Glce = 5.0
Glye = 0.0
HXKgHXKgKeq = 759.0
HXKgHXKgKmADP = 0.126
HXKgHXKgKmATP = 0.116
HXKgHXKgKmGlc = 0.1
HXKgHXKgKmGlc6P = 12.0
HXKgHXKgVmax = 1774.68
P3GATgP3GATgk = 250.0
PFKgPFKgKeq = 1035.0
PFKgPFKgKi1 = 15.8
PFKgPFKgKi2 = 10.7
PFKgPFKgKmADP = 1.0
PFKgPFKgKmATP = 0.0648
PFKgPFKgKmFru6P = 0.999
PFKgPFKgKsATP = 0.0393
PFKgPFKgVmax = 1708.0
PGAMcPGAMcKeq = 0.17
PGAMcPGAMcKm2PGA = 0.16
PGAMcPGAMcKm3PGA = 0.15
PGAMcPGAMcVmax = 225.0
PGIg6PGg = 0.0841958
PGIgPGIgKeq = 0.457
PGIgPGIgKi6PG = 0.14
PGIgPGIgKmFru6P = 0.12
PGIgPGIgKmGlc6P = 0.4
PGIgPGIgVmax = 1305.0
PGKgPGKgKeq = 3377.0
PGKgPGKgKm13BPGA = 0.003
PGKgPGKgKm3PGA = 1.62
PGKgPGKgKmADP = 0.1
PGKgPGKgKmATP = 0.29
PGKgPGKgVmax = 2862.0
PYKcPYKcKeq = 10800.0
PYKcPYKcKiADP = 0.64
PYKcPYKcKiATP = 0.57
PYKcPYKcKmADP = 0.114
PYKcPYKcKmATP = 15.0
PYKcPYKcKmPEP = 0.34
PYKcPYKcKmPyr = 50.0
PYKcPYKcVmax = 1020.0
PYKcPYKcn = 2.5
Pig = 0.0
PyrTcPyrTcKmPyr = 1.96
PyrTcPyrTcVmax = 230.0
Pyre = 0.0
TPIgTPIgKeq = 0.046
TPIgTPIgKmDHAP = 1.2
TPIgTPIgKmGA3P = 0.25
TPIgTPIgVmax = 999.3
cytosol = 5.4549
glycosome = 0.2451
default = 1.0


#initial values
ADPc = 1.3165
ADPg = 1.519
AMPc = 2.2418
AMPg = 4.2405
ATPc = 0.3417
ATPg = 0.2405
B13PGAg = 0.5
DHAPc = 2.23132912
DHAPg = 8.483130623
Fru16BPg = 10.0
Fru6Pg = 0.5
GA3Pg = 2.5
Glc6Pg = 0.5
Glcc = 0.1
Glcg = 0.1
Gly3Pc = 2.76867088
Gly3Pg = 10.51686938
NADHg = 2.0
NADg = 2.0
P2GAc = 0.1
P3GAc = 0.1
P3GAg = 0.1
PEPc = 1.0
Pyrc = 10.0


#assignment rules