cameo uses the model data structures defined by cobrapy, our favorite COnstraints-Based Reconstruction and Analysis tool for Python. cameo is thus 100% compatible with cobrapy. For efficiency reasons though cameo implements its own analysis methods that take advantage of a more advanced solver interface.
In [1]:
from cameo import models
model = models.bigg.e_coli_core
In [2]:
from cameo import flux_variability_analysis
In [3]:
fva_result = flux_variability_analysis(model)
fva_result.data_frame
Out[3]:
In [4]:
fva_result.plot(index=fva_result.data_frame.index[:25])
One very useful application of FVA is determining if alternative optimal solution exist.
In [5]:
fva_result2 = flux_variability_analysis(model,fraction_of_optimum=0.5)
fva_result2.data_frame
Out[5]:
In [6]:
fva_result2.plot()
In [7]:
from cameo.visualization import plotting
The phenotypic phase plane is a modeling technique was developed to do a theoretical accessement of what cell can or cannot do in terms of the stoichiometric constraints [Edawards et al. 2001].
The phenotipic phase plane between growth and a product of interest yields the production envelope: a representation between the trade of between production of the desired product and growth.
In [8]:
from cameo import phenotypic_phase_plane
In [9]:
model.reactions.EX_o2_e.lower_bound = -10
result = phenotypic_phase_plane(model,
variables=[model.reactions.BIOMASS_Ecoli_core_w_GAM],
objective=model.reactions.EX_succ_e,
points=10)
In [10]:
result.plot()
The production envelope allows is a quick way to inspect the limitations of the system to design and how the production relates for with growth. In the previous example, succinate prudction is completly decoupled from growth and by decreasing the growth rate it is theoretically possible to produce up to 15 times more succinate.
In [11]:
result.plot(points=[(0.52, 0), (0.23, 12.2)], points_colors=["green", "red"])
The production envelope can show the coupling between growth and production. There is no stoichiometric couple between growth and production for succinate under aerobic conditions, but that is not the case for acetate under anaerobic conditions.
In [12]:
result = phenotypic_phase_plane(model,
variables=[model.reactions.BIOMASS_Ecoli_core_w_GAM],
objective=model.reactions.EX_ac_e)
result.plot()
In [13]:
result.data_frame
Out[13]:
One can immediately see if a design is feasible within the new defined constraints.
In [14]:
result.plot(points=[(0.2, 8), (0.2, 2)], points_colors=["green", "red"])
The computed data can be inspected in the format of a pandas data frame by calling result.data_frame
In [15]:
result.data_frame
Out[15]:
In [16]:
model.reactions.EX_o2_e.lower_bound = 0
result2 = phenotypic_phase_plane(model,
variables=[model.reactions.BIOMASS_Ecoli_core_w_GAM],
objective=model.reactions.EX_ac_e,
points=10)
result2.plot()
In [17]:
from cameo.flux_analysis.analysis import flux_balance_impact_degree
In [18]:
model.reactions.EX_o2_e.lower_bound = -10
In [19]:
%time fbid = flux_balance_impact_degree(model, ["EX_o2_e"])
In [20]:
fbid
Out[20]:
In [ ]: