(* 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.402845161595828, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[0] == 6.28632882916418, IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[0] == 0.000157166190926943, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[0] == 7.18140328980903, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[0] == 7.18140328980903, JNK[0] == 0.980466198972063, JNK\[LetterSpace]P[0] == 0.019533801027937, LPS[0] == 5.0, MSK1[0] == 14.9999306731002, MSK1\[LetterSpace]P[0] == 6.93268997811892*^-05, TNFR1\[LetterSpace]EL[0] == 0.00122834645669288, TNFR1\[LetterSpace]cytoplasm[0] == 5.52755905511796*^-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.929650735059334, p38\[LetterSpace]P[0] == 0.0703492649406657, p65\[LetterSpace]2P[0] == 2.189217886387*^-07, p65\[LetterSpace]IkBa\[LetterSpace]nucleus[0] == 0.995652059221173, p65\[LetterSpace]cytoplasm[0] == 0.00082222617042193, p65\[LetterSpace]nucleus[0] == 0.0411030139624838, s11[0] == 0.805690176476432, s14[0] == 0.299989869269092, s15[0] == 1.01307309084452*^-05, s16[0] == 0.0, s19[0] == 0.0497846885277208, s3[0] == 0.00376885551788092, s5[0] == 0.0566777858683295, s8[0] == 0.188414573904168 }; 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.243146348449513, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kM -> 5.0, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> 277.103899203292, IKKb\[LetterSpace]basal\[LetterSpace]phosphoylation\[LetterSpace]k -> 1.02344651825526*^-07, IKKb\[LetterSpace]dephopshorylation\[LetterSpace]k -> 0.00303061634929536, 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 -> 4.9999999999778*^-05, IkBa\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k -> 1*^-05, IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k -> 0.000192317445092793, 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.0153889708315641, IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.000100004543264477, IkBa\[LetterSpace]p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k -> 0.00199999999999978, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.372292128739572, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 24.0559815899914, IkBa\[LetterSpace]p\[LetterSpace]active\[LetterSpace]degradation\[LetterSpace]k -> 0.00395159188864045, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.10004489890377, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal -> 0.00273564663219183, 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.000100392930213726, IkBa\[LetterSpace]transcription\[LetterSpace]kA\[LetterSpace]p65 -> 0.477481779382945, IkBa\[LetterSpace]translation\[LetterSpace]Km -> 10.0, IkBa\[LetterSpace]translation\[LetterSpace]vmax -> 0.00063658986783369, JNK\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 5.53631315145271, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 42478.015487462, JNK\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 1.40201022425048*^-05, JNK\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00070371538725355, LPS\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.000281375430692567, LPS\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.000133176015283507, MSK1\[LetterSpace]activation\[LetterSpace]kM -> 0.332175, MSK1\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00999999999999978, MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 1.00729056845474*^-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.02438212695228, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 1.19438296662267*^-05, TNFa\[LetterSpace]LSEC\[LetterSpace]translation\[LetterSpace]k -> 2.12090661254614*^-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.777847391047605, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 5.6070630385483*^-06, TNFa\[LetterSpace]MC\[LetterSpace]translation\[LetterSpace]k -> 2.05455299505843*^-06, TNFa\[LetterSpace]degradation\[LetterSpace]k -> 0.000125263672194682, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.605694690908452, TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.0541021024827527, k\[LetterSpace]IkBa -> 159472.010610349, k\[LetterSpace]IkBa\[LetterSpace]P -> 666711.077600919, k\[LetterSpace]IkBa\[LetterSpace]mRNA -> 4459.19353465518, k\[LetterSpace]JNK\[LetterSpace]P -> 2685090.23853192, k\[LetterSpace]MSK\[LetterSpace]P -> 11744112207.0905, 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 -> 6.9007617563905*^-05, k\[LetterSpace]p38\[LetterSpace]P -> 899937.831283299, k\[LetterSpace]p65 -> 0.102372088824604, k\[LetterSpace]p65\[LetterSpace]P -> 2803253.51262869, p38\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 9.99601266623388, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 284561.328834224, p38\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 8.04034531971613*^-05, p38\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00106251471751828, p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.00550375653893175, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> 3.18759281450581*^-06, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00185771613858199, 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), Catalytic\[LetterSpace]activation\[LetterSpace]\[LetterSpace]irrev[substrate_,Activator_,Kms_,V_,Ka_] -> (Activator*substrate*V)/((Activator + Ka)*(Kms + substrate)), 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, p65\[LetterSpace]2P\[LetterSpace]dephosphorylation -> p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k, scaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> k\[LetterSpace]IkBa\[LetterSpace]mRNA*unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA, unscaled\[LetterSpace]p38\[LetterSpace]P -> default*p38\[LetterSpace]P[t], k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]dissociation -> 0.000019*k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat, IkBa\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], unscaled\[LetterSpace]IkBa\[LetterSpace]P -> default*s5[t], scaled\[LetterSpace]p65 -> k\[LetterSpace]p65*unscaled\[LetterSpace]p65, unscaled\[LetterSpace]p65\[LetterSpace]P -> default*s3[t] + c1*s8[t], unscaled\[LetterSpace]MSK\[LetterSpace]P -> c1*MSK1\[LetterSpace]P[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, scaled\[LetterSpace]p65\[LetterSpace]P -> k\[LetterSpace]p65\[LetterSpace]P*unscaled\[LetterSpace]p65\[LetterSpace]P, p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.02*p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, scaled\[LetterSpace]JNK\[LetterSpace]P -> k\[LetterSpace]JNK\[LetterSpace]P*unscaled\[LetterSpace]JNK\[LetterSpace]P, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], MSK1\[LetterSpace]activation\[LetterSpace]vmax -> MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat*p38\[LetterSpace]P[t], IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k -> 0.003*IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k, IkBa\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.5*IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, unscaled\[LetterSpace]JNK\[LetterSpace]P -> default*JNK\[LetterSpace]P[t], IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat, scaled\[LetterSpace]IkBa -> k\[LetterSpace]IkBa*unscaled\[LetterSpace]IkBa, scaled\[LetterSpace]IkBa\[LetterSpace]P -> k\[LetterSpace]IkBa\[LetterSpace]P*unscaled\[LetterSpace]IkBa\[LetterSpace]P, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], scaled\[LetterSpace]MSK\[LetterSpace]P -> k\[LetterSpace]MSK\[LetterSpace]P*unscaled\[LetterSpace]MSK\[LetterSpace]P, 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]), unscaled\[LetterSpace]IkBa -> c1*(p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] + s11[t]) + default*(IkBa\[LetterSpace]cytoplasm[t] + s19[t] + s5[t]), p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> default*IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t] + c1*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t] }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { IkBa\[LetterSpace]cytoplasm'[t] == 1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c +1.0*IkBa\[LetterSpace]translation -1.0*Complex\[LetterSpace]association\[LetterSpace]c -1.0*IkBa\[LetterSpace]degradation\[LetterSpace]c -1.0*IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*IkBa\[LetterSpace]phosphorylation, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm'[t] == 1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]transport +1.0*IkBa\[LetterSpace]translation -1.0*IkBa\[LetterSpace]mRNA\[LetterSpace]degradation\[LetterSpace]c -1.0*IkBa\[LetterSpace]translation, 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]active\[LetterSpace]phosphorylation +1.0*JNK\[LetterSpace]basal\[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]active\[LetterSpace]phosphorylation -1.0*p38\[LetterSpace]basal\[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]degradation\[LetterSpace]IkBa\[LetterSpace]n -1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n -1.0*Complex\[LetterSpace]nuclear\[LetterSpace]shuttle, p65\[LetterSpace]cytoplasm'[t] == 1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c +1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n +1.0*p65\[LetterSpace]translation -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]2P\[LetterSpace]dephosphorylation\[LetterSpace]0 +1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c +1.0*p65\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*p65\[LetterSpace]phosphorylation\[LetterSpace]MSK1 -1.0*re19, s11'[t] == 1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]n +1.0*IkBa\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*IkBa\[LetterSpace]degradation\[LetterSpace]n -1.0*re19, 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]TNFR1\[LetterSpace]association -1.0*TNFa\[LetterSpace]degradation, s19'[t] == 1.0*Complex\[LetterSpace]association\[LetterSpace]c +1.0*Complex\[LetterSpace]nuclear\[LetterSpace]shuttle -1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]c -1.0*Complex\[LetterSpace]degradation\[LetterSpace]p65\[LetterSpace]c -1.0*re15, s3'[t] == 1.0*re15 -1.0*p65\[LetterSpace]P\[LetterSpace]degradation\[LetterSpace]c -1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]n -1.0*re3, s5'[t] == 1.0*IkBa\[LetterSpace]phosphorylation +1.0*re15 -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]}]