Systematic applications of metabolomics in metabolic engineering

Precision fermentation with mass spectrometry-based spent media analysis

Optimization and monitoring of bioprocesses requires the measurement of several process parameters and quality attributes. Mass spectrometry (MS)-based techniques such as those coupled to gas chromatography (GCMS) and liquid Chromatography (LCMS) enable the simultaneous measurement of hundreds of metabolites with high sensitivity. When applied to spent media, such metabolome analysis can help determine the sequence of substrate uptake and metabolite secretion, consequently facilitating better design of initial media and feeding strategy. Furthermore, the analysis of metabolite diversity and abundance from spent media will aid the determination of metabolic phases of the culture and the identification of metabolites as surrogate markers for product titer and quality. This review covers the recent advances in metabolomics analysis applied to the development and monitoring of bioprocesses. In this regard, we recommend a stepwise workflow and guidelines that a bioprocesses engineer can adopt to develop and optimize a fermentation process using spent media analysis. Finally, we show examples of how the use of MS can revolutionize the design and monitoring of bioprocesses.

Biomarker enrichment medium: A defined medium for metabolomic analysis of microbial pathogens

Microbes have diverse metabolic capabilities and differences in these phenotypes are critical for differentiating strains, species, and broader taxa of microorganisms. Recent advances in liquid chromatography-mass spectrometry (LC-MS) allow researchers to track the complex combinations of molecules that are taken up by each cell type and to quantify the rates that individual metabolites enter or exit the cells. This metabolomics-based approach allows complex metabolic phenotypes to be captured in a single assay, enables computational models of microbial metabolism to be constructed, and can serve as a diagnostic approach for clinical microbiology. Unfortunately, metabolic phenotypes are directly affected by the molecular composition of the culture medium and many traditional media are subject to molecular-level heterogeneity. Herein, we show that commercially sourced Mueller Hinton (MH) medium, a Clinical and Laboratory Standards Institute (CLSI) approved medium for clinical microbiology, has significant lot-to-lot and supplier-to-supplier variability in the concentrations of individual nutrients. We show that this variability does not affect microbial growth rates but does affect the metabolic phenotypes observed in vitro—including metabolic phenotypes that distinguish six common pathogens. To address this, we used a combination of isotope-labeling, substrate exclusion, and nutritional supplementation experiments using Roswell Park Memorial Institute (RPMI) medium to identify the specific nutrients used by the microbes to produce diagnostic biomarkers, and to formulate a Biomarker Enrichment Medium (BEM) as an alternative to complex undefined media for metabolomics research, clinical diagnostics, antibiotic susceptibility testing, and other applications where the analysis of stable microbial metabolic phenotypes is important.

Chemical-damage MINE: A database of curated and predicted spontaneous metabolic reactions

Spontaneous reactions between metabolites are often neglected in favor of emphasizing enzyme-catalyzed chemistry because spontaneous reaction rates are assumed to be insignificant under physiological conditions. However, synthetic biology and engineering efforts can raise natural metabolites' levels or introduce unnatural ones, so that previously innocuous or nonexistent spontaneous reactions become an issue. Problems arise when spontaneous reaction rates exceed the capacity of a platform organism to dispose of toxic or chemically active reaction products. While various reliable sources list competing or toxic enzymatic pathways' side-reactions, no corresponding compilation of spontaneous side-reactions exists, nor is it possible to predict their occurrence. We addressed this deficiency by creating the Chemical Damage (CD)-MINE resource. First, we used literature data to construct a comprehensive database of metabolite reactions that occur spontaneously in physiological conditions. We then leveraged this data to construct 148 reaction rules describing the known spontaneous chemistry in a substrate-generic way. We applied these rules to all compounds in the ModelSEED database, predicting 180,891 spontaneous reactions. The resulting (CD)-MINE is available at https://minedatabase.mcs.anl.gov/cdmine/#/home and through developer tools. We also demonstrate how damage-prone intermediates and end products are widely distributed among metabolic pathways, and how predicting spontaneous chemical damage helps rationalize toxicity and carbon loss using examples from published pathways to commercial products. We explain how analyzing damage-prone areas in metabolism helps design effective engineering strategies. Finally, we use the CD-MINE toolset to predict the formation of the novel damage product N-carbamoyl proline, and present mass spectrometric evidence for its presence in Escherichia coli.

Synergism between SLC6A14 blockade and gemcitabine in pancreactic cancer: a 1H-NMR-based metabolomic study in pancreatic cancer cells

Gemcitabine is the first-line chemotherapy for pancreatic cancer. To overcome the often-acquired gemcitabine resistance, other drugs are used in combination with gemcitabine. It is well-known that cancer cells reprogram cellular metabolism, coupled with the up-regulation of selective nutrient transporters to feed into the altered metabolic pathways. Our previous studies have demonstrated that the amino acid transporter SLC6A14 is markedly up-regulated in pancreatic cancer and that it is a viable therapeutic target. α-Methyltryptophan (α-MT) is a blocker of SLC6A14 and is effective against pancreatic cancer in vitro and in vivo. In the present study, we tested the hypothesis that α-MT could synergize with gemcitabine in the treatment of pancreatic cancer. We investigated the effects of combination of α-MT and gemcitabine on proliferation, migration, and apoptosis in a human pancreatic cancer cell line, and examined the underlying mechanisms using 1H-NMR-based metabolomic analysis. These studies examined the intracellular metabolite profile and the extracellular metabolite profile separately. Combination of α-MT with gemcitabine elicited marked changes in a wide variety of metabolic pathways, particularly amino acid metabolism with notable alterations in pathways involving tryptophan, branched-chain amino acids, ketone bodies, and membrane phospholipids. The metabolomic profiles of untreated control cells and cells treated with gemcitabine or α-MT were distinctly separable, and the combination regimen showed a certain extent of overlap with the individual α-MT and gemcitabine groups. This represents the first study detailing the metabolomic basis of the anticancer efficacy of gemcitabine, α-MT and their combination.

NS-kNN: a modified k-nearest neighbors approach for imputing metabolomics data

INTRODUCTION

A common problem in metabolomics data analysis is the existence of a substantial number of missing values, which can complicate, bias, or even prevent certain downstream analyses. One of the most widely-used solutions to this problem is imputation of missing values using a k-nearest neighbors (kNN) algorithm to estimate missing metabolite abundances. kNN implicitly assumes that missing values are uniformly distributed at random in the dataset, but this is typically not true in metabolomics, where many values are missing because they are below the limit of detection of the analytical instrumentation.

OBJECTIVES

Here, we explore the impact of nonuniformly distributed missing values (missing not at random, or MNAR) on imputation performance. We present a new model for generating synthetic missing data and a new algorithm, No-Skip kNN (NS-kNN), that accounts for MNAR values to provide more accurate imputations.

METHODS

We compare the imputation errors of the original kNN algorithm using two distance metrics, NS-kNN, and a recently developed algorithm KNN-TN, when applied to multiple experimental datasets with different types and levels of missing data.

RESULTS

Our results show that NS-kNN typically outperforms kNN when at least 20-30% of missing values in a dataset are MNAR. NS-kNN also has lower imputation errors than KNN-TN on realistic datasets when at least 50% of missing values are MNAR.

CONCLUSION

Accounting for the nonuniform distribution of missing values in metabolomics data can significantly improve the results of imputation algorithms. The NS-kNN method imputes missing metabolomics data more accurately than existing kNN-based approaches when used on realistic datasets.

Metabolomics Analysis of the Toxic Effects of the Production of Lycopene and Its Precursors

Using cells as microbial factories enables highly specific production of chemicals with many advantages over chemical syntheses. A number of exciting new applications of this approach are in the area of precision metabolic engineering, which focuses on improving the specificity of target production. In recent work, we have used precision metabolic engineering to design lycopene-producing Escherichia coli for use as a low-cost diagnostic biosensor. To increase precursor availability and thus the rate of lycopene production, we heterologously expressed the mevalonate pathway. We found that simultaneous induction of these pathways increases lycopene production, but induction of the mevalonate pathway before induction of the lycopene pathway decreases both lycopene production and growth rate. Here, we aim to characterize the metabolic changes the cells may be undergoing during expression of either or both of these heterologous pathways. After establishing an improved method for quenching E. coli for metabolomics analysis, we used two-dimensional gas chromatography coupled to mass spectrometry (GCxGC-MS) to characterize the metabolomic profile of our lycopene-producing strains in growth conditions characteristic of our biosensor application. We found that the metabolic impacts of producing low, non-toxic levels of lycopene are of much smaller magnitude than the typical metabolic changes inherent to batch growth. We then used metabolomics to study differences in metabolism caused by the time of mevalonate pathway induction and the presence of the lycopene biosynthesis genes. We found that overnight induction of the mevalonate pathway was toxic to cells, but that the cells could recover if the lycopene pathway was not also heterologously expressed. The two pathways appeared to have an antagonistic metabolic effect that was clearly reflected in the cells’ metabolic profiles. The metabolites homocysteine and homoserine exhibited particularly interesting behaviors and may be linked to the growth inhibition seen when the mevalonate pathway is induced overnight, suggesting potential future work that may be useful in engineering increased lycopene biosynthesis.

High-Performance Chemical Isotope Labeling Liquid Chromatography-Mass Spectrometry for Profiling the Metabolomic Reprogramming Elicited by Ammonium Limitation in Yeast

Information about how yeast metabolism is rewired in response to internal and external cues can inform the development of metabolic engineering strategies for food, fuel, and chemical production in this organism. We report a new metabolomics workflow for the characterization of such metabolic rewiring. The workflow combines efficient cell lysis without using chemicals that may interfere with downstream sample analysis and differential chemical isotope labeling liquid chromatography mass spectrometry (CIL LC-MS) for in-depth yeast metabolome profiling. Using (12)C- and (13)C-dansylation (Dns) labeling to analyze the amine/phenol submetabolome, we detected and quantified a total of 5719 peak pairs or metabolites. Among them, 120 metabolites were positively identified using a library of 275 Dns-metabolite standards, and 2980 metabolites were putatively identified based on accurate mass matches to metabolome databases. We also applied (12)C- and (13)C-dimethylaminophenacyl (DmPA) labeling to profile the carboxylic acid submetabolome and detected over 2286 peak pairs, from which 33 metabolites were positively identified using a library of 188 DmPA-metabolite standards, and 1595 metabolites were putatively identified. Using this workflow for metabolomic profiling of cells challenged by ammonium limitation revealed unexpected links between ammonium assimilation and pantothenate accumulation that might be amenable to engineering for better acetyl-CoA production in yeast. We anticipate that efforts to improve other schemes of metabolic engineering will benefit from application of this workflow to multiple cell types.

Metabolomics, Standards, and Metabolic Modeling for Synthetic Biology in Plants

Life on earth depends on dynamic chemical transformations that enable cellular functions, including electron transfer reactions, as well as synthesis and degradation of biomolecules. Biochemical reactions are coordinated in metabolic pathways that interact in a complex way to allow adequate regulation. Biotechnology, food, biofuel, agricultural, and pharmaceutical industries are highly interested in metabolic engineering as an enabling technology of synthetic biology to exploit cells for the controlled production of metabolites of interest. These approaches have only recently been extended to plants due to their greater metabolic complexity (such as primary and secondary metabolism) and highly compartmentalized cellular structures and functions (including plant-specific organelles) compared with bacteria and other microorganisms. Technological advances in analytical instrumentation in combination with advances in data analysis and modeling have opened up new approaches to engineer plant metabolic pathways and allow the impact of modifications to be predicted more accurately. In this article, we review challenges in the integration and analysis of large-scale metabolic data, present an overview of current bioinformatics methods for the modeling and visualization of metabolic networks, and discuss approaches for interfacing bioinformatics approaches with metabolic models of cellular processes and flux distributions in order to predict phenotypes derived from specific genetic modifications or subjected to different environmental conditions.

MINEs: open access databases of computationally predicted enzyme promiscuity products for untargeted metabolomics

BACKGROUND

In spite of its great promise, metabolomics has proven difficult to execute in an untargeted and generalizable manner. Liquid chromatography-mass spectrometry (LC-MS) has made it possible to gather data on thousands of cellular metabolites. However, matching metabolites to their spectral features continues to be a bottleneck, meaning that much of the collected information remains uninterpreted and that new metabolites are seldom discovered in untargeted studies. These challenges require new approaches that consider compounds beyond those available in curated biochemistry databases.

DESCRIPTION

Here we present Metabolic In silico Network Expansions (MINEs), an extension of known metabolite databases to include molecules that have not been observed, but are likely to occur based on known metabolites and common biochemical reactions. We utilize an algorithm called the Biochemical Network Integrated Computational Explorer (BNICE) and expert-curated reaction rules based on the Enzyme Commission classification system to propose the novel chemical structures and reactions that comprise MINE databases. Starting from the Kyoto Encyclopedia of Genes and Genomes (KEGG) COMPOUND database, the MINE contains over 571,000 compounds, of which 93% are not present in the PubChem database. However, these MINE compounds have on average higher structural similarity to natural products than compounds from KEGG or PubChem. MINE databases were able to propose annotations for 98.6% of a set of 667 MassBank spectra, 14% more than KEGG alone and equivalent to PubChem while returning far fewer candidates per spectra than PubChem (46 vs. 1715 median candidates). Application of MINEs to LC-MS accurate mass data enabled the identity of an unknown peak to be confidently predicted.

CONCLUSIONS

MINE databases are freely accessible for non-commercial use via user-friendly web-tools at http://minedatabase.mcs.anl.gov and developer-friendly APIs. MINEs improve metabolomics peak identification as compared to general chemical databases whose results include irrelevant synthetic compounds. Furthermore, MINEs complement and expand on previous in silico generated compound databases that focus on human metabolism. We are actively developing the database; future versions of this resource will incorporate transformation rules for spontaneous chemical reactions and more advanced filtering and prioritization of candidate structures. Graphical abstractMINE database construction and access methods. The process of constructing a MINE database from the curated source databases is depicted on the left. The methods for accessing the database are shown on the right.