(* 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 = { A1[t], A1\[LetterSpace]2[t], A2[t], A2\[LetterSpace]2[t], A3[t], A3\[LetterSpace]2[t], CREB\[LetterSpace]n[t], DUSP\[LetterSpace]c[t], DUSP\[LetterSpace]n[t], DUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[t], DUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t], DUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t], DUSPmRNA[t], ERK\[LetterSpace]c[t], ERK\[LetterSpace]n[t], Elk1\[LetterSpace]n[t], F[t], FOSn[t], FOSn\[LetterSpace]2[t], FmRNA[t], Fn[t], Kin[t], Kin\[LetterSpace]2[t], MEK[t], PreDUSPmRNA[t], PreFOSmRNA[t], PreFmRNA[t], RSK\[LetterSpace]c[t], RsD[t], RsT[t], c\[LetterSpace]FOS\[LetterSpace]c[t], c\[LetterSpace]FOSmRNA[t], pCREB\[LetterSpace]n[t], pDUSP\[LetterSpace]c[t], pDUSP\[LetterSpace]n[t], pDUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[t], pDUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t], pDUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t], pERK\[LetterSpace]c[t], pERK\[LetterSpace]n[t], pElk1\[LetterSpace]n[t], pMEK[t], pRSK\[LetterSpace]c[t], pRSK\[LetterSpace]n[t], pc\[LetterSpace]FOS\[LetterSpace]c[t], ppERK\[LetterSpace]c[t], ppERK\[LetterSpace]n[t] }; initialValues = { A1[0] == 182.3524, A1\[LetterSpace]2[0] == 0.0, A2[0] == 25.38702, A2\[LetterSpace]2[0] == 0.0, A3[0] == 13.09262, A3\[LetterSpace]2[0] == 0.0, CREB\[LetterSpace]n[0] == 1000.0, DUSP\[LetterSpace]c[0] == 0.0, DUSP\[LetterSpace]n[0] == 0.0, DUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[0] == 0.0, DUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[0] == 0.0, DUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[0] == 0.0, DUSPmRNA[0] == 0.0, ERK\[LetterSpace]c[0] == 570.4179, ERK\[LetterSpace]n[0] == 1624.9, Elk1\[LetterSpace]n[0] == 1510.0, F[0] == 0.0, FOSn[0] == 0.0, FOSn\[LetterSpace]2[0] == 0.0, FmRNA[0] == 0.0, Fn[0] == 0.0, Kin[0] == 82.66574, Kin\[LetterSpace]2[0] == 0.0, MEK[0] == 637.3212, PreDUSPmRNA[0] == 0.0, PreFOSmRNA[0] == 0.0, PreFmRNA[0] == 0.0, RSK\[LetterSpace]c[0] == 353.0, RsD[0] == 247.4035, RsT[0] == 0.0, c\[LetterSpace]FOS\[LetterSpace]c[0] == 0.0, c\[LetterSpace]FOSmRNA[0] == 0.0, pCREB\[LetterSpace]n[0] == 0.0, pDUSP\[LetterSpace]c[0] == 0.0, pDUSP\[LetterSpace]n[0] == 0.0, pDUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[0] == 0.0, pDUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[0] == 0.0, pDUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[0] == 0.0, pERK\[LetterSpace]c[0] == 0.0, pERK\[LetterSpace]n[0] == 0.0, pElk1\[LetterSpace]n[0] == 0.0, pMEK[0] == 0.0, pRSK\[LetterSpace]c[0] == 0.0, pRSK\[LetterSpace]n[0] == 0.0, pc\[LetterSpace]FOS\[LetterSpace]c[0] == 0.0, ppERK\[LetterSpace]c[0] == 0.0, ppERK\[LetterSpace]n[0] == 0.0 }; rates = { re106, re107, re108, re109, re110, re111, re112, re113, re115, re116, re117, re118, re119, re85, re86, re87, re88, re89, re90, re92, re94, re95, re96, re99, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]14, reaction\[LetterSpace]15, reaction\[LetterSpace]17, reaction\[LetterSpace]21, reaction\[LetterSpace]22, reaction\[LetterSpace]23, reaction\[LetterSpace]24, reaction\[LetterSpace]25, reaction\[LetterSpace]26, reaction\[LetterSpace]27, reaction\[LetterSpace]28, reaction\[LetterSpace]29, reaction\[LetterSpace]30, reaction\[LetterSpace]31, reaction\[LetterSpace]32, reaction\[LetterSpace]33, reaction\[LetterSpace]34, reaction\[LetterSpace]35, reaction\[LetterSpace]36, reaction\[LetterSpace]41, reaction\[LetterSpace]42, reaction\[LetterSpace]45, reaction\[LetterSpace]46, reaction\[LetterSpace]47, reaction\[LetterSpace]50, reaction\[LetterSpace]52, reaction\[LetterSpace]53, reaction\[LetterSpace]54, reaction\[LetterSpace]55, reaction\[LetterSpace]56, reaction\[LetterSpace]57, reaction\[LetterSpace]58, reaction\[LetterSpace]59, reaction\[LetterSpace]6, reaction\[LetterSpace]60, reaction\[LetterSpace]61, reaction\[LetterSpace]62, reaction\[LetterSpace]67, reaction\[LetterSpace]68, reaction\[LetterSpace]69, reaction\[LetterSpace]70, reaction\[LetterSpace]71, reaction\[LetterSpace]72, reaction\[LetterSpace]73, reaction\[LetterSpace]74, reaction\[LetterSpace]75, reaction\[LetterSpace]76, reaction\[LetterSpace]77, reaction\[LetterSpace]78, reaction\[LetterSpace]79, reaction\[LetterSpace]80, reaction\[LetterSpace]81 }; rateEquations = { re106 -> function\[LetterSpace]4\[LetterSpace]v007[KexERK, KimERK, Vc, Vn, ERK\[LetterSpace]c[t], ERK\[LetterSpace]n[t]], re107 -> function\[LetterSpace]4\[LetterSpace]v008[KexERKP, KimERKP, Vc, Vn, pERK\[LetterSpace]c[t], pERK\[LetterSpace]n[t]], re108 -> function\[LetterSpace]4\[LetterSpace]v009[KexERKPP, KimERKPP, Vc, Vn, ppERK\[LetterSpace]c[t], ppERK\[LetterSpace]n[t]], re109 -> nucleus*function\[LetterSpace]4\[LetterSpace]v010[K10, V10, n10, ppERK\[LetterSpace]n[t]], re110 -> function\[LetterSpace]4\[LetterSpace]v011[Vn, p11, PreDUSPmRNA[t]], re111 -> function\[LetterSpace]4\[LetterSpace]v018[KexDUSP, KimDUSP, Vc, Vn, DUSP\[LetterSpace]n[t], DUSP\[LetterSpace]c[t]], re112 -> function\[LetterSpace]4\[LetterSpace]v019[KexDUSPP, KimDUSPP, Vc, Vn, pDUSP\[LetterSpace]c[t], pDUSP\[LetterSpace]n[t]], re113 -> function\[LetterSpace]4\[LetterSpace]v026[KexRSKP, KimRSKP, Vc, Vn, pRSK\[LetterSpace]n[t], pRSK\[LetterSpace]c[t]], re115 -> function\[LetterSpace]4\[LetterSpace]v032[Vn, p32, PreFOSmRNA[t]], re116 -> function\[LetterSpace]4\[LetterSpace]v040[KexFOS, KimFOS, Vc, Vn, c\[LetterSpace]FOS\[LetterSpace]c[t], FOSn[t]], re117 -> function\[LetterSpace]4\[LetterSpace]v041[KexFOSP, KimFOSP, Vc, Vn, pc\[LetterSpace]FOS\[LetterSpace]c[t], FOSn\[LetterSpace]2[t]], re118 -> function\[LetterSpace]4\[LetterSpace]v058[Vn, p58, PreFmRNA[t]], re119 -> function\[LetterSpace]4\[LetterSpace]v062[KexF, KimF, Vc, Vn, F[t], Fn[t]], re85 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v003[K3, K4, V3, pERK\[LetterSpace]c[t], ppERK\[LetterSpace]c[t]], re86 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v004[K3, K4, V4, pERK\[LetterSpace]c[t], ppERK\[LetterSpace]c[t]], re87 -> nucleus*function\[LetterSpace]4\[LetterSpace]v005[K5, K6, V5, pERK\[LetterSpace]n[t], ppERK\[LetterSpace]n[t]], re88 -> nucleus*function\[LetterSpace]4\[LetterSpace]v006[K5, K6, V6, pERK\[LetterSpace]n[t], ppERK\[LetterSpace]n[t]], re89 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v013[p13, DUSPmRNA[t]], re90 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v034[p34, c\[LetterSpace]FOSmRNA[t]], re92 -> nucleus*function\[LetterSpace]4\[LetterSpace]v057[K57, V57, n57, FOSn\[LetterSpace]2[t]], re94 -> cytoplasm*p59*FmRNA[t], re95 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v060[p60, FmRNA[t]], re96 -> cytoplasm*p61*F[t], re99 -> nucleus*function\[LetterSpace]4\[LetterSpace]v063[cytoplasm, nucleus, p63, Fn[t]], reaction\[LetterSpace]12 -> nucleus*function\[LetterSpace]4\[LetterSpace]v020[K20, V20, ppERK\[LetterSpace]n[t], DUSP\[LetterSpace]n[t]], reaction\[LetterSpace]13 -> nucleus*function\[LetterSpace]4\[LetterSpace]v021[K21, V21, pDUSP\[LetterSpace]n[t]], reaction\[LetterSpace]14 -> nucleus*p22*DUSP\[LetterSpace]n[t], reaction\[LetterSpace]15 -> nucleus*p23*pDUSP\[LetterSpace]n[t], reaction\[LetterSpace]17 -> cytoplasm*p12*DUSPmRNA[t], reaction\[LetterSpace]21 -> nucleus*(-(m47*DUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t]) + p47*DUSP\[LetterSpace]n[t]*ppERK\[LetterSpace]n[t]), reaction\[LetterSpace]22 -> nucleus*p48*DUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t], reaction\[LetterSpace]23 -> nucleus*(-(m49*DUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t]) + p49*DUSP\[LetterSpace]n[t]*pERK\[LetterSpace]n[t]), reaction\[LetterSpace]24 -> nucleus*p50*DUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t], reaction\[LetterSpace]25 -> nucleus*(-(m52*pDUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t]) + p52*pDUSP\[LetterSpace]n[t]*ppERK\[LetterSpace]n[t]), reaction\[LetterSpace]26 -> nucleus*p53*pDUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n[t], reaction\[LetterSpace]27 -> nucleus*(-(m54*pDUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t]) + p54*pDUSP\[LetterSpace]n[t]*pERK\[LetterSpace]n[t]), reaction\[LetterSpace]28 -> nucleus*p55*pDUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n[t], reaction\[LetterSpace]29 -> nucleus*(-(m51*DUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[t]) + p51*DUSP\[LetterSpace]n[t]*ERK\[LetterSpace]n[t]), reaction\[LetterSpace]30 -> nucleus*(p56*ERK\[LetterSpace]n[t]*pDUSP\[LetterSpace]n[t] - m56*pDUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n[t]), reaction\[LetterSpace]31 -> cytoplasm*EGFm*function\[LetterSpace]4\[LetterSpace]v101[K101, V101, A1[t], EGF], reaction\[LetterSpace]32 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v001[Fct, K1, K2, V1, ERK\[LetterSpace]c[t], pERK\[LetterSpace]c[t], pMEK[t]], reaction\[LetterSpace]33 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v014[K14, V14, ppERK\[LetterSpace]c[t], DUSP\[LetterSpace]c[t]], reaction\[LetterSpace]34 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v015[K15, V15, pDUSP\[LetterSpace]c[t]], reaction\[LetterSpace]35 -> cytoplasm*p16*DUSP\[LetterSpace]c[t], reaction\[LetterSpace]36 -> cytoplasm*p17*pDUSP\[LetterSpace]c[t], reaction\[LetterSpace]41 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v024[K24, V24, ppERK\[LetterSpace]c[t], RSK\[LetterSpace]c[t]], reaction\[LetterSpace]42 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v025[K25, V25, pRSK\[LetterSpace]c[t]], reaction\[LetterSpace]45 -> nucleus*function\[LetterSpace]4\[LetterSpace]v027[K27, V27, CREB\[LetterSpace]n[t], pRSK\[LetterSpace]n[t]], reaction\[LetterSpace]46 -> nucleus*function\[LetterSpace]4\[LetterSpace]v028[K28, V28, pCREB\[LetterSpace]n[t]], reaction\[LetterSpace]47 -> nucleus*function\[LetterSpace]4\[LetterSpace]v031[K31, KF31, V31, n31, nF31, Fn[t], pCREB\[LetterSpace]n[t], pElk1\[LetterSpace]n[t]], reaction\[LetterSpace]50 -> cytoplasm*p33*c\[LetterSpace]FOSmRNA[t], reaction\[LetterSpace]52 -> nucleus*function\[LetterSpace]4\[LetterSpace]v042[K42, V42, ppERK\[LetterSpace]n[t], FOSn[t]], reaction\[LetterSpace]53 -> nucleus*function\[LetterSpace]4\[LetterSpace]v043[K43, V43, FOSn[t], pRSK\[LetterSpace]n[t]], reaction\[LetterSpace]54 -> nucleus*function\[LetterSpace]4\[LetterSpace]v044[K44, V44, FOSn\[LetterSpace]2[t]], reaction\[LetterSpace]55 -> nucleus*p45*FOSn[t], reaction\[LetterSpace]56 -> nucleus*p46*FOSn\[LetterSpace]2[t], reaction\[LetterSpace]57 -> nucleus*function\[LetterSpace]4\[LetterSpace]v029[K29, V29, ppERK\[LetterSpace]n[t], Elk1\[LetterSpace]n[t]], reaction\[LetterSpace]58 -> nucleus*function\[LetterSpace]4\[LetterSpace]v030[K30, V30, pElk1\[LetterSpace]n[t]], reaction\[LetterSpace]59 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v035[K35, V35, c\[LetterSpace]FOS\[LetterSpace]c[t], ppERK\[LetterSpace]c[t]], reaction\[LetterSpace]6 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v002[Fct, K1, K2, V2, ERK\[LetterSpace]c[t], pERK\[LetterSpace]c[t], pMEK[t]], reaction\[LetterSpace]60 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v037[K37, V37, pc\[LetterSpace]FOS\[LetterSpace]c[t]], reaction\[LetterSpace]61 -> cytoplasm*p38*c\[LetterSpace]FOS\[LetterSpace]c[t], reaction\[LetterSpace]62 -> cytoplasm*p39*pc\[LetterSpace]FOS\[LetterSpace]c[t], reaction\[LetterSpace]67 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v036[K36, V36, c\[LetterSpace]FOS\[LetterSpace]c[t], pRSK\[LetterSpace]c[t]], reaction\[LetterSpace]68 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v102[K102, V102, A1\[LetterSpace]2[t]], reaction\[LetterSpace]69 -> cytoplasm*HRGm*function\[LetterSpace]4\[LetterSpace]v103[K103, V103, A2[t], HRG], reaction\[LetterSpace]70 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v104[K104, V104, A2\[LetterSpace]2[t]], reaction\[LetterSpace]71 -> cytoplasm*EGFm*function\[LetterSpace]4\[LetterSpace]v105[K105, V105, EGF, RsD[t]], reaction\[LetterSpace]72 -> cytoplasm*HRGm*function\[LetterSpace]4\[LetterSpace]v106[K106, V106, HRG, RsD[t]], reaction\[LetterSpace]73 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v107[K107, V107, A1\[LetterSpace]2[t], RsT[t]], reaction\[LetterSpace]74 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v108[K108, V108, A2\[LetterSpace]2[t], RsT[t]], reaction\[LetterSpace]75 -> cytoplasm*HRGm*function\[LetterSpace]4\[LetterSpace]v109[K109, V109, HRG, A3[t]], reaction\[LetterSpace]76 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v110[K110, V110, A3\[LetterSpace]2[t]], reaction\[LetterSpace]77 -> cytoplasm*HRGm*function\[LetterSpace]4\[LetterSpace]v111[K111, V111, HRG, Kin[t]], reaction\[LetterSpace]78 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v112[K112, V112, A3\[LetterSpace]2[t], Kin\[LetterSpace]2[t]], reaction\[LetterSpace]79 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v113[K113, V113, RsT[t], MEK[t]], reaction\[LetterSpace]80 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v114[K114, V114, Kin\[LetterSpace]2[t], MEK[t]], reaction\[LetterSpace]81 -> cytoplasm*function\[LetterSpace]4\[LetterSpace]v115[K115, V115, pMEK[t]] }; parameters = { EGFm -> 1.0, Fct -> 0.7485, HRGm -> 1.0, K1 -> 307.041525298866, K10 -> 169.0473748, K101 -> 3475.168, K102 -> 237.2001, K103 -> 1334.132, K104 -> 4046.71, K105 -> 1.027895, K106 -> 606.871, K107 -> 424.6884, K108 -> 11.5048, K109 -> 7424.816, K110 -> 425.5268, K111 -> 858.3423, K112 -> 4665.217, K113 -> 20.50809, K114 -> 7.774197, K115 -> 2122.045, K14 -> 34180.48, K15 -> 0.001172165, K2 -> 350.0, K20 -> 735598.6967, K21 -> 387.8377182, K24 -> 29516.06587, K25 -> 0.913939859, K27 -> 441.5834425, K28 -> 14.99180922, K29 -> 21312.69109, K3 -> 160.0, K30 -> 15.04396629, K31 -> 185.9760682, K35 -> 8562.744184, K36 -> 528.552427, K37 -> 0.070379236, K4 -> 60.0, K42 -> 3992.061328, K43 -> 1157.116021, K44 -> 0.051168202, K5 -> 29.94073716, K57 -> 0.637490056, K6 -> 29.9407371620698, KF31 -> 0.013844393, KexDUSP -> 0.070467899, KexDUSPP -> 0.070467899, KexERK -> 0.018, KexERKP -> 0.018, KexERKPP -> 0.013, KexF -> 0.396950616, KexFOS -> 0.133249762, KexFOSP -> 0.133249762, KexRSKP -> 0.129803956, KimDUSP -> 0.024269764, KimDUSPP -> 0.024269764, KimERK -> 0.012, KimERKP -> 0.012, KimERKPP -> 0.011, KimF -> 0.019898797, KimFOS -> 0.54528521, KimFOSP -> 0.54528521, KimRSKP -> 0.025925065, V1 -> 0.342848369838443, V10 -> 29.24109258, V101 -> 0.01807448, V102 -> 0.09858154, V103 -> 0.3573399, V104 -> 4.635749, V105 -> 0.05393704, V106 -> 0.109304, V107 -> 5.291093, V108 -> 0.03436149, V109 -> 0.1374307, V110 -> 0.08258693, V111 -> 0.02487469, V112 -> 0.8850982, V113 -> 0.05377297, V114 -> 0.03957055, V115 -> 13.74244, V14 -> 5.636949216, V15 -> 2.992346912, V2 -> 0.22, V20 -> 0.157678678, V21 -> 0.005648117, V24 -> 0.550346114, V25 -> 10.09063736, V27 -> 19.23118154, V28 -> 6.574759504, V29 -> 0.518529841, V3 -> 0.72, V30 -> 13.79479021, V31 -> 0.655214248, V35 -> 8.907637012, V36 -> 0.000597315, V37 -> 1.745848179, V4 -> 0.648, V42 -> 0.909968714, V43 -> 0.076717457, V44 -> 0.078344305, V5 -> 19.49872346, V57 -> 1.026834758, V6 -> 19.4987234631759, Vc -> 940.0, Vn -> 220.0, m47 -> 15.80783969, m49 -> 2.335459127, m51 -> 9.544308421, m52 -> 15.80783969, m54 -> 2.335459127, m56 -> 9.544308421, n10 -> 3.970849295, n31 -> 1.988003164, n57 -> 3.584464176, nF31 -> 2.800340453, p11 -> 0.000126129, p12 -> 0.007875765, p13 -> 0.001245747, p16 -> 0.000257, p17 -> 4.81*^-05, p22 -> 0.000257, p23 -> 4.81*^-05, p32 -> 0.003284434, p33 -> 0.000601234209304622, p34 -> 7.64816282169636*^-05, p38 -> 0.000257, p39 -> 4.81*^-05, p45 -> 0.000257, p46 -> 4.81*^-05, p47 -> 0.001670815, p48 -> 0.686020478, p49 -> 0.314470502, p50 -> 26.59483436, p51 -> 0.01646825, p52 -> 0.001670815, p53 -> 0.686020478, p54 -> 0.314470502, p55 -> 26.59483436, p56 -> 0.01646825, p58 -> 0.000270488, p59 -> 0.001443889, p60 -> 0.002448164, p61 -> 3.49860901414122*^-05, p63 -> 4.13466150826031*^-05, EGF -> 1.0, HRG -> 1.0, cytoplasm -> 940.0, default -> 1.0, nucleus -> 220.0 }; assignments = { function\[LetterSpace]4\[LetterSpace]v062[KexF_,KimF_,Vc_,Vn_,F_,Fn_] -> F*KimF*Vc - Fn*KexF*Vn, function\[LetterSpace]4\[LetterSpace]v109[K109_,V109_,HRG_,A3_] -> (A3*HRG*V109)/(A3 + K109), function\[LetterSpace]4\[LetterSpace]v043[K43_,V43_,FOSn_,pRSK\[LetterSpace]n_] -> (FOSn*pRSK\[LetterSpace]n*V43)/(FOSn + K43), function\[LetterSpace]4\[LetterSpace]v104[K104_,V104_,A2\[LetterSpace]2_] -> (A2\[LetterSpace]2*V104)/(A2\[LetterSpace]2 + K104), function\[LetterSpace]4\[LetterSpace]v115[K115_,V115_,MEKc_] -> (MEKc*V115)/(K115 + MEKc), function\[LetterSpace]4\[LetterSpace]v035[K35_,V35_,c\[LetterSpace]FOS\[LetterSpace]c_,ppERK\[LetterSpace]c_] -> (c\[LetterSpace]FOS\[LetterSpace]c*ppERK\[LetterSpace]c*V35)/(c\[LetterSpace]FOS\[LetterSpace]c + K35), function\[LetterSpace]4\[LetterSpace]v041[KexFOSP_,KimFOSP_,Vc_,Vn_,pc\[LetterSpace]FOS\[LetterSpace]c_,FOSn\[LetterSpace]2_] -> KimFOSP*pc\[LetterSpace]FOS\[LetterSpace]c*Vc - FOSn\[LetterSpace]2*KexFOSP*Vn, function\[LetterSpace]4\[LetterSpace]v008[KexERKP_,KimERKP_,Vc_,Vn_,pERK\[LetterSpace]c_,pERK\[LetterSpace]n_] -> KimERKP*pERK\[LetterSpace]c*Vc - KexERKP*pERK\[LetterSpace]n*Vn, function\[LetterSpace]4\[LetterSpace]v015[K15_,V15_,pDUSP\[LetterSpace]c_] -> (pDUSP\[LetterSpace]c*V15)/(K15 + pDUSP\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v018[KexDUSP_,KimDUSP_,Vc_,Vn_,DUSP\[LetterSpace]n_,DUSP\[LetterSpace]c_] -> DUSP\[LetterSpace]c*KimDUSP*Vc - DUSP\[LetterSpace]n*KexDUSP*Vn, function\[LetterSpace]4\[LetterSpace]v003[K3_,K4_,V3_,pERK\[LetterSpace]c_,ppERK\[LetterSpace]c_] -> (pERK\[LetterSpace]c*V3)/(pERK\[LetterSpace]c + K3*(1 + ppERK\[LetterSpace]c/K4)), function\[LetterSpace]4\[LetterSpace]v005[K5_,K6_,V5_,pERK\[LetterSpace]n_,ppERK\[LetterSpace]n_] -> (pERK\[LetterSpace]n*V5)/(pERK\[LetterSpace]n + K5*(1 + ppERK\[LetterSpace]n/K6)), function\[LetterSpace]4\[LetterSpace]v006[K5_,K6_,V6_,pERK\[LetterSpace]n_,ppERK\[LetterSpace]n_] -> (ppERK\[LetterSpace]n*V6)/(K6*(1 + pERK\[LetterSpace]n/K5) + ppERK\[LetterSpace]n), function\[LetterSpace]4\[LetterSpace]v020[K20_,V20_,ppERK\[LetterSpace]n_,DUSP\[LetterSpace]n_] -> (DUSP\[LetterSpace]n*ppERK\[LetterSpace]n*V20)/(DUSP\[LetterSpace]n + K20), function\[LetterSpace]4\[LetterSpace]v021[K21_,V21_,pDUSP\[LetterSpace]n_] -> (pDUSP\[LetterSpace]n*V21)/(K21 + pDUSP\[LetterSpace]n), function\[LetterSpace]4\[LetterSpace]v024[K24_,V24_,ppERK\[LetterSpace]c_,RSK\[LetterSpace]c_] -> (ppERK\[LetterSpace]c*RSK\[LetterSpace]c*V24)/(K24 + RSK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v002[Fct_,K1_,K2_,V2_,ERK\[LetterSpace]c_,pERK\[LetterSpace]c_,MEKc_] -> (Fct*MEKc*pERK\[LetterSpace]c*V2)/((1 + ERK\[LetterSpace]c/K1)*K2 + pERK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v036[K36_,V36_,c\[LetterSpace]FOS\[LetterSpace]c_,pRSK\[LetterSpace]c_] -> (c\[LetterSpace]FOS\[LetterSpace]c*pRSK\[LetterSpace]c*V36)/(c\[LetterSpace]FOS\[LetterSpace]c + K36), function\[LetterSpace]4\[LetterSpace]v037[K37_,V37_,pc\[LetterSpace]FOS\[LetterSpace]c_] -> (pc\[LetterSpace]FOS\[LetterSpace]c*V37)/(K37 + pc\[LetterSpace]FOS\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v040[KexFOS_,KimFOS_,Vc_,Vn_,c\[LetterSpace]FOS\[LetterSpace]c_,FOSn_] -> c\[LetterSpace]FOS\[LetterSpace]c*KimFOS*Vc - FOSn*KexFOS*Vn, function\[LetterSpace]4\[LetterSpace]v107[K107_,V107_,A1\[LetterSpace]2_,RsT_] -> (A1\[LetterSpace]2*RsT*V107)/(K107 + RsT), function\[LetterSpace]4\[LetterSpace]v014[K14_,V14_,ppERK\[LetterSpace]c_,DUSP\[LetterSpace]c_] -> (DUSP\[LetterSpace]c*ppERK\[LetterSpace]c*V14)/(DUSP\[LetterSpace]c + K14), function\[LetterSpace]4\[LetterSpace]v110[K110_,V110_,A3\[LetterSpace]2_] -> (A3\[LetterSpace]2*V110)/(A3\[LetterSpace]2 + K110), function\[LetterSpace]4\[LetterSpace]v114[K114_,V114_,Kin\[LetterSpace]2_,pMEK\[LetterSpace]c_] -> (Kin\[LetterSpace]2*pMEK\[LetterSpace]c*V114)/(K114 + pMEK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v031[K31_,KF31_,V31_,n31_,nF31_,Fn_,pCREB\[LetterSpace]n_,pElk1\[LetterSpace]n_] -> ((pCREB\[LetterSpace]n*pElk1\[LetterSpace]n)^n31*V31)/(K31^n31 + (Fn/KF31)^nF31 + (pCREB\[LetterSpace]n*pElk1\[LetterSpace]n)^n31), function\[LetterSpace]4\[LetterSpace]v007[KexERK_,KimERK_,Vc_,Vn_,ERK\[LetterSpace]c_,ERK\[LetterSpace]n_] -> ERK\[LetterSpace]c*KimERK*Vc - ERK\[LetterSpace]n*KexERK*Vn, function\[LetterSpace]4\[LetterSpace]v004[K3_,K4_,V4_,pERK\[LetterSpace]c_,ppERK\[LetterSpace]c_] -> (ppERK\[LetterSpace]c*V4)/(K4*(1 + pERK\[LetterSpace]c/K3) + ppERK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v001[Fct_,K1_,K2_,V1_,ERK\[LetterSpace]c_,pERK\[LetterSpace]c_,MEKc_] -> (ERK\[LetterSpace]c*Fct*MEKc*V1)/(ERK\[LetterSpace]c + K1*(1 + pERK\[LetterSpace]c/K2)), function\[LetterSpace]4\[LetterSpace]v010[K10_,V10_,n10_,ppERK\[LetterSpace]n_] -> (ppERK\[LetterSpace]n^n10*V10)/(K10^n10 + ppERK\[LetterSpace]n^n10), function\[LetterSpace]4\[LetterSpace]v113[K113_,V113_,RsT_,pMEK\[LetterSpace]c_] -> (pMEK\[LetterSpace]c*RsT*V113)/(K113 + pMEK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v111[K111_,V111_,HRG_,Kin_] -> (HRG*Kin*V111)/(K111 + Kin), function\[LetterSpace]4\[LetterSpace]v112[K112_,V112_,A3\[LetterSpace]2_,Kin\[LetterSpace]2_] -> (A3\[LetterSpace]2*Kin\[LetterSpace]2*V112)/(K112 + Kin\[LetterSpace]2), function\[LetterSpace]4\[LetterSpace]v108[K108_,V108_,A2\[LetterSpace]2_,RsT_] -> (A2\[LetterSpace]2*RsT*V108)/(K108 + RsT), function\[LetterSpace]4\[LetterSpace]v106[K106_,V106_,HRG_,RsD_] -> (HRG*RsD*V106)/(K106 + RsD), function\[LetterSpace]4\[LetterSpace]v105[K105_,V105_,EGF_,RsD_] -> (EGF*RsD*V105)/(K105 + RsD), function\[LetterSpace]4\[LetterSpace]v103[K103_,V103_,A2_,HRG_] -> (A2*HRG*V103)/(A2 + K103), function\[LetterSpace]4\[LetterSpace]v102[K102_,V102_,A1\[LetterSpace]2_] -> (A1\[LetterSpace]2*V102)/(A1\[LetterSpace]2 + K102), function\[LetterSpace]4\[LetterSpace]v101[K101_,V101_,A1_,EGF_] -> (A1*EGF*V101)/(A1 + K101), function\[LetterSpace]4\[LetterSpace]v063[cytoplasm_,nucleus_,p63_,Fn_] -> (cytoplasm*Fn*p63)/nucleus, function\[LetterSpace]4\[LetterSpace]v060[p60_,FmRNA_] -> FmRNA*p60, function\[LetterSpace]4\[LetterSpace]v057[K57_,V57_,n57_,FOSn\[LetterSpace]2_] -> (FOSn\[LetterSpace]2^n57*V57)/(FOSn\[LetterSpace]2^n57 + K57^n57), function\[LetterSpace]4\[LetterSpace]v044[K44_,V44_,FOSn\[LetterSpace]2_] -> (FOSn\[LetterSpace]2*V44)/(FOSn\[LetterSpace]2 + K44), function\[LetterSpace]4\[LetterSpace]v042[K42_,V42_,ppERK\[LetterSpace]n_,FOSn_] -> (FOSn*ppERK\[LetterSpace]n*V42)/(FOSn + K42), function\[LetterSpace]4\[LetterSpace]v030[K30_,V30_,pElk1\[LetterSpace]n_] -> (pElk1\[LetterSpace]n*V30)/(K30 + pElk1\[LetterSpace]n), function\[LetterSpace]4\[LetterSpace]v028[K28_,V28_,pCREB\[LetterSpace]n_] -> (pCREB\[LetterSpace]n*V28)/(K28 + pCREB\[LetterSpace]n), function\[LetterSpace]4\[LetterSpace]v027[K27_,V27_,CREB\[LetterSpace]n_,pRSK\[LetterSpace]n_] -> (CREB\[LetterSpace]n*pRSK\[LetterSpace]n*V27)/(CREB\[LetterSpace]n + K27), function\[LetterSpace]4\[LetterSpace]v058[Vn_,p58_,PreFmRNA_] -> p58*PreFmRNA*Vn, function\[LetterSpace]4\[LetterSpace]v034[p34_,c\[LetterSpace]FOSmRNA_] -> c\[LetterSpace]FOSmRNA*p34, function\[LetterSpace]4\[LetterSpace]v026[KexRSKP_,KimRSKP_,Vc_,Vn_,pRSK\[LetterSpace]n_,pRSK\[LetterSpace]c_] -> KimRSKP*pRSK\[LetterSpace]c*Vc - KexRSKP*pRSK\[LetterSpace]n*Vn, function\[LetterSpace]4\[LetterSpace]v025[K25_,V25_,pRSK\[LetterSpace]c_] -> (pRSK\[LetterSpace]c*V25)/(K25 + pRSK\[LetterSpace]c), function\[LetterSpace]4\[LetterSpace]v032[Vn_,p32_,PreFOSmRNA_] -> p32*PreFOSmRNA*Vn, function\[LetterSpace]4\[LetterSpace]v013[p13_,DUSPmRNA_] -> DUSPmRNA*p13, function\[LetterSpace]4\[LetterSpace]v011[Vn_,p11_,PreDUSPmRNA_] -> p11*PreDUSPmRNA*Vn, function\[LetterSpace]4\[LetterSpace]v009[KexERKPP_,KimERKPP_,Vc_,Vn_,ppERK\[LetterSpace]c_,ppERK\[LetterSpace]n_] -> KimERKPP*ppERK\[LetterSpace]c*Vc - KexERKPP*ppERK\[LetterSpace]n*Vn, function\[LetterSpace]4\[LetterSpace]v029[K29_,V29_,ppERK\[LetterSpace]n_,Elk1\[LetterSpace]n_] -> (Elk1\[LetterSpace]n*ppERK\[LetterSpace]n*V29)/(Elk1\[LetterSpace]n + K29), function\[LetterSpace]4\[LetterSpace]v019[KexDUSPP_,KimDUSPP_,Vc_,Vn_,pDUSP\[LetterSpace]c_,pDUSP\[LetterSpace]n_] -> KimDUSPP*pDUSP\[LetterSpace]c*Vc - KexDUSPP*pDUSP\[LetterSpace]n*Vn }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { A1'[t] == 1.0*reaction\[LetterSpace]68 -1.0*reaction\[LetterSpace]31, A1\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]31 -1.0*reaction\[LetterSpace]68, A2'[t] == 1.0*reaction\[LetterSpace]70 -1.0*reaction\[LetterSpace]69, A2\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]69 -1.0*reaction\[LetterSpace]70, A3'[t] == 1.0*reaction\[LetterSpace]76 -1.0*reaction\[LetterSpace]75, A3\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]75 -1.0*reaction\[LetterSpace]76, CREB\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]45, DUSP\[LetterSpace]c'[t] == 1.0*re89 +1.0*reaction\[LetterSpace]34 -1.0*re111 -1.0*reaction\[LetterSpace]33 -1.0*reaction\[LetterSpace]35, DUSP\[LetterSpace]n'[t] == 1.0*re111 +1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]24 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]29, DUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]29 , DUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]24, DUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]22, DUSPmRNA'[t] == 1.0*re110 -1.0*reaction\[LetterSpace]17, ERK\[LetterSpace]c'[t] == 1.0*re85 -1.0*re106 -1.0*reaction\[LetterSpace]32, ERK\[LetterSpace]n'[t] == 1.0*re106 +1.0*re87 +1.0*reaction\[LetterSpace]24 +1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]30, Elk1\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]57, F'[t] == 1.0*re95 -1.0*re119 -1.0*re96, FOSn'[t] == 1.0*re116 +1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]55, FOSn\[LetterSpace]2'[t] == 1.0*re117 +1.0*reaction\[LetterSpace]52 +1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]56, FmRNA'[t] == 1.0*re118 -1.0*re94, Fn'[t] == 1.0*re119 -1.0*re99, Kin'[t] == 1.0*reaction\[LetterSpace]78 -1.0*reaction\[LetterSpace]77, Kin\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]77 -1.0*reaction\[LetterSpace]78, MEK'[t] == 1.0*reaction\[LetterSpace]81 -1.0*reaction\[LetterSpace]79 -1.0*reaction\[LetterSpace]80, PreDUSPmRNA'[t] == 1.0*re109 -1.0*re110, PreFOSmRNA'[t] == 1.0*reaction\[LetterSpace]47 -1.0*re115, PreFmRNA'[t] == 1.0*re92 -1.0*re118, RSK\[LetterSpace]c'[t] == 1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]41, RsD'[t] == 1.0*reaction\[LetterSpace]73 +1.0*reaction\[LetterSpace]74 -1.0*reaction\[LetterSpace]71 -1.0*reaction\[LetterSpace]72, RsT'[t] == 1.0*reaction\[LetterSpace]71 +1.0*reaction\[LetterSpace]72 -1.0*reaction\[LetterSpace]73 -1.0*reaction\[LetterSpace]74, c\[LetterSpace]FOS\[LetterSpace]c'[t] == 1.0*re90 +1.0*reaction\[LetterSpace]60 -1.0*re116 -1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]61 -1.0*reaction\[LetterSpace]67, c\[LetterSpace]FOSmRNA'[t] == 1.0*re115 -1.0*reaction\[LetterSpace]50, pCREB\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]46, pDUSP\[LetterSpace]c'[t] == 1.0*reaction\[LetterSpace]33 -1.0*re112 -1.0*reaction\[LetterSpace]34 -1.0*reaction\[LetterSpace]36, pDUSP\[LetterSpace]n'[t] == 1.0*re112 +1.0*reaction\[LetterSpace]12 +1.0*reaction\[LetterSpace]26 +1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]25 -1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]30, pDUSP\[LetterSpace]n\[LetterSpace]ERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]30 , pDUSP\[LetterSpace]n\[LetterSpace]pERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]28, pDUSP\[LetterSpace]n\[LetterSpace]ppERK\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]25 -1.0*reaction\[LetterSpace]26, pERK\[LetterSpace]c'[t] == 1.0*re86 +1.0*reaction\[LetterSpace]32 -1.0*re107 -1.0*re85 -1.0*reaction\[LetterSpace]6, pERK\[LetterSpace]n'[t] == 1.0*re107 +1.0*re88 +1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]26 -1.0*re87 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]27, pElk1\[LetterSpace]n'[t] == 1.0*reaction\[LetterSpace]57 -1.0*reaction\[LetterSpace]58, pMEK'[t] == 1.0*reaction\[LetterSpace]79 +1.0*reaction\[LetterSpace]80 -1.0*reaction\[LetterSpace]81, pRSK\[LetterSpace]c'[t] == 1.0*reaction\[LetterSpace]41 -1.0*re113 -1.0*reaction\[LetterSpace]42, pRSK\[LetterSpace]n'[t] == 1.0*re113 , pc\[LetterSpace]FOS\[LetterSpace]c'[t] == 1.0*reaction\[LetterSpace]59 +1.0*reaction\[LetterSpace]67 -1.0*re117 -1.0*reaction\[LetterSpace]60 -1.0*reaction\[LetterSpace]62, ppERK\[LetterSpace]c'[t] == 1.0*reaction\[LetterSpace]6 -1.0*re108 -1.0*re86, ppERK\[LetterSpace]n'[t] == 1.0*re108 -1.0*re88 -1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]25 }; 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]}]