(* 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 = { AggP1[t], AggP2[t], AggP3[t], AggP4[t], AggP5[t], AggP\[LetterSpace]Proteasome[t], AggPolyQ1[t], AggPolyQ2[t], AggPolyQ3[t], AggPolyQ4[t], AggPolyQ5[t], MisP[t], MisP\[LetterSpace]Proteasome[t], NatP[t], PIdeath[t], PolyQ[t], PolyQ\[LetterSpace]Proteasome[t], Proteasome[t], ROS[t], SeqAggP[t], mRFPu[t], mRFPu\[LetterSpace]Proteasome[t], p38[t], p38\[LetterSpace]P[t], p38death[t] }; initialValues = { AggP1[0] == 0.0, AggP2[0] == 0.0, AggP3[0] == 0.0, AggP4[0] == 0.0, AggP5[0] == 0.0, AggP\[LetterSpace]Proteasome[0] == 0.0, AggPolyQ1[0] == 0.0, AggPolyQ2[0] == 0.0, AggPolyQ3[0] == 0.0, AggPolyQ4[0] == 0.0, AggPolyQ5[0] == 0.0, MisP[0] == 0.0, MisP\[LetterSpace]Proteasome[0] == 0.0, NatP[0] == 19500.0, PIdeath[0] == 0.0, PolyQ[0] == 1000.0, PolyQ\[LetterSpace]Proteasome[0] == 0.0, Proteasome[0] == 1000.0, ROS[0] == 10.0, SeqAggP[0] == 0.0, mRFPu[0] == 300.0, mRFPu\[LetterSpace]Proteasome[0] == 0.0, p38[0] == 100.0, p38\[LetterSpace]P[0] == 0.0, p38death[0] == 0.0 }; rates = { AggP\[LetterSpace]ProteasomeSequestering, MisPInclusionGrowth, MisP\[LetterSpace]Aggregation1, MisP\[LetterSpace]Aggregation2, MisP\[LetterSpace]Aggregation3, MisP\[LetterSpace]Aggregation4, MisP\[LetterSpace]Aggregation5, MisP\[LetterSpace]Degradation, MisP\[LetterSpace]Disaggregation1, MisP\[LetterSpace]Disaggregation2, MisP\[LetterSpace]Disaggregation3, MisP\[LetterSpace]Disaggregation4, MisP\[LetterSpace]Disaggregation5, MisP\[LetterSpace]InclusionFormation, MisP\[LetterSpace]ProteasomeBinding, MisP\[LetterSpace]ProteasomeRelease, MisP\[LetterSpace]ProteasomeSequestering, Misfolding, P38DeathPathway, PIDeathPathway, PolyQAggregation1, PolyQAggregation2, PolyQAggregation3, PolyQAggregation4, PolyQAggregation5, PolyQDegradation, PolyQDisaggregation1, PolyQDisaggregation2, PolyQDisaggregation3, PolyQDisaggregation4, PolyQDisaggregation5, PolyQInclusionFormation, PolyQInclusionGrowth, PolyQ\[LetterSpace]ProteasomeSequestering, ProteasomeInhibition1, ProteasomeInhibition2, ProteasomeInhibition3, ProteasomeInhibition4, ProteasomeInhibition5, ProteasomeInhibitionAggP1, ProteasomeInhibitionAggP2, ProteasomeInhibitionAggP3, ProteasomeInhibitionAggP4, ProteasomeInhibitionAggP5, ProteinSynthesis, ROSgenerationAggPProteasome, ROSgenerationBasal, ROSgenerationSeqAggP, ROSgenerationSmallAggP1, ROSgenerationSmallAggP2, ROSgenerationSmallAggP3, ROSgenerationSmallAggP4, ROSgenerationSmallAggP5, ROSgenerationSmallAggPolyQ1, ROSgenerationSmallAggPolyQ2, ROSgenerationSmallAggPolyQ3, ROSgenerationSmallAggPolyQ4, ROSgenerationSmallAggPolyQ5, ROSremoval, Refolding, mRFPuDegradation, mRFPuProteasomeBinding, mRFPuProteasomeRelease, mRFPuProteasomeSequestering, mRFPuSequestering, mRFPuSynthesis, p38\[LetterSpace]P\[LetterSpace]ROS\[LetterSpace]Generation, p38activation, p38inactivation, polyQSynthesis, polyqProteasomeBinding, polyqProteasomeRelease }; rateEquations = { AggP\[LetterSpace]ProteasomeSequestering -> kalive*kseqAggPProt*AggP\[LetterSpace]Proteasome[t]*SeqAggP[t], MisPInclusionGrowth -> kalive*kseqMisP*MisP[t]*SeqAggP[t], MisP\[LetterSpace]Aggregation1 -> 0.5*kaggMisP*kalive*(-1 + MisP[t])*MisP[t], MisP\[LetterSpace]Aggregation2 -> kagg2MisP*kalive*AggP1[t]*MisP[t], MisP\[LetterSpace]Aggregation3 -> kagg2MisP*kalive*AggP2[t]*MisP[t], MisP\[LetterSpace]Aggregation4 -> kagg2MisP*kalive*AggP3[t]*MisP[t], MisP\[LetterSpace]Aggregation5 -> kagg2MisP*kalive*AggP4[t]*MisP[t], MisP\[LetterSpace]Degradation -> kalive*kdegMisP*kproteff*MisP\[LetterSpace]Proteasome[t], MisP\[LetterSpace]Disaggregation1 -> kalive*kdisaggMisP1*AggP1[t], MisP\[LetterSpace]Disaggregation2 -> kalive*kdisaggMisP2*AggP2[t], MisP\[LetterSpace]Disaggregation3 -> kalive*kdisaggMisP3*AggP3[t], MisP\[LetterSpace]Disaggregation4 -> kalive*kdisaggMisP4*AggP4[t], MisP\[LetterSpace]Disaggregation5 -> kalive*kdisaggMisP5*AggP5[t], MisP\[LetterSpace]InclusionFormation -> kagg2MisP*kalive*AggP5[t]*MisP[t], MisP\[LetterSpace]ProteasomeBinding -> kalive*kbinMisPProt*MisP[t]*Proteasome[t], MisP\[LetterSpace]ProteasomeRelease -> kalive*krelMisPProt*MisP\[LetterSpace]Proteasome[t], MisP\[LetterSpace]ProteasomeSequestering -> kalive*kseqMisPProt*MisP\[LetterSpace]Proteasome[t]*SeqAggP[t], Misfolding -> kalive*kmisfold*NatP[t]*ROS[t], P38DeathPathway -> kalive*kp38act*kp38death*p38\[LetterSpace]P[t], PIDeathPathway -> kalive*kPIdeath*AggP\[LetterSpace]Proteasome[t], PolyQAggregation1 -> (0.5*kaggPolyQ*kalive*(-1 + PolyQ[t])*PolyQ[t]*ROS[t]^2)/(100 + ROS[t]^2), PolyQAggregation2 -> (kaggPolyQ*kalive*AggPolyQ1[t]*PolyQ[t]*ROS[t]^2)/(100 + ROS[t]^2), PolyQAggregation3 -> (kaggPolyQ*kalive*AggPolyQ2[t]*PolyQ[t]*ROS[t]^2)/(100 + ROS[t]^2), PolyQAggregation4 -> (kaggPolyQ*kalive*AggPolyQ3[t]*PolyQ[t]*ROS[t]^2)/(100 + ROS[t]^2), PolyQAggregation5 -> (kaggPolyQ*kalive*AggPolyQ4[t]*PolyQ[t]*ROS[t]^2)/(100 + ROS[t]^2), PolyQDegradation -> kalive*kdegPolyQ*kproteff*PolyQ\[LetterSpace]Proteasome[t], PolyQDisaggregation1 -> kalive*kdisaggPolyQ1*AggPolyQ1[t], PolyQDisaggregation2 -> kalive*kdisaggPolyQ2*AggPolyQ2[t], PolyQDisaggregation3 -> kalive*kdisaggPolyQ3*AggPolyQ3[t], PolyQDisaggregation4 -> kalive*kdisaggPolyQ4*AggPolyQ4[t], PolyQDisaggregation5 -> kalive*kdisaggPolyQ5*AggPolyQ5[t], PolyQInclusionFormation -> kaggPolyQ*kalive*AggPolyQ5[t]*PolyQ[t], PolyQInclusionGrowth -> kalive*kseqPolyQ*PolyQ[t]*SeqAggP[t], PolyQ\[LetterSpace]ProteasomeSequestering -> kalive*kseqPolyQProt*PolyQ\[LetterSpace]Proteasome[t]*SeqAggP[t], ProteasomeInhibition1 -> kalive*kinhprot*AggPolyQ1[t]*Proteasome[t], ProteasomeInhibition2 -> kalive*kinhprot*AggPolyQ2[t]*Proteasome[t], ProteasomeInhibition3 -> kalive*kinhprot*AggPolyQ3[t]*Proteasome[t], ProteasomeInhibition4 -> kalive*kinhprot*AggPolyQ4[t]*Proteasome[t], ProteasomeInhibition5 -> kalive*kinhprot*AggPolyQ5[t]*Proteasome[t], ProteasomeInhibitionAggP1 -> kalive*kinhprot*AggP1[t]*Proteasome[t], ProteasomeInhibitionAggP2 -> kalive*kinhprot*AggP2[t]*Proteasome[t], ProteasomeInhibitionAggP3 -> kalive*kinhprot*AggP3[t]*Proteasome[t], ProteasomeInhibitionAggP4 -> kalive*kinhprot*AggP4[t]*Proteasome[t], ProteasomeInhibitionAggP5 -> kalive*kinhprot*AggP5[t]*Proteasome[t], ProteinSynthesis -> kalive*ksynNatP*Source, ROSgenerationAggPProteasome -> kalive*kgenROSAggP*AggP\[LetterSpace]Proteasome[t], ROSgenerationBasal -> kalive*kgenROS*Source, ROSgenerationSeqAggP -> kalive*kgenROSSeqAggP*SeqAggP[t], ROSgenerationSmallAggP1 -> kalive*kgenROSAggP*AggP1[t], ROSgenerationSmallAggP2 -> kalive*kgenROSAggP*AggP2[t], ROSgenerationSmallAggP3 -> kalive*kgenROSAggP*AggP3[t], ROSgenerationSmallAggP4 -> kalive*kgenROSAggP*AggP4[t], ROSgenerationSmallAggP5 -> kalive*kgenROSAggP*AggP5[t], ROSgenerationSmallAggPolyQ1 -> kalive*kgenROSAggP*AggPolyQ1[t], ROSgenerationSmallAggPolyQ2 -> kalive*kgenROSAggP*AggPolyQ2[t], ROSgenerationSmallAggPolyQ3 -> kalive*kgenROSAggP*AggPolyQ3[t], ROSgenerationSmallAggPolyQ4 -> kalive*kgenROSAggP*AggPolyQ4[t], ROSgenerationSmallAggPolyQ5 -> kalive*kgenROSAggP*AggPolyQ5[t], ROSremoval -> kalive*kremROS*ROS[t], Refolding -> kalive*krefold*MisP[t], mRFPuDegradation -> kalive*kdegmRFPu*kproteff*mRFPu\[LetterSpace]Proteasome[t], mRFPuProteasomeBinding -> kalive*kbinmRFPu*mRFPu[t]*Proteasome[t], mRFPuProteasomeRelease -> kalive*krelmRFPu*mRFPu\[LetterSpace]Proteasome[t], mRFPuProteasomeSequestering -> kalive*kseqmRFPuProt*mRFPu\[LetterSpace]Proteasome[t]*SeqAggP[t], mRFPuSequestering -> kalive*kseqmRFPu*mRFPu[t]*SeqAggP[t], mRFPuSynthesis -> kalive*ksynmRFPu*Source, p38\[LetterSpace]P\[LetterSpace]ROS\[LetterSpace]Generation -> kalive*kgenROSp38*kp38act*p38\[LetterSpace]P[t], p38activation -> kactp38*kalive*p38[t]*ROS[t], p38inactivation -> kalive*kinactp38*p38\[LetterSpace]P[t], polyQSynthesis -> kalive*ksynPolyQ*Source, polyqProteasomeBinding -> kalive*kbinPolyQ*PolyQ[t]*Proteasome[t], polyqProteasomeRelease -> kalive*krelPolyQ*PolyQ\[LetterSpace]Proteasome[t] }; parameters = { kPIdeath -> 2.5*^-08, kactp38 -> 5*^-06, kagg2MisP -> 1*^-10, kaggMisP -> 1*^-11, kaggPolyQ -> 5*^-08, kalive -> 1.0, kbinMisPProt -> 5*^-08, kbinPolyQ -> 5*^-08, kbinmRFPu -> 5*^-07, kdegMisP -> 0.01, kdegPolyQ -> 0.0025, kdegmRFPu -> 0.005, kdisaggMisP1 -> 5*^-07, kdisaggMisP2 -> 4*^-07, kdisaggMisP3 -> 3*^-07, kdisaggMisP4 -> 2*^-07, kdisaggMisP5 -> 1*^-07, kdisaggPolyQ1 -> 5*^-07, kdisaggPolyQ2 -> 4*^-07, kdisaggPolyQ3 -> 3*^-07, kdisaggPolyQ4 -> 2*^-07, kdisaggPolyQ5 -> 1*^-07, kgenROS -> 0.0017, kgenROSAggP -> 5*^-06, kgenROSSeqAggP -> 1*^-07, kgenROSp38 -> 0.0007, kinactp38 -> 0.002, kinhprot -> 5*^-09, kmisfold -> 2*^-06, kp38act -> 1.0, kp38death -> 9*^-08, kproteff -> 1.0, krefold -> 8*^-05, krelMisPProt -> 1*^-08, krelPolyQ -> 1*^-09, krelmRFPu -> 1*^-08, kremROS -> 0.0002, kseqAggPProt -> 5*^-07, kseqMisP -> 1*^-09, kseqMisPProt -> 5*^-07, kseqPolyQ -> 8*^-07, kseqPolyQProt -> 5*^-07, kseqmRFPu -> 1*^-10, kseqmRFPuProt -> 5*^-07, ksynNatP -> 2.4, ksynPolyQ -> 0.007, ksynmRFPu -> 0.138, Sink -> 1.0, Source -> 1.0, cytosol -> 1.0 }; assignments = { oligomers -> AggPolyQ1[t] + AggPolyQ2[t] + AggPolyQ3[t] + AggPolyQ4[t] + AggPolyQ5[t] }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { AggP1'[t] == 1.0*MisP\[LetterSpace]Aggregation1 +1.0*MisP\[LetterSpace]Disaggregation2 +1.0*ROSgenerationSmallAggP1 -1.0*MisP\[LetterSpace]Aggregation2 -1.0*MisP\[LetterSpace]Disaggregation1 -1.0*ProteasomeInhibitionAggP1 -1.0*ROSgenerationSmallAggP1, AggP2'[t] == 1.0*MisP\[LetterSpace]Aggregation2 +1.0*MisP\[LetterSpace]Disaggregation3 +1.0*ROSgenerationSmallAggP2 -1.0*MisP\[LetterSpace]Aggregation3 -1.0*MisP\[LetterSpace]Disaggregation2 -1.0*ProteasomeInhibitionAggP2 -1.0*ROSgenerationSmallAggP2, AggP3'[t] == 1.0*MisP\[LetterSpace]Aggregation3 +1.0*MisP\[LetterSpace]Disaggregation4 +1.0*ROSgenerationSmallAggP3 -1.0*MisP\[LetterSpace]Aggregation4 -1.0*MisP\[LetterSpace]Disaggregation3 -1.0*ProteasomeInhibitionAggP3 -1.0*ROSgenerationSmallAggP3, AggP4'[t] == 1.0*MisP\[LetterSpace]Aggregation4 +1.0*MisP\[LetterSpace]Disaggregation5 +1.0*ROSgenerationSmallAggP4 -1.0*MisP\[LetterSpace]Aggregation5 -1.0*MisP\[LetterSpace]Disaggregation4 -1.0*ProteasomeInhibitionAggP4 -1.0*ROSgenerationSmallAggP4, AggP5'[t] == 1.0*MisP\[LetterSpace]Aggregation5 +1.0*ROSgenerationSmallAggP5 -1.0*MisP\[LetterSpace]Disaggregation5 -1.0*MisP\[LetterSpace]InclusionFormation -1.0*ProteasomeInhibitionAggP5 -1.0*ROSgenerationSmallAggP5, AggP\[LetterSpace]Proteasome'[t] == 1.0*ProteasomeInhibition1 +1.0*ProteasomeInhibition2 +1.0*ProteasomeInhibition3 +1.0*ProteasomeInhibition4 +1.0*ProteasomeInhibition5 +1.0*ROSgenerationAggPProteasome +1.0*ProteasomeInhibitionAggP1 +1.0*ProteasomeInhibitionAggP2 +1.0*ProteasomeInhibitionAggP3 +1.0*ProteasomeInhibitionAggP4 +1.0*ProteasomeInhibitionAggP5 +1.0*PIDeathPathway -1.0*ROSgenerationAggPProteasome -1.0*AggP\[LetterSpace]ProteasomeSequestering -1.0*PIDeathPathway, AggPolyQ1'[t] == 1.0*PolyQAggregation1 +1.0*PolyQDisaggregation2 +1.0*ROSgenerationSmallAggPolyQ1 -1.0*PolyQAggregation2 -1.0*PolyQDisaggregation1 -1.0*ProteasomeInhibition1 -1.0*ROSgenerationSmallAggPolyQ1, AggPolyQ2'[t] == 1.0*PolyQAggregation2 +1.0*PolyQDisaggregation3 +1.0*ROSgenerationSmallAggPolyQ2 -1.0*PolyQAggregation3 -1.0*PolyQDisaggregation2 -1.0*ProteasomeInhibition2 -1.0*ROSgenerationSmallAggPolyQ2, AggPolyQ3'[t] == 1.0*PolyQAggregation3 +1.0*PolyQDisaggregation4 +1.0*ROSgenerationSmallAggPolyQ3 -1.0*PolyQAggregation4 -1.0*PolyQDisaggregation3 -1.0*ProteasomeInhibition3 -1.0*ROSgenerationSmallAggPolyQ3, AggPolyQ4'[t] == 1.0*PolyQAggregation4 +1.0*PolyQDisaggregation5 +1.0*ROSgenerationSmallAggPolyQ4 -1.0*PolyQAggregation5 -1.0*PolyQDisaggregation4 -1.0*ProteasomeInhibition4 -1.0*ROSgenerationSmallAggPolyQ4, AggPolyQ5'[t] == 1.0*PolyQAggregation5 +1.0*ROSgenerationSmallAggPolyQ5 -1.0*PolyQDisaggregation5 -1.0*PolyQInclusionFormation -1.0*ProteasomeInhibition5 -1.0*ROSgenerationSmallAggPolyQ5, MisP'[t] == 1.0*Misfolding +1.0*MisP\[LetterSpace]ProteasomeRelease +2.0*MisP\[LetterSpace]Disaggregation1 +1.0*MisP\[LetterSpace]Disaggregation2 +1.0*MisP\[LetterSpace]Disaggregation3 +1.0*MisP\[LetterSpace]Disaggregation4 +1.0*MisP\[LetterSpace]Disaggregation5 -1.0*Refolding -1.0*MisP\[LetterSpace]ProteasomeBinding -2.0*MisP\[LetterSpace]Aggregation1 -1.0*MisP\[LetterSpace]Aggregation2 -1.0*MisP\[LetterSpace]Aggregation3 -1.0*MisP\[LetterSpace]Aggregation4 -1.0*MisP\[LetterSpace]Aggregation5 -1.0*MisP\[LetterSpace]InclusionFormation -1.0*MisPInclusionGrowth, MisP\[LetterSpace]Proteasome'[t] == 1.0*MisP\[LetterSpace]ProteasomeBinding -1.0*MisP\[LetterSpace]ProteasomeSequestering -1.0*MisP\[LetterSpace]ProteasomeRelease -1.0*MisP\[LetterSpace]Degradation, NatP'[t] == 1.0*ProteinSynthesis +1.0*Refolding -1.0*Misfolding, PIdeath'[t] == 1.0*PIDeathPathway , PolyQ'[t] == 1.0*polyQSynthesis +1.0*polyqProteasomeRelease +1.0*PolyQDisaggregation5 +1.0*PolyQDisaggregation4 +1.0*PolyQDisaggregation3 +1.0*PolyQDisaggregation2 +2.0*PolyQDisaggregation1 -1.0*polyqProteasomeBinding -2.0*PolyQAggregation1 -1.0*PolyQAggregation2 -1.0*PolyQAggregation3 -1.0*PolyQAggregation4 -1.0*PolyQAggregation5 -1.0*PolyQInclusionFormation -1.0*PolyQInclusionGrowth, PolyQ\[LetterSpace]Proteasome'[t] == 1.0*polyqProteasomeBinding -1.0*polyqProteasomeRelease -1.0*PolyQDegradation -1.0*PolyQ\[LetterSpace]ProteasomeSequestering, Proteasome'[t] == 1.0*polyqProteasomeRelease +1.0*PolyQDegradation +1.0*mRFPuProteasomeRelease +1.0*mRFPuDegradation +1.0*MisP\[LetterSpace]ProteasomeRelease +1.0*MisP\[LetterSpace]Degradation -1.0*polyqProteasomeBinding -1.0*mRFPuProteasomeBinding -1.0*ProteasomeInhibition1 -1.0*ProteasomeInhibition2 -1.0*ProteasomeInhibition3 -1.0*ProteasomeInhibition4 -1.0*ProteasomeInhibition5 -1.0*MisP\[LetterSpace]ProteasomeBinding -1.0*ProteasomeInhibitionAggP1 -1.0*ProteasomeInhibitionAggP2 -1.0*ProteasomeInhibitionAggP3 -1.0*ProteasomeInhibitionAggP4 -1.0*ProteasomeInhibitionAggP5, ROS'[t] == 1.0*PolyQAggregation1 +1.0*PolyQAggregation2 +1.0*PolyQAggregation3 +1.0*PolyQAggregation4 +1.0*PolyQAggregation5 +1.0*ROSgenerationBasal +1.0*ROSgenerationSmallAggPolyQ1 +1.0*ROSgenerationSmallAggPolyQ2 +1.0*ROSgenerationSmallAggPolyQ3 +1.0*ROSgenerationSmallAggPolyQ4 +1.0*ROSgenerationSmallAggPolyQ5 +1.0*ROSgenerationAggPProteasome +1.0*p38activation +1.0*Misfolding +1.0*ROSgenerationSmallAggP1 +1.0*ROSgenerationSmallAggP2 +1.0*ROSgenerationSmallAggP3 +1.0*ROSgenerationSmallAggP4 +1.0*ROSgenerationSmallAggP5 +1.0*p38\[LetterSpace]P\[LetterSpace]ROS\[LetterSpace]Generation +1.0*ROSgenerationSeqAggP -1.0*PolyQAggregation1 -1.0*PolyQAggregation2 -1.0*PolyQAggregation3 -1.0*PolyQAggregation4 -1.0*PolyQAggregation5 -1.0*ROSremoval -1.0*p38activation -1.0*Misfolding, SeqAggP'[t] == 7.0*PolyQInclusionFormation +2.0*PolyQInclusionGrowth +2.0*mRFPuProteasomeSequestering +2.0*mRFPuSequestering +2.0*AggP\[LetterSpace]ProteasomeSequestering +2.0*PolyQ\[LetterSpace]ProteasomeSequestering +2.0*MisP\[LetterSpace]ProteasomeSequestering +7.0*MisP\[LetterSpace]InclusionFormation +2.0*MisPInclusionGrowth +1.0*ROSgenerationSeqAggP -1.0*PolyQInclusionGrowth -1.0*mRFPuProteasomeSequestering -1.0*mRFPuSequestering -1.0*AggP\[LetterSpace]ProteasomeSequestering -1.0*PolyQ\[LetterSpace]ProteasomeSequestering -1.0*MisP\[LetterSpace]ProteasomeSequestering -1.0*MisPInclusionGrowth -1.0*ROSgenerationSeqAggP, mRFPu'[t] == 1.0*mRFPuSynthesis +1.0*mRFPuProteasomeRelease -1.0*mRFPuProteasomeBinding -1.0*mRFPuSequestering, mRFPu\[LetterSpace]Proteasome'[t] == 1.0*mRFPuProteasomeBinding -1.0*mRFPuProteasomeRelease -1.0*mRFPuDegradation -1.0*mRFPuProteasomeSequestering, p38'[t] == 1.0*p38inactivation -1.0*p38activation, p38\[LetterSpace]P'[t] == 1.0*p38activation +1.0*p38\[LetterSpace]P\[LetterSpace]ROS\[LetterSpace]Generation +1.0*P38DeathPathway -1.0*p38inactivation -1.0*p38\[LetterSpace]P\[LetterSpace]ROS\[LetterSpace]Generation -1.0*P38DeathPathway, p38death'[t] == 1.0*P38DeathPathway }; 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]}]