(* 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 = { s1[t], s107[t], s108[t], s109[t], s110[t], s111[t], s112[t], s113[t], s17[t], s18[t], s2[t], s20[t], s200[t], s21[t], s211[t], s213[t], s22[t], s23[t], s24[t], s25[t], s26[t], s27[t], s29[t], s30[t], s31[t], s32[t], s33[t], s35[t], s381[t], s383[t], s385[t], s387[t], s389[t], s391[t], s393[t], s445[t], s446[t], s447[t], s448[t], s451[t], s473[t], s474[t], s475[t], s476[t], s477[t], s478[t], s479[t], s482[t], s483[t], s484[t], s489[t], s490[t], s491[t], s492[t], s493[t], s494[t], s495[t], s496[t], s498[t], s499[t], s5[t], s500[t], s501[t], s502[t], s51[t], s517[t], s518[t], s519[t], s52[t], s520[t], s521[t], s522[t], s523[t], s524[t], s525[t], s526[t], s527[t], s528[t], s529[t], s53[t], s530[t], s531[t], s533[t], s535[t], s536[t], s6[t], s7[t], s78[t], s85[t] }; initialValues = { s1[0] == 8.0, s107[0] == 0.0, s108[0] == 0.0, s109[0] == 0.0, s110[0] == 0.0, s111[0] == 0.0, s112[0] == 0.0, s113[0] == 0.0, s17[0] == 1650.0, s18[0] == 22.0, s2[0] == 8.0, s20[0] == 8.0, s200[0] == 0.0, s21[0] == 8.0, s211[0] == 0.0, s213[0] == 0.0, s22[0] == 750.0, s23[0] == 8.0, s24[0] == 8.0, s25[0] == 8.0, s26[0] == 8.0, s27[0] == 8.0, s29[0] == 8.0, s30[0] == 8.0, s31[0] == 8.0, s32[0] == 8.0, s33[0] == 0.0, s35[0] == 1500.0, s381[0] == 0.0, s383[0] == 0.0, s385[0] == 0.0, s387[0] == 0.0, s389[0] == 0.0, s391[0] == 0.0, s393[0] == 0.0, s445[0] == 0.0, s446[0] == 0.0, s447[0] == 0.0, s448[0] == 0.0, s451[0] == 0.0, s473[0] == 0.0, s474[0] == 0.0, s475[0] == 0.0, s476[0] == 0.0, s477[0] == 0.0, s478[0] == 0.0, s479[0] == 0.0, s482[0] == 0.0, s483[0] == 0.0, s484[0] == 0.0, s489[0] == 0.0, s490[0] == 0.0, s491[0] == 0.0, s492[0] == 0.0, s493[0] == 0.0, s494[0] == 0.0, s495[0] == 0.0, s496[0] == 0.0, s498[0] == 0.0, s499[0] == 0.0, s5[0] == 8.0, s500[0] == 0.0, s501[0] == 0.0, s502[0] == 0.0, s51[0] == 200.0, s517[0] == 0.0, s518[0] == 0.0, s519[0] == 0.0, s52[0] == 0.0, s520[0] == 0.0, s521[0] == 0.0, s522[0] == 0.0, s523[0] == 0.0, s524[0] == 0.0, s525[0] == 0.0, s526[0] == 0.0, s527[0] == 0.0, s528[0] == 0.0, s529[0] == 0.0, s53[0] == 0.0, s530[0] == 0.0, s531[0] == 0.0, s533[0] == 0.0, s535[0] == 0.0, s536[0] == 0.0, s6[0] == 0.0, s7[0] == 953.0, s78[0] == 0.0, s85[0] == 0.0 }; rates = { re1, re10, re100, re101, re102, re103, re104, re105, re106, re107, re108, re109, re11, re110, re111, re112, re113, re114, re115, re116, re117, re118, re119, re12, re120, re121, re122, re123, re124, re125, re126, re127, re128, re129, re13, re130, re131, re132, re133, re134, re135, re136, re137, re14, re15, re16, re17, re18, re19, re2, re20, re21, re22, re23, re24, re25, re26, re27, re28, re29, re3, re30, re31, re32, re33, re34, re35, re36, re37, re38, re4, re40, re41, re42, re43, re44, re45, re46, re47, re48, re49, re5, re50, re51, re52, re53, re54, re55, re56, re57, re58, re59, re6, re60, re61, re62, re63, re64, re65, re66, re67, re68, re69, re7, re70, re71, re72, re73, re74, re75, re76, re77, re78, re79, re8, re80, re81, re82, re83, re84, re85, re86, re87, re88, re89, re9, re90, re91, re92, re93, re94, re95, re96, re97, re98, re99 }; rateEquations = { re1 -> c1*re1\[LetterSpace]k1*s3, re10 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s20[t], re100 -> k\[LetterSpace]OligomerForm*s17[t]*s499[t], re101 -> k\[LetterSpace]OligomerForm*s17[t]*s500[t], re102 -> c1*k\[LetterSpace]ProtOligDegr*s381[t], re103 -> c1*k\[LetterSpace]ProtOligDegr*s383[t], re104 -> c1*k\[LetterSpace]ProtOligDegr*s385[t], re105 -> c1*k\[LetterSpace]ProtOligDegr*s387[t], re106 -> c1*k\[LetterSpace]ProtOligDegr*s389[t], re107 -> c1*k\[LetterSpace]ProtOligDegr*s391[t], re108 -> c1*k\[LetterSpace]ProtOligDegr*s393[t], re109 -> c1*k\[LetterSpace]ProtOligDegr*s473[t], re11 -> c1*k\[LetterSpace]ProteasomeBind*s20[t]*s35[t], re110 -> c1*k\[LetterSpace]ProtOligDegr*s474[t], re111 -> c1*k\[LetterSpace]ProtOligDegr*s475[t], re112 -> c1*k\[LetterSpace]ProtOligDegr*s476[t], re113 -> c1*k\[LetterSpace]ProtOligDegr*s477[t], re114 -> c1*k\[LetterSpace]ProtOligDegr*s478[t], re115 -> c1*k\[LetterSpace]ProtOligDegr*s479[t], re116 -> c1*k\[LetterSpace]ProteasomeBind*s33[t]*s35[t], re117 -> c1*k\[LetterSpace]DisRate*s27[t], re118 -> c1*k\[LetterSpace]DisRate*s26[t], re119 -> c1*k\[LetterSpace]DisRate*s25[t], re12 -> k\[LetterSpace]OligAutophagUptake*s24[t], re120 -> c1*k\[LetterSpace]DisRate*s21[t], re121 -> c1*k\[LetterSpace]DisRate*s2[t], re122 -> c1*k\[LetterSpace]DisRate*s1[t], re123 -> c1*k\[LetterSpace]DisRate*s5[t], re124 -> c1*k\[LetterSpace]DisRate*s6[t], re125 -> c1*k\[LetterSpace]DisRate*s29[t], re126 -> c1*k\[LetterSpace]DisRate*s30[t], re127 -> c1*k\[LetterSpace]DisRate*s31[t], re128 -> c1*k\[LetterSpace]DisRate*s32[t], re129 -> c1*k\[LetterSpace]DisRate*s23[t], re13 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s24[t], re130 -> c1*k\[LetterSpace]DisRate*s24[t], re131 -> c1*k\[LetterSpace]DisRate*s20[t], re132 -> c1*k\[LetterSpace]DisRate*s18[t], re133 -> c1*re133\[LetterSpace]k1*s17[t]*s33[t], re134 -> k\[LetterSpace]OligAutophagUptake*s17[t], re135 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s533[t], re136 -> k\[LetterSpace]OligAutophagUptake*s7[t], re137 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s535[t], re14 -> c1*k\[LetterSpace]ProteasomeBind*s24[t]*s35[t], re15 -> k\[LetterSpace]OligAutophagUptake*s23[t], re16 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s23[t], re17 -> c1*k\[LetterSpace]ProteasomeBind*s23[t]*s35[t], re18 -> k\[LetterSpace]WTOligBindOnLamp*s23[t]*s51[t], re19 -> k\[LetterSpace]OligAutophagUptake*s32[t], re2 -> c1*re2\[LetterSpace]k1*s3, re20 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s32[t], re21 -> c1*k\[LetterSpace]ProteasomeBind*s32[t]*s35[t], re22 -> k\[LetterSpace]OligAutophagUptake*s31[t], re23 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s31[t], re24 -> c1*k\[LetterSpace]ProteasomeBind*s31[t]*s35[t], re25 -> k\[LetterSpace]WTOligBindOnLamp*s31[t]*s51[t], re26 -> k\[LetterSpace]OligAutophagUptake*s30[t], re27 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s30[t], re28 -> c1*k\[LetterSpace]ProteasomeBind*s30[t]*s35[t], re29 -> c1*k\[LetterSpace]ProteasomeBind*s29[t]*s35[t], re3 -> c1*k\[LetterSpace]2merForm*s17[t]^2, re30 -> k\[LetterSpace]WTOligBindOnLamp*s29[t]*s51[t], re31 -> c1*re31\[LetterSpace]k1*s22[t], re32 -> c1*k\[LetterSpace]2merForm*s7[t]^2, re33 -> k\[LetterSpace]2merForm*s536[t]*s7[t], re34 -> k\[LetterSpace]OligomerForm*s490[t]*s7[t], re35 -> k\[LetterSpace]OligomerForm*s489[t]*s7[t], re36 -> k\[LetterSpace]OligomerForm*s492[t]*s7[t], re37 -> c3*re37\[LetterSpace]k1*s78[t], re38 -> c3*re38\[LetterSpace]k1*s85[t], re4 -> c1*re4\[LetterSpace]k1*s17[t]*s22[t], re40 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s517[t], re41 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s523[t], re42 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s520[t], re43 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s521[t], re44 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s522[t], re45 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s518[t], re46 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s519[t], re47 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s29[t], re48 -> k\[LetterSpace]OligAutophagUptake*s6[t], re49 -> c1*k\[LetterSpace]OligomerForm*s6[t]*s7[t], re5 -> k\[LetterSpace]WTasyn1\[LetterSpace]2merBindOnLamp*s17[t]*s51[t], re50 -> k\[LetterSpace]OligAutophagUptake*s5[t], re51 -> c1*k\[LetterSpace]OligomerForm*s5[t]*s7[t], re52 -> c1*k\[LetterSpace]ProteasomeBind*s35[t]*s5[t], re53 -> k\[LetterSpace]OligAutophagUptake*s2[t], re54 -> c1*k\[LetterSpace]OligomerForm*s2[t]*s7[t], re55 -> c1*k\[LetterSpace]ProteasomeBind*s2[t]*s35[t], re56 -> k\[LetterSpace]OligAutophagUptake*s1[t], re57 -> c1*k\[LetterSpace]OligomerForm*s1[t]*s7[t], re58 -> c1*k\[LetterSpace]ProteasomeBind*s1[t]*s35[t], re59 -> k\[LetterSpace]OligAutophagUptake*s21[t], re6 -> k\[LetterSpace]OligAutophagUptake*s18[t], re60 -> c1*k\[LetterSpace]OligomerForm*s21[t]*s7[t], re61 -> c1*k\[LetterSpace]ProteasomeBind*s21[t]*s35[t], re62 -> k\[LetterSpace]OligAutophagUptake*s25[t], re63 -> c1*k\[LetterSpace]OligomerForm*s25[t]*s7[t], re64 -> c1*k\[LetterSpace]ProteasomeBind*s25[t]*s35[t], re65 -> k\[LetterSpace]OligAutophagUptake*s26[t], re66 -> c1*k\[LetterSpace]OligomerForm*s26[t]*s7[t], re67 -> c1*k\[LetterSpace]ProteasomeBind*s26[t]*s35[t], re68 -> c1*k\[LetterSpace]ProteasomeBind*s27[t]*s35[t], re69 -> c3*re69\[LetterSpace]k1*s53[t], re7 -> c1*k\[LetterSpace]OligomerForm*s17[t]*s18[t], re70 -> c3*re70\[LetterSpace]k1*s52[t], re71 -> k\[LetterSpace]LampFreeWTasyn*s501[t], re72 -> k\[LetterSpace]OligomerForm*s482[t]*s7[t], re73 -> k\[LetterSpace]OligomerForm*s483[t]*s7[t], re74 -> k\[LetterSpace]OligomerForm*s484[t]*s7[t], re75 -> k\[LetterSpace]OligomerForm*s491[t]*s7[t], re76 -> k\[LetterSpace]LampFreeWTasyn*s494[t], re77 -> k\[LetterSpace]LampFreeWTasyn*s495[t], re78 -> k\[LetterSpace]LampFreeWTasyn*s496[t], re79 -> k\[LetterSpace]LampFreeWTasyn*s498[t], re8 -> k\[LetterSpace]WTasyn1\[LetterSpace]2merBindOnLamp*s18[t]*s51[t], re80 -> k\[LetterSpace]LampFreeWTasyn*s499[t], re81 -> k\[LetterSpace]LampFreeWTasyn*s500[t], re82 -> k\[LetterSpace]WTOligBindOnLamp*s30[t]*s500[t], re83 -> k\[LetterSpace]WTOligBindOnLamp*s20[t]*s51[t], re84 -> k\[LetterSpace]WTOligBindOnLamp*s24[t]*s51[t], re85 -> k\[LetterSpace]WTOligBindOnLamp*s32[t]*s51[t], re86 -> k\[LetterSpace]DopModWTasynLampBind*s51[t]*s7[t], re87 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s530[t], re88 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s531[t], re89 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s527[t], re9 -> k\[LetterSpace]OligAutophagUptake*s20[t], re90 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s529[t], re91 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s528[t], re92 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s526[t], re93 -> c2*k\[LetterSpace]M\[LetterSpace]autophagyDegr*s525[t], re94 -> k\[LetterSpace]2merForm*s17[t]*s78[t], re95 -> k\[LetterSpace]OligomerForm*s17[t]*s85[t], re96 -> k\[LetterSpace]OligomerForm*s17[t]*s494[t], re97 -> k\[LetterSpace]OligomerForm*s17[t]*s495[t], re98 -> k\[LetterSpace]OligomerForm*s17[t]*s496[t], re99 -> k\[LetterSpace]OligomerForm*s17[t]*s498[t] }; parameters = { k\[LetterSpace]2merForm -> 1.462941015*^-09, k\[LetterSpace]DisRate -> 4.999533748*^-07, k\[LetterSpace]DopModWTasynLampBind -> 7.6715997*^-09, k\[LetterSpace]LampFreeWTasyn -> 0.0003044571674, k\[LetterSpace]M\[LetterSpace]autophagyDegr -> 0.1, k\[LetterSpace]OligAutophagUptake -> 2.39034347*^-08, k\[LetterSpace]OligomerForm -> 3.350497192*^-08, k\[LetterSpace]ProtOligDegr -> 0.000370096, k\[LetterSpace]ProteasomeBind -> 3.424693672*^-09, k\[LetterSpace]WTOligBindOnLamp -> 4*^-06, k\[LetterSpace]WTasyn1\[LetterSpace]2merBindOnLamp -> 6.865455081*^-07, s3 -> 1.0, re37\[LetterSpace]k1 -> 0.00999558, re38\[LetterSpace]k1 -> 0.00995312, re69\[LetterSpace]k1 -> 0.1, re70\[LetterSpace]k1 -> 0.1, re133\[LetterSpace]k1 -> 4.90556*^-07, re1\[LetterSpace]k1 -> 0.0294219, re2\[LetterSpace]k1 -> 0.0791823, re4\[LetterSpace]k1 -> 6.74768*^-07, re31\[LetterSpace]k1 -> 0.00679501, c1 -> 1.0, c2 -> 1.0, c3 -> 1.0 }; assignments = { Total\[LetterSpace]Cytosolic\[LetterSpace]WTASYN\[LetterSpace]Oligomers -> c1*s1[t] + c1*s2[t] + c1*s20[t] + c1*s21[t] + c1*s23[t] + c1*s24[t] + c1*s25[t] + c1*s26[t] + c1*s27[t] + c1*s29[t] + c1*s30[t] + c1*s31[t] + c1*s32[t] + c1*s5[t], Total\[LetterSpace]Cytosolic\[LetterSpace]WTASYN\[LetterSpace]Monomer -> c1*s17[t] + c1*s7[t], Total\[LetterSpace]Cytosolic\[LetterSpace]WTASYN\[LetterSpace]Dimer -> c1*s18[t] + c1*s6[t] }; events = { }; speciesAnnotations = { s17[t]->"http://identifiers.org/uniprot/P37840", s22[t]->"http://identifiers.org/chebi/CHEBI:18243", s51[t]->"http://identifiers.org/uniprot/P13473", s517[t]->"http://identifiers.org/uniprot/P37840", s518[t]->"http://identifiers.org/uniprot/P37840", s519[t]->"http://identifiers.org/uniprot/P37840", s520[t]->"http://identifiers.org/uniprot/P37840", s521[t]->"http://identifiers.org/uniprot/P37840", s522[t]->"http://identifiers.org/uniprot/P37840", s523[t]->"http://identifiers.org/uniprot/P37840", s533[t]->"http://identifiers.org/uniprot/P37840" }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { s1'[t] == 1.0*re54 +1.0*re120 -1.0*re56 -1.0*re57 -1.0*re58 -1.0*re122, s107'[t] == 1.0*re43 +1.0*re135 , s108'[t] == 1.0*re42 , s109'[t] == 1.0*re40 , s110'[t] == 1.0*re45 , s111'[t] == 1.0*re46 , s112'[t] == 1.0*re44 , s113'[t] == 1.0*re41 , s17'[t] == 1.0*re1 +1.0*re125 +1.0*re126 +1.0*re127 +1.0*re128 +1.0*re129 +1.0*re130 +1.0*re131 +2.0*re132 -2.0*re3 -1.0*re4 -1.0*re5 -1.0*re7 -1.0*re10 -1.0*re13 -1.0*re16 -1.0*re20 -1.0*re23 -1.0*re27 -1.0*re47 -1.0*re94 -1.0*re95 -1.0*re96 -1.0*re97 -1.0*re98 -1.0*re99 -1.0*re100 -1.0*re101 -1.0*re133 -1.0*re134, s18'[t] == 1.0*re3 +1.0*re131 -1.0*re6 -1.0*re7 -1.0*re8 -1.0*re132, s2'[t] == 1.0*re51 +1.0*re122 -1.0*re53 -1.0*re54 -1.0*re55 -1.0*re121, s20'[t] == 1.0*re7 +1.0*re76 +1.0*re130 -1.0*re9 -1.0*re10 -1.0*re11 -1.0*re83 -1.0*re131, s200'[t] == 1.0*re31 , s21'[t] == 1.0*re57 +1.0*re119 -1.0*re59 -1.0*re60 -1.0*re61 -1.0*re120, s211'[t] == 1.0*re69 , s213'[t] == 1.0*re70 , s22'[t] == 1.0*re2 -1.0*re4 -1.0*re31, s23'[t] == 1.0*re13 +1.0*re78 +1.0*re128 -1.0*re15 -1.0*re16 -1.0*re17 -1.0*re18 -1.0*re129, s24'[t] == 1.0*re10 +1.0*re77 +1.0*re129 -1.0*re12 -1.0*re13 -1.0*re14 -1.0*re84 -1.0*re130, s25'[t] == 1.0*re60 +1.0*re118 -1.0*re62 -1.0*re63 -1.0*re64 -1.0*re119, s26'[t] == 1.0*re63 +1.0*re117 -1.0*re65 -1.0*re66 -1.0*re67 -1.0*re118, s27'[t] == 1.0*re66 -1.0*re68 -1.0*re117, s29'[t] == 1.0*re27 +1.0*re71 -1.0*re29 -1.0*re30 -1.0*re47 -1.0*re125, s30'[t] == 1.0*re23 +1.0*re81 +1.0*re125 -1.0*re26 -1.0*re27 -1.0*re28 -1.0*re82 -1.0*re126, s31'[t] == 1.0*re20 +1.0*re80 +1.0*re126 -1.0*re22 -1.0*re23 -1.0*re24 -1.0*re25 -1.0*re127, s32'[t] == 1.0*re16 +1.0*re79 +1.0*re127 -1.0*re19 -1.0*re20 -1.0*re21 -1.0*re85 -1.0*re128, s33'[t] == 1.0*re47 +1.0*re133 -1.0*re116 -1.0*re133, s35'[t] == 1.0*re102 +1.0*re103 +1.0*re104 +1.0*re105 +1.0*re106 +1.0*re107 +1.0*re108 +1.0*re109 +1.0*re110 +1.0*re111 +1.0*re112 +1.0*re113 +1.0*re114 +1.0*re115 -1.0*re11 -1.0*re14 -1.0*re17 -1.0*re21 -1.0*re24 -1.0*re28 -1.0*re29 -1.0*re52 -1.0*re55 -1.0*re58 -1.0*re61 -1.0*re64 -1.0*re67 -1.0*re68 -1.0*re116, s381'[t] == 1.0*re11 -1.0*re102, s383'[t] == 1.0*re14 -1.0*re103, s385'[t] == 1.0*re17 -1.0*re104, s387'[t] == 1.0*re21 -1.0*re105, s389'[t] == 1.0*re24 -1.0*re106, s391'[t] == 1.0*re28 -1.0*re107, s393'[t] == 1.0*re29 -1.0*re108, s445'[t] == 1.0*re89 +1.0*re137 , s446'[t] == 1.0*re88 , s447'[t] == 1.0*re87 , s448'[t] == 1.0*re90 , s451'[t] == 1.0*re93 , s473'[t] == 1.0*re52 -1.0*re109, s474'[t] == 1.0*re55 -1.0*re110, s475'[t] == 1.0*re58 -1.0*re111, s476'[t] == 1.0*re61 -1.0*re112, s477'[t] == 1.0*re64 -1.0*re113, s478'[t] == 1.0*re67 -1.0*re114, s479'[t] == 1.0*re68 -1.0*re115, s482'[t] == 1.0*re33 -1.0*re72, s483'[t] == 1.0*re72 -1.0*re73, s484'[t] == 1.0*re73 -1.0*re74, s489'[t] == 1.0*re34 -1.0*re35, s490'[t] == 1.0*re75 -1.0*re34, s491'[t] == 1.0*re74 -1.0*re75, s492'[t] == 1.0*re35 -1.0*re36, s493'[t] == 1.0*re36 , s494'[t] == 1.0*re83 +1.0*re95 -1.0*re76 -1.0*re96, s495'[t] == 1.0*re84 +1.0*re96 -1.0*re77 -1.0*re97, s496'[t] == 1.0*re18 +1.0*re97 -1.0*re78 -1.0*re98, s498'[t] == 1.0*re85 +1.0*re98 -1.0*re79 -1.0*re99, s499'[t] == 1.0*re25 +1.0*re99 -1.0*re80 -1.0*re100, s5'[t] == 1.0*re49 +1.0*re121 -1.0*re50 -1.0*re51 -1.0*re52 -1.0*re123, s500'[t] == 1.0*re100 -1.0*re81 -1.0*re82 -1.0*re101, s501'[t] == 1.0*re30 +1.0*re101 -1.0*re71, s502'[t] == 1.0*re116 , s51'[t] == 1.0*re37 +1.0*re38 +1.0*re71 +1.0*re76 +1.0*re77 +1.0*re78 +1.0*re79 +1.0*re80 +1.0*re81 +1.0*re82 -1.0*re5 -1.0*re8 -1.0*re18 -1.0*re25 -1.0*re30 -1.0*re83 -1.0*re84 -1.0*re85 -1.0*re86, s517'[t] == 1.0*re12 -1.0*re40, s518'[t] == 1.0*re15 -1.0*re45, s519'[t] == 1.0*re19 -1.0*re46, s52'[t] == 1.0*re37 -1.0*re70, s520'[t] == 1.0*re9 -1.0*re42, s521'[t] == 1.0*re6 -1.0*re43, s522'[t] == 1.0*re22 -1.0*re44, s523'[t] == 1.0*re26 -1.0*re41, s524'[t] == 1.0*re91 +1.0*re92 , s525'[t] == 1.0*re65 -1.0*re93, s526'[t] == 1.0*re62 -1.0*re92, s527'[t] == 1.0*re48 -1.0*re89, s528'[t] == 1.0*re59 -1.0*re91, s529'[t] == 1.0*re56 -1.0*re90, s53'[t] == 1.0*re38 -1.0*re69, s530'[t] == 1.0*re53 -1.0*re87, s531'[t] == 1.0*re50 -1.0*re88, s533'[t] == 1.0*re134 -1.0*re135, s535'[t] == 1.0*re136 -1.0*re137, s536'[t] == 1.0*re86 -1.0*re33, s6'[t] == 1.0*re32 +1.0*re123 -1.0*re48 -1.0*re49 -1.0*re124, s7'[t] == 1.0*re4 +1.0*re117 +1.0*re118 +1.0*re119 +1.0*re120 +1.0*re121 +1.0*re122 +1.0*re123 +2.0*re124 -2.0*re32 -1.0*re33 -1.0*re34 -1.0*re35 -1.0*re36 -1.0*re49 -1.0*re51 -1.0*re54 -1.0*re57 -1.0*re60 -1.0*re63 -1.0*re66 -1.0*re72 -1.0*re73 -1.0*re74 -1.0*re75 -1.0*re86 -1.0*re136, s78'[t] == 1.0*re5 -1.0*re37 -1.0*re94, s85'[t] == 1.0*re8 +1.0*re94 -1.0*re38 -1.0*re95 }; 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]}]