Evaluating E. coli genome-scale metabolic model accuracy with high-throughput mutant fitness data.

The Escherichia coli genome-scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision-recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene-protein-reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high-throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

Construction and Modeling of a Coculture Microplate for Real-Time Measurement of Microbial Interactions.

The dynamic structures of microbial communities emerge from the complex network of interactions between their constituent microorganisms. Quantitative measurements of these interactions are important for understanding and engineering ecosystem structure. Here, we present the development and application of the BioMe plate, a redesigned microplate device in which pairs of wells are separated by porous membranes. BioMe facilitates the measurement of dynamic microbial interactions and integrates easily with standard laboratory equipment. We first applied BioMe to recapitulate recently characterized, natural symbiotic interactions between bacteria isolated from the Drosophila melanogaster gut microbiome. Specifically, the BioMe plate allowed us to observe the benefit provided by two Lactobacillus strains to an Acetobacter strain. We next explored the use of BioMe to gain quantitative insight into the engineered obligate syntrophic interaction between a pair of Escherichia coli amino acid auxotrophs. We integrated experimental observations with a mechanistic computational model to quantify key parameters associated with this syntrophic interaction, including metabolite secretion and diffusion rates. This model also allowed us to explain the slow growth observed for auxotrophs growing in adjacent wells by demonstrating that, under the relevant range of parameters, local exchange between auxotrophs is essential for efficient growth. The BioMe plate provides a scalable and flexible approach for the study of dynamic microbial interactions. IMPORTANCE Microbial communities participate in many essential processes from biogeochemical cycles to the maintenance of human health. The structure and functions of these communities are dynamic properties that depend on poorly understood interactions among different species. Unraveling these interactions is therefore a crucial step toward understanding natural microbiota and engineering artificial ones. Microbial interactions have been difficult to measure directly, largely due to limitations of existing methods to disentangle the contribution of different organisms in mixed cocultures. To overcome these limitations, we developed the BioMe plate, a custom microplate-based device that enables direct measurement of microbial interactions, by detecting the abundance of segregated populations of microbes that can exchange small molecules through a membrane. We demonstrated the possible application of the BioMe plate for studying both natural and artificial consortia. BioMe is a scalable and accessible platform that can be used to broadly characterize microbial interactions mediated by diffusible molecules.

Enhancing Microbiome Research through Genome-Scale Metabolic Modeling.

Construction and analysis of genome-scale metabolic models (GEMs) is a well-established systems biology approach that can be used to predict metabolic and growth phenotypes. The ability of GEMs to produce mechanistic insight into microbial ecological processes makes them appealing tools that can open a range of exciting opportunities in microbiome research. Here, we briefly outline these opportunities, present current rate-limiting challenges for the trustworthy application of GEMs to microbiome research, and suggest approaches for moving the field forward.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS).

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure.

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis.

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

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.

Coordinated regulation of acid resistance in Escherichia coli.

BACKGROUND: Enteric Escherichia coli survives the highly acidic environment of the stomach through multiple acid resistance (AR) mechanisms. The most effective system, AR2, decarboxylates externally-derived glutamate to remove cytoplasmic protons and excrete GABA. The first described system, AR1, does not require an external amino acid. Its mechanism has not been determined. The regulation of the multiple AR systems and their coordination with broader cellular metabolism has not been fully explored. RESULTS: We utilized a combination of ChIP-Seq and gene expression analysis to experimentally map the regulatory interactions of four TFs: nac, ntrC, ompR, and csiR. Our data identified all previously in vivo confirmed direct interactions and revealed several others previously inferred from gene expression data. Our data demonstrate that nac and csiR directly modulate AR, and leads to a regulatory network model in which all four TFs participate in coordinating acid resistance, glutamate metabolism, and nitrogen metabolism. This model predicts a novel mechanism for AR1 by which the decarboxylation enzymes of AR2 are used with internally derived glutamate. This hypothesis makes several testable predictions that we confirmed experimentally. CONCLUSIONS: Our data suggest that the regulatory network underlying AR is complex and deeply interconnected with the regulation of GABA and glutamate metabolism, nitrogen metabolism. These connections underlie and experimentally validated model of AR1 in which the decarboxylation enzymes of AR2 are used with internally derived glutamate.

Controlling mechanical ventilation in acute respiratory distress syndrome with fuzzy logic.

PURPOSE: The current ventilatory care goal for acute respiratory distress syndrome (ARDS) and the only evidence-based approach for managing ARDS is to ventilate with a tidal volume (VT) of 6 mL/kg predicted body weight (PBW). However, it is not uncommon for some caregivers to feel inclined to deviate from this strategy for one reason or another. To accommodate this inclination in a rationalized manner, we previously developed an algorithm that allows for VT to depart from 6 mL/kg PBW based on physiological criteria. The goal of the present study was to test the feasibility of this algorithm in a small retrospective study. MATERIALS AND METHODS: Current values of peak airway pressure, positive end-expiratory pressure (PEEP), and arterial oxygen saturation are used in a fuzzy logic algorithm to decide how much VT should differ from 6 mL/kg PBW and how much PEEP should change from its current setting. We retrospectively tested the predictions of the algorithm against 26 cases of decision making in 17 patients with ARDS. RESULTS: Differences between algorithm and physician VT decisions were within 2.5 mL/kg PBW, except in 1 of 26 cases, and differences between PEEP decisions were within 2.5 cm H2O, except in 3 of 26 cases. The algorithm was consistently more conservative than physicians in changing VT but was slightly less conservative when changing PEEP. CONCLUSIONS: Within the limits imposed by a small retrospective study, we conclude that our fuzzy logic algorithm makes sensible decisions while at the same time keeping practice close to the current ventilatory care goal.

Elucidating the fuzziness in physician decision making in ARDS.

The current standard of care for patients suffering from acute respiratory distress syndrome (ARDS) is ventilation with a tidal volume of 6 ml/kg predicted body weight (PBW), but variability remains in the tidal volumes that are actually used. This study aims to identify patient scenarios for which there is discordance between physicians in choice of tidal volume and positive end-expiratory pressure (PEEP) in ARDS patients. We developed an algorithm based on fuzzy logic for encapsulating the expertise of individual physicians regarding their use of tidal volume and PEEP in ARDS patients. The algorithm uses three input measurements: (1) peak airway pressure (PAP), (2) PEEP, and (3) arterial oxygen saturation (SaO2). It then generates two output parameters: (1) the deviation of tidal volume from 6 ml/kg PBW, and (2) the change in PEEP from its current value. We captured 6 realizations of intensivist expertise in this algorithm and assessed their degree of concordance using a Monte Carlo simulation. Variability in the tidal volume recommended by the algorithm increased for PAP > 30 cmH2O and PEEP > 5 cmH2O. Tidal volume variability decreased for SaO2 > 90 %. Variability in the recommended change in PEEP increased for PEEP > 5 cmH2O and for SaO2 near 90 %. Intensivists vary in their management of ARDS patients when peak airway pressures and PEEP are high, suggesting that the current goal of 6 ml/kg PBW may need to be revisited under these conditions.