begin name
achcar2
end name

begin reactions
AKc {2.0}ADPc = {1.0}AMPc + {1.0}ATPc
AKg {44.5116}ADPg = {22.2558}AMPg + {22.2558}ATPg
ALDg {22.2558}Fru16BPg = {22.2558}DHAPg + {22.2558}GA3Pg
ATPuc {1.0}ATPc = {1.0}ADPc
DHAPtg {1.0}DHAPc = {22.2558}DHAPg
ENOc {1.0}PGA2c = {1.0}PEPc
G3PDHg {22.2558}DHAPg + {22.2558}NADHg = {22.2558}Gly3Pg + {22.2558}NADg
GAPDHg {1.0}$Pig + {22.2558}GA3Pg + {22.2558}NADg = {22.2558}BPGA13g + {22.2558}NADHg
GKg {22.2558}ADPg + {22.2558}Gly3Pg = {22.2558}ATPg + {22.2558}Glyg
GPOc {1.0}Gly3Pc = {1.0}DHAPc
GlcTc {1.0}$Glce = {1.0}Glcc
GlcTg {1.0}Glcc = {22.2558}Glcg
Gly3Ptg {22.2558}Gly3Pg = {1.0}Gly3Pc
GlyTc {1.0}Glyc = {1.0}$Glye
Glytg {22.2558}Glyg = {1.0}Glyc
HXKg {22.2558}ATPg + {22.2558}Glcg = {22.2558}ADPg + {22.2558}Glc6Pg
PFKg {22.2558}ATPg + {22.2558}Fru6Pg = {22.2558}ADPg + {22.2558}Fru16BPg
PGA3Tg {22.2558}PGA3g = {1.0}PGA3c
PGAMc {1.0}PGA3c = {1.0}PGA2c
PGIg {22.2558}Glc6Pg = {22.2558}Fru6Pg
PGKg {22.2558}BPGA13g + {22.2558}ADPg = {22.2558}PGA3g + {22.2558}ATPg
PYKc {1.0}ADPc + {1.0}PEPc = {1.0}ATPc + {1.0}Pyrc
PyrTc {1.0}Pyrc = {1.0}$Pyre
TPIg {22.2558}DHAPg = {22.2558}GA3Pg
end reactions

begin rate equations
AKc = (AKcAKck1*ADPc[t]^2 - AKcAKck2*AMPc[t]*ATPc[t])/cytosol
AKg = (AKgAKgk1*ADPg[t]^2 - AKgAKgk2*AMPg[t]*ATPg[t])/cytosol
ALDg = (ALDgALDgVmax*Fru16BPg[t]*(1 - (DHAPg[t]*GA3Pg[t])/(ALDgALDgKeq*Fru16BPg[t])))/(ALDgALDgKmFru16BP*cytosol*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)*(1 + DHAPg[t]/ALDgALDgKmDHAP + Fru16BPg[t]/(ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)) + GA3Pg[t]/ALDgALDgKmGA3P + (DHAPg[t]*GA3Pg[t])/(ALDgALDgKmDHAP*ALDgALDgKmGA3P) + (Fru16BPg[t]*GA3Pg[t])/(ALDgALDgKiGA3P*ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP))))
ATPuc = (ATPucATPuck*ATPc[t])/(cytosol*ADPc[t])
DHAPtg = (DHAPtgDHAPtgk*(DHAPc[t] - DHAPg[t]))/cytosol
ENOc = (ENOcENOcVmax*(1 - PEPc[t]/(ENOcENOcKeq*PGA2c[t]))*PGA2c[t])/(cytosol*ENOcENOcKmPGA2*(1 + PEPc[t]/ENOcENOcKmPEP + PGA2c[t]/ENOcENOcKmPGA2))
G3PDHg = (G3PDHgG3PDHgVmax*DHAPg[t]*(1 - (Gly3Pg[t]*NADg[t])/(G3PDHgG3PDHgKeq*DHAPg[t]*NADHg[t]))*NADHg[t])/(cytosol*G3PDHgG3PDHgKmDHAP*G3PDHgG3PDHgKmNADH*(1 + DHAPg[t]/G3PDHgG3PDHgKmDHAP + Gly3Pg[t]/G3PDHgG3PDHgKmGly3P)*(1 + NADg[t]/G3PDHgG3PDHgKmNAD + NADHg[t]/G3PDHgG3PDHgKmNADH))
GAPDHg = (GAPDHgGAPDHgVmax*GA3Pg[t]*NADg[t]*(1 - (BPGA13g[t]*NADHg[t])/(GAPDHgGAPDHgKeq*GA3Pg[t]*NADg[t])))/(cytosol*GAPDHgGAPDHgKmGA3P*GAPDHgGAPDHgKmNAD*(1 + BPGA13g[t]/GAPDHgGAPDHgKmBPGA13 + GA3Pg[t]/GAPDHgGAPDHgKmGA3P)*(1 + NADg[t]/GAPDHgGAPDHgKmNAD + NADHg[t]/GAPDHgGAPDHgKmNADH))
GKg = (GKgGKgVmax*ADPg[t]*Gly3Pg[t]*(1 - (ATPg[t]*Glyg[t])/(GKgGKgKeq*ADPg[t]*Gly3Pg[t])))/(cytosol*GKgGKgKmADP*GKgGKgKmGly3P*(1 + ADPg[t]/GKgGKgKmADP + ATPg[t]/GKgGKgKmATP)*(1 + Gly3Pg[t]/GKgGKgKmGly3P + Glyg[t]/GKgGKgKmGly))
GPOc = (GPOcGPOcVmax*Gly3Pc[t])/(cytosol*GPOcGPOcKmGly3P*(1 + SHAM/GPOcGPOcKiSHAM)*(1 + Gly3Pc[t]/GPOcGPOcKmGly3P))
GlcTc = (GlcTcGlcTcVmax*(Glce - Glcc[t]))/(cytosol*(Glce + GlcTcGlcTcKmGlc + Glcc[t] + (Glce*GlcTcGlcTcalpha*Glcc[t])/GlcTcGlcTcKmGlc))
GlcTg = (GlcTgGlcTgk*Glcc[t] - GlcTgGlcTgk*Glcg[t])/cytosol
Gly3Ptg = (Gly3PtgGly3Ptgk*(-Gly3Pc[t] + Gly3Pg[t]))/cytosol
GlyTc = (GlyTcGlyTck*(-Glye + Glyc[t]) + (GlyTcGlyTcVrmax*GlyTcGlyTcy*(-Glye + Glyc[t]))/(Glye*GlyTcGlyTcy + GlyTcGlyTcKmGlye*GlyTcGlyTcy + Glyc[t] + (Glye*GlyTcGlyTcx*(1 - GlyTcGlyTcx + GlyTcGlyTcy)*Glyc[t])/(GlyTcGlyTcKmGlye*GlyTcGlyTcy)))/cytosol
Glytg = (GlytgGlytgk*(-Glyc[t] + Glyg[t]))/cytosol
HXKg = (HXKgHXKgVmax*ATPg[t]*(1 - (ADPg[t]*Glc6Pg[t])/(HXKgHXKgKeq*ATPg[t]*Glcg[t]))*Glcg[t])/(cytosol*HXKgHXKgKmATP*HXKgHXKgKmGlc*(1 + ADPg[t]/HXKgHXKgKmADP + ATPg[t]/HXKgHXKgKmATP)*(1 + Glc6Pg[t]/HXKgHXKgKmGlc6P + Glcg[t]/HXKgHXKgKmGlc))
PFKg = (PFKgPFKgKi1*PFKgPFKgVmax*ATPg[t]*(1 - (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKeq*ATPg[t]*Fru6Pg[t]))*Fru6Pg[t])/(cytosol*PFKgPFKgKmATP*PFKgPFKgKmFru6P*(PFKgPFKgKi1 + Fru16BPg[t])*(PFKgPFKgKsATP/PFKgPFKgKmATP + ADPg[t]/PFKgPFKgKmADP + ATPg[t]/PFKgPFKgKmATP + (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKi2*PFKgPFKgKmADP) + Fru6Pg[t]/PFKgPFKgKmFru6P + (ATPg[t]*Fru6Pg[t])/(PFKgPFKgKmATP*PFKgPFKgKmFru6P)))
PGA3Tg = (-(PGA3TgPGA3Tgk*PGA3c[t]) + PGA3TgPGA3Tgk*PGA3g[t])/cytosol
PGAMc = (PGAMcPGAMcVmax*(1 - PGA2c[t]/(PGAMcPGAMcKeq*PGA3c[t]))*PGA3c[t])/(cytosol*PGAMcPGAMcKmPGA3*(1 + PGA2c[t]/PGAMcPGAMcKmPGA2 + PGA3c[t]/PGAMcPGAMcKmPGA3))
PGIg = (PGIgPGIgVmax*(1 - Fru6Pg[t]/(PGIgPGIgKeq*Glc6Pg[t]))*Glc6Pg[t])/(cytosol*PGIgPGIgKmGlc6P*(1 + Fru6Pg[t]/PGIgPGIgKmFru6P + Glc6Pg[t]/PGIgPGIgKmGlc6P))
PGKg = (PGKgPGKgVmax*ADPg[t]*BPGA13g[t]*(1 - (ATPg[t]*PGA3g[t])/(PGKgPGKgKeq*ADPg[t]*BPGA13g[t])))/(cytosol*PGKgPGKgKmADP*PGKgPGKgKmBPGA13*(1 + ADPg[t]/PGKgPGKgKmADP + ATPg[t]/PGKgPGKgKmATP)*(1 + BPGA13g[t]/PGKgPGKgKmBPGA13 + PGA3g[t]/PGKgPGKgKmPGA3))
PYKc = (PYKcPYKcVmax*ADPc[t]*(PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn*(1 - (ATPc[t]*Pyrc[t])/(PYKcPYKcKeq*ADPc[t]*PEPc[t])))/(cytosol*PYKcPYKcKmADP*(1 + ADPc[t]/PYKcPYKcKmADP + ATPc[t]/PYKcPYKcKmATP)*(1 + (PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn + Pyrc[t]/PYKcPYKcKmPyr))
PyrTc = (PyrTcPyrTcVmax*Pyrc[t])/(cytosol*PyrTcPyrTcKmPyr*(1 + Pyrc[t]/PyrTcPyrTcKmPyr))
TPIg = (TPIgTPIgVmax*DHAPg[t]*(1 - GA3Pg[t]/(TPIgTPIgKeq*DHAPg[t])))/(cytosol*TPIgTPIgKmDHAP*(1 + DHAPg[t]/TPIgTPIgKmDHAP + GA3Pg[t]/TPIgTPIgKmGA3P))
end rate equations

begin parameters
AKcAKck1 = 442.0
AKcAKck2 = 1000.0
AKgAKgk1 = 442.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
DHAPtgDHAPtgk = 1000.0
ENOcENOcKeq = 4.12
ENOcENOcKmPEP = 0.24
ENOcENOcKmPGA2 = 0.054
ENOcENOcVmax = 598.0
EXTERNAL = 0.0
G3PDHgG3PDHgKeq = 17085.0
G3PDHgG3PDHgKmDHAP = 0.1
G3PDHgG3PDHgKmGly3P = 2.0
G3PDHgG3PDHgKmNAD = 0.4
G3PDHgG3PDHgKmNADH = 0.01
G3PDHgG3PDHgVmax = 465.0
GAPDHgGAPDHgKeq = 0.066
GAPDHgGAPDHgKmBPGA13 = 0.1
GAPDHgGAPDHgKmGA3P = 0.15
GAPDHgGAPDHgKmNAD = 0.45
GAPDHgGAPDHgKmNADH = 0.02
GAPDHgGAPDHgVmax = 720.9
GKgGKgKeq = 0.000837
GKgGKgKmADP = 0.56
GKgGKgKmATP = 0.24
GKgGKgKmGly = 0.44
GKgGKgKmGly3P = 3.83
GKgGKgVmax = 200.0
GPOcGPOcKiSHAM = 0.015
GPOcGPOcKmGly3P = 1.7
GPOcGPOcVmax = 368.0
GlcTcGlcTcKmGlc = 1.1
GlcTcGlcTcVmax = 111.7
GlcTcGlcTcalpha = 0.75
GlcTgGlcTgk = 250000.0
Gly3PtgGly3Ptgk = 1000.0
GlyTcGlyTcKmGlye = 0.17
GlyTcGlyTcVrmax = 85.0
GlyTcGlyTck = 8.5
GlyTcGlyTcx = 1.0
GlyTcGlyTcy = 1.0
GlytgGlytgk = 1000.0
HXKgHXKgKeq = 759.0
HXKgHXKgKmADP = 0.126
HXKgHXKgKmATP = 0.116
HXKgHXKgKmGlc = 0.1
HXKgHXKgKmGlc6P = 12.0
HXKgHXKgVmax = 1929.0
O2c = 1.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
PGA3TgPGA3Tgk = 250.0
PGAMcPGAMcKeq = 0.17
PGAMcPGAMcKmPGA2 = 0.16
PGAMcPGAMcKmPGA3 = 0.15
PGAMcPGAMcVmax = 225.0
PGIgPGIgKeq = 0.46
PGIgPGIgKmFru6P = 0.12
PGIgPGIgKmGlc6P = 0.4
PGIgPGIgVmax = 1305.0
PGKgPGKgKeq = 3376.0
PGKgPGKgKmADP = 0.1
PGKgPGKgKmATP = 0.29
PGKgPGKgKmBPGA13 = 0.003
PGKgPGKgKmPGA3 = 1.62
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
PyrTcPyrTcKmPyr = 1.96
PyrTcPyrTcVmax = 230.0
SHAM = 0.0
TPIgTPIgKeq = 0.046
TPIgTPIgKmDHAP = 1.2
TPIgTPIgKmGA3P = 0.25
TPIgTPIgVmax = 999.3
cytosol = 5.4549
default = 1.0
glycosome = 0.2451
Glce = 5.0
Glye = "1e-06"
Pig = 0.0
Pyre = 0.0
defaultcompartment = 1.0
end parameters

begin initial conditions
ADPc[0] = ADPci
ADPg[0] = ADPgi
AMPc[0] = AMPci
AMPg[0] = AMPgi
ATPc[0] = ATPci
ATPg[0] = ATPgi
BPGA13g[0] = BPGA13gi
DHAPc[0] = DHAPci
DHAPg[0] = DHAPgi
Fru16BPg[0] = Fru16BPgi
Fru6Pg[0] = Fru6Pgi
GA3Pg[0] = GA3Pgi
Glc6Pg[0] = Glc6Pgi
Glcc[0] = Glcci
Glcg[0] = Glcgi
Gly3Pc[0] = Gly3Pci
Gly3Pg[0] = Gly3Pgi
Glyc[0] = Glyci
Glyg[0] = Glygi
NADHg[0] = NADHgi
NADg[0] = NADgi
PEPc[0] = PEPci
PGA2c[0] = PGA2ci
PGA3c[0] = PGA3ci
PGA3g[0] = PGA3gi
Pyrc[0] = Pyrci
end initial conditions

begin initial values
ADPci = 1.0
ADPgi = 1.0
AMPci = 0.1
AMPgi = 0.1
ATPci = 4.9
ATPgi = 4.9
BPGA13gi = 0.4
DHAPci = 4.6
DHAPgi = 4.6
Fru16BPgi = 13.0
Fru6Pgi = 0.2
GA3Pgi = 0.14
Glc6Pgi = 0.8
Glcci = 0.01
Glcgi = 0.01
Gly3Pci = 1.5
Gly3Pgi = 1.5
Glyci = 1e-06
Glygi = 1e-06
NADHgi = 0.1
NADgi = 3.9
PEPci = 1.4
PGA2ci = 0.5
PGA3ci = 17.0
PGA3gi = 17.0
Pyrci = 14.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units