begin name
arnold9
end name

begin reactions
ATPS {1.0}$EXTERNAL = {1.0}ATP
F6PTK {1.0}GAP + {1.0}F6P = {1.0}X5P + {1.0}E4P
FBPA {1.0}DHAP + {1.0}GAP = {1.0}FBP
FBPase {1.0}FBP = {1.0}F6P
GAPDH {1.0}DPGA + {1.0}$NADPH + {1.0}$H = {1.0}GAP + {1.0}$NADP
PGAK {1.0}PGA + {1.0}ATP = {1.0}DPGA
PGI {1.0}F6P = {1.0}G6P
PGM {1.0}G6P = {1.0}G1P
R5PI {1.0}R5P = {1.0}Ru5P
Ru5PE {1.0}X5P = {1.0}Ru5P
Ru5PK {1.0}Ru5P + {1.0}ATP = {1.0}RuBP
RuBisCO {1.0}RuBP = {2.0}PGA
S7PTK {1.0}GAP + {1.0}S7P = {1.0}X5P + {1.0}R5P
SBPA {1.0}DHAP + {1.0}E4P = {1.0}SBP
SBPase {1.0}SBP = {1.0}S7P
StarchP {1.0}$EXTERNAL = {1.0}G1P
StarchS {1.0}G1P + {1.0}ATP = {1.0}$EXTERNAL
TPTDHAP {1.0}DHAP = {1.0}$EXTERNAL
TPTGAP {1.0}GAP = {1.0}$EXTERNAL
TPTPGA {1.0}PGA = {1.0}$EXTERNAL
TPI {1.0}GAP = {1.0}DHAP
end reactions

begin rate equations
ATPS = chloroplast*function6[ATPSVm, ADP, Pi, ATPSK1, ATPSK2]
F6PTK = chloroplast*(F6PTKk1*F6P[t]*GAP[t] - F6PTKk2*E4P[t]*X5P[t])
FBPA = chloroplast*(-(FBPAk2*FBP[t]) + FBPAk1*DHAP[t]*GAP[t])
FBPase = chloroplast*function2[FBPaseVm, FBP[t], FBPaseK, F6P[t], FBPaseKR1, Pi, FBPaseKR2]
GAPDH = chloroplast*(GAPDHk1*H*NADPH*DPGA[t] - GAPDHk2*NADP*Pi*GAP[t])
PGAK = chloroplast*(-(ADP*PGAKk2*DPGA[t]) + PGAKk1*ATP[t]*PGA[t])
PGI = chloroplast*(PGIk1*F6P[t] - PGIk2*G6P[t])
PGM = chloroplast*(-(PGMk2*G1P[t]) + PGMk1*G6P[t])
R5PI = chloroplast*(R5PIk1*R5P[t] - R5PIk2*Ru5P[t])
Ru5PE = chloroplast*(-(Ru5PEk2*Ru5P[t]) + Ru5PEk1*X5P[t])
Ru5PK = chloroplast*function4[Ru5PKVm, Ru5P[t], ATP[t], Ru5PKK1, PGA[t], Ru5PKKR1, RuBP[t], Ru5PKKR2, Pi, Ru5PKKR3, ADP, Ru5PKKR41, Ru5PKK2, Ru5PKKR42]
RuBisCO = chloroplast*function1[RuBisCOVm, RuBP[t], RuBisCOK, PGA[t], RuBisCOKR1, FBP[t], RuBisCOKR2, SBP[t], RuBisCOKR3, Pi, RuBisCOKR4, NADPH, RuBisCOKR5]
S7PTK = chloroplast*(S7PTKk1*GAP[t]*S7P[t] - S7PTKk2*R5P[t]*X5P[t])
SBPA = chloroplast*(SBPAk1*DHAP[t]*E4P[t] - SBPAk2*SBP[t])
SBPase = chloroplast*function3[SBPaseVm, SBP[t], SBPaseK, Pi, SBPaseKR1]
StarchP = chloroplast*function3[StarchPVm, Pi, StarchPK, G1P[t], StarchPKR1]
StarchS = chloroplast*function5[StarchSVm, G1P[t], ATP[t], StarchSK1, StarchSK2, ADP, StarchSKR1, Pi, StarchSKA1, PGA[t], StarchSKA2, F6P[t], StarchSKA3, FBP[t]]
TPTDHAP = chloroplast*function7[TPTDHAPVm, DHAP[t], TPTDHAPKA, Pext, TPTDHAPK, Pi, TPTDHAPKR1, PGA[t], TPTDHAPKR2, GAP[t], TPTDHAPKR3]
TPTGAP = chloroplast*function7[TPTGAPVm, GAP[t], TPTGAPKA, Pext, TPTGAPK, Pi, TPTGAPKR1, PGA[t], TPTGAPKR2, DHAP[t], TPTGAPKR3]
TPTPGA = chloroplast*function7[TPTPGAVm, PGA[t], TPTPGAKA, Pext, TPTPGAK, Pi, TPTPGAKR1, GAP[t], TPTPGAKR2, DHAP[t], TPTPGAKR3]
TPI = chloroplast*(-(TPIk2*DHAP[t]) + TPIk1*GAP[t])
end rate equations

begin parameters
H = "1.25892541179417e-05"
NADP = 0.29
NADPH = 0.21
Pext = 0.5
RuBisCOVm = 340.0
RuBisCOK = 0.02
RuBisCOKR1 = 0.84
RuBisCOKR2 = 0.04
RuBisCOKR3 = 0.0075
RuBisCOKR4 = 0.9
RuBisCOKR5 = 0.07
PGAKk1 = 500000000.0
PGAKk2 = 1612900000000.0
GAPDHk1 = 500000000.0
GAPDHk2 = 31.25
TPIk1 = 500000000.0
TPIk2 = 22727000.0
FBPAk1 = 500000000.0
FBPAk2 = 70423000.0
FBPaseVm = 200.0
FBPaseK = 0.03
FBPaseKR1 = 0.7
FBPaseKR2 = 12.0
F6PTKk1 = 500000000.0
F6PTKk2 = 5952400000.0
SBPAk1 = 500000000.0
SBPAk2 = 38461500.0
SBPaseVm = 40.0
SBPaseK = 0.02
SBPaseKR1 = 12.0
S7PTKk1 = 500000000.0
S7PTKk2 = 588240000.0
R5PIk1 = 500000000.0
R5PIk2 = 1250000000.0
Ru5PEk1 = 500000000.0
Ru5PEk2 = 746269000.0
Ru5PKVm = 1000.0
Ru5PKK1 = 0.05
Ru5PKKR1 = 2.0
Ru5PKKR2 = 0.7
Ru5PKKR3 = 4.0
Ru5PKKR41 = 2.5
Ru5PKK2 = 0.05
Ru5PKKR42 = 0.4
PGIk1 = 500000000.0
PGIk2 = 217400000.0
PGMk1 = 500000000.0
PGMk2 = 8621000000.0
StarchSVm = 40.0
StarchSK1 = 0.08
StarchSK2 = 0.08
StarchSKR1 = 10.0
StarchSKA1 = 0.1
StarchSKA2 = 0.02
StarchSKA3 = 0.02
StarchPVm = 40.0
StarchPK = 0.1
StarchPKR1 = 0.05
ATPSVm = 350.0
ATPSK1 = 0.014
ATPSK2 = 0.3
TPTPGAVm = 250.0
TPTPGAKA = 0.74
TPTPGAK = 0.25
TPTPGAKR1 = 0.63
TPTPGAKR2 = 0.075
TPTPGAKR3 = 0.077
TPTGAPVm = 250.0
TPTGAPKA = 0.74
TPTGAPK = 0.075
TPTGAPKR1 = 0.63
TPTGAPKR2 = 0.25
TPTGAPKR3 = 0.077
TPTDHAPVm = 250.0
TPTDHAPKA = 0.74
TPTDHAPK = 0.077
TPTDHAPKR1 = 0.63
TPTDHAPKR2 = 0.25
TPTDHAPKR3 = 0.075
chloroplast = 1.0
cytosol = 1.0
end parameters

begin initial conditions
ATP[0] = ATPi
DHAP[0] = DHAPi
DPGA[0] = DPGAi
E4P[0] = E4Pi
F6P[0] = F6Pi
FBP[0] = FBPi
G1P[0] = G1Pi
G6P[0] = G6Pi
GAP[0] = GAPi
PGA[0] = PGAi
R5P[0] = R5Pi
Ru5P[0] = Ru5Pi
RuBP[0] = RuBPi
S7P[0] = S7Pi
SBP[0] = SBPi
X5P[0] = X5Pi
end initial conditions

begin initial values
ATPi = 0.68
DHAPi = 0.48
DPGAi = 0.0011
E4Pi = 0.05
F6Pi = 0.640764257004719
FBPi = 0.669999999999999
G1Pi = 0.0854779518844294
G6Pi = 1.47375779111085
GAPi = 0.02
PGAi = 2.4
R5Pi = 0.125186846038864
Ru5Pi = 0.0500747384155456
RuBPi = 2.0
S7Pi = 2.0
SBPi = 0.3
X5Pi = 0.0747384155455904
end initial values

begin assignment rules
ADP := 1.5 - ATP[t]
Pi := 15 - ATP[t] - DHAP[t] - E4P[t] - F6P[t] - G1P[t] - G6P[t] - GAP[t] - PGA[t] - R5P[t] - Ru5P[t] - S7P[t] - 2*(DPGA[t] + FBP[t] + RuBP[t] + SBP[t]) - X5P[t]
end assignment rules

begin function definitions
function5[Vm_,S1_,S2_,K1_,K2_,R1_,KR1_,R2_,KA1_,A1_,KA2_,A2_,KA3_,A3_] := (S1*S2*Vm)/((1 + R1/KR1)*(K1 + S1)*(K2*(1 + (K2*R2)/(A1*KA1 + A2*KA2 + A3*KA3)) + S2))
function3[Vm_,S_,K_,R1_,KR1_] := (S*Vm)/(K*(1 + R1/KR1) + S)
function1[Vm_,S_,K_,R1_,KR1_,R2_,KR2_,R3_,KR3_,R4_,KR4_,R5_,KR5_] := (S*Vm)/(K*(1 + R1/KR1 + R2/KR2 + R3/KR3 + R4/KR4 + R5/KR5) + S)
function6[Vm_,s1_,s2_,K1_,K2_] := (s1*s2*Vm)/((K1 + s1)*(K2 + s2))
function4[Vm_,S1_,S2_,K1_,R1_,KR1_,R2_,KR2_,R3_,KR3_,R4_,KR41_,K2_,KR42_] := (S1*S2*Vm)/((K1*(1 + R1/KR1 + R2/KR2 + R3/KR3) + S1)*(K2*(1 + R4/KR42) + (1 + R4/KR41)*S2))
function7[Vm_,S_,KA_,A_,K_,R1_,KR1_,R2_,KR2_,R3_,KR3_] := (S*Vm)/(K*(1 + (1 + KA/A)*(R1/KR1 + R2/KR2 + R3/KR3)) + (1 + KA/A)*S)
function2[Vm_,S_,K_,R1_,KR1_,R2_,KR2_] := (S*Vm)/(K*(1 + R1/KR1 + R2/KR2) + S)
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
G6P http://identifiers.org/pubchem.compound/439958
NADP http://identifiers.org/chebi/CHEBI:25523
NADPH http://identifiers.org/chebi/CHEBI:16474
end species annotations

begin units
end units