{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Standard FBA example" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "cobrakbase 0.2.8\n" ] } ], "source": [ "import cobra\n", "from cobra.flux_analysis import (\n", " single_gene_deletion, single_reaction_deletion, double_gene_deletion,\n", " double_reaction_deletion)\n", "from cobra.flux_analysis.parsimonious import pfba\n", "import cobrakbase\n", "import escher\n", "from escher import Builder\n", "from modelseedpy import MSPackageManager,FBAHelper" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "

Objective

1.0 bio1 = 0.1680011696126881

Uptake

\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
MetaboliteReactionFluxC-NumberC-Flux
cpd00009_e0EX_cpd00009_e00.162100.00%
cpd00013_e0EX_cpd00013_e01.81400.00%
cpd00027_e0EX_cpd00027_e056100.00%
cpd00030_e0EX_cpd00030_e00.000116100.00%
cpd00034_e0EX_cpd00034_e05.729E-0500.00%
cpd00048_e0EX_cpd00048_e00.0423100.00%
cpd00058_e0EX_cpd00058_e00.000119100.00%
cpd00063_e0EX_cpd00063_e00.000874400.00%
cpd00099_e0EX_cpd00099_e00.000874400.00%
cpd00149_e0EX_cpd00149_e04.2E-0600.00%
cpd00205_e0EX_cpd00205_e00.0327900.00%
cpd00244_e0EX_cpd00244_e05.426E-0500.00%
cpd00254_e0EX_cpd00254_e00.00145700.00%
cpd10515_e0EX_cpd10515_e00.00138700.00%
cpd10516_e0EX_cpd10516_e00.00131200.00%
cpd11574_e0EX_cpd11574_e01.176E-0600.00%

Secretion

\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
MetaboliteReactionFluxC-NumberC-Flux
cpd01042_c0DM_cpd01042_c0-3.746E-0570.00%
cpd15380_c0DM_cpd15380_c0-0.000112750.00%
cpd00001_e0EX_cpd00001_e0-0.798500.00%
cpd00029_e0EX_cpd00029_e0-3.722232.22%
cpd00036_e0EX_cpd00036_e0-0.0560540.97%
cpd00047_e0EX_cpd00047_e0-8.048134.83%
cpd00067_e0EX_cpd00067_e0-13.4300.00%
cpd00116_e0EX_cpd00116_e0-3.36E-0710.00%
cpd00363_e0EX_cpd00363_e0-3.693231.97%
cpd02701_c0SK_cpd02701_c0-3.36E-07150.00%
" ], "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "kbase_api = cobrakbase.KBaseAPI()\n", "model = kbase_api.get_from_ws(\"E_iML1515.kb\",40576)\n", "media = kbase_api.get_from_ws(\"Btheta_Ecoli_minimal_media\",40576)\n", "pkgmgr = MSPackageManager.get_pkg_mgr(model)\n", "pkgmgr.getpkg(\"KBaseMediaPkg\").build_package(media)\n", "FBAHelper.set_objective_from_target_reaction(model, \"bio1\")\n", "solution=model.optimize()\n", "model.summary()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
growthstatus
ids
(b1123)0.168001optimal
(b1134)0.168001optimal
(b0158)0.168001optimal
(b0775)0.000000optimal
(b1683)0.168001optimal
.........
(b4386)0.168001optimal
(b0867)0.168001optimal
(b1463)0.168001optimal
(b0688)0.168001optimal
(b1475)0.168001optimal
\n", "

1515 rows × 2 columns

\n", "
" ], "text/plain": [ " growth status\n", "ids \n", "(b1123) 0.168001 optimal\n", "(b1134) 0.168001 optimal\n", "(b0158) 0.168001 optimal\n", "(b0775) 0.000000 optimal\n", "(b1683) 0.168001 optimal\n", "... ... ...\n", "(b4386) 0.168001 optimal\n", "(b0867) 0.168001 optimal\n", "(b1463) 0.168001 optimal\n", "(b0688) 0.168001 optimal\n", "(b1475) 0.168001 optimal\n", "\n", "[1515 rows x 2 columns]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "deletion_results = cobra.flsingle_gene_deletion(model)\n", "deletion_results" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "d30ab0c86c3544b890993baee65b5198", "version_major": 2, "version_minor": 0 }, "text/plain": [ "Builder(reaction_data={'34dhpactex_e0': 0.0, 'GUAtex_e0': 0.0, 'rxn01256_c0': 0.05432889023402956, 'rxn00411_c…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "file = open('CentralCarbon.json',)\n", "ccmap = file.read()\n", "builder = Builder(map_json=ccmap)\n", "builder.reaction_scale=[\n", " { \"type\": 'value', \"value\": 0, \"color\": '#dcdcdc', \"size\": 10 },\n", " { \"type\": 'value', \"value\": 0.000001, \"color\": '#9696ff', \"size\": 25 },\n", " { \"type\": 'value', \"value\": 1, \"color\": '#ff0000', \"size\": 35 },\n", " { \"type\": 'value', \"value\": 20, \"color\": '#209123', \"size\": 50 },\n", " ]\n", "builder.model = model\n", "builder.reaction_data = solution.fluxes\n", "builder" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.9.12" } }, "nbformat": 4, "nbformat_minor": 4 }