(* 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.564988584462963, IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[0] == 2.79615255717589, IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[0] == 8.66748867250049*^-05, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]1[0] == 5.33067431606035, IkBa\[LetterSpace]pre\[LetterSpace]mRNA\[LetterSpace]2[0] == 5.33067431606035, JNK[0] == 0.610494254930111, JNK\[LetterSpace]P[0] == 0.38950574506989, LPS[0] == 5.0, MSK1[0] == 14.9534172138408, MSK1\[LetterSpace]P[0] == 0.0465827861591864, TNFR1\[LetterSpace]EL[0] == 0.00122834645669316, TNFR1\[LetterSpace]cytoplasm[0] == 5.52755905511923*^-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.573949218303294, p38\[LetterSpace]P[0] == 0.426050781696706, p65\[LetterSpace]2P[0] == 0.228927278358454, p65\[LetterSpace]IkBa\[LetterSpace]nucleus[0] == 0.84799128781868, p65\[LetterSpace]cytoplasm[0] == 0.000334685105954136, p65\[LetterSpace]nucleus[0] == 0.0167298646806817, s11[0] == 1.12997707322809, s14[0] == 0.299976297560879, s15[0] == 2.37024391214021*^-05, s16[0] == 0.0, s19[0] == 0.042399452680612, s3[0] == 0.00435353355955181, s5[0] == 0.073948097360867, s8[0] == 0.217657350208359 }; 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.00100000751793317, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kM -> 5.0, IKKb\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> 999.994251892825, IKKb\[LetterSpace]basal\[LetterSpace]phosphoylation\[LetterSpace]k -> 7.50988514561613*^-05, IKKb\[LetterSpace]dephopshorylation\[LetterSpace]k -> 0.950445449749497, 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.00023844546542126, 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.0182773778651534, 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 -> 9.35331903485181, IkBa\[LetterSpace]p\[LetterSpace]active\[LetterSpace]degradation\[LetterSpace]k -> 0.0027815203293068, IkBa\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.1, IkBa\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]kbasal -> 0.00203245678351198, 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.644877105549519, IkBa\[LetterSpace]translation\[LetterSpace]Km -> 8.62022804518077, IkBa\[LetterSpace]translation\[LetterSpace]vmax -> 0.000869570552810696, JNK\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 0.00100000000000022, JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 128223.281718797, JNK\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 0.000297281217739997, JNK\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000465945567751045, LPS\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.000237613611887944, LPS\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.000119938547034168, MSK1\[LetterSpace]activation\[LetterSpace]kM -> 0.332175, MSK1\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.00894728276588594, MSK1\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 0.000999993266243514, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]elongation\[LetterSpace]k -> 1.00000000022204*^-06, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kA\[LetterSpace]LPS -> 1.06727068122671, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]LSEC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 1.21773968583006*^-05, TNFa\[LetterSpace]LSEC\[LetterSpace]translation\[LetterSpace]k -> 2.16264815881264*^-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.820922917690253, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]kbasal -> 0.0, TNFa\[LetterSpace]MC\[LetterSpace]transcription\[LetterSpace]initiation\[LetterSpace]vA\[LetterSpace]LPS -> 5.85088364603215*^-06, TNFa\[LetterSpace]MC\[LetterSpace]translation\[LetterSpace]k -> 2.14389156706575*^-06, TNFa\[LetterSpace]degradation\[LetterSpace]k -> 0.000155294844602142, TNFa\[LetterSpace]mRNA\[LetterSpace]LSEC\[LetterSpace]degradation\[LetterSpace]k -> 0.595159854597121, TNFa\[LetterSpace]mRNA\[LetterSpace]MC\[LetterSpace]degradation\[LetterSpace]k -> 0.0518402782192837, k\[LetterSpace]IkBa -> 123151.961018769, k\[LetterSpace]IkBa\[LetterSpace]P -> 511003.109490479, k\[LetterSpace]IkBa\[LetterSpace]mRNA -> 10025.184947558, k\[LetterSpace]JNK\[LetterSpace]P -> 134657.881495759, k\[LetterSpace]MSK\[LetterSpace]P -> 17478192.2064024, 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.99671747509247*^-05, k\[LetterSpace]p38\[LetterSpace]P -> 148597.227473568, k\[LetterSpace]p65 -> 0.102372088824604, k\[LetterSpace]p65\[LetterSpace]P -> 2426666.054545, p38\[LetterSpace]active\[LetterSpace]phoshorylation\[LetterSpace]Km -> 0.10616767037577, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat -> 166798.143477709, p38\[LetterSpace]basal\[LetterSpace]phosphorylation\[LetterSpace]k -> 0.000370463577422213, p38\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000499065580456255, p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM -> 0.945158930722195, p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> 0.0135882089748922, p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k -> 0.000552351266083683, 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)), scaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> k\[LetterSpace]IkBa\[LetterSpace]mRNA*unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA, 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]), k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]dissociation -> 0.000019*k\[LetterSpace]TNFa\[LetterSpace]TNFR1\[LetterSpace]association, scaled\[LetterSpace]IkBa -> k\[LetterSpace]IkBa*unscaled\[LetterSpace]IkBa, unscaled\[LetterSpace]IkBa -> c1*(p65\[LetterSpace]IkBa\[LetterSpace]nucleus[t] + s11[t]) + default*(IkBa\[LetterSpace]cytoplasm[t] + s19[t] + s5[t]), p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kM -> p65\[LetterSpace]MSK1\[LetterSpace]phosphorylation\[LetterSpace]kM, scaled\[LetterSpace]p38\[LetterSpace]P -> k\[LetterSpace]p38\[LetterSpace]P*unscaled\[LetterSpace]p38\[LetterSpace]P, p38\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> p38\[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], p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]JNK\[LetterSpace]P -> default*JNK\[LetterSpace]P[t], IkBa\[LetterSpace]phosphorylation\[LetterSpace]kcat -> IkBa\[LetterSpace]p65\[LetterSpace]phosphorylation\[LetterSpace]kcat, scaled\[LetterSpace]JNK\[LetterSpace]P -> k\[LetterSpace]JNK\[LetterSpace]P*unscaled\[LetterSpace]JNK\[LetterSpace]P, p65\[LetterSpace]P\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat, p65\[LetterSpace]2P\[LetterSpace]dephosphorylation -> p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]k, 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]p38\[LetterSpace]P -> default*p38\[LetterSpace]P[t], JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]vmax -> JNK\[LetterSpace]active\[LetterSpace]phosphorylation\[LetterSpace]kcat*TNFa\[LetterSpace]TNFR1\[LetterSpace]EL[t], p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]vmax -> p65\[LetterSpace]MSK1\[LetterSpace]phosphrylation\[LetterSpace]kcat*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]MSK\[LetterSpace]P -> c1*MSK1\[LetterSpace]P[t], unscaled\[LetterSpace]IkBa\[LetterSpace]P -> default*s5[t], scaled\[LetterSpace]IkBa\[LetterSpace]P -> k\[LetterSpace]IkBa\[LetterSpace]P*unscaled\[LetterSpace]IkBa\[LetterSpace]P, IkBa\[LetterSpace]nuclear\[LetterSpace]export\[LetterSpace]k -> 0.5*IkBa\[LetterSpace]nuclear\[LetterSpace]import\[LetterSpace]k, IkBa\[LetterSpace]p65\[LetterSpace]dissociation\[LetterSpace]k -> 0.003*IkBa\[LetterSpace]p65\[LetterSpace]association\[LetterSpace]k, unscaled\[LetterSpace]IkBa\[LetterSpace]mRNA -> default*IkBa\[LetterSpace]mRNA\[LetterSpace]cytoplasm[t] + c1*IkBa\[LetterSpace]mRNA\[LetterSpace]nucleus[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], unscaled\[LetterSpace]p65\[LetterSpace]P -> default*s3[t] + c1*s8[t], scaled\[LetterSpace]p65\[LetterSpace]P -> k\[LetterSpace]p65\[LetterSpace]P*unscaled\[LetterSpace]p65\[LetterSpace]P }; 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]2P\[LetterSpace]dephosphorylation\[LetterSpace]0 +1.0*p65\[LetterSpace]nuclear\[LetterSpace]shuttle +1.0*Complex\[LetterSpace]degradation\[LetterSpace]IkBa\[LetterSpace]n +1.0*p65\[LetterSpace]P\[LetterSpace]dephosphorylation\[LetterSpace]c -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]}]