Bactericidal antibiotics alter microbial metabolism as part of their lethality and can damage mitochondria in mammalian cells. In addition, antibiotic susceptibility is sensitive to extracellular ...metabolites, but it remains unknown whether metabolites present at an infection site can affect either treatment efficacy or immune function. Here, we quantify local metabolic changes in the host microenvironment following antibiotic treatment for a peritoneal Escherichia coli infection. Antibiotic treatment elicits microbiome-independent changes in local metabolites, but not those distal to the infection site, by acting directly on host cells. The metabolites induced during treatment, such as AMP, reduce antibiotic efficacy and enhance phagocytic killing. Moreover, antibiotic treatment impairs immune function by inhibiting respiratory activity in immune cells. Collectively, these results highlight the immunomodulatory potential of antibiotics and reveal the local metabolic microenvironment to be an important determinant of infection resolution.
Display omitted
•Antibiotic treatment depletes central metabolism intermediates in the peritoneum•Antibiotic treatment elicits microbiome-independent changes in host metabolites•Metabolites altered by antibiotic treatment during infection inhibit drug efficacy•Antibiotic treatment impairs phagocytic killing by inhibiting respiratory activity
Antibiotic susceptibility is sensitive to metabolites, but how this affects in vivo treatment efficacy remains unexplored. Yang, Bhargava et al. characterize antibiotic-induced changes to the metabolic environment during infection and find that direct actions of antibiotics on host cells induce metabolites that impair drug efficacy and enhance phagocytic activity.
The genome‐scale model (GEM) of metabolism in the bacterium Escherichia coli K‐12 has been in development for over a decade and is now in wide use. GEM‐enabled studies of E. coli have been primarily ...focused on six applications: (1) metabolic engineering, (2) model‐driven discovery, (3) prediction of cellular phenotypes, (4) analysis of biological network properties, (5) studies of evolutionary processes, and (6) models of interspecies interactions. In this review, we provide an overview of these applications along with a critical assessment of their successes and limitations, and a perspective on likely future developments in the field. Taken together, the studies performed over the past decade have established a genome‐scale mechanistic understanding of genotype–phenotype relationships in E. coli metabolism that forms the basis for similar efforts for other microbial species. Future challenges include the expansion of GEMs by integrating additional cellular processes beyond metabolism, the identification of key constraints based on emerging data types, and the development of computational methods able to handle such large‐scale network models with sufficient accuracy.
This review summarizes the applications enabled by genome‐scale models of metabolism for the bacterium E. coli. It provides an overview of the applications along with a critical assessment of their successes and limitations, and a perspective on likely future developments in the field.
Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop ...an integrated “white-box” biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.
Display omitted
•A white-box machine learning approach is developed for antibiotics research•Network modeling is coupled to a biochemical screen to identify pathway mechanisms•Antibiotic-induced adenine limitation increases purine biosynthesis and ATP demand•Increased ATP demand drives central carbon metabolism and oxygen consumption
Causal metabolic pathways underlying antibiotic lethality in bacteria are illuminated by a network model-driven machine learning approach, overcoming limitations of existing “black-box” approaches that cannot reveal causal relationships from large biological datasets.
Generative modeling and representation learning of tandem mass spectrometry data aim to learn an interpretable and instrument-agnostic digital representation of metabolites directly from MS/MS ...spectra. Interpretable and instrument-agnostic digital representations would facilitate comparisons of MS/MS spectra between instrument vendors and enable better and more accurate queries of large MS/MS spectra databases for metabolite identification. In this study, we apply generative modeling and representation learning using variational autoencoders to understand the extent to which tandem mass spectra can be disentangled into their factors of generation (e.g., collision energy, ionization mode, instrument type, etc.) with minimal prior knowledge of the factors. We find that variational autoencoders can disentangle tandem mass spectra data with the proper choice of hyperparameters into meaningful latent representations aligned with known factors of variation. We develop a two-step approach to facilitate the selection of models that are disentangled, which could be applied to other complex and high-dimensional data sets.
Eukaryotic cells compartmentalize biochemical processes in different organelles, often relying on metabolic cycles to shuttle reducing equivalents across intracellular membranes. NADPH serves as the ...electron carrier for the maintenance of redox homeostasis and reductive biosynthesis, with separate cytosolic and mitochondrial pools providing reducing power in each respective location. This cellular organization is critical for numerous functions but complicates analysis of metabolic pathways using available methods. Here we develop an approach to resolve NADP(H)-dependent pathways present within both the cytosol and the mitochondria. By tracing hydrogen in compartmentalized reactions that use NADPH as a cofactor, including the production of 2-hydroxyglutarate by mutant isocitrate dehydrogenase enzymes, we can observe metabolic pathway activity in these distinct cellular compartments. Using this system we determine the direction of serine/glycine interconversion within the mitochondria and cytosol, highlighting the ability of this approach to resolve compartmentalized reactions in intact cells.
Display omitted
•Use of 2H stable isotopes to trace NAD(P)H metabolism in cells•Quantification of pentose phosphate pathway contribution to cytosolic NADPH•Reporter system to distinguish cytosolic and mitochondrial NADPH•Resolve direction of otherwise identical compartmentalized redox reactions
NAD(P)H cofactors are not transported across cell membranes, with metabolic shuttles required for transport of reducing equivalents between the mitochondria and cytosol. Lewis et al. develop an approach to resolve NADPH-dependent pathways in these distinct metabolic compartments and show that directionality of serine/glycine interconversions can be determined within intact cells.
Machine learning has greatly advanced over the past decade, owing to advances in algorithmic innovations, hardware acceleration, and benchmark datasets to train on domains such as computer vision, ...natural-language processing, and more recently the life sciences. In particular, the subfield of machine learning known as deep learning has found applications in genomics, proteomics, and metabolomics. However, a thorough assessment of how the data preprocessing methods required for the analysis of life science data affect the performance of deep learning is lacking. This work contributes to filling that gap by assessing the impact of commonly used as well as newly developed methods employed in data preprocessing workflows for metabolomics that span from raw data to processed data. The results from these analyses are summarized into a set of best practices that can be used by researchers as a starting point for downstream classification and reconstruction tasks using deep learning.
Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not ...just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.
The analytical challenges to acquire accurate isotopic data of intracellular metabolic intermediates for stationary, nonstationary, and dynamic metabolic flux analysis (MFA) are numerous. This work ...presents MID Max, a novel LC–MS/MS workflow, acquisition, and isotopomer deconvolution method for MFA that takes advantage of additional scan types that maximizes the number of mass isotopomer distributions (MIDs) that can be acquired in a given experiment. The analytical method was found to measure the MIDs of 97 metabolites, corresponding to 74 unique metabolite-fragment pairs (32 precursor spectra and 42 product spectra) with accuracy and precision. The compounds measured included metabolic intermediates in central carbohydrate metabolism and cofactors of peripheral metabolism (e.g., ATP). Using only a subset of the acquired MIDs, the method was found to improve the precision of flux estimations and number of resolved exchange fluxes for wild-type E. coli compared to traditional methods and previously published data sets.
Understanding individual variation is fundamental to personalized medicine. Yet interpreting complex phenotype data, such as multi-compartment metabolomic profiles, in the context of genotype data ...for an individual is complicated by interactions within and between cells and remains an unresolved challenge. Here, we constructed multi-omic, data-driven, personalized whole-cell kinetic models of erythrocyte metabolism for 24 healthy individuals based on fasting-state plasma and erythrocyte metabolomics and whole-genome genotyping. We show that personalized kinetic rate constants, rather than metabolite levels, better represent the genotype. Additionally, changes in erythrocyte dynamics between individuals occur on timescales of circulation, suggesting detected differences play a role in physiology. Finally, we use the models to identify individuals at risk for a drug side effect (ribavirin-induced anemia) and how genetic variation (inosine triphosphatase deficiency) may protect against this side effect. This study demonstrates the feasibility of personalized kinetic models, and we anticipate their use will accelerate discoveries in characterizing individual metabolic variation.
Display omitted
•Personalized kinetic models of erythrocytes are built using metabolomic data•Calculated kinetic parameters better represent the genotype than metabolite levels•Individual differences in dynamics occur on physiologically relevant timescales•Model simulations identify a potential mechanism for ribavirin-induced anemia
Bordbar et al. integrate plasma and erythrocyte metabolite measurements of 24 individuals with a metabolic network to construct personalized kinetic models. Estimated kinetic parameters better represent the individual’s genotype than metabolite levels. Subsequently, the models are used to predict drug side effect susceptibility.
Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational ...predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach.