(* 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 = { ATP[t], DPGA[t], E4P[t], F26BPc[t], FBP[t], FBPc[t], GCA[t], GCAc[t], GCEA[t], GCEAc[t], GLYc[t], GOAc[t], HPRc[t], HeP[t], HePc[t], PGA[t], PGAc[t], PGCA[t], PeP[t], RuBP[t], S7P[t], SBP[t], SERc[t], SucPc[t], Succ[t], TP[t], TPc[t], UDPGc[t] }; initialValues = { ATP[0] == 0.68, DPGA[0] == 0.0011, E4P[0] == 0.05, F26BPc[0] == 7.8*^-06, FBP[0] == 0.669999999999999, FBPc[0] == 2.0, GCA[0] == 0.36, GCAc[0] == 0.36, GCEA[0] == 0.1812, GCEAc[0] == 0.1812, GLYc[0] == 1.8, GOAc[0] == 0.028, HPRc[0] == 0.0035, HeP[0] == 2.2, HePc[0] == 5.8, PGA[0] == 2.4, PGAc[0] == 0.0, PGCA[0] == 0.0029, PeP[0] == 0.25, RuBP[0] == 2.0, S7P[0] == 2.0, SBP[0] == 0.3, SERc[0] == 7.5, SucPc[0] == 0.0, Succ[0] == 0.0, TP[0] == 0.5, TPc[0] == 2.3, UDPGc[0] == 0.57 }; rates = { AGPase, ATP\[LetterSpace]S, F26BPc\[LetterSpace]P, F6P\[LetterSpace]TK, F6Pc\[LetterSpace]K, FBP\[LetterSpace]A, FBPase, FBPc\[LetterSpace]A, FBPcase, GAP\[LetterSpace]DH, GCAc\[LetterSpace]Ox, GCEA\[LetterSpace]DH, GCEA\[LetterSpace]K, GLUcGOAc\[LetterSpace]AT, GLYc\[LetterSpace]DC, GPT\[LetterSpace]GCA, GPT\[LetterSpace]GCAc, GPT\[LetterSpace]GCEA, GPT\[LetterSpace]GCEAc, PGA\[LetterSpace]K, PGAc\[LetterSpace]Deg, PGCA\[LetterSpace]P, Ru5P\[LetterSpace]K, RuBisCO\[LetterSpace]CO2, RuBisCO\[LetterSpace]O2, S7P\[LetterSpace]TK, SBP\[LetterSpace]A, SBPase, SERcGOAc\[LetterSpace]AT, SucPc\[LetterSpace]P, SucPc\[LetterSpace]S, Succ\[LetterSpace]Deg, TPT\[LetterSpace]DHAP, TPT\[LetterSpace]GAP, TPT\[LetterSpace]PGA, UGPase }; rateEquations = { AGPase -> chloroplast*function\[LetterSpace]10[AGPase\[LetterSpace]Vm, G1P, ATP[t], AGPase\[LetterSpace]K1, AGPase\[LetterSpace]K2, ADP, AGPase\[LetterSpace]KR1, Pi, AGPase\[LetterSpace]KA1, PGA[t], AGPase\[LetterSpace]KA2, F6P, AGPase\[LetterSpace]KA3, FBP[t]], ATP\[LetterSpace]S -> chloroplast*function\[LetterSpace]4[ADP, Pi, ATP[t], ATP\[LetterSpace]S\[LetterSpace]q, ATP\[LetterSpace]S\[LetterSpace]Ks1, ATP\[LetterSpace]S\[LetterSpace]Ks2, ATP\[LetterSpace]S\[LetterSpace]Kp1, ATP\[LetterSpace]S\[LetterSpace]Vm], F26BPc\[LetterSpace]P -> cytosol*function\[LetterSpace]16[F26BPc\[LetterSpace]P\[LetterSpace]Vm, F26BPc[t], Pic, F26BPc\[LetterSpace]P\[LetterSpace]Kr1, F6Pc, F26BPc\[LetterSpace]P\[LetterSpace]Kr2, F26BPc\[LetterSpace]P\[LetterSpace]Ks1], F6P\[LetterSpace]TK -> chloroplast*function\[LetterSpace]6[F6P\[LetterSpace]TK\[LetterSpace]Vm, F6P, GAP, X5P, E4P[t], F6P\[LetterSpace]TK\[LetterSpace]Ks1, F6P\[LetterSpace]TK\[LetterSpace]Ks2, X5P, F6P\[LetterSpace]TK\[LetterSpace]Kr1, E4P[t], F6P\[LetterSpace]TK\[LetterSpace]Kr2, F6P\[LetterSpace]TK\[LetterSpace]q], F6Pc\[LetterSpace]K -> cytosol*function\[LetterSpace]15[F6Pc\[LetterSpace]K\[LetterSpace]Vm, F6Pc, ATPc, F26BPc[t], ADPc, F6Pc\[LetterSpace]K\[LetterSpace]q, F6Pc\[LetterSpace]K\[LetterSpace]Ks1, F26BPc[t], F6Pc\[LetterSpace]K\[LetterSpace]Kr1, DHAPc, F6Pc\[LetterSpace]K\[LetterSpace]Kr2, F6Pc\[LetterSpace]K\[LetterSpace]Ks2, ADPc, F6Pc\[LetterSpace]K\[LetterSpace]Kr3], FBP\[LetterSpace]A -> chloroplast*function\[LetterSpace]4[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]5[FBPase\[LetterSpace]Vm, FBP[t], F6P, Pi, FBPase\[LetterSpace]q, FBPase\[LetterSpace]Ks1, F6P, FBPase\[LetterSpace]Kr1, Pi, FBPase\[LetterSpace]Kr2], FBPc\[LetterSpace]A -> cytosol*function\[LetterSpace]4[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]12[FBPc[t], F6Pc, Pic, FBPcase\[LetterSpace]q, K52a, FBPcase\[LetterSpace]Kp1, FBPcase\[LetterSpace]Kp2, FBPcase\[LetterSpace]Vm], GAP\[LetterSpace]DH -> chloroplast*function\[LetterSpace]3[GAP\[LetterSpace]DH\[LetterSpace]Vm, DPGA[t], NADPH, GAP\[LetterSpace]DH\[LetterSpace]K1, GAP\[LetterSpace]DH\[LetterSpace]K2], GCAc\[LetterSpace]Ox -> cytosol*function\[LetterSpace]17[GCAc[t], GCAc\[LetterSpace]Ox\[LetterSpace]Km, GCAc\[LetterSpace]Ox\[LetterSpace]V], GCEA\[LetterSpace]DH -> function\[LetterSpace]21[GCEA\[LetterSpace]DH\[LetterSpace]Vm, HPRc[t], NADH, GCEAc[t], NAD, GCEA\[LetterSpace]DH\[LetterSpace]q, GCEA\[LetterSpace]DH\[LetterSpace]Ks1, HPRc[t], GCEA\[LetterSpace]DH\[LetterSpace]Kr1], GCEA\[LetterSpace]K -> chloroplast*function\[LetterSpace]20[GCEA\[LetterSpace]K\[LetterSpace]Vm, ATP[t], GCEA[t], PGA[t], ADP, GCEA\[LetterSpace]K\[LetterSpace]q, GCEA\[LetterSpace]K\[LetterSpace]Ks1, PGA[t], GCEA\[LetterSpace]K\[LetterSpace]Kr1, GCEA\[LetterSpace]K\[LetterSpace]Ks2], GLUcGOAc\[LetterSpace]AT -> cytosol*function\[LetterSpace]20[GLUcGOAc\[LetterSpace]AT\[LetterSpace]Vm, GLUc, GOAc[t], KGc, GLYc[t], GLUcGOAc\[LetterSpace]AT\[LetterSpace]q, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Ks1, GLYc[t], GLUcGOAc\[LetterSpace]AT\[LetterSpace]Kr1, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Ks2], GLYc\[LetterSpace]DC -> cytosol*function\[LetterSpace]22[GLYc\[LetterSpace]DC\[LetterSpace]Vm, GLYc[t], GLYc\[LetterSpace]DC\[LetterSpace]Ks1, SERc[t], GLYc\[LetterSpace]DC\[LetterSpace]Kr1], GPT\[LetterSpace]GCA -> function\[LetterSpace]18[GCA[t], GCEA[t], GPT\[LetterSpace]GCA\[LetterSpace]Km, GPT\[LetterSpace]GCA\[LetterSpace]V, GPT\[LetterSpace]GCA\[LetterSpace]Ki], GPT\[LetterSpace]GCAc -> function\[LetterSpace]18[GCAc[t], GCEAc[t], GPT\[LetterSpace]GCAc\[LetterSpace]Km, GPT\[LetterSpace]GCAc\[LetterSpace]V, GPT\[LetterSpace]GCAc\[LetterSpace]Ki], GPT\[LetterSpace]GCEA -> function\[LetterSpace]18[GCEA[t], GCA[t], GPT\[LetterSpace]GCEA\[LetterSpace]Km, GPT\[LetterSpace]GCEA\[LetterSpace]V, GPT\[LetterSpace]GCEA\[LetterSpace]Ki], GPT\[LetterSpace]GCEAc -> function\[LetterSpace]18[GCEAc[t], GCAc[t], GPT\[LetterSpace]GCEAc\[LetterSpace]Km, GPT\[LetterSpace]GCEAc\[LetterSpace]V, GPT\[LetterSpace]GCEAc\[LetterSpace]Ki], PGA\[LetterSpace]K -> chloroplast*function\[LetterSpace]2[PGA\[LetterSpace]K\[LetterSpace]Vm, PGA[t], ATP[t], PGA\[LetterSpace]K\[LetterSpace]Ks1, PGA\[LetterSpace]K\[LetterSpace]Ks2, ADP, PGA\[LetterSpace]K\[LetterSpace]Kr1], PGAc\[LetterSpace]Deg -> cytosol*function\[LetterSpace]17[PGAc[t], PGAc\[LetterSpace]Deg\[LetterSpace]Km, PGAc\[LetterSpace]Deg\[LetterSpace]V], PGCA\[LetterSpace]P -> chloroplast*function\[LetterSpace]19[PGCA\[LetterSpace]P\[LetterSpace]Vm, PGCA[t], PGCA\[LetterSpace]P\[LetterSpace]Ks1, GCA[t], PGCA\[LetterSpace]P\[LetterSpace]Kr1, Pi, PGCA\[LetterSpace]P\[LetterSpace]Kr2], Ru5P\[LetterSpace]K -> chloroplast*function\[LetterSpace]9[Ru5P\[LetterSpace]K\[LetterSpace]Vm, Ru5P, ATP[t], RuBP[t], ADP, Ru5P\[LetterSpace]K\[LetterSpace]q, Ru5P\[LetterSpace]K\[LetterSpace]Ks1, PGA[t], Ru5P\[LetterSpace]K\[LetterSpace]Kr1, RuBP[t], Ru5P\[LetterSpace]K\[LetterSpace]Kr2, Pi, Ru5P\[LetterSpace]K\[LetterSpace]Kr3, ADP, Ru5P\[LetterSpace]K\[LetterSpace]Kr41, Ru5P\[LetterSpace]K\[LetterSpace]Ks2, Ru5P\[LetterSpace]K\[LetterSpace]Kr42], RuBisCO\[LetterSpace]CO2 -> chloroplast*function\[LetterSpace]1[Wc\[LetterSpace]min, RuBP[t], RuBisCO\[LetterSpace]CO2\[LetterSpace]K, PGA[t], RuBisCO\[LetterSpace]CO2\[LetterSpace]KR1, FBP[t], RuBisCO\[LetterSpace]CO2\[LetterSpace]KR2, SBP[t], RuBisCO\[LetterSpace]CO2\[LetterSpace]KR3, Pi, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR4, NADPH, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR5], RuBisCO\[LetterSpace]O2 -> chloroplast*function\[LetterSpace]1[Wo\[LetterSpace]min, RuBP[t], RuBisCO\[LetterSpace]O2\[LetterSpace]K, PGA[t], RuBisCO\[LetterSpace]O2\[LetterSpace]KR1, FBP[t], RuBisCO\[LetterSpace]O2\[LetterSpace]KR2, SBP[t], RuBisCO\[LetterSpace]O2\[LetterSpace]KR3, Pi, RuBisCO\[LetterSpace]O2\[LetterSpace]KR4, NADPH, RuBisCO\[LetterSpace]O2\[LetterSpace]KR5], S7P\[LetterSpace]TK -> chloroplast*function\[LetterSpace]6[S7P\[LetterSpace]TK\[LetterSpace]Vm, GAP, S7P[t], X5P, R5P, S7P\[LetterSpace]TK\[LetterSpace]Ks1, S7P\[LetterSpace]TK\[LetterSpace]Ks2, X5P, S7P\[LetterSpace]TK\[LetterSpace]Kr1, R5P, S7P\[LetterSpace]TK\[LetterSpace]Kr2, S7P\[LetterSpace]TK\[LetterSpace]q], SBP\[LetterSpace]A -> chloroplast*function\[LetterSpace]7[SBP\[LetterSpace]A\[LetterSpace]Vm, DHAP, E4P[t], SBP[t], SBP\[LetterSpace]A\[LetterSpace]q, SBP\[LetterSpace]A\[LetterSpace]Ks1, SBP\[LetterSpace]A\[LetterSpace]Ks2], SBPase -> chloroplast*function\[LetterSpace]8[SBPase\[LetterSpace]Vm, SBP[t], S7P[t], Pi, SBPase\[LetterSpace]q, SBPase\[LetterSpace]Ks1, Pi, SBPase\[LetterSpace]Kr1], SERcGOAc\[LetterSpace]AT -> function\[LetterSpace]20[SERcGOAc\[LetterSpace]AT\[LetterSpace]Vm, SERc[t], GOAc[t], HPRc[t], GLYc[t], SERcGOAc\[LetterSpace]AT\[LetterSpace]q, SERcGOAc\[LetterSpace]AT\[LetterSpace]Ks1, GLYc[t], SERcGOAc\[LetterSpace]AT\[LetterSpace]Kr1, SERcGOAc\[LetterSpace]AT\[LetterSpace]Ks2], SucPc\[LetterSpace]P -> cytosol*function\[LetterSpace]8[SucPc\[LetterSpace]P\[LetterSpace]Vm, SucPc[t], Succ[t], Pic, SucPc\[LetterSpace]P\[LetterSpace]q, SucPc\[LetterSpace]P\[LetterSpace]Ks1, Succ[t], SucPc\[LetterSpace]P\[LetterSpace]Kr1], SucPc\[LetterSpace]S -> cytosol*function\[LetterSpace]14[SucPc\[LetterSpace]S\[LetterSpace]Vm, F6Pc, UDPGc[t], SucPc[t], UDPc, SucPc\[LetterSpace]S\[LetterSpace]q, SucPc\[LetterSpace]S\[LetterSpace]Ks1, FBPc[t], SucPc\[LetterSpace]S\[LetterSpace]Kr1, SucPc\[LetterSpace]S\[LetterSpace]Ks2, UDPc, SucPc\[LetterSpace]S\[LetterSpace]Kr2, SucPc[t], SucPc\[LetterSpace]S\[LetterSpace]Kr3, Succ[t], SucPc\[LetterSpace]S\[LetterSpace]Kr4, Pic, SucPc\[LetterSpace]S\[LetterSpace]Kr5], Succ\[LetterSpace]Deg -> cytosol*function\[LetterSpace]17[Succ[t], Succ\[LetterSpace]Deg\[LetterSpace]Km, Succ\[LetterSpace]Deg\[LetterSpace]V], TPT\[LetterSpace]DHAP -> function\[LetterSpace]11[TPT\[LetterSpace]DHAP\[LetterSpace]Vm, DHAP, TPT\[LetterSpace]DHAP\[LetterSpace]KA, Pic, TPT\[LetterSpace]DHAP\[LetterSpace]K, Pi, TPT\[LetterSpace]DHAP\[LetterSpace]KR1, PGA[t], TPT\[LetterSpace]DHAP\[LetterSpace]KR2, GAP, TPT\[LetterSpace]DHAP\[LetterSpace]KR3], TPT\[LetterSpace]GAP -> function\[LetterSpace]11[TPT\[LetterSpace]GAP\[LetterSpace]Vm, GAP, TPT\[LetterSpace]GAP\[LetterSpace]KA, Pic, TPT\[LetterSpace]GAP\[LetterSpace]K, Pi, TPT\[LetterSpace]GAP\[LetterSpace]KR1, PGA[t], TPT\[LetterSpace]GAP\[LetterSpace]KR2, DHAP, TPT\[LetterSpace]GAP\[LetterSpace]KR3], TPT\[LetterSpace]PGA -> function\[LetterSpace]11[TPT\[LetterSpace]PGA\[LetterSpace]Vm, PGA[t], TPT\[LetterSpace]PGA\[LetterSpace]KA, Pic, TPT\[LetterSpace]PGA\[LetterSpace]K, Pi, TPT\[LetterSpace]PGA\[LetterSpace]KR1, GAP, TPT\[LetterSpace]PGA\[LetterSpace]KR2, DHAP, TPT\[LetterSpace]PGA\[LetterSpace]KR3], UGPase -> cytosol*function\[LetterSpace]13[G1Pc, UTPc, UDPGc[t], PiPic, UGPase\[LetterSpace]q, UGPase\[LetterSpace]Ks1, UGPase\[LetterSpace]Ks2, UGPase\[LetterSpace]Kp1, UGPase\[LetterSpace]Kp2, UGPase\[LetterSpace]Vm] }; parameters = { E -> 1.456965457, cA -> 1.5, cAc -> 1.0, cNc -> 0.87, cP -> 15.0, cPc -> 15.0, cUc -> 1.5, ATPc -> 0.35, CO2 -> 0.009, GLUc -> 24.0, KGc -> 0.4, NAD -> 0.4, NADH -> 0.47, NADP -> 0.29, NADPH -> 0.21, O2 -> 0.2646, UTPc -> 0.75, RuBisCO\[LetterSpace]CO2\[LetterSpace]K -> 0.02, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR1 -> 0.84, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR2 -> 0.04, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR3 -> 0.075, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR4 -> 0.9, RuBisCO\[LetterSpace]CO2\[LetterSpace]KR5 -> 0.07, PGA\[LetterSpace]K\[LetterSpace]Vm -> 30.1408, PGA\[LetterSpace]K\[LetterSpace]Ks1 -> 0.24, PGA\[LetterSpace]K\[LetterSpace]Ks2 -> 0.39, PGA\[LetterSpace]K\[LetterSpace]Kr1 -> 0.23, GAP\[LetterSpace]DH\[LetterSpace]Vm -> 4.03948, GAP\[LetterSpace]DH\[LetterSpace]K1 -> 0.004, GAP\[LetterSpace]DH\[LetterSpace]K2 -> 0.1, FBP\[LetterSpace]A\[LetterSpace]q -> 7.1, FBP\[LetterSpace]A\[LetterSpace]Ks1 -> 0.3, FBP\[LetterSpace]A\[LetterSpace]Ks2 -> 0.4, FBP\[LetterSpace]A\[LetterSpace]Kp1 -> 0.02, FBP\[LetterSpace]A\[LetterSpace]Vm -> 1.21889, FBPase\[LetterSpace]Vm -> 0.72626, FBPase\[LetterSpace]q -> 666000.0, FBPase\[LetterSpace]Ks1 -> 0.033, FBPase\[LetterSpace]Kr1 -> 0.7, FBPase\[LetterSpace]Kr2 -> 12.0, F6P\[LetterSpace]TK\[LetterSpace]Vm -> 3.12207, F6P\[LetterSpace]TK\[LetterSpace]Ks1 -> 0.1, F6P\[LetterSpace]TK\[LetterSpace]Ks2 -> 0.1, F6P\[LetterSpace]TK\[LetterSpace]Kr1 -> 0.1, F6P\[LetterSpace]TK\[LetterSpace]Kr2 -> 0.1, F6P\[LetterSpace]TK\[LetterSpace]q -> 10.0, SBP\[LetterSpace]A\[LetterSpace]Vm -> 1.21889, SBP\[LetterSpace]A\[LetterSpace]q -> 1.017, SBP\[LetterSpace]A\[LetterSpace]Ks1 -> 0.4, SBP\[LetterSpace]A\[LetterSpace]Ks2 -> 0.2, SBPase\[LetterSpace]Vm -> 0.324191, SBPase\[LetterSpace]q -> 666000.0, SBPase\[LetterSpace]Ks1 -> 0.05, SBPase\[LetterSpace]Kr1 -> 12.0, S7P\[LetterSpace]TK\[LetterSpace]Vm -> 3.12207, S7P\[LetterSpace]TK\[LetterSpace]Ks1 -> 0.072, S7P\[LetterSpace]TK\[LetterSpace]Ks2 -> 0.46, S7P\[LetterSpace]TK\[LetterSpace]Kr1 -> 0.1, S7P\[LetterSpace]TK\[LetterSpace]Kr2 -> 1.5, S7P\[LetterSpace]TK\[LetterSpace]q -> 1.17647, Ru5P\[LetterSpace]K\[LetterSpace]Vm -> 10.8348, Ru5P\[LetterSpace]K\[LetterSpace]q -> 6846.0, Ru5P\[LetterSpace]K\[LetterSpace]Ks1 -> 0.05, Ru5P\[LetterSpace]K\[LetterSpace]Kr1 -> 2.0, Ru5P\[LetterSpace]K\[LetterSpace]Kr2 -> 0.7, Ru5P\[LetterSpace]K\[LetterSpace]Kr3 -> 4.0, Ru5P\[LetterSpace]K\[LetterSpace]Kr41 -> 2.5, Ru5P\[LetterSpace]K\[LetterSpace]Ks2 -> 0.059, Ru5P\[LetterSpace]K\[LetterSpace]Kr42 -> 0.4, ATP\[LetterSpace]S\[LetterSpace]q -> 5.734, ATP\[LetterSpace]S\[LetterSpace]Ks1 -> 0.014, ATP\[LetterSpace]S\[LetterSpace]Ks2 -> 0.3, ATP\[LetterSpace]S\[LetterSpace]Kp1 -> 0.3, ATP\[LetterSpace]S\[LetterSpace]Vm -> 15.0, AGPase\[LetterSpace]Vm -> 0.266843, AGPase\[LetterSpace]K1 -> 0.08, AGPase\[LetterSpace]K2 -> 0.08, AGPase\[LetterSpace]KR1 -> 10.0, AGPase\[LetterSpace]KA1 -> 0.1, AGPase\[LetterSpace]KA2 -> 0.02, AGPase\[LetterSpace]KA3 -> 0.02, TPT\[LetterSpace]PGA\[LetterSpace]Vm -> 1.24333, TPT\[LetterSpace]PGA\[LetterSpace]KA -> 0.74, TPT\[LetterSpace]PGA\[LetterSpace]K -> 0.25, TPT\[LetterSpace]PGA\[LetterSpace]KR1 -> 0.63, TPT\[LetterSpace]PGA\[LetterSpace]KR2 -> 0.075, TPT\[LetterSpace]PGA\[LetterSpace]KR3 -> 0.077, TPT\[LetterSpace]GAP\[LetterSpace]Vm -> 1.24333, TPT\[LetterSpace]GAP\[LetterSpace]KA -> 0.74, TPT\[LetterSpace]GAP\[LetterSpace]K -> 0.075, TPT\[LetterSpace]GAP\[LetterSpace]KR1 -> 0.63, TPT\[LetterSpace]GAP\[LetterSpace]KR2 -> 0.25, TPT\[LetterSpace]GAP\[LetterSpace]KR3 -> 0.077, TPT\[LetterSpace]DHAP\[LetterSpace]Vm -> 1.24333, TPT\[LetterSpace]DHAP\[LetterSpace]KA -> 0.74, TPT\[LetterSpace]DHAP\[LetterSpace]K -> 0.077, TPT\[LetterSpace]DHAP\[LetterSpace]KR1 -> 0.63, TPT\[LetterSpace]DHAP\[LetterSpace]KR2 -> 0.25, TPT\[LetterSpace]DHAP\[LetterSpace]KR3 -> 0.075, FBPc\[LetterSpace]A\[LetterSpace]q -> 12.0, FBPc\[LetterSpace]A\[LetterSpace]Ks1 -> 0.3, FBPc\[LetterSpace]A\[LetterSpace]Ks2 -> 0.4, FBPc\[LetterSpace]A\[LetterSpace]Kp1 -> 0.02, FBPc\[LetterSpace]A\[LetterSpace]Vm -> 0.107377, FBPcase\[LetterSpace]q -> 6663.0, FBPcase\[LetterSpace]Kp1 -> 0.7, FBPcase\[LetterSpace]Kp2 -> 12.0, FBPcase\[LetterSpace]Vm -> 0.063979, UGPase\[LetterSpace]q -> 0.31, UGPase\[LetterSpace]Ks1 -> 0.14, UGPase\[LetterSpace]Ks2 -> 0.1, UGPase\[LetterSpace]Kp1 -> 0.12, UGPase\[LetterSpace]Kp2 -> 0.11, UGPase\[LetterSpace]Vm -> 0.115403, SucPc\[LetterSpace]S\[LetterSpace]Vm -> 0.0555034, SucPc\[LetterSpace]S\[LetterSpace]q -> 10.0, SucPc\[LetterSpace]S\[LetterSpace]Ks1 -> 0.8, SucPc\[LetterSpace]S\[LetterSpace]Kr1 -> 0.8, SucPc\[LetterSpace]S\[LetterSpace]Ks2 -> 2.4, SucPc\[LetterSpace]S\[LetterSpace]Kr2 -> 0.7, SucPc\[LetterSpace]S\[LetterSpace]Kr3 -> 0.4, SucPc\[LetterSpace]S\[LetterSpace]Kr4 -> 50.0, SucPc\[LetterSpace]S\[LetterSpace]Kr5 -> 11.0, SucPc\[LetterSpace]P\[LetterSpace]Vm -> 0.555034, SucPc\[LetterSpace]P\[LetterSpace]q -> 780.0, SucPc\[LetterSpace]P\[LetterSpace]Ks1 -> 0.35, SucPc\[LetterSpace]P\[LetterSpace]Kr1 -> 80.0, F6Pc\[LetterSpace]K\[LetterSpace]Vm -> 0.100915, F6Pc\[LetterSpace]K\[LetterSpace]q -> 590.0, F6Pc\[LetterSpace]K\[LetterSpace]Ks1 -> 0.5, F6Pc\[LetterSpace]K\[LetterSpace]Kr1 -> 0.021, F6Pc\[LetterSpace]K\[LetterSpace]Kr2 -> 0.7, F6Pc\[LetterSpace]K\[LetterSpace]Ks2 -> 0.5, F6Pc\[LetterSpace]K\[LetterSpace]Kr3 -> 0.16, F26BPc\[LetterSpace]P\[LetterSpace]Vm -> 0.0168192, F26BPc\[LetterSpace]P\[LetterSpace]Kr1 -> 0.5, F26BPc\[LetterSpace]P\[LetterSpace]Kr2 -> 0.1, F26BPc\[LetterSpace]P\[LetterSpace]Ks1 -> 0.032, Succ\[LetterSpace]Deg\[LetterSpace]Km -> 5.0, Succ\[LetterSpace]Deg\[LetterSpace]V -> 2.0, PGAc\[LetterSpace]Deg\[LetterSpace]Km -> 1.0, PGAc\[LetterSpace]Deg\[LetterSpace]V -> 0.5, GPT\[LetterSpace]GCEAc\[LetterSpace]Km -> 0.39, GPT\[LetterSpace]GCEAc\[LetterSpace]V -> 5.0, GPT\[LetterSpace]GCEAc\[LetterSpace]Ki -> 0.28, GPT\[LetterSpace]GCEA\[LetterSpace]Km -> 0.39, GPT\[LetterSpace]GCEA\[LetterSpace]V -> 5.0, GPT\[LetterSpace]GCEA\[LetterSpace]Ki -> 0.28, GPT\[LetterSpace]GCA\[LetterSpace]Km -> 0.2, GPT\[LetterSpace]GCA\[LetterSpace]V -> 6.0, GPT\[LetterSpace]GCA\[LetterSpace]Ki -> 0.22, GPT\[LetterSpace]GCAc\[LetterSpace]Km -> 0.2, GPT\[LetterSpace]GCAc\[LetterSpace]V -> 6.0, GPT\[LetterSpace]GCAc\[LetterSpace]Ki -> 0.22, RuBisCO\[LetterSpace]O2\[LetterSpace]K -> 0.02, RuBisCO\[LetterSpace]O2\[LetterSpace]KR1 -> 0.84, RuBisCO\[LetterSpace]O2\[LetterSpace]KR2 -> 0.04, RuBisCO\[LetterSpace]O2\[LetterSpace]KR3 -> 0.075, RuBisCO\[LetterSpace]O2\[LetterSpace]KR4 -> 0.9, RuBisCO\[LetterSpace]O2\[LetterSpace]KR5 -> 0.07, PGCA\[LetterSpace]P\[LetterSpace]Vm -> 52.4199, PGCA\[LetterSpace]P\[LetterSpace]Ks1 -> 0.026, PGCA\[LetterSpace]P\[LetterSpace]Kr1 -> 94.0, PGCA\[LetterSpace]P\[LetterSpace]Kr2 -> 2.55, GCEA\[LetterSpace]K\[LetterSpace]Vm -> 5.71579, GCEA\[LetterSpace]K\[LetterSpace]q -> 300.0, GCEA\[LetterSpace]K\[LetterSpace]Ks1 -> 0.21, GCEA\[LetterSpace]K\[LetterSpace]Kr1 -> 0.36, GCEA\[LetterSpace]K\[LetterSpace]Ks2 -> 0.25, GCAc\[LetterSpace]Ox\[LetterSpace]Km -> 0.1, GCAc\[LetterSpace]Ox\[LetterSpace]V -> 1.45611, SERcGOAc\[LetterSpace]AT\[LetterSpace]Vm -> 3.30619, SERcGOAc\[LetterSpace]AT\[LetterSpace]q -> 0.24, SERcGOAc\[LetterSpace]AT\[LetterSpace]Ks1 -> 2.7, SERcGOAc\[LetterSpace]AT\[LetterSpace]Kr1 -> 33.0, SERcGOAc\[LetterSpace]AT\[LetterSpace]Ks2 -> 0.15, GCEA\[LetterSpace]DH\[LetterSpace]Vm -> 10.0098, GCEA\[LetterSpace]DH\[LetterSpace]q -> 250000.0, GCEA\[LetterSpace]DH\[LetterSpace]Ks1 -> 0.09, GCEA\[LetterSpace]DH\[LetterSpace]Kr1 -> 12.0, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Vm -> 2.74582, GLUcGOAc\[LetterSpace]AT\[LetterSpace]q -> 607.0, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Ks1 -> 1.7, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Kr1 -> 2.0, GLUcGOAc\[LetterSpace]AT\[LetterSpace]Ks2 -> 0.15, GLYc\[LetterSpace]DC\[LetterSpace]Vm -> 2.49475, GLYc\[LetterSpace]DC\[LetterSpace]Ks1 -> 6.0, GLYc\[LetterSpace]DC\[LetterSpace]Kr1 -> 4.0, chloroplast -> 1.0, cytosol -> 1.0 }; assignments = { function\[LetterSpace]4[s1_,s2_,p1_,q_,Ks1_,Ks2_,Kp1_,Vm_] -> ((-(p1/q) + s1*s2)*Vm)/(Ks1*Ks2*(p1/Kp1 + (1 + s1/Ks1)*(1 + s2/Ks2))), function\[LetterSpace]6[Vm_,s1_,s2_,p1_,p2_,Ks1_,Ks2_,r1_,Kr1_,r2_,Kr2_,q_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/((Ks1*(1 + r1/Kr1 + r2/Kr2) + s1)*(Ks2 + s2)), function\[LetterSpace]16[Vm_,s1_,r1_,Kr1_,r2_,Kr2_,Ks1_] -> (s1*Vm)/((1 + r1/Kr1)*(1 + r2/Kr2)*(Ks1 + s1)), function\[LetterSpace]5[Vm_,s1_,p1_,p2_,q_,Ks1_,r1_,Kr1_,r2_,Kr2_] -> ((-((p1*p2)/q) + s1)*Vm)/(Ks1*(1 + r1/Kr1 + r2/Kr2) + s1), function\[LetterSpace]17[substrate_,Km_,V_] -> (substrate*V)/(Km + substrate), function\[LetterSpace]18[substrate_,Inhibitor_,Km_,V_,Ki_] -> (substrate*V)/(Km + (Inhibitor*Km)/Ki + substrate), function\[LetterSpace]9[Vm_,s1_,s2_,p1_,p2_,q_,Ks1_,r1_,Kr1_,r2_,Kr2_,r3_,Kr3_,r4_,Kr41_,Ks2_,Kr42_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/((Ks1*(1 + r1/Kr1 + r2/Kr2 + r3/Kr3) + s1)*(Ks2*(1 + r4/Kr42) + (1 + r4/Kr41)*s2)), function\[LetterSpace]12[s1_,p1_,p2_,q_,Ks1_,Kp1_,Kp2_,Vm_] -> ((-((p1*p2)/q) + s1)*Vm)/(Ks1*((1 + p1/Kp1)*(1 + p2/Kp2) + s1/Ks1)), function\[LetterSpace]8[Vm_,s1_,p1_,p2_,q_,Ks1_,r1_,Kr1_] -> ((-((p1*p2)/q) + s1)*Vm)/(Ks1*(1 + r1/Kr1) + s1), function\[LetterSpace]15[Vm_,s1_,s2_,p1_,p2_,q_,Ks1_,r1_,Kr1_,r2_,Kr2_,Ks2_,r3_,Kr3_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/((Ks1*(1 + r1/Kr1)*(1 + r2/Kr2) + s1)*(Ks2*(1 + r3/Kr3) + s2)), function\[LetterSpace]21[Vm_,s1_,s2_,p1_,p2_,q_,Ks1_,r1_,Kr1_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/(Ks1*(1 + r1/Kr1) + s1), function\[LetterSpace]14[Vm_,s1_,s2_,p1_,p2_,q_,Ks1_,r1_,Kr1_,Ks2_,r2_,Kr2_,r3_,Kr3_,r4_,Kr4_,r5_,Kr5_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/((Ks1*(1 + r1/Kr1) + s1)*(Ks2*(1 + r2/Kr2)*(1 + r3/Kr3)*(1 + r4/Kr4)*(1 + r5/Kr5) + s2)), function\[LetterSpace]22[Vm_,s1_,Ks1_,r1_,Kr1_] -> (s1*Vm)/(Ks1*(1 + r1/Kr1) + s1), function\[LetterSpace]3[Vm_,s1_,s2_,K1_,K2_] -> (s1*s2*Vm)/((K1 + s1)*(K2 + s2)), function\[LetterSpace]1[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), function\[LetterSpace]20[Vm_,s1_,s2_,p1_,p2_,q_,Ks1_,r1_,Kr1_,Ks2_] -> ((-((p1*p2)/q) + s1*s2)*Vm)/((Ks1*(1 + r1/Kr1) + s1)*(Ks2 + s2)), function\[LetterSpace]11[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), function\[LetterSpace]7[Vm_,s1_,s2_,p1_,q_,Ks1_,Ks2_] -> ((-(p1/q) + s1*s2)*Vm)/((Ks1 + s1)*(Ks2 + s2)), function\[LetterSpace]19[Vm_,s1_,Ks1_,r1_,Kr1_,r2_,Kr2_] -> (s1*Vm)/(Ks1*(1 + r1/Kr1)*(1 + r2/Kr2) + s1), function\[LetterSpace]10[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)), function\[LetterSpace]2[Vm_,s1_,s2_,Ks1_,Ks2_,r1_,Kr1_] -> (s1*s2*Vm)/((Ks1 + s1)*(Ks2*(1 + r1/Kr1) + s2)), function\[LetterSpace]13[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))), Ru5P -> 0.20029895366218237*PeP[t], Pi -> cP - ATP[t] - E4P[t] - HeP[t] - PeP[t] - PGA[t] - PGCA[t] - S7P[t] - 2*(DPGA[t] + FBP[t] + RuBP[t] + SBP[t]) - TP[t], X5P -> 0.2989536621823617*PeP[t], G1P -> 0.03885361449292246*HeP[t], GAP -> 0.047619047619047616*TP[t], DHAP -> 0.9523809523809523*TP[t], G6P -> 0.6698899050503873*HeP[t], R5P -> 0.5007473841554559*PeP[t], K52a -> 0.0025*(1 + (100000*F26BPc[t])/7), ADPc -> -ATPc + cAc, Wo\[LetterSpace]min -> (0.34967170967999994*O2*(1 - Abs[1 - RuBP[t]/E] + RuBP[t]/E))/(0.222*(1 + 86.95652173913044*CO2) + O2), UDPc -> cUc - UTPc - UDPGc[t], DHAPc -> 0.9523809523809523*TPc[t], G1Pc -> 0.039111090405087465*HePc[t], Wc\[LetterSpace]min -> (1.456965457*CO2*(1 - Abs[1 - RuBP[t]/E] + RuBP[t]/E))/(CO2 + 0.0115*(1 + 4.504504504504505*O2)), PiPic -> -Pic + PiTc, Pic -> 6000000*(-1 + Sqrt[1 + PiTc/3000000]), F6Pc -> 0.29117845745300375*HePc[t], PiTc -> -ATPc + cPc - UTPc - 2*(F26BPc[t] + FBPc[t]) - HePc[t] - PGAc[t] - SucPc[t] - TPc[t], GAPc -> 0.047619047619047616*TPc[t], G6Pc -> 0.6697104521419086*HePc[t], ADP -> cA - ATP[t], F6P -> 0.2912564804566901*HeP[t] }; events = { }; speciesAnnotations = { ADPc[t]->"http://identifiers.org/chebi/CHEBI:16761", ATPc[t]->"http://identifiers.org/chebi/CHEBI:15422", Pic[t]->"http://identifiers.org/chebi/CHEBI:43474", SucPc[t]->"http://identifiers.org/chebi/CHEBI:16308" }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { ATP'[t] == 1.0*ATP\[LetterSpace]S -1.0*PGA\[LetterSpace]K -1.0*Ru5P\[LetterSpace]K -1.0*AGPase -1.0*GCEA\[LetterSpace]K, DPGA'[t] == 1.0*PGA\[LetterSpace]K -1.0*GAP\[LetterSpace]DH, E4P'[t] == 1.0*F6P\[LetterSpace]TK -1.0*SBP\[LetterSpace]A, 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, GCA'[t] == 1.0*GPT\[LetterSpace]GCAc +1.0*PGCA\[LetterSpace]P -1.0*GPT\[LetterSpace]GCA, GCAc'[t] == 1.0*GPT\[LetterSpace]GCA -1.0*GPT\[LetterSpace]GCAc -1.0*GCAc\[LetterSpace]Ox, GCEA'[t] == 1.0*GPT\[LetterSpace]GCEAc -1.0*GPT\[LetterSpace]GCEA -1.0*GCEA\[LetterSpace]K, GCEAc'[t] == 1.0*GPT\[LetterSpace]GCEA +1.0*GCEA\[LetterSpace]DH -1.0*GPT\[LetterSpace]GCEAc, GLYc'[t] == 1.0*SERcGOAc\[LetterSpace]AT +1.0*GLUcGOAc\[LetterSpace]AT -2.0*GLYc\[LetterSpace]DC, GOAc'[t] == 1.0*GCAc\[LetterSpace]Ox -1.0*SERcGOAc\[LetterSpace]AT -1.0*GLUcGOAc\[LetterSpace]AT, HPRc'[t] == 1.0*SERcGOAc\[LetterSpace]AT -1.0*GCEA\[LetterSpace]DH, 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] == 2.0*RuBisCO\[LetterSpace]CO2 +1.0*RuBisCO\[LetterSpace]O2 +1.0*GCEA\[LetterSpace]K -1.0*PGA\[LetterSpace]K -1.0*TPT\[LetterSpace]PGA, PGAc'[t] == 1.0*TPT\[LetterSpace]PGA -1.0*PGAc\[LetterSpace]Deg, PGCA'[t] == 1.0*RuBisCO\[LetterSpace]O2 -1.0*PGCA\[LetterSpace]P, 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]CO2 -1.0*RuBisCO\[LetterSpace]O2, S7P'[t] == 1.0*SBPase -1.0*S7P\[LetterSpace]TK, SBP'[t] == 1.0*SBP\[LetterSpace]A -1.0*SBPase, SERc'[t] == 1.0*GLYc\[LetterSpace]DC -1.0*SERcGOAc\[LetterSpace]AT, SucPc'[t] == 1.0*SucPc\[LetterSpace]S -1.0*SucPc\[LetterSpace]P, Succ'[t] == 1.0*SucPc\[LetterSpace]P -1.0*Succ\[LetterSpace]Deg, TP'[t] == 1.0*GAP\[LetterSpace]DH -2.0*FBP\[LetterSpace]A -1.0*F6P\[LetterSpace]TK -1.0*SBP\[LetterSpace]A -1.0*S7P\[LetterSpace]TK -1.0*TPT\[LetterSpace]GAP -1.0*TPT\[LetterSpace]DHAP, TPc'[t] == 1.0*TPT\[LetterSpace]GAP +1.0*TPT\[LetterSpace]DHAP -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]}]