(* 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 = { amp[t], ampE[t], arg[t], argE[t], argQ[t], argQE[t], argR[t], argRE[t], as[t], asE[t], asR[t], asRE[t], asp[t], aspE[t], atp[t], atpE[t], cit[t], citE[t], citR[t], citRE[t], citims[t], citimsR[t], citmat[t], citmatR[t], cp[t], cpR[t], cpe[t], cpeE[t], cpeR[t], cpeRE[t], cpims[t], cpimsR[t], fum[t], fumE[t], orn[t], ornE[t], ornims[t], ornmat[t], pi[t], ppi[t], ppiE[t], ure[t], ureE[t], ureR[t], ureRE[t] }; initialValues = { amp[0] == 1*^-06, ampE[0] == 1*^-06, arg[0] == 1*^-06, argE[0] == 1*^-06, argQ[0] == 0.0, argQE[0] == 0.0, argR[0] == 0.0, argRE[0] == 0.0, as[0] == 1*^-06, asE[0] == 1*^-06, asR[0] == 0.0, asRE[0] == 0.0, asp[0] == 0.0001, aspE[0] == 0.0001, atp[0] == 0.0015, atpE[0] == 0.0015, cit[0] == 1*^-06, citE[0] == 1*^-06, citR[0] == 0.0, citRE[0] == 0.0, citims[0] == 1*^-06, citimsR[0] == 0.0, citmat[0] == 1*^-06, citmatR[0] == 0.0, cp[0] == 1*^-05, cpR[0] == 0.0, cpe[0] == 1*^-05, cpeE[0] == 0.0005, cpeR[0] == 0.0, cpeRE[0] == 2.26667*^-06, cpims[0] == 1*^-05, cpimsR[0] == 0.0, fum[0] == 1*^-06, fumE[0] == 1*^-06, orn[0] == 0.0002, ornE[0] == 0.0002, ornims[0] == 0.00045, ornmat[0] == 1*^-06, pi[0] == 0.0001, ppi[0] == 1*^-06, ppiE[0] == 1*^-06, ure[0] == 1*^-06, ureE[0] == 1*^-06, ureR[0] == 0.0, ureRE[0] == 0.0 }; rates = { v\[LetterSpace]1, v\[LetterSpace]10, v\[LetterSpace]11, v\[LetterSpace]12, v\[LetterSpace]13, v\[LetterSpace]14, v\[LetterSpace]15, v\[LetterSpace]16, v\[LetterSpace]17, v\[LetterSpace]18, v\[LetterSpace]19, v\[LetterSpace]2, v\[LetterSpace]20, v\[LetterSpace]21, v\[LetterSpace]22, v\[LetterSpace]23, v\[LetterSpace]24, v\[LetterSpace]25, v\[LetterSpace]26, v\[LetterSpace]27, v\[LetterSpace]28, v\[LetterSpace]29, v\[LetterSpace]3, v\[LetterSpace]30, v\[LetterSpace]31, v\[LetterSpace]32, v\[LetterSpace]33, v\[LetterSpace]34, v\[LetterSpace]35, v\[LetterSpace]36, v\[LetterSpace]37, v\[LetterSpace]38, v\[LetterSpace]39, v\[LetterSpace]4, v\[LetterSpace]40, v\[LetterSpace]41, v\[LetterSpace]42, v\[LetterSpace]43, v\[LetterSpace]44, v\[LetterSpace]45, v\[LetterSpace]46, v\[LetterSpace]47, v\[LetterSpace]48, v\[LetterSpace]5, v\[LetterSpace]6, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9 }; rateEquations = { v\[LetterSpace]1 -> (e10*k1e1*k3e1*k5e1*k7e1*matV*otcAdj*cp[t]*ornmat[t])/(k2e1*(k4e1 + k5e1)*k7e1 + k2e1*k4e1*k6e1*(citmat[t] + citmatR[t]) + k1e1*(k4e1 + k5e1)*k7e1*(cp[t] + cpR[t]) + k1e1*k4e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t]) + k3e1*k5e1*k7e1*ornmat[t] + k1e1*k3e1*(k5e1 + k7e1)*(cp[t] + cpR[t])*ornmat[t] + k1e1*k3e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t])*ornmat[t] + k2e1*(k4e1 + k5e1)*k8e1*pi[t] + (k2e1 + k4e1)*k6e1*k8e1*(citmat[t] + citmatR[t])*pi[t] + k3e1*k5e1*k8e1*ornmat[t]*pi[t] + k3e1*k6e1*k8e1*(citmat[t] + citmatR[t])*ornmat[t]*pi[t]) - (e10*k2e1*k4e1*k6e1*k8e1*matV*otcAdj*citmat[t]*pi[t])/(k2e1*(k4e1 + k5e1)*k7e1 + k2e1*k4e1*k6e1*(citmat[t] + citmatR[t]) + k1e1*(k4e1 + k5e1)*k7e1*(cp[t] + cpR[t]) + k1e1*k4e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t]) + k3e1*k5e1*k7e1*ornmat[t] + k1e1*k3e1*(k5e1 + k7e1)*(cp[t] + cpR[t])*ornmat[t] + k1e1*k3e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t])*ornmat[t] + k2e1*(k4e1 + k5e1)*k8e1*pi[t] + (k2e1 + k4e1)*k6e1*k8e1*(citmat[t] + citmatR[t])*pi[t] + k3e1*k5e1*k8e1*ornmat[t]*pi[t] + k3e1*k6e1*k8e1*(citmat[t] + citmatR[t])*ornmat[t]*pi[t]), v\[LetterSpace]10 -> (aslAdj*cytV*e30*k1e3*k3e3*k5e3*asR[t])/((k2e3 + k3e3)*k5e3 + (k2e3 + k3e3)*k6e3*(arg[t] + fm*argQ[t] + argR[t]) + k1e3*(k3e3 + k5e3)*(as[t] + asR[t]) + k2e3*k4e3*fum[t] + k4e3*k6e3*(arg[t] + fm*argQ[t] + argR[t])*fum[t] + k1e3*k4e3*(fm*argQ[t] + as[t] + asR[t])*fum[t]) - (aslAdj*cytV*e30*k2e3*k4e3*k6e3*argR[t]*fum[t])/((k2e3 + k3e3)*k5e3 + (k2e3 + k3e3)*k6e3*(arg[t] + fm*argQ[t] + argR[t]) + k1e3*(k3e3 + k5e3)*(as[t] + asR[t]) + k2e3*k4e3*fum[t] + k4e3*k6e3*(arg[t] + fm*argQ[t] + argR[t])*fum[t] + k1e3*k4e3*(fm*argQ[t] + as[t] + asR[t])*fum[t]), v\[LetterSpace]11 -> (argfac*cytV*e40*k1e4*k3e4*k4e4*arg[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(arg[t] + fm*argQ[t] + argR[t]) + (k2e4 + k3e4)*k5e4*orn[t]), v\[LetterSpace]12 -> (argfac*cytV*e40*k1e4*k3e4*k4e4*argR[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(arg[t] + fm*argQ[t] + argR[t]) + (k2e4 + k3e4)*k5e4*orn[t]), v\[LetterSpace]13 -> (argfac*cytV*e40*fm*k1e4*k3e4*k4e4*argQ[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(arg[t] + fm*argQ[t] + argR[t]) + (k2e4 + k3e4)*k5e4*orn[t]), v\[LetterSpace]14 -> (argfac*e40E*eV*k1e4*k3e4*k4e4*argE[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(argE[t] + fmE*argQE[t] + argRE[t]) + (k2e4 + k3e4)*k5e4*ornE[t]), v\[LetterSpace]15 -> (argfac*e40E*eV*k1e4*k3e4*k4e4*argRE[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(argE[t] + fmE*argQE[t] + argRE[t]) + (k2e4 + k3e4)*k5e4*ornE[t]), v\[LetterSpace]16 -> (argfac*e40E*eV*fm*k1e4*k3e4*k4e4*argQE[t])/((k2e4 + k3e4)*k4e4 + k1e4*(k3e4 + k4e4)*(argE[t] + fmE*argQE[t] + argRE[t]) + (k2e4 + k3e4)*k5e4*ornE[t]), v\[LetterSpace]17 -> cytV*kornci*orn[t] - imsV*kornic*ornims[t], v\[LetterSpace]18 -> imsV*kornim*ornims[t] - kornmi*matV*ornmat[t], v\[LetterSpace]19 -> cppool*eV*kcp, v\[LetterSpace]2 -> (e10*k1e1*k3e1*k5e1*k7e1*matV*otcAdj*cpR[t]*ornmat[t])/(k2e1*(k4e1 + k5e1)*k7e1 + k2e1*k4e1*k6e1*(citmat[t] + citmatR[t]) + k1e1*(k4e1 + k5e1)*k7e1*(cp[t] + cpR[t]) + k1e1*k4e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t]) + k3e1*k5e1*k7e1*ornmat[t] + k1e1*k3e1*(k5e1 + k7e1)*(cp[t] + cpR[t])*ornmat[t] + k1e1*k3e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t])*ornmat[t] + k2e1*(k4e1 + k5e1)*k8e1*pi[t] + (k2e1 + k4e1)*k6e1*k8e1*(citmat[t] + citmatR[t])*pi[t] + k3e1*k5e1*k8e1*ornmat[t]*pi[t] + k3e1*k6e1*k8e1*(citmat[t] + citmatR[t])*ornmat[t]*pi[t]) - (e10*k2e1*k4e1*k6e1*k8e1*matV*otcAdj*citmatR[t]*pi[t])/(k2e1*(k4e1 + k5e1)*k7e1 + k2e1*k4e1*k6e1*(citmat[t] + citmatR[t]) + k1e1*(k4e1 + k5e1)*k7e1*(cp[t] + cpR[t]) + k1e1*k4e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t]) + k3e1*k5e1*k7e1*ornmat[t] + k1e1*k3e1*(k5e1 + k7e1)*(cp[t] + cpR[t])*ornmat[t] + k1e1*k3e1*k6e1*(citmat[t] + citmatR[t])*(cp[t] + cpR[t])*ornmat[t] + k2e1*(k4e1 + k5e1)*k8e1*pi[t] + (k2e1 + k4e1)*k6e1*k8e1*(citmat[t] + citmatR[t])*pi[t] + k3e1*k5e1*k8e1*ornmat[t]*pi[t] + k3e1*k6e1*k8e1*(citmat[t] + citmatR[t])*ornmat[t]*pi[t]), v\[LetterSpace]20 -> cytV*kcpei*cpe[t] - imsV*kcpie*cpims[t], v\[LetterSpace]21 -> -(kcpmi*matV*cp[t]) + imsV*kcpim*cpims[t], v\[LetterSpace]22 -> cytV*kcpei*cpeR[t] - imsV*kcpie*cpimsR[t], v\[LetterSpace]23 -> imsV*kcpim*cpimsR[t] - kcpmi*matV*cpR[t], v\[LetterSpace]24 -> atppool*cytV*katp, v\[LetterSpace]25 -> asppool*cytV*kasp, v\[LetterSpace]26 -> cytV*kppi*ppi[t], v\[LetterSpace]27 -> cytV*kfum*fum[t], v\[LetterSpace]28 -> kpi*matV*pi[t], v\[LetterSpace]29 -> cytV*kamp*amp[t], v\[LetterSpace]3 -> -(imsV*kcitim*citims[t]) + kcitmi*matV*citmat[t], v\[LetterSpace]30 -> cytV*kure*ure[t], v\[LetterSpace]31 -> 0, v\[LetterSpace]32 -> cytV*kcitCE*cit[t] - eV*kcitCEREV*citE[t], v\[LetterSpace]33 -> cytV*kcitRCE*citR[t] - eV*kcitRCEREV*citRE[t], v\[LetterSpace]34 -> cytV*katpCE*atp[t] - eV*katpCEREV*atpE[t], v\[LetterSpace]35 -> cytV*kaspCE*asp[t] - eV*kaspCEREV*aspE[t], v\[LetterSpace]36 -> cytV*kampCE*amp[t] - eV*kampCEREV*ampE[t], v\[LetterSpace]37 -> cytV*kppiCE*ppi[t] - eV*kppiCEREV*ppiE[t], v\[LetterSpace]38 -> cytV*kasCE*as[t] - eV*kasCEREV*asE[t], v\[LetterSpace]39 -> cytV*kasRCE*asR[t] - eV*kasRCEREV*asRE[t], v\[LetterSpace]4 -> -(imsV*kcitim*citimsR[t]) + kcitmi*matV*citmatR[t], v\[LetterSpace]40 -> cytV*kfumCE*fum[t] - eV*kfumCEREV*fumE[t], v\[LetterSpace]41 -> cytV*kargQCE*argQ[t] - eV*kargQCEREV*argQE[t], v\[LetterSpace]42 -> cytV*kargCE*arg[t] - eV*kargCEREV*argE[t], v\[LetterSpace]43 -> cytV*kargRCE*argR[t] - eV*kargRCEREV*argRE[t], v\[LetterSpace]44 -> cytV*kureCE*ure[t] - eV*kureCEREV*ureE[t], v\[LetterSpace]45 -> cytV*kureRCE*ureR[t] - eV*kureRCEREV*ureRE[t], v\[LetterSpace]46 -> cytV*kornCE*orn[t] - eV*kornCEREV*ornE[t], v\[LetterSpace]47 -> cytV*kcpeCE*cpe[t] - eV*kcpeCEREV*cpeE[t], v\[LetterSpace]48 -> cytV*kcpeRCE*cpeR[t] - eV*kcpeRCEREV*cpeRE[t], v\[LetterSpace]5 -> -(cytV*kcitci*cit[t]) + imsV*kcitic*citims[t], v\[LetterSpace]6 -> imsV*kcitic*citimsR[t] - cytV*kcitci*citR[t], v\[LetterSpace]7 -> (assAdj*cytV*e20*k11e2*k1e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t]*cit[t])/((1 + (arg[t] + argQ[t] + argR[t])/Kiarg)*(k11e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2 + k11e2*k2e2*k5e2*k7e2*k9e2*asp[t] + k12e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k4e2*(k6e2 + k7e2)*amp[t]*(as[t] + asR[t]) + k12e2*k2e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t]) + k11e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t] + k12e2*k3e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t])*atp[t] + k10e2*k12e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t] + k11e2*k1e2*k4e2*(k6e2 + k7e2)*k9e2*(cit[t] + citR[t]) + k11e2*k1e2*k5e2*k7e2*k9e2*asp[t]*(cit[t] + citR[t]) + k11e2*k1e2*k3e2*(k6e2 + k7e2)*k9e2*atp[t]*(cit[t] + citR[t]) + k1e2*k3e2*k5e2*(k11e2*k7e2 + k11e2*k9e2 + k7e2*k9e2)*asp[t]*atp[t]*(cit[t] + citR[t]) + k10e2*k1e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t]) + k11e2*k2e2*k4e2*k6e2*k8e2*ppi[t] + k10e2*k2e2*k4e2*k6e2*k8e2*amp[t]*ppi[t] + k12e2*k2e2*k4e2*k6e2*k8e2*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*(k2e2*k4e2 + k2e2*k6e2 + k4e2*k6e2)*k8e2*amp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k2e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k3e2*k6e2*k8e2*amp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k10e2*k12e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k11e2*k1e2*k4e2*k6e2*k8e2*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k4e2*k6e2*k8e2*amp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k6e2*k8e2*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k6e2*k8e2*amp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k5e2*k8e2*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t])) - (assAdj*cytV*e20*k10e2*k12e2*k2e2*k4e2*k6e2*k8e2*amp[t]*as[t]*ppi[t])/((1 + (arg[t] + argQ[t] + argR[t])/Kiarg)*(k11e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2 + k11e2*k2e2*k5e2*k7e2*k9e2*asp[t] + k12e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k4e2*(k6e2 + k7e2)*amp[t]*(as[t] + asR[t]) + k12e2*k2e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t]) + k11e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t] + k12e2*k3e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t])*atp[t] + k10e2*k12e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t] + k11e2*k1e2*k4e2*(k6e2 + k7e2)*k9e2*(cit[t] + citR[t]) + k11e2*k1e2*k5e2*k7e2*k9e2*asp[t]*(cit[t] + citR[t]) + k11e2*k1e2*k3e2*(k6e2 + k7e2)*k9e2*atp[t]*(cit[t] + citR[t]) + k1e2*k3e2*k5e2*(k11e2*k7e2 + k11e2*k9e2 + k7e2*k9e2)*asp[t]*atp[t]*(cit[t] + citR[t]) + k10e2*k1e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t]) + k11e2*k2e2*k4e2*k6e2*k8e2*ppi[t] + k10e2*k2e2*k4e2*k6e2*k8e2*amp[t]*ppi[t] + k12e2*k2e2*k4e2*k6e2*k8e2*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*(k2e2*k4e2 + k2e2*k6e2 + k4e2*k6e2)*k8e2*amp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k2e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k3e2*k6e2*k8e2*amp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k10e2*k12e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k11e2*k1e2*k4e2*k6e2*k8e2*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k4e2*k6e2*k8e2*amp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k6e2*k8e2*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k6e2*k8e2*amp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k5e2*k8e2*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t])), v\[LetterSpace]8 -> (assAdj*cytV*e20*k11e2*k1e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t]*citR[t])/((1 + (arg[t] + argQ[t] + argR[t])/Kiarg)*(k11e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2 + k11e2*k2e2*k5e2*k7e2*k9e2*asp[t] + k12e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k4e2*(k6e2 + k7e2)*amp[t]*(as[t] + asR[t]) + k12e2*k2e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t]) + k11e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t] + k12e2*k3e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t])*atp[t] + k10e2*k12e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t] + k11e2*k1e2*k4e2*(k6e2 + k7e2)*k9e2*(cit[t] + citR[t]) + k11e2*k1e2*k5e2*k7e2*k9e2*asp[t]*(cit[t] + citR[t]) + k11e2*k1e2*k3e2*(k6e2 + k7e2)*k9e2*atp[t]*(cit[t] + citR[t]) + k1e2*k3e2*k5e2*(k11e2*k7e2 + k11e2*k9e2 + k7e2*k9e2)*asp[t]*atp[t]*(cit[t] + citR[t]) + k10e2*k1e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t]) + k11e2*k2e2*k4e2*k6e2*k8e2*ppi[t] + k10e2*k2e2*k4e2*k6e2*k8e2*amp[t]*ppi[t] + k12e2*k2e2*k4e2*k6e2*k8e2*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*(k2e2*k4e2 + k2e2*k6e2 + k4e2*k6e2)*k8e2*amp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k2e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k3e2*k6e2*k8e2*amp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k10e2*k12e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k11e2*k1e2*k4e2*k6e2*k8e2*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k4e2*k6e2*k8e2*amp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k6e2*k8e2*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k6e2*k8e2*amp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k5e2*k8e2*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t])) - (assAdj*cytV*e20*k10e2*k12e2*k2e2*k4e2*k6e2*k8e2*amp[t]*asR[t]*ppi[t])/((1 + (arg[t] + argQ[t] + argR[t])/Kiarg)*(k11e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2 + k11e2*k2e2*k5e2*k7e2*k9e2*asp[t] + k12e2*k2e2*k4e2*(k6e2 + k7e2)*k9e2*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k4e2*(k6e2 + k7e2)*amp[t]*(as[t] + asR[t]) + k12e2*k2e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t]) + k10e2*k12e2*k2e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t]) + k11e2*k3e2*k5e2*k7e2*k9e2*asp[t]*atp[t] + k12e2*k3e2*k5e2*k7e2*k9e2*asp[t]*(as[t] + asR[t])*atp[t] + k10e2*k12e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t] + k11e2*k1e2*k4e2*(k6e2 + k7e2)*k9e2*(cit[t] + citR[t]) + k11e2*k1e2*k5e2*k7e2*k9e2*asp[t]*(cit[t] + citR[t]) + k11e2*k1e2*k3e2*(k6e2 + k7e2)*k9e2*atp[t]*(cit[t] + citR[t]) + k1e2*k3e2*k5e2*(k11e2*k7e2 + k11e2*k9e2 + k7e2*k9e2)*asp[t]*atp[t]*(cit[t] + citR[t]) + k10e2*k1e2*k3e2*k5e2*k7e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t]) + k11e2*k2e2*k4e2*k6e2*k8e2*ppi[t] + k10e2*k2e2*k4e2*k6e2*k8e2*amp[t]*ppi[t] + k12e2*k2e2*k4e2*k6e2*k8e2*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*(k2e2*k4e2 + k2e2*k6e2 + k4e2*k6e2)*k8e2*amp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k2e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*ppi[t] + k10e2*k12e2*k3e2*k6e2*k8e2*amp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k10e2*k12e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*(as[t] + asR[t])*atp[t]*ppi[t] + k11e2*k1e2*k4e2*k6e2*k8e2*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k4e2*k6e2*k8e2*amp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k6e2*k8e2*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k6e2*k8e2*amp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k11e2*k1e2*k3e2*k5e2*k8e2*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t] + k10e2*k1e2*k3e2*k5e2*k8e2*amp[t]*asp[t]*atp[t]*(cit[t] + citR[t])*ppi[t])), v\[LetterSpace]9 -> (aslAdj*cytV*e30*k1e3*k3e3*k5e3*as[t])/((k2e3 + k3e3)*k5e3 + (k2e3 + k3e3)*k6e3*(arg[t] + fm*argQ[t] + argR[t]) + k1e3*(k3e3 + k5e3)*(as[t] + asR[t]) + k2e3*k4e3*fum[t] + k4e3*k6e3*(arg[t] + fm*argQ[t] + argR[t])*fum[t] + k1e3*k4e3*(fm*argQ[t] + as[t] + asR[t])*fum[t]) - (aslAdj*cytV*e30*k2e3*k4e3*k6e3*arg[t]*fum[t])/((k2e3 + k3e3)*k5e3 + (k2e3 + k3e3)*k6e3*(arg[t] + fm*argQ[t] + argR[t]) + k1e3*(k3e3 + k5e3)*(as[t] + asR[t]) + k2e3*k4e3*fum[t] + k4e3*k6e3*(arg[t] + fm*argQ[t] + argR[t])*fum[t] + k1e3*k4e3*(fm*argQ[t] + as[t] + asR[t])*fum[t]) }; parameters = { Kiarg -> 0.0007, ampfinal -> 5*^-05, argfac -> 1.0, aslAdj -> 1.0, assAdj -> 2.0, cellV -> 3.868*^-05, cytV -> 2.47552*^-05, e10 -> 1.49288*^-05, e20 -> 6.25*^-06, e30 -> 3.4375*^-06, e40 -> 1.20984*^-05, e40E -> 2.28177*^-08, eV -> 0.00196132, fm -> 1.0, fmE -> 1.0, fumfinal -> 1*^-05, imsV -> 4.09389*^-06, k10e2 -> 640000.0, k11e2 -> 50.0, k12e2 -> 170000.0, k1e1 -> 17000000.0, k1e2 -> 240000.0, k1e3 -> 2700000.0, k1e4 -> 10000000.0, k2e1 -> 63.0, k2e2 -> 2.3, k2e3 -> 70.0, k2e4 -> 54000.0, k3e1 -> 2100000.0, k3e2 -> 350000.0, k3e3 -> 75.0, k3e4 -> 5300.0, k4e1 -> 1000.0, k4e2 -> 10.0, k4e3 -> 1500000.0, k4e4 -> 30000.0, k5e1 -> 3000.0, k5e2 -> 480000.0, k5e3 -> 1100.0, k5e4 -> 10000000.0, k6e1 -> 90000.0, k6e2 -> 10.0, k6e3 -> 700000.0, k7e1 -> 2600.0, k7e2 -> 20.0, k8e1 -> 500000.0, k8e2 -> 890000.0, k9e2 -> 50.0, kamp -> 0.00630208, kampCE -> 10000000.0, kampCEREV -> 126217.0, kargCE -> 10000000.0, kargCEREV -> 126217.0, kargQCE -> 10000000.0, kargQCEREV -> 126217.0, kargRCE -> 10000000.0, kargRCEREV -> 126217.0, kasCE -> 10000000.0, kasCEREV -> 126217.0, kasRCE -> 10000000.0, kasRCEREV -> 126217.0, kasp -> 0.000315104, kaspCE -> 10000000.0, kaspCEREV -> 126217.0, katp -> 0.000157552, katpCE -> 10000000.0, katpCEREV -> 126217.0, kcitCE -> 10000000.0, kcitCEREV -> 126217.0, kcitRCE -> 10000000.0, kcitRCEREV -> 126217.0, kcitci -> 10000.0, kcitic -> 60468.6, kcitim -> 100.0, kcitmi -> 60.7717, kcp -> 2.65143*^-07, kcpeCE -> 10000000.0, kcpeCEREV -> 126217.0, kcpeRCE -> 10000000.0, kcpeRCEREV -> 126217.0, kcpei -> 10000.0, kcpie -> 60468.6, kcpim -> 0.822751, kcpmi -> 0.5, kfum -> 0.0315104, kfumCE -> 10000000.0, kfumCEREV -> 126217.0, kornCE -> 10000000.0, kornCEREV -> 126217.0, kornci -> 10000.0, kornic -> 60468.6, kornim -> 110.0, kornmi -> 66.8489, kpi -> 0.0231588, kppi -> 0.0315104, kppiCE -> 10000000.0, kppiCEREV -> 126217.0, kure -> 0.315104, kureCE -> 10000000.0, kureCEREV -> 126217.0, kureRCE -> 10000000.0, kureRCEREV -> 126217.0, matV -> 6.73651*^-06, mitoV -> 1.08304*^-05, nucV -> 3.0944*^-06, otcAdj -> 1.0, pifinal -> 5*^-05, ppifinal -> 1*^-05, sampleV -> 0.002, steadyStateUreaProdn -> 7.80047*^-12, urefinal -> 1*^-06, XX -> 0.0, asppool -> 0.001, atppool -> 0.002, cppool -> 0.015, default -> 1.0 }; assignments = { }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { amp'[t] == 40395.6*v\[LetterSpace]7 +40395.6*v\[LetterSpace]8 -40395.6*v\[LetterSpace]29 -40395.6*v\[LetterSpace]36, ampE'[t] == 509.861*v\[LetterSpace]36 , arg'[t] == 40395.6*v\[LetterSpace]9 -40395.6*v\[LetterSpace]42 -40395.6*v\[LetterSpace]11, argE'[t] == 509.861*v\[LetterSpace]42 -509.861*v\[LetterSpace]14, argQ'[t] == -40395.6*v\[LetterSpace]41 -40395.6*v\[LetterSpace]13, argQE'[t] == 509.861*v\[LetterSpace]41 -509.861*v\[LetterSpace]16, argR'[t] == 40395.6*v\[LetterSpace]10 -40395.6*v\[LetterSpace]43 -40395.6*v\[LetterSpace]12, argRE'[t] == 509.861*v\[LetterSpace]43 -509.861*v\[LetterSpace]15, as'[t] == 40395.6*v\[LetterSpace]7 -40395.6*v\[LetterSpace]38 -40395.6*v\[LetterSpace]9, asE'[t] == 509.861*v\[LetterSpace]38 , asR'[t] == 40395.6*v\[LetterSpace]8 -40395.6*v\[LetterSpace]39 -40395.6*v\[LetterSpace]10, asRE'[t] == 509.861*v\[LetterSpace]39 , asp'[t] == 40395.6*v\[LetterSpace]25 -40395.6*v\[LetterSpace]35 -40395.6*v\[LetterSpace]8 -40395.6*v\[LetterSpace]7, aspE'[t] == 509.861*v\[LetterSpace]35 , atp'[t] == 40395.6*v\[LetterSpace]24 -40395.6*v\[LetterSpace]34 -40395.6*v\[LetterSpace]8 -40395.6*v\[LetterSpace]7, atpE'[t] == 509.861*v\[LetterSpace]34 , cit'[t] == 40395.6*v\[LetterSpace]5 -40395.6*v\[LetterSpace]32 -40395.6*v\[LetterSpace]7, citE'[t] == 509.861*v\[LetterSpace]32 , citR'[t] == 40395.6*v\[LetterSpace]6 -40395.6*v\[LetterSpace]33 -40395.6*v\[LetterSpace]8, citRE'[t] == 509.861*v\[LetterSpace]33 , citims'[t] == 244266.0*v\[LetterSpace]3 -244266.0*v\[LetterSpace]5, citimsR'[t] == 244266.0*v\[LetterSpace]4 -244266.0*v\[LetterSpace]6, citmat'[t] == 148445.0*v\[LetterSpace]1 -148445.0*v\[LetterSpace]3, citmatR'[t] == 148445.0*v\[LetterSpace]2 -148445.0*v\[LetterSpace]4, cp'[t] == 148445.0*v\[LetterSpace]21 -148445.0*v\[LetterSpace]1, cpR'[t] == 148445.0*v\[LetterSpace]23 -148445.0*v\[LetterSpace]2, cpe'[t] == -40395.6*v\[LetterSpace]47 -40395.6*v\[LetterSpace]20, cpeE'[t] == 509.861*v\[LetterSpace]47 +509.861*v\[LetterSpace]19 , cpeR'[t] == -40395.6*v\[LetterSpace]48 -40395.6*v\[LetterSpace]22, cpeRE'[t] == 509.861*v\[LetterSpace]48 , cpims'[t] == 244266.0*v\[LetterSpace]20 -244266.0*v\[LetterSpace]21, cpimsR'[t] == 244266.0*v\[LetterSpace]22 -244266.0*v\[LetterSpace]23, fum'[t] == 40395.6*v\[LetterSpace]9 +40395.6*v\[LetterSpace]10 -40395.6*v\[LetterSpace]27 -40395.6*v\[LetterSpace]40, fumE'[t] == 509.861*v\[LetterSpace]40 , orn'[t] == 40395.6*v\[LetterSpace]11 +40395.6*v\[LetterSpace]12 +40395.6*v\[LetterSpace]13 -40395.6*v\[LetterSpace]17 -40395.6*v\[LetterSpace]46, ornE'[t] == 509.861*v\[LetterSpace]14 +509.861*v\[LetterSpace]46 +509.861*v\[LetterSpace]16 +509.861*v\[LetterSpace]15 , ornims'[t] == 244266.0*v\[LetterSpace]17 -244266.0*v\[LetterSpace]18, ornmat'[t] == 148445.0*v\[LetterSpace]18 -148445.0*v\[LetterSpace]1 -148445.0*v\[LetterSpace]2, pi'[t] == 148445.0*v\[LetterSpace]1 +148445.0*v\[LetterSpace]2 -148445.0*v\[LetterSpace]28, ppi'[t] == 40395.6*v\[LetterSpace]7 +40395.6*v\[LetterSpace]8 -40395.6*v\[LetterSpace]37 -40395.6*v\[LetterSpace]26, ppiE'[t] == 509.861*v\[LetterSpace]37 , ure'[t] == 40395.6*v\[LetterSpace]11 +40395.6*v\[LetterSpace]13 -40395.6*v\[LetterSpace]44 -40395.6*v\[LetterSpace]30, ureE'[t] == 509.861*v\[LetterSpace]44 +509.861*v\[LetterSpace]14 +509.861*v\[LetterSpace]16 , ureR'[t] == 40395.6*v\[LetterSpace]12 -40395.6*v\[LetterSpace]31 -40395.6*v\[LetterSpace]45, ureRE'[t] == 509.861*v\[LetterSpace]45 +509.861*v\[LetterSpace]15 }; 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]}]