

In [1]:

    
from cameo import load_model
from cameo.flux_analysis.analysis import phenotypic_phase_plane
from cameo.strain_design.deterministic import DifferentialFVA

E. coli model and succinate production

Load the E. coli core model.



In [2]:

    
model = load_model('iJO1366')

The production envelope looks like this.



In [3]:

    
production_envelope = phenotypic_phase_plane(model, 
                                             variables=[model.reactions.Ec_biomass_iJO1366_core_53p95M],
                                             objective=model.reactions.EX_succ_lp_e_rp_)
production_envelope.plot(height=300)

Set up a model that represents a reference state (in this case a model with a constrained growth rate).



In [4]:

    
reference_model = model.copy()
biomass_rxn = reference_model.reactions.Ec_biomass_iJO1366_core_53p95M
biomass_rxn.lower_bound = 0.3
target = reference_model.reactions.EX_succ_lp_e_rp_
target.lower_bound = 2

Set up the differential flux variability analysis strain design method.



In [5]:

    
diffFVA = DifferentialFVA(design_space_model=model,
                          reference_model=reference_model,
                          objective=target,
                          variables=[biomass_rxn],
                          normalize_ranges_by=biomass_rxn,
                          points=10)

Run differential flux variability analysis (only on the surface of the production envelope)



In [6]:

    
result = diffFVA.run(surface_only=True)
result.plot(height=300)









    




Scanning grid points  0%



In [7]:

    
result.display_on_map("iJO1366.Central metabolism")









    





        
            Ec_biomass_iJO1366_core_53p95M 0.000000
            EX_succ_lp_e_rp_ 17.096429
        
        






    














    







  
  
  

  
  



  
Tip: Hover over an option to see more details about it.
View and build options
Identifiers: ID'sDescriptive names
Scroll to zoom (instead of scroll to pan)Hide secondary metabolitesShow gene reaction rulesHide reaction, gene, and metabolite labelsAllow duplicate reactionsHighlight reactions not in model
Tip: To increase map performance, turn off text boxes (i.e. labels and gene reaction rules).
Reactions
Reaction and gene data not loaded
Value:
Color:
Size:
Styles for reactions with no dataColor:
Size:
Reaction or Gene data
Options: Absolute valueSizeColorText (Show data in label)
Comparison: Fold ChangeLog2(Fold Change)Difference
Method for evaluating AND: MeanMin
Metabolites
Metabolite data not loaded
Value:
Color:
Size:
Styles for metabolites with no dataColor:
Size:
Metabolite data
Options: Absolute valueSizeColorText (Show data in label)
Comparison: Fold ChangeLog2(Fold Change)Difference



In [ ]:

Options:	Absolute valueSizeColorText (Show data in label)
Comparison:	Fold ChangeLog2(Fold Change)Difference

Differential flux variability analysis

E. coli model and succinate production