{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simple EGFR model\n",
    "\n",
    "See http://bionetgen.org/index.php/Simple_EGFR_model for details."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from ecell4.prelude import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/markdown": [
       "Blinov ML, Faeder JR, Goldstein B, Hlavacek WS, A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity. *Bio Systems*, **2-3**(83), 136-51, 2006. [10.1016/j.biosystems.2005.06.014](https://doi.org/10.1016/j.biosystems.2005.06.014). PubMed PMID: [16233948](https://www.ncbi.nlm.nih.gov/pubmed/16233948)."
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "citation(16233948)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "NA = 6.02e23  # Avogadro's number (molecules/mol)\n",
    "f = 1  # Fraction of the cell to simulate\n",
    "Vo = f * 1.0e-10  # Extracellular volume=1/cell_density (L)\n",
    "V = f * 3.0e-12  # Cytoplasmic volume (L)\n",
    "\n",
    "EGF_init = 20 * 1e-9 * NA * Vo  # Initial amount of ligand (20 nM) converted to copies per cell\n",
    "\n",
    "# Initial amounts of cellular components (copies per cell)\n",
    "EGFR_init = f * 1.8e5\n",
    "Grb2_init = f * 1.5e5\n",
    "Sos1_init = f * 6.2e4\n",
    "\n",
    "# Rate constants\n",
    "# Divide by NA*V to convert bimolecular rate constants\n",
    "# from /M/sec to /(molecule/cell)/sec\n",
    "kp1 = 9.0e7 / (NA * Vo)  # ligand-monomer binding\n",
    "km1 = 0.06  # ligand-monomer dissociation\n",
    "kp2 = 1.0e7 / (NA * V)  # aggregation of bound monomers\n",
    "km2 = 0.1  # dissociation of bound monomers\n",
    "kp3 = 0.5  # dimer transphosphorylation\n",
    "km3 = 4.505  # dimer dephosphorylation\n",
    "kp4 = 1.5e6 / (NA * V)  # binding of Grb2 to receptor\n",
    "km4 = 0.05  # dissociation of Grb2 from receptor\n",
    "kp5 = 1.0e7 / (NA * V)  # binding of Grb2 to Sos1\n",
    "km5 = 0.06  # dissociation of Grb2 from Sos1\n",
    "deg = 0.01  # degradation of receptor dimers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "with reaction_rules():\n",
    "  # R1: Ligand-receptor binding\n",
    "  EGFR(L, CR1) + EGF(R) == EGFR(L^1, CR1).EGF(R^1) | (kp1, km1)\n",
    "  \n",
    "  # R2: Receptor-aggregation\n",
    "  EGFR(L^_, CR1) + EGFR(L^_, CR1) == EGFR(L^_,CR1^1).EGFR(L^_,CR1^1) | (kp2, km2)\n",
    "  \n",
    "  # R3: Transphosphorylation of EGFR by RTK\n",
    "  EGFR(CR1^_, Y1068=U) > EGFR(CR1^_, Y1068=P) | kp3\n",
    "  \n",
    "  # R4: Dephosphorylation\n",
    "  EGFR(Y1068=P) > EGFR(Y1068=U) | km3\n",
    "  \n",
    "  # R5: Grb2 binding to pY1068\n",
    "  EGFR(Y1068=P) + Grb2(SH2) == EGFR(Y1068=P^1).Grb2(SH2^1) | (kp4, km4)\n",
    "  \n",
    "  # R6: Grb2 binding to Sos1\n",
    "  Grb2(SH3) + Sos1(PxxP) == Grb2(SH3^1).Sos1(PxxP^1) | (kp5, km5)\n",
    "  \n",
    "  # R7: Receptor dimer internalization/degradation\n",
    "  (EGF(R^1).EGF(R^2).EGFR(L^1,CR1^3).EGFR(L^2,CR1^3) > ~EmptySet | deg\n",
    "       | ReactionRule.STRICT | ReactionRule.DESTROY)\n",
    "\n",
    "m1 = get_model(is_netfree=True, effective=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "EGFR(L, CR1) + EGF(R) > EGFR(L ^ 1, CR1).EGF(R ^ 1) | 1.495016611295681e-06\n",
      "EGFR(L ^ 1, CR1).EGF(R ^ 1) > EGFR(L, CR1) + EGF(R) | 0.06\n",
      "EGFR(L ^ _, CR1) + EGFR(L ^ _, CR1) > EGFR(L ^ _, CR1 ^ 1).EGFR(L ^ _, CR1 ^ 1) | 5.537098560354375e-06\n",
      "EGFR(L ^ _, CR1 ^ 1).EGFR(L ^ _, CR1 ^ 1) > EGFR(L ^ _, CR1) + EGFR(L ^ _, CR1) | 0.1\n",
      "EGFR(CR1 ^ _, Y1068=U) > EGFR(CR1 ^ _, Y1068=P) | 0.5\n",
      "EGFR(Y1068=P) > EGFR(Y1068=U) | 4.505\n",
      "EGFR(Y1068=P) + Grb2(SH2) > EGFR(Y1068=P ^ 1).Grb2(SH2 ^ 1) | 8.305647840531562e-07\n",
      "EGFR(Y1068=P ^ 1).Grb2(SH2 ^ 1) > EGFR(Y1068=P) + Grb2(SH2) | 0.05\n",
      "Grb2(SH3) + Sos1(PxxP) > Grb2(SH3 ^ 1).Sos1(PxxP ^ 1) | 5.537098560354375e-06\n",
      "Grb2(SH3 ^ 1).Sos1(PxxP ^ 1) > Grb2(SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R ^ 1).EGF(R ^ 2).EGFR(L ^ 1, CR1 ^ 3).EGFR(L ^ 2, CR1 ^ 3) >  | 0.01\n"
     ]
    }
   ],
   "source": [
    "show(m1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "y0 = {\"EGF(R)\": EGF_init, \"EGFR(L,CR1,Y1068=U)\": EGFR_init, \"Grb2(SH2,SH3)\": Grb2_init, \"Sos1(PxxP)\": Sos1_init}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "m2 = m1.expand([Species(serial) for serial in y0.keys()])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "22 species and 86 reactions were generated.\n"
     ]
    }
   ],
   "source": [
    "print(\"{} species and {} reactions were generated.\".format(len(m2.list_species()), len(m2.reaction_rules())))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) | {}\n",
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) | {}\n",
      "EGFR(CR1,L,Y1068=P) | {}\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | {}\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) | {}\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | {}\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | {}\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=U) | {}\n",
      "Grb2(SH2,SH3^1).Sos1(PxxP^1) | {}\n",
      "EGF(R) | {}\n",
      "EGFR(CR1,L,Y1068=U) | {}\n",
      "Grb2(SH2,SH3) | {}\n",
      "Sos1(PxxP) | {}\n",
      "\n",
      "EGF(R) + EGFR(CR1,L,Y1068=U) > EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 1.495016611295681e-06\n",
      "Grb2(SH2,SH3) + Sos1(PxxP) > Grb2(SH2,SH3^1).Sos1(PxxP^1) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=U) > EGF(R) + EGFR(CR1,L,Y1068=U) | 0.06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=U) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | 2.7685492801771873e-06\n",
      "Grb2(SH2,SH3^1).Sos1(PxxP^1) > Grb2(SH2,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) > EGF(R^1).EGFR(CR1,L^1,Y1068=U) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | 1.0\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) | 0.5\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | 4.505\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) + Grb2(SH2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) > EGF(R) + EGFR(CR1,L,Y1068=P) | 0.06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) | 2.7685492801771873e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) | 0.5\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) | 0.5\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) > EGF(R^1).EGFR(CR1,L^1,Y1068=U) | 4.505\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) | 9.01\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + Grb2(SH2,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) | 1.6611295681063123e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) + Grb2(SH2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) | 1.6611295681063123e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3) + Grb2(SH2,SH3) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) + Sos1(PxxP) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) >  | 0.01\n",
      "EGF(R) + EGFR(CR1,L,Y1068=P) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) | 1.495016611295681e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R) + EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) | 0.06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) > EGF(R) + EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) | 0.06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) | 2.7685492801771873e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) | 2.7685492801771873e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=U) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) | 0.1\n",
      "EGFR(CR1,L,Y1068=P) > EGFR(CR1,L,Y1068=U) | 4.505\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5) | 4.505\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3) | 4.505\n",
      "EGFR(CR1,L,Y1068=P) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) | 8.305647840531562e-07\n",
      "EGFR(CR1,L,Y1068=P) + Grb2(SH2,SH3) > EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) + Grb2(SH2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Grb2(SH2,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) | 8.305647840531562e-07\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.05\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P) + Grb2(SH2,SH3) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3) + Grb2(SH2,SH3) | 0.05\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + Sos1(PxxP) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Sos1(PxxP) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) >  | 0.01\n",
      "EGF(R) + EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 1.495016611295681e-06\n",
      "EGF(R) + EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) | 1.495016611295681e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) > EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) + EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3) | 0.1\n",
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) > EGFR(CR1,L,Y1068=P) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.05\n",
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) > EGFR(CR1,L,Y1068=P) + Grb2(SH2,SH3) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.1\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Grb2(SH2,SH3^1).Sos1(PxxP^1) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5) + Grb2(SH2,SH3) | 0.05\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3) + Grb2(SH2,SH3) | 0.1\n",
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) + Sos1(PxxP) > EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) + Sos1(PxxP) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) | 5.537098560354375e-06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) + Sos1(PxxP) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) | 1.107419712070875e-05\n",
      "EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2) > EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) + Sos1(PxxP) | 0.12\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) > EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) + Sos1(PxxP) | 0.06\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6) >  | 0.01\n",
      "EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3) >  | 0.01\n"
     ]
    }
   ],
   "source": [
    "show(m2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "ret = run_simulation(120, model=m2, y0=y0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEGCAYAAABVSfMhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de5xkZX3n8c+vbn2f6QEGgRlwUEcNKN5awOjLdcXLYHwxJBKDLzcSw8qa1YgbdwPEZFkTk5BdV4JZZYNAxKzLRSAL8RJERKOJCIMmXEVGQGa4Dsy1r3X77R/Pc6pO9VR3V/d0VU9Xf9++jvWc33nOeZ4zRdevzjlPnWPujoiISKdklroDIiKysijxiIhIRynxiIhIRynxiIhIRynxiIhIR+WWugMHu8MOO8w3bNiw1N0QEVlW7r777ufcfW2zZUo8c9iwYQNbtmxZ6m6IiCwrZvaLmZbpVJuIiHSUEo+IiHSUEo+IiHSUEo+IiHSUEo+IiHSUEo+IiHRU2xKPmV1pZs+a2X2p2P8ws5+a2T1m9ndmNpxadoGZbTWzh8zsnan4phjbambnp+LHmtmPzOxhM7vWzAox3hPnt8blG+ZqQ0REOqedRzxfAjZNi90KvMLdTwB+BlwAYGbHAWcCx8d1vmBmWTPLAp8HTgWOA94X6wL8BXCxu28EdgFnx/jZwC53fwlwcaw3YxuLvdM1z9wP3/4UTOxuWxMiIstR2xKPu/8jsHNa7FvuXo6zdwDrY3kzcI27T7n7o8BW4MQ4bXX3R9y9CFwDbDYzA94KXB/Xvwo4PbWtq2L5euCUWH+mNtpj12Pwg8/Czkfa1oSIyHK0lNd4fhv4ZiyvA7allm2PsZnihwK7U0ksiTdsKy7fE+vPtK39mNk5ZrbFzLbs2LFjQTvH6qPD655ts9cTEVlhliTxmNkngTLwlSTUpJovIL6Qbe0fdL/M3UfcfWTt2qa3GprbcEw8u5V4RETSOn6vNjM7C3g3cIrXn7u9HTg6VW098GQsN4s/BwybWS4e1aTrJ9vabmY5YDXhlN9sbSy+3mEoDMHux9vWhIjIctTRIx4z2wScB5zm7uOpRTcDZ8YRaccCG4E7gbuAjXEEW4EwOODmmLBuB86I658F3JTa1lmxfAbwnVh/pjbawywc9ehUm4hIg7Yd8ZjZ1cBbgMPMbDtwIWEUWw9wa7jezx3u/mF3v9/MrgMeIJyC+4i7V+J2PgrcAmSBK939/tjEecA1ZvZp4CfAFTF+BfC3ZraVcKRzJsBsbbTN8DE61SYiMo3Vz3ZJMyMjI77gxyJ8/T/DvdfB+TrdJiIri5nd7e4jzZbpzgXtNHw0TO4Jk4iIAEo87bVaI9tERKZT4mmn4ReGVw0wEBGpUeJpJ/2WR0RkP0o87TSwFnK9sHvGR4+LiKw4SjztZAar1+tUm4hIihJPu+m3PCIiDZR42m217l4gIpKmxNNuw0fD2A4ojs9dV0RkBVDiabfVx4TXPduXth8iIgcJJZ52S4ZU79Ftc0REQImn/YbjEY8GGIiIAEo87Td0JGRyei6PiEikxNNumSysOkoj20REIiWeTlit3/KIiCSUeDph+Bgd8YiIREo8nTB8NOx9EsrFpe6JiMiSU+LphNVHA66jHhERlHg6Y91rwTLwzfOgUlrq3oiILCklnk54wfHw7oth663w9x8H96XukYjIksktdQdWjNf9Fux9Cr53EQwdAf/2k5BR3heRlUeJp5Pecj7sexK+/xm484tw1KvgyFfDwGFQGAxTNh9+cJrNh9NzWHg1YtmIM4Glyul4ms0Qn6n+rOu0uP5ibOdg6EM7tzWv6h3ez3lts13bbVdfO7X9TrTR5u33rQmfT4tMiaeTzOBXLoYXvgm23QFP/BjuuBSquu4jIgehN34c3v6pRd+sEk+nZXPwqt8IE0C1CqVxKI5CcQwqRaiW4yAED9eDvBqvC/m060Op8ozXjWaIz3qdqYVrUC1dpzqQ7XSyDy1spuVtzff63cHwb30A22zXdtvV105tvxNtdOJa8WEvbctm25Z4zOxK4N3As+7+ihg7BLgW2AA8BrzX3XeZmQGXAO8CxoHfcvcfx3XOAv4wbvbT7n5VjL8O+BLQB3wDONfdfSFtLKlMBnoGwyQisgK08+r2l4BN02LnA7e5+0bgtjgPcCqwMU7nAJdCLVFdCJwEnAhcaGZr4jqXxrrJepsW0oaIiHRW2xKPu/8jsHNaeDNwVSxfBZyein/ZgzuAYTM7EngncKu773T3XcCtwKa4bJW7/9DdHfjytG3Npw0REemgTo/nfYG7PwUQXw+P8XVA+mf922Nstvj2JvGFtLEfMzvHzLaY2ZYdO3bMawdFRGR2B8sPSZqN8fMFxBfSxv5B98vcfcTdR9auXTvHZkVEZD46nXieSU5vxddnY3w7cHSq3nrgyTni65vEF9KGiIh0UKcTz83AWbF8FnBTKv4BC04G9sTTZLcA7zCzNXFQwTuAW+KyfWZ2chyt9oFp25pPGyIi0kFzJh4z+3UzG4rlPzSzG83stS2sdzXwQ+BlZrbdzM4GLgLebmYPA2+P8xCGQz8CbAW+CPxHAHffCfwJcFec/jjGAH4HuDyu83PgmzE+rzZERKSzzOf4EZKZ3ePuJ5jZm4A/Bz4D/IG7n9SJDi61kZER37Jly1J3Q0RkWTGzu919pNmyVk61VeLrrwCXuvtNQGGxOiciIitLK4nnCTP7a+C9wDfMrKfF9URERPbTSgJ5L+Ei/yZ33w0cAvyXtvZKRES61pyJx93HCUOS3xRDZeDhdnZKRES6Vyuj2i4EzgMuiKE88H/a2SkREelerZxq+1XgNGAMwN2fBIba2SkREelerSSeYrwRpwOY2UB7uyQiIt2slcRzXRzVNmxmHwK+TfgBpoiIyLzN+SA4d/+Mmb0d2Au8DPiv7n5r23smIiJdqaUnkMZEo2QjIiIHbMbEY2b7qD+CIH1fHQPc3Ve1uW8iItKFZkw87q6RayIisujmPNVmZsc0i7v744vfHRER6XatXOP5eqrcCxwLPAQc35YeiYhIV2tlVNsr0/PxWTz/oW09EhGRrjbvu0y7+4+B17ehLyIisgK0co3n91KzGeC1wI629UhERLpaK9d40qPbyoRrPje0pzsiItLtWrnG86lOdERERFaGVh6LcKuZDafm15jZLe3tloiIdKtWBhesjU8eBcDddwGHt69LIiLSzVpJPJX0j0jN7IU03kJHRESkZa0MLvgk8AMz+16cfzNwTvu6JCIi3WzOIx53/wfCEOprgeuA17n7AV3jMbP/ZGb3m9l9Zna1mfWa2bFm9iMze9jMrjWzQqzbE+e3xuUbUtu5IMYfMrN3puKbYmyrmZ2fijdtQ0REOqeVwQUGbAJe6+5/D/Sb2YkLbdDM1gEfA0bc/RVAFjgT+AvgYnffCOwCzo6rnA3scveXABfHepjZcXG942P/vmBmWTPLAp8HTgWOA94X6zJLGyIi0iGtXOP5AvAG4H1xfh/hg/1A5IA+M8sB/cBTwFuB6+Pyq4DTY3lznCcuPyUmw83ANe4+5e6PAluBE+O01d0fcfcicA2wOa4zUxsiItIhrSSek9z9I8Ak1Ea1LfgUlbs/AXwGeJyQcPYAdwO73b0cq20H1sXyOmBbXLcc6x+ajk9bZ6b4obO0ISIiHdJK4inF01cOYGZrgepCGzSzNYSjlWOBo4ABwmmx6ZKRczbDssWKN+vjOWa2xcy27NihuwOJiCymVhLP54C/Aw43sz8FfgD82QG0+TbgUXff4e4l4Ebgl4HheOoNYD3wZCxvB44GiMtXAzvT8WnrzBR/bpY2Grj7Ze4+4u4ja9euPYBdFRGR6VoZ1fYV4PeBPyecGjvd3b96AG0+DpxsZv3xusspwAPA7cAZsc5ZwE2xfHOcJy7/jrt7jJ8ZR70dC2wE7gTuAjbGEWwFwgCEm+M6M7UhIiIdMuPveMzskNTss8DV6WXuvnMhDbr7j8zseuDHhJuO/gS4jHDz0WvM7NMxdkVc5Qrgb81sK+FI58y4nfvN7DpC0ioDH3H3SuzfR4FbCCPmrnT3++O2zpuhDRER6RALBwJNFpg9yizXRtz9Re3s2MFiZGTEt2zZstTdEBFZVszsbncfabZsxiMedz+2fV0SEZGVqpVb5mBmpxFulQPwXXf/Wvu6JCIi3ayVOxdcBJxLuJbyAHCumf15uzsmIiLdqZUjnncBr3b3KoCZXUW4MH9BOzsmIiLdqZXf8QAMp8qr29ERERFZGVo54vlz4CdmdjthhNub0dGOiIgs0JyJx92vNrPvAq8nJJ7z3P3pdndMRES602w/IH3ttND2+HqUmR3l7j9uX7dERKRbzXbEswW4H0jukpn+IakTHjEgIiIyL7Mlnk8A7wEmCM+0+Tt3H+1Ir0REpGvNOKrN3S929zcBHyXc7fk2M7vOzF7dsd6JiEjXaeXu1I8S7uL8LcLTPV/a7k6JiEj3mm1wwYsId4LeTHii5zXAn7r7ZIf6JiIiXWi2azxbgXsIRzt7gWOA/xgeoQPu/tm2905ERLrObInnj6k/GnqwA30REZEVYLbHIvy3DvZDRERWiFbv1SYiIrIolHhERKSjZkw8ZnZufH1j57ojIiLdbrYjng/G17/qREdERGRlmG1U24Nm9hiw1szuScUNcHc/oa09ExGRrjTbqLb3mdkRwC3AaZ3rkoiIdLNZn8cTn7vzKjMrUL9VzkPuXmp7z0REpCvN+SA4M/s3wJeBxwin2Y42s7Pc/R/b3DcREelCrQyn/izwDnf/N+7+ZuCdwMUH0qiZDZvZ9Wb2UzN70MzeYGaHmNmtZvZwfF0T65qZfc7MtprZPekH1JnZWbH+w2Z2Vir+OjO7N67zOYv3+ZmpDRER6ZxWEk/e3R9KZtz9Z0D+ANu9BPgHd3858CrgQeB84DZ33wjcFucBTgU2xukc4FIISQS4EDiJcNfsC1OJ5NJYN1lvU4zP1IaIiHRIK4lni5ldYWZvidMXgbsX2qCZrQLeDFwB4O5Fd99NuAv2VbHaVcDpsbwZ+LIHdwDDZnYk4cjrVnff6e67gFuBTXHZKnf/obs74TRhelvN2hARkQ5pJfH8DuER2B8DzgUeAD58AG2+iPA47b8xs5+Y2eVmNgC8wN2fAoivh8f66wiPZUhsj7HZ4tubxJmljQZmdo6ZbTGzLTt27GhWRUREFqiVB8FNuftn3f3X3P1X45NJpw6gzRzwWuBSd38NMMbsp7ysWbcWEG+Zu1/m7iPuPrJ27dr5rCoiInNYinu1bQe2u/uP4vz1hET0TDxNRnx9NlX/6NT664En54ivbxJnljZERKRDOp544m+DtpnZy2LoFMLpu5uBZGTaWYQH0BHjH4ij204G9sTTZLcA7zCzNXFQwTuAW+KyfWZ2chzN9oFp22rWhoiIdMisv+Mxsyxwkbv/l0Vu93eBr8Qfpj5CuC9cBrjOzM4GHgd+Pdb9BvAuwhNRx2Nd3H2nmf0JcFes98fuvjOWfwf4EtAHfDNOABfN0IaIiHSIhYFfs1Qw+w5wis9VsUuNjIz4li1blrobIiLLipnd7e4jzZbNeecC4CfATWb2VcJAAADc/cZF6p+IiKwgrSSeQ4DngbemYg4o8YiIyLzNmXjc/YNz1REREWnVnKPazOylZnabmd0X508wsz9sf9dERKQbtTKc+ovABUAJwN3vAc5sZ6dERKR7tZJ4+t39zmmxcjs6IyIi3a+VxPOcmb2YeNsZMzsDeKqtvRIRka7Vyqi2jwCXAS83syeAR4H3t7VXIiLStVoZ1fYI8LZ4B+mMu+9rf7dERKRbtTKq7VAz+xzwfeC7ZnaJmR3a/q6JiEg3auUazzWE5+e8Bzgjlq9tZ6dERKR7tXTnAnf/k9T8p81MT+4UEZEFaeWI53YzO9PMMnF6L/D1dndMRES604xHPGa2j/oTPX8P+D9xUQYYBS5se+9ERKTrzJh43H2okx0REZGVoZVrPJjZCcCGdH09FkFERBZizsRjZlcCJwD3A9UY1mMRRERkQVo54jnZ3Y9re09ERGRFaGVU2w/NTIlHREQWRStHPFcRks/TwBRhlJu7+wlt7ZmIiHSlVhLPlcBvAvdSv8YjIiKyIK0knsfd/ea290RERFaEVhLPT83s/wJ/TzjVBmg4tYiILEwriaePkHDekYppOLWIiCxIK8/j+WA7GjazLLAFeMLd321mxxLuhH0I8GPgN929aGY9wJeB1wHPA7/h7o/FbVwAnA1UgI+5+y0xvgm4BMgCl7v7RTHetI127J+IiDTXyvN4/sbMrpw+LULb5wIPpub/ArjY3TcCuwgJhfi6y91fAlwc6xGHeJ8JHA9sAr5gZtmY0D4PnAocB7wvNRx8pjZERKRDWvkdz9cId6P+OnAbsIpwk9AFM7P1wK8Al8d5A94KXB+rXAUkj17YHOeJy0+J9TcD17j7lLs/CmwFTozTVnd/JB7NXANsnqMNERHpkFZOtd2Qnjezq4FvH2C7fwn8PpDciPRQYLe7l+P8dmBdLK8DtsW+lM1sT6y/Drgjtc30OtumxU+ao40GZnYOcA7AMcccs4DdExGRmbRyxDPdRmDBn8Zm9m7gWXe/Ox1uUtXnWLZY8f2D7pe5+4i7j6xdu7ZZFRERWaBWbhKaPJcn8TRw3gG0+UbgNDN7F9BLOHX3l8CwmeXiEcl64MlYfztwNLDdzHLAamBnKp5Ir9Ms/twsbYiISIfMecTj7kPuvio1vXT66bf5cPcL3H29u28gDA74jru/H7gdOCNWOwu4KZZvjvPE5d9xd4/xM82sJ45W2wjcCdwFbDSzY82sENu4Oa4zUxsiItIhsz2BdNbTae7++CL35TzgGjP7NPAT4IoYvwL4WzPbSjjSOTO2f7+ZXQc8AJSBj7h7Jfb9o8AthOHUV7r7/XO0ISIiHWLhQKDJArN72f/aiANrgcPdPdv+7i29kZER37Jly1J3Q0RkWTGzu919pNmy2R59/cppG9lAOGJ4G/Bni9g/ERFZQVr5AelGM/sS8E3gbuA4d/+rdndMRES602zXeF4BfJJwZ4D/DpydXEMRERFZqNmGU/8r4YeYXyfeESD8+D9w94+1t2siItKNZks8v92xXoiIyIox2+CCq2ZaJiIislALuWWOiIjIginxiIhIRynxiIhIR7XyO56XmtltZnZfnD/BzP6w/V0TEZFu1MoRzxeBC4ASgLvfQ7xfmoiIyHy1knj63f3OabFy05oiIiJzaCXxPGdmLyY+k8fMzgCeamuvRESka835IDjgI8BlwMvN7AngUeDftbVXIiLSteZMPO7+CPA2MxsAMu6+r/3dEhGRbtXKo697gPcAG4Bccr82d//jtvZMRES6Uiun2m4C9hAeiTDV3u6IiEi3ayXxrHf3TW3viYiIrAitjGr7ZzN75dzVRERE5jbbg+DuA6qxzgfN7BHCqTYD3N1P6EwXRUSkm8x2qm0d8OpOdURERFaG2RLPo+7+i471REREVoTZEs/hZvZ7My1098+2oT8iItLlZks8WWCQcE1HRERkUcyWeJ5qx49Ezexo4MvAEYTBC5e5+yVmdghwLeGHqo8B73X3XRZ+sXoJ8C5gHPgtd/9x3NZZQPKIhk8nj+s2s9cBXwL6gG8A57q7z9TGYu+jiIjMbLbh1O060ikDn3D3XwJOBj5iZscB5wO3uftG4LY4D3AqsDFO5wCXAsQkciFwEnAicKGZrYnrXBrrJuslv0OaqQ0REemQ2RLPKe1o0N2fSo5Y4n3fHiSMoNsMXBWrXQWcHsubgS97cAcwbGZHAu8EbnX3nfGo5VZgU1y2yt1/6O5OOLpKb6tZGyIi0iEzJh5339nuxs1sA/Aa4EfAC9z9qdj2U8Dhsdo6YFtqte0xNlt8e5M4s7QxvV/nmNkWM9uyY8eOhe6eiIg00cqdC9rCzAaBG4CPu/ve2ao2ifkC4i1z98vcfcTdR9auXTufVUVEZA6t3Ktt0ZlZnpB0vuLuN8bwM2Z2pLs/FU+XPRvj24GjU6uvB56M8bdMi383xtc3qT9bG4tu285x/vnnz5ExI5sxctkM/fks/T1ZBgo5VvflWd2XZ1VfnmxGAwdFZOXoeOKJo9SuAB6c9lugm4GzgIvi602p+EfN7BrCQII9MXHcAvxZakDBO4AL3H2nme0zs5MJp/A+APzVHG0sun/dvpvzbri3pbqr+/Ks6c+zZqDAIf2F8DpQYE1/gUMG8vE1xA8dKLCqN09GyUpElikL19872KDZm4DvA/cShlMD/AEhSVwHHAM8Dvx6TCIG/C/CyLRx4IPuviVu67fjugB/6u5/E+Mj1IdTfxP43Tic+tBmbczW35GREd+yZcu893OiWGHXeJFK1am6U6o4E8UKo1NlRqfK7J0osWeixO6JEnvGi+wcL7FrrMiu8SK7x0s8PzbFZKnadNsZgzVJguovsGYgn0pUBYZjwhruj7H+AkO9OSUrEekYM7vb3UeaLut04lluFpp4FsNEscLO8SI7R0NC2jVe5PlY3jkWpl3jRXaNldg5XmT3eJFSpfn7mc0Ya/qTZFR/XdNfaIgNx1gyX8gt2WVAEVnGZks8S3KNR1rTV8iyrtDHuuG+luq7O/umyvHIKRxB7UwdRSXJaddYiW07x7lne6hXLDc/sgLoL2RZ018IpwMH8gz3FVjdH04NJuXhvnrSGu7Ls7o/T08uu1j/DCLSZZR4uoiZsao3z6rePC88tLV13J2JUoVd4yV2xwQVjq7CKcAQD8t2jRd5es/eMD9RolKd+Wi5L59ldV+e4f4wgGJ1X0xKMba6L8/qmNDS06reHLmsjrJEupkSzwpnZvQXcvQXci0fWUFIWKNT5ZiUkutVIVHtnagnsT3xWta2nePcNxHqTpQqs257sCfHqt5cLWElo/8aElRfjlW99WWrekOsL58lXBYUkYOVEo8siJkx1JtnqDfP0YfMb92pcoW9E2X2TDQmpz0TpRhPzU+WeHzneG1+vDh70splLCaiXENCCq95hnpCfKg3xIZ6U/N9eQYLGoQh0m5KPNJxPbksa4eyrB3qmfe6pUqVvRMl9k6GBLVvsp6k9k3WRwvumyyzNy57eu9kXKc040jBhBkMFnIM9eZiYq0npqHeHIM9SdKqLx/sqZeTeZ0uFJmZEo8sK/lshkMHezh0cP5JC6BYrrJvMiSufZP1ZJUkqmbxZ/dN8vMdZfbFZTONHEzry2cZjIloKCamwZ4cg7VElWsyH+uk4hpVKN1Iiecgsbe4l+9t+x6P7nmUx/c9zpOjTzJRnmCqMkWxUiRjGfKZPLlMjt5cL325Pvpz/QzkBxjIDzCYH2SgMMBQfojBwiBDhSFWFVYxmA/locIQg/lBspmVPdqskDuwxOXuTJWrtSQUXsuMTiVJq8xoatlosT6/Y99U+A3XZImxqTKzjM1o6O9QT46BOIVyloEkQSXx3nqdwXh3jKROf0+WwR5d/5KDhxLPErvr6bu44eEb+PYvvs1UZYqc5Thy8EjWD67niIEjyGfyFLIFql6lXC1TqpaYLE8yUZ7g+cnn2bZvG2OlMUZLo0yUJ+ZsbyA/UEtEQ/mYnGKiGswPNswnSWywMMhQPqzTk+1Z0R9eZkZvPktvfmGnChPuznixEpNWmPbFhJTExqbK7JsKiWtsqszoVIXRqRLPjRb5xfPj7It15rrulcgYtYSUJK96gsrSnySqQja+xnqFXO1WTwM92RjPMVDI6pSiLIgSzxK69F8u5Qv/+gWG8kOc/pLTOf0lp/PyQ15OLrOwt6VcLTNWGmNvcS+jxVFGS6PsndrLvtI+9hX3n0ZLozwz/gwP7364Nl/12a+B5DK5ekLK1xNWOnkl8w2vqXJ/vp+MrewPLDOrHaEcqErVGSuGJJQkqLGYzEYny4wXG2PjxTJjU5Va+YndEzEWls91HSytkMswUMjWklRfIdcw3x/L/Q2xZL5e7iuExNYX43kltK6mxLME3J3P/8vn+et7/prTXnwaf3TyH9Gb6z3g7eYyOVb3rGZ1z+oF92u8PN6QmGrl4mgtgSXl0eIoY6Uxto1uC4kuJjtv4Wbg6VOESVKqnTKcdvpwIDfAYGGQ/lx/qJcbCPH8AIVMYUUfgUG4K0Xy+63FUKl6LTmNFcuMxyQ1UYqxqTJjxQrj8TU56hovhvmJYpknd0+EdWO98VKF+dwkpZDN1JJQ8tqfz+0fK4RTiEki683HWCFDXz43LR7W681lNXJxiSnxdJi787mffI7L772cX9v4a1z4hgsPmm//Zlb70D9i4IgFbaPqVSbKE7UENVqKU0xSzeaT16fHnmasNFabWklgOcvRn++vHUkl/R/IDzRcA2tYlhugL99XKyfL+nP95LOL8+G9nGUz9aHyi8XdmSxVGU+SUbGe1MaLZSZKFcaScrFSS2DjxQrjpQoTMbHtGi/y5O6w/kQpxOZzhJbozWdqSas3H5NcPkdvIUtfPkNfPiSpvnySxELyCrFsbX6/WKqujtpmpsTTYd9/4vtcfu/lnPHSM/ijk//ooEk6iyVjmdoHPAML306SwJIENVYaY6w81jA/Xh4PSWtanX3FfbUkNl4aZ6w8NucpxEQ+kw+JKCak/nw//bk45fd/7cv17T+fa4yv9OtiEL7U9MUjjhZvqtGyajXcfWOilqDqSWmiVq7UyunXUDckr4lihT0TJZ7ZU19nqhQS32x36ZhJLmOp5JShN1c/4uotZOnNZcKymPx681l6kkRYS2xhvd70fCz35DINdZbTUZwST4dd/dOrWdu3lj846Q+6LukspoYEdoDcncnKJGOlMSZKE4yVxxqS0nhpPJRj4pooTTBeDrHkdefkzrB+eaI2zWdf+nJ9tSmdmPaL5/saYnNNSmqQySze9bJm3J1ipcpkqcpkTFiT5ZjAauVqLflNlSpMxsQ1WQrxyTgl5b0TJZ6txevrznbfxLkUshl68vVk1JNrTFxJourJZeiZNj89kfXkwvyxawd48drBRfzXDJR4Omjb3m380xP/xIdf9WHyGZ3S6RQzq31Q0/pdgWZVqVaYrEzWEtZEuZ6sZisn88nR3LPjzzYks/kkNADDasPr+3J99GZDOYlNXzY9ngz0QSkAABBKSURBVMSmz6fLKz25mVn8IA73H2ynajUM158shYSWHImFcoWpUuOyJHFNliq19aYaloX4eLHMzrEqU+WkXiyXqhQrMye733nLizlv08sXfT+VeDroqz/7KhnL8J6N71nqrsgBymayDGTCEdlaFu/x6FWv1obLzzRNlidryStdd/p6z08+X4snyyYrkwvqVzoJJUdavbnehsTVbFmtHJf3ZnvpyfXUlqXLSZ2FjursBplM/ZRkpyTJLklY6deF/t5tLiv3He6wyfIkN269kbce81ZeMPCCpe6OHKQylqldW2oHd2eqMlVPVJWQpKbKUw3JabI82fA6UQrxqcpULZEl83vG9zBZTi2L8Vavq02Xsxw9uZ5aoipkC7Wk1DDlevaL9eZ6KWRC/UK2QE+2p/aanprFcpncijyyW4pkp8TTId/6xbfYM7WH33jZbyx1V2QFM7PaEUg7uXv4sXNlspbU0olrqjJVW5aUk3i6Tm0qT9XWHy2NUqwUmSxPhtdKve6BMKwhKRWyhf2SVyFboCdTX1aLZwr7rZfE0nXy2fz+86l189k8Oev+BKjE0yHX/vRaNqzawIlHnFiLuTt7n5tkbPcU43uLTI4WKZeqVMpVKmUn/cMHyxhmRiZrWMbIZEI5mc9mDcsamUymFs9kkti0cmobZqFsmfDNp9ZOjFktVu+DGV3/hyEHxsxqH7oUOtOmu1OshoRUqpZqCSudwIqVYi3RlSql/eLNysl8sVpk79Te/dYtVUoUq8UDTnwJw5omr+QuJumElU5cM9bJ5Mln8/NalswnozUXmxJPBzz4/IPc89w9nH/i+ZSLVR7856d44me7eGrrbib2lZa6ewsTk4+lXsk0zptZrJcux7oW/sCS5Zhhqe1CjMdCUrZapdRymq/TMJ9et7FKOtSwrcZYs7pNtjWLhjozrNBSOl+8SvPSju8a8/lR6YExoDdOQS5OjeMmD6BDDlUc9yrVZKJK1VMxnGq1GuZxqp6U969X3443qZ8sd5xqjDuT7kyktgXFOC3MEScM8PH3/9bC/01moMTTATc+fCM92R5ePf4m/u+n7mB05xSrDuvlmOMP5cgXr2bo0F76V/XQN5QnX8iSzWVqRyUQvsm5g1edatXDayVMXnUqlWqtXIt7Uq5SrYJXQj33dNnxauN2p7fj1aT9xuXujf0ima+CMy1OUg4/CfX4f+E11k8uB8Q61A746ttI2iDZBrOsU1+1XqfJZ0pD/WnFdDueanu/yjN8VjVW9xni89PSup37NO9Oi5JhDchiZGl25aThi0wbTx6Ev8WYsGLZY9Jy2C+eJLDkf8esXthdUOaixNNmU5UpvvXwt3nv47/Hd3/wc9YcOcCv/ufjOeolwy1vI31EsbLvLS0i3UCJp81u33Y7Gx99AwNPHsFJm1/Ea95+DFk9Y0VEVjAlnjb72r23cMLTm9h44uGMnLph0bbrxSKVsTF8YoLqxATV8Ql8coLq5BTVyQl8qohPTeHFKapTU/X5UilMxWJ4LZfrr+USlCuhXCmHcqUClfDq1QpUqlCthFNnlUo4j1x1iOetiafUqIZD+Nq5NffG81a1HZnltFD6wk7TcrxOlMnU49NjGYvXkuoxyxhh5ESIWbIsKdfmDQujLpqXs9kZ61g2A/uVDctkIJNt2EatP9lM4/JMNqyfsVQsM608bZ1sFiy1TrJNi/1Il9Ox2nYat1lfJ2674d8r9j1ZJ1me3o5Z+Hcy04AUqVHiaaNnxp6Buw4jQ5aTT3tx0zruTmX3bspPP03p6aepPP885Z27qOzcSWXPHip791Lds4fK6CjVffuojo5SGR+H0sIGJVg+X5vI57FCHsvlsVyuNpHN1ucLhfCBks2ED8Jstv6Bl40feFb/EKx92NdGENQTQuOHz35X/af/wySF+rUVT120qSWyuLzq0xJcuCYVLnCl5uNyr1bqSbFSCUkyJtOGupVKSMIlD8k3Vc9jgqVaqfchFfdqpZ6Uk35UKk3LXq2GtqsLv2XKQS8moYZEv1+ySxL0tCSXzTRP6tOT3PSEN98vCcmXgGZfNtJfBJK4sf+XiVr9aV84Wipba/UyGcAalyVfqpIvF5npX7iarJdelrRd26cM2aEhssOtXxZo1YpMPGa2CbgEyAKXu/tF7Wjnprv/gZc9exIveuMaVh0W7tVSfOwxRn/wT0w99FMmH/oZxa1bqY6P79/H/n6yq1eTXbWK7KpV5I86iuzQIJmBgfrU30+mvx/r6yMTJ+vpJdPbg/X2kunpCYmjtze85vPhPyg5qHlMSuEoMklQDl5tTGyVSmMCi1NS9nh0intjOSa5ZJteiQnS69sMR7rVEGvY9rSkmqyT1E2vM61cPzpO+uv1pF2J22wo05ioU/tWays50q7EOul4pQzFsE7t3zQOlKmVk4Q/Pd6sPNt2al9CuuuLw6Ef+vcc/olPLPp2V1ziMbMs8Hng7cB24C4zu9ndH1jMdtydbd+eZG22zOsP38ez//Oz7Lv9OxS3/hyA7PAwPS97Gavf8x7y644ic/jhZA47FFatgsFBqrkslVKJSrlMpVSiWi5TqpSplstUymWqleS1QqVcxiuTVHaP4dUK1UqcqhWq8Y84zIcPiWr8A/ZkWGe1/sfp6bJX4wi2ODQzWZ6Mcms46iDWoWFZ7dEGnoxO82kjzNJHNKl/v4aZNo3SmnFI8/Tx1elVrOmC6UOla9uYMT593WT4d7oxZqmfVGiyXnreUnsz/XRXwzbrR6Lp7de30+zfJdV+anvWbLtxcIwlo2MseQx3Ptl4k/1I9WWGf9OGI+iGIfTN+5D+92z8J0vHm/270HRf9+tDIiZBzJJhnWF5Eif2JRnOWasTt1OLp09Rh3mD2t9X0qJXHav9bXncT6/93aWGkNbLEJJ/uuwediPejbv3+F+iHVZc4gFOBLa6+yMAZnYNsBlY1MTz9c/8T9Y8fA9WuoUbvjAa/gMc7MFOfiWez1N1p1Ieo/LAnVTuKS9m0/uxTIZMPEWWyWTJxFMU6akWs3C6oF62WA6H4Jnk0Dw5XI9/8EndUCeb+r1Nkz/k2m92mn3YzvBHv+hjTpsns/2GTM+Q9Brqpcqemp++rYZThqk16mcV09upJ2qfVr9WNWmnPo48XaUhnh427ulTlkm9pFJtmzPvR30b+9ed9YvEHH3Z78vGAfcltV5qbL2nP3gbup7+otS4nfS/UWMXpvVrhn1erl6//iiObcN2V2LiWQdsS81vB05KVzCzc4BzAI455pgFNZIZXk0u08vhR/TQs+YlZFavJtfTExNAlmw+TzaXi695csl8oVCfz+fJ5HIhngvLM9kwn8nlyGSzZHN5MtlsbcrmcjHR5MjE8+K6qCty8GhMlI1fGDyVwOv1maH+LEl7v+K0LzBNvpzU+paSK7TnthMrMfE0+xRu/A7jfhlwGcDIyMiCvrK860Mfgg99aCGrikgXq42cZPGP45eLlXileTtwdGp+PfDkEvVFRGTFWYmJ5y5go5kda2YF4Ezg5iXuk4jIirHiTrW5e9nMPgrcQhhOfaW737/E3RIRWTFWXOIBcPdvAN9Y6n6IiKxEK/FUm4iILCElHhER6SglHhER6SglHhER6ShrdgsKqTOzHcAvFrj6YcBzi9idpaR9OTh1y750y36A9iXxQndf22yBEk8bmdkWdx9Z6n4sBu3Lwalb9qVb9gO0L63QqTYREekoJR4REekoJZ72umypO7CItC8Hp27Zl27ZD9C+zEnXeEREpKN0xCMiIh2lxCMiIh2lxNMmZrbJzB4ys61mdv5S92c+zOxoM7vdzB40s/vN7NwYP8TMbjWzh+PrmqXuayvMLGtmPzGzr8X5Y83sR3E/ro2PxzjomdmwmV1vZj+N780blvF78p/if1v3mdnVZta7XN4XM7vSzJ41s/tSsabvgwWfi58D95jZa5eu541m2I//Ef/7usfM/s7MhlPLLoj78ZCZvfNA2lbiaQMzywKfB04FjgPeZ2bHLW2v5qUMfMLdfwk4GfhI7P/5wG3uvhG4Lc4vB+cCD6bm/wK4OO7HLuDsJenV/F0C/IO7vxx4FWGflt17YmbrgI8BI+7+CsLjSc5k+bwvXwI2TYvN9D6cCmyM0znApR3qYyu+xP77cSvwCnc/AfgZcAFA/Ps/Ezg+rvOF+Dm3IEo87XEisNXdH3H3InANsHmJ+9Qyd3/K3X8cy/sIH3DrCPtwVax2FXD60vSwdWa2HvgV4PI4b8BbgetjleWyH6uANwNXALh70d13swzfkygH9JlZDugHnmKZvC/u/o/Azmnhmd6HzcCXPbgDGDazIzvT09k12w93/5a7l+PsHYQnNEPYj2vcfcrdHwW2Ej7nFkSJpz3WAdtS89tjbNkxsw3Aa4AfAS9w96cgJCfg8KXrWcv+Evh9oBrnDwV2p/64lst78yJgB/A38bTh5WY2wDJ8T9z9CeAzwOOEhLMHuJvl+b4kZnoflvNnwW8D34zlRd0PJZ72sCaxZTdu3cwGgRuAj7v73qXuz3yZ2buBZ9397nS4SdXl8N7kgNcCl7r7a4AxlsFptWbi9Y/NwLHAUcAA4ZTUdMvhfZnLsvzvzcw+STjl/pUk1KTagvdDiac9tgNHp+bXA08uUV8WxMzyhKTzFXe/MYafSU4TxNdnl6p/LXojcJqZPUY43flWwhHQcDzFA8vnvdkObHf3H8X56wmJaLm9JwBvAx519x3uXgJuBH6Z5fm+JGZ6H5bdZ4GZnQW8G3i/13/ouaj7ocTTHncBG+MonQLhotzNS9ynlsXrIFcAD7r7Z1OLbgbOiuWzgJs63bf5cPcL3H29u28gvAffcff3A7cDZ8RqB/1+ALj708A2M3tZDJ0CPMAye0+ix4GTzaw//reW7Muye19SZnofbgY+EEe3nQzsSU7JHYzMbBNwHnCau4+nFt0MnGlmPWZ2LGGwxJ0LbsjdNbVhAt5FGBXyc+CTS92fefb9TYTD6HuAf4nTuwjXR24DHo6vhyx1X+exT28BvhbLL4p/NFuBrwI9S92/Fvfh1cCW+L78P2DNcn1PgE8BPwXuA/4W6Fku7wtwNeHaVIlwJHD2TO8D4RTV5+PnwL2EkXxLvg+z7MdWwrWc5O/+f6fqfzLux0PAqQfStm6ZIyIiHaVTbSIi0lFKPCIi0lFKPCIi0lFKPCIi0lFKPCIi0lG5uauISCeYWTIkF+AIoEK4TQ7AuLv/8pJ0TGSRaTi1yEHIzP4bMOrun1nqvogsNp1qE1kGzGw0vr7FzL5nZteZ2c/M7CIze7+Z3Wlm95rZi2O9tWZ2g5ndFac3Lu0eiNQp8YgsP68iPGPolcBvAi919xMJj3743VjnEsKzbV4PvCcuEzko6BqPyPJzl8f7fZnZz4Fvxfi9wL+N5bcBx4VboQGwysyGPDxfSWRJKfGILD9TqXI1NV+l/jedAd7g7hOd7JhIK3SqTaQ7fQv4aDJjZq9ewr6INFDiEelOHwNGzOweM3sA+PBSd0gkoeHUIiLSUTriERGRjlLiERGRjlLiERGRjlLiERGRjlLiERGRjlLiERGRjlLiERGRjvr/aRsPEHnbo+cAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ret.plot(y=[\n",
    "    \"EGFR\", \"EGF(R)\", \"EGFR(CR1^_)\", \"EGFR(Y1068=P^_0)\", \"Grb2(SH2, SH3^1).Sos1(PxxP^1)\",\n",
    "    \"EGFR(Y1068^1).Grb2(SH2^1, SH3^2).Sos1(PxxP^2)\"])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}