(* 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 = { ADP[t], AMP[t], ATP[t], AcAld[t], BPG[t], DHAP[t], F16bP[t], F6P[t], G1P[t], G3P[t], G6P[t], GAP[t], GLC[t], NAD[t], NADH[t], P2G[t], P3G[t], PEP[t], PYR[t], T6P[t], UDG[t], UDP[t], UTP[t] }; initialValues = { ADP[0] == 1.29, AMP[0] == 0.44, ATP[0] == 4.29, AcAld[0] == 0.178140579850657, BPG[0] == 0.000736873499865602, DHAP[0] == 0.290344213186674, F16bP[0] == 1.14580464751733, F6P[0] == 0.0588603054728053, G1P[0] == 0.13481212658623, G3P[0] == 0.068500732297821, G6P[0] == 0.193120800911304, GAP[0] == 0.0789727571926259, GLC[0] == 1.5700004483456, NAD[0] == 1.50329030201531, NADH[0] == 0.0867096979846952, P2G[0] == 0.0169344770274836, P3G[0] == 0.117456252783611, PEP[0] == 0.152501353339511, PYR[0] == 0.527117851793548, T6P[0] == 0.02, UDG[0] == 0.467246194874247, UDP[0] == 0.2815, UTP[0] == 0.6491 }; rates = { ADHADH1, ADHADH5, AK, ATPase, ENOENO1, ENOENO2, FBA, GPD, GPM, GPP, HXKGLK1, HXKHXK2, HXT, PDCPDC1, PDCPDC5, PDCPDC6, PFK, PGI, PGK, PGM, PYK, TDHTDH1, TDHTDH2, TDHTDH3, TPI, TPP, TPS, UGP, acetatebranch, succinatebranch, udptoutp }; rateEquations = { ADHADH1 -> (ADH1*ADHADH1kcat*(-((EtOH*NAD[t])/(ADHADH1Kacald*ADHADH1Kinadh*KeqADH)) + (AcAld[t]*NADH[t])/(ADHADH1Kacald*ADHADH1Kinadh)))/(1 + (ADHADH1Knad*EtOH)/(ADHADH1Ketoh*ADHADH1Kinad) + (ADHADH1Knadh*AcAld[t])/(ADHADH1Kacald*ADHADH1Kinadh) + NAD[t]/ADHADH1Kinad + (EtOH*NAD[t])/(ADHADH1Ketoh*ADHADH1Kinad) + (ADHADH1Knadh*AcAld[t]*NAD[t])/(ADHADH1Kacald*ADHADH1Kinad*ADHADH1Kinadh) + (EtOH*AcAld[t]*NAD[t])/(ADHADH1Ketoh*ADHADH1Kiacald*ADHADH1Kinad) + NADH[t]/ADHADH1Kinadh + (ADHADH1Knad*EtOH*NADH[t])/(ADHADH1Ketoh*ADHADH1Kinad*ADHADH1Kinadh) + (AcAld[t]*NADH[t])/(ADHADH1Kacald*ADHADH1Kinadh) + (EtOH*AcAld[t]*NADH[t])/(ADHADH1Kacald*ADHADH1Kietoh*ADHADH1Kinadh)), ADHADH5 -> (ADH5*ADHADH5kcat*(-((EtOH*NAD[t])/(ADHADH5Kacald*ADHADH5Kinadh*KeqADH)) + (AcAld[t]*NADH[t])/(ADHADH5Kacald*ADHADH5Kinadh)))/(1 + (ADHADH5Knad*EtOH)/(ADHADH5Ketoh*ADHADH5Kinad) + (ADHADH5Knadh*AcAld[t])/(ADHADH5Kacald*ADHADH5Kinadh) + NAD[t]/ADHADH5Kinad + (EtOH*NAD[t])/(ADHADH5Ketoh*ADHADH5Kinad) + (ADHADH5Knadh*AcAld[t]*NAD[t])/(ADHADH5Kacald*ADHADH5Kinad*ADHADH5Kinadh) + (EtOH*AcAld[t]*NAD[t])/(ADHADH5Ketoh*ADHADH5Kiacald*ADHADH5Kinad) + NADH[t]/ADHADH5Kinadh + (ADHADH5Knad*EtOH*NADH[t])/(ADHADH5Ketoh*ADHADH5Kinad*ADHADH5Kinadh) + (AcAld[t]*NADH[t])/(ADHADH5Kacald*ADHADH5Kinadh) + (EtOH*AcAld[t]*NADH[t])/(ADHADH5Kacald*ADHADH5Kietoh*ADHADH5Kinadh)), AK -> AKk*(ADP[t]^2 - (AMP[t]*ATP[t])/AKKeq), ATPase -> (ATPaseVmax*ATP[t])/(ATPaseKatp*(1 + ATP[t]/ATPaseKatp)), ENOENO1 -> (ENO1*ENOENO1kcat*(P2G[t]/ENOENO1Kp2g - PEP[t]/(ENOENO1Kp2g*KeqENO)))/(1 + P2G[t]/ENOENO1Kp2g + PEP[t]/ENOENO1Kpep), ENOENO2 -> (ENO2*ENOENO2kcat*(P2G[t]/ENOENO2Kp2g - PEP[t]/(ENOENO2Kp2g*KeqENO)))/(1 + P2G[t]/ENOENO2Kp2g + PEP[t]/ENOENO2Kpep), FBA -> (FBA1*FBAkcat*(F16bP[t]/FBAKf16bp - (DHAP[t]*GAP[t])/(FBAKeq*FBAKf16bp)))/(1 + DHAP[t]/FBAKdhap + F16bP[t]/FBAKf16bp + GAP[t]/FBAKgap + (DHAP[t]*GAP[t])/(FBAKdhap*FBAKgap) + (F16bP[t]*GAP[t])/(FBAKf16bp*FBAKigap)), GPD -> (GPDVmax*(-((G3P[t]*NAD[t])/GPDKeq) + DHAP[t]*NADH[t]))/(GPDKdhap*GPDKnadh*(1 + ADP[t]/GPDKadp + ATP[t]/GPDKatp + F16bP[t]/GPDKfbp)*(1 + DHAP[t]/GPDKdhap + G3P[t]/GPDKg3p)*(1 + NAD[t]/GPDKnad + NADH[t]/GPDKnadh)), GPM -> (GPMVmax*(-(P2G[t]/(GPMKeq*GPMKp3g)) + P3G[t]/GPMKp3g))/(1 + P2G[t]/GPMKp2g + P3G[t]/GPMKp3g), GPP -> (GPPVmax*G3P[t])/(GPPKg3p*(1 + G3P[t]/GPPKg3p)), HXKGLK1 -> (GLK1*HXKGLK1kcat*(-((ADP[t]*G6P[t])/(HXKGLK1Katp*HXKGLK1Kglc*KeqHXK)) + (ATP[t]*GLC[t])/(HXKGLK1Katp*HXKGLK1Kglc)))/((1 + ADP[t]/HXKGLK1Kadp + ATP[t]/HXKGLK1Katp)*(1 + G6P[t]/HXKGLK1Kg6p + GLC[t]/HXKGLK1Kglc)), HXKHXK2 -> (HXK2*HXKHXK2kcat*(-((ADP[t]*G6P[t])/(HXKHXK2Katp*HXKHXK2Kglc*KeqHXK)) + (ATP[t]*GLC[t])/(HXKHXK2Katp*HXKHXK2Kglc)))/((1 + ADP[t]/HXKHXK2Kadp + ATP[t]/HXKHXK2Katp)*(1 + G6P[t]/HXKHXK2Kg6p + GLC[t]/HXKHXK2Kglc + T6P[t]/HXKHXK2Kit6p)), HXT -> (HXTVmax*(GLCx - GLC[t]))/(HXTKglc*(1 + GLCx/HXTKglc + GLC[t]/HXTKglc + (GLCx*HXTKi*GLC[t])/HXTKglc^2)), PDCPDC1 -> (PDC1*PDCPDC1kcat*PYR[t])/(PDCPDC1Kpyr*(1 + PYR[t]/PDCPDC1Kpyr)), PDCPDC5 -> (PDC5*PDCPDC5kcat*PYR[t])/(PDCPDC5Kpyr*(1 + PYR[t]/PDCPDC5Kpyr)), PDCPDC6 -> (PDC6*PDCPDC6kcat*PYR[t])/(PDCPDC6Kpyr*(1 + PYR[t]/PDCPDC6Kpyr)), PFK -> (PFKgR*PFKkcat*(PFK1 + PFK2 - Abs[PFK1 - PFK2])*ATP[t]*(1 - (ADP[t]*F16bP[t])/(PFKKeq*ATP[t]*F6P[t]))*F6P[t]*(1 + ADP[t]/PFKKadp + ATP[t]/PFKKatp + F16bP[t]/PFKKf16 + (PFKgR*ADP[t]*F16bP[t])/(PFKKadp*PFKKf16) + F6P[t]/PFKKf6p + (PFKgR*ATP[t]*F6P[t])/(PFKKatp*PFKKf6p)))/(2*PFKKatp*PFKKf6p*((PFKL0*(1 + (PFKCamp*AMP[t])/PFKKamp)^2*(1 + (PFKCatp*ATP[t])/PFKKatp)^2*(1 + (PFKCiatp*ATP[t])/PFKKiatp)^2*(1 + (F26bP*PFKCf26)/PFKKf26 + (PFKCf16*F16bP[t])/PFKKf16)^2)/((1 + AMP[t]/PFKKamp)^2*(1 + ATP[t]/PFKKiatp)^2*(1 + F26bP/PFKKf26 + F16bP[t]/PFKKf16)^2) + (1 + ADP[t]/PFKKadp + ATP[t]/PFKKatp + F16bP[t]/PFKKf16 + (PFKgR*ADP[t]*F16bP[t])/(PFKKadp*PFKKf16) + F6P[t]/PFKKf6p + (PFKgR*ATP[t]*F6P[t])/(PFKKatp*PFKKf6p))^2)), PGI -> (PGI1*PGIkcat*(-(F6P[t]/(PGIKeq*PGIKg6p)) + G6P[t]/PGIKg6p))/(1 + F6P[t]/PGIKf6p + G6P[t]/PGIKg6p), PGK -> (PGK1*PGKkcat*(ADP[t]/PGKKadp)^(-1 + PGKnHadp)*((ADP[t]*BPG[t])/(PGKKadp*PGKKbpg) - (ATP[t]*P3G[t])/(PGKKadp*PGKKbpg*PGKKeq)))/((1 + (ADP[t]/PGKKadp)^PGKnHadp + ATP[t]/PGKKatp)*(1 + BPG[t]/PGKKbpg + P3G[t]/PGKKp3g)), PGM -> (PGMVmax*rescaleTRH*(-(G1P[t]/(PGMKeq*PGMKg6p)) + G6P[t]/PGMKg6p))/(1 + G1P[t]/PGMKg1p + G6P[t]/PGMKg6p), PYK -> (PYKVmax*((ADP[t]*PEP[t])/(PYKKadp*PYKKpep) - (ATP[t]*PYR[t])/(PYKKadp*PYKKeq*PYKKpep)))/((1 + ADP[t]/PYKKadp + ATP[t]/PYKKatp)*(1 + PEP[t]/PYKKpep + PYR[t]/PYKKpyr)), TDHTDH1 -> (TDH1*TDHTDH1kcat*((GAP[t]*NAD[t])/(TDHTDH1Kgap*TDHTDH1Knad) - (BPG[t]*NADH[t])/(KeqTDH*TDHTDH1Kgap*TDHTDH1Knad)))/((1 + BPG[t]/TDHTDH1Kbpg + GAP[t]/TDHTDH1Kgap)*(1 + NAD[t]/TDHTDH1Knad + NADH[t]/TDHTDH1Knadh)), TDHTDH2 -> (TDH2*TDHTDH2kcat*((GAP[t]*NAD[t])/(TDHTDH2Kgap*TDHTDH2Knad) - (BPG[t]*NADH[t])/(KeqTDH*TDHTDH2Kgap*TDHTDH2Knad)))/((1 + BPG[t]/TDHTDH2Kbpg + GAP[t]/TDHTDH2Kgap)*(1 + NAD[t]/TDHTDH2Knad + NADH[t]/TDHTDH2Knadh)), TDHTDH3 -> (TDH3*TDHTDH3kcat*((GAP[t]*NAD[t])/(TDHTDH3Kgap*TDHTDH3Knad) - (BPG[t]*NADH[t])/(KeqTDH*TDHTDH3Kgap*TDHTDH3Knad)))/((1 + BPG[t]/TDHTDH3Kbpg + GAP[t]/TDHTDH3Kgap)*(1 + NAD[t]/TDHTDH3Knad + NADH[t]/TDHTDH3Knadh)), TPI -> TPIk*(DHAP[t] - GAP[t]/TPIKeq), TPP -> (rescaleTRH*TPPVmax*T6P[t])/(TPPKt6p*(1 + T6P[t]/TPPKt6p)), TPS -> (rescaleTRH*TPSVmax*G6P[t]*UDG[t])/(TPSKg6p*TPSKudg*(1 + G6P[t]/TPSKg6p)*(1 + UDG[t]/TPSKudg)), UGP -> (rescaleTRH*UGPVmax*G1P[t]*UTP[t])/(UGPKg1p*UGPKutp*(UGPKiutp/UGPKutp + G1P[t]/UGPKg1p + (UGPKiutp*UDG[t])/(UGPKiudg*UGPKutp) + (G1P[t]*UDG[t])/(UGPKg1p*UGPKiudg) + UTP[t]/UGPKutp + (G1P[t]*UTP[t])/(UGPKg1p*UGPKutp))), acetatebranch -> acetatebranchk*AcAld[t]*NAD[t], succinatebranch -> succinatebranchk*NAD[t]*PYR[t], udptoutp -> rescaleTRH*udptoutpk*ATP[t]*UDP[t] }; parameters = { ADH1 -> 0.0409771277320022, ADH5 -> 0.00106249605621922, ADHADH1Kacald -> 0.4622, ADHADH1Ketoh -> 17.0, ADHADH1Kiacald -> 1.1, ADHADH1Kietoh -> 90.0, ADHADH1Kinad -> 0.92, ADHADH1Kinadh -> 0.031, ADHADH1Knad -> 0.17, ADHADH1Knadh -> 0.11, ADHADH1kcat -> 176.0, ADHADH5Kacald -> 1.11, ADHADH5Ketoh -> 17.0, ADHADH5Kiacald -> 1.1, ADHADH5Kietoh -> 90.0, ADHADH5Kinad -> 0.92, ADHADH5Kinadh -> 0.031, ADHADH5Knad -> 0.17, ADHADH5Knadh -> 0.11, ADHADH5kcat -> 0.0, AKKeq -> 0.45, AKk -> 0.75, ATPaseKatp -> 3.0, ATPaseVmax -> 1.13333333333333, CDC19 -> 0.512097526792801, ENO1 -> 0.171592988538958, ENO2 -> 0.493611573294543, ENOENO1Kp2g -> 0.043, ENOENO1Kpep -> 0.5, ENOENO1kcat -> 7.6, ENOENO2Kp2g -> 0.104, ENOENO2Kpep -> 0.5, ENOENO2kcat -> 19.87, EXTERNAL -> 0.0, F26bP -> 0.003, FBA1 -> 0.334598664262206, FBAKdhap -> 2.0, FBAKeq -> 0.069, FBAKf16bp -> 0.4507, FBAKgap -> 2.4, FBAKigap -> 10.0, FBAkcat -> 4.139, GLK1 -> 0.0112717406104807, GPD1 -> 0.00170877794272468, GPD2 -> 0.000198351416606057, GPDKadp -> 2.0, GPDKatp -> 0.73, GPDKdhap -> 0.54, GPDKeq -> 10000.0, GPDKfbp -> 4.8, GPDKg3p -> 1.2, GPDKnad -> 0.93, GPDKnadh -> 0.023, GPDVmax -> 0.783333333333333, GPM1 -> 0.182500074724267, GPMKeq -> 0.19, GPMKp2g -> 0.08, GPMKp3g -> 1.2, GPMVmax -> 43.0833333333333, GPPKg3p -> 3.5, GPPVmax -> 0.883333333333333, HOR2 -> 0.0013683673909939, HXK1 -> 0.0041951864287446, HXK2 -> 0.0153328384926289, HXKGLK1Kadp -> 0.23, HXKGLK1Katp -> 0.865, HXKGLK1Kg6p -> 30.0, HXKGLK1Kglc -> 0.0106, HXKGLK1kcat -> 0.0721, HXKHXK1Kadp -> 0.23, HXKHXK1Katp -> 0.293, HXKHXK1Kg6p -> 30.0, HXKHXK1Kglc -> 0.15, HXKHXK1Kit6p -> 0.2, HXKHXK1kcat -> 10.2, HXKHXK2Kadp -> 0.23, HXKHXK2Katp -> 0.195, HXKHXK2Kg6p -> 30.0, HXKHXK2Kglc -> 0.2, HXKHXK2Kit6p -> 0.04, HXKHXK2kcat -> 63.1, HXTKglc -> 0.9, HXTKi -> 0.91, HXTVmax -> 3.35, KeqADH -> 14492.7536231884, KeqENO -> 6.7, KeqHXK -> 2000.0, KeqTDH -> 0.00533412710224736, NA -> 6.02214*^+20, PDC1 -> 0.266952694557084, PDC5 -> 0.00308868608169189, PDC6 -> 0.0016352160527653, PDCPDC1Kpyr -> 8.5, PDCPDC1kcat -> 12.14, PDCPDC5Kpyr -> 7.08, PDCPDC5kcat -> 10.32, PDCPDC6Kpyr -> 2.92, PDCPDC6kcat -> 9.21, PFK1 -> 0.0116962574765781, PFK2 -> 0.00975915538330228, PFKCamp -> 0.0845, PFKCatp -> 3.0, PFKCf16 -> 0.397, PFKCf26 -> 0.0174, PFKCiatp -> 100.0, PFKKadp -> 1.0, PFKKamp -> 0.0995, PFKKatp -> 0.71, PFKKeq -> 800.0, PFKKf16 -> 0.111, PFKKf26 -> 0.000682, PFKKf6p -> 0.1, PFKKiatp -> 0.65, PFKL0 -> 0.66, PFKgR -> 5.12, PFKkcat -> 209.6, PGI1 -> 0.0345726768225249, PGIKeq -> 0.29, PGIKf6p -> 0.307, PGIKg6p -> 1.0257, PGIkcat -> 487.36, PGK1 -> 0.0644142281647388, PGKKadp -> 0.2, PGKKatp -> 1.99, PGKKbpg -> 0.003, PGKKeq -> 3200.0, PGKKp3g -> 4.58, PGKkcat -> 58.6, PGKnHadp -> 2.0, PGM1 -> 0.000815573865768647, PGM2 -> 0.00031467219294138, PGMKeq -> 0.1667, PGMKg1p -> 0.023, PGMKg6p -> 0.05, PGMVmax -> 0.00590833333333333, PYK2 -> 0.00151748381804475, PYKKadp -> 0.53, PYKKatp -> 1.5, PYKKeq -> 6500.0, PYKKpep -> 0.14, PYKKpyr -> 21.0, PYKVmax -> 16.6666666666667, RHR2 -> 0.0127951193429578, TDH1 -> 0.087716160700349, TDH2 -> 0.0, TDH3 -> 1.05110118662137, TDHTDH1Kbpg -> 0.0098, TDHTDH1Kgap -> 0.495, TDHTDH1Knad -> 0.09, TDHTDH1Knadh -> 0.06, TDHTDH1kcat -> 19.12, TDHTDH2Kbpg -> 0.0098, TDHTDH2Kgap -> 0.77, TDHTDH2Knad -> 0.09, TDHTDH2Knadh -> 0.06, TDHTDH2kcat -> 8.633, TDHTDH3Kbpg -> 0.909, TDHTDH3Kgap -> 0.423, TDHTDH3Knad -> 0.09, TDHTDH3Knadh -> 0.06, TDHTDH3kcat -> 18.162, TPI1 -> 0.073589454911377, TPIKeq -> 0.045, TPIk -> 7500.0, TPPKt6p -> 0.5, TPPVmax -> 0.108333333333333, TPS1 -> 0.000848120435592663, TPS2 -> 0.000664962953368736, TPSKg6p -> 3.8, TPSKudg -> 0.886, TPSVmax -> 0.02285, UGP1 -> 0.00155052854965178, UGPKg1p -> 0.32, UGPKiudg -> 0.0035, UGPKiutp -> 0.11, UGPKutp -> 0.11, UGPVmax -> 0.613666666666667, acetatebranchk -> 0.00554339592436782, cell -> 1.0, extracellular -> 1.0, rescaleTRH -> 22.0, succinatebranchk -> 0.00549663630823988, sumAXP -> 6.02, sumNAD -> 1.59, sumUXP -> 1.39784619487425, udptoutpk -> 0.00345026988011002, volume -> 2*^-14, ACE -> 55.7500633495734, EtOH -> 55.472577854384, GLCx -> 74.0, GLY -> 0.15, SUC -> 0.0, TRH -> 0.00384697134241316, default\[LetterSpace]compartment -> 1.0 }; assignments = { energycharge -> (ADP[t]/2 + ATP[t])/sumAXP, sumPXG -> P2G[t] + P3G[t], fitconc -> Sqrt[(1 - (NA*sumPXG*volume)/1618640)^2 + (1 - (NA*volume*DHAP[t])/3496987)^2 + (1 - (NA*volume*F16bP[t])/13800392)^2 + (1 - (NA*volume*F6P[t])/708930)^2 + (1 - (NA*volume*G6P[t])/2326001)^2 + (1 - (NA*volume*GAP[t])/951170)^2 + (1 - (NA*volume*GLC[t])/18909525)^2 + (1 - (NA*volume*PEP[t])/1836769)^2 + (1 - (NA*volume*PYR[t])/6348755)^2]/3 }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { ADP'[t] == 1.0*HXKGLK1 +1.0*PFK +1.0*HXKHXK2 +1.0*ATPase +1.0*udptoutp -2.0*AK -1.0*PYK -1.0*PGK, AMP'[t] == 1.0*AK , ATP'[t] == 1.0*AK +1.0*PYK +1.0*PGK -1.0*HXKGLK1 -1.0*PFK -1.0*HXKHXK2 -1.0*ATPase -1.0*udptoutp, AcAld'[t] == 1.0*PDCPDC1 +1.0*PDCPDC5 +1.0*PDCPDC6 -1.0*ADHADH5 -1.0*ADHADH1 -1.0*acetatebranch, BPG'[t] == 1.0*TDHTDH1 +1.0*TDHTDH3 +1.0*TDHTDH2 -1.0*PGK, DHAP'[t] == 1.0*FBA -1.0*GPD -1.0*TPI, F16bP'[t] == 1.0*PFK -1.0*FBA, F6P'[t] == 1.0*PGI -1.0*PFK, G1P'[t] == 1.0*PGM -1.0*UGP, G3P'[t] == 1.0*GPD -1.0*GPP, G6P'[t] == 1.0*HXKGLK1 +1.0*HXKHXK2 -1.0*TPS -1.0*PGM -1.0*PGI, GAP'[t] == 1.0*FBA +1.0*TPI -1.0*TDHTDH1 -1.0*TDHTDH3 -1.0*TDHTDH2, GLC'[t] == 1.0*HXT -1.0*HXKGLK1 -1.0*HXKHXK2, NAD'[t] == 1.0*ADHADH5 +1.0*ADHADH1 +1.0*GPD -1.0*TDHTDH1 -3.0*succinatebranch -1.0*acetatebranch -1.0*TDHTDH3 -1.0*TDHTDH2, NADH'[t] == 1.0*TDHTDH1 +3.0*succinatebranch +1.0*acetatebranch +1.0*TDHTDH3 +1.0*TDHTDH2 -1.0*GPD -1.0*ADHADH5 -1.0*ADHADH1, P2G'[t] == 1.0*GPM -1.0*ENOENO2 -1.0*ENOENO1, P3G'[t] == 1.0*PGK -1.0*GPM, PEP'[t] == 1.0*ENOENO2 +1.0*ENOENO1 -1.0*PYK, PYR'[t] == 1.0*PYK -1.0*PDCPDC1 -1.0*PDCPDC5 -1.0*PDCPDC6 -1.0*succinatebranch, T6P'[t] == 1.0*TPS -1.0*TPP, UDG'[t] == 1.0*UGP -1.0*TPS, UDP'[t] == 1.0*TPS -1.0*udptoutp, UTP'[t] == 1.0*udptoutp -1.0*UGP }; 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]}]