begin name
jamshidi
end name

begin reactions
v1 {1.0}ADO = {1.0}INO
v10 {1.0}ATP = {1.0}ADP
v11 {1.0}DPG23 = {1.0}PG3
v12 {1.0}DPG13 = {1.0}DPG23
v13 {1.0}PG2 = {1.0}PEP
v14 {1.0}G6P = {1.0}GL6P + {1.0}NADPH
v15 {1.0}GAP = {1.0}DPG13 + {1.0}NADH
v16 {1.0}GO6P = {1.0}NADPH + {1.0}RU5P
v17 {1.0}GSH = {1.0}$GSHex
v18 {1.0}NADPH = {2.0}GSH
v19 {1.0}HX + {1.0}PRPP = {1.0}IMP
v2 {1.0}ADE = {1.0}$adeex
v20 {1.0}ATP = {1.0}ADP + {1.0}G6P
v21 {1.0}HX = {1.0}$HXex
v22 {1.0}IMP = {1.0}INO
v23 {1.0}INO = {1.0}$inoex
v24 {1.0}$KLEak = {1.0}KI
v25 {1.0}NADH + {1.0}PYR = {1.0}LAC
v26 {1.0}LAC = {1.0}$LACEXT
v27 {1.0}$NALEak = {1.0}NAI
v28 {1.0}ADP + {1.0}MG = {1.0}MGADP
v29 {1.0}AMP + {1.0}MG = {1.0}MGAMP
v3 {2.0}ADP = {1.0}AMP + {1.0}ATP
v30 {1.0}ATP + {1.0}MG = {1.0}MGATP
v31 {1.0}DPG23 + {1.0}MG = {1.0}MGDPG23
v32 {1.0}PYR = {1.0}$PYREXT
v33 {1.0}ATP + {1.0}F6P = {1.0}ADP + {1.0}FDP
v34 {1.0}G6P = {1.0}F6P
v35 {1.0}ADP + {1.0}DPG13 = {1.0}ATP + {1.0}PG3
v36 {1.0}GL6P = {1.0}GO6P
v37 {1.0}PG3 = {1.0}PG2
v38 {1.0}ADP + {1.0}PEP = {1.0}PYR + {1.0}ATP
v39 {1.0}INO = {1.0}HX + {1.0}R1P
v4 {1.0}ADO = {1.0}$adoex
v40 {1.0}R1P = {1.0}R5P
v41 {1.0}R5P + {1.0}ATP = {1.0}AMP + {1.0}PRPP
v42 {3.0}NAI + {1.0}ATP = {1.0}ADP + {2.0}KI
v43 {1.0}RU5P = {1.0}R5P
v44 {1.0}GAP + {1.0}S7P = {1.0}E4P + {1.0}F6P
v45 {1.0}X5P + {1.0}R5P = {1.0}GAP + {1.0}S7P
v46 {1.0}X5P + {1.0}E4P = {1.0}F6P + {1.0}GAP
v47 {1.0}DHAP = {1.0}GAP
v48 {1.0}RU5P = {1.0}X5P
v5 {1.0}PRPP + {1.0}ADE = {1.0}AMP
v6 {1.0}ADO + {1.0}ATP = {1.0}ADP + {1.0}AMP
v7 {1.0}FDP = {1.0}DHAP + {1.0}GAP
v8 {1.0}AMP = {1.0}ADO
v9 {1.0}AMP = {1.0}IMP
end reactions

begin rate equations
v1 = (adaVM*scaleada*ADO[t])/(adaKA + ADO[t])
v10 = katpase*scaleatpase*ATP[t]
v11 = (dpgaseVM*scaledpgase*DPG23[t])/(dpgasek + DPG23[t])
v12 = (KDPGM*scaledpgm*DPG13[t])/(1 + DPG23[t]/dpgmk)
v13 = E^(19.876 - 6152.2/TEMP)*scaleen*(-(KEN2*PEP[t]) + KEN1*PG2[t])
v14 = (3600000*NUMG6PDH*scaleg6pd)/DENG6PDH
v15 = E^(16.203 - 5021.3/TEMP)*scalegapdh*(KGA1*PI*GAP[t]*(NADHT - NADH[t]) - HI*KGA2*DPG13[t]*NADH[t])
v16 = ((160.84 - 39.173*HI + 4.1322*HI^2)*NUMG16PD*scalegl6pd)/DENG16PD
v17 = gshconst*(-2880 + 902.33*PHI - 65.272*PHI^2)*scalegshr*GSH[t]
v18 = (0.00429046444277593*gssgNUM*(-2880 + 902.33*PHI - 65.272*PHI^2)*scalegssgr)/gssgDEN
v19 = (E^(14.543 - 4501.5/TEMP)*hgprtVM*scalehgprt*HX[t]*PRPP[t])/((hgprtKHX + HX[t])*(hgprtKPRPP + PRPP[t]))
v2 = scaleadeex*vmade*(-(adeex/(adeex + kmade)) + ADE[t]/(kmade + ADE[t]))
v20 = (GLC*scalehk*VM1HK*(1 + (2.1187*MG[t])/KMGAMG)*MGATP[t])/(KMAGLC*KMGATP*N1HK)
v21 = E^(14.543 - 4501.5/TEMP)*scalehxex*(hxexPM*HX[t] + (hxexVM*HX[t])/(hxexKM + HX[t]))
v22 = (E^(14.543 - 4501.5/TEMP)*KIASE*scaleimpase*IMP[t])/(1 + INO[t]/KINO)
v23 = vmino*(-(inoex/(inoex + kmino)) + INO[t]/(kmino + INO[t]))
v24 = E^(22.7283 - 7045.7976/TEMP)*(kV1 + kV2)*(-1.16 + 0.3*PHI)*scalekleak
v25 = E^(20.237 - 6267.2/TEMP)*scaleldh*(-(KLD2*LAC[t]*(NADHT - NADH[t])) + KLD1*NADH[t]*PYR[t])
v26 = E^(14.543 - 4501.5/TEMP)*scalelex*VLAC
v27 = E^(29.2222 - 9058.88/TEMP)*(naV1 + naV2)*(-1.16 + 0.3*PHI)*scalenaleak
v28 = kAss*((ADP[t]*MG[t])/(mgadpk + MG[t]) - MGADP[t])
v29 = kAss*((AMP[t]*MG[t])/(mgampk + MG[t]) - MGAMP[t])
v3 = E^(19.876 - 6152.2/TEMP)*scaleadk*(adkKAP1*ADP[t]^2 - adkKAP2*AMP[t]*ATP[t])
v30 = kAss*((ATP[t]*MG[t])/(mgatpk + MG[t]) - MGATP[t])
v31 = kAss*((DPG23[t]*MG[t])/(mgdpgk + MG[t]) - MGDPG23[t])
v32 = -(E^(14.52096 - 4501.5/TEMP)*KPYR*scalepex*(PYREXT - PYR[t]))
v33 = (pfkVM*scalepfk*F6P[t]*MGATP[t])/((1 + LPFK)*(pfkKF6P + F6P[t])*(pfkKMGATP + MGATP[t]))
v34 = E^(13.042 - 4040.7/TEMP)*scalepgi*(-(kpgi2*F6P[t]) + kpgi1*G6P[t])
v35 = E^(12.56 - 3895.3/TEMP)*scalepgk*(kpgk1*ADP[t]*DPG13[t] - kpgk2*ATP[t]*PG3[t])
v36 = (E^(14.167 - 4380.4/TEMP)*pglaseVM*scalepglase*GL6P[t])/(pglaseKM + GL6P[t])
v37 = E^(16.907 - 5236.7/TEMP)*scalepgm*(-(KPGM2*PG2[t]) + KPGM1*PG3[t])
v38 = (pkDEN1*scalepk*pkVM[PHI, TEMP])/(1 + (pkDEN2*pkL0)/pkDEN3)
v39 = E^(14.543 - 4501.5/TEMP)*scalepnpase*(PI*pnpaseKA1*INO[t] - pnpaseKA2*HX[t]*R1P[t])
v4 = scaleadoex*vmaex*(-(adoex/(adoex + kmaex)) + ADO[t]/(kmaex + ADO[t]))
v40 = E^(14.543 - 4501.5/TEMP)*scaleprm*(prmKA1*R1P[t] - prmKA2*R5P[t])
v41 = (E^(14.543 - 4501.5/TEMP)*scaleprppsyn*VMPRPPSYN*VRPRPPSYN*(-((AMP[t]*PRPP[t])/KEQPRPPSYN) + ATP[t]*R5P[t]))/((KPRPPPRPPSYN*VMPRPPSYN*AMP[t])/KEQPRPPSYN + KR5PPRPPSYN*VRPRPPSYN*ATP[t] + (KAMPPRPPSYN*VMPRPPSYN*PRPP[t])/KEQPRPPSYN + (VMPRPPSYN*AMP[t]*PRPP[t])/KEQPRPPSYN + KATPPRPPSYN*VRPRPPSYN*R5P[t] + VRPRPPSYN*ATP[t]*R5P[t])
v42 = (0.5*E^(40.7298 - 12626.2626/TEMP)*FATP*NUMPUMP*scalepump)/DENPUMP
v43 = (ru5iVF*scaleru5pi*(-(R5P[t]/ru5iKEQ) + RU5P[t]))/(ru5iK1 + (ru5iK1*R5P[t])/ru5iK2 + RU5P[t])
v44 = (ETALD*scaletald*(-(NALD2*E4P[t]*F6P[t]) + NALD1*GAP[t]*S7P[t]))/DENTALD
v45 = (scaletki*tkiET*(-(tkiN2*GAP[t]*S7P[t]) + tkiN1*R5P[t]*X5P[t]))/tkiD
v46 = (ETTKII*scaletkii*(-(NTKII2*F6P[t]*GAP[t]) + NTKII1*E4P[t]*X5P[t]))/DENTKII
v47 = E^(18.475 - 5726.8/TEMP)*scaletpi*(tpiKTP1*DHAP[t] - tpiKTP2*GAP[t])
v48 = (scalexu5pe*xu5eVF*(RU5P[t] - X5P[t]/xu5eKEQ))/(xu5eK1 + RU5P[t] + (xu5eK1*X5P[t])/xu5eK2)
v5 = (adprtVM*scaleadprt*ADE[t]*PRPP[t])/((adprtKADE + ADE[t])*(adprtKPRPP + PRPP[t]))
v6 = (akVM*scaleak*ADO[t]*ATP[t])/((akKA + ADO[t])*(akKMGATP + ATP[t]))
v7 = (aldET*scaleald*(aldN1*FDP[t] - aldN2*DHAP[t]*GAP[t]))/(aldD1 + aldD4*DHAP[t] + aldD2*FDP[t] + aldD3*GAP[t] + aldD6*DHAP[t]*GAP[t] + aldD5*FDP[t]*GAP[t])
v8 = E^(19.876 - 6152.2/TEMP)*KAMPAS*scaleampase*AMP[t]
v9 = (ampdaVM*E^(19.876 - 6152.2/TEMP)*scaleampda*AMP[t])/(ampdaKAMP + AMP[t])
end rate equations

begin parameters
CO2 = 1.2
CONSTHK = 1.55
ETALD = 0.00069
ETG16PD = "2.1e-06"
ETG6PD = "9.3e-08"
ETTKII = 0.00033
FMAXN = 2.8166
G16P = 0.16
GLC = 5.0
GSSGT = 3.32
KALD1 = 2088000.0
KALD2 = 163080.0
KALD3 = 58680.0
KALD4 = 3636000.0
KALD5 = 1764000.0
KALD6 = 216000.0
KALD7 = 61200.0
KALD8 = 284400.0
KAMPPRPPSYN = 0.275
KATPPRPPSYN = 0.17
KE = 5.0
KEN = 1.708543057
KEN1 = 50000.0
KEQALD = 1.05
KEQPRPPSYN = 28.6
KEQTKII = 10.3
KG16P = 0.334
KG16PD1 = 1200000.0
KG16PD10 = 450000.0
KG16PD11 = 300.0
KG16PD12 = 9900000.0
KG16PD2 = 410.0
KG16PD3 = 1000000000.0
KG16PD4 = 26000.0
KG16PD5 = 48.0
KG16PD6 = 30.0
KG16PD7 = 630.0
KG16PD8 = 36000.0
KG16PD9 = 800.0
KG6P = 0.334
KG6PD1 = 110000000.0
KG6PD10 = 1400000000.0
KG6PD2 = 870.0
KG6PD3 = 26000000.0
KG6PD4 = 300.0
KG6PD5 = 750.0
KG6PD6 = 2000.0
KG6PD7 = 220000.0
KG6PD8 = 1100000000.0
KG6PD9 = 10000.0
KGA1 = 115975.0
KGAPDH = "1.1294e-06"
KGLC = 0.04
KINO = 25.0
KLAC = 12.8
KLD1 = 308642.2
KLDH = 44.64
KMAGLC = 0.051
KMG = 1.03
KMGAMG = 1.14
KMGATP = 1.13
KMGPG = 3.44
KNAPUMP = 6.2672
KPG = 1800.0
KPG23 = 2.7
KPGI = 0.41
KPGM = 0.1470588235
KPGM1 = 5000.0
KPRPPPRPPSYN = 0.09
KPUMP1 = 0.0082
KPUMP2 = 0.0501
KPYR = 63.8
KR5PPRPPSYN = 0.65
KSPPUMP = 0.1328
KSPUMP = 0.0617
KTKII1 = 777600.0
KTKII2 = 136800.0
KTKII3 = 122400.0
KTKII4 = 561600.0
KTKII5 = 806400.0
KTKII6 = 630000.0
KTKII7 = 144000.0
KTKII8 = 76680.0
Kglcg6p = 0.069
MGT = 2.7
NADHT = 0.089
NADPHT = 0.066
NAE = 140.0
PHE = 7.4
PI = 1.2
PNA = "3.96e-07"
PSIPUMP = 0.7114
RT = 0.6
S = "1.55e-06"
TEMP = 310.0
V = "8.7e-11"
adaKA = 0.052
adaKII = 1.3
adkK1 = 0.081
adkK2 = 0.81
adkK3 = 22.2
adkk0 = 0.22933781
adprtKADE = 0.0023
adprtKPRPP = 0.0195
akKA = 0.0004
akKMGATP = 0.8
aldET = 0.00037
aldK1 = 38500000.0
aldK2 = 838800.0
aldK3 = 6840000.0
aldK4 = 40320000.0
aldK5 = 252000.0
aldK6 = 23076000.0
aldkeq = 0.081
ampdaKAMP = 0.8
ampdaKIMP = 2.5
ampdaVM = 0.01
dpgaseVM = 0.52
gssgET = "1.25e-07"
gssgK1 = 85000000.0
gssgK10 = 50000000.0
gssgK11 = 7000.0
gssgK12 = 100000000.0
gssgK2 = 510.0
gssgK3 = 100000000.0
gssgK4 = 7200.0
gssgK5 = 810.0
gssgK6 = 1000.0
gssgK7 = 1000000.0
gssgK8 = 50000000.0
gssgK9 = 1000000.0
hgprtKHX = 0.22
hgprtKPRPP = 0.005
hgprtVM = 0.2011
hxexKM = 0.4
hxexPM = 37.8
hxexVM = 151.6
kAss = 10000.0
kFMAXK = 3.115
kKM = 4.0
kPK = "3.564e-07"
kmade = 2.6
kmaex = 0.13
kmino = 0.13
kpgi1 = 60290.0
kpgk1 = 8633310.0
mgadpk = 0.81
mgampk = 22.2
mgatpk = 0.081
mgdpgk = 1.67
naKM = 21.0
pfkKAMP = 0.033
pfkKF6P = 0.1
pfkKMG = 0.44
pfkKMGATP = 0.068
pfkLO = 0.001072
pglaseKM = 0.08
pglaseVM = 9.93
pkKADP = 0.474
pkKATP = 3.39
pkKFDP = 0.005
pkKPEP = 0.225
pkL0 = 19.0
pnpaseKA1 = 1950.0
pnpaseKeq = 0.09
prmKA1 = 427.59
prmKeq = 13.3
ru5iK1 = 0.78
ru5iK2 = 2.2
ru5iKEQ = 2.57
ru5iVF = 1702.296
scaleada = 1.0
scaleadeex = 1.0
scaleadk = 1.0
scaleadoex = 1.0
scaleadprt = 1.0
scaleak = 1.0
scaleald = 1.0
scaleampase = 1.0
scaleampda = 1.0
scaleatpase = 1.0
scaledpgase = 1.0
scaledpgm = 1.0
scaleen = 1.0
scaleg6pd = 1.0
scalegapdh = 1.0
scalegl6pd = 1.0
scalegshr = 1.0
scalegssgr = 1.0
scalehgprt = 1.0
scalehk = 1.0
scalehxex = 1.0
scaleimpase = 1.0
scalekleak = 1.0
scaleldh = 1.0
scalelex = 1.0
scalenaleak = 1.0
scalepex = 1.0
scalepfk = 1.0
scalepgi = 1.0
scalepgk = 1.0
scalepglase = 1.0
scalepgm = 1.0
scalepk = 1.0
scalepnpase = 1.0
scaleprm = 1.0
scaleprppsyn = 1.0
scalepump = 1.0
scaleru5pi = 1.0
scaletald = 1.0
scaletki = 1.0
scaletkii = 1.0
scaletpi = 1.0
scalexu5pe = 1.0
tkiET = 0.00033
tkiK1 = 777600.0
tkiK2 = 136800.0
tkiK3 = 122400.0
tkiK4 = 561600.0
tkiK5 = 1184400.0
tkiK6 = 630000.0
tkiK7 = 144000.0
tkiK8 = 161280.0
tkiKEQ = 1.2
tpiEPSI = 200.0
tpiKTP1 = 31056.0
tpiKTPI = 22.0
vmade = 90.0
vmaex = 61.2
vmino = 61.2
xu5eK1 = 0.19
xu5eK2 = 1.5
xu5eKEQ = 3.0
xu5eVF = 4633.56
GSHex = 0.0
HXex = 0.0
KLEak = 0.0
LACEXT = 1.2
NALEak = 0.0
PYREXT = 0.06
adeex = 0.0013
adoex = "8e-05"
inoex = 0.0002
defaultcompartment = 1.0
end parameters

begin initial conditions
ADE[0] = ADEi
ADO[0] = ADOi
ADP[0] = ADPi
AMP[0] = AMPi
ATP[0] = ATPi
DHAP[0] = DHAPi
DPG13[0] = DPG13i
DPG23[0] = DPG23i
E4P[0] = E4Pi
F6P[0] = F6Pi
FDP[0] = FDPi
G6P[0] = G6Pi
GAP[0] = GAPi
GL6P[0] = GL6Pi
GO6P[0] = GO6Pi
GSH[0] = GSHi
HX[0] = HXi
IMP[0] = IMPi
INO[0] = INOi
KI[0] = KIi
LAC[0] = LACi
MG[0] = MGi
MGADP[0] = MGADPi
MGAMP[0] = MGAMPi
MGATP[0] = MGATPi
MGDPG23[0] = MGDPG23i
NADH[0] = NADHi
NADPH[0] = NADPHi
NAI[0] = NAIi
PEP[0] = PEPi
PG2[0] = PG2i
PG3[0] = PG3i
PRPP[0] = PRPPi
PYR[0] = PYRi
R1P[0] = R1Pi
R5P[0] = R5Pi
RU5P[0] = RU5Pi
S7P[0] = S7Pi
X5P[0] = X5Pi
end initial conditions

begin initial values
ADEi = 0.00103
ADOi = 5.19e-05
ADPi = 0.492
AMPi = 0.148
ATPi = 2.401
DHAPi = 0.147
DPG13i = 0.000228
DPG23i = 4.65
E4Pi = 0.00507
F6Pi = 0.0272
FDPi = 0.01386
G6Pi = 0.066
GAPi = 0.00669
GL6Pi = 0.001759
GO6Pi = 0.0374
GSHi = 3.312
HXi = 0.000226
IMPi = 0.0192
INOi = 4.27e-05
KIi = 139.882
LACi = 1.398
MGi = 0.233836
MGADPi = 0.110216
MGAMPi = 0.00154265
MGATPi = 1.78328
MGDPG23i = 0.571129
NADHi = 0.030592
NADPHi = 0.0659
NAIi = 14.135
PEPi = 0.02087
PG2i = 0.01226
PG3i = 0.0839
PRPPi = 0.0136
PYRi = 0.062
R1Pi = 0.0117
R5Pi = 0.01266
RU5Pi = 0.00493
S7Pi = 0.02398
X5Pi = 0.01478
end initial values

begin assignment rules
NADP := NADPHT - NADPH[t]
HE := 10^(3 - PHE)
NAD := NADHT - NADH[t]
HI := (25*HE)/(23*RT)
pfkKATP := E^(-0.058693 + 0.32796*PHI)/1000
kpgi2 := kpgi1/KPGI
aldN1 := aldK1*aldK3*aldK5
DTKII5 := KTKII1*KTKII5*(KTKII3 + KTKII7)
DTKII6 := KTKII4*(KTKII2 + KTKII6)*KTKII8
aldD1 := (aldK2 + aldK3)*aldK5
aldD5 := aldK1*aldK4
aldD6 := aldK4*aldK6
tpiKTP2 := tpiKTP1*tpiKTPI
aldD3 := aldK2*aldK4
DALD3 := KALD2*KALD4*(KALD6 + KALD7)
DALD5 := KALD1*KALD5*(KALD3 + KALD7)
DALD7 := (KALD2 + KALD3)*KALD5*KALD8
DALD2 := (KALD2 + KALD3)*KALD5*KALD7
dpgasek := 0.5208*E^(-8.0287 + 0.98217*PHI)
DENTKII := DTKII2*E4P[t] + DTKII4*F6P[t] + DTKII7*E4P[t]*F6P[t] + DTKII3*GAP[t] + DTKII6*F6P[t]*GAP[t] + DTKII1*X5P[t] + DTKII5*E4P[t]*X5P[t] + DTKII8*GAP[t]*X5P[t]
NALD1 := KALD1*KALD3*KALD5*KALD7
KEN2 := KEN1/KEN
prmKA2 := prmKA1/prmKeq
pkDEN3 := (1 + FDP[t]/pkKFDP)^4*(1 + PEP[t]/pkKPEP)^4
KLD2 := KLD1/KLDH
akVM := (8.6229*E^(19.876 - 6152.2/TEMP))/(1 + 0.0001636/HI)
adkKAP1 := 54.400000000000006
adkKAP2 := adkKAP1/adkKADK
NALD2 := NALD1/KEQALD
DENG6PD1 := DG6PD1 + (DG6PD2*NADP + DG6PD3*G6P[t] + DG6PD4*GL6P[t] + DG6PD5*NADPH[t])/1000
KPGM2 := KPGM1/KPGM
katpase := (0.4149938*E^(14.543 - 4501.5/TEMP))/(1 + 8.271399999999999*^-7/HI + 267.7949761662471*HI)
MG6PD2 := (1.1631*E^(14.167 - 4380.4/TEMP)*ETG6PD*KG6PD10*KG6PD2*KG6PD4*KG6PD6*KG6PD8)/(1 + 2585.983966899405*HI)
DG6PD1 := KG6PD2*(KG6PD4*KG6PD6 + KG6PD4*KG6PD7 + KG6PD5*KG6PD7)*KG6PD9
DG6PD5 := KG6PD10*KG6PD2*(KG6PD4*KG6PD6 + KG6PD4*KG6PD7 + KG6PD5*KG6PD7)
DG6PD3 := KG6PD3*KG6PD5*KG6PD7*KG6PD9
DG6PD4 := KG6PD2*KG6PD4*KG6PD6*KG6PD8
DG6PD7 := KG6PD1*KG6PD4*KG6PD6*KG6PD8
DG6PD8 := KG6PD10*KG6PD3*KG6PD5*KG6PD7
pkDEN2 := (1 + ATP[t]/pkKATP)^4
DG6PD11 := KG6PD10*KG6PD3*(KG6PD5 + KG6PD6)*KG6PD8
DENG6PD3 := (DG6PD10*NADP*G6P[t]*GL6P[t])/1000000000 + (DG6PD9*GL6P[t]*NADPH[t])/1000000
DENG6PDH := DENG6PD1 + DENG6PD2 + DENG6PD3 + DENG6PD4
MG16PD2 := E^(17.896 - 5546.8/TEMP)*ETG16PD*KG16PD10*KG16PD12*KG16PD2*KG16PD4*KG16PD6*KG16PD8
DG6PD10 := KG6PD1*KG6PD3*(KG6PD5 + KG6PD6)*KG6PD8
DG16PD3 := KG16PD11*KG16PD3*KG16PD5*KG16PD7*KG16PD9
gssgDEN3 := ((gssgD10 + gssgD11)*NADP*GSH[t])/1000000 + (gssgD9*GSH[t]^2)/1000000
gssgDEN5 := (GSSG*gssgD15*NADP*GSH[t] + gssgD14*GSH[t]^2*NADPH[t])/1000000000
gssgDEN7 := (GSSG*gssgD18*NADP*GSH[t]^2)/1000000000000
DG16PD4 := KG16PD11*KG16PD2*KG16PD4*KG16PD6*KG16PD8
AG16PD1 := KG16PD1*KG16PD3*(KG16PD11*KG16PD5*KG16PD9 + KG16PD11*KG16PD6*KG16PD9)
DG16PD6 := AG16PD1 + BG16PD1 + CG16PD1
DG16PD18 := KG16PD10*KG16PD12*KG16PD3*(KG16PD5 + KG16PD6)*KG16PD8
gssgD12 := gssgK1*gssgK11*gssgK3*(gssgK5 + gssgK6)*gssgK8
gssgD13 := gssgK1*gssgK10*gssgK3*gssgK5*gssgK7
gshconst := 0.0011429*E^(5.0476 - 1567.2/TEMP)
tkiN1 := tkiK1*tkiK3*tkiK5*tkiK7
tkiN2 := tkiN1/tkiKEQ
tkiD1 := tkiK1*tkiK3*(tkiK6 + tkiK7)
tkiD3 := tkiK2*tkiK4*(tkiK6 + tkiK7)
tkiD2 := (tkiK2 + tkiK3)*tkiK5*tkiK7
tkiD4 := (tkiK2 + tkiK3)*tkiK6*tkiK8
tkiD6 := tkiK4*(tkiK2 + tkiK6)*tkiK8
tkiD8 := tkiK1*tkiK4*(tkiK6 + tkiK7)
NTKII1 := KTKII1*KTKII3*KTKII5*KTKII7
DG16PD17 := KG16PD1*KG16PD10*KG16PD3*(KG16PD5 + KG16PD6)*KG16PD8
NTKII2 := NTKII1/KEQTKII
DTKII1 := KTKII1*KTKII3*(KTKII6 + KTKII7)
DTKII3 := KTKII2*KTKII4*(KTKII6 + KTKII7)
DTKII4 := (KTKII2 + KTKII3)*KTKII6*KTKII8
DENG16PD4 := (CO2*DG16PD12*NADP*GO6P[t] + DG16PD13*NADP*GO6P[t]*RU5P[t])/1000000000
DENG16PD6A := (CO2*DG16PD16*NADPH[t]*RU5P[t])/1000000000
DENG16PD6B := (CO2*DG16PD17*NADP*GO6P[t]*RU5P[t])/1000000000000
DENG16PD7 := (CO2*DG16PD18*GO6P[t]*NADPH[t]*RU5P[t])/1000000000000
gssgM2 := 15455.6*E^(15.136 - 4690.8/TEMP)*gssgET*gssgK10*gssgK12*gssgK2*gssgK4*gssgK6*gssgK8
gssgD2 := gssgK1*gssgK11*(gssgK4*gssgK6 + gssgK4*gssgK7 + gssgK5*gssgK7)*gssgK9
gssgD3 := gssgK11*gssgK3*gssgK5*gssgK7*gssgK9
gssgD4 := gssgK11*gssgK2*gssgK4*gssgK6*gssgK8
DG16PD15 := KG16PD10*KG16PD12*KG16PD3*KG16PD5*KG16PD7
DTKII7 := (KTKII2 + KTKII3)*KTKII5*KTKII8
DTKII8 := KTKII1*KTKII4*(KTKII6 + KTKII7)
naV2 := FMAXN*(NAE/(NAE + naKM) - (RT*NAI[t])/(naKM + RT*NAI[t]))
gssgD7 := gssgK1*gssgK11*gssgK4*gssgK6*gssgK8
gssgD8 := gssgK12*gssgK3*gssgK5*gssgK7*gssgK9
VMPUMP := 2.318/(1 + 8.271399999999999*^-7/HI + 267.7949761662471*HI)
NUMPUMP := (KE^2 + 0.5*KE*KSPPUMP*PSIPUMP)*VMPUMP
DENPUMP1 := KE^2 + 2*KE*KSPPUMP + KSPPUMP*KSPUMP
DENPUMP2 := (1 + KNAPUMP/NAI[t])^3
DENPUMP3 := KPUMP1*KSPPUMP*KSPUMP + KPUMP2*(KE^2 + KE*KSPPUMP*PSIPUMP)
DENPUMP := DENPUMP1 + DENPUMP2*DENPUMP3
FATP := ATP[t]/(0.04 + ATP[t])
DALD1 := KALD1*KALD3*(KALD6 + KALD7)
DG16PD10 := KG16PD12*KG16PD2*KG16PD4*KG16PD6*KG16PD8
gssgD10 := gssgK12*gssgK2*gssgK4*gssgK6*gssgK8
pnpaseKA2 := pnpaseKA1/pnpaseKeq
DG16PD12 := KG16PD1*KG16PD11*KG16PD3*(KG16PD5 + KG16PD6)*KG16PD8
DALD4 := (KALD2 + KALD3)*KALD6*KALD8
DALD6 := KALD4*(KALD2 + KALD6)*KALD8
DENG16PD1 := DG16PD1 + (CO2*DG16PD4 + DG16PD2*NADP + DG16PD3*GO6P[t] + DG16PD5*NADPH[t])/1000
DENG16PD3 := (CO2*DG16PD10*NADPH[t] + CO2*DG16PD9*RU5P[t] + DG16PD11*NADPH[t]*RU5P[t])/1000000
gssgD15 := gssgK10*gssgK12*gssgK3*gssgK5*gssgK7
gssgD17 := gssgK1*gssgK10*gssgK3*(gssgK5 + gssgK6)*gssgK8
gssgD18 := gssgK10*gssgK12*gssgK3*(gssgK5 + gssgK6)*gssgK8
gssgDEN2 := (GSSG*gssgD8*NADP + GSSG*gssgD6*NADPH[t] + gssgD7*GSH[t]*NADPH[t])/1000000
DALD8 := KALD1*KALD4*(KALD6 + KALD7)
DENTALD := DALD3*E4P[t] + DALD4*F6P[t] + DALD6*E4P[t]*F6P[t] + DALD2*GAP[t] + DALD7*F6P[t]*GAP[t] + DALD1*S7P[t] + DALD8*E4P[t]*S7P[t] + DALD5*GAP[t]*S7P[t]
KIASE := 0.09387899999999999/(1 + 1.2158999999999999*^-6/HI + 389.9851805631386*HI)
VMPRPPSYN := 0.83633592/(1 + 1.8834999999999998*^-6/HI + 9901.970492127934*HI)
VRPRPPSYN := 0.13333333333333333*VMPRPPSYN
kV1 := (kPK*SV*(KE - RT*KI[t])*Log[RT])/(-1 + RT)
kV2 := kFMAXK*(KE/(KE + kKM) - (RT*KI[t])/(kKM + RT*KI[t]))
naV1 := (PNA*SV*Log[RT]*(NAE - RT*NAI[t]))/(-1 + RT)
SV := S/V
GSSG := 0.5*(GSSGT - GSH[t])
ATPF := ATP[t] - MGATP[t]
VM1HK := (7.5095*E^(18.21 - 5637.9/TEMP))/(1 + 4.6015*^-6/HI + 4401.795932740558*HI)
N1HK := 1 + GLC/KGLC + ((1 + GLC/KGLC)*MG[t])/KMG + (G16P/KG16P + (CONSTHK*GLC)/KGLC + G6P[t]/KG6P + (GLC*G6P[t])/(KGLC*Kglcg6p))*(1 + MG[t]/KMG) + (1 + GLC/KGLC)*(DPG23[t]/KPG23 + (DPG23[t]*MG[t])/(KMG*KMGPG)) + ((1 + GLC/KMAGLC)*MGATP[t])/KMGATP + ((1 + GLC/KMAGLC)*MG[t]*MGATP[t])/(KMGAMG*KMGATP)
pfkVM := 250*E^(13.474 - 4174.8/TEMP)
LPFK := ((1 + ATPF/pfkKATP)^4*pfkLO*(1 + MG[t]/pfkKMG)^4)/((1 + AMP[t]/pfkKAMP)^4*(1 + F6P[t]/pfkKF6P)^4)
aldN2 := aldN1/aldkeq
aldD2 := aldK1*(aldK3 + aldK5)
KGA2 := KGA1/KGAPDH
kpgk2 := kpgk1/KPG
dpgmk := (8889*E^(-144503/10000 + (37*PHI)/25))/5000
adaVM := (36.348*E^(19.876 - 6152.2/TEMP))/(1 + 0.00003162/HI + 5012.531328320802*HI)
DG6PD9 := KG6PD10*(KG6PD2*KG6PD4 + KG6PD2*KG6PD5 + KG6PD2*KG6PD6 + KG6PD4*KG6PD6)*KG6PD8
DG16PD2 := KG16PD1*KG16PD11*(KG16PD4*KG16PD6 + KG16PD4*KG16PD7 + KG16PD5*KG16PD7)*KG16PD9
CG16PD1 := KG16PD1*KG16PD11*KG16PD3*KG16PD5*KG16PD7
DG16PD8 := KG16PD12*KG16PD3*KG16PD5*KG16PD7*KG16PD9
DG16PD11 := KG16PD10*KG16PD12*KG16PD2*(KG16PD4*KG16PD6 + KG16PD4*KG16PD7 + KG16PD5*KG16PD7)
DG16PD13 := KG16PD1*KG16PD10*KG16PD3*KG16PD5*KG16PD7
DG16PD14 := KG16PD1*KG16PD10*KG16PD4*KG16PD6*KG16PD8
gssgD6 := gssgK1*gssgK3*(gssgK11*gssgK5*gssgK7 + gssgK11*gssgK5*gssgK9 + gssgK11*gssgK6*gssgK9 + gssgK11*gssgK7*gssgK9 + gssgK5*gssgK7*gssgK9)
PHI := 3 - Log[HI]/Log[10]
aldD4 := (aldK2 + aldK3)*aldK6
pkDEN1 := (MGADP[t]*PEP[t])/((pkKADP + MGADP[t])*(pkKPEP + PEP[t]))
VLAC := -(KLAC*(LACEXT - LAC[t]))
KAMPAS := 1.648/(1 + 1.2158999999999999*^-6/HI + 389.9851805631386*HI)
adkKADK := (adkk0*(adkK1 + MG[t])*(adkK3 + MG[t]))/(adkK2 + MG[t])^2
adprtVM := (0.202*E^(23.356 - 7240.2/TEMP))/(1 + 3.162*^-7/HI + 25119.316754584277*HI)
MG6PD1 := (1.1631*E^(14.167 - 4380.4/TEMP)*ETG6PD*KG6PD1*KG6PD3*KG6PD5*KG6PD7*KG6PD9)/(1 + 2585.983966899405*HI)
DG6PD2 := KG6PD1*(KG6PD4*KG6PD6 + KG6PD4*KG6PD7 + KG6PD5*KG6PD7)*KG6PD9
DG6PD6 := KG6PD1*KG6PD3*(KG6PD5*KG6PD7 + KG6PD5*KG6PD9 + KG6PD6*KG6PD9 + KG6PD7*KG6PD9)
NUMG6PDH := (MG6PD1*NADP*G6P[t] - MG6PD2*GL6P[t]*NADPH[t])/1000000
DENG6PD2 := (DG6PD6*NADP*G6P[t] + DG6PD7*NADP*GL6P[t] + DG6PD8*G6P[t]*NADPH[t])/1000000
DENG6PD4 := (DG6PD11*G6P[t]*GL6P[t]*NADPH[t])/1000000000
MG16PD1 := E^(17.896 - 5546.8/TEMP)*ETG16PD*KG16PD1*KG16PD11*KG16PD3*KG16PD5*KG16PD7*KG16PD9
DG16PD1 := KG16PD11*KG16PD2*(KG16PD4*KG16PD6 + KG16PD4*KG16PD7 + KG16PD5*KG16PD7)*KG16PD9
DG16PD5 := KG16PD12*KG16PD2*(KG16PD4*KG16PD6 + KG16PD4*KG16PD7 + KG16PD5*KG16PD7)*KG16PD9
BG16PD1 := KG16PD1*KG16PD3*(KG16PD11*KG16PD7*KG16PD9 + KG16PD5*KG16PD7*KG16PD9)
DG16PD7 := KG16PD1*KG16PD11*KG16PD4*KG16PD6*KG16PD8
DG16PD9 := KG16PD10*KG16PD2*KG16PD4*KG16PD6*KG16PD8
DG16PD16 := KG16PD10*KG16PD12*(KG16PD2*KG16PD4 + KG16PD2*KG16PD5 + KG16PD2*KG16PD6 + KG16PD4*KG16PD6)*KG16PD8
NUMG16PD := 0.036706499999999996*MG16PD1*NADP*GO6P[t] - 0.0000367065*CO2*MG16PD2*NADPH[t]*RU5P[t]
DENG16PD2 := (CO2*DG16PD7*NADP + DG16PD6*NADP*GO6P[t] + DG16PD8*GO6P[t]*NADPH[t])/1000000
DENG16PD5 := (CO2*DG16PD14*NADP*RU5P[t] + DG16PD15*GO6P[t]*NADPH[t]*RU5P[t])/1000000000
DENG16PD6 := DENG16PD6A + DENG16PD6B
DENG16PD := DENG16PD1 + DENG16PD2 + DENG16PD3 + DENG16PD4 + DENG16PD5 + DENG16PD6 + DENG16PD7
gssgM1 := 15455.6*E^(15.136 - 4690.8/TEMP)*gssgET*gssgK1*gssgK11*gssgK3*gssgK5*gssgK7*gssgK9
gssgD1 := gssgK11*gssgK2*(gssgK4*gssgK6 + gssgK4*gssgK7 + gssgK5*gssgK7)*gssgK9
gssgD5 := gssgK12*gssgK2*(gssgK4*gssgK6 + gssgK4*gssgK7 + gssgK5*gssgK7)*gssgK9
gssgD9 := gssgK10*gssgK2*gssgK4*gssgK6*gssgK8
gssgD11 := gssgK10*gssgK12*gssgK2*(gssgK4*gssgK6 + gssgK4*gssgK7 + gssgK5*gssgK7)
gssgD14 := gssgK1*gssgK10*gssgK4*gssgK6*gssgK8
gssgD16 := gssgK10*gssgK12*(gssgK2*gssgK4 + gssgK2*gssgK5 + gssgK2*gssgK6 + gssgK4*gssgK6)*gssgK8
gssgNUM := -(gssgM2*NADP*GSH[t]^2) + GSSG*gssgM1*NADPH[t]
gssgDEN1 := gssgD1 + (GSSG*gssgD3 + gssgD5*NADP + gssgD4*GSH[t] + gssgD2*NADPH[t])/1000
gssgDEN4 := (GSSG*(gssgD12 + gssgD13)*GSH[t]*NADPH[t])/1000000000
gssgDEN6 := (gssgD16*NADP*GSH[t]^2)/1000000000 + (GSSG*gssgD17*GSH[t]^2*NADPH[t])/1000000000000
gssgDEN := gssgDEN1 + gssgDEN2 + gssgDEN3 + gssgDEN4 + gssgDEN5 + gssgDEN6 + gssgDEN7
tkiD5 := tkiK1*tkiK5*(tkiK3 + tkiK7)
tkiD7 := (tkiK2 + tkiK3)*tkiK5*tkiK8
tkiD := tkiD3*GAP[t] + tkiD2*R5P[t] + tkiD4*S7P[t] + tkiD6*GAP[t]*S7P[t] + tkiD7*R5P[t]*S7P[t] + tkiD1*X5P[t] + tkiD8*GAP[t]*X5P[t] + tkiD5*R5P[t]*X5P[t]
DTKII2 := (KTKII2 + KTKII3)*KTKII5*KTKII7
KDPGM := Piecewise[{{0, PHI < 6.85}}, 275000.*E^(16.907 - 5236.7/TEMP)]
end assignment rules

begin function definitions
pkVM[PHI_,TEMP_] := Piecewise[{{(5*E^(20.61 - 6379.6/TEMP)*(300 - 90.9*(-7.2 + PHI)))/6, PHI > 7.21}}, (5*E^(20.61 - 6379.6/TEMP)*(300 - 113.6*(7.2 - PHI)))/6]
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units