community_model_shuttle_RM_14.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import cobra\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from scipy.optimize import linprog\n",
    "import scipy"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Community modeling shuttle\n",
    "\n",
    "In this notebook we will implement a method to create community models of two or more species specific metabolic models using cobrapy."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Growth:  0.2823653925322476\n",
      "Growth:  2.558694612613393\n"
     ]
    }
   ],
   "source": [
    "model_DP = cobra.io.read_sbml_model(\"models/consistent_DP_SNM.xml\")\n",
    "model_SA = cobra.io.read_sbml_model(\"models/consistent_iYS854_SNM.xml\")\n",
    "print(\"Growth: \", model_DP.slim_optimize())\n",
    "print(\"Growth: \", model_SA.slim_optimize())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "for rec in model_SA.reactions:\n",
    "    rec.lower_bound = max(rec.lower_bound, -1000)\n",
    "    rec.upper_bound = min(rec.upper_bound, 1000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Compound</th>\n",
       "      <th>BiGG</th>\n",
       "      <th>ModelSeed</th>\n",
       "      <th>KEGG</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Alanine</td>\n",
       "      <td>ala__L</td>\n",
       "      <td>cpd00035</td>\n",
       "      <td>C00041</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Arginine</td>\n",
       "      <td>arg__L</td>\n",
       "      <td>cpd00051</td>\n",
       "      <td>C00062</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Cysteine</td>\n",
       "      <td>cys__L</td>\n",
       "      <td>cpd00084</td>\n",
       "      <td>C00097</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Glutamic acid</td>\n",
       "      <td>glu__L</td>\n",
       "      <td>cpd00023</td>\n",
       "      <td>C00025</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Glycine</td>\n",
       "      <td>gly</td>\n",
       "      <td>cpd00033</td>\n",
       "      <td>C00037</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Compound    BiGG ModelSeed    KEGG\n",
       "0        Alanine  ala__L  cpd00035  C00041\n",
       "1       Arginine  arg__L  cpd00051  C00062\n",
       "2       Cysteine  cys__L  cpd00084  C00097\n",
       "3  Glutamic acid  glu__L  cpd00023  C00025\n",
       "4        Glycine     gly  cpd00033  C00037"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "snm3 = pd.read_csv(\"SNM3.csv\", sep =\";\")\n",
    "snm3.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "BIOMASS_DP = \"Growth\" \n",
    "BIOMASS_SA = \"BIOMASS_iYS_wild_type\"\n",
    "models = [model_DP.copy(), model_SA.copy()]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1) Constructing of community model explicitely"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "def create_stoichiometry_matrix(model):\n",
    "    metabolites = model.metabolites \n",
    "    reactions = model.reactions \n",
    "    S = np.zeros((len(metabolites), len(reactions)))\n",
    "    \n",
    "    met_id = dict()\n",
    "    rec_id = dict()\n",
    "    for i,reaction in enumerate(model.reactions):\n",
    "        rec_id[reaction.id] = i\n",
    "        for metabolite, stoich in reaction.metabolites.items():\n",
    "            met_id[metabolite.id] = int(metabolites.index(metabolite))\n",
    "            S[metabolites.index(metabolite), i] = stoich\n",
    "    return S, met_id, rec_id "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "class Model():\n",
    "    def __init__(self, model, biomass_function):\n",
    "        \"\"\" This is a new class of metabolic model, capable of flux balance analysis\n",
    "        Attributes:\n",
    "        models (list): CobraPy models of single organisms which will be used in construction\n",
    "        biomass_reactions (list): List of strings containing the ids for the growth reactions\n",
    "        \"\"\"\n",
    "        self.biomass_function = biomass_function\n",
    "        self.model = model\n",
    "        self.id = model.id\n",
    "        # Compute stoichimetry_matrix\n",
    "        S, met_id, rec_id = create_stoichiometry_matrix(model)\n",
    "        self.num_reactions = S.shape[1]\n",
    "        self.num_metabolites = S.shape[0]\n",
    "        self.stoichiometry_matrix = scipy.sparse.csr_matrix(S)\n",
    "        self.met_id = met_id\n",
    "        self.rec_id = rec_id \n",
    "        # Set objective\n",
    "        idx = self.rec_id[biomass_function]\n",
    "        c = np.zeros(self.num_reactions)\n",
    "        c[idx] = 1\n",
    "        self.objective_c = c\n",
    "        # Set bounds\n",
    "        self._reset_bounds()\n",
    "    \n",
    "    @property\n",
    "    def reactions(self):\n",
    "        return self.model.reactions\n",
    "    @property\n",
    "    def exchanges(self):\n",
    "        return self.model.exchanges\n",
    "    @property\n",
    "    def metabolites(self):\n",
    "        return self.model.metabolites\n",
    "    @property\n",
    "    def medium(self):\n",
    "        return self.model.medium\n",
    "\n",
    "    def set_medium(self, medium):\n",
    "        ex_ids = [ex.id for ex in self.exchanges]\n",
    "        new_med = {}\n",
    "        for key,val in medium.items():\n",
    "            if key in ex_ids:\n",
    "                new_med[key] = val\n",
    "        self.model.medium = new_med\n",
    "        self._reset_bounds()\n",
    "        \n",
    "    def optimize(self, disp=False):\n",
    "        sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
    "        sol[\"fun\"] = -sol[\"fun\"] # As we have to minimize\n",
    "        return sol \n",
    "    \n",
    "    def slim_optimize(self, disp=False):\n",
    "        sol = self.optimize(disp=disp)\n",
    "        return sol[\"fun\"]\n",
    "\n",
    "    def summary(self):\n",
    "        sol = self.optimize()\n",
    "        flux = sol[\"x\"]\n",
    "        ex_ids = [ex.id for ex in self.exchanges]\n",
    "        fluxes = []\n",
    "        for ex in ex_ids:\n",
    "            idx = self.rec_id[ex]\n",
    "            fluxes.append(flux[idx])\n",
    "        summary_df = pd.DataFrame({\"Exchange reaction\": ex_ids, \"Flux\": fluxes})\n",
    "        summary_df.sort_values([\"Flux\"], inplace=True)\n",
    "        return summary_df\n",
    "\n",
    "    def _reset_bounds(self):\n",
    "        self.bounds = []\n",
    "        for rec in self.model.reactions:\n",
    "            self.bounds.append((rec.lower_bound, rec.upper_bound))\n",
    "\n",
    "    def __add__(self, model2):\n",
    "        \"\"\" Adding another model creates a community model \"\"\"\n",
    "        return CommunityModel([self,model2], [1.,1.])\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "class CommunityModel(Model):\n",
    "    def __init__(self, models, weight):\n",
    "        self.models = models\n",
    "        self.id = \"|\".join([model.id for model in models])\n",
    "        self.shared_exchanges = []\n",
    "        self.weight = weight\n",
    "        for model in models:\n",
    "            for rec in model.exchanges:\n",
    "                if rec.id not in self.shared_exchanges:\n",
    "                    self.shared_exchanges.append(rec.id)\n",
    "        # Create community stoichimetry matrix with shuttel reactions!\n",
    "        self._shifts = [len(self.shared_exchanges)]\n",
    "        for i,model in enumerate(models):\n",
    "            self._shifts.append(self._shifts[i] + model.num_reactions)\n",
    "        S_EX = -np.eye(self._shifts[0])\n",
    "        matrices = [S_EX] + [m.stoichiometry_matrix.todense() for m in models]\n",
    "        S = scipy.linalg.block_diag(*matrices)\n",
    "        self.num_reactions = S.shape[1]\n",
    "        self.num_metabolites = S.shape[0]\n",
    "        for i, id in enumerate(self.shared_exchanges):\n",
    "            for j,model in enumerate(models):\n",
    "                if id in model.rec_id:\n",
    "                    S[i,self._shifts[j] + model.rec_id[id]] = 1\n",
    "        self.stoichiometry_matrix = scipy.sparse.csr_matrix(S)\n",
    "        # Cretae objective:\n",
    "        self._weighted_objective(weight)\n",
    "        # Create bounds\n",
    "        self._reset_bounds()\n",
    "    \n",
    "    @property\n",
    "    def reactions(self):\n",
    "        return [model.reactions for model in self.models]\n",
    "    @property\n",
    "    def exchanges(self):\n",
    "        return self.shared_exchanges\n",
    "    @property\n",
    "    def metabolites(self):\n",
    "        return [model.metabolites for model in self.models]\n",
    "    @property\n",
    "    def medium(self):\n",
    "        medium = {}\n",
    "        for model in self.models:\n",
    "            for key,val in model.medium:\n",
    "                medium[key] = val\n",
    "        return medium\n",
    "\n",
    "    def _weighted_objective(self, weight):\n",
    "        self.weight = weight\n",
    "        self.objective_c = np.zeros(self._shifts[0])\n",
    "        for i,model in enumerate(self.models):\n",
    "            self.objective_c = np.append(self.objective_c, weight[i]*model.objective_c)\n",
    "    def _reset_bounds(self):\n",
    "        self.bounds = []\n",
    "        for id in self.shared_exchanges:\n",
    "            min_lower_bound = 0\n",
    "            for model in self.models:\n",
    "                if id in model.rec_id:\n",
    "                    rec = model.reactions.get_by_id(id)\n",
    "                    if rec.lower_bound < min_lower_bound:\n",
    "                        min_lower_bound = rec.lower_bound \n",
    "            self.bounds.append((min_lower_bound, 1000))\n",
    "        for model in self.models:\n",
    "            self.bounds += model.bounds\n",
    "\n",
    "    def set_medium(self, medium):\n",
    "        for model in self.models:\n",
    "            model.set_medium(medium)\n",
    "        self._reset_bounds()\n",
    "\n",
    "    def summary(self):\n",
    "        sol = self.optimize()\n",
    "        flux = sol[\"x\"]\n",
    "        ex_ids = self.shared_exchanges\n",
    "        ex_flux = flux[:len(ex_ids)]\n",
    "        df_ex = pd.DataFrame({\"Exchanges\": ex_ids, \"Flux\": ex_flux})\n",
    "        df_ex.sort_values([\"Flux\"], inplace=True)\n",
    "        for i,model in enumerate(self.models):\n",
    "            shuttel_ids = [ex.id for ex in model.exchanges]\n",
    "            id = str(model.id) + \" Shuttel Flux\"\n",
    "            df_ex[id] = 0.\n",
    "            for sh in shuttel_ids:\n",
    "                idx = model.rec_id[sh]\n",
    "                df_ex[id][ex_ids.index(sh)] = flux[self._shifts[i] +idx]\n",
    "        \n",
    "        return df_ex\n",
    "    \n",
    "    def set_weights(self, weight):\n",
    "        self._weighted_objective(weight)\n",
    "    \n",
    "    def get_model_growths(self):\n",
    "        mask = self.objective_c != 0\n",
    "        sol = self.optimize()\n",
    "        flux = sol[\"x\"]\n",
    "        return flux[mask]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Growth:  0.28236539253224774\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchange reaction</th>\n",
       "      <th>Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>57</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>EX_glc__D_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>EX_arg__L_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>EX_fum_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Exchange reaction  Flux\n",
       "32           EX_o2_e -20.0\n",
       "57       EX_cys__L_e -10.0\n",
       "1        EX_glc__D_e -10.0\n",
       "2           EX_h2o_e -10.0\n",
       "16       EX_arg__L_e -10.0\n",
       "39          EX_fum_e -10.0"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model_DP = Model(model_DP, BIOMASS_DP)\n",
    "print(\"Growth: \",model_DP.slim_optimize())\n",
    "model_DP.summary().head(6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Growth:  2.5586946126118373\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchange reaction</th>\n",
       "      <th>Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>EX_glc__D_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>EX_orn_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td>EX_pyr_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54</th>\n",
       "      <td>EX_glu__L_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>EX_ser__L_e</td>\n",
       "      <td>-10.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Exchange reaction  Flux\n",
       "42           EX_o2_e -20.0\n",
       "51       EX_glc__D_e -10.0\n",
       "43          EX_orn_e -10.0\n",
       "52          EX_pyr_e -10.0\n",
       "54       EX_glu__L_e -10.0\n",
       "56       EX_ser__L_e -10.0"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model_SA = Model(model_SA, BIOMASS_SA)\n",
    "print(\"Growth: \",model_SA.slim_optimize())\n",
    "model_SA.summary().head(6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Weights 1:1:  2.5586946126136296\n",
      "SA growth:  2.5586946126136296\n",
      "DP growth:  -0.0\n",
      "weights 1:10  4.294950548440156\n",
      "SA growth:  1.4712966231176634\n",
      "DP growth:  0.28236539253224924\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n"
     ]
    }
   ],
   "source": [
    "# Community model\n",
    "model = model_SA + model_DP\n",
    "print(\"Weights 1:1: \", model.slim_optimize())\n",
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])\n",
    "model.set_weights([1,10])\n",
    "print(\"weights 1:10 \", model.slim_optimize())\n",
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-8-65ea33eaf871>:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_ex[id][ex_ids.index(sh)] = flux[self._shifts[i] +idx]\n",
      "<ipython-input-8-65ea33eaf871>:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_ex[id][ex_ids.index(sh)] = flux[self._shifts[i] +idx]\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchanges</th>\n",
       "      <th>Flux</th>\n",
       "      <th>iYS854 Shuttel Flux</th>\n",
       "      <th>DP_83VPs_KB5 Shuttel Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-20.000000</td>\n",
       "      <td>-19.874221</td>\n",
       "      <td>-0.125779</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>EX_ser__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54</th>\n",
       "      <td>EX_glu__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-9.828551</td>\n",
       "      <td>-0.171449</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-0.265291</td>\n",
       "      <td>-9.734709</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td>EX_pyr_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>13.604625</td>\n",
       "      <td>23.604625</td>\n",
       "      <td>-10.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>EX_nh4_e</td>\n",
       "      <td>19.684068</td>\n",
       "      <td>19.684068</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>EX_lac__L_e</td>\n",
       "      <td>21.124497</td>\n",
       "      <td>21.124497</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>EX_ser__D_e</td>\n",
       "      <td>25.506917</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>25.506917</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EX_co2_e</td>\n",
       "      <td>31.933467</td>\n",
       "      <td>26.807083</td>\n",
       "      <td>5.126384</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>111 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      Exchanges       Flux  iYS854 Shuttel Flux  DP_83VPs_KB5 Shuttel Flux\n",
       "42      EX_o2_e -20.000000           -19.874221                  -0.125779\n",
       "56  EX_ser__L_e -10.000000           -10.000000                   0.000000\n",
       "54  EX_glu__L_e -10.000000            -9.828551                  -0.171449\n",
       "24  EX_cys__L_e -10.000000            -0.265291                  -9.734709\n",
       "52     EX_pyr_e -10.000000           -10.000000                   0.000000\n",
       "..          ...        ...                  ...                        ...\n",
       "10     EX_h2o_e  13.604625            23.604625                 -10.000000\n",
       "39     EX_nh4_e  19.684068            19.684068                   0.000000\n",
       "19  EX_lac__L_e  21.124497            21.124497                   0.000000\n",
       "98  EX_ser__D_e  25.506917             0.000000                  25.506917\n",
       "23     EX_co2_e  31.933467            26.807083                   5.126384\n",
       "\n",
       "[111 rows x 4 columns]"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = model.summary()\n",
    "df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "df.to_csv(\"exchange_summary_shuttels_RM.csv\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Investigate COMPM\n",
    "\n",
    "The medium should induce competition!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "import json\n",
    "compm_SA = json.loads(open(\"compm_SA.json\").read())\n",
    "compm_DP = json.loads(open(\"compm_DP.json\").read())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "compm = dict()\n",
    "for key, val in compm_SA.items():\n",
    "    compm[key] = val \n",
    "for key, val in compm_DP.items():\n",
    "    compm[key] = val "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "model.set_medium(compm)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SA growth:  2.175506751508334\n",
      "DP growth:  0.07290151281852739\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n"
     ]
    }
   ],
   "source": [
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-7-5ceee4aecab6>:49: OptimizeWarning: Unknown solver options: sparse\n",
      "  sol = linprog(-self.objective_c, A_eq=self.stoichiometry_matrix, b_eq=np.zeros(self.num_metabolites), bounds=self.bounds, method=\"highs\", options={\"disp\":disp})\n",
      "<ipython-input-8-65ea33eaf871>:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_ex[id][ex_ids.index(sh)] = flux[self._shifts[i] +idx]\n",
      "<ipython-input-8-65ea33eaf871>:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_ex[id][ex_ids.index(sh)] = flux[self._shifts[i] +idx]\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchanges</th>\n",
       "      <th>Flux</th>\n",
       "      <th>iYS854 Shuttel Flux</th>\n",
       "      <th>DP_83VPs_KB5 Shuttel Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-20.000000</td>\n",
       "      <td>-19.967526</td>\n",
       "      <td>-0.032474</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>EX_ser__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>EX_orn_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>EX_glc__D_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-7.418185</td>\n",
       "      <td>-2.581815</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-0.361048</td>\n",
       "      <td>-9.638952</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>EX_nh4_e</td>\n",
       "      <td>15.292785</td>\n",
       "      <td>12.282427</td>\n",
       "      <td>3.010359</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>EX_ser__D_e</td>\n",
       "      <td>15.845526</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>15.845526</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>EX_lac__L_e</td>\n",
       "      <td>26.201712</td>\n",
       "      <td>26.201712</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EX_co2_e</td>\n",
       "      <td>29.082032</td>\n",
       "      <td>29.136054</td>\n",
       "      <td>-0.054022</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>31.159853</td>\n",
       "      <td>41.159853</td>\n",
       "      <td>-10.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>111 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      Exchanges       Flux  iYS854 Shuttel Flux  DP_83VPs_KB5 Shuttel Flux\n",
       "42      EX_o2_e -20.000000           -19.967526                  -0.032474\n",
       "56  EX_ser__L_e -10.000000           -10.000000                   0.000000\n",
       "43     EX_orn_e -10.000000           -10.000000                  -0.000000\n",
       "51  EX_glc__D_e -10.000000            -7.418185                  -2.581815\n",
       "24  EX_cys__L_e -10.000000            -0.361048                  -9.638952\n",
       "..          ...        ...                  ...                        ...\n",
       "39     EX_nh4_e  15.292785            12.282427                   3.010359\n",
       "98  EX_ser__D_e  15.845526             0.000000                  15.845526\n",
       "19  EX_lac__L_e  26.201712            26.201712                   0.000000\n",
       "23     EX_co2_e  29.082032            29.136054                  -0.054022\n",
       "10     EX_h2o_e  31.159853            41.159853                 -10.000000\n",
       "\n",
       "[111 rows x 4 columns]"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.summary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "df.to_csv(\"exchange_summary_shuttels_compm_RM.csv\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Investigating COOPM\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "coopm_SA = json.loads(open(\"coopm_SA.json\").read())\n",
    "coopm_DP = json.loads(open(\"coopm_DP.json\").read())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "model.set_medium(coopm_DP)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SA growth:  -0.0\n",
      "DP growth:  0.028236539253225203\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n"
     ]
    }
   ],
   "source": [
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchanges</th>\n",
       "      <th>Flux</th>\n",
       "      <th>iYS854 Shuttel Flux</th>\n",
       "      <th>DP_83VPs_KB5 Shuttel Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>EX_glc__D_e</td>\n",
       "      <td>-1.031864</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-1.031864</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>EX_val__L_e</td>\n",
       "      <td>-0.109437</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.109437</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>EX_thr__L_e</td>\n",
       "      <td>-0.035130</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-0.035130</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>46</th>\n",
       "      <td>EX_pi_e</td>\n",
       "      <td>-0.024136</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.024136</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-0.021831</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.021831</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>108</th>\n",
       "      <td>EX_s_e</td>\n",
       "      <td>0.019177</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.019177</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>EX_12ppd__S_e</td>\n",
       "      <td>0.723979</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.723979</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>0.969719</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.969719</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EX_co2_e</td>\n",
       "      <td>1.302146</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.302146</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>EX_acald_e</td>\n",
       "      <td>1.369230</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.369230</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>111 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         Exchanges      Flux  iYS854 Shuttel Flux  DP_83VPs_KB5 Shuttel Flux\n",
       "51     EX_glc__D_e -1.031864                 -0.0                  -1.031864\n",
       "63     EX_val__L_e -0.109437                 -0.0                  -0.109437\n",
       "89     EX_thr__L_e -0.035130                  0.0                  -0.035130\n",
       "46         EX_pi_e -0.024136                 -0.0                  -0.024136\n",
       "24     EX_cys__L_e -0.021831                 -0.0                  -0.021831\n",
       "..             ...       ...                  ...                        ...\n",
       "108         EX_s_e  0.019177                  0.0                   0.019177\n",
       "83   EX_12ppd__S_e  0.723979                  0.0                   0.723979\n",
       "10        EX_h2o_e  0.969719                  0.0                   0.969719\n",
       "23        EX_co2_e  1.302146                  0.0                   1.302146\n",
       "5       EX_acald_e  1.369230                  0.0                   1.369230\n",
       "\n",
       "[111 rows x 4 columns]"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.summary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "model.set_medium(coopm_SA)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SA growth:  0.25586946126113846\n",
      "DP growth:  -0.0\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n"
     ]
    }
   ],
   "source": [
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchanges</th>\n",
       "      <th>Flux</th>\n",
       "      <th>iYS854 Shuttel Flux</th>\n",
       "      <th>DP_83VPs_KB5 Shuttel Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-15.182832</td>\n",
       "      <td>-15.182832</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td>EX_pyr_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>EX_na1_e</td>\n",
       "      <td>-0.509052</td>\n",
       "      <td>-0.509052</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>EX_k_e</td>\n",
       "      <td>-0.047512</td>\n",
       "      <td>-0.047512</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-0.046136</td>\n",
       "      <td>-0.046136</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>110</th>\n",
       "      <td>EX_26dap__M_e</td>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>EX_glyclt_e</td>\n",
       "      <td>0.000114</td>\n",
       "      <td>0.000114</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>EX_ac_e</td>\n",
       "      <td>1.575224</td>\n",
       "      <td>1.575224</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>15.623422</td>\n",
       "      <td>15.623422</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EX_co2_e</td>\n",
       "      <td>18.033212</td>\n",
       "      <td>18.033212</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>111 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         Exchanges       Flux  iYS854 Shuttel Flux  DP_83VPs_KB5 Shuttel Flux\n",
       "42         EX_o2_e -15.182832           -15.182832                       -0.0\n",
       "52        EX_pyr_e -10.000000           -10.000000                        0.0\n",
       "35        EX_na1_e  -0.509052            -0.509052                       -0.0\n",
       "16          EX_k_e  -0.047512            -0.047512                       -0.0\n",
       "24     EX_cys__L_e  -0.046136            -0.046136                       -0.0\n",
       "..             ...        ...                  ...                        ...\n",
       "110  EX_26dap__M_e  -0.000000             0.000000                        0.0\n",
       "82     EX_glyclt_e   0.000114             0.000114                       -0.0\n",
       "4          EX_ac_e   1.575224             1.575224                        0.0\n",
       "10        EX_h2o_e  15.623422            15.623422                        0.0\n",
       "23        EX_co2_e  18.033212            18.033212                        0.0\n",
       "\n",
       "[111 rows x 4 columns]"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.summary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "coopm = dict()\n",
    "for key, val in compm_SA.items():\n",
    "    coopm[key] = min(coopm.get(key, 11),val)\n",
    "for key, val in compm_DP.items():\n",
    "    coopm[key] = min(coopm.get(key, 11),val) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "model.set_medium(coopm) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SA growth:  1.4982832805635986\n",
      "DP growth:  0.1356402884180246\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n"
     ]
    }
   ],
   "source": [
    "growths = model.get_model_growths()\n",
    "print(\"SA growth: \", growths[0])\n",
    "print(\"DP growth: \", growths[1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Sparse constraint matrix detected; setting 'sparse':True.\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:49: OptimizeWarning: Unknown solver options: sparse\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "C:\\Users\\Manuel Glöckler\\AppData\\Roaming\\Python\\Python37\\site-packages\\ipykernel_launcher.py:83: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Exchanges</th>\n",
       "      <th>Flux</th>\n",
       "      <th>iYS854 Shuttel Flux</th>\n",
       "      <th>DP_83VPs_KB5 Shuttel Flux</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>EX_o2_e</td>\n",
       "      <td>-11.000000</td>\n",
       "      <td>-10.939580</td>\n",
       "      <td>-0.060420</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>EX_ser__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>EX_glc__D_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-5.196285</td>\n",
       "      <td>-4.803715</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EX_cys__L_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-0.248288</td>\n",
       "      <td>-9.751712</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td>EX_pyr_e</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>-10.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>EX_nh4_e</td>\n",
       "      <td>11.639244</td>\n",
       "      <td>11.639244</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>EX_h2o_e</td>\n",
       "      <td>16.842553</td>\n",
       "      <td>26.842553</td>\n",
       "      <td>-10.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>EX_ser__D_e</td>\n",
       "      <td>19.986799</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>19.986799</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EX_co2_e</td>\n",
       "      <td>22.174276</td>\n",
       "      <td>21.831561</td>\n",
       "      <td>0.342715</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>EX_lac__L_e</td>\n",
       "      <td>24.579253</td>\n",
       "      <td>24.579253</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>111 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      Exchanges       Flux  iYS854 Shuttel Flux  DP_83VPs_KB5 Shuttel Flux\n",
       "42      EX_o2_e -11.000000           -10.939580                  -0.060420\n",
       "56  EX_ser__L_e -10.000000           -10.000000                   0.000000\n",
       "51  EX_glc__D_e -10.000000            -5.196285                  -4.803715\n",
       "24  EX_cys__L_e -10.000000            -0.248288                  -9.751712\n",
       "52     EX_pyr_e -10.000000           -10.000000                   0.000000\n",
       "..          ...        ...                  ...                        ...\n",
       "39     EX_nh4_e  11.639244            11.639244                   0.000000\n",
       "10     EX_h2o_e  16.842553            26.842553                 -10.000000\n",
       "98  EX_ser__D_e  19.986799             0.000000                  19.986799\n",
       "23     EX_co2_e  22.174276            21.831561                   0.342715\n",
       "19  EX_lac__L_e  24.579253            24.579253                   0.000000\n",
       "\n",
       "[111 rows x 4 columns]"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.summary()"
   ]
  }
 ],
 "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.8.5"
  },
  "metadata": {
   "interpreter": {
    "hash": "a8f61be024eba58adef938c9aa1e29e02cb3dece83a5348b1a2dafd16a070453"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}