(* Generated by JWS Online *) (* This is an experimental feature of JWS Online. Please report any mistakes.*) (* Note that the following notable SBML entities or features are not supported in notebook outputyet: *) (* Events *) (* Constraints *) (* Units and UnitDefinitions *) (* AlgebraicRules *) (* conversionFactors *) variables = { ADPG[t], ATP[t], E4P[t], EOP[t], EP[t], EPG[t], EPP[t], ER[t], F26BPc[t], FBP[t], FBPc[t], HeP[t], HePc[t], PGA[t], PGAc[t], PeP[t], RuBP[t], S7P[t], SBP[t], SucPc[t], TP[t], TPc[t], UDPGc[t] }; initialValues = { ADPG[0] == 6.295*^-06, ATP[0] == 0.00068, E4P[0] == 5*^-05, EOP[0] == 0.0002966, EP[0] == 0.0007045, EPG[0] == 0.0003206, EPP[0] == 0.0002711, ER[0] == 0.001201, F26BPc[0] == 7.8*^-09, FBP[0] == 0.00067, FBPc[0] == 0.002, HeP[0] == 0.0022, HePc[0] == 0.0058, PGA[0] == 0.0024, PGAc[0] == 0.0, PeP[0] == 0.00025, RuBP[0] == 0.002, S7P[0] == 0.002, SBP[0] == 0.0003, SucPc[0] == 0.0, TP[0] == 0.0005, TPc[0] == 0.0023, UDPGc[0] == 0.00057 }; rates = { AGPase, ATP\[LetterSpace]S, F26BPc\[LetterSpace]P, F6P\[LetterSpace]TK, F6Pc\[LetterSpace]K, FBP\[LetterSpace]A, FBPase, FBPc\[LetterSpace]A, FBPcase, PGA\[LetterSpace]K, Ru5P\[LetterSpace]K, RuBisCO\[LetterSpace]1, RuBisCO\[LetterSpace]2\[LetterSpace]CO2, RuBisCO\[LetterSpace]4, RuBisCO\[LetterSpace]5\[LetterSpace]EOP, RuBisCO\[LetterSpace]5\[LetterSpace]EP, RuBisCO\[LetterSpace]6\[LetterSpace]O2, RuBisCO\[LetterSpace]7, S7P\[LetterSpace]TK, SBP\[LetterSpace]A, SBPase, StS, SucPc\[LetterSpace]P, SucPc\[LetterSpace]S, TPT\[LetterSpace]PGA, TPT\[LetterSpace]TP, UGPase }; rateEquations = { AGPase -> chloroplast*function\[LetterSpace]28[AGPase\[LetterSpace]Vm, PGA[t], Pi, G1P, ATP[t], ADPG[t], PiPi, AGPase\[LetterSpace]q, AGPase\[LetterSpace]Ks1, AGPase\[LetterSpace]Ks2, AGPase\[LetterSpace]Kp1, AGPase\[LetterSpace]Kp2], ATP\[LetterSpace]S -> chloroplast*function\[LetterSpace]24[ADP, Pi, ATP[t], q12, ATP\[LetterSpace]S\[LetterSpace]Ks1, ATP\[LetterSpace]S\[LetterSpace]Ks2, Kp12, ATP\[LetterSpace]S\[LetterSpace]Vm], F26BPc\[LetterSpace]P -> cytosol*function\[LetterSpace]33[F26BPc\[LetterSpace]P\[LetterSpace]Vm, F26BPc[t], F26BPc\[LetterSpace]P\[LetterSpace]Ks1, TPc[t], PGAc[t], F26BPc\[LetterSpace]P\[LetterSpace]Kr1, Pic, F26BPc\[LetterSpace]P\[LetterSpace]Kr3, HePc[t], F26BPc\[LetterSpace]P\[LetterSpace]Kr4], F6P\[LetterSpace]TK -> chloroplast*Transketolase[F6P\[LetterSpace]TK\[LetterSpace]Vm, F6P\[LetterSpace]TK\[LetterSpace]q, F6P, GAP, E4P[t], X5P, F6P\[LetterSpace]TK\[LetterSpace]K1, F6P\[LetterSpace]TK\[LetterSpace]K2, F6P\[LetterSpace]TK\[LetterSpace]K1s2, F6P\[LetterSpace]TK\[LetterSpace]K2s1, S7P[t], F6P\[LetterSpace]TK\[LetterSpace]K2r1, F6P\[LetterSpace]TK\[LetterSpace]K2s2, R5P], F6Pc\[LetterSpace]K -> cytosol*function\[LetterSpace]32[F6Pc\[LetterSpace]K\[LetterSpace]Vm, F6Pc, F6Pc\[LetterSpace]K\[LetterSpace]Ks1, Pic, F6Pc\[LetterSpace]K\[LetterSpace]Kr1, TPc[t], PGAc[t], F6Pc\[LetterSpace]K\[LetterSpace]Kr2], FBP\[LetterSpace]A -> chloroplast*function\[LetterSpace]24[GAP, DHAP, FBP[t], FBP\[LetterSpace]A\[LetterSpace]q, FBP\[LetterSpace]A\[LetterSpace]Ks1, FBP\[LetterSpace]A\[LetterSpace]Ks2, FBP\[LetterSpace]A\[LetterSpace]Kp1, FBP\[LetterSpace]A\[LetterSpace]Vm], FBPase -> chloroplast*function\[LetterSpace]25[FBPase\[LetterSpace]Vm, FBP[t], F6P, Pi, FBPase\[LetterSpace]q, FBPase\[LetterSpace]Ks1, FBPase\[LetterSpace]Kp1, FBPase\[LetterSpace]Kp2], FBPc\[LetterSpace]A -> cytosol*function\[LetterSpace]24[GAPc, DHAPc, FBPc[t], FBPc\[LetterSpace]A\[LetterSpace]q, FBPc\[LetterSpace]A\[LetterSpace]Ks1, FBPc\[LetterSpace]A\[LetterSpace]Ks2, FBPc\[LetterSpace]A\[LetterSpace]Kp1, FBPc\[LetterSpace]A\[LetterSpace]Vm], FBPcase -> cytosol*function\[LetterSpace]30[FBPcase\[LetterSpace]Vm, FBPc[t], F6Pc, Pic, FBPcase\[LetterSpace]q, FBPcase\[LetterSpace]Ks1, F26BPc[t], FBPcase\[LetterSpace]Kr1, FBPcase\[LetterSpace]Kp1, FBPcase\[LetterSpace]Kp2], PGA\[LetterSpace]K -> chloroplast*function\[LetterSpace]23[PGA\[LetterSpace]K\[LetterSpace]Vm, PGA[t], ATP[t], GAP, ADP, Pi, q1, PGA\[LetterSpace]K\[LetterSpace]Ks1, PGA\[LetterSpace]K\[LetterSpace]Ks2, PGA\[LetterSpace]K\[LetterSpace]Kp1, PGA\[LetterSpace]K\[LetterSpace]Kp2, PGA\[LetterSpace]K\[LetterSpace]Kp3], Ru5P\[LetterSpace]K -> chloroplast*function\[LetterSpace]27[Ru5P, ATP[t], RuBP[t], ADP, Ru5P\[LetterSpace]K\[LetterSpace]q, Ru5P\[LetterSpace]K\[LetterSpace]Ks1, Ru5P\[LetterSpace]K\[LetterSpace]Ks2, Ru5P\[LetterSpace]K\[LetterSpace]Kp1, Ru5P\[LetterSpace]K\[LetterSpace]Kp2, Ru5P\[LetterSpace]K\[LetterSpace]Vm], RuBisCO\[LetterSpace]1 -> chloroplast*(-(RuBisCO\[LetterSpace]1\[LetterSpace]k2*ER[t]) + E*RuBisCO\[LetterSpace]1\[LetterSpace]k1*RuBP[t]), RuBisCO\[LetterSpace]2\[LetterSpace]CO2 -> chloroplast*CO2*RuBisCO\[LetterSpace]2\[LetterSpace]CO2\[LetterSpace]k1*ER[t], RuBisCO\[LetterSpace]4 -> chloroplast*(RuBisCO\[LetterSpace]4\[LetterSpace]k1*EPP[t] - RuBisCO\[LetterSpace]4\[LetterSpace]k2*EP[t]*PGA[t]), RuBisCO\[LetterSpace]5\[LetterSpace]EOP -> chloroplast*function\[LetterSpace]34[W4], RuBisCO\[LetterSpace]5\[LetterSpace]EP -> chloroplast*(RuBisCO\[LetterSpace]5\[LetterSpace]EP\[LetterSpace]k1*EP[t] - E*RuBisCO\[LetterSpace]5\[LetterSpace]EP\[LetterSpace]k2*PGA[t]), RuBisCO\[LetterSpace]6\[LetterSpace]O2 -> chloroplast*O2*RuBisCO\[LetterSpace]6\[LetterSpace]O2\[LetterSpace]k1*ER[t], RuBisCO\[LetterSpace]7 -> chloroplast*RuBisCO\[LetterSpace]7\[LetterSpace]k1*EPG[t], S7P\[LetterSpace]TK -> chloroplast*Transketolase[S7P\[LetterSpace]TK\[LetterSpace]Vm, S7P\[LetterSpace]TK\[LetterSpace]q, S7P[t], GAP, R5P, X5P, S7P\[LetterSpace]TK\[LetterSpace]K1, S7P\[LetterSpace]TK\[LetterSpace]K2, S7P\[LetterSpace]TK\[LetterSpace]K1s2, S7P\[LetterSpace]TK\[LetterSpace]K2s1, F6P, S7P\[LetterSpace]TK\[LetterSpace]K2r1, S7P\[LetterSpace]TK\[LetterSpace]K2s2, E4P[t]], SBP\[LetterSpace]A -> chloroplast*function\[LetterSpace]24[E4P[t], DHAP, SBP[t], SBP\[LetterSpace]A\[LetterSpace]q, SBP\[LetterSpace]A\[LetterSpace]Ks1, SBP\[LetterSpace]A\[LetterSpace]Ks2, SBP\[LetterSpace]A\[LetterSpace]Kp1, SBP\[LetterSpace]A\[LetterSpace]Vm], SBPase -> chloroplast*function\[LetterSpace]26[SBPase\[LetterSpace]Vm, SBP[t], S7P[t], Pi, SBPase\[LetterSpace]q, SBPase\[LetterSpace]Ks1, SBPase\[LetterSpace]Kp1, SBPase\[LetterSpace]Kp2], StS -> chloroplast*function\[LetterSpace]29[StS\[LetterSpace]Vm, ADPG[t], ADP, StS\[LetterSpace]q, StS\[LetterSpace]Ks1, StS\[LetterSpace]Kp1], SucPc\[LetterSpace]P -> cytosol*function\[LetterSpace]25[SucPc\[LetterSpace]P\[LetterSpace]Vm, SucPc[t], Succ, Pic, SucPc\[LetterSpace]P\[LetterSpace]q, SucPc\[LetterSpace]P\[LetterSpace]Ks1, SucPc\[LetterSpace]P\[LetterSpace]Kp1, SucPc\[LetterSpace]P\[LetterSpace]Kp2], SucPc\[LetterSpace]S -> cytosol*function\[LetterSpace]31[SucPc\[LetterSpace]S\[LetterSpace]Vm, F6Pc, UDPGc[t], UDPc, SucPc[t], Hc, SucPc\[LetterSpace]S\[LetterSpace]q, SucPc\[LetterSpace]S\[LetterSpace]Ks1, Pic, SucPc\[LetterSpace]S\[LetterSpace]Kr11, SucPc\[LetterSpace]S\[LetterSpace]Ks2, SucPc\[LetterSpace]S\[LetterSpace]Kp1, SucPc\[LetterSpace]S\[LetterSpace]Kp2, SucPc\[LetterSpace]S\[LetterSpace]Kr12], TPT\[LetterSpace]PGA -> TPTout[TPT\[LetterSpace]PGA\[LetterSpace]Vm, PGA[t], TPT\[LetterSpace]PGA\[LetterSpace]Ks, TP[t], TPT\[LetterSpace]PGA\[LetterSpace]Kr1, Pi, TPT\[LetterSpace]PGA\[LetterSpace]Kr2, PGAc[t], TPT\[LetterSpace]PGA\[LetterSpace]Kp, TPc[t], TPT\[LetterSpace]PGA\[LetterSpace]Kr3, Pic, TPT\[LetterSpace]PGA\[LetterSpace]Kr4], TPT\[LetterSpace]TP -> TPTout[TPT\[LetterSpace]TP\[LetterSpace]Vm, TP[t], TPT\[LetterSpace]TP\[LetterSpace]Ks, PGA[t], TPT\[LetterSpace]TP\[LetterSpace]Kr1, Pi, TPT\[LetterSpace]TP\[LetterSpace]Kr2, TPc[t], TPT\[LetterSpace]TP\[LetterSpace]Kp, PGAc[t], TPT\[LetterSpace]TP\[LetterSpace]Kr3, Pic, TPT\[LetterSpace]TP\[LetterSpace]Kr4], UGPase -> cytosol*function\[LetterSpace]27[G1Pc, UTPc, UDPGc[t], PiPic, UGPase\[LetterSpace]q, UGPase\[LetterSpace]Ks1, UGPase\[LetterSpace]Ks2, UGPase\[LetterSpace]Kp1, UGPase\[LetterSpace]Kp2, UGPase\[LetterSpace]Vm] }; 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.386364*^-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 -> 1*^-06, PiPic -> 4*^-05, Succ -> 0.0, UDPc -> 0.00064, UTPc -> 0.00036, FBPcase\[LetterSpace]Vm -> 0.00113636, RuBisCO\[LetterSpace]1\[LetterSpace]k1 -> 50000.0, RuBisCO\[LetterSpace]1\[LetterSpace]k2 -> 0.9, RuBisCO\[LetterSpace]2\[LetterSpace]CO2\[LetterSpace]k1 -> 300000.0, RuBisCO\[LetterSpace]4\[LetterSpace]k1 -> 6.0, RuBisCO\[LetterSpace]4\[LetterSpace]k2 -> 0.0, RuBisCO\[LetterSpace]5\[LetterSpace]EP\[LetterSpace]k1 -> 6.0, RuBisCO\[LetterSpace]5\[LetterSpace]EP\[LetterSpace]k2 -> 70000.0, RuBisCO\[LetterSpace]6\[LetterSpace]O2\[LetterSpace]k1 -> 3030.3, RuBisCO\[LetterSpace]7\[LetterSpace]k1 -> 3.0, PGA\[LetterSpace]K\[LetterSpace]Vm -> 0.0170455, PGA\[LetterSpace]K\[LetterSpace]Ks1 -> 0.0011122, PGA\[LetterSpace]K\[LetterSpace]Ks2 -> 0.0003307, PGA\[LetterSpace]K\[LetterSpace]Kp1 -> 0.00027035, PGA\[LetterSpace]K\[LetterSpace]Kp2 -> 0.00053013, PGA\[LetterSpace]K\[LetterSpace]Kp3 -> 0.0027397, FBP\[LetterSpace]A\[LetterSpace]q -> 1.18815, FBP\[LetterSpace]A\[LetterSpace]Ks1 -> 0.00027035, FBP\[LetterSpace]A\[LetterSpace]Ks2 -> 0.00036393, FBP\[LetterSpace]A\[LetterSpace]Kp1 -> 2.0129*^-05, FBP\[LetterSpace]A\[LetterSpace]Vm -> 0.022727, FBPase\[LetterSpace]Vm -> 0.011364, FBPase\[LetterSpace]q -> 0.77294, FBPase\[LetterSpace]Ks1 -> 3.2842*^-05, FBPase\[LetterSpace]Kp1 -> 6.3429*^-05, FBPase\[LetterSpace]Kp2 -> 0.0017914, F6P\[LetterSpace]TK\[LetterSpace]Vm -> 0.170455, F6P\[LetterSpace]TK\[LetterSpace]q -> 0.99943, F6P\[LetterSpace]TK\[LetterSpace]K1 -> 0.00061349, F6P\[LetterSpace]TK\[LetterSpace]K2 -> 0.00011438, F6P\[LetterSpace]TK\[LetterSpace]K1s2 -> 0.00027035, F6P\[LetterSpace]TK\[LetterSpace]K2s1 -> 0.0005407, F6P\[LetterSpace]TK\[LetterSpace]K2r1 -> 0.00017677, F6P\[LetterSpace]TK\[LetterSpace]K2s2 -> 9.0464*^-05, SBP\[LetterSpace]A\[LetterSpace]q -> 1.18815, SBP\[LetterSpace]A\[LetterSpace]Ks1 -> 0.00017677, SBP\[LetterSpace]A\[LetterSpace]Ks2 -> 0.00036393, SBP\[LetterSpace]A\[LetterSpace]Kp1 -> 2.0129*^-05, SBP\[LetterSpace]A\[LetterSpace]Vm -> 0.011364, SBPase\[LetterSpace]Vm -> 0.00568182, SBPase\[LetterSpace]q -> 0.77294, SBPase\[LetterSpace]Ks1 -> 1.2713*^-05, SBPase\[LetterSpace]Kp1 -> 1.5597*^-05, SBPase\[LetterSpace]Kp2 -> 0.006744, S7P\[LetterSpace]TK\[LetterSpace]Vm -> 0.0821023, S7P\[LetterSpace]TK\[LetterSpace]q -> 0.99996, S7P\[LetterSpace]TK\[LetterSpace]K1 -> 0.00061349, S7P\[LetterSpace]TK\[LetterSpace]K2 -> 0.00011438, S7P\[LetterSpace]TK\[LetterSpace]K1s2 -> 0.00027035, S7P\[LetterSpace]TK\[LetterSpace]K2s1 -> 0.00017677, S7P\[LetterSpace]TK\[LetterSpace]K2r1 -> 0.0005407, S7P\[LetterSpace]TK\[LetterSpace]K2s2 -> 9.0464*^-05, Ru5P\[LetterSpace]K\[LetterSpace]q -> 1.05289, Ru5P\[LetterSpace]K\[LetterSpace]Ks1 -> 3.63934*^-05, Ru5P\[LetterSpace]K\[LetterSpace]Ks2 -> 0.00055117, Ru5P\[LetterSpace]K\[LetterSpace]Kp1 -> 9.95868*^-05, Ru5P\[LetterSpace]K\[LetterSpace]Kp2 -> 9.11825*^-05, Ru5P\[LetterSpace]K\[LetterSpace]Vm -> 0.568182, ATP\[LetterSpace]S\[LetterSpace]Ks1 -> 0.00031808, ATP\[LetterSpace]S\[LetterSpace]Ks2 -> 0.00031612, ATP\[LetterSpace]S\[LetterSpace]Vm -> 0.0284091, AGPase\[LetterSpace]Vm -> 0.00113636, AGPase\[LetterSpace]q -> 0.11059, AGPase\[LetterSpace]Ks1 -> 0.0010398, AGPase\[LetterSpace]Ks2 -> 0.00011023, AGPase\[LetterSpace]Kp1 -> 0.00053013, AGPase\[LetterSpace]Kp2 -> 0.01951, StS\[LetterSpace]Vm -> 0.00284091, StS\[LetterSpace]q -> 1.00326, StS\[LetterSpace]Ks1 -> 0.000212052, StS\[LetterSpace]Kp1 -> 0.000636157, TPT\[LetterSpace]TP\[LetterSpace]Vm -> 0.0568182, TPT\[LetterSpace]TP\[LetterSpace]Ks -> 9.3583*^-05, TPT\[LetterSpace]TP\[LetterSpace]Kr1 -> 0.00089213, TPT\[LetterSpace]TP\[LetterSpace]Kr2 -> 9.8597*^-05, TPT\[LetterSpace]TP\[LetterSpace]Kp -> 9.6372*^-05, TPT\[LetterSpace]TP\[LetterSpace]Kr3 -> 0.00054107, TPT\[LetterSpace]TP\[LetterSpace]Kr4 -> 9.4837*^-05, TPT\[LetterSpace]PGA\[LetterSpace]Vm -> 0.0568182, TPT\[LetterSpace]PGA\[LetterSpace]Ks -> 0.00089213, TPT\[LetterSpace]PGA\[LetterSpace]Kr1 -> 9.3583*^-05, TPT\[LetterSpace]PGA\[LetterSpace]Kr2 -> 9.8597*^-05, TPT\[LetterSpace]PGA\[LetterSpace]Kp -> 0.00054107, TPT\[LetterSpace]PGA\[LetterSpace]Kr3 -> 9.6372*^-05, TPT\[LetterSpace]PGA\[LetterSpace]Kr4 -> 9.4837*^-05, FBPc\[LetterSpace]A\[LetterSpace]q -> 1.00224, FBPc\[LetterSpace]A\[LetterSpace]Ks1 -> 0.000278407, FBPc\[LetterSpace]A\[LetterSpace]Ks2 -> 0.000374778, FBPc\[LetterSpace]A\[LetterSpace]Kp1 -> 2.10226*^-05, FBPc\[LetterSpace]A\[LetterSpace]Vm -> 0.00568182, FBPcase\[LetterSpace]q -> 0.792367, FBPcase\[LetterSpace]Ks1 -> 2.2129*^-05, FBPcase\[LetterSpace]Kr1 -> 1.1065*^-06, FBPcase\[LetterSpace]Kp1 -> 6.5319*^-05, FBPcase\[LetterSpace]Kp2 -> 0.0018624, UGPase\[LetterSpace]q -> 1.6219, UGPase\[LetterSpace]Ks1 -> 3.2124*^-05, UGPase\[LetterSpace]Ks2 -> 0.0002364, UGPase\[LetterSpace]Kp1 -> 0.00014393, UGPase\[LetterSpace]Kp2 -> 0.0013192, UGPase\[LetterSpace]Vm -> 0.00410568, SucPc\[LetterSpace]S\[LetterSpace]Vm -> 7.38636*^-05, SucPc\[LetterSpace]S\[LetterSpace]q -> 1.00012, SucPc\[LetterSpace]S\[LetterSpace]Ks1 -> 0.000278407, SucPc\[LetterSpace]S\[LetterSpace]Kr11 -> 0.00920241, SucPc\[LetterSpace]S\[LetterSpace]Ks2 -> 0.000110717, SucPc\[LetterSpace]S\[LetterSpace]Kp1 -> 0.000642157, SucPc\[LetterSpace]S\[LetterSpace]Kp2 -> 0.000374778, SucPc\[LetterSpace]S\[LetterSpace]Kr12 -> 0.00164329, SucPc\[LetterSpace]P\[LetterSpace]Vm -> 0.0010267, SucPc\[LetterSpace]P\[LetterSpace]q -> 1.35286, SucPc\[LetterSpace]P\[LetterSpace]Ks1 -> 5.354*^-05, SucPc\[LetterSpace]P\[LetterSpace]Kp1 -> 0.01, SucPc\[LetterSpace]P\[LetterSpace]Kp2 -> 0.002191, F6Pc\[LetterSpace]K\[LetterSpace]Vm -> 1.02614*^-07, F6Pc\[LetterSpace]K\[LetterSpace]Ks1 -> 0.001, F6Pc\[LetterSpace]K\[LetterSpace]Kr1 -> 0.001, F6Pc\[LetterSpace]K\[LetterSpace]Kr2 -> 0.0015, F26BPc\[LetterSpace]P\[LetterSpace]Vm -> 2.05284*^-10, F26BPc\[LetterSpace]P\[LetterSpace]Ks1 -> 1*^-09, F26BPc\[LetterSpace]P\[LetterSpace]Kr1 -> 0.002, F26BPc\[LetterSpace]P\[LetterSpace]Kr3 -> 0.001, F26BPc\[LetterSpace]P\[LetterSpace]Kr4 -> 0.0001, chloroplast -> 1.0, cytosol -> 1.0 }; assignments = { function\[LetterSpace]28[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))), function\[LetterSpace]25[Vm_,s1_,p1_,p2_,q_,Ks1_,Kp1_,Kp2_] -> ((-((p1*p2)/q) + s1)*Vm)/(Ks1*(1 + p1/Kp1 + (p1*p2)/(Kp1*Kp2) + s1/Ks1)), function\[LetterSpace]31[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))), function\[LetterSpace]33[Vm_,s1_,Ks1_,r1_,r2_,Kr1_,r3_,Kr3_,r4_,Kr4_] -> ((1 + (r1 + r2)/Kr1)*s1*Vm)/(Ks1*(1 + r3/Kr3 + r4/Kr4)), function\[LetterSpace]26[Vm_,s1_,p1_,p2_,q_,Ks1_,Kp1_,Kp2_] -> ((-((p1*p2)/q) + s1)*Vm)/(Ks1*((1 + p1/Kp1)*(1 + p2/Kp2) + s1/Ks1)), function\[LetterSpace]24[s1_,s2_,p1_,q_,Ks1_,Ks2_,Kp1_,Vm_] -> ((-(p1/q) + s1*s2)*Vm)/(Ks1*Ks2*(p1/Kp1 + (1 + s1/Ks1)*(1 + s2/Ks2))), 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))), function\[LetterSpace]23[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))), function\[LetterSpace]27[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))), function\[LetterSpace]32[Vm_,s1_,Ks1_,r1_,Kr1_,r2_,r3_,Kr2_] -> ((1 + r1/Kr1)*s1*Vm)/(Ks1*(1 + (r2 + r3)/Kr2)), function\[LetterSpace]29[Vm_,s1_,p1_,q_,Ks1_,Kp1_] -> ((-(p1/q) + s1)*Vm)/(Ks1*(1 + p1/Kp1 + s1/Ks1)), 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)), function\[LetterSpace]30[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))), function\[LetterSpace]34[v_] -> v, Pic -> -ADPc + PiTc - UDPc - 2*(ATPc + PiPic + UTPc + FBPc[t]) - HePc[t] - PGAc[t] - SucPc[t] - TPc[t] - UDPGc[t], F6P -> 0.33337401159330404*HeP[t], Ru5P -> 0.33331032709614467*PeP[t], GAP -> 0.49981684211820576*TP[t], ADPc -> ADTc - ATPc, W4 -> -70000*E*Pi + 6*EOP[t], E -> Et - EOP[t] - EP[t] - EPG[t] - EPP[t] - ER[t], CO2 -> (0.000030379746835443035 + 26.666640000000005*O2*ER[t])/(2.531645569620253 + 5280.*ER[t]), DHAPc -> 0.5001831578817942*TPc[t], GAPc -> 0.49981684211820576*TPc[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], ADP -> ADT - ATP[t], R5P -> 0.33339701031882757*PeP[t], F6Pc -> 0.33337401159330404*HePc[t], G6Pc -> 0.3334283604160519*HePc[t], DHAP -> 0.5001831578817942*TP[t], G1Pc -> 0.33319762799064395*HePc[t], G1P -> 0.33319762799064395*HeP[t], G6P -> 0.3334283604160519*HeP[t], X5P -> 0.33329266258502765*PeP[t] }; events = { }; speciesAnnotations = { DHAP[t]->"http://identifiers.org/chebi/CHEBI:16108", E4P[t]->"http://identifiers.org/pubchem.compound/122357", F6P[t]->"http://identifiers.org/chebi/CHEBI:15946", G6P[t]->"http://identifiers.org/pubchem.compound/439958", Pi[t]->"http://identifiers.org/chebi/CHEBI:43474", RuBP[t]->"http://identifiers.org/chebi/CHEBI:16710" }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { ADPG'[t] == 1.0*AGPase -1.0*StS, ATP'[t] == 1.0*ATP\[LetterSpace]S -0.5*RuBisCO\[LetterSpace]6\[LetterSpace]O2 -1.0*PGA\[LetterSpace]K -1.0*Ru5P\[LetterSpace]K -1.0*AGPase, E4P'[t] == 1.0*F6P\[LetterSpace]TK -1.0*SBP\[LetterSpace]A, EOP'[t] == -1.0*RuBisCO\[LetterSpace]5\[LetterSpace]EOP, EP'[t] == 1.0*RuBisCO\[LetterSpace]4 +1.0*RuBisCO\[LetterSpace]7 -1.0*RuBisCO\[LetterSpace]5\[LetterSpace]EP, EPG'[t] == 1.0*RuBisCO\[LetterSpace]6\[LetterSpace]O2 -1.0*RuBisCO\[LetterSpace]7, EPP'[t] == 1.0*RuBisCO\[LetterSpace]2\[LetterSpace]CO2 -1.0*RuBisCO\[LetterSpace]4, ER'[t] == 1.0*RuBisCO\[LetterSpace]1 -1.0*RuBisCO\[LetterSpace]2\[LetterSpace]CO2 -1.0*RuBisCO\[LetterSpace]6\[LetterSpace]O2, F26BPc'[t] == 1.0*F6Pc\[LetterSpace]K -1.0*F26BPc\[LetterSpace]P, FBP'[t] == 1.0*FBP\[LetterSpace]A -1.0*FBPase, FBPc'[t] == 1.0*FBPc\[LetterSpace]A -1.0*FBPcase, HeP'[t] == 1.0*FBPase -1.0*F6P\[LetterSpace]TK -1.0*AGPase, HePc'[t] == 1.0*FBPcase +1.0*F26BPc\[LetterSpace]P -1.0*UGPase -1.0*SucPc\[LetterSpace]S -1.0*F6Pc\[LetterSpace]K, PGA'[t] == 1.0*RuBisCO\[LetterSpace]4 +1.0*RuBisCO\[LetterSpace]5\[LetterSpace]EP +0.5*RuBisCO\[LetterSpace]6\[LetterSpace]O2 -1.0*PGA\[LetterSpace]K -1.0*TPT\[LetterSpace]PGA, PGAc'[t] == 1.0*TPT\[LetterSpace]PGA , PeP'[t] == 1.0*F6P\[LetterSpace]TK +2.0*S7P\[LetterSpace]TK -1.0*Ru5P\[LetterSpace]K, RuBP'[t] == 1.0*Ru5P\[LetterSpace]K -1.0*RuBisCO\[LetterSpace]1, S7P'[t] == 1.0*SBPase -1.0*S7P\[LetterSpace]TK, SBP'[t] == 1.0*SBP\[LetterSpace]A -1.0*SBPase, SucPc'[t] == 1.0*SucPc\[LetterSpace]S -1.0*SucPc\[LetterSpace]P, TP'[t] == 1.0*PGA\[LetterSpace]K -2.0*FBP\[LetterSpace]A -1.0*F6P\[LetterSpace]TK -1.0*SBP\[LetterSpace]A -1.0*S7P\[LetterSpace]TK -1.0*TPT\[LetterSpace]TP, TPc'[t] == 1.0*TPT\[LetterSpace]TP -2.0*FBPc\[LetterSpace]A, UDPGc'[t] == 1.0*UGPase -1.0*SucPc\[LetterSpace]S }; timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}]; (* Steady-state solution initialized with result of time evolution *) findRootEquations = odes /.D[_[t],t]->0; findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2]; steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100] fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations (* Plot the time evolution of the variables *) plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]