Enumeration of minimal stoichiometric precursor sets in metabolic networks

Background: What an organism needs at least from its environment to produce a set of metabolites, e.g. target(s) of interest and/or biomass, has been called a minimal precursor set. Early approaches to enumerate all minimal precursor sets took into account only the topology of the metabolic network (topological precursor sets). Due to cycles and the stoichiometric values of the reactions, it is often not possible to produce the target(s) from a topological precursor set in the sense that there is no feasible flux. Although considering the stoichiometry makes the problem harder, it enables to obtain biologically reasonable precursor sets that we call stoichiometric. Recently a method to enumerate all minimal stoichiometric precursor sets was proposed in the literature. The relationship between topological and stoichiometric precursor sets had however not yet been studied. Results: Such relationship between topological and stoichiometric precursor sets is highlighted. We also present two algorithms that enumerate all minimal stoichiometric precursor sets. The first one is of theoretical interest only and is based on the above mentioned relationship. The second approach solves a series of mixed integer linear programming problems. We compared the computed minimal precursor sets to experimentally obtained growth media of several Escherichia coli strains using genome-scale metabolic networks. Conclusions: The results show that the second approach efficiently enumerates minimal precursor sets taking stoichiometry into account, and allows for broad in silico studies of strains or species interactions that may help to understand e.g. pathotype and niche-specific metabolic capabilities. sasita is written in Java, uses cplex as LP solver and can be downloaded together with all networks and input files used in this paper at http://www.sasita.gforge.inria.fr

Genome-scale gene/reaction essentiality and synthetic lethality analysis

Synthetic lethals are to pairs of non-essential genes whose simultaneous deletion prohibits growth. One can extend the concept of synthetic lethality by considering gene groups of increasing size where only the simultaneous elimination of all genes is lethal, whereas individual gene deletions are not. We developed optimization-based procedures for the exhaustive and targeted enumeration of multi-gene (and by extension multi-reaction) lethals for genome-scale metabolic models. Specifically, these approaches are applied to iAF1260, the latest model of Escherichia coli, leading to the complete identification of all double and triple gene and reaction synthetic lethals as well as the targeted identification of quadruples and some higher-order ones. Graph representations of these synthetic lethals reveal a variety of motifs ranging from hub-like to highly connected subgraphs providing a birds-eye view of the avenues available for redirecting metabolism and uncovering complex patterns of gene utilization and interdependence. The procedure also enables the use of falsely predicted synthetic lethals for metabolic model curation. By analyzing the functional classifications of the genes involved in synthetic lethals, we reveal surprising connections within and across clusters of orthologous group functional classifications

Metabolic network percolation quantifies biosynthetic capabilities across the human oral microbiome

The biosynthetic capabilities of microbes underlie their growth and interactions, playing a prominent role in microbial community structure. For large, diverse microbial communities, prediction of these capabilities is limited by uncertainty about metabolic functions and environmental conditions. To address this challenge, we propose a probabilistic method, inspired by percolation theory, to computationally quantify how robustly a genome-derived metabolic network produces a given set of metabolites under an ensemble of variable environments. We used this method to compile an atlas of predicted biosynthetic capabilities for 97 metabolites across 456 human oral microbes. This atlas captures taxonomically-related trends in biomass composition, and makes it possible to estimate inter-microbial metabolic distances that correlate with microbial co-occurrences. We also found a distinct cluster of fastidious/uncultivated taxa, including several Saccharibacteria (TM7) species, characterized by their abundant metabolic deficiencies. By embracing uncertainty, our approach can be broadly applied to understanding metabolic interactions in complex microbial ecosystems.T32GM008764 - NIGMS NIH HHS; T32 GM008764 - NIGMS NIH HHS; R01 DE024468 - NIDCR NIH HHS; R01 GM121950 - NIGMS NIH HHS; DE-SC0012627 - Biological and Environmental Research; RGP0020/2016 - Human Frontier Science Program; NSFOCE-BSF 1635070 - National Science Foundation; HR0011-15-C-0091 - Defense Advanced Research Projects Agency; R37DE016937 - NIDCR NIH HHS; R37 DE016937 - NIDCR NIH HHS; R01GM121950 - NIGMS NIH HHS; R01DE024468 - NIDCR NIH HHS; 1457695 - National Science FoundationPublished versio

Quantifying biosynthetic network robustness across the human oral microbiome

Metabolic interactions, such as cross-feeding, play a prominent role in microbial communitystructure. For example, they may underlie the ubiquity of uncultivated microorganisms. We investigated this phenomenon in the human oral microbiome, by analyzing microbial metabolic networks derived from sequenced genomes. Specifically, we devised a probabilistic biosynthetic network robustness metric that describes the chance that an organism could produce a given metabolite, and used it to assemble a comprehensive atlas of biosynthetic capabilities for 88 metabolites across 456 human oral microbiome strains. A cluster of organisms characterized by reduced biosynthetic capabilities stood out within this atlas. This cluster included several uncultivated taxa and three recently co-cultured Saccharibacteria (TM7) phylum species. Comparison across strains also allowed us to systematically identify specific putative metabolic interdependences between organisms. Our method, which provides a new way of converting annotated genomes into metabolic predictions, is easily extendible to other microbial communities and metabolic products.https://www.biorxiv.org/content/10.1101/392621v1First author draf

Elementary approaches to microbial growth rate maximisation

This thesis, called Elementary approaches to microbial growth rate maximisation, reports on a theoretical search for principles underlying single cell growth, in particular for microbial species that are selected for fast growth rates. First, the optimally growing cell is characterised in terms of its elementary modes. We prove an extremum principle: a cell that maximises a metabolic rate uses few Elementary Flux Modes (EFMs, the minimal pathways that support steady-state metabolism). The number of active EFMs is bounded by the number of growth-limiting constraints. Later, this extremum principle is extended in a theory that explicitly accounts for self-fabrication. For this, we had to define the elementary modes that underlie balanced self-fabrication: minimal self-supporting sets of expressed enzymes that we call Elementary Growth Modes (EGMs). It turns out that many of the results for EFMs can be extended to their more general self-fabrication analogue. Where the above extremum principles tell us that few elementary modes are used by a rate-maximising cell, it does not tell us how the cell can find them. Therefore, we also search for an elementary adaptation method. It turns out that stochastic phenotype switching with growth rate dependent switching rates provides an adaptation mechanism that is often competitive with more conventional regulatory-circuitry based mechanisms. The derived theory is applied in two ways. First, the extremum principles are used to review the mathematical fundaments of all optimisation-based explanations of overflow metabolism. Second, a computational tool is presented that enumerates Elementary Conversion Modes. These elementary modes can be computed for larger networks than EFMs and EGMs, and still provide an overview of the metabolic capabilities of an organism

Flux analysis in central carbon metabolism in plants: 13C NMR experiments and analysis

Metabolic flux analysis is crucial in metabolic engineering. This research concentrated on improvements in 13C labeling-based flux analysis, a powerful flux quantification method, particularly oriented toward application to plants. Furthermore, systemic 13C flux analyses were performed on two model plant systems: Glycine max (soybean) embryos, and Catharanthus roseus hairy roots.;The concepts \u27bond integrity\u27, \u27bondomer\u27 and the algorithm \u27Boolean function mapping\u27 were introduced, to facilitate efficient flux evaluation from carbon bond labeling experiments, and easier flux identifiability analysis.;13C labeling experiments were performed on developing soybean (Glycine max) embryos and C. roseus hairy roots. A computer program, NMR2Flux, was developed to automatically calculate fluxes from the labeling data. This program accepts a user-defined metabolic network model, and incorporates recent mathematical advances toward accurate and efficient evaluation of fluxes and their standard deviations. Several physiological insights were obtained from the flux results. For instance, in soybean embryos, the reductive pentose phosphate pathway was active in the plastid and negligible in the cytosol. Also, unknown fluxes (such as plastidic fructose-1,6-bisphosphatase) could be identified and quantified. To the best of the author\u27s knowledge, this is the most comprehensive flux analysis of a plant system to date.;Investigations on flux identifiability were carried out for the soybean embryo system. Using these, optimal labeling experiments were designed, that utilize judicious combinations of labeled varieties of two substrates (sucrose and glutamine), to maximize the statistical quality of the evaluated fluxes.;The identity of four intense peaks observed in the 2-D [13C, 1H] spectra of protein isolated from soybean embryos, was investigated. These peaks were identified as levulinic acid and 5-hydroxymethyl furfural, and were degradation products of glycosylating sugars associated with soybean embryo protein. A 2-D NMR study was conducted on them, and it was shown that the metabolic information in the degradation products can be used toward metabolic flux or pathway analysis.;In addition, the elemental make-up and composition of the biomass of C. roseus hairy roots (crucial toward flux analysis) is reported. 89.2% (+/-9.7%) of the biomass was accounted for.*;*This dissertation is a compound document (contains both a paper copy and a CD as part of the dissertation). The CD requires the following system requirements: Adobe Acrobat

A retrosynthetic biology approach to metabolic pathway design for therapeutic production

<p>Abstract</p> <p>Background</p> <p>Synthetic biology is used to develop cell factories for production of chemicals by constructively importing heterologous pathways into industrial microorganisms. In this work we present a retrosynthetic approach to the production of therapeutics with the goal of developing an <it>in situ </it>drug delivery device in host cells. Retrosynthesis, a concept originally proposed for synthetic chemistry, iteratively applies reversed chemical transformations (reversed enzyme-catalyzed reactions in the metabolic space) starting from a target product to reach precursors that are endogenous to the chassis. So far, a wider adoption of retrosynthesis into the manufacturing pipeline has been hindered by the complexity of enumerating all feasible biosynthetic pathways for a given compound.</p> <p>Results</p> <p>In our method, we efficiently address the complexity problem by coding substrates, products and reactions into molecular signatures. Metabolic maps are represented using hypergraphs and the complexity is controlled by varying the specificity of the molecular signature. Furthermore, our method enables candidate pathways to be ranked to determine which ones are best to engineer. The proposed ranking function can integrate data from different sources such as host compatibility for inserted genes, the estimation of steady-state fluxes from the genome-wide reconstruction of the organism's metabolism, or the estimation of metabolite toxicity from experimental assays. We use several machine-learning tools in order to estimate enzyme activity and reaction efficiency at each step of the identified pathways. Examples of production in bacteria and yeast for two antibiotics and for one antitumor agent, as well as for several essential metabolites are outlined.</p> <p>Conclusions</p> <p>We present here a unified framework that integrates diverse techniques involved in the design of heterologous biosynthetic pathways through a retrosynthetic approach in the reaction signature space. Our engineering methodology enables the flexible design of industrial microorganisms for the efficient on-demand production of chemical compounds with therapeutic applications.</p