In [1]:
from marsi.cobra.strain_design.post_processing import replace_strain_design_results
In [2]:
from cameo import models
from cameo.strain_design import OptGene
from cameo.util import TimeMachine
from pandas import DataFrame
In [3]:
model = models.bigg.iJO1366
In [4]:
model.reactions.EX_o2_e.lower_bound = 0
We search for reactions, because they can be directly translated into metabolite targets
In [ ]:
optgene = OptGene(model, manipulation_type="reactions", essential_reactions=["ATPM"])
In [ ]:
results = optgene.run(target="EX_lac__D_e",
biomass="BIOMASS_Ec_iJO1366_core_53p95M",
substrate="EX_glc__D_e",
max_knockouts=15,
max_evaluations=15000)
In [ ]:
df = results.data_frame
In [ ]:
df
In [ ]:
replacements = replace_strain_design_results(model, results, results._objective_function,
results._simulation_method, results._simulation_kwargs)
In [ ]:
replacements.style
In [ ]: