(* 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 = { species\[LetterSpace]11[t], species\[LetterSpace]13[t], species\[LetterSpace]14[t], species\[LetterSpace]15[t], species\[LetterSpace]19[t], species\[LetterSpace]20[t], species\[LetterSpace]21[t], species\[LetterSpace]24[t], species\[LetterSpace]28[t], species\[LetterSpace]29[t], species\[LetterSpace]30[t], species\[LetterSpace]31[t], species\[LetterSpace]32[t], species\[LetterSpace]33[t], species\[LetterSpace]34[t], species\[LetterSpace]35[t], species\[LetterSpace]36[t], species\[LetterSpace]37[t], species\[LetterSpace]38[t], species\[LetterSpace]39[t], species\[LetterSpace]40[t], species\[LetterSpace]41[t], species\[LetterSpace]42[t], species\[LetterSpace]43[t], species\[LetterSpace]44[t], species\[LetterSpace]45[t], species\[LetterSpace]46[t], species\[LetterSpace]47[t], species\[LetterSpace]48[t], species\[LetterSpace]49[t], species\[LetterSpace]5[t], species\[LetterSpace]50[t], species\[LetterSpace]51[t], species\[LetterSpace]52[t], species\[LetterSpace]9[t] }; initialValues = { species\[LetterSpace]11[0] == 0.0, species\[LetterSpace]13[0] == 0.0, species\[LetterSpace]14[0] == 0.1, species\[LetterSpace]15[0] == 0.0, species\[LetterSpace]19[0] == 1.0, species\[LetterSpace]20[0] == 0.0, species\[LetterSpace]21[0] == 0.0, species\[LetterSpace]24[0] == 0.0, species\[LetterSpace]28[0] == 0.0, species\[LetterSpace]29[0] == 0.0, species\[LetterSpace]30[0] == 0.0, species\[LetterSpace]31[0] == 0.1, species\[LetterSpace]32[0] == 0.0, species\[LetterSpace]33[0] == 0.0, species\[LetterSpace]34[0] == 0.0, species\[LetterSpace]35[0] == 0.0, species\[LetterSpace]36[0] == 0.0, species\[LetterSpace]37[0] == 0.0, species\[LetterSpace]38[0] == 0.0, species\[LetterSpace]39[0] == 0.0, species\[LetterSpace]40[0] == 0.0, species\[LetterSpace]41[0] == 0.0, species\[LetterSpace]42[0] == 0.0, species\[LetterSpace]43[0] == 0.0, species\[LetterSpace]44[0] == 0.0, species\[LetterSpace]45[0] == 0.0, species\[LetterSpace]46[0] == 0.0, species\[LetterSpace]47[0] == 0.0, species\[LetterSpace]48[0] == 0.0, species\[LetterSpace]49[0] == 0.0, species\[LetterSpace]5[0] == 0.0, species\[LetterSpace]50[0] == 0.0, species\[LetterSpace]51[0] == 0.0, species\[LetterSpace]52[0] == 0.0, species\[LetterSpace]9[0] == 1.0 }; rates = { reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]11, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]14, reaction\[LetterSpace]15, reaction\[LetterSpace]16, reaction\[LetterSpace]17, reaction\[LetterSpace]18, reaction\[LetterSpace]19, reaction\[LetterSpace]2, reaction\[LetterSpace]20, reaction\[LetterSpace]21, reaction\[LetterSpace]22, reaction\[LetterSpace]23, reaction\[LetterSpace]24, reaction\[LetterSpace]25, reaction\[LetterSpace]26, reaction\[LetterSpace]27, reaction\[LetterSpace]28, reaction\[LetterSpace]29, reaction\[LetterSpace]3, reaction\[LetterSpace]30, reaction\[LetterSpace]31, reaction\[LetterSpace]32, reaction\[LetterSpace]33, reaction\[LetterSpace]34, reaction\[LetterSpace]35, reaction\[LetterSpace]36, reaction\[LetterSpace]37, reaction\[LetterSpace]38, reaction\[LetterSpace]39, reaction\[LetterSpace]4, reaction\[LetterSpace]40, reaction\[LetterSpace]41, reaction\[LetterSpace]42, reaction\[LetterSpace]43, reaction\[LetterSpace]44, reaction\[LetterSpace]45, reaction\[LetterSpace]46, reaction\[LetterSpace]47, reaction\[LetterSpace]48, reaction\[LetterSpace]49, reaction\[LetterSpace]5, reaction\[LetterSpace]50, reaction\[LetterSpace]51, reaction\[LetterSpace]52, reaction\[LetterSpace]53, reaction\[LetterSpace]54, reaction\[LetterSpace]55, reaction\[LetterSpace]56, reaction\[LetterSpace]57, reaction\[LetterSpace]58, reaction\[LetterSpace]59, reaction\[LetterSpace]6, reaction\[LetterSpace]7, reaction\[LetterSpace]8, reaction\[LetterSpace]9 }; rateEquations = { reaction\[LetterSpace]1 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]1\[LetterSpace]k1*species\[LetterSpace]1*species\[LetterSpace]2*species\[LetterSpace]3*species\[LetterSpace]4*species\[LetterSpace]6*species\[LetterSpace]7*species\[LetterSpace]8 - reaction\[LetterSpace]1\[LetterSpace]k2*species\[LetterSpace]5[t]), reaction\[LetterSpace]10 -> compartment\[LetterSpace]1*reaction\[LetterSpace]10\[LetterSpace]k1*species\[LetterSpace]21[t]^2, reaction\[LetterSpace]11 -> compartment\[LetterSpace]1*reaction\[LetterSpace]11\[LetterSpace]k1*species\[LetterSpace]22, reaction\[LetterSpace]12 -> compartment\[LetterSpace]1*reaction\[LetterSpace]12\[LetterSpace]k1*species\[LetterSpace]24[t], reaction\[LetterSpace]13 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]13\[LetterSpace]k1*species\[LetterSpace]16*species\[LetterSpace]25*species\[LetterSpace]26*species\[LetterSpace]27*species\[LetterSpace]24[t] - reaction\[LetterSpace]13\[LetterSpace]k2*species\[LetterSpace]28[t]), reaction\[LetterSpace]14 -> compartment\[LetterSpace]1*function\[LetterSpace]2[reaction\[LetterSpace]14\[LetterSpace]k, species\[LetterSpace]28[t], species\[LetterSpace]14[t], reaction\[LetterSpace]14\[LetterSpace]Km], reaction\[LetterSpace]15 -> compartment\[LetterSpace]1*(-(reaction\[LetterSpace]15\[LetterSpace]k2*species\[LetterSpace]30[t]) + reaction\[LetterSpace]15\[LetterSpace]k1*species\[LetterSpace]29[t]*species\[LetterSpace]31[t]), reaction\[LetterSpace]16 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]16\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]30[t] - reaction\[LetterSpace]16\[LetterSpace]k2*species\[LetterSpace]32[t]), reaction\[LetterSpace]17 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]17\[LetterSpace]k1*species\[LetterSpace]31[t]*species\[LetterSpace]33[t] - reaction\[LetterSpace]17\[LetterSpace]k2*species\[LetterSpace]34[t]), reaction\[LetterSpace]18 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]18\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]34[t] - reaction\[LetterSpace]18\[LetterSpace]k2*species\[LetterSpace]35[t]), reaction\[LetterSpace]19 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]19\[LetterSpace]k1*species\[LetterSpace]31[t]*species\[LetterSpace]36[t] - reaction\[LetterSpace]19\[LetterSpace]k2*species\[LetterSpace]37[t]), reaction\[LetterSpace]2 -> compartment\[LetterSpace]1*function\[LetterSpace]1[reaction\[LetterSpace]2\[LetterSpace]k1, species\[LetterSpace]5[t], species\[LetterSpace]9[t], species\[LetterSpace]10, reaction\[LetterSpace]2\[LetterSpace]k2, species\[LetterSpace]11[t]], reaction\[LetterSpace]20 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]20\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]37[t] - reaction\[LetterSpace]20\[LetterSpace]k2*species\[LetterSpace]38[t]), reaction\[LetterSpace]21 -> compartment\[LetterSpace]1*reaction\[LetterSpace]21\[LetterSpace]k1*species\[LetterSpace]32[t], reaction\[LetterSpace]22 -> compartment\[LetterSpace]1*reaction\[LetterSpace]22\[LetterSpace]k1*species\[LetterSpace]35[t], reaction\[LetterSpace]23 -> compartment\[LetterSpace]1*reaction\[LetterSpace]23\[LetterSpace]k1*species\[LetterSpace]38[t], reaction\[LetterSpace]24 -> compartment\[LetterSpace]1*reaction\[LetterSpace]24\[LetterSpace]k1*species\[LetterSpace]30[t], reaction\[LetterSpace]25 -> compartment\[LetterSpace]1*reaction\[LetterSpace]25\[LetterSpace]k1*species\[LetterSpace]34[t], reaction\[LetterSpace]26 -> compartment\[LetterSpace]1*reaction\[LetterSpace]26\[LetterSpace]k1*species\[LetterSpace]37[t], reaction\[LetterSpace]27 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]27\[LetterSpace]k1*species\[LetterSpace]31[t] - reaction\[LetterSpace]27\[LetterSpace]k2*species\[LetterSpace]39[t]), reaction\[LetterSpace]28 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]28\[LetterSpace]k1*species\[LetterSpace]39[t]*species\[LetterSpace]40[t] - reaction\[LetterSpace]28\[LetterSpace]k2*species\[LetterSpace]41[t]), reaction\[LetterSpace]29 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]29\[LetterSpace]k1*species\[LetterSpace]39[t]*species\[LetterSpace]42[t] - reaction\[LetterSpace]29\[LetterSpace]k2*species\[LetterSpace]43[t]), reaction\[LetterSpace]3 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]3\[LetterSpace]k1*species\[LetterSpace]12*species\[LetterSpace]11[t] - reaction\[LetterSpace]3\[LetterSpace]k2*species\[LetterSpace]13[t]), reaction\[LetterSpace]30 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]30\[LetterSpace]k1*species\[LetterSpace]39[t]*species\[LetterSpace]44[t] - reaction\[LetterSpace]30\[LetterSpace]k2*species\[LetterSpace]45[t]), reaction\[LetterSpace]31 -> compartment\[LetterSpace]1*reaction\[LetterSpace]31\[LetterSpace]k1*species\[LetterSpace]22, reaction\[LetterSpace]32 -> compartment\[LetterSpace]1*reaction\[LetterSpace]32\[LetterSpace]k1*species\[LetterSpace]39[t]^2, reaction\[LetterSpace]33 -> compartment\[LetterSpace]1*reaction\[LetterSpace]33\[LetterSpace]k1*species\[LetterSpace]46[t], reaction\[LetterSpace]34 -> compartment\[LetterSpace]1*reaction\[LetterSpace]34\[LetterSpace]k1*species\[LetterSpace]22, reaction\[LetterSpace]35 -> compartment\[LetterSpace]1*reaction\[LetterSpace]35\[LetterSpace]k1*species\[LetterSpace]47[t], reaction\[LetterSpace]36 -> compartment\[LetterSpace]1*reaction\[LetterSpace]36\[LetterSpace]k1*species\[LetterSpace]22, reaction\[LetterSpace]37 -> compartment\[LetterSpace]1*reaction\[LetterSpace]37\[LetterSpace]k1*species\[LetterSpace]48[t], reaction\[LetterSpace]38 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]38\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]29[t] - reaction\[LetterSpace]38\[LetterSpace]k2*species\[LetterSpace]49[t]), reaction\[LetterSpace]39 -> compartment\[LetterSpace]1*reaction\[LetterSpace]39\[LetterSpace]k1*species\[LetterSpace]46[t], reaction\[LetterSpace]4 -> compartment\[LetterSpace]1*function\[LetterSpace]2[reaction\[LetterSpace]4\[LetterSpace]k, species\[LetterSpace]13[t], species\[LetterSpace]14[t], reaction\[LetterSpace]4\[LetterSpace]Km], reaction\[LetterSpace]40 -> compartment\[LetterSpace]1*reaction\[LetterSpace]40\[LetterSpace]k1*species\[LetterSpace]29[t], reaction\[LetterSpace]41 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]41\[LetterSpace]k1*species\[LetterSpace]29[t] - reaction\[LetterSpace]41\[LetterSpace]k2*species\[LetterSpace]40[t]), reaction\[LetterSpace]42 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]42\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]33[t] - reaction\[LetterSpace]42\[LetterSpace]k2*species\[LetterSpace]50[t]), reaction\[LetterSpace]43 -> compartment\[LetterSpace]1*reaction\[LetterSpace]43\[LetterSpace]k1*species\[LetterSpace]47[t], reaction\[LetterSpace]44 -> compartment\[LetterSpace]1*reaction\[LetterSpace]44\[LetterSpace]k1*species\[LetterSpace]33[t], reaction\[LetterSpace]45 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]45\[LetterSpace]k1*species\[LetterSpace]33[t] - reaction\[LetterSpace]45\[LetterSpace]k2*species\[LetterSpace]42[t]), reaction\[LetterSpace]46 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]46\[LetterSpace]k1*species\[LetterSpace]15[t]*species\[LetterSpace]36[t] - reaction\[LetterSpace]46\[LetterSpace]k2*species\[LetterSpace]51[t]), reaction\[LetterSpace]47 -> compartment\[LetterSpace]1*reaction\[LetterSpace]47\[LetterSpace]k1*species\[LetterSpace]48[t], reaction\[LetterSpace]48 -> compartment\[LetterSpace]1*reaction\[LetterSpace]48\[LetterSpace]k1*species\[LetterSpace]36[t], reaction\[LetterSpace]49 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]49\[LetterSpace]k1*species\[LetterSpace]36[t] - reaction\[LetterSpace]49\[LetterSpace]k2*species\[LetterSpace]44[t]), reaction\[LetterSpace]5 -> compartment\[LetterSpace]1*reaction\[LetterSpace]5\[LetterSpace]k1*species\[LetterSpace]15[t], reaction\[LetterSpace]50 -> compartment\[LetterSpace]1*(-(reaction\[LetterSpace]50\[LetterSpace]k2*species\[LetterSpace]32[t]) + reaction\[LetterSpace]50\[LetterSpace]k1*species\[LetterSpace]31[t]*species\[LetterSpace]49[t]), reaction\[LetterSpace]51 -> compartment\[LetterSpace]1*reaction\[LetterSpace]51\[LetterSpace]k1*species\[LetterSpace]41[t], reaction\[LetterSpace]52 -> compartment\[LetterSpace]1*(-(reaction\[LetterSpace]52\[LetterSpace]k2*species\[LetterSpace]35[t]) + reaction\[LetterSpace]52\[LetterSpace]k1*species\[LetterSpace]31[t]*species\[LetterSpace]50[t]), reaction\[LetterSpace]53 -> compartment\[LetterSpace]1*reaction\[LetterSpace]53\[LetterSpace]k1*species\[LetterSpace]43[t], reaction\[LetterSpace]54 -> compartment\[LetterSpace]1*(-(reaction\[LetterSpace]54\[LetterSpace]k2*species\[LetterSpace]38[t]) + reaction\[LetterSpace]54\[LetterSpace]k1*species\[LetterSpace]31[t]*species\[LetterSpace]51[t]), reaction\[LetterSpace]55 -> compartment\[LetterSpace]1*reaction\[LetterSpace]55\[LetterSpace]k1*species\[LetterSpace]45[t], reaction\[LetterSpace]56 -> compartment\[LetterSpace]1*reaction\[LetterSpace]56\[LetterSpace]k1*species\[LetterSpace]49[t], reaction\[LetterSpace]57 -> compartment\[LetterSpace]1*reaction\[LetterSpace]57\[LetterSpace]k1*species\[LetterSpace]50[t], reaction\[LetterSpace]58 -> compartment\[LetterSpace]1*reaction\[LetterSpace]58\[LetterSpace]k1*species\[LetterSpace]51[t], reaction\[LetterSpace]59 -> compartment\[LetterSpace]1*reaction\[LetterSpace]59\[LetterSpace]k1*species\[LetterSpace]39[t]^2, reaction\[LetterSpace]6 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]6\[LetterSpace]k1*species\[LetterSpace]1*species\[LetterSpace]16*species\[LetterSpace]17*species\[LetterSpace]18*species\[LetterSpace]3*species\[LetterSpace]4*species\[LetterSpace]53*species\[LetterSpace]8 - reaction\[LetterSpace]6\[LetterSpace]k2*species\[LetterSpace]52[t]), reaction\[LetterSpace]7 -> compartment\[LetterSpace]1*function\[LetterSpace]2[reaction\[LetterSpace]7\[LetterSpace]k, species\[LetterSpace]52[t], species\[LetterSpace]19[t], reaction\[LetterSpace]7\[LetterSpace]Km], reaction\[LetterSpace]8 -> compartment\[LetterSpace]1*reaction\[LetterSpace]8\[LetterSpace]k1*species\[LetterSpace]20[t], reaction\[LetterSpace]9 -> compartment\[LetterSpace]1*(reaction\[LetterSpace]9\[LetterSpace]k1*species\[LetterSpace]20[t] - reaction\[LetterSpace]9\[LetterSpace]k2*species\[LetterSpace]21[t]) }; parameters = { species\[LetterSpace]1 -> 1.0, species\[LetterSpace]10 -> 1.0, species\[LetterSpace]12 -> 1.0, species\[LetterSpace]16 -> 1.0, species\[LetterSpace]17 -> 1.0, species\[LetterSpace]18 -> 1.0, species\[LetterSpace]2 -> 1.0, species\[LetterSpace]22 -> 1.0, species\[LetterSpace]23 -> 0.0, species\[LetterSpace]25 -> 1.0, species\[LetterSpace]26 -> 1.0, species\[LetterSpace]27 -> 1.0, species\[LetterSpace]3 -> 1.0, species\[LetterSpace]4 -> 1.0, species\[LetterSpace]53 -> 1.0, species\[LetterSpace]6 -> 1.0, species\[LetterSpace]7 -> 1.0, species\[LetterSpace]8 -> 1.0, reaction\[LetterSpace]27\[LetterSpace]k2 -> 8*^-05, reaction\[LetterSpace]28\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]28\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]29\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]29\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]30\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]30\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]31\[LetterSpace]k1 -> 1.54*^-06, reaction\[LetterSpace]32\[LetterSpace]k1 -> 0.0165, reaction\[LetterSpace]33\[LetterSpace]k1 -> 0.00028, reaction\[LetterSpace]34\[LetterSpace]k1 -> 1.78*^-07, reaction\[LetterSpace]35\[LetterSpace]k1 -> 0.00028, reaction\[LetterSpace]36\[LetterSpace]k1 -> 1.27*^-07, reaction\[LetterSpace]37\[LetterSpace]k1 -> 0.00028, reaction\[LetterSpace]38\[LetterSpace]k1 -> 0.0225, reaction\[LetterSpace]38\[LetterSpace]k2 -> 0.00125, reaction\[LetterSpace]39\[LetterSpace]k1 -> 0.00408, reaction\[LetterSpace]40\[LetterSpace]k1 -> 0.000113, reaction\[LetterSpace]41\[LetterSpace]k1 -> 0.0003, reaction\[LetterSpace]41\[LetterSpace]k2 -> 0.0002, reaction\[LetterSpace]42\[LetterSpace]k1 -> 0.006, reaction\[LetterSpace]42\[LetterSpace]k2 -> 0.00175, reaction\[LetterSpace]43\[LetterSpace]k1 -> 0.00408, reaction\[LetterSpace]44\[LetterSpace]k1 -> 0.000113, reaction\[LetterSpace]45\[LetterSpace]k1 -> 0.00015, reaction\[LetterSpace]45\[LetterSpace]k2 -> 0.0001, reaction\[LetterSpace]46\[LetterSpace]k1 -> 0.009, reaction\[LetterSpace]46\[LetterSpace]k2 -> 0.00175, reaction\[LetterSpace]47\[LetterSpace]k1 -> 0.00408, reaction\[LetterSpace]48\[LetterSpace]k1 -> 0.000113, reaction\[LetterSpace]49\[LetterSpace]k1 -> 0.00015, reaction\[LetterSpace]49\[LetterSpace]k2 -> 0.0001, reaction\[LetterSpace]50\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]50\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]51\[LetterSpace]k1 -> 0.0138, reaction\[LetterSpace]52\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]52\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]53\[LetterSpace]k1 -> 0.0052, reaction\[LetterSpace]54\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]54\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]55\[LetterSpace]k1 -> 0.0052, reaction\[LetterSpace]56\[LetterSpace]k1 -> 0.00407, reaction\[LetterSpace]57\[LetterSpace]k1 -> 0.0015, reaction\[LetterSpace]58\[LetterSpace]k1 -> 0.0022, reaction\[LetterSpace]59\[LetterSpace]k1 -> 0.001, reaction\[LetterSpace]20\[LetterSpace]k1 -> 0.07, reaction\[LetterSpace]20\[LetterSpace]k2 -> 0.00175, reaction\[LetterSpace]21\[LetterSpace]k1 -> 0.0204, reaction\[LetterSpace]22\[LetterSpace]k1 -> 0.0075, reaction\[LetterSpace]23\[LetterSpace]k1 -> 0.011, reaction\[LetterSpace]24\[LetterSpace]k1 -> 2.25*^-05, reaction\[LetterSpace]25\[LetterSpace]k1 -> 2.25*^-05, reaction\[LetterSpace]26\[LetterSpace]k1 -> 2.25*^-05, reaction\[LetterSpace]27\[LetterSpace]k1 -> 0.09, reaction\[LetterSpace]1\[LetterSpace]k1 -> 0.0001, reaction\[LetterSpace]1\[LetterSpace]k2 -> 0.0001, reaction\[LetterSpace]2\[LetterSpace]k1 -> 0.001, reaction\[LetterSpace]2\[LetterSpace]k2 -> 0.001, reaction\[LetterSpace]3\[LetterSpace]k1 -> 0.003, reaction\[LetterSpace]3\[LetterSpace]k2 -> 0.01, reaction\[LetterSpace]4\[LetterSpace]k -> 0.1, reaction\[LetterSpace]4\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]5\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]6\[LetterSpace]k1 -> 0.0001, reaction\[LetterSpace]6\[LetterSpace]k2 -> 0.0001, reaction\[LetterSpace]7\[LetterSpace]k -> 0.1, reaction\[LetterSpace]7\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]8\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]9\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]9\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]10\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]11\[LetterSpace]k1 -> 0.001, reaction\[LetterSpace]12\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]13\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]13\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]14\[LetterSpace]k -> 0.1, reaction\[LetterSpace]14\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]15\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]15\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]16\[LetterSpace]k1 -> 0.185, reaction\[LetterSpace]16\[LetterSpace]k2 -> 0.0125, reaction\[LetterSpace]17\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]17\[LetterSpace]k2 -> 0.0005, reaction\[LetterSpace]18\[LetterSpace]k1 -> 0.048, reaction\[LetterSpace]18\[LetterSpace]k2 -> 0.00175, reaction\[LetterSpace]19\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]19\[LetterSpace]k2 -> 0.0005, compartment\[LetterSpace]1 -> 1.0 }; assignments = { function\[LetterSpace]2[k_,Enzyme_,Substrate_,Km_] -> (Enzyme*k*Substrate)/(Km + Substrate), function\[LetterSpace]1[k1_,Enzyme_,Substrate1_,Substrate2_,k2_,Product_] -> -(k2*Product) + Enzyme*k1*Substrate1*Substrate2 }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { species\[LetterSpace]11'[t] == 1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]3, species\[LetterSpace]13'[t] == 1.0*reaction\[LetterSpace]3 , species\[LetterSpace]14'[t] == 1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]14, species\[LetterSpace]15'[t] == 1.0*reaction\[LetterSpace]4 +1.0*reaction\[LetterSpace]14 +1.0*reaction\[LetterSpace]21 +1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]23 +1.0*reaction\[LetterSpace]56 +1.0*reaction\[LetterSpace]57 +1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]18 -1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]38 -1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]46, species\[LetterSpace]19'[t] == 1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]7, species\[LetterSpace]20'[t] == 1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]9, species\[LetterSpace]21'[t] == 1.0*reaction\[LetterSpace]9 +2.0*reaction\[LetterSpace]10 -2.0*reaction\[LetterSpace]10, species\[LetterSpace]24'[t] == 1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]11 +1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]13, species\[LetterSpace]28'[t] == 1.0*reaction\[LetterSpace]13 , species\[LetterSpace]29'[t] == 1.0*reaction\[LetterSpace]39 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]38 -1.0*reaction\[LetterSpace]40 -1.0*reaction\[LetterSpace]41, species\[LetterSpace]30'[t] == 1.0*reaction\[LetterSpace]15 +1.0*reaction\[LetterSpace]51 -1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]24, species\[LetterSpace]31'[t] == 1.0*reaction\[LetterSpace]21 +1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]23 +1.0*reaction\[LetterSpace]24 +1.0*reaction\[LetterSpace]25 +1.0*reaction\[LetterSpace]26 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]17 -1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]54, species\[LetterSpace]32'[t] == 1.0*reaction\[LetterSpace]16 +1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]21, species\[LetterSpace]33'[t] == 1.0*reaction\[LetterSpace]43 -1.0*reaction\[LetterSpace]17 -1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]44 -1.0*reaction\[LetterSpace]45, species\[LetterSpace]34'[t] == 1.0*reaction\[LetterSpace]17 +1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]18 -1.0*reaction\[LetterSpace]25, species\[LetterSpace]35'[t] == 1.0*reaction\[LetterSpace]18 +1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]22, species\[LetterSpace]36'[t] == 1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]49, species\[LetterSpace]37'[t] == 1.0*reaction\[LetterSpace]19 +1.0*reaction\[LetterSpace]55 -1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]26, species\[LetterSpace]38'[t] == 1.0*reaction\[LetterSpace]20 +1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]23, species\[LetterSpace]39'[t] == 1.0*reaction\[LetterSpace]27 +2.0*reaction\[LetterSpace]32 +2.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]30 -2.0*reaction\[LetterSpace]32 -2.0*reaction\[LetterSpace]59, species\[LetterSpace]40'[t] == 1.0*reaction\[LetterSpace]41 -1.0*reaction\[LetterSpace]28, species\[LetterSpace]41'[t] == 1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]51, species\[LetterSpace]42'[t] == 1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]29, species\[LetterSpace]43'[t] == 1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]53, species\[LetterSpace]44'[t] == 1.0*reaction\[LetterSpace]49 -1.0*reaction\[LetterSpace]30, species\[LetterSpace]45'[t] == 1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]55, species\[LetterSpace]46'[t] == 1.0*reaction\[LetterSpace]31 +1.0*reaction\[LetterSpace]32 +1.0*reaction\[LetterSpace]39 -1.0*reaction\[LetterSpace]33 -1.0*reaction\[LetterSpace]39, species\[LetterSpace]47'[t] == 1.0*reaction\[LetterSpace]34 +1.0*reaction\[LetterSpace]43 -1.0*reaction\[LetterSpace]35 -1.0*reaction\[LetterSpace]43, species\[LetterSpace]48'[t] == 1.0*reaction\[LetterSpace]36 +1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]37 -1.0*reaction\[LetterSpace]47, species\[LetterSpace]49'[t] == 1.0*reaction\[LetterSpace]38 -1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]56, species\[LetterSpace]5'[t] == 1.0*reaction\[LetterSpace]1 , species\[LetterSpace]50'[t] == 1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]57, species\[LetterSpace]51'[t] == 1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]58, species\[LetterSpace]52'[t] == 1.0*reaction\[LetterSpace]6 , species\[LetterSpace]9'[t] == -1.0*reaction\[LetterSpace]2 }; 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]}]