(* 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 = { IkBa\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t], IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[t], IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[t], JNK[t], JNK\[LetterSpace]P[t], LPS[t], MSK1[t], MSK1\[LetterSpace]P[t], TNFR1\[LetterSpace]EL[t], TNFR1\[LetterSpace]cytoplasm[t], TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], TNFa\[LetterSpace]TNFR1\[LetterSpace]cytoplasm[t], TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC[t], TNFa\[LetterSpace]mRNA\[LetterSpace]MC[t], TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]LSEC[t], TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]MC[t], p38[t], p38\[LetterSpace]P[t], p65\[LetterSpace]2P[t], p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t], p65\[LetterSpace]cytoplasm[t], p65\[LetterSpace]nucleus[t], s11[t], s14[t], s15[t], s16[t], s19[t], s3[t], s5[t], s8[t] }; initialValues = { IkBa\[LetterSpace]cytoplasm[0] == 0.424204224794365, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[0] == 1.62849759450558, IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[0] == 0.000211704687285678, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[0] == 28.4610527359882, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[0] == 28.4610527359882, JNK[0] == 0.983803981324523, JNK\[LetterSpace]P[0] == 0.0161960186754769, LPS[0] == 5.0, MSK1[0] == 14.9991737488977, MSK1\[LetterSpace]P[0] == 0.000826251102268073, TNFR1\[LetterSpace]EL[0] == 0.00122834645669292, TNFR1\[LetterSpace]cytoplasm[0] == 5.52755905511812*^-05, TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[0] == 0.0, TNFa\[LetterSpace]TNFR1\[LetterSpace]cytoplasm[0] == 0.0, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC[0] == 0.0, TNFa\[LetterSpace]mRNA\[LetterSpace]MC[0] == 0.0, TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]LSEC[0] == 0.0, TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]MC[0] == 0.0, p38[0] == 0.826703614971449, p38\[LetterSpace]P[0] == 0.173296385028551, p65\[LetterSpace]2P[0] == 0.153433031574853, p65\[LetterSpace]IkBa\[LetterSpace]nucleus[0] == 0.915071201909209, p65\[LetterSpace]cytoplasm[0] == 0.00026669174188777, p65\[LetterSpace]nucleus[0] == 0.0133264213400785, s11[0] == 0.848408327993857, s14[0] == 0.299987669185545, s15[0] == 1.23308144552886*^-05, s16[0] == 0.0, s19[0] == 0.045753415155894, s3[0] == 0.00365791495212158, s5[0] == 0.0636533313200277, s8[0] == 0.182873595048844 }; rates = { Complex\[LetterSpace]association\[LetterSpace]c, Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c, Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n, Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c, Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n, Complex\[LetterSpace]nuclear\[LetterSpace]shuttle, IKK\[LetterSpace]phosphorylation, IkBa\[LetterSpace]degradation\[LetterSpace]c, IkBa\[LetterSpace]degradation\[LetterSpace]n, IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]c, IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]n, IkBa\[LetterSpace]mRNA\[LetterSpace]transport, IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle, IkBa\[LetterSpace]phosphorylation, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]1, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]2, IkBa\[LetterSpace]translation, JNK\[LetterSpace]active\[LetterSpace]phosphorylation, JNK\[LetterSpace]basal\[LetterSpace]phosphorylation, JNK\[LetterSpace]dephosphorylation, LPS\[LetterSpace]LSEC\[LetterSpace]degradation, LPS\[LetterSpace]MC\[LetterSpace]degradation, MSK1\[LetterSpace]activation, MSK1\[LetterSpace]deactivation, TNFR1\[LetterSpace]shuttle, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation, TNFa\[LetterSpace]LSEC\[LetterSpace]translation, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation, TNFa\[LetterSpace]MC\[LetterSpace]translation, TNFa\[LetterSpace]TNFR1\[LetterSpace]association, TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation, TNFa\[LetterSpace]degradation, TNFa\[LetterSpace]internal\[LetterSpace]degradation, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation, TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation, basal\[LetterSpace]IKKb\[LetterSpace]phosphorylation, p38\[LetterSpace]active\[LetterSpace]phosphorylation, p38\[LetterSpace]basal\[LetterSpace]phosphorylation, p38\[LetterSpace]dephosphorylation, p65\[LetterSpace]2P\[LetterSpace]degradation, p65\[LetterSpace]2P\[LetterSpace]dephosphorylation\[LetterSpace]0, p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]c, p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]n, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n, p65\[LetterSpace]P\[LetterSpace]phosphorylation\[LetterSpace]MSK1, p65\[LetterSpace]nuclear\[LetterSpace]shuttle, p65\[LetterSpace]phosphorylation\[LetterSpace]MSK1, p65\[LetterSpace]translation, re10, re14, re15, re17, re19, re3 }; rateEquations = { Complex\[LetterSpace]association\[LetterSpace]c -> default*(IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k*IkBa\[LetterSpace]cytoplasm[t]*p65\[LetterSpace]cytoplasm[t] - IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k*s19[t]), Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c -> default*IkBa\[LetterSpace]deg\[LetterSpace]complex\[LetterSpace]cyt\[LetterSpace]k*s19[t], Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n -> c1*IkBa\[LetterSpace]deg\[LetterSpace]complex\[LetterSpace]nuc\[LetterSpace]k*p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t], Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c -> default*p65\[LetterSpace]degradation\[LetterSpace]k*s19[t], Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n -> c1*p65\[LetterSpace]degradation\[LetterSpace]k*p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t], Complex\[LetterSpace]nuclear\[LetterSpace]shuttle -> IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k*p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] - IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k*s19[t], IKK\[LetterSpace]phosphorylation -> default*Catalytic\[LetterSpace]activation\[LetterSpace]\[LetterSpace]irrev[s14[t], TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kM, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kA], IkBa\[LetterSpace]degradation\[LetterSpace]c -> default*IkBa\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k*IkBa\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]degradation\[LetterSpace]n -> c1*IkBa\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k*s11[t], IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]c -> default*IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k*IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]n -> c1*IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t], IkBa\[LetterSpace]mRNA\[LetterSpace]transport -> IkBa\[LetterSpace]mRNA\[LetterSpace]transport\[LetterSpace]k*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t], IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle -> IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k*IkBa\[LetterSpace]cytoplasm[t] - IkBa\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k*s11[t], IkBa\[LetterSpace]phosphorylation -> default*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[IkBa\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]phosphorylation\[LetterSpace]kM, IkBa\[LetterSpace]phosphorylation\[LetterSpace]vmax], IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]1 -> c1*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal*IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[t], IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]2 -> c1*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal*IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[t], IkBa\[LetterSpace]translation -> default*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]translation\[LetterSpace]Km, IkBa\[LetterSpace]translation\[LetterSpace]vmax], JNK\[LetterSpace]active\[LetterSpace]phosphorylation -> default*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[JNK[t], JNK\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax], JNK\[LetterSpace]basal\[LetterSpace]phosphorylation -> default*JNK\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k*JNK[t], JNK\[LetterSpace]dephosphorylation -> default*JNK\[LetterSpace]dephosphorylation\[LetterSpace]k*JNK\[LetterSpace]P[t], LPS\[LetterSpace]LSEC\[LetterSpace]degradation -> EL*LPS\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k*LPS[t], LPS\[LetterSpace]MC\[LetterSpace]degradation -> EL*LPS\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k*LPS[t], MSK1\[LetterSpace]activation -> c1*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[MSK1[t], MSK1\[LetterSpace]activation\[LetterSpace]kM, MSK1\[LetterSpace]activation\[LetterSpace]vmax], MSK1\[LetterSpace]deactivation -> c1*MSK1\[LetterSpace]dephosphorylation\[LetterSpace]k*MSK1\[LetterSpace]P[t], TNFR1\[LetterSpace]shuttle -> k\[LetterSpace]TNFR1\[LetterSpace]vessicle2outer\[LetterSpace]membrane\[LetterSpace]shuttle*TNFR1\[LetterSpace]cytoplasm[t] - k\[LetterSpace]TNFR1\[LetterSpace]outer\[LetterSpace]membrane2vessicle\[LetterSpace]shuttle*TNFR1\[LetterSpace]EL[t], TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation -> LSEC*TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k*TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]LSEC[t], TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation -> LSEC*TNFa\[LetterSpace]transcription\[LetterSpace]initiation[TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS, LPS[t], TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS], TNFa\[LetterSpace]LSEC\[LetterSpace]translation -> TNFa\[LetterSpace]LSEC\[LetterSpace]translation\[LetterSpace]k*TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC[t], TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation -> MC*TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k*TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]MC[t], TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation -> MC*TNFa\[LetterSpace]transcription\[LetterSpace]initiation[TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS, LPS[t], TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS], TNFa\[LetterSpace]MC\[LetterSpace]translation -> TNFa\[LetterSpace]MC\[LetterSpace]translation\[LetterSpace]k*TNFa\[LetterSpace]mRNA\[LetterSpace]MC[t], TNFa\[LetterSpace]TNFR1\[LetterSpace]association -> EL*(-(k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]dissociation*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t]) + k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association*s16[t]*TNFR1\[LetterSpace]EL[t]), TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation -> k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], TNFa\[LetterSpace]degradation -> EL*TNFa\[LetterSpace]degradation\[LetterSpace]k*s16[t], TNFa\[LetterSpace]internal\[LetterSpace]degradation -> default*k\[LetterSpace]TNFa\[LetterSpace]internal\[LetterSpace]degradation*TNFa\[LetterSpace]TNFR1\[LetterSpace]cytoplasm[t], TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation -> LSEC*TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k*TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC[t], TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation -> MC*TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k*TNFa\[LetterSpace]mRNA\[LetterSpace]MC[t], basal\[LetterSpace]IKKb\[LetterSpace]phosphorylation -> default*IKKb\[LetterSpace]basal\[LetterSpace]phosphoylation\[LetterSpace]k*s14[t], p38\[LetterSpace]active\[LetterSpace]phosphorylation -> default*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[p38[t], p38\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax], p38\[LetterSpace]basal\[LetterSpace]phosphorylation -> default*p38\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k*p38[t], p38\[LetterSpace]dephosphorylation -> default*p38\[LetterSpace]dephosphorylation\[LetterSpace]k*p38\[LetterSpace]P[t], p65\[LetterSpace]2P\[LetterSpace]degradation -> c1*p65\[LetterSpace]degradation\[LetterSpace]k*p65\[LetterSpace]2P[t], p65\[LetterSpace]2P\[LetterSpace]dephosphorylation\[LetterSpace]0 -> c1*p65\[LetterSpace]2P\[LetterSpace]dephosphorylation*p65\[LetterSpace]2P[t], p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]c -> default*p65\[LetterSpace]degradation\[LetterSpace]k*s3[t], p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]n -> c1*p65\[LetterSpace]degradation\[LetterSpace]k*s8[t], p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c -> c1*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k*s8[t], p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n -> default*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k*s3[t], p65\[LetterSpace]P\[LetterSpace]phosphorylation\[LetterSpace]MSK1 -> c1*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[s8[t], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kM, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax], p65\[LetterSpace]nuclear\[LetterSpace]shuttle -> p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k*p65\[LetterSpace]cytoplasm[t] - p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k*p65\[LetterSpace]nucleus[t], p65\[LetterSpace]phosphorylation\[LetterSpace]MSK1 -> c1*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[p65\[LetterSpace]nucleus[t], p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax], p65\[LetterSpace]translation -> default*p65\[LetterSpace]mRNA*p65\[LetterSpace]translation\[LetterSpace]k, re10 -> c1*activated\[LetterSpace]Transcription\[LetterSpace]\[LetterSpace]1[IkBa\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal, IkBa\[LetterSpace]transcription\[LetterSpace]kA2\[LetterSpace]p65\[LetterSpace]P, s8[t], IkBa\[LetterSpace]transcription\[LetterSpace]kA\[LetterSpace]p65, p65\[LetterSpace]nucleus[t], p65\[LetterSpace]2P[t], IkBa\[LetterSpace]transcription\[LetterSpace]kA3\[LetterSpace]p65\[LetterSpace]2P], re14 -> default*IkBa\[LetterSpace]p\[LetterSpace]active\[LetterSpace]degradation\[LetterSpace]k*s5[t], re15 -> default*Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[s19[t], IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kM, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]vmax], re17 -> default*IKKb\[LetterSpace]dephopshorylation\[LetterSpace]k*s15[t], re19 -> c1*(-(IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k*p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t]) + IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k*p65\[LetterSpace]nucleus[t]*s11[t]), re3 -> p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k*s3[t] - p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k*s8[t] }; parameters = { IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kA -> 0.477886214802612, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kM -> 5.0, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> 974.780691842994, IKKb\[LetterSpace]basal\[LetterSpace]phosphoylation\[LetterSpace]k -> 1.38817875620378*^-07, IKKb\[LetterSpace]dephopshorylation\[LetterSpace]k -> 0.00337720197636294, IkBa\[LetterSpace]deg\[LetterSpace]complex\[LetterSpace]cyt\[LetterSpace]k -> 1*^-05, IkBa\[LetterSpace]deg\[LetterSpace]complex\[LetterSpace]nuc\[LetterSpace]k -> 1.00000000002221*^-05, IkBa\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k -> 1.39705394692627*^-05, IkBa\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k -> 1*^-05, IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k -> 0.000999999999999778, IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k -> 5*^-05, IkBa\[LetterSpace]mRNA\[LetterSpace]transport\[LetterSpace]k -> 0.0001, IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k -> 0.005, IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k -> 0.0448147292565213, IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.000100000000000222, IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k -> 0.00199999999999978, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.5, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 19.4934850553468, IkBa\[LetterSpace]p\[LetterSpace]active\[LetterSpace]degradation\[LetterSpace]k -> 0.00337245067021992, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.100000000019363, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal -> 0.000929800339482254, IkBa\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 1*^-06, IkBa\[LetterSpace]transcription\[LetterSpace]kA2\[LetterSpace]p65\[LetterSpace]P -> 0.000100000000000222, IkBa\[LetterSpace]transcription\[LetterSpace]kA3\[LetterSpace]p65\[LetterSpace]2P -> 0.101069836420617, IkBa\[LetterSpace]transcription\[LetterSpace]kA\[LetterSpace]p65 -> 0.820652255737923, IkBa\[LetterSpace]translation\[LetterSpace]Km -> 4.74886881333832, IkBa\[LetterSpace]translation\[LetterSpace]vmax -> 0.000869911007708551, JNK\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 2.90744335946462, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 21203.3919242898, JNK\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 1.22484274931687*^-05, JNK\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000744013203133041, LPS\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.000197681956972319, LPS\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.000106908917668639, MSK1\[LetterSpace]activation\[LetterSpace]kM -> 0.332175, MSK1\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00999999999999978, MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 4.87343945878311*^-05, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k -> 1.00000000022204*^-06, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS -> 1.13302040932107, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 1.26292824593068*^-05, TNFa\[LetterSpace]LSEC\[LetterSpace]translation\[LetterSpace]k -> 2.24199564035721*^-05, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k -> 1.00000000022204*^-06, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS -> 0.855913113601739, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 6.27636799940596*^-06, TNFa\[LetterSpace]MC\[LetterSpace]translation\[LetterSpace]k -> 2.29979473005635*^-06, TNFa\[LetterSpace]degradation\[LetterSpace]k -> 0.000213041462325339, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.55139637922104, TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.0479810150640656, k\[LetterSpace]IkBa -> 154014.992996784, k\[LetterSpace]IkBa\[LetterSpace]P -> 593648.547667112, k\[LetterSpace]IkBa\[LetterSpace]mRNA -> 17213.2742721703, k\[LetterSpace]JNK\[LetterSpace]P -> 3238451.34489474, k\[LetterSpace]MSK\[LetterSpace]P -> 985394013.714541, k\[LetterSpace]TNFR1\[LetterSpace]outer\[LetterSpace]membrane2vessicle\[LetterSpace]shuttle -> 1.8*^-05, k\[LetterSpace]TNFR1\[LetterSpace]vessicle2outer\[LetterSpace]membrane\[LetterSpace]shuttle -> 0.0004, k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association -> 16.833, k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation -> 0.00077, k\[LetterSpace]TNFa\[LetterSpace]internal\[LetterSpace]degradation -> 8.80079235349194*^-05, k\[LetterSpace]p38\[LetterSpace]P -> 365327.672084138, k\[LetterSpace]p65 -> 0.102372088824604, k\[LetterSpace]p65\[LetterSpace]P -> 2888210.91907808, p38\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 2.90903981195396, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 231523.367249302, p38\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 0.00011778990894447, p38\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000561912145573657, p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM -> 5.29277823905396, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> 4.39829577360994, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000837226665087327, p65\[LetterSpace]degradation\[LetterSpace]k -> 1.328912101*^-05, p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k -> 0.0005, p65\[LetterSpace]translation\[LetterSpace]k -> 1.68*^-06, p65\[LetterSpace]mRNA -> 1.0, EL -> 2.59*^-06, LSEC -> 4.7*^-07, MC -> 3.5*^-07, c1 -> 8*^-07, default -> 1.3*^-05 }; assignments = { TNFa\[LetterSpace]transcription\[LetterSpace]initiation[kbasal_,vA_,A_,kA_] -> kbasal + (A*vA)/(A + kA), Henri\[LetterSpace]Michaelis\[LetterSpace]Menten\[LetterSpace]\[LetterSpace]irreversible[substrate_,Km_,V_] -> (substrate*V)/(Km + substrate), activated\[LetterSpace]Transcription\[LetterSpace]\[LetterSpace]1[k1_,k2_,Modifier_,k3_,Modifier2_,Modifier3_,k4_] -> k1 + k2*Modifier + k3*Modifier2 + k4*Modifier3, Catalytic\[LetterSpace]activation\[LetterSpace]\[LetterSpace]irrev[substrate_,Activator_,Kms_,V_,Ka_] -> (Activator*substrate*V)/((Activator + Ka)*(Kms + substrate)), p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]p65\[LetterSpace]P -> default*s3[t] + c1*s8[t], p65\[LetterSpace]2P\[LetterSpace]dephosphorylation -> p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k, unscaled\[LetterSpace]JNK\[LetterSpace]P -> default*JNK\[LetterSpace]P[t], MSK1\[LetterSpace]activation\[LetterSpace]vmax -> MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat*p38\[LetterSpace]P[t], unscaled\[LetterSpace]MSK\[LetterSpace]P -> c1*MSK1\[LetterSpace]P[t], scaled\[LetterSpace]MSK\[LetterSpace]P -> k\[LetterSpace]MSK\[LetterSpace]P*unscaled\[LetterSpace]MSK\[LetterSpace]P, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]p38\[LetterSpace]P -> default*p38\[LetterSpace]P[t], scaled\[LetterSpace]p65 -> k\[LetterSpace]p65*unscaled\[LetterSpace]p65, IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k -> 0.003*IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k, scaled\[LetterSpace]p65\[LetterSpace]P -> k\[LetterSpace]p65\[LetterSpace]P*unscaled\[LetterSpace]p65\[LetterSpace]P, IkBa\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.5*IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat, IkBa\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], unscaled\[LetterSpace]IkBa\[LetterSpace]P -> default*s5[t], scaled\[LetterSpace]IkBa -> k\[LetterSpace]IkBa*unscaled\[LetterSpace]IkBa, scaled\[LetterSpace]IkBa\[LetterSpace]P -> k\[LetterSpace]IkBa\[LetterSpace]P*unscaled\[LetterSpace]IkBa\[LetterSpace]P, scaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> k\[LetterSpace]IkBa\[LetterSpace]mRNA*unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA, unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> default*IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t] + c1*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t], unscaled\[LetterSpace]IkBa -> c1*(p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] + s11[t]) + default*(IkBa\[LetterSpace]cytoplasm[t] + s19[t] + s5[t]), scaled\[LetterSpace]p38\[LetterSpace]P -> k\[LetterSpace]p38\[LetterSpace]P*unscaled\[LetterSpace]p38\[LetterSpace]P, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kM -> p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat, k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]dissociation -> 0.000019*k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association, scaled\[LetterSpace]JNK\[LetterSpace]P -> k\[LetterSpace]JNK\[LetterSpace]P*unscaled\[LetterSpace]JNK\[LetterSpace]P, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.02*p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], unscaled\[LetterSpace]p65 -> default*(p65\[LetterSpace]cytoplasm[t] + s19[t] + s3[t]) + c1*(p65\[LetterSpace]2P[t] + p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] + p65\[LetterSpace]nucleus[t] + s8[t]) }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { IkBa\[LetterSpace]cytoplasm'[t] == 1.0*IkBa\[LetterSpace]translation +1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c -1.0*IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*IkBa\[LetterSpace]degradation\[LetterSpace]c -1.0*IkBa\[LetterSpace]phosphorylation -1.0*Complex\[LetterSpace]association\[LetterSpace]c, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm'[t] == 1.0*IkBa\[LetterSpace]translation +1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]transport -1.0*IkBa\[LetterSpace]translation -1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]c, IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus'[t] == 1.0*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]2 -1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]n -1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]transport, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1'[t] == 1.0*re10 -1.0*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]1, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2'[t] == 1.0*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]1 -1.0*IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]2, JNK'[t] == 1.0*JNK\[LetterSpace]dephosphorylation -1.0*JNK\[LetterSpace]active\[LetterSpace]phosphorylation -1.0*JNK\[LetterSpace]basal\[LetterSpace]phosphorylation, JNK\[LetterSpace]P'[t] == 1.0*JNK\[LetterSpace]basal\[LetterSpace]phosphorylation +1.0*JNK\[LetterSpace]active\[LetterSpace]phosphorylation -1.0*JNK\[LetterSpace]dephosphorylation, LPS'[t] == -1.0*LPS\[LetterSpace]LSEC\[LetterSpace]degradation -1.0*LPS\[LetterSpace]MC\[LetterSpace]degradation, MSK1'[t] == 1.0*MSK1\[LetterSpace]deactivation -1.0*MSK1\[LetterSpace]activation, MSK1\[LetterSpace]P'[t] == 1.0*MSK1\[LetterSpace]activation -1.0*MSK1\[LetterSpace]deactivation, TNFR1\[LetterSpace]EL'[t] == 1.0*TNFR1\[LetterSpace]shuttle -1.0*TNFa\[LetterSpace]TNFR1\[LetterSpace]association, TNFR1\[LetterSpace]cytoplasm'[t] == 1.0*TNFa\[LetterSpace]internal\[LetterSpace]degradation -1.0*TNFR1\[LetterSpace]shuttle, TNFa\[LetterSpace]TNFR1\[LetterSpace]EL'[t] == 1.0*TNFa\[LetterSpace]TNFR1\[LetterSpace]association -1.0*TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation, TNFa\[LetterSpace]TNFR1\[LetterSpace]cytoplasm'[t] == 1.0*TNFa\[LetterSpace]TNFR1\[LetterSpace]internalisation -1.0*TNFa\[LetterSpace]internal\[LetterSpace]degradation, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC'[t] == 1.0*TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation +1.0*TNFa\[LetterSpace]LSEC\[LetterSpace]translation -1.0*TNFa\[LetterSpace]LSEC\[LetterSpace]translation -1.0*TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation, TNFa\[LetterSpace]mRNA\[LetterSpace]MC'[t] == 1.0*TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation +1.0*TNFa\[LetterSpace]MC\[LetterSpace]translation -1.0*TNFa\[LetterSpace]MC\[LetterSpace]translation -1.0*TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation, TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]LSEC'[t] == 2.0*TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation -1.0*TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation, TNFa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]MC'[t] == 2.0*TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation -1.0*TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]elongation, p38'[t] == 1.0*p38\[LetterSpace]dephosphorylation -1.0*p38\[LetterSpace]basal\[LetterSpace]phosphorylation -1.0*p38\[LetterSpace]active\[LetterSpace]phosphorylation, p38\[LetterSpace]P'[t] == 1.0*p38\[LetterSpace]active\[LetterSpace]phosphorylation +1.0*p38\[LetterSpace]basal\[LetterSpace]phosphorylation -1.0*p38\[LetterSpace]dephosphorylation, p65\[LetterSpace]2P'[t] == 1.0*p65\[LetterSpace]P\[LetterSpace]phosphorylation\[LetterSpace]MSK1 +1.0*p65\[LetterSpace]phosphorylation\[LetterSpace]MSK1 -1.0*p65\[LetterSpace]2P\[LetterSpace]degradation -1.0*p65\[LetterSpace]2P\[LetterSpace]dephosphorylation\[LetterSpace]0, p65\[LetterSpace]IkBa\[LetterSpace]nucleus'[t] == 1.0*re19 -1.0*Complex\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n -1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n, p65\[LetterSpace]cytoplasm'[t] == 1.0*p65\[LetterSpace]translation +1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c +1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n -1.0*Complex\[LetterSpace]association\[LetterSpace]c -1.0*p65\[LetterSpace]nuclear\[LetterSpace]shuttle, p65\[LetterSpace]nucleus'[t] == 1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n +1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c +1.0*p65\[LetterSpace]2P\[LetterSpace]dephosphorylation\[LetterSpace]0 +1.0*p65\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*re19 -1.0*p65\[LetterSpace]phosphorylation\[LetterSpace]MSK1, s11'[t] == 1.0*IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle +1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n -1.0*re19 -1.0*IkBa\[LetterSpace]degradation\[LetterSpace]n, s14'[t] == 1.0*re17 -1.0*IKK\[LetterSpace]phosphorylation -1.0*basal\[LetterSpace]IKKb\[LetterSpace]phosphorylation, s15'[t] == 1.0*IKK\[LetterSpace]phosphorylation +1.0*basal\[LetterSpace]IKKb\[LetterSpace]phosphorylation -1.0*re17, s16'[t] == 1.0*TNFa\[LetterSpace]LSEC\[LetterSpace]translation +1.0*TNFa\[LetterSpace]MC\[LetterSpace]translation -1.0*TNFa\[LetterSpace]degradation -1.0*TNFa\[LetterSpace]TNFR1\[LetterSpace]association, s19'[t] == 1.0*Complex\[LetterSpace]nuclear\[LetterSpace]shuttle +1.0*Complex\[LetterSpace]association\[LetterSpace]c -1.0*re15 -1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c -1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c, s3'[t] == 1.0*re15 -1.0*re3 -1.0*p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]c -1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n, s5'[t] == 1.0*re15 +1.0*IkBa\[LetterSpace]phosphorylation -1.0*re14, s8'[t] == 1.0*re3 -1.0*p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]n -1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c -1.0*p65\[LetterSpace]P\[LetterSpace]phosphorylation\[LetterSpace]MSK1 }; 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]}]