begin name
messiha2
end name

begin reactions
vE4Psink {1.0}E4P = {1.0}$EXTERNAL
vGND {1.0}P6G = {1.0}NADPH + {1.0}Ru5P
vNADPHoxidase {1.0}NADPH = {1.0}$EXTERNAL
vR5Psink {1.0}R5P = {1.0}$EXTERNAL
vRKI {1.0}Ru5P = {1.0}R5P
vRPE {1.0}Ru5P = {1.0}X5P
vSOL {1.0}G6L = {1.0}P6G
vTAL {1.0}S7P = {1.0}E4P
vTKLE4P {1.0}E4P + {1.0}X5P = {1.0}$EXTERNAL
vTKLR5P {1.0}R5P + {1.0}X5P = {1.0}S7P
vZWF {1.0}$EXTERNAL = {1.0}G6L + {1.0}NADPH
end reactions

begin rate equations
vE4Psink = E4Psinkk*E4P[t]
vGND = (GND1*kcatGND1*NADP*P6G[t])/(KnadpGND1*Kp6gGND1*(1 + NADP/KnadpGND1 + NADPH[t]/KnadphGND1)*(1 + P6G[t]/Kp6gGND1 + Ru5P[t]/Kru5pGND1)) + (GND2*kcatGND2*NADP*P6G[t])/((1 + NADP/KnadpGND2 + NADPH[t]/KnadphGND2)*(1 + P6G[t]/Kp6gGND2 + Ru5P[t]/Kru5pGND2))
vNADPHoxidase = NADPHoxidasek*NADPH[t]
vR5Psink = R5Psinkk*R5P[t]
vRKI = (RKI1*RKIkcat*(-(R5P[t]/RKIKeq) + Ru5P[t]))/(RKIKru5p*(1 + R5P[t]/RKIKr5p + Ru5P[t]/RKIKru5p))
vRPE = (RPE1*RPEkcat*(Ru5P[t] - X5P[t]/RPEKeq))/(RPEKru5p*(1 + Ru5P[t]/RPEKru5p + X5P[t]/RPEKx5p))
vSOL = (SOL3*SOLkcat*G6L[t])/(SOLKg6l*(1 + G6L[t]/SOLKg6l + P6G[t]/SOLKp6g))
vTAL = (kcatTAL1*TAL1*(-((F6P*E4P[t])/TALKeq) + GAP*S7P[t]))/((1 + F6P/Kf6pTAL1 + GAP/KgapTAL1)*KgapTAL1*Ks7pTAL1*(1 + E4P[t]/Ke4pTAL1 + S7P[t]/Ks7pTAL1)) + (NQM1*TALkcatNQM1*(-((F6P*E4P[t])/TALKeq) + GAP*S7P[t]))/((1 + F6P/TALKf6pNQM1 + GAP/TALKgapNQM1)*TALKgapNQM1*TALKs7pNQM1*(1 + E4P[t]/TALKe4pNQM1 + S7P[t]/TALKs7pNQM1))
vTKLE4P = (TKL1*TKLE4Pkcat*(-((F6P*GAP)/TKLE4PKeq) + E4P[t]*X5P[t]))/(Ke4pTAL*Kx5pTAL*(1 + F6P/Kf6pTAL + E4P[t]/Ke4pTAL + R5P[t]/Kr5pTAL + S7P[t]/Ks7pTAL)*(1 + GAP/KgapTAL + X5P[t]/Kx5pTAL))
vTKLR5P = (TKL1*TKLR5Pkcat*(-((GAP*S7P[t])/TKLR5PKeq) + R5P[t]*X5P[t]))/(Kr5pTAL*Kx5pTAL*(1 + F6P/Kf6pTAL + E4P[t]/Ke4pTAL + R5P[t]/Kr5pTAL + S7P[t]/Ks7pTAL)*(1 + GAP/KgapTAL + X5P[t]/Kx5pTAL))
vZWF = (G6P*NADP*ZWF1*ZWFkcat)/(ZWFKg6p*ZWFKnadp*(1 + G6P/ZWFKg6p + G6L[t]/ZWFKg6l)*(1 + NADP/ZWFKnadp + NADPH[t]/ZWFKnadph))
end rate equations

begin parameters
E4Psinkk = 1.0
EXTERNAL = 0.0
F6P = 0.325
GAP = 0.067
GND1 = 0.013
GND2 = 0.003
Ke4pTAL = 0.946
Ke4pTAL1 = 0.362
Kf6pTAL = 1.1
Kf6pTAL1 = 1.44
KgapTAL = 0.1
KgapTAL1 = 0.272
KnadpGND1 = 0.094
KnadpGND2 = 0.094
KnadphGND1 = 0.055
KnadphGND2 = 0.055
Kp6gGND1 = 0.062
Kp6gGND2 = 0.115
Kr5pTAL = 0.235
Kru5pGND1 = 0.1
Kru5pGND2 = 0.1
Ks7pTAL = 0.15
Ks7pTAL1 = 0.786
Kx5pTAL = 0.67
NADPHoxidasek = 1.0
NQM1 = 0.02
R5Psinkk = 1.0
RKI1 = 0.05
RKIKeq = 4.0
RKIKiru5p = 9.88
RKIKr5p = 5.7
RKIKru5p = 2.47
RKIkcat = 335.0
RPE1 = 0.03
RPEKeq = 1.4
RPEKru5p = 5.97
RPEKx5p = 7.7
RPEkcat = 4020.0
SOL3 = 0.0296
SOLKg6l = 0.8
SOLKp6g = 0.5
SOLkcat = 4.3
TAL1 = 0.144
TALKe4pNQM1 = 0.305
TALKeq = 1.05
TALKf6pNQM1 = 1.04
TALKgapNQM1 = 0.272
TALKs7pNQM1 = 0.786
TALkcatNQM1 = 0.694
TKL1 = 0.455
TKLE4PKeq = 10.0
TKLE4Pkcat = 47.1
TKLR5PKeq = 1.2
TKLR5Pkcat = 40.5
ZWF1 = 0.02
ZWFKg6l = 0.01
ZWFKg6p = 0.042
ZWFKnadp = 0.045
ZWFKnadph = 0.017
ZWFkcat = 189.0
cell = 1.0
kcatGND1 = 28.0
kcatGND2 = 27.3
kcatTAL1 = 0.694
sumNADP = 0.33
default = 1.0
end parameters

begin initial conditions
E4P[0] = E4Pi
G6L[0] = G6Li
NADPH[0] = NADPHi
P6G[0] = P6Gi
R5P[0] = R5Pi
Ru5P[0] = Ru5Pi
S7P[0] = S7Pi
X5P[0] = X5Pi
end initial conditions

begin initial values
E4Pi = 0.029
G6Li = 0.1
NADPHi = 0.16
P6Gi = 0.25
R5Pi = 0.118
Ru5Pi = 0.033
S7Pi = 0.082
X5Pi = 0.041
end initial values

begin assignment rules
G6P := 0.9
NADP := sumNADP - NADPH[t]
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