On the effects of alternative optima in context-specific metabolic model predictions

Recent methodological developments have facilitated the integration of high-throughput data into genome-scale models to obtain context-specific metabolic reconstructions. A unique solution to this data integration problem often may not be guaranteed, leading to a multitude of context-specific predictions equally concordant with the integrated data. Yet, little attention has been paid to the alternative optima resulting from the integration of context-specific data. Here we present computational approaches to analyze alternative optima for different context-specific data integration instances. By using these approaches on metabolic reconstructions for the leaf of Arabidopsis thaliana and the human liver, we show that the analysis of alternative optima is key to adequately evaluating the specificity of the predictions in particular cellular contexts. While we provide several ways to reduce the ambiguity in the context-specific predictions, our findings indicate that the existence of alternative optimal solutions warrant caution in detailed context-specific analyses of metabolism

Experiments and simulations on short chain fatty acid production in a colonic bacterial community

Understanding how production of specific metabolites by gut microbes is modulated by interactions with surrounding species and by environmental nutrient availability is an important open challenge in microbiome research. As part of this endeavor, we explore interactions between F. prausnitzii, a major butyrate producer, and B. thetaiotaomicron, an acetate producer, under three different in vitro media conditions in monoculture and coculture. In silico Genome-scale dynamic flux balance analysis (dFBA) models of metabolism in the system using COMETS (Computation of Microbial Ecosystems in Time and Space) are also tested for explanatory, predictive and inferential power. Experimental findings indicate enhancement of butyrate production in coculture relative to F. prausnitzii monoculture but defy a simple model of monotonic increases in butyrate production as a function of acetate availability in the medium. Simulations recapitulate biomass production curves for monocultures and accurately predict the growth curve of coculture total biomass, using parameters learned from monocultures, suggesting that the model captures some aspects of how the two bacteria interact. However, a comparison of data and simulations for environmental acetate and butyrate changes suggest that the organisms adopt one of many possible metabolic strategies equivalent in terms of growth efficiency. Furthermore, the model seems not to capture subsequent shifts in metabolic activities observed experimentally under low-nutrient regimes. Some discrepancies can be explained by the multiplicity of possible fermentative states for F. prausnitzii. In general, these results provide valuable guidelines for design of future experiments aimed at better determining the mechanisms leading to enhanced butyrate in this ecosystem.https://www.biorxiv.org/content/10.1101/444760v1https://www.biorxiv.org/content/10.1101/444760v1Othe

Dispensability of Escherichia coli's latent pathways

Gene-knockout experiments on single-cell organisms have established that expression of a substantial fraction of genes is not needed for optimal growth. This problem acquired a new dimension with the recent discovery that environmental and genetic perturbations of the bacterium Escherichia coli are followed by the temporary activation of a large number of latent metabolic pathways, which suggests the hypothesis that temporarily activated reactions impact growth and hence facilitate adaptation in the presence of perturbations. Here we test this hypothesis computationally and find, surprisingly, that the availability of latent pathways consistently offers no growth advantage, and tends in fact to inhibit growth after genetic perturbations. This is shown to be true even for latent pathways with a known function in alternate conditions, thus extending the significance of this adverse effect beyond apparently nonessential genes. These findings raise the possibility that latent pathway activation is in fact derivative of another, potentially suboptimal, adaptive response

Impact of stoichiometry representation on simulation of genotype-phenotype relationships in metabolic networks.

<div><p>Genome-scale metabolic networks provide a comprehensive structural framework for modeling genotype-phenotype relationships through flux simulations. The solution space for the metabolic flux state of the cell is typically very large and optimization-based approaches are often necessary for predicting the active metabolic state under specific environmental conditions. The objective function to be used in such optimization algorithms is directly linked with the biological hypothesis underlying the model and therefore it is one of the most relevant parameters for successful modeling. Although linear combination of selected fluxes is widely used for formulating metabolic objective functions, we show that the resulting optimization problem is sensitive towards stoichiometry representation of the metabolic network. This undesirable sensitivity leads to different simulation results when using numerically different but biochemically equivalent stoichiometry representations and thereby makes biological interpretation intrinsically subjective and ambiguous. We hereby propose a new method, Minimization of Metabolites Balance (MiMBl), which decouples the artifacts of stoichiometry representation from the formulation of the desired objective functions, by casting objective functions using metabolite turnovers rather than fluxes. By simulating perturbed metabolic networks, we demonstrate that the use of stoichiometry representation independent algorithms is fundamental for unambiguously linking modeling results with biological interpretation. For example, MiMBl allowed us to expand the scope of metabolic modeling in elucidating the mechanistic basis of several genetic interactions in <em>Saccharomyces cerevisiae</em>.</p> </div

A robust and efficient method for estimating enzyme complex abundance and metabolic flux from expression data

A major theme in constraint-based modeling is unifying experimental data, such as biochemical information about the reactions that can occur in a system or the composition and localization of enzyme complexes, with highthroughput data including expression data, metabolomics, or DNA sequencing. The desired result is to increase predictive capability resulting in improved understanding of metabolism. The approach typically employed when only gene (or protein) intensities are available is the creation of tissue-specific models, which reduces the available reactions in an organism model, and does not provide an objective function for the estimation of fluxes, which is an important limitation in many modeling applications. We develop a method, flux assignment with LAD (least absolute deviation) convex objectives and normalization (FALCON), that employs metabolic network reconstructions along with expression data to estimate fluxes. In order to use such a method, accurate measures of enzyme complex abundance are needed, so we first present a new algorithm that addresses quantification of complex abundance. Our extensions to prior techniques include the capability to work with large models and significantly improved run-time performance even for smaller models, an improved analysis of enzyme complex formation logic, the ability to handle very large enzyme complex rules that may incorporate multiple isoforms, and depending on the model constraints, either maintained or significantly improved correlation with experimentally measured fluxes. FALCON has been implemented in MATLAB and ATS, and can be downloaded from: https://github.com/bbarker/FALCON. ATS is not required to compile the software, as intermediate C source code is available, and binaries are provided for Linux x86-64 systems. FALCON requires use of the COBRA Toolbox, also implemented in MATLAB.Comment: 30 pages, 12 figures, 4 table

Allometric Trajectories and \u201cStress\u201d: A Quantitative Approach

The term \u201cstress\u201d is an important but vague term in plant biology. We show situations in which thinking in terms of \u201cstress\u201d is profitably replaced by quantifying distance from functionally optimal scaling relationships between plant parts. These relationships include, for example, the often-cited one between leaf area and sapwood area, which presumably reflects mutual dependence between source and sink tissues and which scales positively within individuals and across species. These relationships seem to be so basic to plant functioning that they are favored by selection across nearly all plant lineages. Within a species or population, individuals that are far from the common scaling patterns are thus expected to perform negatively. For instance, \u201ctoo little\u201d leaf area (e.g. due to herbivory or disease) per unit of active stem mass would be expected to incur to low carbon income per respiratory cost and thus lead to lower growth. We present a framework that allows quantitative study of phenomena traditionally assigned to \u201cstress,\u201d without need for recourse to this term. Our approach contrasts with traditional approaches for studying \u201cstress,\u201d e.g. revealing that small \u201cstressed\u201d plants likely are in fact well suited to local conditions. We thus offer a quantitative perspective to the study of phenomena often referred to under such terms as \u201cstress,\u201d plasticity, adaptation, and acclimation

Mapping the landscape of metabolic goals of a cell

Genome-scale flux balance models of metabolism provide testable predictions of all metabolic rates in an organism, by assuming that the cell is optimizing a metabolic goal known as the objective function. We introduce an efficient inverse flux balance analysis (invFBA) approach, based on linear programming duality, to characterize the space of possible objective functions compatible with measured fluxes. After testing our algorithm on simulated E. coli data and time-dependent S. oneidensis fluxes inferred from gene expression, we apply our inverse approach to flux measurements in long-term evolved E. coli strains, revealing objective functions that provide insight into metabolic adaptation trajectories.MURI W911NF-12-1-0390 - Army Research Office (US); MURI W911NF-12-1-0390 - Army Research Office (US); 5R01GM089978-02 - National Institutes of Health (US); IIS-1237022 - National Science Foundation (US); DE-SC0012627 - U.S. Department of Energy; HR0011-15-C-0091 - Defense Sciences Office, DARPA; National Institutes of Health; R01GM103502; 5R01DE024468; 1457695 - National Science Foundatio

Dynamic estimations of metabolic fluxes with constraint-based models and possibility theory

Living cells can be modelled by successively imposing known constraints that limit their behaviour, such as mass balances, thermodynamic laws or enzyme capacities. The resulting constraint-based models enclose all the functional states that the modelled cells may exhibit. Then, predictions can be obtained from the models in two main ways: adding experimental data to determine the state of cells at given conditions (MFA) or invoking an assumption of evolved optimal behaviour (FBA). Both MFA and FBA predictions are typically performed at steady state. However, it is easy to take extracellular dynamics into account. This work explores the benefits of using possibility theory to get these dynamic predictions. It will be shown that the possibilistic methods (a) provide rich estimates for time-varying fluxes and metabolite concentrations, (b) account for uncertainty and data scarcity, and (c) give predictions relaxing the optimality assumption of FBA. On the other hand, these methods could serve as basis for monitoring and fault detection systems in industrial bioprocesses.This research has been partially supported by the Spanish Government MINECO (1st and 3rd authors are grateful to grant CICYT DPI2011-28112-C04-01, and A. Sala is grateful to grant DPI2011-27845-C02-01).Llaneras Estrada, F.; Sala, A.; Picó Marco, JA. (2012). Dynamic estimations of metabolic fluxes with constraint-based models and possibility theory. Journal of Process Control. 22(10):1946-1955. https://doi.org/10.1016/j.jprocont.2012.09.00119461955221

Integrated human-virus metabolic stoichiometric modelling predicts host-based antiviral targets against Chikungunya, Dengue and Zika viruses

Current and reoccurring viral epidemic outbreaks such as those caused by the Zika virus illustrate the need for rapid development of antivirals. Such development would be facilitated by computational approaches that can provide experimentally testable predictions for possible antiviral strategies. To this end, we focus here on the fact that viruses are directly dependent on their host metabolism for reproduction. We develop a stoichiometric, genome-scale metabolic model that integrates human macrophage cell metabolism with the biochemical demands arising from virus production and use it to determine the virus impact on host metabolism and vice versa. While this approach applies to any host–virus pair, we first apply it to currently epidemic viruses Chikungunya, Dengue and Zika in this study. We find that each of these viruses causes specific alterations in the host metabolic flux towards fulfilling their biochemical demands as predicted by their genome and capsid structure. Subsequent analysis of this integrated model allows us to predict a set of host reactions, which, when constrained, inhibit virus production. We show that this prediction recovers known targets of existing antiviral drugs, specifically those targeting nucleotide production, while highlighting a set of hitherto unexplored reactions involving both amino acid and nucleotide metabolic pathways, with either broad or virus-specific antiviral potential. Thus, this computational approach allows rapid generation of experimentally testable hypotheses for novel antiviral targets within a host