begin name
arnold10
end name

begin reactions
AGPase {1.0}HeP + {1.0}ATP = {1.0}ADPG + {1.0}$PiPi
ATPS {1.0}$EXTERNAL = {1.0}ATP
F26BPcP {1.0}F26BPc = {1.0}HePc
F6PTK {1.0}HeP + {1.0}TP = {1.0}E4P + {1.0}PeP
F6PcK {1.0}HePc = {1.0}F26BPc
FBPA {2.0}TP = {1.0}FBP
FBPase {1.0}FBP = {1.0}HeP
FBPcA {2.0}TPc = {1.0}FBPc
FBPcase {1.0}FBPc = {1.0}HePc
PGAK {1.0}PGA + {1.0}ATP = {1.0}TP
Ru5PK {1.0}PeP + {1.0}ATP = {1.0}RuBP
RuBisCO1 {1.0}RuBP = {1.0}ER
RuBisCO2CO2 {1.0}ER = {1.0}EPP
RuBisCO4 {1.0}EPP = {1.0}PGA + {1.0}EP
RuBisCO5EOP {1.0}EOP = {1.0}$EXTERNAL
RuBisCO5EP {1.0}EP = {1.0}PGA
RuBisCO6O2 {1.0}ER + {1.0}$O2 + {0.5}ATP = {1.0}EPG + {0.5}PGA
RuBisCO7 {1.0}EPG = {1.0}EP
S7PTK {1.0}S7P + {1.0}TP = {2.0}PeP
SBPA {1.0}E4P + {1.0}TP = {1.0}SBP
SBPase {1.0}SBP = {1.0}S7P
StS {1.0}ADPG = {1.0}$EXTERNAL
SucPcP {1.0}SucPc = {1.0}$Succ
SucPcS {1.0}HePc + {1.0}UDPGc = {1.0}$UDPc + {1.0}SucPc + {1.0}$Hc
TPTPGA {1.0}PGA = {1.0}PGAc
TPTTP {1.0}TP = {1.0}TPc
UGPase {1.0}HePc + {1.0}$UTPc = {1.0}UDPGc + {1.0}$PiPic
end reactions

begin rate equations
AGPase = chloroplast*function28[AGPaseVm, PGA[t], Pi, G1P, ATP[t], ADPG[t], PiPi, AGPaseq, AGPaseKs1, AGPaseKs2, AGPaseKp1, AGPaseKp2]
ATPS = chloroplast*function24[ADP, Pi, ATP[t], q12, ATPSKs1, ATPSKs2, Kp12, ATPSVm]
F26BPcP = cytosol*function33[F26BPcPVm, F26BPc[t], F26BPcPKs1, TPc[t], PGAc[t], F26BPcPKr1, Pic, F26BPcPKr3, HePc[t], F26BPcPKr4]
F6PTK = chloroplast*Transketolase[F6PTKVm, F6PTKq, F6P, GAP, E4P[t], X5P, F6PTKK1, F6PTKK2, F6PTKK1s2, F6PTKK2s1, S7P[t], F6PTKK2r1, F6PTKK2s2, R5P]
F6PcK = cytosol*function32[F6PcKVm, F6Pc, F6PcKKs1, Pic, F6PcKKr1, TPc[t], PGAc[t], F6PcKKr2]
FBPA = chloroplast*function24[GAP, DHAP, FBP[t], FBPAq, FBPAKs1, FBPAKs2, FBPAKp1, FBPAVm]
FBPase = chloroplast*function25[FBPaseVm, FBP[t], F6P, Pi, FBPaseq, FBPaseKs1, FBPaseKp1, FBPaseKp2]
FBPcA = cytosol*function24[GAPc, DHAPc, FBPc[t], FBPcAq, FBPcAKs1, FBPcAKs2, FBPcAKp1, FBPcAVm]
FBPcase = cytosol*function30[FBPcaseVm, FBPc[t], F6Pc, Pic, FBPcaseq, FBPcaseKs1, F26BPc[t], FBPcaseKr1, FBPcaseKp1, FBPcaseKp2]
PGAK = chloroplast*function23[PGAKVm, PGA[t], ATP[t], GAP, ADP, Pi, q1, PGAKKs1, PGAKKs2, PGAKKp1, PGAKKp2, PGAKKp3]
Ru5PK = chloroplast*function27[Ru5P, ATP[t], RuBP[t], ADP, Ru5PKq, Ru5PKKs1, Ru5PKKs2, Ru5PKKp1, Ru5PKKp2, Ru5PKVm]
RuBisCO1 = chloroplast*(-(RuBisCO1k2*ER[t]) + E*RuBisCO1k1*RuBP[t])
RuBisCO2CO2 = chloroplast*CO2*RuBisCO2CO2k1*ER[t]
RuBisCO4 = chloroplast*(RuBisCO4k1*EPP[t] - RuBisCO4k2*EP[t]*PGA[t])
RuBisCO5EOP = chloroplast*function34[W4]
RuBisCO5EP = chloroplast*(RuBisCO5EPk1*EP[t] - E*RuBisCO5EPk2*PGA[t])
RuBisCO6O2 = chloroplast*O2*RuBisCO6O2k1*ER[t]
RuBisCO7 = chloroplast*RuBisCO7k1*EPG[t]
S7PTK = chloroplast*Transketolase[S7PTKVm, S7PTKq, S7P[t], GAP, R5P, X5P, S7PTKK1, S7PTKK2, S7PTKK1s2, S7PTKK2s1, F6P, S7PTKK2r1, S7PTKK2s2, E4P[t]]
SBPA = chloroplast*function24[E4P[t], DHAP, SBP[t], SBPAq, SBPAKs1, SBPAKs2, SBPAKp1, SBPAVm]
SBPase = chloroplast*function26[SBPaseVm, SBP[t], S7P[t], Pi, SBPaseq, SBPaseKs1, SBPaseKp1, SBPaseKp2]
StS = chloroplast*function29[StSVm, ADPG[t], ADP, StSq, StSKs1, StSKp1]
SucPcP = cytosol*function25[SucPcPVm, SucPc[t], Succ, Pic, SucPcPq, SucPcPKs1, SucPcPKp1, SucPcPKp2]
SucPcS = cytosol*function31[SucPcSVm, F6Pc, UDPGc[t], UDPc, SucPc[t], Hc, SucPcSq, SucPcSKs1, Pic, SucPcSKr11, SucPcSKs2, SucPcSKp1, SucPcSKp2, SucPcSKr12]
TPTPGA = TPTout[TPTPGAVm, PGA[t], TPTPGAKs, TP[t], TPTPGAKr1, Pi, TPTPGAKr2, PGAc[t], TPTPGAKp, TPc[t], TPTPGAKr3, Pic, TPTPGAKr4]
TPTTP = TPTout[TPTTPVm, TP[t], TPTTPKs, PGA[t], TPTTPKr1, Pi, TPTTPKr2, TPc[t], TPTTPKp, PGAc[t], TPTTPKr3, Pic, TPTTPKr4]
UGPase = cytosol*function27[G1Pc, UTPc, UDPGc[t], PiPic, UGPaseq, UGPaseKs1, UGPaseKs2, UGPaseKp1, UGPaseKp2, UGPaseVm]
end rate equations

begin parameters
ADT = 0.0015
ADTc = 0.001
Et = 0.0028030303030303
Kp12 = 224014.808032967
NADPT = 0.0005
PiT = 0.0284090909090909
PiTc = 0.0170454545454545
UDTc = 0.001
V28 = "7.386364e-05"
q1 = 0.129053067280279
q12 = 2227862541257.35
ATPc = 0.00036
H = 0.0891250931577478
Hc = 0.158489318357816
NADP = 0.00029
NADPH = 0.00021
O2 = 0.00026
PiPi = "1e-06"
PiPic = "4e-05"
Succ = 0.0
UDPc = 0.00064
UTPc = 0.00036
RuBisCO1k1 = 50000.0
RuBisCO1k2 = 0.9
RuBisCO2CO2k1 = 300000.0
RuBisCO4k1 = 6.0
RuBisCO4k2 = 0.0
RuBisCO5EPk1 = 6.0
RuBisCO5EPk2 = 70000.0
RuBisCO6O2k1 = 3030.3
RuBisCO7k1 = 3.0
PGAKVm = 0.0170455
PGAKKs1 = 0.0011122
PGAKKs2 = 0.0003307
PGAKKp1 = 0.00027035
PGAKKp2 = 0.00053013
PGAKKp3 = 0.0027397
FBPAq = 1.18815
FBPAKs1 = 0.00027035
FBPAKs2 = 0.00036393
FBPAKp1 = "2.0129e-05"
FBPAVm = 0.022727
FBPaseVm = 0.011364
FBPaseq = 0.77294
FBPaseKs1 = "3.2842e-05"
FBPaseKp1 = "6.3429e-05"
FBPaseKp2 = 0.0017914
F6PTKVm = 0.170455
F6PTKq = 0.99943
F6PTKK1 = 0.00061349
F6PTKK2 = 0.00011438
F6PTKK1s2 = 0.00027035
F6PTKK2s1 = 0.0005407
F6PTKK2r1 = 0.00017677
F6PTKK2s2 = "9.0464e-05"
SBPAq = 1.18815
SBPAKs1 = 0.00017677
SBPAKs2 = 0.00036393
SBPAKp1 = "2.0129e-05"
SBPAVm = 0.011364
SBPaseVm = 0.00568182
SBPaseq = 0.77294
SBPaseKs1 = "1.2713e-05"
SBPaseKp1 = "1.5597e-05"
SBPaseKp2 = 0.006744
S7PTKVm = 0.0821023
S7PTKq = 0.99996
S7PTKK1 = 0.00061349
S7PTKK2 = 0.00011438
S7PTKK1s2 = 0.00027035
S7PTKK2s1 = 0.00017677
S7PTKK2r1 = 0.0005407
S7PTKK2s2 = "9.0464e-05"
Ru5PKq = 1.05289
Ru5PKKs1 = "3.63934e-05"
Ru5PKKs2 = 0.00055117
Ru5PKKp1 = "9.95868e-05"
Ru5PKKp2 = "9.11825e-05"
Ru5PKVm = 0.568182
ATPSKs1 = 0.00031808
ATPSKs2 = 0.00031612
ATPSVm = 0.0284091
AGPaseVm = 0.00113636
AGPaseq = 0.11059
AGPaseKs1 = 0.0010398
AGPaseKs2 = 0.00011023
AGPaseKp1 = 0.00053013
AGPaseKp2 = 0.01951
StSVm = 0.00284091
StSq = 1.00326
StSKs1 = 0.000212052
StSKp1 = 0.000636157
TPTTPVm = 0.0568182
TPTTPKs = "9.3583e-05"
TPTTPKr1 = 0.00089213
TPTTPKr2 = "9.8597e-05"
TPTTPKp = "9.6372e-05"
TPTTPKr3 = 0.00054107
TPTTPKr4 = "9.4837e-05"
TPTPGAVm = 0.0568182
TPTPGAKs = 0.00089213
TPTPGAKr1 = "9.3583e-05"
TPTPGAKr2 = "9.8597e-05"
TPTPGAKp = 0.00054107
TPTPGAKr3 = "9.6372e-05"
TPTPGAKr4 = "9.4837e-05"
FBPcAq = 1.00224
FBPcAKs1 = 0.000278407
FBPcAKs2 = 0.000374778
FBPcAKp1 = "2.10226e-05"
FBPcAVm = 0.00568182
FBPcaseVm = 0.00113636
FBPcaseq = 0.792367
FBPcaseKs1 = "2.2129e-05"
FBPcaseKr1 = "1.1065e-06"
FBPcaseKp1 = "6.5319e-05"
FBPcaseKp2 = 0.0018624
UGPaseq = 1.6219
UGPaseKs1 = "3.2124e-05"
UGPaseKs2 = 0.0002364
UGPaseKp1 = 0.00014393
UGPaseKp2 = 0.0013192
UGPaseVm = 0.00410568
SucPcSVm = "7.38636e-05"
SucPcSq = 1.00012
SucPcSKs1 = 0.000278407
SucPcSKr11 = 0.00920241
SucPcSKs2 = 0.000110717
SucPcSKp1 = 0.000642157
SucPcSKp2 = 0.000374778
SucPcSKr12 = 0.00164329
SucPcPVm = 0.0010267
SucPcPq = 1.35286
SucPcPKs1 = "5.354e-05"
SucPcPKp1 = 0.01
SucPcPKp2 = 0.002191
F6PcKVm = "1.02614e-07"
F6PcKKs1 = 0.001
F6PcKKr1 = 0.001
F6PcKKr2 = 0.0015
F26BPcPVm = "2.05284e-10"
F26BPcPKs1 = "1e-09"
F26BPcPKr1 = 0.002
F26BPcPKr3 = 0.001
F26BPcPKr4 = 0.0001
chloroplast = 1.0
cytosol = 1.0
end parameters

begin initial conditions
ADPG[0] = ADPGi
ATP[0] = ATPi
E4P[0] = E4Pi
EOP[0] = EOPi
EP[0] = EPi
EPG[0] = EPGi
EPP[0] = EPPi
ER[0] = ERi
F26BPc[0] = F26BPci
FBP[0] = FBPi
FBPc[0] = FBPci
HeP[0] = HePi
HePc[0] = HePci
PGA[0] = PGAi
PGAc[0] = PGAci
PeP[0] = PePi
RuBP[0] = RuBPi
S7P[0] = S7Pi
SBP[0] = SBPi
SucPc[0] = SucPci
TP[0] = TPi
TPc[0] = TPci
UDPGc[0] = UDPGci
end initial conditions

begin initial values
ADPGi = 6.295e-06
ATPi = 0.00068
E4Pi = 5e-05
EOPi = 0.0002966
EPi = 0.0007045
EPGi = 0.0003206
EPPi = 0.0002711
ERi = 0.001201
F26BPci = 7.8e-09
FBPi = 0.00067
FBPci = 0.002
HePi = 0.0022
HePci = 0.0058
PGAi = 0.0024
PGAci = 0.0
PePi = 0.00025
RuBPi = 0.002
S7Pi = 0.002
SBPi = 0.0003
SucPci = 0.0
TPi = 0.0005
TPci = 0.0023
UDPGci = 0.00057
end initial values

begin assignment rules
GAP := 0.49981684211820576*TP[t]
F6P := 0.33337401159330404*HeP[t]
X5P := 0.33329266258502765*PeP[t]
G6P := 0.3334283604160519*HeP[t]
G1P := 0.33319762799064395*HeP[t]
G1Pc := 0.33319762799064395*HePc[t]
DHAP := 0.5001831578817942*TP[t]
G6Pc := 0.3334283604160519*HePc[t]
R5P := 0.33339701031882757*PeP[t]
ADP := ADT - ATP[t]
GAPc := 0.49981684211820576*TPc[t]
DHAPc := 0.5001831578817942*TPc[t]
E := Et - EOP[t] - EP[t] - EPG[t] - EPP[t] - ER[t]
W4 := -70000*E*Pi + 6*EOP[t]
ADPc := ADTc - ATPc
Ru5P := 0.33331032709614467*PeP[t]
Pic := -ADPc + PiTc - UDPc - 2*(ATPc + PiPic + UTPc + FBPc[t]) - HePc[t] - PGAc[t] - SucPc[t] - TPc[t] - UDPGc[t]
F6Pc := 0.33337401159330404*HePc[t]
Pi := -ADP + PiT - ADPG[t] - E4P[t] - EP[t] - HeP[t] - PeP[t] - PGA[t] - S7P[t] - 2*(PiPi + ATP[t] + EPG[t] + EPP[t] + FBP[t] + RuBP[t] + SBP[t]) - TP[t]
CO2 := (0.000030379746835443035 + 26.666640000000005*O2*ER[t])/(2.531645569620253 + 5280.*ER[t])
end assignment rules

begin function definitions
function27[s1_,s2_,p1_,p2_,q_,Ks1_,Ks2_,Kp1_,Kp2_,Vm_] := ((-((p1*p2)/q) + s1*s2)*Vm)/(Ks1*Ks2*(-1 + (1 + p1/Kp1)*(1 + p2/Kp2) + (1 + s1/Ks1)*(1 + s2/Ks2)))
function23[Vm_,s1_,s2_,p1_,p2_,p3_,q_,Ks1_,Ks2_,Kp1_,Kp2_,Kp3_] := ((-((p1*p2*p3)/q) + s1*s2)*Vm)/(Ks1*Ks2*(p1/Kp1 + p2/Kp2 + p3/Kp3 + (p1*p2*p3)/(Kp1*Kp2*Kp3) + (1 + s1/Ks1)*(1 + s2/Ks2)))
function33[Vm_,s1_,Ks1_,r1_,r2_,Kr1_,r3_,Kr3_,r4_,Kr4_] := ((1 + (r1 + r2)/Kr1)*s1*Vm)/(Ks1*(1 + r3/Kr3 + r4/Kr4))
function25[Vm_,s1_,p1_,p2_,q_,Ks1_,Kp1_,Kp2_] := ((-((p1*p2)/q) + s1)*Vm)/(Ks1*(1 + p1/Kp1 + (p1*p2)/(Kp1*Kp2) + s1/Ks1))
Transketolase[Vm_,q_,s1_,s2_,p1_,p2_,K1_,K2_,K1s2_,K2s1_,r1_,K2r1_,K2s2_,r2_] := ((-(p1*p2) + q*s1*s2)*Vm)/(K1*K2*(1 + (p1 + r2 + p2*(1 + (p1*r2)/K1))/K2 + s2/K2s2 + (r1/K2r1 + s1/K2s1)*(1 + s2/K1s2)))
function24[s1_,s2_,p1_,q_,Ks1_,Ks2_,Kp1_,Vm_] := ((-(p1/q) + s1*s2)*Vm)/(Ks1*Ks2*(p1/Kp1 + (1 + s1/Ks1)*(1 + s2/Ks2)))
function26[Vm_,s1_,p1_,p2_,q_,Ks1_,Kp1_,Kp2_] := ((-((p1*p2)/q) + s1)*Vm)/(Ks1*((1 + p1/Kp1)*(1 + p2/Kp2) + s1/Ks1))
function34[v_] := v
function30[Vm_,s1_,p1_,p2_,q_,Ks1_,r1_,Kr1_,Kp1_,Kp2_] := (s1*(-((p1*p2)/q) + s1)*Vm)/(Ks1^2*(1 + r1/Kr1)^2*((1 + p1/Kp1)*(1 + p2/Kp2) + s1^2/(Ks1^2*(1 + r1/Kr1)^2)))
TPTout[Vm_,s_,Ks_,r1_,Kr1_,r2_,Kr2_,p_,Kp_,r3_,Kr3_,r4_,Kr4_] := ((((p/Kp + r3/Kr3 + r4/Kr4)*s)/Ks - (p*(r1/Kr1 + r2/Kr2 + s/Ks))/Kp)*Vm)/(p/Kp + r1/Kr1 + r2/Kr2 + r3/Kr3 + r4/Kr4 + s/Ks + (p/Kp + r3/Kr3 + r4/Kr4)*(r1/Kr1 + r2/Kr2 + s/Ks))
function29[Vm_,s1_,p1_,q_,Ks1_,Kp1_] := ((-(p1/q) + s1)*Vm)/(Ks1*(1 + p1/Kp1 + s1/Ks1))
function32[Vm_,s1_,Ks1_,r1_,Kr1_,r2_,r3_,Kr2_] := ((1 + r1/Kr1)*s1*Vm)/(Ks1*(1 + (r2 + r3)/Kr2))
function28[Vm_,r1_,r2_,s1_,s2_,p1_,p2_,q_,Ks1_,Ks2_,Kp1_,Kp2_] := (r1^2*(-((p1*p2)/q) + s1*s2)*Vm)/(Ks1*Ks2*r2^2*(-1 + (1 + p1/Kp1)*(1 + p2/Kp2) + (1 + s1/Ks1)*(1 + s2/Ks2)))
function31[Vm_,s1_,s2_,p1_,p2_,p3_,q_,Ks1_,r1_,Kr11_,Ks2_,Kp1_,Kp2_,Kr12_] := (s1*(-((p1*p2*p3)/q) + s1*s2)*Vm)/(Ks1^2*Ks2*(1 + r1/Kr11)^2*(-1 + (1 + p1/Kp1)*(1 + p2/Kp2) + r1/Kr12 + (1 + s1^2/(Ks1^2*(1 + r1/Kr11)^2))*(1 + s2/Ks2)))
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
DHAP http://identifiers.org/chebi/CHEBI:16108
E4P http://identifiers.org/pubchem.compound/122357
F6P http://identifiers.org/chebi/CHEBI:15946
G6P http://identifiers.org/pubchem.compound/439958
Pi http://identifiers.org/chebi/CHEBI:43474
RuBP http://identifiers.org/chebi/CHEBI:16710
end species annotations

begin units
end units