begin name
lambeth
end name

begin reactions
v1 {1.0}$Ph + {1.0}$GLY = {1.0}G1P
v10 {1.0}P2G = {1.0}PEP
v11 {1.0}PEP + {1.0}$adp = {1.0}PYR + {1.0}$atp
v12 {1.0}PYR + {1.0}NADH = {1.0}$LAC + {1.0}NAD
v2 {1.0}G1P = {1.0}G6P
v3 {1.0}G6P = {1.0}F6P
v4 {1.0}F6P + {1.0}$atp = {1.0}FDP + {1.0}$adp
v5 {1.0}FDP = {1.0}DHAP + {1.0}GAP
v6 {1.0}GAP = {1.0}DHAP
v7 {1.0}GAP + {1.0}NAD + {1.0}$Ph = {1.0}DPG + {1.0}NADH
v8 {1.0}DPG + {1.0}$adp = {1.0}P3G + {1.0}$atp
v9 {1.0}P3G = {1.0}P2G
end reactions

begin rate equations
v1 = fracA*((GLY*Ph*Vfgly)/(KgpAigly*KgpApi*(1 + GLY/KgpAglyb + GLY/KgpAglyf + (GLY*Ph)/(KgpAglyf*KgpAipi) + Ph/KgpApi + G1P[t]/KgpAg1p + (GLY*G1P[t])/(KgpAglyb*KgpAig1p))) - (3.2258064516129035*GLY*Vfgly*G1P[t])/(KgpAigly*KgpApi*(1 + GLY/KgpAglyb + GLY/KgpAglyf + (GLY*Ph)/(KgpAglyf*KgpAipi) + Ph/KgpApi + G1P[t]/KgpAg1p + (GLY*G1P[t])/(KgpAglyb*KgpAig1p)))) + fracB*((50.*amp^nH*GLY*Ph*Vfgly)/(Kgpamp^nH*(1 + (50.*amp^nH)/Kgpamp^nH)*KgpBiglyf*KgpBpi*(1 + GLY/KgpBiglyb + GLY/KgpBipi + Ph/KgpBiglyf + (GLY*Ph)/(KgpBiglyf*KgpBpi) + G1P[t]/KgpBig1p + (GLY*G1P[t])/(KgpBg1p*KgpBiglyb))) - (161.29032258064518*amp^nH*GLY*Vfgly*G1P[t])/(Kgpamp^nH*(1 + (50.*amp^nH)/Kgpamp^nH)*KgpBiglyf*KgpBpi*(1 + GLY/KgpBiglyb + GLY/KgpBipi + Ph/KgpBiglyf + (GLY*Ph)/(KgpBiglyf*KgpBpi) + G1P[t]/KgpBig1p + (GLY*G1P[t])/(KgpBg1p*KgpBiglyb))))
v10 = ((Vfen*P2G[t])/Ken2pg - (2.0408163265306123*Vfen*PEP[t])/Ken2pg)/(1 + P2G[t]/Ken2pg + PEP[t]/Kenpep)
v11 = ((adp*Vfpk*PEP[t])/(Kpkadp*Kpkpep) - (atp*Vfpk*PYR[t])/(10304*Kpkadp*Kpkpep))/(1 + adp/Kpkadp + atp/Kpkatp + PEP[t]/Kpkpep + (adp*PEP[t])/(Kpkadp*Kpkpep) + PYR[t]/Kpkpyr + (atp*PYR[t])/(Kpkatp*Kpkpyr))
v12 = (-(LAC*Vfldh*NAD[t])/(16198*Kldhnadh*Kldhpyr) + (Vfldh*NADH[t]*PYR[t])/(Kldhnadh*Kldhpyr))/(1 + LAC/Kldhlac + NAD[t]/Kldhnad + (LAC*NAD[t])/(Kldhlac*Kldhnad) + NADH[t]/Kldhnadh + PYR[t]/Kldhpyr + (NADH[t]*PYR[t])/(Kldhnadh*Kldhpyr))
v2 = ((Vfpglm*G1P[t])/Kpglmg1p - (0.06016847172081829*Vfpglm*G6P[t])/Kpglmg1p)/(1 + G1P[t]/Kpglmg1p + G6P[t]/Kpglmg6p)
v3 = (-((Vbpgi*F6P[t])/Kpgif6p) + (0.45*Vbpgi*G6P[t])/Kpgif6p)/(1 + F6P[t]/Kpgif6p + G6P[t]/Kpgig6p)
v4 = (atp*Vfpfk*F6P[t]*(1 + ((1 + (amp*en)/Kpfkamp)^4*Kpfkatp*Kpfkf6p*(1 + atp/Kpfkiatp)^4*Lo*((1 + atp/KpfkatpT)*(1 + F6P[t]/Kpfkf6pT) + FDP[t]/KpfkfdpT + (adp*(1 + FDP[t]/KpfkfdpT))/KpfkadpT)^3)/((1 + amp/Kpfkamp)^4*KpfkatpT*Kpfkf6pT*(1 + (atp*dn)/Kpfkiatp)^4*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)^3)))/(Kpfkatp*Kpfkf6p*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)*(1 + ((1 + (amp*en)/Kpfkamp)^4*(1 + atp/Kpfkiatp)^4*Lo*((1 + atp/KpfkatpT)*(1 + F6P[t]/Kpfkf6pT) + FDP[t]/KpfkfdpT + (adp*(1 + FDP[t]/KpfkfdpT))/KpfkadpT)^4)/((1 + amp/Kpfkamp)^4*(1 + (atp*dn)/Kpfkiatp)^4*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)^4))) - (0.00411742907728414*adp*Vfpfk*FDP[t]*(1 + ((1 + (amp*en)/Kpfkamp)^4*Kpfkatp*Kpfkf6p*(1 + atp/Kpfkiatp)^4*Lo*((1 + atp/KpfkatpT)*(1 + F6P[t]/Kpfkf6pT) + FDP[t]/KpfkfdpT + (adp*(1 + FDP[t]/KpfkfdpT))/KpfkadpT)^3)/((1 + amp/Kpfkamp)^4*KpfkatpT*Kpfkf6pT*(1 + (atp*dn)/Kpfkiatp)^4*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)^3)))/(Kpfkatp*Kpfkf6p*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)*(1 + ((1 + (amp*en)/Kpfkamp)^4*(1 + atp/Kpfkiatp)^4*Lo*((1 + atp/KpfkatpT)*(1 + F6P[t]/Kpfkf6pT) + FDP[t]/KpfkfdpT + (adp*(1 + FDP[t]/KpfkfdpT))/KpfkadpT)^4)/((1 + amp/Kpfkamp)^4*(1 + (atp*dn)/Kpfkiatp)^4*((1 + atp/Kpfkatp)*(1 + F6P[t]/Kpfkf6p) + FDP[t]/Kpfkfdp + (adp*(1 + FDP[t]/Kpfkfdp))/Kpfkadp)^4)))
v5 = ((Vfald*FDP[t])/Kaldfdp - (10526.315789473683*Vfald*DHAP[t]*GAP[t])/Kaldfdp)/(1 + DHAP[t]/Kalddhap + FDP[t]/Kaldfdp + GAP[t]/Kaldgap)
v6 = ((-0.052083333333333336*Vftpi*DHAP[t])/Ktpigap + (Vftpi*GAP[t])/Ktpigap)/(1 + DHAP[t]/Ktpidhap + GAP[t]/Ktpigap)
v7 = ((Ph*Vfgad*GAP[t]*NAD[t])/(Kgapdhgap*Kgapdhnad*Kgapdhpi) - (11.2359550561798*Vfgad*DPG[t]*NADH[t])/(Kgapdhgap*Kgapdhnad*Kgapdhpi))/(1 + Ph/Kgapdhpi + DPG[t]/Kgapdh13dpg + GAP[t]/Kgapdhgap + NAD[t]/Kgapdhnad + (GAP[t]*NAD[t])/(Kgapdhgap*Kgapdhnad) + (Ph*GAP[t]*NAD[t])/(Kgapdhgap*Kgapdhnad*Kgapdhpi) + NADH[t]/Kgapdhnadh + (DPG[t]*NADH[t])/(Kgapdh13dpg*Kgapdhnadh))
v8 = ((57109*adp*Vbpgk*DPG[t])/(Kpgk3pg*Kpgkatp) - (atp*Vbpgk*P3G[t])/(Kpgk3pg*Kpgkatp))/(1 + adp/Kpgkadp + atp/Kpgkatp + DPG[t]/Kpgk13dpg + (adp*DPG[t])/(Kpgk13dpg*Kpgkadp) + P3G[t]/Kpgk3pg + (atp*P3G[t])/(Kpgk3pg*Kpgkatp))
v9 = ((-5.555555555555555*Vfpgm*P2G[t])/Kpgm3pg + (Vfpgm*P3G[t])/Kpgm3pg)/(1 + P2G[t]/Kpgm2pg + P3G[t]/Kpgm3pg)
end rate equations

begin parameters
Kadkadp = 0.00035
Kadkamp = 0.00032
Kadkatp = 0.00027
Kalddhap = 0.002
Kaldfdp = 0.0005
Kaldgap = 0.001
KckPCr = 0.00111
Kckcr = 0.0038
KckiPCr = 0.0039
Kckiadp = 0.000135
Kckiatp = 0.0035
Ken2pg = 0.0001
Kenpep = 0.00037
KeqCK = 233.0
Kgapdh13dpg = "8e-07"
Kgapdhgap = "2.5e-06"
Kgapdhnad = "9e-05"
Kgapdhnadh = "3.3e-06"
Kgapdhpi = 0.00029
KgpAg1p = 0.0027
KgpAglyb = 0.00015
KgpAglyf = 0.0017
KgpAig1p = 0.0101
KgpAigly = 0.002
KgpAipi = 0.0047
KgpApi = 0.004
KgpBg1p = 0.0015
KgpBig1p = 0.0074
KgpBiglyb = 0.0044
KgpBiglyf = 0.015
KgpBipi = 0.0046
KgpBpi = 0.0002
Kgpamp = "9.7e-05"
Kldhlac = 0.01717
Kldhnad = 0.000849
Kldhnadh = "2.167e-06"
Kldhpyr = 0.000335
Kpfkadp = 0.00271
KpfkadpT = 0.00271
Kpfkamp = "6e-05"
Kpfkatp = "8e-05"
KpfkatpT = 0.00025
Kpfkf6p = 0.00018
Kpfkf6pT = 0.02
Kpfkfdp = 0.0042
KpfkfdpT = 0.0042
Kpfkiatp = 0.00087
Kpgif6p = 0.000119
Kpgig6p = 0.00048
Kpgk13dpg = "1.9e-06"
Kpgk3pg = 0.0012
Kpgkadp = "8.3e-05"
Kpgkatp = 0.00035
Kpglmg1p = "6.3e-05"
Kpglmg6p = "3e-05"
Kpgm2pg = "1.4e-05"
Kpgm3pg = 0.0002
Kpkadp = 0.0003
Kpkatp = 0.00113
Kpkpep = "8e-05"
Kpkpyr = 0.00705
Ktpidhap = 0.00061
Ktpigap = 0.00032
Lo = 13.0
Vbpgi = 0.88
Vbpgk = 1.12
Vfadk = 0.88
Vfald = 0.104
Vfen = 0.5
Vfgad = 1.65
Vfgly = 0.15
Vfldh = 1.92
Vfpfk = 0.056
Vfpglm = 0.48
Vfpgm = 1.12
Vfpk = 3.0
Vftpi = 12.0
VrevCK = 0.5
amp = "9.52098310973979e-06"
dn = 0.01
en = 0.01
fracA = 0.4
fracB = 0.6
k = 0.075
kout = 0.2
nH = 1.75
GLY = 0.112
LAC = 0.00205129800531363
Ph = 0.0307416000593277
adp = 0.000405079560143532
atp = 0.00779841945674673
defaultcompartment = 1.0
end parameters

begin initial conditions
DHAP[0] = DHAPi
DPG[0] = DPGi
F6P[0] = F6Pi
FDP[0] = FDPi
G1P[0] = G1Pi
G6P[0] = G6Pi
GAP[0] = GAPi
NAD[0] = NADi
NADH[0] = NADHi
P2G[0] = P2Gi
P3G[0] = P3Gi
PEP[0] = PEPi
PYR[0] = PYRi
end initial conditions

begin initial values
DHAPi = 7e-05
DPGi = 6.5e-05
F6Pi = 0.0002
FDPi = 7e-05
G1Pi = 5e-05
G6Pi = 0.00075
GAPi = 3e-05
NADi = 0.0005
NADHi = 5e-07
P2Gi = 5e-06
P3Gi = 5e-05
PEPi = 1e-05
PYRi = 0.000144
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