(* 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], EGF[t], Erk[t], GAB[t], GABp[t], GABp\[LetterSpace]GS[t], GABp\[LetterSpace]PI3K[t], GABp\[LetterSpace]RasGAP[t], GABp\[LetterSpace]SHP2[t], GABp\[LetterSpace]pSHP2[t], GABp\[LetterSpace]pSHP2\[LetterSpace]GS[t], GS[t], I[t], IR[t], IRL[t], IRS[t], IRSp[t], IRSp\[LetterSpace]GS[t], IRSp\[LetterSpace]PI3K[t], IRSp\[LetterSpace]SHP2[t], IRi[t], IRp[t], IRp\[LetterSpace]IRS[t], IRp\[LetterSpace]IRSp[t], IRp\[LetterSpace]PI3K[t], IRp\[LetterSpace]RasGAP[t], Mek[t], Null[t], PDK1[t], PI3K[t], PIP3[t], PKA[t], R[t], RE[t], Raf[t], RasGAP[t], Rd[t], Ri[t], Rp[t], Rp\[LetterSpace]GS[t], Rp\[LetterSpace]PI3K[t], Rp\[LetterSpace]RasGAP[t], Rp\[LetterSpace]Shc[t], Rp\[LetterSpace]pShc[t], Rp\[LetterSpace]pShc\[LetterSpace]GS[t], SHP2[t], Shc[t], aPX[t], aPX\[LetterSpace]GS[t], aRaf[t], aSrc[t], aaRaf[t], amTOR[t], bRasGAP[t], dRas[t], iGS[t], iPX[t], iSrc[t], imGAB[t], imGABp[t], imIRS[t], mGAB[t], mGABp[t], mGABp\[LetterSpace]GS[t], mGABp\[LetterSpace]PI3K[t], mGABp\[LetterSpace]RasGAP[t], mGABp\[LetterSpace]SHP2[t], mGABp\[LetterSpace]pSHP2[t], mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t], mIRS[t], mIRSp[t], mIRSp\[LetterSpace]GS[t], mIRSp\[LetterSpace]PI3K[t], mIRSp\[LetterSpace]SHP2[t], mPDK1[t], mTOR[t], pAkt[t], pErk[t], pShc[t], pShc\[LetterSpace]GS[t], ppAkt[t], ppErk[t], ppMek[t], tRas[t], tRas\[LetterSpace]PI3K[t]
};

initialValues = {
Akt[0] == 100.0, EGF[0] == 1.0, Erk[0] == 400.0, GAB[0] == 225.0, GABp[0] == 0.0, GABp\[LetterSpace]GS[0] == 0.0, GABp\[LetterSpace]PI3K[0] == 0.0, GABp\[LetterSpace]RasGAP[0] == 0.0, GABp\[LetterSpace]SHP2[0] == 0.0, GABp\[LetterSpace]pSHP2[0] == 0.0, GABp\[LetterSpace]pSHP2\[LetterSpace]GS[0] == 0.0, GS[0] == 200.0, I[0] == 0.0, IR[0] == 150.0, IRL[0] == 0.0, IRS[0] == 300.0, IRSp[0] == 0.0, IRSp\[LetterSpace]GS[0] == 0.0, IRSp\[LetterSpace]PI3K[0] == 0.0, IRSp\[LetterSpace]SHP2[0] == 0.0, IRi[0] == 0.0, IRp[0] == 0.0, IRp\[LetterSpace]IRS[0] == 0.0, IRp\[LetterSpace]IRSp[0] == 0.0, IRp\[LetterSpace]PI3K[0] == 0.0, IRp\[LetterSpace]RasGAP[0] == 0.0, Mek[0] == 200.0, Null[0] == 0.0, PDK1[0] == 100.0, PI3K[0] == 200.0, PIP3[0] == 0.0, PKA[0] == 100.0, R[0] == 100.0, RE[0] == 0.0, Raf[0] == 100.0, RasGAP[0] == 50.0, Rd[0] == 0.0, Ri[0] == 0.0, Rp[0] == 0.0, Rp\[LetterSpace]GS[0] == 0.0, Rp\[LetterSpace]PI3K[0] == 0.0, Rp\[LetterSpace]RasGAP[0] == 0.0, Rp\[LetterSpace]Shc[0] == 0.0, Rp\[LetterSpace]pShc[0] == 0.0, Rp\[LetterSpace]pShc\[LetterSpace]GS[0] == 0.0, SHP2[0] == 300.0, Shc[0] == 270.0, aPX[0] == 0.0, aPX\[LetterSpace]GS[0] == 0.0, aRaf[0] == 0.0, aSrc[0] == 0.0, aaRaf[0] == 0.0, amTOR[0] == 0.0, bRasGAP[0] == 1*^-05, dRas[0] == 150.0, iGS[0] == 0.0, iPX[0] == 200.0, iSrc[0] == 518.0, imGAB[0] == 0.0, imGABp[0] == 0.0, imIRS[0] == 0.0, mGAB[0] == 0.0, mGABp[0] == 0.0, mGABp\[LetterSpace]GS[0] == 0.0, mGABp\[LetterSpace]PI3K[0] == 0.0, mGABp\[LetterSpace]RasGAP[0] == 0.0, mGABp\[LetterSpace]SHP2[0] == 0.0, mGABp\[LetterSpace]pSHP2[0] == 0.0, mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[0] == 0.0, mIRS[0] == 0.0, mIRSp[0] == 0.0, mIRSp\[LetterSpace]GS[0] == 0.0, mIRSp\[LetterSpace]PI3K[0] == 0.0, mIRSp\[LetterSpace]SHP2[0] == 0.0, mPDK1[0] == 0.0, mTOR[0] == 100.0, pAkt[0] == 0.0, pErk[0] == 0.0, pShc[0] == 0.0, pShc\[LetterSpace]GS[0] == 0.0, ppAkt[0] == 0.0, ppErk[0] == 0.0, ppMek[0] == 0.0, tRas[0] == 0.0, tRas\[LetterSpace]PI3K[0] == 0.0
};

rates = {
reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]100, reaction\[LetterSpace]101, reaction\[LetterSpace]102, reaction\[LetterSpace]103, reaction\[LetterSpace]104, reaction\[LetterSpace]105, reaction\[LetterSpace]106, reaction\[LetterSpace]107, reaction\[LetterSpace]108, reaction\[LetterSpace]109, reaction\[LetterSpace]11, reaction\[LetterSpace]110, reaction\[LetterSpace]111, reaction\[LetterSpace]112, reaction\[LetterSpace]113, reaction\[LetterSpace]114, reaction\[LetterSpace]115, reaction\[LetterSpace]117, reaction\[LetterSpace]118, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]17, reaction\[LetterSpace]18, reaction\[LetterSpace]19, reaction\[LetterSpace]2, reaction\[LetterSpace]20, 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]3, reaction\[LetterSpace]30, reaction\[LetterSpace]31, reaction\[LetterSpace]34, reaction\[LetterSpace]35, reaction\[LetterSpace]36, reaction\[LetterSpace]37, reaction\[LetterSpace]38, reaction\[LetterSpace]4, reaction\[LetterSpace]40, reaction\[LetterSpace]41, reaction\[LetterSpace]42, reaction\[LetterSpace]43, reaction\[LetterSpace]44, reaction\[LetterSpace]45, reaction\[LetterSpace]46, reaction\[LetterSpace]47, reaction\[LetterSpace]48, reaction\[LetterSpace]49, reaction\[LetterSpace]5, reaction\[LetterSpace]50, reaction\[LetterSpace]51, 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]63, reaction\[LetterSpace]64, reaction\[LetterSpace]65, reaction\[LetterSpace]66, reaction\[LetterSpace]67, reaction\[LetterSpace]68, reaction\[LetterSpace]69, reaction\[LetterSpace]7, 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]8, reaction\[LetterSpace]80, reaction\[LetterSpace]81, reaction\[LetterSpace]82, reaction\[LetterSpace]83, reaction\[LetterSpace]84, reaction\[LetterSpace]85, reaction\[LetterSpace]88, reaction\[LetterSpace]89, reaction\[LetterSpace]9, reaction\[LetterSpace]90, reaction\[LetterSpace]91, reaction\[LetterSpace]92, reaction\[LetterSpace]93, reaction\[LetterSpace]94, reaction\[LetterSpace]95, reaction\[LetterSpace]96, reaction\[LetterSpace]97, reaction\[LetterSpace]98, reaction\[LetterSpace]99
};

rateEquations = {
reaction\[LetterSpace]1 -> cell*(k1*EGF[t]*R[t] - k\[LetterSpace]1*RE[t]), reaction\[LetterSpace]10 -> cell*(-(k10*pShc\[LetterSpace]GS[t]*Rp[t]) + k\[LetterSpace]10*Rp\[LetterSpace]pShc\[LetterSpace]GS[t]), reaction\[LetterSpace]100 -> (cell*V31*IRSp\[LetterSpace]GS[t])/(Km31 + IRSp\[LetterSpace]GS[t]), reaction\[LetterSpace]101 -> cell*k48*IRSp\[LetterSpace]SHP2[t], reaction\[LetterSpace]102 -> cell*k56*mGABp\[LetterSpace]pSHP2[t], reaction\[LetterSpace]103 -> cell*k56*mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t], reaction\[LetterSpace]104 -> (cell*V51*GABp[t])/(Km51 + GABp[t]), reaction\[LetterSpace]105 -> (cell*V51*GABp\[LetterSpace]PI3K[t])/(Km51 + GABp\[LetterSpace]PI3K[t]), reaction\[LetterSpace]106 -> (cell*V51*GABp\[LetterSpace]GS[t])/(Km51 + GABp\[LetterSpace]GS[t]), reaction\[LetterSpace]107 -> (cell*V51*GABp\[LetterSpace]RasGAP[t])/(Km51 + GABp\[LetterSpace]RasGAP[t]), reaction\[LetterSpace]108 -> cell*k56*GABp\[LetterSpace]SHP2[t], reaction\[LetterSpace]109 -> cell*k56*GABp\[LetterSpace]pSHP2[t], reaction\[LetterSpace]11 -> cell*(-(k11*GS[t]*pShc[t]) + k\[LetterSpace]11*pShc\[LetterSpace]GS[t]), reaction\[LetterSpace]110 -> cell*k56*GABp\[LetterSpace]pSHP2\[LetterSpace]GS[t], reaction\[LetterSpace]111 -> cell*k111*mGABp\[LetterSpace]RasGAP[t]*(mGABp\[LetterSpace]pSHP2[t] + mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t] + mGABp\[LetterSpace]SHP2[t] + mIRSp\[LetterSpace]SHP2[t]), reaction\[LetterSpace]112 -> cell*k111*(mGABp\[LetterSpace]pSHP2[t] + mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t] + mGABp\[LetterSpace]SHP2[t])*Rp\[LetterSpace]RasGAP[t], reaction\[LetterSpace]113 -> cell*k111*IRp\[LetterSpace]RasGAP[t]*(mGABp\[LetterSpace]pSHP2[t] + mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t] + mGABp\[LetterSpace]SHP2[t]), reaction\[LetterSpace]114 -> cell*k111*mIRSp\[LetterSpace]SHP2[t]*Rp\[LetterSpace]RasGAP[t], reaction\[LetterSpace]115 -> cell*k111*IRp\[LetterSpace]RasGAP[t]*mIRSp\[LetterSpace]SHP2[t], reaction\[LetterSpace]117 -> cell*(-(k\[LetterSpace]80*imGABp[t]) + (2*kcat80*mGABp[t]*ppErk[t])/(Km80 + mGABp[t])), reaction\[LetterSpace]118 -> cell*k118*imGABp[t], reaction\[LetterSpace]12 -> cell*(k12*PI3K[t]*Rp[t] - k\[LetterSpace]12*Rp\[LetterSpace]PI3K[t]), reaction\[LetterSpace]13 -> cell*(k13*RasGAP[t]*Rp[t] - k\[LetterSpace]13*Rp\[LetterSpace]RasGAP[t]), reaction\[LetterSpace]17 -> cell*k17*Rp[t], reaction\[LetterSpace]18 -> cell*k17*Rp\[LetterSpace]GS[t], reaction\[LetterSpace]19 -> cell*k17*Rp\[LetterSpace]Shc[t], reaction\[LetterSpace]2 -> cell*(-(k\[LetterSpace]2*Rd[t]) + k2*RE[t]^2), reaction\[LetterSpace]20 -> cell*k17*Rp\[LetterSpace]pShc[t], reaction\[LetterSpace]21 -> cell*k17*Rp\[LetterSpace]pShc\[LetterSpace]GS[t], reaction\[LetterSpace]22 -> cell*k17*Rp\[LetterSpace]PI3K[t], reaction\[LetterSpace]23 -> cell*k17*Rp\[LetterSpace]RasGAP[t], reaction\[LetterSpace]24 -> cell*(k24*I[t]*IR[t] - k\[LetterSpace]24*IRL[t]), reaction\[LetterSpace]25 -> cell*k25*IRL[t], reaction\[LetterSpace]26 -> cell*(-(k\[LetterSpace]26*IRp\[LetterSpace]PI3K[t]) + k26*IRp[t]*PI3K[t]), reaction\[LetterSpace]27 -> cell*(-(k\[LetterSpace]27*IRp\[LetterSpace]RasGAP[t]) + k27*IRp[t]*RasGAP[t]), reaction\[LetterSpace]28 -> cell*(-(k\[LetterSpace]28*IRp\[LetterSpace]IRS[t]) + k28*IRp[t]*IRS[t]), reaction\[LetterSpace]29 -> cell*k29*IRp\[LetterSpace]IRS[t], reaction\[LetterSpace]3 -> cell*k3*Rd[t], reaction\[LetterSpace]30 -> cell*(k\[LetterSpace]30*IRp\[LetterSpace]IRSp[t] - k30*IRp[t]*IRSp[t]), reaction\[LetterSpace]31 -> (cell*V31*IRSp[t])/(Km31 + IRSp[t]), reaction\[LetterSpace]34 -> cell*k17*IRp[t], reaction\[LetterSpace]35 -> cell*k17*IRp\[LetterSpace]PI3K[t], reaction\[LetterSpace]36 -> cell*k17*IRp\[LetterSpace]RasGAP[t], reaction\[LetterSpace]37 -> cell*k17*IRp\[LetterSpace]IRS[t], reaction\[LetterSpace]38 -> cell*k17*IRp\[LetterSpace]IRSp[t], reaction\[LetterSpace]4 -> cell*(k4*GS[t]*Rp[t] - k\[LetterSpace]4*Rp\[LetterSpace]GS[t]), reaction\[LetterSpace]40 -> (cell*kcat40*iSrc[t]*(alpha40*IRp[t] + Rp[t]))/(Km40 + iSrc[t]), reaction\[LetterSpace]41 -> (cell*V41*aSrc[t])/(Km41 + aSrc[t]), reaction\[LetterSpace]42 -> cell*(-(k\[LetterSpace]42*mIRS[t]) + k42*IRS[t]*PIP3[t]), reaction\[LetterSpace]43 -> (cell*kcat43*mIRS[t]*(IRp[t] + alpha43*Rp[t]))/(Km43 + mIRS[t]), reaction\[LetterSpace]44 -> (cell*V31*mIRSp[t])/(Km31 + mIRSp[t]), reaction\[LetterSpace]45 -> cell*(k45*GS[t]*mIRSp[t] - k\[LetterSpace]45*mIRSp\[LetterSpace]GS[t]), reaction\[LetterSpace]46 -> cell*(-(k\[LetterSpace]46*mIRSp\[LetterSpace]PI3K[t]) + k46*mIRSp[t]*PI3K[t]), reaction\[LetterSpace]47 -> cell*(-(k\[LetterSpace]47*mIRSp\[LetterSpace]SHP2[t]) + k47*mIRSp[t]*SHP2[t]), reaction\[LetterSpace]48 -> cell*k48*mIRSp\[LetterSpace]SHP2[t], reaction\[LetterSpace]49 -> cell*(-(k\[LetterSpace]49*mGAB[t]) + k49*GAB[t]*PIP3[t]), reaction\[LetterSpace]5 -> cell*(-(k\[LetterSpace]5*Rp\[LetterSpace]Shc[t]) + k5*Rp[t]*Shc[t]), reaction\[LetterSpace]50 -> (cell*kcat50*mGAB[t]*(alpha50*aSrc[t] + Rp[t]))/(Km50 + mGAB[t]), reaction\[LetterSpace]51 -> (cell*V51*mGABp[t])/(Km51 + mGABp[t]), reaction\[LetterSpace]52 -> cell*(k52*GS[t]*mGABp[t] - k\[LetterSpace]52*mGABp\[LetterSpace]GS[t]), reaction\[LetterSpace]53 -> cell*(-(k\[LetterSpace]53*mGABp\[LetterSpace]PI3K[t]) + k53*mGABp[t]*PI3K[t]), reaction\[LetterSpace]54 -> cell*(-(k\[LetterSpace]54*mGABp\[LetterSpace]RasGAP[t]) + k54*mGABp[t]*RasGAP[t]), reaction\[LetterSpace]55 -> cell*(-(k\[LetterSpace]55*mGABp\[LetterSpace]SHP2[t]) + k55*mGABp[t]*SHP2[t]), reaction\[LetterSpace]56 -> cell*k56*mGABp\[LetterSpace]SHP2[t], reaction\[LetterSpace]57 -> (cell*kcat57*mGABp\[LetterSpace]SHP2[t]*(aSrc[t] + Rp[t]))/(Km57 + mGABp\[LetterSpace]SHP2[t]), reaction\[LetterSpace]58 -> (cell*V58*mGABp\[LetterSpace]pSHP2[t])/(Km58 + mGABp\[LetterSpace]pSHP2[t]), reaction\[LetterSpace]59 -> cell*(k59*GS[t]*mGABp\[LetterSpace]pSHP2[t] - k\[LetterSpace]59*mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t]), reaction\[LetterSpace]6 -> cell*k6*Rp\[LetterSpace]Shc[t], reaction\[LetterSpace]60 -> cell*k60*(IRp\[LetterSpace]PI3K[t] + mGABp\[LetterSpace]PI3K[t] + mIRSp\[LetterSpace]PI3K[t] + Rp\[LetterSpace]PI3K[t] + tRas\[LetterSpace]PI3K[t]), reaction\[LetterSpace]61 -> cell*k61*PIP3[t], reaction\[LetterSpace]62 -> (cell*kcat62*dRas[t]*(mGABp\[LetterSpace]GS[t] + mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t] + mIRSp\[LetterSpace]GS[t] + Rp\[LetterSpace]GS[t] + Rp\[LetterSpace]pShc\[LetterSpace]GS[t]))/(Km62 + dRas[t]), reaction\[LetterSpace]63 -> (cell*kcat63*(bRasGAP[t] + IRp\[LetterSpace]RasGAP[t] + mGABp\[LetterSpace]RasGAP[t] + Rp\[LetterSpace]RasGAP[t])*tRas[t])/(Km63 + tRas[t]), reaction\[LetterSpace]64 -> cell*(k64*PI3K[t]*tRas[t] - k\[LetterSpace]64*tRas\[LetterSpace]PI3K[t]), reaction\[LetterSpace]65 -> (cell*kcat65*Raf[t]*tRas[t])/(Km65 + Raf[t]), reaction\[LetterSpace]66 -> (cell*kcat66*aRaf[t]*aSrc[t])/(Km66 + aRaf[t]), reaction\[LetterSpace]67 -> cell*((kcat67*aaRaf[t]*PKA[t])/(Km67 + aaRaf[t]) + alpha67*aaRaf[t]*(pAkt[t] + beta67*ppAkt[t])), reaction\[LetterSpace]68 -> (cell*kcat68*aaRaf[t]*Mek[t])/(Km68 + Mek[t]), reaction\[LetterSpace]69 -> (cell*V69*ppMek[t])/(Km69 + ppMek[t]), reaction\[LetterSpace]7 -> cell*(-(k7*pShc[t]*Rp[t]) + k\[LetterSpace]7*Rp\[LetterSpace]pShc[t]), reaction\[LetterSpace]70 -> (cell*kcat70*Erk[t]*ppMek[t])/(Km70 + Erk[t] + (Km70*pErk[t])/Km71), reaction\[LetterSpace]71 -> (cell*kcat71*pErk[t]*ppMek[t])/(Km71 + (Km71*Erk[t])/Km70 + pErk[t]), reaction\[LetterSpace]72 -> (cell*V72*ppErk[t])/(Km72 + (Km72*pErk[t])/Km73 + ppErk[t]), reaction\[LetterSpace]73 -> (cell*V73*pErk[t])/(Km73 + pErk[t] + (Km73*ppErk[t])/Km72), reaction\[LetterSpace]74 -> cell*(-(k\[LetterSpace]74*mPDK1[t]) + k74*PDK1[t]*PIP3[t]), reaction\[LetterSpace]75 -> (cell*kcat75*Akt[t]*mPDK1[t])/(Km75 + Akt[t]), reaction\[LetterSpace]76 -> (cell*V76*pAkt[t])/(Km76 + pAkt[t]), reaction\[LetterSpace]77 -> cell*(-(k\[LetterSpace]77*amTOR[t]) + (kcat77*mTOR[t]*(alpha77*pAkt[t] + ppAkt[t]))/(Km77 + mTOR[t])), reaction\[LetterSpace]78 -> cell*((kcat78*amTOR[t]*pAkt[t])/(Km78 + pAkt[t]) - k\[LetterSpace]78*ppAkt[t]), reaction\[LetterSpace]79 -> cell*(-(k\[LetterSpace]79*iGS[t]) + (kcat79*GS[t]*ppErk[t])/(Km79 + GS[t])), reaction\[LetterSpace]8 -> (cell*V8*pShc[t])/(Km8 + pShc[t]), reaction\[LetterSpace]80 -> cell*(-(k\[LetterSpace]80*imGAB[t]) + (kcat80*mGAB[t]*ppErk[t])/(Km80 + mGAB[t])), reaction\[LetterSpace]81 -> cell*(-(k\[LetterSpace]81*imIRS[t]) + (kcat81*amTOR[t]*mIRS[t])/(Km81 + mIRS[t])), reaction\[LetterSpace]82 -> (cell*V\[LetterSpace]82*Rp[t])/(Km82 + Rp[t]), reaction\[LetterSpace]83 -> cell*k83*Ri[t], reaction\[LetterSpace]84 -> (cell*V\[LetterSpace]84*IRp[t])/(Km84 + IRp[t]), reaction\[LetterSpace]85 -> cell*k85*IRi[t], reaction\[LetterSpace]88 -> cell*(k\[LetterSpace]42*mIRSp[t] - k42*IRSp[t]*PIP3[t]), reaction\[LetterSpace]89 -> cell*(k\[LetterSpace]42*mIRSp\[LetterSpace]PI3K[t] - k42*IRSp\[LetterSpace]PI3K[t]*PIP3[t]), reaction\[LetterSpace]9 -> cell*(k9*GS[t]*Rp\[LetterSpace]pShc[t] - k\[LetterSpace]9*Rp\[LetterSpace]pShc\[LetterSpace]GS[t]), reaction\[LetterSpace]90 -> cell*(k\[LetterSpace]42*mIRSp\[LetterSpace]GS[t] - k42*IRSp\[LetterSpace]GS[t]*PIP3[t]), reaction\[LetterSpace]91 -> cell*(k\[LetterSpace]42*mIRSp\[LetterSpace]SHP2[t] - k42*IRSp\[LetterSpace]SHP2[t]*PIP3[t]), reaction\[LetterSpace]92 -> cell*(k\[LetterSpace]42*mGABp[t] - k42*GABp[t]*PIP3[t]), reaction\[LetterSpace]93 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]PI3K[t] - k42*GABp\[LetterSpace]PI3K[t]*PIP3[t]), reaction\[LetterSpace]94 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]GS[t] - k42*GABp\[LetterSpace]GS[t]*PIP3[t]), reaction\[LetterSpace]95 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]RasGAP[t] - k49*GABp\[LetterSpace]RasGAP[t]*PIP3[t]), reaction\[LetterSpace]96 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]SHP2[t] - k42*GABp\[LetterSpace]SHP2[t]*PIP3[t]), reaction\[LetterSpace]97 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]pSHP2[t] - k42*GABp\[LetterSpace]pSHP2[t]*PIP3[t]), reaction\[LetterSpace]98 -> cell*(k\[LetterSpace]42*mGABp\[LetterSpace]pSHP2\[LetterSpace]GS[t] - k42*GABp\[LetterSpace]pSHP2\[LetterSpace]GS[t]*PIP3[t]), reaction\[LetterSpace]99 -> (cell*V31*IRSp\[LetterSpace]PI3K[t])/(Km31 + IRSp\[LetterSpace]PI3K[t])
};

parameters = {
Kd1 -> 0.58824, Kd10 -> 400.0, Kd12 -> 12.45, Kd13 -> 200.0, Kd2 -> 15.0, Kd24 -> 0.029412, Kd26 -> 124.5, Kd27 -> 2000000.0, Kd28 -> 3.75, Kd30 -> 10.0, Kd4 -> 10.0, Kd42 -> 10.0, Kd45 -> 100000.0, Kd46 -> 1.0, Kd47 -> 1000.0, Kd49 -> 1.0, Kd5 -> 10.0, Kd52 -> 1.0, Kd53 -> 2.5, Kd54 -> 66666.0, Kd55 -> 100.0, Kd59 -> 20.0, Kd7 -> 400.0, Kd74 -> 100.0, Kd9 -> 10.0, Km31 -> 143.3, Km40 -> 110.0, Km41 -> 50.0, Km43 -> 150.0, Km50 -> 150.0, Km51 -> 130.0, Km57 -> 150.0, Km58 -> 130.0, Km62 -> 50.0, Km63 -> 50.0, Km65 -> 400.0, Km66 -> 10.0, Km67 -> 10000.0, Km68 -> 50.0, Km69 -> 675.299, Km70 -> 500.0, Km71 -> 500.0, Km72 -> 500.0, Km73 -> 500.0, Km75 -> 500.0, Km76 -> 1.0, Km77 -> 100.0, Km78 -> 100.0, Km79 -> 5000.0, Km8 -> 100.0, Km80 -> 700.0, Km81 -> 300.0, Km82 -> 50.0, Km84 -> 266.0, V31 -> 333.0, V41 -> 6.66, V51 -> 333.0, V58 -> 2.0, V69 -> 16.6, V72 -> 33.3, V73 -> 23.33, V76 -> 16.66, V8 -> 200.0, V\[LetterSpace]82 -> 133.0, V\[LetterSpace]84 -> 333.0, alpha40 -> 0.00025, alpha43 -> 0.05, alpha50 -> 0.0001, alpha67 -> 1*^-06, alpha77 -> 0.5, beta67 -> 2.0, k1 -> 0.068, k10 -> 0.0004, k111 -> 0.0133, k118 -> 0.001, k12 -> 0.00933, k13 -> 6.66*^-06, k17 -> 0.000185, k2 -> 0.033, k24 -> 0.011322, k25 -> 1.66, k26 -> 0.00933, k27 -> 6.66*^-08, k28 -> 0.1066, k29 -> 0.66, k3 -> 0.4, k30 -> 0.0066, k4 -> 0.000666, k42 -> 0.00666, k45 -> 0.000666, k46 -> 0.00666, k47 -> 0.000666, k48 -> 0.666, k49 -> 0.000666, k5 -> 0.0133, k52 -> 0.002, k53 -> 0.0133, k54 -> 1*^-05, k55 -> 0.000666, k56 -> 0.666, k59 -> 0.01, k6 -> 0.333, k60 -> 4.66, k61 -> 3.33, k64 -> 0.0, k7 -> 0.000666, k74 -> 0.00666, k83 -> 0.0166, k85 -> 0.0166, k9 -> 0.00666, k\[LetterSpace]64 -> 2.5, k\[LetterSpace]77 -> 0.666, k\[LetterSpace]78 -> 0.666, k\[LetterSpace]79 -> 6.66*^-05, k\[LetterSpace]80 -> 6.66*^-05, k\[LetterSpace]81 -> 6.66*^-05, kcat40 -> 6.6, kcat43 -> 33.3, kcat50 -> 3333.0, kcat57 -> 0.133, kcat62 -> 5.33, kcat63 -> 20000.0, kcat65 -> 0.1, kcat66 -> 3.33, kcat67 -> 0.666, kcat68 -> 0.133, kcat70 -> 0.333, kcat71 -> 0.666, kcat75 -> 4.66, kcat77 -> 0.666, kcat78 -> 0.666, kcat79 -> 0.0466, kcat80 -> 0.04, kcat81 -> 0.166, cell -> 1.0, extra -> 34.0
};

assignments = {
k\[LetterSpace]5 -> k5*Kd5, k\[LetterSpace]4 -> k4*Kd4, k\[LetterSpace]2 -> k2*Kd2, k\[LetterSpace]1 -> k1*Kd1, k\[LetterSpace]10 -> k10*Kd10, k\[LetterSpace]9 -> k9*Kd9, k\[LetterSpace]7 -> k7*Kd7, k\[LetterSpace]74 -> k74*Kd74, k\[LetterSpace]59 -> k59*Kd59, k\[LetterSpace]55 -> k55*Kd55, k\[LetterSpace]54 -> k54*Kd54, k\[LetterSpace]53 -> k53*Kd53, k\[LetterSpace]52 -> k52*Kd52, k\[LetterSpace]49 -> k49*Kd49, k\[LetterSpace]47 -> k47*Kd47, k\[LetterSpace]46 -> k46*Kd46, k\[LetterSpace]45 -> k45*Kd45, k\[LetterSpace]42 -> k42*Kd42, k\[LetterSpace]30 -> k30*Kd30, k\[LetterSpace]28 -> k28*Kd28, k\[LetterSpace]27 -> k27*Kd27, k\[LetterSpace]26 -> k26*Kd26, k\[LetterSpace]24 -> k24*Kd24, k\[LetterSpace]13 -> k13*Kd13, k\[LetterSpace]12 -> k12*Kd12, k11 -> k9, k\[LetterSpace]11 -> k\[LetterSpace]9, EGF\[LetterSpace]tot -> EGF[t] + (cell*(RE[t] + 2*(Rd[t] + Ri[t] + Rp[t] + Rp\[LetterSpace]GS[t] + Rp\[LetterSpace]PI3K[t] + Rp\[LetterSpace]pShc[t] + Rp\[LetterSpace]pShc\[LetterSpace]GS[t] + Rp\[LetterSpace]RasGAP[t] + Rp\[LetterSpace]Shc[t])))/extra, phosphorylated\[LetterSpace]Akt -> pAkt[t] + ppAkt[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Akt'[t] == 1.0*reaction\[LetterSpace]76 -1.0*reaction\[LetterSpace]75, EGF'[t] ==  -1.0*reaction\[LetterSpace]1, Erk'[t] == 1.0*reaction\[LetterSpace]73 -1.0*reaction\[LetterSpace]70, GAB'[t] == 1.0*reaction\[LetterSpace]104 +1.0*reaction\[LetterSpace]105 +1.0*reaction\[LetterSpace]106 +1.0*reaction\[LetterSpace]107 +1.0*reaction\[LetterSpace]108 +1.0*reaction\[LetterSpace]109 +1.0*reaction\[LetterSpace]110 -1.0*reaction\[LetterSpace]49, GABp'[t] == 1.0*reaction\[LetterSpace]92 -1.0*reaction\[LetterSpace]104, GABp\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]94 -1.0*reaction\[LetterSpace]106, GABp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]93 -1.0*reaction\[LetterSpace]105, GABp\[LetterSpace]RasGAP'[t] == 1.0*reaction\[LetterSpace]95 -1.0*reaction\[LetterSpace]107, GABp\[LetterSpace]SHP2'[t] == 1.0*reaction\[LetterSpace]96 -1.0*reaction\[LetterSpace]108, GABp\[LetterSpace]pSHP2'[t] == 1.0*reaction\[LetterSpace]97 -1.0*reaction\[LetterSpace]109, GABp\[LetterSpace]pSHP2\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]98 -1.0*reaction\[LetterSpace]110, GS'[t] == 1.0*reaction\[LetterSpace]11 +1.0*reaction\[LetterSpace]18 +1.0*reaction\[LetterSpace]100 +1.0*reaction\[LetterSpace]103 +1.0*reaction\[LetterSpace]106 +1.0*reaction\[LetterSpace]110 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]79, I'[t] ==  -1.0*reaction\[LetterSpace]24, IR'[t] ==  -1.0*reaction\[LetterSpace]24, IRL'[t] == 1.0*reaction\[LetterSpace]24 +1.0*reaction\[LetterSpace]85 -1.0*reaction\[LetterSpace]25, IRS'[t] == 1.0*reaction\[LetterSpace]31 +1.0*reaction\[LetterSpace]37 +1.0*reaction\[LetterSpace]99 +1.0*reaction\[LetterSpace]100 +1.0*reaction\[LetterSpace]101 -1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]42, IRSp'[t] == 1.0*reaction\[LetterSpace]30 +1.0*reaction\[LetterSpace]38 +1.0*reaction\[LetterSpace]88 -1.0*reaction\[LetterSpace]31, IRSp\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]90 -1.0*reaction\[LetterSpace]100, IRSp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]89 -1.0*reaction\[LetterSpace]99, IRSp\[LetterSpace]SHP2'[t] == 1.0*reaction\[LetterSpace]91 -1.0*reaction\[LetterSpace]101, IRi'[t] == 1.0*reaction\[LetterSpace]84 -1.0*reaction\[LetterSpace]85, IRp'[t] == 1.0*reaction\[LetterSpace]25 +1.0*reaction\[LetterSpace]30 +1.0*reaction\[LetterSpace]113 +1.0*reaction\[LetterSpace]115 -1.0*reaction\[LetterSpace]26 -1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]34 -1.0*reaction\[LetterSpace]84, IRp\[LetterSpace]IRS'[t] == 1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]37, IRp\[LetterSpace]IRSp'[t] == 1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]38, IRp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]26 -1.0*reaction\[LetterSpace]35, IRp\[LetterSpace]RasGAP'[t] == 1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]36 -1.0*reaction\[LetterSpace]113 -1.0*reaction\[LetterSpace]115, Mek'[t] == 1.0*reaction\[LetterSpace]69 -1.0*reaction\[LetterSpace]68, Null'[t] == 1.0*reaction\[LetterSpace]17 +1.0*reaction\[LetterSpace]34 +1.0*reaction\[LetterSpace]61 , PDK1'[t] ==  -1.0*reaction\[LetterSpace]74, PI3K'[t] == 1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]35 +1.0*reaction\[LetterSpace]99 +1.0*reaction\[LetterSpace]105 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]26 -1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]64, PIP3'[t] == 1.0*reaction\[LetterSpace]60 +1.0*reaction\[LetterSpace]88 +1.0*reaction\[LetterSpace]89 +1.0*reaction\[LetterSpace]90 +1.0*reaction\[LetterSpace]91 +1.0*reaction\[LetterSpace]92 +1.0*reaction\[LetterSpace]93 +1.0*reaction\[LetterSpace]94 +1.0*reaction\[LetterSpace]95 +1.0*reaction\[LetterSpace]96 +1.0*reaction\[LetterSpace]97 +1.0*reaction\[LetterSpace]98 -1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]49 -1.0*reaction\[LetterSpace]61 -1.0*reaction\[LetterSpace]74, PKA'[t] == 0.0 , R'[t] ==  -1.0*reaction\[LetterSpace]1, RE'[t] == 1.0*reaction\[LetterSpace]1 -2.0*reaction\[LetterSpace]2, Raf'[t] == 1.0*reaction\[LetterSpace]67 -1.0*reaction\[LetterSpace]65, RasGAP'[t] == 1.0*reaction\[LetterSpace]23 +1.0*reaction\[LetterSpace]36 +1.0*reaction\[LetterSpace]107 +1.0*reaction\[LetterSpace]111 +1.0*reaction\[LetterSpace]112 +1.0*reaction\[LetterSpace]113 +1.0*reaction\[LetterSpace]114 +1.0*reaction\[LetterSpace]115 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]54, Rd'[t] == 1.0*reaction\[LetterSpace]2 +1.0*reaction\[LetterSpace]83 -1.0*reaction\[LetterSpace]3, Ri'[t] == 1.0*reaction\[LetterSpace]82 -1.0*reaction\[LetterSpace]83, Rp'[t] == 1.0*reaction\[LetterSpace]3 +1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]112 +1.0*reaction\[LetterSpace]114 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]17 -1.0*reaction\[LetterSpace]82, Rp\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]18, Rp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]22, Rp\[LetterSpace]RasGAP'[t] == 1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]112 -1.0*reaction\[LetterSpace]114, Rp\[LetterSpace]Shc'[t] == 1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]19, Rp\[LetterSpace]pShc'[t] == 1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]20, Rp\[LetterSpace]pShc\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]21, SHP2'[t] == 1.0*reaction\[LetterSpace]48 +1.0*reaction\[LetterSpace]56 +1.0*reaction\[LetterSpace]101 +1.0*reaction\[LetterSpace]102 +1.0*reaction\[LetterSpace]103 +1.0*reaction\[LetterSpace]108 +1.0*reaction\[LetterSpace]109 +1.0*reaction\[LetterSpace]110 -1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]55, Shc'[t] == 1.0*reaction\[LetterSpace]8 +1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]5, aPX'[t] == 0.0 , aPX\[LetterSpace]GS'[t] == 0.0 , aRaf'[t] == 1.0*reaction\[LetterSpace]65 -1.0*reaction\[LetterSpace]66, aSrc'[t] == 1.0*reaction\[LetterSpace]40 -1.0*reaction\[LetterSpace]41, aaRaf'[t] == 1.0*reaction\[LetterSpace]66 -1.0*reaction\[LetterSpace]67, amTOR'[t] == 1.0*reaction\[LetterSpace]77 , bRasGAP'[t] == 0.0 , dRas'[t] == 1.0*reaction\[LetterSpace]63 -1.0*reaction\[LetterSpace]62, iGS'[t] == 1.0*reaction\[LetterSpace]79 , iPX'[t] == 0.0 , iSrc'[t] == 1.0*reaction\[LetterSpace]41 -1.0*reaction\[LetterSpace]40, imGAB'[t] == 1.0*reaction\[LetterSpace]80 +1.0*reaction\[LetterSpace]118 , imGABp'[t] == 1.0*reaction\[LetterSpace]117 -1.0*reaction\[LetterSpace]118, imIRS'[t] == 1.0*reaction\[LetterSpace]81 , mGAB'[t] == 1.0*reaction\[LetterSpace]49 +1.0*reaction\[LetterSpace]51 +1.0*reaction\[LetterSpace]56 +1.0*reaction\[LetterSpace]102 +1.0*reaction\[LetterSpace]103 +1.0*reaction\[LetterSpace]111 -1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]80, mGABp'[t] == 1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]51 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]55 -1.0*reaction\[LetterSpace]92 -1.0*reaction\[LetterSpace]117, mGABp\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]94, mGABp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]93, mGABp\[LetterSpace]RasGAP'[t] == 1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]95 -1.0*reaction\[LetterSpace]111, mGABp\[LetterSpace]SHP2'[t] == 1.0*reaction\[LetterSpace]55 +1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]56 -1.0*reaction\[LetterSpace]57 -1.0*reaction\[LetterSpace]96, mGABp\[LetterSpace]pSHP2'[t] == 1.0*reaction\[LetterSpace]57 -1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]97 -1.0*reaction\[LetterSpace]102, mGABp\[LetterSpace]pSHP2\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]98 -1.0*reaction\[LetterSpace]103, mIRS'[t] == 1.0*reaction\[LetterSpace]42 +1.0*reaction\[LetterSpace]44 +1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]43 -1.0*reaction\[LetterSpace]81, mIRSp'[t] == 1.0*reaction\[LetterSpace]43 -1.0*reaction\[LetterSpace]44 -1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]88, mIRSp\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]90, mIRSp\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]89, mIRSp\[LetterSpace]SHP2'[t] == 1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]91, mPDK1'[t] == 1.0*reaction\[LetterSpace]74 , mTOR'[t] ==  -1.0*reaction\[LetterSpace]77, pAkt'[t] == 1.0*reaction\[LetterSpace]75 -1.0*reaction\[LetterSpace]76 -1.0*reaction\[LetterSpace]78, pErk'[t] == 1.0*reaction\[LetterSpace]70 +1.0*reaction\[LetterSpace]72 -1.0*reaction\[LetterSpace]71 -1.0*reaction\[LetterSpace]73, pShc'[t] == 1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]11 +1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]8, pShc\[LetterSpace]GS'[t] == 1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]11, ppAkt'[t] == 1.0*reaction\[LetterSpace]78 , ppErk'[t] == 1.0*reaction\[LetterSpace]71 -1.0*reaction\[LetterSpace]72, ppMek'[t] == 1.0*reaction\[LetterSpace]68 -1.0*reaction\[LetterSpace]69, tRas'[t] == 1.0*reaction\[LetterSpace]62 -1.0*reaction\[LetterSpace]63 -1.0*reaction\[LetterSpace]64, tRas\[LetterSpace]PI3K'[t] == 1.0*reaction\[LetterSpace]64 
};
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]}]