(* 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.554236210497145, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[0] == 1.20735759082122, IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[0] == 3.7519558085674*^-05, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[0] == 2.28074483990426, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[0] == 2.28074483990426, JNK[0] == 0.60623758410987, JNK\[LetterSpace]P[0] == 0.39376241589013, LPS[0] == 5.0, MSK1[0] == 14.99559259493, MSK1\[LetterSpace]P[0] == 0.0044074050699931, TNFR1\[LetterSpace]EL[0] == 0.00122834645669292, TNFR1\[LetterSpace]cytoplasm[0] == 5.52755905511814*^-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.647950363084063, p38\[LetterSpace]P[0] == 0.352049636915937, p65\[LetterSpace]2P[0] == 0.211938917601913, p65\[LetterSpace]IkBa\[LetterSpace]nucleus[0] == 0.857207842972257, p65\[LetterSpace]cytoplasm[0] == 0.00034883688586427, p65\[LetterSpace]nucleus[0] == 0.0174373858516052, s11[0] == 1.1084723259792, s14[0] == 0.299994079312355, s15[0] == 5.92068764529143*^-06, s16[0] == 0.0, s19[0] == 0.0428602812425719, s3[0] == 0.00435780683888014, s5[0] == 0.0725590666809302, s8[0] == 0.217871250890394 }; 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.00100000000000022, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kM -> 5.0, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> 1000.0, IKKb\[LetterSpace]basal\[LetterSpace]phosphoylation\[LetterSpace]k -> 7.52966991388947*^-05, IKKb\[LetterSpace]dephopshorylation\[LetterSpace]k -> 3.81519264091726, 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.00000000002221*^-05, IkBa\[LetterSpace]deg\[LetterSpace]nuc\[LetterSpace]k -> 1*^-05, IkBa\[LetterSpace]mRNA\[LetterSpace]deg\[LetterSpace]cyt\[LetterSpace]k -> 0.000239044329094007, 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.0177373504282182, 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 -> 36.8226975109864, IkBa\[LetterSpace]p\[LetterSpace]active\[LetterSpace]degradation\[LetterSpace]k -> 0.00278262839781987, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.1, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal -> 0.00205632237093388, 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.000100000000000222, IkBa\[LetterSpace]transcription\[LetterSpace]kA\[LetterSpace]p65 -> 0.26643704849888, IkBa\[LetterSpace]translation\[LetterSpace]Km -> 3.67448875058377, IkBa\[LetterSpace]translation\[LetterSpace]vmax -> 0.000845419331920816, JNK\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 0.00100000000000022, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 129582.438480298, JNK\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 0.000299993806270227, JNK\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000461871202080221, LPS\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.000238245550264726, LPS\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.000120112331975526, MSK1\[LetterSpace]activation\[LetterSpace]kM -> 0.332175, MSK1\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00998972839132086, MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 0.000127834483364289, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k -> 1.00000000022204*^-06, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS -> 1.06602156258362, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 1.21692029414147*^-05, TNFa\[LetterSpace]LSEC\[LetterSpace]translation\[LetterSpace]k -> 2.16196526283029*^-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.82029823696587, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 5.84583032546336*^-06, TNFa\[LetterSpace]MC\[LetterSpace]translation\[LetterSpace]k -> 2.14203984011005*^-06, TNFa\[LetterSpace]degradation\[LetterSpace]k -> 0.00015479362017954, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.595075132442134, TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.0518587991406785, k\[LetterSpace]IkBa -> 125089.189993148, k\[LetterSpace]IkBa\[LetterSpace]P -> 520785.470663157, k\[LetterSpace]IkBa\[LetterSpace]mRNA -> 23217.6006562486, k\[LetterSpace]JNK\[LetterSpace]P -> 133202.195905293, k\[LetterSpace]MSK\[LetterSpace]P -> 184730669.650311, 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.90167499840012*^-05, k\[LetterSpace]p38\[LetterSpace]P -> 179832.496001677, k\[LetterSpace]p65 -> 0.102372088824604, k\[LetterSpace]p65\[LetterSpace]P -> 2424284.31087238, p38\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 0.159221354746047, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 138790.20176403, p38\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 0.000304993336054619, p38\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000561342839509866, p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM -> 1.06554415971021, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> 0.15005984082624, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000566704472713779, 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 = { 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, 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)), unscaled\[LetterSpace]p65\[LetterSpace]P -> default*s3[t] + c1*s8[t], unscaled\[LetterSpace]IkBa\[LetterSpace]P -> default*s5[t], scaled\[LetterSpace]IkBa\[LetterSpace]P -> k\[LetterSpace]IkBa\[LetterSpace]P*unscaled\[LetterSpace]IkBa\[LetterSpace]P, MSK1\[LetterSpace]activation\[LetterSpace]vmax -> MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat*p38\[LetterSpace]P[t], unscaled\[LetterSpace]p38\[LetterSpace]P -> default*p38\[LetterSpace]P[t], scaled\[LetterSpace]p65\[LetterSpace]P -> k\[LetterSpace]p65\[LetterSpace]P*unscaled\[LetterSpace]p65\[LetterSpace]P, scaled\[LetterSpace]IkBa -> k\[LetterSpace]IkBa*unscaled\[LetterSpace]IkBa, p65\[LetterSpace]2P\[LetterSpace]dephosphorylation -> p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k, k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]dissociation -> 0.000019*k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association, IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], IkBa\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.5*IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, p65\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.02*p65\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, scaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> k\[LetterSpace]IkBa\[LetterSpace]mRNA*unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA, scaled\[LetterSpace]p38\[LetterSpace]P -> k\[LetterSpace]p38\[LetterSpace]P*unscaled\[LetterSpace]p38\[LetterSpace]P, 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, IkBa\[LetterSpace]phosphorylation\[LetterSpace]vmax -> IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat*s15[t], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], scaled\[LetterSpace]JNK\[LetterSpace]P -> k\[LetterSpace]JNK\[LetterSpace]P*unscaled\[LetterSpace]JNK\[LetterSpace]P, unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> default*IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t] + c1*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[t], scaled\[LetterSpace]p65 -> k\[LetterSpace]p65*unscaled\[LetterSpace]p65, 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]), JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k -> 0.003*IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k, unscaled\[LetterSpace]IkBa -> c1*(p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] + s11[t]) + default*(IkBa\[LetterSpace]cytoplasm[t] + s19[t] + s5[t]), p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]JNK\[LetterSpace]P -> default*JNK\[LetterSpace]P[t], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kM -> p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM }; 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]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]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*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c +1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n +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*basal\[LetterSpace]IKKb\[LetterSpace]phosphorylation +1.0*IKK\[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]}]