(* 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 = {
Akt[t], BRaf[t], BRaf\[LetterSpace]inhibitor[t], EGFR[t], EGFR\[LetterSpace]Sos[t], EGFR\[LetterSpace]inhibitor[t], EGF\[LetterSpace]EGFR[t], ERKcyt[t], ERKnuc[t], MEK\[LetterSpace]ERK\[LetterSpace]cyt[t], MEK\[LetterSpace]inhibitor[t], MEKcyt[t], PI3K[t], PI3K\[LetterSpace]inhibitor[t], RSK[t], RasEffector[t], RasGAP[t], RasGEF[t], Ras\[LetterSpace]GDP[t], Ras\[LetterSpace]GTP[t], Ras\[LetterSpace]Raf[t], Ras\[LetterSpace]Shoc2[t], Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf[t], RassEffector[t], Rheb\[LetterSpace]GDP[t], Rheb\[LetterSpace]GTP[t], S6K[t], Shoc2[t], Sos1\[LetterSpace]2[t], TSC[t], active\[LetterSpace]PI3K[t], cRaf[t], degra\[LetterSpace]EGF\[LetterSpace]pEGFR[t], inactive\[LetterSpace]cRaf[t], mutated\[LetterSpace]PI3K[t], mutated\[LetterSpace]Ras[t], pAkt[t], pMEK\[LetterSpace]ERK\[LetterSpace]cyt[t], pMEK\[LetterSpace]pYERK[t], pMEKcyt[t], pRSK[t], pS6K[t], pSos[t], pTERK\[LetterSpace]cyt[t], pTERK\[LetterSpace]nuc[t], pTSC[t], pTpYERK\[LetterSpace]cyt[t], pTpYERK\[LetterSpace]nuc[t], pYERK\[LetterSpace]cyt[t], pYERK\[LetterSpace]nuc[t], ppSos[t], pppSos[t], ppppSos[t]
};

initialValues = {
Akt[0] == 1.0, BRaf[0] == 0.0, BRaf\[LetterSpace]inhibitor[0] == 1.0, EGFR[0] == 0.42, EGFR\[LetterSpace]Sos[0] == 0.0, EGFR\[LetterSpace]inhibitor[0] == 1.0, EGF\[LetterSpace]EGFR[0] == 0.0, ERKcyt[0] == 0.96, ERKnuc[0] == 0.0, MEK\[LetterSpace]ERK\[LetterSpace]cyt[0] == 0.0, MEK\[LetterSpace]inhibitor[0] == 1.0, MEKcyt[0] == 1.4, PI3K[0] == 1.0, PI3K\[LetterSpace]inhibitor[0] == 1.0, RSK[0] == 1.0, RasEffector[0] == 2.0, RasGAP[0] == 1.0, RasGEF[0] == 1.0, Ras\[LetterSpace]GDP[0] == 0.43, Ras\[LetterSpace]GTP[0] == 0.0, Ras\[LetterSpace]Raf[0] == 0.0, Ras\[LetterSpace]Shoc2[0] == 0.0, Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf[0] == 0.0, RassEffector[0] == 0.0, Rheb\[LetterSpace]GDP[0] == 1.0, Rheb\[LetterSpace]GTP[0] == 0.0, S6K[0] == 1.0, Shoc2[0] == 0.67, Sos1\[LetterSpace]2[0] == 0.12, TSC[0] == 1.0, active\[LetterSpace]PI3K[0] == 0.0, cRaf[0] == 0.013, degra\[LetterSpace]EGF\[LetterSpace]pEGFR[0] == 0.0, inactive\[LetterSpace]cRaf[0] == 0.0, mutated\[LetterSpace]PI3K[0] == 0.0, mutated\[LetterSpace]Ras[0] == 0.0, pAkt[0] == 0.0, pMEK\[LetterSpace]ERK\[LetterSpace]cyt[0] == 0.0, pMEK\[LetterSpace]pYERK[0] == 0.0, pMEKcyt[0] == 0.0, pRSK[0] == 0.0, pS6K[0] == 0.0, pSos[0] == 0.0, pTERK\[LetterSpace]cyt[0] == 0.0, pTERK\[LetterSpace]nuc[0] == 0.0, pTSC[0] == 0.0, pTpYERK\[LetterSpace]cyt[0] == 0.0, pTpYERK\[LetterSpace]nuc[0] == 0.0, pYERK\[LetterSpace]cyt[0] == 0.0, pYERK\[LetterSpace]nuc[0] == 0.0, ppSos[0] == 0.0, pppSos[0] == 0.0, ppppSos[0] == 0.0
};

rates = {
re1, re12, re15, re16, re2, re22, re27, re28, re3, re30, re31, re32, re4, re42, re43, re44, re45, re46, re47, re48, re49, re5, re50, re51, re52, re53, re54, re55, re56, re57, re58, re6, re62, re63, re64, re65, re66, re67, re68, re69, re70, re9
};

rateEquations = {
re1 -> function\[LetterSpace]1[EGF, EGFR[t], kf\[LetterSpace]EGFEGFR, kphos\[LetterSpace]EGFR\[LetterSpace]mutant], re12 -> c1*function\[LetterSpace]8[ERKcyt[t], c1, kb\[LetterSpace]MEKERK, kf\[LetterSpace]MEKERK, pMEK\[LetterSpace]ERK\[LetterSpace]cyt[t], pMEKcyt[t]], re15 -> function\[LetterSpace]9[c2, kexport\[LetterSpace]pTpYERK, kimport\[LetterSpace]pTpYERK, pTpYERK\[LetterSpace]cyt[t], pTpYERK\[LetterSpace]nuc[t]], re16 -> function\[LetterSpace]10[ERKcyt[t], ERKnuc[t], c2, kexport\[LetterSpace]ERK, kimport\[LetterSpace]ERK], re2 -> c1*function\[LetterSpace]2[EGF\[LetterSpace]EGFR[t], c1, degra\[LetterSpace]EGF\[LetterSpace]pEGFR[t], kb\[LetterSpace]degra, kf\[LetterSpace]degra], re22 -> c1*function\[LetterSpace]11[c1, kf\[LetterSpace]dpMEK, pMEKcyt[t]], re27 -> c1*function\[LetterSpace]12[RasEffector[t], RassEffector[t], Ras\[LetterSpace]GTP[t], c1, kb\[LetterSpace]RasEffector, kf\[LetterSpace]RasEffector], re28 -> c1*function\[LetterSpace]13[Km\[LetterSpace]pERKcyt, Sos1\[LetterSpace]2[t], c1, kb\[LetterSpace]pSos, kcat\[LetterSpace]pERKcyt, kf\[LetterSpace]pSos, pSos[t], pTpYERK\[LetterSpace]cyt[t]], re3 -> c1*function\[LetterSpace]3[EGFR\[LetterSpace]Sos[t], EGFR\[LetterSpace]inhibitor[t], EGF\[LetterSpace]EGFR[t], Sos1\[LetterSpace]2[t], c1, kb\[LetterSpace]Sos, kf\[LetterSpace]Sos], re30 -> c1*function\[LetterSpace]14[Km\[LetterSpace]pERKcyt, c1, kb\[LetterSpace]pSos, kcat\[LetterSpace]pERKcyt, kfb\[LetterSpace]ppSos, pSos[t], pTpYERK\[LetterSpace]cyt[t], ppSos[t]], re31 -> c1*function\[LetterSpace]15[Km\[LetterSpace]pERKcyt, c1, kb\[LetterSpace]pSos, kcat\[LetterSpace]pERKcyt, kfb\[LetterSpace]pppSos, pTpYERK\[LetterSpace]cyt[t], ppSos[t], pppSos[t]], re32 -> c1*function\[LetterSpace]16[Km\[LetterSpace]pERKcyt, c1, kb\[LetterSpace]pSos, kcat\[LetterSpace]pERKcyt, kfb\[LetterSpace]ppppSos, pTpYERK\[LetterSpace]cyt[t], pppSos[t], ppppSos[t]], re4 -> c1*function\[LetterSpace]4[EGFR\[LetterSpace]Sos[t], RasGEF[t], Ras\[LetterSpace]GDP[t], c1, kf\[LetterSpace]RasGEF, kf\[LetterSpace]Sos\[LetterSpace]Ras], re42 -> c1*function\[LetterSpace]17[MEK\[LetterSpace]inhibitor[t], c1, kphos\[LetterSpace]np\[LetterSpace]pY, pMEK\[LetterSpace]ERK\[LetterSpace]cyt[t]], re43 -> c1*function\[LetterSpace]18[MEK\[LetterSpace]inhibitor[t], c1, kphos\[LetterSpace]pY\[LetterSpace]pTpY\[LetterSpace]MEKERK, pMEK\[LetterSpace]pYERK[t]], re44 -> c1*function\[LetterSpace]19[c1, kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]cyt, pTpYERK\[LetterSpace]cyt[t]], re45 -> c1*function\[LetterSpace]20[c1, kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]cyt, pTpYERK\[LetterSpace]cyt[t]], re46 -> c1*function\[LetterSpace]21[c1, kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]cyt, pYERK\[LetterSpace]cyt[t]], re47 -> c1*function\[LetterSpace]22[c1, kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]cyt, pTERK\[LetterSpace]cyt[t]], re48 -> c1*function\[LetterSpace]23[MEK\[LetterSpace]inhibitor[t], c1, kphos\[LetterSpace]pY\[LetterSpace]pTpY, pMEKcyt[t], pYERK\[LetterSpace]cyt[t]], re49 -> c1*function\[LetterSpace]24[MEK\[LetterSpace]inhibitor[t], c1, kphos\[LetterSpace]pT\[LetterSpace]pTpY, pMEKcyt[t], pTERK\[LetterSpace]cyt[t]], re5 -> c1*function\[LetterSpace]5[RasGAP[t], Ras\[LetterSpace]GTP[t], c1, kf\[LetterSpace]RasGAP, mutated\[LetterSpace]Ras[t]], re50 -> function\[LetterSpace]25[c2, kexport\[LetterSpace]pYERK, kimport\[LetterSpace]pYERK, pYERK\[LetterSpace]cyt[t], pYERK\[LetterSpace]nuc[t]], re51 -> function\[LetterSpace]26[c2, kexport\[LetterSpace]pTERK, kimport\[LetterSpace]pTERK, pTERK\[LetterSpace]cyt[t], pTERK\[LetterSpace]nuc[t]], re52 -> c2*kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]nuc*pTpYERK\[LetterSpace]nuc[t], re53 -> c2*kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]nuc*pTpYERK\[LetterSpace]nuc[t], re54 -> c2*kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]nuc*pYERK\[LetterSpace]nuc[t], re55 -> c2*kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]nuc*pTERK\[LetterSpace]nuc[t], re56 -> c1*function\[LetterSpace]27[BRaf[t], BRaf\[LetterSpace]inhibitor[t], MEKcyt[t], Ras\[LetterSpace]Raf[t], Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf[t], c1, kf\[LetterSpace]pMEK, kphos\[LetterSpace]BRafMEK], re57 -> c1*function\[LetterSpace]28[Ras\[LetterSpace]GTP[t], Ras\[LetterSpace]Shoc2[t], Shoc2[t], c1, kb\[LetterSpace]Ras\[LetterSpace]Shoc2, kf\[LetterSpace]Ras\[LetterSpace]Shoc2], re58 -> c1*function\[LetterSpace]29[Ras\[LetterSpace]Shoc2[t], Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf[t], c1, cRaf[t], kb\[LetterSpace]RasRaf\[LetterSpace]Shoc2, kf\[LetterSpace]RasRaf\[LetterSpace]Shoc2], re6 -> c1*function\[LetterSpace]6[Ras\[LetterSpace]GTP[t], Ras\[LetterSpace]Raf[t], c1, cRaf[t], kb\[LetterSpace]RasRaf, kf\[LetterSpace]RasRaf], re62 -> c1*function\[LetterSpace]30[EGFR\[LetterSpace]inhibitor[t], EGF\[LetterSpace]EGFR[t], PI3K[t], Ras\[LetterSpace]GTP[t], active\[LetterSpace]PI3K[t], c1, kact\[LetterSpace]EGFRRasPI3K, kdact\[LetterSpace]PI3K, mutated\[LetterSpace]PI3K[t]], re63 -> c1*function\[LetterSpace]31[Akt[t], PI3K\[LetterSpace]inhibitor[t], active\[LetterSpace]PI3K[t], c1, kdphos\[LetterSpace]pAkt, kphos\[LetterSpace]PI3KAkt, pAkt[t]], re64 -> c1*function\[LetterSpace]32[TSC[t], c1, kdphos\[LetterSpace]pTSC, kphos\[LetterSpace]AktTSC, kphos\[LetterSpace]RSKTSC, pAkt[t], pRSK[t], pTSC[t]], re65 -> c1*function\[LetterSpace]33[Rheb\[LetterSpace]GDP[t], Rheb\[LetterSpace]GTP[t], TSC[t], c1, kact\[LetterSpace]Rheb, kdact\[LetterSpace]TSCRheb], re66 -> c1*function\[LetterSpace]34[Rheb\[LetterSpace]GTP[t], S6K[t], c1, kdphos\[LetterSpace]pS6K, kphos\[LetterSpace]RhebS6K, pS6K[t]], re67 -> c1*function\[LetterSpace]35[RSK[t], c1, kdphos\[LetterSpace]pRSK, kphos\[LetterSpace]ERKRSK, pRSK[t], pTpYERK\[LetterSpace]cyt[t]], re68 -> c1*function\[LetterSpace]36[Rheb\[LetterSpace]GTP[t], c1, cRaf[t], inactive\[LetterSpace]cRaf[t], kact\[LetterSpace]RhebRaf, kdact\[LetterSpace]RhebRaf], re69 -> c1*function\[LetterSpace]37[EGF\[LetterSpace]EGFR[t], c1, kf\[LetterSpace]degra\[LetterSpace]ERKEGFR, pTpYERK\[LetterSpace]cyt[t]], re70 -> c1*function\[LetterSpace]38[PI3K\[LetterSpace]inhibitor[t], active\[LetterSpace]PI3K[t], c1, inactive\[LetterSpace]cRaf[t], kact\[LetterSpace]PI3KRaf], re9 -> c1*function\[LetterSpace]7[ERKcyt[t], MEK\[LetterSpace]ERK\[LetterSpace]cyt[t], MEKcyt[t], c1, kb\[LetterSpace]MEKERK, kf\[LetterSpace]MEKERK]
};

parameters = {
Km\[LetterSpace]pERKcyt -> 3.0, kact\[LetterSpace]EGFRRasPI3K -> 3.0, kact\[LetterSpace]PI3KRaf -> 1.8, kact\[LetterSpace]Rheb -> 0.0021, kact\[LetterSpace]RhebRaf -> 0.04, kb\[LetterSpace]MEKERK -> 0.27, kb\[LetterSpace]RafMEK -> 0.065, kb\[LetterSpace]RasEffector -> 0.1, kb\[LetterSpace]RasRaf -> 0.0049, kb\[LetterSpace]RasRaf\[LetterSpace]Shoc2 -> 0.049, kb\[LetterSpace]Ras\[LetterSpace]Shoc2 -> 0.1, kb\[LetterSpace]Sos -> 0.1, kb\[LetterSpace]degra -> 0.001, kb\[LetterSpace]pSos -> 0.0025, kcat\[LetterSpace]pERKcyt -> 0.075, kdact\[LetterSpace]PI3K -> 0.028, kdact\[LetterSpace]RhebRaf -> 4.0, kdact\[LetterSpace]TSCRheb -> 0.097, kdphos\[LetterSpace]pAkt -> 0.048, kdphos\[LetterSpace]pRSK -> 0.3, kdphos\[LetterSpace]pS6K -> 0.001, kdphos\[LetterSpace]pTSC -> 0.0088, kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]cyt -> 0.0068, kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]nuc -> 0.0058, kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]cyt -> 0.0055, kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]nuc -> 0.0038, kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]cyt -> 0.004, kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]nuc -> 0.0032, kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]cyt -> 0.0067, kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]nuc -> 0.0077, kexport\[LetterSpace]ERK -> 0.013, kexport\[LetterSpace]pTERK -> 0.049, kexport\[LetterSpace]pTpYERK -> 0.0076, kexport\[LetterSpace]pYERK -> 0.017, kf\[LetterSpace]EGFEGFR -> 10.0, kf\[LetterSpace]MEKERK -> 0.18, kf\[LetterSpace]RafMEK -> 0.65, kf\[LetterSpace]RasEffector -> 1.0, kf\[LetterSpace]RasGAP -> 0.05, kf\[LetterSpace]RasGEF -> 9*^-05, kf\[LetterSpace]RasRaf -> 0.049, kf\[LetterSpace]RasRaf\[LetterSpace]Shoc2 -> 0.49, kf\[LetterSpace]Ras\[LetterSpace]Shoc2 -> 1.0, kf\[LetterSpace]Sos -> 0.1, kf\[LetterSpace]Sos\[LetterSpace]Ras -> 1.0, kf\[LetterSpace]degra -> 0.0014, kf\[LetterSpace]degra\[LetterSpace]ERKEGFR -> 0.02, kf\[LetterSpace]dpMEK -> 0.0096, kf\[LetterSpace]pMEK -> 2.0, kf\[LetterSpace]pSos -> 1.0, kfb\[LetterSpace]ppSos -> 1.0, kfb\[LetterSpace]pppSos -> 1.0, kfb\[LetterSpace]ppppSos -> 1.0, kimport\[LetterSpace]ERK -> 0.0017, kimport\[LetterSpace]pTERK -> 0.0022, kimport\[LetterSpace]pTpYERK -> 0.0082, kimport\[LetterSpace]pYERK -> 0.0025, kphos\[LetterSpace]AktTSC -> 1.2, kphos\[LetterSpace]BRafMEK -> 0.01, kphos\[LetterSpace]EGFR\[LetterSpace]mutant -> 0.0, kphos\[LetterSpace]ERKRSK -> 0.7, kphos\[LetterSpace]PI3KAkt -> 0.1, kphos\[LetterSpace]RSKTSC -> 0.6, kphos\[LetterSpace]RhebS6K -> 0.002, kphos\[LetterSpace]np\[LetterSpace]pY -> 0.073, kphos\[LetterSpace]pT\[LetterSpace]pTpY -> 0.02, kphos\[LetterSpace]pY\[LetterSpace]pTpY -> 0.02, kphos\[LetterSpace]pY\[LetterSpace]pTpY\[LetterSpace]MEKERK -> 0.05, EGF -> 0.0, c1 -> 1.0, c2 -> 0.2, default -> 1.0
};

assignments = {
function\[LetterSpace]31[Akt_,PI3K\[LetterSpace]inhibitor_,active\[LetterSpace]PI3K_,c1_,kdphos\[LetterSpace]pAkt_,kphos\[LetterSpace]PI3KAkt_,pAkt_] -> (-(kdphos\[LetterSpace]pAkt*pAkt) + active\[LetterSpace]PI3K*Akt*c1*kphos\[LetterSpace]PI3KAkt*PI3K\[LetterSpace]inhibitor)/c1, function\[LetterSpace]11[c1_,kf\[LetterSpace]dpMEK_,pMEKcyt_] -> (kf\[LetterSpace]dpMEK*pMEKcyt)/c1, function\[LetterSpace]32[TSC_,c1_,kdphos\[LetterSpace]pTSC_,kphos\[LetterSpace]AktTSC_,kphos\[LetterSpace]RSKTSC_,pAkt_,pRSK_,pTSC_] -> (-(kdphos\[LetterSpace]pTSC*pTSC) + kphos\[LetterSpace]AktTSC*pAkt*TSC + kphos\[LetterSpace]RSKTSC*pRSK*TSC)/c1, function\[LetterSpace]28[Ras\[LetterSpace]GTP_,Ras\[LetterSpace]Shoc2_,Shoc2_,c1_,kb\[LetterSpace]Ras\[LetterSpace]Shoc2_,kf\[LetterSpace]Ras\[LetterSpace]Shoc2_] -> (-(kb\[LetterSpace]Ras\[LetterSpace]Shoc2*Ras\[LetterSpace]Shoc2) + kf\[LetterSpace]Ras\[LetterSpace]Shoc2*Ras\[LetterSpace]GTP*Shoc2)/c1, function\[LetterSpace]25[c2_,kexport\[LetterSpace]pYERK_,kimport\[LetterSpace]pYERK_,pYERK\[LetterSpace]cyt_,pYERK\[LetterSpace]nuc_] -> kimport\[LetterSpace]pYERK*pYERK\[LetterSpace]cyt - c2*kexport\[LetterSpace]pYERK*pYERK\[LetterSpace]nuc, function\[LetterSpace]21[c1_,kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]cyt_,pYERK\[LetterSpace]cyt_] -> (kdphos\[LetterSpace]pY\[LetterSpace]np\[LetterSpace]cyt*pYERK\[LetterSpace]cyt)/c1, function\[LetterSpace]19[c1_,kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]cyt_,pTpYERK\[LetterSpace]cyt_] -> (kdphos\[LetterSpace]pTpY\[LetterSpace]pT\[LetterSpace]cyt*pTpYERK\[LetterSpace]cyt)/c1, function\[LetterSpace]17[MEK\[LetterSpace]inhibitor_,c1_,kphos\[LetterSpace]np\[LetterSpace]pY_,pMEK\[LetterSpace]ERK\[LetterSpace]cyt_] -> kphos\[LetterSpace]np\[LetterSpace]pY*MEK\[LetterSpace]inhibitor*pMEK\[LetterSpace]ERK\[LetterSpace]cyt, function\[LetterSpace]16[Km\[LetterSpace]pERKcyt_,c1_,kb\[LetterSpace]pSos_,kcat\[LetterSpace]pERKcyt_,kfb\[LetterSpace]ppppSos_,pTpYERK\[LetterSpace]cyt_,pppSos_,ppppSos_] -> (kfb\[LetterSpace]ppppSos*(-(kb\[LetterSpace]pSos*ppppSos) + (kcat\[LetterSpace]pERKcyt*pppSos*pTpYERK\[LetterSpace]cyt)/(Km\[LetterSpace]pERKcyt + pppSos)))/c1, function\[LetterSpace]13[Km\[LetterSpace]pERKcyt_,Sos1\[LetterSpace]2_,c1_,kb\[LetterSpace]pSos_,kcat\[LetterSpace]pERKcyt_,kf\[LetterSpace]pSos_,pSos_,pTpYERK\[LetterSpace]cyt_] -> (-(kb\[LetterSpace]pSos*pSos) + (kcat\[LetterSpace]pERKcyt*kf\[LetterSpace]pSos*pTpYERK\[LetterSpace]cyt*Sos1\[LetterSpace]2)/(Km\[LetterSpace]pERKcyt + Sos1\[LetterSpace]2))/c1, function\[LetterSpace]9[c2_,kexport\[LetterSpace]pTpYERK_,kimport\[LetterSpace]pTpYERK_,pTpYERK\[LetterSpace]cyt_,pTpYERK\[LetterSpace]nuc_] -> kimport\[LetterSpace]pTpYERK*pTpYERK\[LetterSpace]cyt - c2*kexport\[LetterSpace]pTpYERK*pTpYERK\[LetterSpace]nuc, function\[LetterSpace]33[Rheb\[LetterSpace]GDP_,Rheb\[LetterSpace]GTP_,TSC_,c1_,kact\[LetterSpace]Rheb_,kdact\[LetterSpace]TSCRheb_] -> (-(kact\[LetterSpace]Rheb*Rheb\[LetterSpace]GDP) + kdact\[LetterSpace]TSCRheb*Rheb\[LetterSpace]GTP*TSC)/c1, function\[LetterSpace]29[Ras\[LetterSpace]Shoc2_,Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf_,c1_,cRaf_,kb\[LetterSpace]RasRaf\[LetterSpace]Shoc2_,kf\[LetterSpace]RasRaf\[LetterSpace]Shoc2_] -> (cRaf*kf\[LetterSpace]RasRaf\[LetterSpace]Shoc2*Ras\[LetterSpace]Shoc2 - kb\[LetterSpace]RasRaf\[LetterSpace]Shoc2*Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf)/c1, function\[LetterSpace]26[c2_,kexport\[LetterSpace]pTERK_,kimport\[LetterSpace]pTERK_,pTERK\[LetterSpace]cyt_,pTERK\[LetterSpace]nuc_] -> kimport\[LetterSpace]pTERK*pTERK\[LetterSpace]cyt - c2*kexport\[LetterSpace]pTERK*pTERK\[LetterSpace]nuc, function\[LetterSpace]24[MEK\[LetterSpace]inhibitor_,c1_,kphos\[LetterSpace]pT\[LetterSpace]pTpY_,pMEKcyt_,pTERK\[LetterSpace]cyt_] -> kphos\[LetterSpace]pT\[LetterSpace]pTpY*MEK\[LetterSpace]inhibitor*pMEKcyt*pTERK\[LetterSpace]cyt, function\[LetterSpace]23[MEK\[LetterSpace]inhibitor_,c1_,kphos\[LetterSpace]pY\[LetterSpace]pTpY_,pMEKcyt_,pYERK\[LetterSpace]cyt_] -> kphos\[LetterSpace]pY\[LetterSpace]pTpY*MEK\[LetterSpace]inhibitor*pMEKcyt*pYERK\[LetterSpace]cyt, function\[LetterSpace]22[c1_,kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]cyt_,pTERK\[LetterSpace]cyt_] -> (kdphos\[LetterSpace]pT\[LetterSpace]np\[LetterSpace]cyt*pTERK\[LetterSpace]cyt)/c1, function\[LetterSpace]20[c1_,kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]cyt_,pTpYERK\[LetterSpace]cyt_] -> (kdphos\[LetterSpace]pTpY\[LetterSpace]pY\[LetterSpace]cyt*pTpYERK\[LetterSpace]cyt)/c1, function\[LetterSpace]18[MEK\[LetterSpace]inhibitor_,c1_,kphos\[LetterSpace]pY\[LetterSpace]pTpY\[LetterSpace]MEKERK_,pMEK\[LetterSpace]pYERK_] -> kphos\[LetterSpace]pY\[LetterSpace]pTpY\[LetterSpace]MEKERK*MEK\[LetterSpace]inhibitor*pMEK\[LetterSpace]pYERK, function\[LetterSpace]15[Km\[LetterSpace]pERKcyt_,c1_,kb\[LetterSpace]pSos_,kcat\[LetterSpace]pERKcyt_,kfb\[LetterSpace]pppSos_,pTpYERK\[LetterSpace]cyt_,ppSos_,pppSos_] -> (kfb\[LetterSpace]pppSos*(-(kb\[LetterSpace]pSos*pppSos) + (kcat\[LetterSpace]pERKcyt*ppSos*pTpYERK\[LetterSpace]cyt)/(Km\[LetterSpace]pERKcyt + ppSos)))/c1, function\[LetterSpace]14[Km\[LetterSpace]pERKcyt_,c1_,kb\[LetterSpace]pSos_,kcat\[LetterSpace]pERKcyt_,kfb\[LetterSpace]ppSos_,pSos_,pTpYERK\[LetterSpace]cyt_,ppSos_] -> (kfb\[LetterSpace]ppSos*(-(kb\[LetterSpace]pSos*ppSos) + (kcat\[LetterSpace]pERKcyt*pSos*pTpYERK\[LetterSpace]cyt)/(Km\[LetterSpace]pERKcyt + pSos)))/c1, function\[LetterSpace]7[ERKcyt_,MEK\[LetterSpace]ERK\[LetterSpace]cyt_,MEKcyt_,c1_,kb\[LetterSpace]MEKERK_,kf\[LetterSpace]MEKERK_] -> (ERKcyt*kf\[LetterSpace]MEKERK*MEKcyt - kb\[LetterSpace]MEKERK*MEK\[LetterSpace]ERK\[LetterSpace]cyt)/c1, function\[LetterSpace]36[Rheb\[LetterSpace]GTP_,c1_,cRaf_,inactive\[LetterSpace]cRaf_,kact\[LetterSpace]RhebRaf_,kdact\[LetterSpace]RhebRaf_] -> (-(inactive\[LetterSpace]cRaf*kact\[LetterSpace]RhebRaf) + cRaf*kdact\[LetterSpace]RhebRaf*Rheb\[LetterSpace]GTP)/c1, function\[LetterSpace]1[EGF_,EGFR_,kf\[LetterSpace]EGFEGFR_,kphos\[LetterSpace]EGFR\[LetterSpace]mutant_] -> EGF*EGFR*kf\[LetterSpace]EGFEGFR + EGFR*kphos\[LetterSpace]EGFR\[LetterSpace]mutant, function\[LetterSpace]5[RasGAP_,Ras\[LetterSpace]GTP_,c1_,kf\[LetterSpace]RasGAP_,mutated\[LetterSpace]Ras_] -> (kf\[LetterSpace]RasGAP*RasGAP*(-mutated\[LetterSpace]Ras + Ras\[LetterSpace]GTP))/c1, function\[LetterSpace]10[ERKcyt_,ERKnuc_,c2_,kexport\[LetterSpace]ERK_,kimport\[LetterSpace]ERK_] -> -(c2*ERKnuc*kexport\[LetterSpace]ERK) + ERKcyt*kimport\[LetterSpace]ERK, function\[LetterSpace]8[ERKcyt_,c1_,kb\[LetterSpace]MEKERK_,kf\[LetterSpace]MEKERK_,pMEK\[LetterSpace]ERK\[LetterSpace]cyt_,pMEKcyt_] -> (ERKcyt*kf\[LetterSpace]MEKERK*pMEKcyt - kb\[LetterSpace]MEKERK*pMEK\[LetterSpace]ERK\[LetterSpace]cyt)/c1, function\[LetterSpace]6[Ras\[LetterSpace]GTP_,Ras\[LetterSpace]Raf_,c1_,cRaf_,kb\[LetterSpace]RasRaf_,kf\[LetterSpace]RasRaf_] -> (cRaf*kf\[LetterSpace]RasRaf*Ras\[LetterSpace]GTP - kb\[LetterSpace]RasRaf*Ras\[LetterSpace]Raf)/c1, function\[LetterSpace]38[PI3K\[LetterSpace]inhibitor_,active\[LetterSpace]PI3K_,c1_,inactive\[LetterSpace]cRaf_,kact\[LetterSpace]PI3KRaf_] -> active\[LetterSpace]PI3K*inactive\[LetterSpace]cRaf*kact\[LetterSpace]PI3KRaf*PI3K\[LetterSpace]inhibitor, function\[LetterSpace]35[RSK_,c1_,kdphos\[LetterSpace]pRSK_,kphos\[LetterSpace]ERKRSK_,pRSK_,pTpYERK\[LetterSpace]cyt_] -> (-(kdphos\[LetterSpace]pRSK*pRSK) + kphos\[LetterSpace]ERKRSK*pTpYERK\[LetterSpace]cyt*RSK)/c1, function\[LetterSpace]37[EGF\[LetterSpace]EGFR_,c1_,kf\[LetterSpace]degra\[LetterSpace]ERKEGFR_,pTpYERK\[LetterSpace]cyt_] -> (EGF\[LetterSpace]EGFR*kf\[LetterSpace]degra\[LetterSpace]ERKEGFR*pTpYERK\[LetterSpace]cyt)/c1, function\[LetterSpace]27[BRaf_,BRaf\[LetterSpace]inhibitor_,MEKcyt_,Ras\[LetterSpace]Raf_,Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf_,c1_,kf\[LetterSpace]pMEK_,kphos\[LetterSpace]BRafMEK_] -> (BRaf*BRaf\[LetterSpace]inhibitor*c1*kphos\[LetterSpace]BRafMEK + kf\[LetterSpace]pMEK*MEKcyt*(Ras\[LetterSpace]Raf + Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf))/c1, function\[LetterSpace]4[EGFR\[LetterSpace]Sos_,RasGEF_,Ras\[LetterSpace]GDP_,c1_,kf\[LetterSpace]RasGEF_,kf\[LetterSpace]Sos\[LetterSpace]Ras_] -> (EGFR\[LetterSpace]Sos*kf\[LetterSpace]Sos\[LetterSpace]Ras*Ras\[LetterSpace]GDP + kf\[LetterSpace]RasGEF*RasGEF*Ras\[LetterSpace]GDP)/c1, function\[LetterSpace]3[EGFR\[LetterSpace]Sos_,EGFR\[LetterSpace]inhibitor_,EGF\[LetterSpace]EGFR_,Sos1\[LetterSpace]2_,c1_,kb\[LetterSpace]Sos_,kf\[LetterSpace]Sos_] -> (-(EGFR\[LetterSpace]Sos*kb\[LetterSpace]Sos) + c1*EGFR\[LetterSpace]inhibitor*EGF\[LetterSpace]EGFR*kf\[LetterSpace]Sos*Sos1\[LetterSpace]2)/c1, function\[LetterSpace]2[EGF\[LetterSpace]EGFR_,c1_,degra\[LetterSpace]EGF\[LetterSpace]pEGFR_,kb\[LetterSpace]degra_,kf\[LetterSpace]degra_] -> (-(degra\[LetterSpace]EGF\[LetterSpace]pEGFR*kb\[LetterSpace]degra) + EGF\[LetterSpace]EGFR*kf\[LetterSpace]degra)/c1, function\[LetterSpace]12[RasEffector_,RassEffector_,Ras\[LetterSpace]GTP_,c1_,kb\[LetterSpace]RasEffector_,kf\[LetterSpace]RasEffector_] -> (-(kb\[LetterSpace]RasEffector*RassEffector) + kf\[LetterSpace]RasEffector*RasEffector*Ras\[LetterSpace]GTP)/c1, function\[LetterSpace]30[EGFR\[LetterSpace]inhibitor_,EGF\[LetterSpace]EGFR_,PI3K_,Ras\[LetterSpace]GTP_,active\[LetterSpace]PI3K_,c1_,kact\[LetterSpace]EGFRRasPI3K_,kdact\[LetterSpace]PI3K_,mutated\[LetterSpace]PI3K_] -> (-(kdact\[LetterSpace]PI3K*(active\[LetterSpace]PI3K - mutated\[LetterSpace]PI3K)) + c1*EGFR\[LetterSpace]inhibitor*EGF\[LetterSpace]EGFR*kact\[LetterSpace]EGFRRasPI3K*PI3K*Ras\[LetterSpace]GTP)/c1, function\[LetterSpace]34[Rheb\[LetterSpace]GTP_,S6K_,c1_,kdphos\[LetterSpace]pS6K_,kphos\[LetterSpace]RhebS6K_,pS6K_] -> (-(kdphos\[LetterSpace]pS6K*pS6K) + kphos\[LetterSpace]RhebS6K*Rheb\[LetterSpace]GTP*S6K)/c1
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Akt'[t] ==  -1.0*re63, BRaf'[t] == 0.0 , BRaf\[LetterSpace]inhibitor'[t] == 0.0 , EGFR'[t] ==  -1.0*re1, EGFR\[LetterSpace]Sos'[t] == 1.0*re3 , EGFR\[LetterSpace]inhibitor'[t] == 0.0 , EGF\[LetterSpace]EGFR'[t] == 1.0*re1 -1.0*re2 -1.0*re3 -1.0*re69, ERKcyt'[t] == 1.0*re46 +1.0*re47 -1.0*re12 -1.0*re16 -1.0*re9, ERKnuc'[t] == 1.0*re16 +1.0*re54 +1.0*re55 , MEK\[LetterSpace]ERK\[LetterSpace]cyt'[t] == 1.0*re9 , MEK\[LetterSpace]inhibitor'[t] == 0.0 , MEKcyt'[t] == 1.0*re22 -1.0*re56 -1.0*re9, PI3K'[t] ==  -1.0*re62, PI3K\[LetterSpace]inhibitor'[t] == 0.0 , RSK'[t] ==  -1.0*re67, RasEffector'[t] ==  -1.0*re27, RasGAP'[t] == 0.0 , RasGEF'[t] == 0.0 , Ras\[LetterSpace]GDP'[t] == 1.0*re5 -1.0*re4, Ras\[LetterSpace]GTP'[t] == 1.0*re4 -1.0*re27 -1.0*re5 -1.0*re57 -1.0*re6, Ras\[LetterSpace]Raf'[t] == 1.0*re6 , Ras\[LetterSpace]Shoc2'[t] == 1.0*re57 -1.0*re58, Ras\[LetterSpace]Shoc2\[LetterSpace]cRaf'[t] == 1.0*re58 , RassEffector'[t] == 1.0*re27 , Rheb\[LetterSpace]GDP'[t] == 1.0*re65 , Rheb\[LetterSpace]GTP'[t] ==  -1.0*re65, S6K'[t] ==  -1.0*re66, Shoc2'[t] ==  -1.0*re57, Sos1\[LetterSpace]2'[t] ==  -1.0*re28 -1.0*re3, TSC'[t] ==  -1.0*re64, active\[LetterSpace]PI3K'[t] == 1.0*re62 , cRaf'[t] == 1.0*re70 -1.0*re58 -1.0*re6 -1.0*re68, degra\[LetterSpace]EGF\[LetterSpace]pEGFR'[t] == 1.0*re2 +1.0*re69 , inactive\[LetterSpace]cRaf'[t] == 1.0*re68 -1.0*re70, mutated\[LetterSpace]PI3K'[t] == 0.0 , mutated\[LetterSpace]Ras'[t] == 0.0 , pAkt'[t] == 1.0*re63 , pMEK\[LetterSpace]ERK\[LetterSpace]cyt'[t] == 1.0*re12 -1.0*re42, pMEK\[LetterSpace]pYERK'[t] == 1.0*re42 -1.0*re43, pMEKcyt'[t] == 1.0*re43 +1.0*re56 -1.0*re12 -1.0*re22, pRSK'[t] == 1.0*re67 , pS6K'[t] == 1.0*re66 , pSos'[t] == 1.0*re28 -1.0*re30, pTERK\[LetterSpace]cyt'[t] == 1.0*re44 -1.0*re47 -1.0*re49 -1.0*re51, pTERK\[LetterSpace]nuc'[t] == 1.0*re51 +1.0*re53 -1.0*re55, pTSC'[t] == 1.0*re64 , pTpYERK\[LetterSpace]cyt'[t] == 1.0*re43 +1.0*re48 +1.0*re49 -1.0*re15 -1.0*re44 -1.0*re45, pTpYERK\[LetterSpace]nuc'[t] == 1.0*re15 -1.0*re52 -1.0*re53, pYERK\[LetterSpace]cyt'[t] == 1.0*re45 -1.0*re46 -1.0*re48 -1.0*re50, pYERK\[LetterSpace]nuc'[t] == 1.0*re50 +1.0*re52 -1.0*re54, ppSos'[t] == 1.0*re30 -1.0*re31, pppSos'[t] == 1.0*re31 -1.0*re32, ppppSos'[t] == 1.0*re32 
};
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]}]