{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Predict expression modulation targets\n",
"\n",
"Cameo provides algorithms to search for genes or reactions that can be over or down regulated in order to achieve a given biological objective."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": []
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/javascript": [
"\n",
"(function(global) {\n",
" function now() {\n",
" return new Date();\n",
" }\n",
"\n",
" if (typeof (window._bokeh_onload_callbacks) === \"undefined\") {\n",
" window._bokeh_onload_callbacks = [];\n",
" }\n",
"\n",
" function run_callbacks() {\n",
" window._bokeh_onload_callbacks.forEach(function(callback) { callback() });\n",
" delete window._bokeh_onload_callbacks\n",
" console.info(\"Bokeh: all callbacks have finished\");\n",
" }\n",
"\n",
" function load_libs(js_urls, callback) {\n",
" window._bokeh_onload_callbacks.push(callback);\n",
" if (window._bokeh_is_loading > 0) {\n",
" console.log(\"Bokeh: BokehJS is being loaded, scheduling callback at\", now());\n",
" return null;\n",
" }\n",
" if (js_urls == null || js_urls.length === 0) {\n",
" run_callbacks();\n",
" return null;\n",
" }\n",
" console.log(\"Bokeh: BokehJS not loaded, scheduling load and callback at\", now());\n",
" window._bokeh_is_loading = js_urls.length;\n",
" for (var i = 0; i < js_urls.length; i++) {\n",
" var url = js_urls[i];\n",
" var s = document.createElement('script');\n",
" s.src = url;\n",
" s.async = false;\n",
" s.onreadystatechange = s.onload = function() {\n",
" window._bokeh_is_loading--;\n",
" if (window._bokeh_is_loading === 0) {\n",
" console.log(\"Bokeh: all BokehJS libraries loaded\");\n",
" run_callbacks()\n",
" }\n",
" };\n",
" s.onerror = function() {\n",
" console.warn(\"failed to load library \" + url);\n",
" };\n",
" console.log(\"Bokeh: injecting script tag for BokehJS library: \", url);\n",
" document.getElementsByTagName(\"head\")[0].appendChild(s);\n",
" }\n",
" };\n",
"\n",
" var js_urls = ['https://cdn.pydata.org/bokeh/release/bokeh-0.11.1.min.js', 'https://cdn.pydata.org/bokeh/release/bokeh-widgets-0.11.1.min.js', 'https://cdn.pydata.org/bokeh/release/bokeh-compiler-0.11.1.min.js'];\n",
"\n",
" var inline_js = [\n",
" function(Bokeh) {\n",
" Bokeh.set_log_level(\"info\");\n",
" },\n",
" \n",
" function(Bokeh) {\n",
" Bokeh.$(\"#012ac511-bf21-4670-bcce-4db334096a9e\").text(\"BokehJS successfully loaded\");\n",
" },\n",
" function(Bokeh) {\n",
" console.log(\"Bokeh: injecting CSS: https://cdn.pydata.org/bokeh/release/bokeh-0.11.1.min.css\");\n",
" Bokeh.embed.inject_css(\"https://cdn.pydata.org/bokeh/release/bokeh-0.11.1.min.css\");\n",
" console.log(\"Bokeh: injecting CSS: https://cdn.pydata.org/bokeh/release/bokeh-widgets-0.11.1.min.css\");\n",
" Bokeh.embed.inject_css(\"https://cdn.pydata.org/bokeh/release/bokeh-widgets-0.11.1.min.css\");\n",
" }\n",
" ];\n",
"\n",
" function run_inline_js() {\n",
" for (var i = 0; i < inline_js.length; i++) {\n",
" inline_js[i](window.Bokeh);\n",
" }\n",
" }\n",
"\n",
" if (window._bokeh_is_loading === 0) {\n",
" console.log(\"Bokeh: BokehJS loaded, going straight to plotting\");\n",
" run_inline_js();\n",
" } else {\n",
" load_libs(js_urls, function() {\n",
" console.log(\"Bokeh: BokehJS plotting callback run at\", now());\n",
" run_inline_js();\n",
" });\n",
" }\n",
"}(this));"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from cameo import models"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Load the E. coli core model."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"model = models.bigg.e_coli_core"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Flux Scanning based on Enforced Objective Flux"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"from cameo.strain_design.deterministic.flux_variability_based import FSEOF"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"fseof = FSEOF(model)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"Model: e_coli_core\n",
"Enforced objective: EX_succ_e\n",
"Primary objective: Maximize\n",
"-1.0*BIOMASS_Ecoli_core_w_GAM_reverse_712e5 + 1.0*BIOMASS_Ecoli_core_w_GAM\n",
"
\n",
"Reaction fluxes
\n",
"
\n", " | 1 | \n", "2 | \n", "3 | \n", "4 | \n", "5 | \n", "6 | \n", "7 | \n", "8 | \n", "9 | \n", "10 | \n", "
---|---|---|---|---|---|---|---|---|---|---|
ADK1 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "-0.000000 | \n", "0.000000 | \n", "-0.000000 | \n", "0.000000 | \n", "0.410338 | \n", "1.527283 | \n", "
ENO | \n", "15.238892 | \n", "15.761645 | \n", "16.212407 | \n", "16.577535 | \n", "16.942662 | \n", "17.307790 | \n", "17.704415 | \n", "18.602292 | \n", "19.293292 | \n", "19.628048 | \n", "
FBA | \n", "7.756941 | \n", "8.036499 | \n", "8.242515 | \n", "8.361048 | \n", "8.479582 | \n", "8.598115 | \n", "8.747793 | \n", "9.393102 | \n", "9.830787 | \n", "9.910940 | \n", "
FRD7 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "1.453235 | \n", "3.175302 | \n", "4.897369 | \n", "6.586703 | \n", "7.755141 | \n", "9.127842 | \n", "10.852285 | \n", "
GAPD | \n", "16.431132 | \n", "16.838737 | \n", "17.173618 | \n", "17.421989 | \n", "17.670359 | \n", "17.918730 | \n", "18.198431 | \n", "18.976723 | \n", "19.547784 | \n", "19.761991 | \n", "
ICL | \n", "0.000000 | \n", "0.000000 | \n", "0.386057 | \n", "1.231342 | \n", "2.076627 | \n", "2.921912 | \n", "3.735322 | \n", "4.041459 | \n", "4.143333 | \n", "3.893465 | \n", "
MALS | \n", "0.000000 | \n", "0.000000 | \n", "0.386057 | \n", "1.231342 | \n", "2.076627 | \n", "2.921912 | \n", "3.735322 | \n", "4.041459 | \n", "4.143333 | \n", "3.893465 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
PGM | \n", "-15.238892 | \n", "-15.761645 | \n", "-16.212407 | \n", "-16.577535 | \n", "-16.942662 | \n", "-17.307790 | \n", "-17.704415 | \n", "-18.602292 | \n", "-19.293292 | \n", "-19.628048 | \n", "
PPC | \n", "3.758320 | \n", "5.012330 | \n", "5.878875 | \n", "6.284516 | \n", "6.690157 | \n", "7.095799 | \n", "7.532996 | \n", "8.472367 | \n", "9.615323 | \n", "11.108854 | \n", "
PPS | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.410338 | \n", "1.527283 | \n", "
SUCCt3 | \n", "1.474575 | \n", "2.949150 | \n", "4.423725 | \n", "5.898300 | \n", "7.372875 | \n", "8.847450 | \n", "10.322025 | \n", "11.796600 | \n", "13.271175 | \n", "14.745750 | \n", "
THD2 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "0.119674 | \n", "2.143848 | \n", "2.916807 | \n", "1.535161 | \n", "
TPI | \n", "7.756941 | \n", "8.036499 | \n", "8.242515 | \n", "8.361048 | \n", "8.479582 | \n", "8.598115 | \n", "8.747793 | \n", "9.393102 | \n", "9.830787 | \n", "9.910940 | \n", "
EX_succ_e | \n", "1.474575 | \n", "2.949150 | \n", "4.423725 | \n", "5.898300 | \n", "7.372875 | \n", "8.847450 | \n", "10.322025 | \n", "11.796600 | \n", "13.271175 | \n", "14.745750 | \n", "
18 rows × 10 columns
\n", "\n", " | upper_bound | \n", "lower_bound | \n", "gaps | \n", "normalized_gaps | \n", "KO | \n", "flux_reversal | \n", "suddenly_essential | \n", "
---|---|---|---|---|---|---|---|
ACALD | \n", "-9.10733e-16 | \n", "-9.10733e-16 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "True | \n", "
ACALDt | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "True | \n", "
ACKr | \n", "-8.69634 | \n", "-8.69634 | \n", "0 | \n", "82.443 | \n", "False | \n", "False | \n", "True | \n", "
ACONTa | \n", "0.101495 | \n", "0.101495 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "False | \n", "
ACONTb | \n", "0.101495 | \n", "0.101495 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "False | \n", "
ACt2r | \n", "-8.69634 | \n", "-8.69634 | \n", "0 | \n", "82.443 | \n", "False | \n", "False | \n", "True | \n", "
ADK1 | \n", "0.583116 | \n", "0.583116 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "False | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
SUCDi | \n", "990.126 | \n", "0 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "False | \n", "
SUCOAS | \n", "0 | \n", "-1.88738e-14 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "True | \n", "
TALA | \n", "-0.0168296 | \n", "-0.0168296 | \n", "0 | \n", "-0.141033 | \n", "False | \n", "False | \n", "True | \n", "
THD2 | \n", "1.61298 | \n", "1.61298 | \n", "0 | \n", "0 | \n", "False | \n", "False | \n", "False | \n", "
TKT1 | \n", "-0.0168296 | \n", "-0.0168296 | \n", "0 | \n", "-0.141033 | \n", "False | \n", "False | \n", "True | \n", "
TKT2 | \n", "-0.0507897 | \n", "-0.0507897 | \n", "0 | \n", "-0.425623 | \n", "False | \n", "False | \n", "True | \n", "
TPI | \n", "9.90643 | \n", "9.90643 | \n", "0 | \n", "95.3064 | \n", "False | \n", "False | \n", "False | \n", "
76 rows × 7 columns
\n", "