Computational tools for strain optimization by adding reactions

This paper introduces a new plug-in for the OptFlux Metabolic Engineering platform, aimed at finding suitable sets of reactions to add to the genomes of microbes (wild type strain), as well as finding complementary sets of deletions, so that the mutant becomes able to overproduce compounds with industrial interest, while preserving their viability. The optimization methods used are Evolutionary Algorithms and Simulated Annealing. The usefulness of this plug-in is demonstrated by a case study, regarding the production of vanillin by the bacterium E. coli

In silico strain optimization by adding reactions to metabolic models

Nowadays, the concerns about the environment and the needs to increase the productivity at low costs, demand for the search of new ways to produce compounds with industrial interest. Based on the increasing knowledge of biological processes, through genome sequencing projects, and high-throughput experimental techniques as well as the available computational tools, the use of microorganisms has been considered as an approach to produce desirable compounds. However, this usually requires to manipulate these organisms by genetic engineering and/ or changing the enviromental conditions to make the production of these compounds possible. In many cases, it is necessary to enrich the genetic material of those microbes with hereologous pathways from other species and consequently adding the potential to produce novel compounds. This paper introduces a new plug-in for the OptFlux Metabolic Engineering platform, aimed at finding suitable sets of reactions to add to the genomes of selected microbes (wild type strain), as well as finding complementary sets of deletions, so that the mutant becomes able to overproduce compounds with industrial interest, while preserving their viability. The necessity of adding reactions to the metabolic model arises from existing gaps in the original model or motivated by the productions of new compounds by the organism. The optimization methods used are metaheuristics such as Evolutionary Algorithms and Simulated Annealing. The usefulness of this plug-in is demonstrated by a case study, regarding the production of vanillin by the bacterium E. coli.This work is supported by project PTDC/EIA-EIA/115176/2009, funded by Portuguese FCT and Programa COMPETE

Visualization of metabolic interaction networks in microbial communities using VisANT 5.0

The complexity of metabolic networks in microbial communities poses an unresolved visualization and interpretation challenge. We address this challenge in the newly expanded version of a software tool for the analysis of biological networks, VisANT 5.0. We focus in particular on facilitating the visual exploration of metabolic interaction between microbes in a community, e.g. as predicted by COMETS (Computation of Microbial Ecosystems in Time and Space), a dynamic stoichiometric modeling framework. Using VisANT's unique metagraph implementation, we show how one can use VisANT 5.0 to explore different time-dependent ecosystem-level metabolic networks. In particular, we analyze the metabolic interaction network between two bacteria previously shown to display an obligate cross-feeding interdependency. In addition, we illustrate how a putative minimal gut microbiome community could be represented in our framework, making it possible to highlight interactions across multiple coexisting species. We envisage that the "symbiotic layout" of VisANT can be employed as a general tool for the analysis of metabolism in complex microbial communities as well as heterogeneous human tissues.This work was supported by the National Institutes of Health, R01GM103502-05 to CD, ZH and DS. Partial support was also provided by grants from the Office of Science (BER), U.S. Department of Energy (DE-SC0004962), the Joslin Diabetes Center (Pilot & Feasibility grant P30 DK036836), the Army Research Office under MURI award W911NF-12-1-0390, National Institutes of Health (1RC2GM092602-01, R01GM089978 and 5R01DE024468), NSF (1457695), and Defense Advanced Research Projects Agency Biological Technologies Office (BTO), Program: Biological Robustness In Complex Settings (BRICS), Purchase Request No. HR0011515303, Program Code: TRS-0 Issued by DARPA/CMO under Contract No. HR0011-15-C-0091. Funding for open access charge: National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. (R01GM103502-05 - National Institutes of Health; 1RC2GM092602-01 - National Institutes of Health; R01GM089978 - National Institutes of Health; 5R01DE024468 - National Institutes of Health; DE-SC0004962 - Office of Science (BER), U.S. Department of Energy; P30 DK036836 - Joslin Diabetes Center; W911NF-12-1-0390 - Army Research Office under MURI; 1457695 - NSF; HR0011515303 - Defense Advanced Research Projects Agency Biological Technologies Office (BTO), Program: Biological Robustness In Complex Settings (BRICS); HR0011-15-C-0091 - DARPA/CMO; National Institutes of Health)Published versio

Coupling metabolic footprinting and flux balance analysis to predict how single gene knockouts perturb microbial metabolism

Tese de mestrado. Biologia (Bioinformática e Biologia Computacional). Universidade de Lisboa, Faculdade de Ciências, 2012The model organisms Caenorhabditis elegans and E. coli form one of the simplest gut microbe host interaction models. Interventions in the microbe that increase the host longevity including inhibition of folate synthesis have been reported previously. To find novel single gene knockouts with an effect on lifespan, a screen of the Keio collection of E. coli was undertaken, and some of the genes found are directly involved in metabolism. The next step in those specific cases is to understand how these mutations perturb metabolism systematically, so that hypotheses can be generated. For that, I employed dynamic Flux Balance Analysis (dFBA), a constraint-based modeling technique capable of simulating the dynamics of metabolism in a batch culture and making predictions about changes in intracellular flux distribution. Since the specificities of the C. elegans lifespan experiments demand us to culture microbes in conditions differing from most of the published literature on E. coli physiology, novel data must be acquired to characterize and make dFBA simulations as realistic as possible. To do this exchange fluxes were measured using quantitative H NMR Time-Resolved Metabolic Footprinting. Furthermore, I also investigate the combination of TReF and dFBA as a tool in microbial metabolism studies. These approaches were tested by comparing wild type E. coli with one of the knockout strains found, ΔmetL, a knockout of the metL gene which encodes a byfunctional enzyme involved in aspartate and threonine metabolism. I found that the strain exhibits a slower growth rate than the wild type. Model simulation results revealed that reduced homoserine and methionine synthesis, as well as impaired sulfur and folate metabolism are the main effects of this knockout and the reasons for the growth deficiency. These results indicate that there are common mechanisms of the lifespan extension between ΔmetL and inhibition of folate biosynthesis and that the flux balance analysis/metabolic footprinting approach can help us understand the nature of these mechanisms.Os organismos modelo Caenorhabditis elegans e E. coli formam um dos modelos mais simples de interacções entre micróbio do tracto digestivo e hospedeiro. Intervenções no micróbio capazes de aumentar a longevidade do hospedeiro, incluindo inibição de síntese de folatos, foram reportadas previamente. Para encontrar novas delecções génicas do micróbio capazes de aumentar a longevidade do hospedeiro, a colecção Keio de deleções génicas de E. coli foi rastreada. Alguns dos genes encontrados participam em processos metabólicos, e nesses casos, o próximpo passo é perceber como as deleções perturbam o metabolismo sistémicamente, para gerar hipóteses. Para isso, utilizo dynamic Flux Balance Analysis (dFBA), uma técnica de modelação metabólica capaz de fazer previsões sobre alterações na distribuição intracelular de fluxos. As especificidades das experiências de tempo de vida em C.elegans obrigam-nos a trabalhar em condições diferentes das usadas na maioria da literatura publicada em fisiologia de E. coli, e para dar o máximo realismo às simulações de dFBA novos dados foram adquiridos, utilizando H NMR Time-Resolved Metabolic Footprinting para medir fluxos de troca de metabolitos entre microorganismo e meio de cultura. A combinação de TReF e dFBA como ferramenta de estudo do metabolism microbiano é também investigada. Estas abordagens foram testadas ao comparar E. coli wild-type com uma das estirpes encontradas no rastreio, ΔmetL, knockout do gene metL, que codifica um enzima bifunctional participante no metabolismo de aspartato e treonina, e que exibe uma taxa de crescimento reduzida comparativamente ao wild-type. Os resultados das simulações revelaram que os principais efeitos da deleção deste gene, e as razões para a menor taxa de crescimento observada, são a produção reduzida de homoserina e metionina e os efeitos que provoca no metabolismo de folatos e enxofre. Estes resultados indicam que há mecanismos comuns na extensão da longevidade causada por esta deleção e inibição de síntese de folatos, e que a combinação metabolic footprinting/flux balance analysis pode ajudar-nos a compreender a natureza desses mecanismos

Automation on the generation of genome scale metabolic models

Background: Nowadays, the reconstruction of genome scale metabolic models is a non-automatized and interactive process based on decision taking. This lengthy process usually requires a full year of one person's work in order to satisfactory collect, analyze and validate the list of all metabolic reactions present in a specific organism. In order to write this list, one manually has to go through a huge amount of genomic, metabolomic and physiological information. Currently, there is no optimal algorithm that allows one to automatically go through all this information and generate the models taking into account probabilistic criteria of unicity and completeness that a biologist would consider. Results: This work presents the automation of a methodology for the reconstruction of genome scale metabolic models for any organism. The methodology that follows is the automatized version of the steps implemented manually for the reconstruction of the genome scale metabolic model of a photosynthetic organism, {\it Synechocystis sp. PCC6803}. The steps for the reconstruction are implemented in a computational platform (COPABI) that generates the models from the probabilistic algorithms that have been developed. Conclusions: For validation of the developed algorithm robustness, the metabolic models of several organisms generated by the platform have been studied together with published models that have been manually curated. Network properties of the models like connectivity and average shortest mean path of the different models have been compared and analyzed.Comment: 24 pages, 2 figures, 2 table

Computational tools for strain optimization by tuning the optimal level of gene expression

In this work, a plug-in for the OptFlux Metabolic Engineering platform is presented, implementing methods that allow the identification of sets of genes to over/under express, relatively to their wild type levels. The optimization methods used are Simulated Annealing and Evolutionary Algorithms, working with a novel representation and operators. This overcomes the limitations of previous approaches based solely on gene knockouts, bringing new avenues for Biotechnology, fostering the discovery of genetic manipulations able to increase the production of certain compounds using a host microbe. The plug-in is made freely available together with appropriate documentation.Support of FCT and Programa COMPETE (ref. PTDC/EIA-EIA/ 115176/2009)

OptForce: An Optimization Procedure for Identifying All Genetic Manipulations Leading to Targeted Overproductions

Computational procedures for predicting metabolic interventions leading to the overproduction of biochemicals in microbial strains are widely in use. However, these methods rely on surrogate biological objectives (e.g., maximize growth rate or minimize metabolic adjustments) and do not make use of flux measurements often available for the wild-type strain. In this work, we introduce the OptForce procedure that identifies all possible engineering interventions by classifying reactions in the metabolic model depending upon whether their flux values must increase, decrease or become equal to zero to meet a pre-specified overproduction target. We hierarchically apply this classification rule for pairs, triples, quadruples, etc. of reactions. This leads to the identification of a sufficient and non-redundant set of fluxes that must change (i.e., MUST set) to meet a pre-specified overproduction target. Starting with this set we subsequently extract a minimal set of fluxes that must actively be forced through genetic manipulations (i.e., FORCE set) to ensure that all fluxes in the network are consistent with the overproduction objective. We demonstrate our OptForce framework for succinate production in Escherichia coli using the most recent in silico E. coli model, iAF1260. The method not only recapitulates existing engineering strategies but also reveals non-intuitive ones that boost succinate production by performing coordinated changes on pathways distant from the last steps of succinate synthesis

OptFlux: an open-source software platform for in silico metabolic engineering

<p>Abstract</p> <p>Background</p> <p>Over the last few years a number of methods have been proposed for the phenotype simulation of microorganisms under different environmental and genetic conditions. These have been used as the basis to support the discovery of successful genetic modifications of the microbial metabolism to address industrial goals. However, the use of these methods has been restricted to bioinformaticians or other expert researchers. The main aim of this work is, therefore, to provide a user-friendly computational tool for Metabolic Engineering applications.</p> <p>Results</p> <p><it>OptFlux </it>is an open-source and modular software aimed at being the reference computational application in the field. It is the first tool to incorporate strain optimization tasks, i.e., the identification of Metabolic Engineering targets, using Evolutionary Algorithms/Simulated Annealing metaheuristics or the previously proposed OptKnock algorithm. It also allows the use of stoichiometric metabolic models for (i) phenotype simulation of both wild-type and mutant organisms, using the methods of Flux Balance Analysis, Minimization of Metabolic Adjustment or Regulatory on/off Minimization of Metabolic flux changes, (ii) Metabolic Flux Analysis, computing the admissible flux space given a set of measured fluxes, and (iii) pathway analysis through the calculation of Elementary Flux Modes.</p> <p><it>OptFlux </it>also contemplates several methods for model simplification and other pre-processing operations aimed at reducing the search space for optimization algorithms.</p> <p>The software supports importing/exporting to several flat file formats and it is compatible with the SBML standard. <it>OptFlux </it>has a visualization module that allows the analysis of the model structure that is compatible with the layout information of <it>Cell Designer</it>, allowing the superimposition of simulation results with the model graph.</p> <p>Conclusions</p> <p>The <it>OptFlux </it>software is freely available, together with documentation and other resources, thus bridging the gap from research in strain optimization algorithms and the final users. It is a valuable platform for researchers in the field that have available a number of useful tools. Its open-source nature invites contributions by all those interested in making their methods available for the community.</p> <p>Given its plug-in based architecture it can be extended with new functionalities. Currently, several plug-ins are being developed, including network topology analysis tools and the integration with Boolean network based regulatory models.</p