(* 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 = {
IKK[t], IkBaIKK[t], IkBaIKKNFkB[t], IkBaNFkB\[LetterSpace]cytoplasm[t], IkBaNFkB\[LetterSpace]nucleus[t], IkBa\[LetterSpace]cytoplasm[t], IkBa\[LetterSpace]mRNA[t], IkBa\[LetterSpace]nucleus[t], IkBbIKK[t], IkBbIKKNFkB[t], IkBbNFkB\[LetterSpace]cytoplasm[t], IkBbNFkB\[LetterSpace]nucleus[t], IkBb\[LetterSpace]cytoplasm[t], IkBb\[LetterSpace]mRNA[t], IkBb\[LetterSpace]nucleus[t], IkBeIKK[t], IkBeIKKNFkB[t], IkBeNFkB\[LetterSpace]cytoplasm[t], IkBeNFkB\[LetterSpace]nucleus[t], IkBe\[LetterSpace]cytoplasm[t], IkBe\[LetterSpace]mRNA[t], IkBe\[LetterSpace]nucleus[t], NFkB\[LetterSpace]cytoplasm[t], NFkB\[LetterSpace]nucleus[t]
};

initialValues = {
IKK[0] == 0.001, IkBaIKK[0] == 0.0, IkBaIKKNFkB[0] == 0.0, IkBaNFkB\[LetterSpace]cytoplasm[0] == 0.06, IkBaNFkB\[LetterSpace]nucleus[0] == 0.0, IkBa\[LetterSpace]cytoplasm[0] == 0.0, IkBa\[LetterSpace]mRNA[0] == 0.0, IkBa\[LetterSpace]nucleus[0] == 0.0, IkBbIKK[0] == 0.0, IkBbIKKNFkB[0] == 0.0, IkBbNFkB\[LetterSpace]cytoplasm[0] == 0.03, IkBbNFkB\[LetterSpace]nucleus[0] == 0.0, IkBb\[LetterSpace]cytoplasm[0] == 0.0, IkBb\[LetterSpace]mRNA[0] == 0.0, IkBb\[LetterSpace]nucleus[0] == 0.0, IkBeIKK[0] == 0.0, IkBeIKKNFkB[0] == 0.0, IkBeNFkB\[LetterSpace]cytoplasm[0] == 0.01, IkBeNFkB\[LetterSpace]nucleus[0] == 0.0, IkBe\[LetterSpace]cytoplasm[0] == 0.0, IkBe\[LetterSpace]mRNA[0] == 0.0, IkBe\[LetterSpace]nucleus[0] == 0.0, NFkB\[LetterSpace]cytoplasm[0] == 0.0, NFkB\[LetterSpace]nucleus[0] == 0.0
};

rates = {
IKK\[LetterSpace]deg, deg\[LetterSpace]a, deg\[LetterSpace]a\[LetterSpace]n, deg\[LetterSpace]ai, deg\[LetterSpace]ain, deg\[LetterSpace]an, deg\[LetterSpace]an\[LetterSpace]n, deg\[LetterSpace]b, deg\[LetterSpace]b\[LetterSpace]n, deg\[LetterSpace]bi, deg\[LetterSpace]bin, deg\[LetterSpace]bn, deg\[LetterSpace]bn\[LetterSpace]n, deg\[LetterSpace]e, deg\[LetterSpace]e\[LetterSpace]n, deg\[LetterSpace]ei, deg\[LetterSpace]ein, deg\[LetterSpace]en, deg\[LetterSpace]en\[LetterSpace]n, int\[LetterSpace]2ain, int\[LetterSpace]2ani, int\[LetterSpace]2bin, int\[LetterSpace]2bni, int\[LetterSpace]2ein, int\[LetterSpace]2eni, int\[LetterSpace]ai, int\[LetterSpace]an, int\[LetterSpace]an\[LetterSpace]n, int\[LetterSpace]bi, int\[LetterSpace]bn, int\[LetterSpace]bn\[LetterSpace]n, int\[LetterSpace]ei, int\[LetterSpace]en, int\[LetterSpace]en\[LetterSpace]n, itxn\[LetterSpace]a, loc\[LetterSpace]a, loc\[LetterSpace]an, loc\[LetterSpace]b, loc\[LetterSpace]bn, loc\[LetterSpace]e, loc\[LetterSpace]en, loc\[LetterSpace]n, mdeg\[LetterSpace]a, mdeg\[LetterSpace]b, mdeg\[LetterSpace]e, tsl\[LetterSpace]a, tsl\[LetterSpace]b, tsl\[LetterSpace]e, txn\[LetterSpace]a, txn\[LetterSpace]b, txn\[LetterSpace]e
};

rateEquations = {
IKK\[LetterSpace]deg -> cytoplasm*IKK\[LetterSpace]deg\[LetterSpace]k\[LetterSpace]IKK\[LetterSpace]deg*IKK[t], deg\[LetterSpace]a -> cytoplasm*deg\[LetterSpace]a\[LetterSpace]deg1\[LetterSpace]c*IkBa\[LetterSpace]cytoplasm[t], deg\[LetterSpace]a\[LetterSpace]n -> deg\[LetterSpace]a\[LetterSpace]n\[LetterSpace]deg1\[LetterSpace]n*nucleus*IkBa\[LetterSpace]nucleus[t], deg\[LetterSpace]ai -> cytoplasm*deg\[LetterSpace]ai\[LetterSpace]r1*IkBaIKK[t], deg\[LetterSpace]ain -> cytoplasm*deg\[LetterSpace]ain\[LetterSpace]r4*IkBaIKKNFkB[t], deg\[LetterSpace]an -> cytoplasm*deg\[LetterSpace]an\[LetterSpace]deg4\[LetterSpace]c*IkBaNFkB\[LetterSpace]cytoplasm[t], deg\[LetterSpace]an\[LetterSpace]n -> deg\[LetterSpace]an\[LetterSpace]n\[LetterSpace]deg4\[LetterSpace]n*nucleus*IkBaNFkB\[LetterSpace]nucleus[t], deg\[LetterSpace]b -> cytoplasm*deg\[LetterSpace]b\[LetterSpace]deg2\[LetterSpace]c*IkBb\[LetterSpace]cytoplasm[t], deg\[LetterSpace]b\[LetterSpace]n -> deg\[LetterSpace]b\[LetterSpace]n\[LetterSpace]deg2\[LetterSpace]n*nucleus*IkBb\[LetterSpace]nucleus[t], deg\[LetterSpace]bi -> cytoplasm*deg\[LetterSpace]bi\[LetterSpace]r2*IkBbIKK[t], deg\[LetterSpace]bin -> cytoplasm*deg\[LetterSpace]bin\[LetterSpace]r5*IkBbIKKNFkB[t], deg\[LetterSpace]bn -> cytoplasm*deg\[LetterSpace]bn\[LetterSpace]deg5\[LetterSpace]c*IkBbNFkB\[LetterSpace]cytoplasm[t], deg\[LetterSpace]bn\[LetterSpace]n -> deg\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]deg5\[LetterSpace]n*nucleus*IkBbNFkB\[LetterSpace]nucleus[t], deg\[LetterSpace]e -> cytoplasm*deg\[LetterSpace]e\[LetterSpace]deg3\[LetterSpace]c*IkBe\[LetterSpace]cytoplasm[t], deg\[LetterSpace]e\[LetterSpace]n -> deg\[LetterSpace]e\[LetterSpace]n\[LetterSpace]deg3\[LetterSpace]n*nucleus*IkBe\[LetterSpace]nucleus[t], deg\[LetterSpace]ei -> cytoplasm*deg\[LetterSpace]ei\[LetterSpace]r3*IkBeIKK[t], deg\[LetterSpace]ein -> cytoplasm*deg\[LetterSpace]ein\[LetterSpace]r6*IkBeIKKNFkB[t], deg\[LetterSpace]en -> cytoplasm*deg\[LetterSpace]en\[LetterSpace]deg6\[LetterSpace]c*IkBeNFkB\[LetterSpace]cytoplasm[t], deg\[LetterSpace]en\[LetterSpace]n -> deg\[LetterSpace]en\[LetterSpace]n\[LetterSpace]deg6\[LetterSpace]n*nucleus*IkBeNFkB\[LetterSpace]nucleus[t], int\[LetterSpace]2ain -> cytoplasm*(-(int\[LetterSpace]2ain\[LetterSpace]d4\[LetterSpace]3*IkBaIKKNFkB[t]) + int\[LetterSpace]2ain\[LetterSpace]a4\[LetterSpace]3*IkBaIKK[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]2ani -> cytoplasm*(-(int\[LetterSpace]2ani\[LetterSpace]d1\[LetterSpace]2*IkBaIKKNFkB[t]) + int\[LetterSpace]2ani\[LetterSpace]a7*IkBaNFkB\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]2bin -> cytoplasm*(-(int\[LetterSpace]2bin\[LetterSpace]d5\[LetterSpace]3*IkBbIKKNFkB[t]) + int\[LetterSpace]2bin\[LetterSpace]a5\[LetterSpace]3*IkBbIKK[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]2bni -> cytoplasm*(-(int\[LetterSpace]2bni\[LetterSpace]d2\[LetterSpace]2*IkBbIKKNFkB[t]) + int\[LetterSpace]2bni\[LetterSpace]a8*IkBbNFkB\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]2ein -> cytoplasm*(-(int\[LetterSpace]2ein\[LetterSpace]d6\[LetterSpace]3*IkBeIKKNFkB[t]) + int\[LetterSpace]2ein\[LetterSpace]a6\[LetterSpace]3*IkBeIKK[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]2eni -> cytoplasm*(-(int\[LetterSpace]2eni\[LetterSpace]d3\[LetterSpace]2*IkBeIKKNFkB[t]) + int\[LetterSpace]2eni\[LetterSpace]a9*IkBeNFkB\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]ai -> cytoplasm*(-(int\[LetterSpace]ai\[LetterSpace]d1\[LetterSpace]1*IkBaIKK[t]) + int\[LetterSpace]ai\[LetterSpace]a1*IkBa\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]an -> cytoplasm*(-(int\[LetterSpace]an\[LetterSpace]d4\[LetterSpace]1*IkBaNFkB\[LetterSpace]cytoplasm[t]) + int\[LetterSpace]an\[LetterSpace]a4\[LetterSpace]1*IkBa\[LetterSpace]cytoplasm[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]an\[LetterSpace]n -> nucleus*(-(int\[LetterSpace]an\[LetterSpace]n\[LetterSpace]d4\[LetterSpace]2*IkBaNFkB\[LetterSpace]nucleus[t]) + int\[LetterSpace]an\[LetterSpace]n\[LetterSpace]a4\[LetterSpace]2*IkBa\[LetterSpace]nucleus[t]*NFkB\[LetterSpace]nucleus[t]), int\[LetterSpace]bi -> cytoplasm*(-(int\[LetterSpace]bi\[LetterSpace]d2\[LetterSpace]1*IkBbIKK[t]) + int\[LetterSpace]bi\[LetterSpace]a2*IkBb\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]bn -> cytoplasm*(-(int\[LetterSpace]bn\[LetterSpace]d5\[LetterSpace]1*IkBbNFkB\[LetterSpace]cytoplasm[t]) + int\[LetterSpace]bn\[LetterSpace]a5\[LetterSpace]1*IkBb\[LetterSpace]cytoplasm[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]bn\[LetterSpace]n -> nucleus*(-(int\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]d5\[LetterSpace]2*IkBbNFkB\[LetterSpace]nucleus[t]) + int\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]a5\[LetterSpace]2*IkBb\[LetterSpace]nucleus[t]*NFkB\[LetterSpace]nucleus[t]), int\[LetterSpace]ei -> cytoplasm*(-(int\[LetterSpace]ei\[LetterSpace]d3\[LetterSpace]1*IkBeIKK[t]) + int\[LetterSpace]ei\[LetterSpace]a3*IkBe\[LetterSpace]cytoplasm[t]*IKK[t]), int\[LetterSpace]en -> cytoplasm*(-(int\[LetterSpace]en\[LetterSpace]d6\[LetterSpace]1*IkBeNFkB\[LetterSpace]cytoplasm[t]) + int\[LetterSpace]en\[LetterSpace]a6\[LetterSpace]1*IkBe\[LetterSpace]cytoplasm[t]*NFkB\[LetterSpace]cytoplasm[t]), int\[LetterSpace]en\[LetterSpace]n -> nucleus*(-(int\[LetterSpace]en\[LetterSpace]n\[LetterSpace]d6\[LetterSpace]2*IkBeNFkB\[LetterSpace]nucleus[t]) + int\[LetterSpace]en\[LetterSpace]n\[LetterSpace]a6\[LetterSpace]2*IkBe\[LetterSpace]nucleus[t]*NFkB\[LetterSpace]nucleus[t]), itxn\[LetterSpace]a -> itxn\[LetterSpace]a\[LetterSpace]tr2a\[LetterSpace]i*nucleus*NFkB\[LetterSpace]nucleus[t]^2, loc\[LetterSpace]a -> cytoplasm*loc\[LetterSpace]a\[LetterSpace]tp1a*IkBa\[LetterSpace]cytoplasm[t] - loc\[LetterSpace]a\[LetterSpace]tp2a*nucleus*IkBa\[LetterSpace]nucleus[t], loc\[LetterSpace]an -> loc\[LetterSpace]an\[LetterSpace]k2\[LetterSpace]a*nucleus*IkBaNFkB\[LetterSpace]nucleus[t], loc\[LetterSpace]b -> cytoplasm*loc\[LetterSpace]b\[LetterSpace]tp1b*IkBb\[LetterSpace]cytoplasm[t] - loc\[LetterSpace]b\[LetterSpace]tp2b*nucleus*IkBb\[LetterSpace]nucleus[t], loc\[LetterSpace]bn -> loc\[LetterSpace]bn\[LetterSpace]k2\[LetterSpace]b*nucleus*IkBbNFkB\[LetterSpace]nucleus[t], loc\[LetterSpace]e -> cytoplasm*loc\[LetterSpace]e\[LetterSpace]tp1e*IkBe\[LetterSpace]cytoplasm[t] - loc\[LetterSpace]e\[LetterSpace]tp2e*nucleus*IkBe\[LetterSpace]nucleus[t], loc\[LetterSpace]en -> loc\[LetterSpace]en\[LetterSpace]k2\[LetterSpace]e*nucleus*IkBeNFkB\[LetterSpace]nucleus[t], loc\[LetterSpace]n -> cytoplasm*loc\[LetterSpace]n\[LetterSpace]k1\[LetterSpace]2*NFkB\[LetterSpace]cytoplasm[t] - loc\[LetterSpace]n\[LetterSpace]k1\[LetterSpace]1*nucleus*NFkB\[LetterSpace]nucleus[t], mdeg\[LetterSpace]a -> mdeg\[LetterSpace]a\[LetterSpace]tr3a*nucleus*IkBa\[LetterSpace]mRNA[t], mdeg\[LetterSpace]b -> mdeg\[LetterSpace]b\[LetterSpace]tr3b*nucleus*IkBb\[LetterSpace]mRNA[t], mdeg\[LetterSpace]e -> mdeg\[LetterSpace]e\[LetterSpace]tr3e*nucleus*IkBe\[LetterSpace]mRNA[t], tsl\[LetterSpace]a -> nucleus*tsl\[LetterSpace]a\[LetterSpace]tr1a*IkBa\[LetterSpace]mRNA[t], tsl\[LetterSpace]b -> nucleus*tsl\[LetterSpace]b\[LetterSpace]tr1b*IkBb\[LetterSpace]mRNA[t], tsl\[LetterSpace]e -> nucleus*tsl\[LetterSpace]e\[LetterSpace]tr1e*IkBe\[LetterSpace]mRNA[t], txn\[LetterSpace]a -> nucleus*txn\[LetterSpace]a\[LetterSpace]tr2a, txn\[LetterSpace]b -> nucleus*txn\[LetterSpace]b\[LetterSpace]tr2b, txn\[LetterSpace]e -> nucleus*txn\[LetterSpace]e\[LetterSpace]tr2e
};

parameters = {
txn\[LetterSpace]a\[LetterSpace]tr2a -> 0.0001848, mdeg\[LetterSpace]a\[LetterSpace]tr3a -> 0.0168, tsl\[LetterSpace]a\[LetterSpace]tr1a -> 0.2448, int\[LetterSpace]ai\[LetterSpace]a1 -> 1.35, int\[LetterSpace]ai\[LetterSpace]d1\[LetterSpace]1 -> 0.075, int\[LetterSpace]an\[LetterSpace]a4\[LetterSpace]1 -> 30.0, int\[LetterSpace]an\[LetterSpace]d4\[LetterSpace]1 -> 6*^-05, int\[LetterSpace]an\[LetterSpace]n\[LetterSpace]a4\[LetterSpace]2 -> 30.0, int\[LetterSpace]an\[LetterSpace]n\[LetterSpace]d4\[LetterSpace]2 -> 6*^-05, int\[LetterSpace]2ani\[LetterSpace]a7 -> 11.1, int\[LetterSpace]2ani\[LetterSpace]d1\[LetterSpace]2 -> 0.075, int\[LetterSpace]2ain\[LetterSpace]a4\[LetterSpace]3 -> 30.0, int\[LetterSpace]2ain\[LetterSpace]d4\[LetterSpace]3 -> 6*^-05, deg\[LetterSpace]a\[LetterSpace]deg1\[LetterSpace]c -> 0.12, deg\[LetterSpace]a\[LetterSpace]n\[LetterSpace]deg1\[LetterSpace]n -> 0.12, deg\[LetterSpace]an\[LetterSpace]n\[LetterSpace]deg4\[LetterSpace]n -> 6*^-05, deg\[LetterSpace]an\[LetterSpace]deg4\[LetterSpace]c -> 6*^-05, deg\[LetterSpace]ai\[LetterSpace]r1 -> 0.072, deg\[LetterSpace]ain\[LetterSpace]r4 -> 0.36, loc\[LetterSpace]a\[LetterSpace]tp1a -> 0.018, loc\[LetterSpace]a\[LetterSpace]tp2a -> 0.012, loc\[LetterSpace]an\[LetterSpace]k2\[LetterSpace]a -> 0.828, loc\[LetterSpace]bn\[LetterSpace]k2\[LetterSpace]b -> 0.414, loc\[LetterSpace]b\[LetterSpace]tp1b -> 0.018, loc\[LetterSpace]b\[LetterSpace]tp2b -> 0.012, deg\[LetterSpace]bin\[LetterSpace]r5 -> 0.12, deg\[LetterSpace]bi\[LetterSpace]r2 -> 0.024, deg\[LetterSpace]bn\[LetterSpace]deg5\[LetterSpace]c -> 6*^-05, deg\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]deg5\[LetterSpace]n -> 6*^-05, deg\[LetterSpace]b\[LetterSpace]n\[LetterSpace]deg2\[LetterSpace]n -> 0.18, deg\[LetterSpace]b\[LetterSpace]deg2\[LetterSpace]c -> 0.18, int\[LetterSpace]2bin\[LetterSpace]a5\[LetterSpace]3 -> 30.0, int\[LetterSpace]2bin\[LetterSpace]d5\[LetterSpace]3 -> 6*^-05, int\[LetterSpace]2bni\[LetterSpace]a8 -> 2.88, int\[LetterSpace]2bni\[LetterSpace]d2\[LetterSpace]2 -> 0.105, int\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]a5\[LetterSpace]2 -> 30.0, int\[LetterSpace]bn\[LetterSpace]n\[LetterSpace]d5\[LetterSpace]2 -> 6*^-05, int\[LetterSpace]bn\[LetterSpace]a5\[LetterSpace]1 -> 30.0, int\[LetterSpace]bn\[LetterSpace]d5\[LetterSpace]1 -> 6*^-05, int\[LetterSpace]bi\[LetterSpace]a2 -> 0.36, int\[LetterSpace]bi\[LetterSpace]d2\[LetterSpace]1 -> 0.105, tsl\[LetterSpace]b\[LetterSpace]tr1b -> 0.2448, mdeg\[LetterSpace]b\[LetterSpace]tr3b -> 0.0168, txn\[LetterSpace]b\[LetterSpace]tr2b -> 4.272*^-05, loc\[LetterSpace]n\[LetterSpace]k1\[LetterSpace]2 -> 5.4, loc\[LetterSpace]n\[LetterSpace]k1\[LetterSpace]1 -> 0.0048, txn\[LetterSpace]e\[LetterSpace]tr2e -> 3.048*^-05, mdeg\[LetterSpace]e\[LetterSpace]tr3e -> 0.0168, tsl\[LetterSpace]e\[LetterSpace]tr1e -> 0.2448, int\[LetterSpace]ei\[LetterSpace]a3 -> 0.54, int\[LetterSpace]ei\[LetterSpace]d3\[LetterSpace]1 -> 0.105, int\[LetterSpace]en\[LetterSpace]a6\[LetterSpace]1 -> 30.0, int\[LetterSpace]en\[LetterSpace]d6\[LetterSpace]1 -> 6*^-05, int\[LetterSpace]en\[LetterSpace]n\[LetterSpace]a6\[LetterSpace]2 -> 30.0, int\[LetterSpace]en\[LetterSpace]n\[LetterSpace]d6\[LetterSpace]2 -> 6*^-05, int\[LetterSpace]2eni\[LetterSpace]a9 -> 4.2, int\[LetterSpace]2eni\[LetterSpace]d3\[LetterSpace]2 -> 0.105, int\[LetterSpace]2ein\[LetterSpace]a6\[LetterSpace]3 -> 30.0, int\[LetterSpace]2ein\[LetterSpace]d6\[LetterSpace]3 -> 6*^-05, deg\[LetterSpace]e\[LetterSpace]deg3\[LetterSpace]c -> 0.18, deg\[LetterSpace]e\[LetterSpace]n\[LetterSpace]deg3\[LetterSpace]n -> 0.18, deg\[LetterSpace]en\[LetterSpace]n\[LetterSpace]deg6\[LetterSpace]n -> 6*^-05, deg\[LetterSpace]en\[LetterSpace]deg6\[LetterSpace]c -> 6*^-05, deg\[LetterSpace]ei\[LetterSpace]r3 -> 0.036, deg\[LetterSpace]ein\[LetterSpace]r6 -> 0.18, loc\[LetterSpace]e\[LetterSpace]tp1e -> 0.018, loc\[LetterSpace]e\[LetterSpace]tp2e -> 0.012, loc\[LetterSpace]en\[LetterSpace]k2\[LetterSpace]e -> 0.414, IKK\[LetterSpace]deg\[LetterSpace]k\[LetterSpace]IKK\[LetterSpace]deg -> 0.0, itxn\[LetterSpace]a\[LetterSpace]tr2a\[LetterSpace]i -> 1.98, cytoplasm -> 1.0, nucleus -> 1.0
};

assignments = {
Total\[LetterSpace]IkBeps -> IkBeIKK[t] + IkBeIKKNFkB[t] + IkBeNFkB\[LetterSpace]cytoplasm[t] + IkBeNFkB\[LetterSpace]nucleus[t] + IkBe\[LetterSpace]cytoplasm[t] + IkBe\[LetterSpace]nucleus[t], Total\[LetterSpace]IkBbeta -> IkBbIKK[t] + IkBbIKKNFkB[t] + IkBbNFkB\[LetterSpace]cytoplasm[t] + IkBbNFkB\[LetterSpace]nucleus[t] + IkBb\[LetterSpace]cytoplasm[t] + IkBb\[LetterSpace]nucleus[t], Total\[LetterSpace]IkBalpha -> IkBaIKK[t] + IkBaIKKNFkB[t] + IkBaNFkB\[LetterSpace]cytoplasm[t] + IkBaNFkB\[LetterSpace]nucleus[t] + IkBa\[LetterSpace]cytoplasm[t] + IkBa\[LetterSpace]nucleus[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {
deg\[LetterSpace]ai->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]ai->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]ain->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]ain->"http://identifiers.org/go/GO:0051092", deg\[LetterSpace]ain->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]bi->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]bi->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]bin->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]bin->"http://identifiers.org/go/GO:0051092", deg\[LetterSpace]bin->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]ei->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]ei->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]ein->"http://identifiers.org/ec-code/2.7.11.10", deg\[LetterSpace]ein->"http://identifiers.org/go/GO:0008384", deg\[LetterSpace]ein->"http://identifiers.org/go/GO:0051092", txn\[LetterSpace]a->"http://identifiers.org/go/GO:0006351", txn\[LetterSpace]b->"http://identifiers.org/go/GO:0006351", txn\[LetterSpace]e->"http://identifiers.org/go/GO:0006351"
};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
IKK'[t] == 1.0*deg\[LetterSpace]ai +1.0*deg\[LetterSpace]ain +1.0*deg\[LetterSpace]bin +1.0*deg\[LetterSpace]bi +1.0*deg\[LetterSpace]ei +1.0*deg\[LetterSpace]ein -1.0*int\[LetterSpace]ai -1.0*int\[LetterSpace]2ani -1.0*int\[LetterSpace]2bni -1.0*int\[LetterSpace]bi -1.0*int\[LetterSpace]ei -1.0*int\[LetterSpace]2eni -1.0*IKK\[LetterSpace]deg, IkBaIKK'[t] == 1.0*int\[LetterSpace]ai -1.0*int\[LetterSpace]2ain -1.0*deg\[LetterSpace]ai, IkBaIKKNFkB'[t] == 1.0*int\[LetterSpace]2ani +1.0*int\[LetterSpace]2ain -1.0*deg\[LetterSpace]ain, IkBaNFkB\[LetterSpace]cytoplasm'[t] == 1.0*int\[LetterSpace]an +1.0*loc\[LetterSpace]an -1.0*int\[LetterSpace]2ani -1.0*deg\[LetterSpace]an, IkBaNFkB\[LetterSpace]nucleus'[t] == 1.0*int\[LetterSpace]an\[LetterSpace]n -1.0*deg\[LetterSpace]an\[LetterSpace]n -1.0*loc\[LetterSpace]an, IkBa\[LetterSpace]cytoplasm'[t] == 1.0*tsl\[LetterSpace]a -1.0*int\[LetterSpace]ai -1.0*int\[LetterSpace]an -1.0*deg\[LetterSpace]a -1.0*loc\[LetterSpace]a, IkBa\[LetterSpace]mRNA'[t] == 1.0*txn\[LetterSpace]a +1.0*tsl\[LetterSpace]a +1.0*itxn\[LetterSpace]a -1.0*mdeg\[LetterSpace]a -1.0*tsl\[LetterSpace]a, IkBa\[LetterSpace]nucleus'[t] == 1.0*loc\[LetterSpace]a -1.0*int\[LetterSpace]an\[LetterSpace]n -1.0*deg\[LetterSpace]a\[LetterSpace]n, IkBbIKK'[t] == 1.0*int\[LetterSpace]bi -1.0*deg\[LetterSpace]bi -1.0*int\[LetterSpace]2bin, IkBbIKKNFkB'[t] == 1.0*int\[LetterSpace]2bin +1.0*int\[LetterSpace]2bni -1.0*deg\[LetterSpace]bin, IkBbNFkB\[LetterSpace]cytoplasm'[t] == 1.0*loc\[LetterSpace]bn +1.0*int\[LetterSpace]bn -1.0*deg\[LetterSpace]bn -1.0*int\[LetterSpace]2bni, IkBbNFkB\[LetterSpace]nucleus'[t] == 1.0*int\[LetterSpace]bn\[LetterSpace]n -1.0*loc\[LetterSpace]bn -1.0*deg\[LetterSpace]bn\[LetterSpace]n, IkBb\[LetterSpace]cytoplasm'[t] == 1.0*tsl\[LetterSpace]b -1.0*loc\[LetterSpace]b -1.0*deg\[LetterSpace]b -1.0*int\[LetterSpace]bn -1.0*int\[LetterSpace]bi, IkBb\[LetterSpace]mRNA'[t] == 1.0*tsl\[LetterSpace]b +1.0*txn\[LetterSpace]b -1.0*tsl\[LetterSpace]b -1.0*mdeg\[LetterSpace]b, IkBb\[LetterSpace]nucleus'[t] == 1.0*loc\[LetterSpace]b -1.0*deg\[LetterSpace]b\[LetterSpace]n -1.0*int\[LetterSpace]bn\[LetterSpace]n, IkBeIKK'[t] == 1.0*int\[LetterSpace]ei -1.0*int\[LetterSpace]2ein -1.0*deg\[LetterSpace]ei, IkBeIKKNFkB'[t] == 1.0*int\[LetterSpace]2eni +1.0*int\[LetterSpace]2ein -1.0*deg\[LetterSpace]ein, IkBeNFkB\[LetterSpace]cytoplasm'[t] == 1.0*int\[LetterSpace]en +1.0*loc\[LetterSpace]en -1.0*int\[LetterSpace]2eni -1.0*deg\[LetterSpace]en, IkBeNFkB\[LetterSpace]nucleus'[t] == 1.0*int\[LetterSpace]en\[LetterSpace]n -1.0*deg\[LetterSpace]en\[LetterSpace]n -1.0*loc\[LetterSpace]en, IkBe\[LetterSpace]cytoplasm'[t] == 1.0*tsl\[LetterSpace]e -1.0*int\[LetterSpace]ei -1.0*int\[LetterSpace]en -1.0*deg\[LetterSpace]e -1.0*loc\[LetterSpace]e, IkBe\[LetterSpace]mRNA'[t] == 1.0*txn\[LetterSpace]e +1.0*tsl\[LetterSpace]e -1.0*mdeg\[LetterSpace]e -1.0*tsl\[LetterSpace]e, IkBe\[LetterSpace]nucleus'[t] == 1.0*loc\[LetterSpace]e -1.0*int\[LetterSpace]en\[LetterSpace]n -1.0*deg\[LetterSpace]e\[LetterSpace]n, NFkB\[LetterSpace]cytoplasm'[t] == 1.0*deg\[LetterSpace]an +1.0*deg\[LetterSpace]ain +1.0*deg\[LetterSpace]bin +1.0*deg\[LetterSpace]bn +1.0*deg\[LetterSpace]en +1.0*deg\[LetterSpace]ein -1.0*int\[LetterSpace]an -1.0*int\[LetterSpace]2ain -1.0*int\[LetterSpace]2bin -1.0*int\[LetterSpace]bn -1.0*loc\[LetterSpace]n -1.0*int\[LetterSpace]en -1.0*int\[LetterSpace]2ein, NFkB\[LetterSpace]nucleus'[t] == 1.0*loc\[LetterSpace]n +1.0*deg\[LetterSpace]an\[LetterSpace]n +1.0*deg\[LetterSpace]bn\[LetterSpace]n +1.0*deg\[LetterSpace]en\[LetterSpace]n -1.0*int\[LetterSpace]an\[LetterSpace]n -1.0*int\[LetterSpace]bn\[LetterSpace]n -1.0*int\[LetterSpace]en\[LetterSpace]n
};
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]}]