MID Max: LC-MS/MS Method for Measuring the Precursor and Product Mass Isotopomer Distributions of Metabolic Intermediates and Cofactors for Metabolic Flux Analysis Applications

Simultaneous in vivo multi-organ fluxomics reveals divergent metabolic adaptations in liver, heart, and skeletal muscle during obesity

We present an isotope-based metabolic flux analysis (MFA) approach to simultaneously quantify metabolic fluxes in the liver, heart, and skeletal muscle of individual mice. The platform was scaled to examine metabolic flux adaptations in age-matched cohorts of mice exhibiting varying levels of chronic obesity. We found that severe obesity increases hepatic gluconeogenesis and citric acid cycle flux, accompanied by elevated glucose oxidation in the heart that compensates for impaired fatty acid oxidation. In contrast, skeletal muscle fluxes exhibit an overall reduction in substrate oxidation. These findings demonstrate the dichotomy in fuel utilization between cardiac and skeletal muscle during worsening metabolic disease and demonstrate the divergent effects of obesity on metabolic fluxes in different organs. This multi-tissue MFA technology can be extended to address important questions about in vivo regulation of metabolism and its dysregulation in disease, which cannot be fully answered through studies of single organs or isolated cells/tissues.

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

Dual strategy for 13C-Metabolic flux analysis of central carbon and energy metabolism in Mammalian cells based on LC-isoMRM-MS

Central carbon and energy metabolism are the most concerned metabolic pathways in 13C-Metabolic flux analysis (13C-MFA). However, some α-keto acids, ribonucleoside triphosphate (NTPs) and deoxyribonucleoside triphosphate (dNTPs) involved in central carbon and energy metabolism pathways were unstable or reactive, leading to inaccurate metabolic flux analysis. To achieve accurate 13C-MFA of central carbon and energy metabolism, we proposed a dual strategy for the detection of 101 metabolites in glucose metabolism pathways. N-Methylphenylethylamine (MPEA) was utilized for derivatization of 4 carboxyl (α-keto acids) and 8 phosphate metabolites (NTPs and dNTPs). After derivatization, the MPEA derivatives were investigated to be stable for 4 weeks under 4 °C and detected with high intensity in ∼104 cells. On the other hand, we analyzed an additional 89 metabolites in central carbon and energy metabolic pathways were directly analyzed by liquid chromatography tandem mass spectrometry (LC-MRM-MS). The limit of detection (LODs) of our method were as low as 0.05 ng/mL and the linear range was at least two orders of magnitude with determination coefficient (R2) > 0.9701. The relative standard divisions (RSDs) of intra- and inter-day of 95% metabolites were below 20%. In addition, the isotope list of 82 detected metabolites in central carbon and energy metabolism were generated according to isotopologues and isotopomers for each metabolite resulting from 13C incorporation. Accurate assessment of mass isotopomer distributions (MIDs) of intracellular 13C-labeled metabolites was achieved in [U-13C]-glucose cultured HepG2 cells by our dual strategy. Finally, we performed MID analysis of 101 metabolites in central carbon and energy metabolism. Overall, this dual method is reproducible and robust for application on 13C-MFA and has a great potential for studying clinical isotope labeled samples.

Isotopically non-stationary 13 C metabolic flux analysis of two closely related fast-growing cyanobacteria, Synechococcus elongatus PCC 11801 and 11802

Synechococcus elongatus PCC 11801 and 11802 are closely related cyanobacterial strains that are fast-growing and tolerant to high light and temperature. These strains hold significant promise as chassis for photosynthetic production of chemicals from carbon dioxide. A detailed quantitative understanding of the central carbon pathways would be a reference for future metabolic engineering studies with these strains. We conducted isotopic non-stationary 13 C metabolic flux analysis to quantitively assess the metabolic potential of these two strains. This study highlights key similarities and differences in the central carbon flux distribution between these and other model/non-model strains. The two strains demonstrated a higher Calvin-Benson-Bassham (CBB) cycle flux coupled with negligible flux through the oxidative pentose phosphate pathway and the photorespiratory pathway and lower anaplerosis fluxes under photoautotrophic conditions. Interestingly, PCC 11802 shows the highest CBB cycle and pyruvate kinase flux values among those reported in cyanobacteria. The unique tricarboxylic acid (TCA) cycle diversion in PCC 11801 makes it ideal for the large-scale production of TCA cycle-derived chemicals. Additionally, dynamic labeling transients were measured for intermediates of amino acid, nucleotide, and nucleotide sugar metabolism. Overall, this study provides the first detailed metabolic flux maps of S. elongatus PCC 11801 and 11802, which may aid metabolic engineering efforts in these strains.

Automating the design-build-test-learn cycle towards next-generation bacterial cell factories

Automation is playing an increasingly significant role in synthetic biology. Groundbreaking technologies, developed over the past 20 years, have enormously accelerated the construction of efficient microbial cell factories. Integrating state-of-the-art tools (e.g. for genome engineering and analytical techniques) into the design-build-test-learn cycle (DBTLc) will shift the metabolic engineering paradigm from an almost artisanal labor towards a fully automated workflow. Here, we provide a perspective on how a fully automated DBTLc could be harnessed to construct the next-generation bacterial cell factories in a fast, high-throughput fashion. Innovative toolsets and approaches that pushed the boundaries in each segment of the cycle are reviewed to this end. We also present the most recent efforts on automation of the DBTLc, which heralds a fully autonomous pipeline for synthetic biology in the near future.

An automated workflow for multi-omics screening of microbial model organisms

Multi-omics datasets are becoming of key importance to drive discovery in fundamental research as much as generating knowledge for applied biotechnology. However, the construction of such large datasets is usually time-consuming and expensive. Automation might enable to overcome these issues by streamlining workflows from sample generation to data analysis. Here, we describe the construction of a complex workflow for the generation of high-throughput microbial multi-omics datasets. The workflow comprises a custom-built platform for automated cultivation and sampling of microbes, sample preparation protocols, analytical methods for sample analysis and automated scripts for raw data processing. We demonstrate possibilities and limitations of such workflow in generating data for three biotechnologically relevant model organisms, namely Escherichia coli, Saccharomyces cerevisiae, and Pseudomonas putida.

Reconstruction of Glutathione Metabolism in the Neuronal Model of Rotenone-Induced Neurodegeneration Using Mass Isotopologue Analysis with Hydrophilic Interaction Liquid Chromatography-Zeno High-Resolution Multiple Reaction Monitoring

Accurate reconstruction of metabolic pathways is an important prerequisite for interpreting metabolomics changes and understanding the diverse biological processes in disease models. A tracer-based metabolomics strategy utilizes stable isotope-labeled precursors to resolve complex pathways by tracing the labeled atom(s) to downstream metabolites through enzymatic reactions. Isotope enrichment analysis is informative and achieved by counting total labeled atoms and acquiring the mass isotopologue distribution (MID) of the intact metabolite. However, quantitative analysis of labeled metabolite substructures/moieties (MS² fragments) can offer more valuable insights into the reaction connections through measuring metabolite transformation. In order to acquire the isotopic labeling information at the intact metabolite and moiety level simultaneously, we developed a method that couples hydrophilic interaction liquid chromatography (HILIC) with Zeno trap-enabled high-resolution multiple reaction monitoring (MRM^HR). The method enabled accurate and reproducible MID quantification for intact metabolites as well as their fragmented moieties, with notably high sensitivity in the MS² fragmentation mode based on the measurement of ¹³C- or ¹⁵N-labeled cellular samples. The method was applied to human-induced pluripotent stem cell-derived neurons to trace the fate of ¹³C/¹⁵N atoms from D-¹³C₆-glucose/L-¹⁵N₂-glutamine added to the media. With the MID analysis of both intact metabolites and fragmented moieties, we validated the pathway reconstruction of de novo glutathione synthesis in mid-brain neurons. We discovered increased glutathione oxidization from both basal and newly synthesized glutathione pools under neuronal oxidative stress. Furthermore, the significantly decreased de novo glutathione synthesis was investigated and associated with altered activities of several key enzymes, as evidenced by suppressed glutamate supply via glucose metabolism and a diminished flux of glutathione synthetic reaction in the neuronal model of rotenone-induced neurodegeneration.

Metabolic flux analysis: a comprehensive review on sample preparation, analytical techniques, data analysis, computational modelling, and main application areas

Metabolic flux analysis (MFA) quantitatively describes cellular fluxes to understand metabolic phenotypes and functional behaviour after environmental and/or genetic perturbations. In the last decade, the application of stable isotopes became extremely important to determine and integrate in vivo measurements of metabolic reactions in systems biology. 13C-MFA is one of the most informative methods used to study central metabolism of biological systems. This review aims to outline the current experimental procedure adopted in 13C-MFA, starting from the preparation of cell cultures and labelled tracers to the quenching and extraction of metabolites and their subsequent analysis performed with very powerful software. Here, the limitations and advantages of nuclear magnetic resonance spectroscopy and mass spectrometry techniques used in carbon labelled experiments are elucidated by reviewing the most recent published papers. Furthermore, we summarise the most successful approaches used for computational modelling in flux analysis and the main application areas with a particular focus in metabolic engineering.

Isotope Calculation Gadgets: A Series of Software for Isotope-Tracing Experiments in Garuda Platform

Isotope tracing is a powerful technique for elucidating intracellular metabolism. Experiments utilizing this technique involve various processes, such as the correction of natural isotopes. Although some previously developed software are available for these procedures, there are still time-consuming steps in isotope tracing including the creation of an isotope measurement method in mass spectrometry (MS) and the interpretation of obtained labeling data. Additionally, these multi-step tasks often require data format conversion, which is also time-consuming. In this study, the Isotope Calculation Gadgets, a series of software that supports an entire workflow of isotope-tracing experiments, was developed in the Garuda platform, an open community. Garuda is a graphical user interface-based platform that allows individual operations to be sequentially performed, without data format conversion, which significantly reduces the required time and effort. The developed software includes new features that construct channels for isotopomer measurements, as well as conventional functions such as natural isotope correction, the calculation of fractional labeling and split ratio, and data mapping, thus facilitating an overall workflow of isotope-tracing experiments through smooth functional integration.

Dynamic flux regulation for high-titer anthranilate production by plasmid-free, conditionally-auxotrophic strains of Pseudomonas putida

Anthranilate, an intermediate of the shikimate pathway, is a high-value aromatic compound widely used as a precursor in the production of dyes, fragrances, plastics and pharmaceuticals. Traditional strategies adopted for microbial anthranilate production rely on the implementation of auxotrophic strains-which requires aromatic amino acids or complex additives to be supplemented in the culture medium, negatively impacting production costs. In this work, we engineered the soil bacterium Pseudomonas putida for high-titer, glucose-dependent anthranilate production by repurposing elements of the Esa quorum sensing (QS) system of Pantoea stewartii. The P_esaS promoter mediated a self-regulated transcriptional response that effectively knocked-down the expression of the trpDC genes. Next, we harnessed the synthetic QS elements to engineer a growth-to-anthranilate production switch. The resulting plasmid-free P. putida strain produced the target compound at 3.8 ± 0.3 mM in shaken-flask cultures after 72 h-a titer >2-fold higher than anthranilate levels reported thus far. Our results highlight the value of dynamic flux regulation for the production of intermediate metabolites within highly-regulated routes (such as the shikimate pathway), thereby circumventing the need of expensive additives.

An efficient LC-MS method for isomer separation and detection of sugars, phosphorylated sugars, and organic acids

Assessing central carbon metabolism in plants can be challenging due to the dynamic range in pool sizes, with low levels of important phosphorylated sugars relative to more abundant sugars and organic acids. Here, we report a sensitive liquid chromatography-mass spectrometry method for analysing central metabolites on a hybrid column, where both anion-exchange and hydrophilic interaction chromatography (HILIC) ligands are embedded in the stationary phase. The liquid chromatography method was developed for enhanced selectivity of 27 central metabolites in a single run with sensitivity at femtomole levels observed for most phosphorylated sugars. The method resolved phosphorylated hexose, pentose, and triose isomers that are otherwise challenging. Compared with a standard HILIC approach, these metabolites had improved peak areas using our approach due to ion enhancement or low ion suppression in the biological sample matrix. The approach was applied to investigate metabolism in high lipid-producing tobacco leaves that exhibited increased levels of acetyl-CoA, a precursor for oil biosynthesis. The application of the method to isotopologue detection and quantification was considered through evaluating 13C-labeled seeds from Camelina sativa. The method provides a means to analyse intermediates more comprehensively in central metabolism of plant tissues.

Model-guided dynamic control of essential metabolic nodes boosts acetyl-coenzyme A-dependent bioproduction in rewired Pseudomonas putida

Pseudomonas putida is evolutionarily endowed with features relevant for bioproduction, especially under harsh operating conditions. The rich metabolic versatility of this species, however, comes at the price of limited formation of acetyl-coenzyme A (CoA) from sugar substrates. Since acetyl-CoA is a key metabolic precursor for a number of added-value products, in this work we deployed an in silico-guided rewiring program of central carbon metabolism for upgrading P. putida as a host for acetyl-CoA-dependent bioproduction. An updated kinetic model, integrating fluxomics and metabolomics datasets in addition to manually-curated information of enzyme mechanisms, identified targets that would lead to increased acetyl-CoA levels. Based on these predictions, a set of plasmids based on clustered regularly interspaced short palindromic repeats (CRISPR) and dead CRISPR-associated protein 9 (dCas9) was constructed to silence genes by CRISPR interference (CRISPRi). Dynamic reduction of gene expression of two key targets (gltA, encoding citrate synthase, and the essential accA gene, encoding subunit A of the acetyl-CoA carboxylase complex) mediated an 8-fold increase in the acetyl-CoA content of rewired P. putida. Poly(3-hydroxybutyrate) (PHB) was adopted as a proxy of acetyl-CoA availability, and two synthetic pathways were engineered for biopolymer accumulation. By including cell morphology as an extra target for the CRISPRi approach, fully rewired P. putida strains programmed for PHB accumulation had a 5-fold increase in PHB titers in bioreactor cultures using glucose. Thus, the strategy described herein allowed for rationally redirecting metabolic fluxes in P. putida from central metabolism towards product biosynthesis-especially relevant when deletion of essential pathways is not an option.

A guide to metabolic flux analysis in metabolic engineering: Methods, tools and applications

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, "Metabolic fluxes and metabolic engineering" (Metabolic Engineering, 1: 1-11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. ¹³C-metabolic flux analysis (¹³C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.

Comprehensive Isotopic Targeted Mass Spectrometry: Reliable Metabolic Flux Analysis with Broad Coverage

Metabolic flux analysis (MFA) is highly relevant to understanding metabolic mechanisms of various biological processes. While the pace of methodology development in MFA has been rapid, a major challenge the field continues to witness is limited metabolite coverage, often restricted to a small to moderate number of well-known compounds. In addition, isotopic peaks from an enriched metabolite tend to have low abundances, which makes liquid chromatography tandem mass spectrometry (LC-MS/MS) highly useful in MFA due to its high sensitivity and specificity. Previously we have built large-scale LC-MS/MS approaches that can be routinely used for measurement of up to ∼1,900 metabolite/feature levels [Gu et al. Anal. Chem. 2015, 87, 12355-12362. Shi et al. Anal. Chem. 2019, 91, 13737-13745.]. In this study, we aim to expand our previous studies focused on metabolite level measurements to flux analysis and establish a novel comprehensive isotopic targeted mass spectrometry (CIT-MS) method for reliable MFA analysis with broad coverage. As a proof-of-principle, we have applied CIT-MS to compare the steady-state enrichment of metabolites between Myc(oncogene)-On and Myc-Off Tet21N human neuroblastoma cells cultured with U-13C6-glucose medium. CIT-MS is operationalized using multiple reaction monitoring (MRM) mode and is able to perform MFA of 310 identified metabolites (142 reliably detected, 46 kinetically profiled) selected from >35 metabolic pathways of strong biological significance. Further, we developed a novel concept of relative flux, which eliminates the requirement of absolute quantitation in traditional MFA and thus enables comparative MFA under the pseudosteady state. As a result, CIT-MS was shown to possess the advantages of broad coverage, easy implementation, fast throughput, and more importantly, high fidelity and accuracy in MFA. In principle, CIT-MS can be easily adapted to track the flux of other labeled tracers (such as 15N-tracers) in any metabolite detectable by LC-MS/MS and in various biological models (such as mice). Therefore, CIT-MS has great potential to bring new insights to both basic and clinical metabolism research.

Enzymes as Parts in Need of Replacement - and How to Extend Their Working Life

Enzymes catalyze reactions in vivo at different rates and each enzyme molecule has a lifetime limit before it is degraded and replaced to enable catalysis to continue. Considering these rates together as a unitless ratio of catalytic cycles until replacement (CCR) provides a new quantitative tool to assess the replacement schedule of and energy investment into enzymes as they relate to function. Here, we outline the challenges of determining CCRs and new approaches to overcome them and then assess the CCRs of selected enzymes in bacteria and plants to reveal a range of seven orders of magnitude for this ratio. Modifying CCRs in plants holds promise to lower cellular costs, to tailor enzymes for particular environments, and to breed enzyme improvements for crop productivity.

Fragmentation of Dicarboxylic and Tricarboxylic Acids in the Krebs Cycle Using GC-EI-MS and GC-EI-MS/MS

Isotope labeling measurements using mass spectrometry can provide informative insights on the metabolic systems of various organisms. The detailed identification of carbon positions included in the fragment ions of dicarboxylic and tricarboxylic acids in central carbon metabolism is needed for precise interpretation of the metabolic states. In this study, fragment ions containing the carbon backbone cleavage of dicarboxylic and tricarboxylic in the Krebs cycle were investigated by using gas chromatography (GC)-electron ionization (EI)-MS and GC-EI-MS/MS. The positions of decarboxylation in the dicarboxylic and tricarboxylic acids were successfully identified by analyses using position-specific ¹³C-labeled standards prepared by in vitro enzymatic reactions. For example, carboxyl groups of C1 and C6 of trimethylsilyl (TMS)- and tert-butyldimethylsilyl (TBDMS)-derivatized malic and citric acids were primarily cleaved by EI. MS/MS analyses were also performed, and fragment ions of TBDMS-citric and α-ketoglutaric acids (αKG) with the loss of two carboxyl groups in collision-induced dissociation (CID) were observed.

SWATH: A Data-Independent Tandem Mass Spectrometry Method to Quantify 13C Enrichment in Cellular Metabolites and Fragments

Recently, the sequential windowed acquisition of all theoretical fragment ion mass spectra (SWATH) method coupled with liquid chromatography has been demonstrated for the quantification of isotopic ¹³C enrichment in a large number of cellular metabolites and fragments. SWATH, a data-independent acquisition (DIA) method, alleviates the need for data deconvolution and shows greater accuracy in the quantification of low abundance isotopologs of fragments thereby resulting in a lower systematic error. Here we provide a detailed protocol for the design of Q1 mass isolation windows and the post-acquisition data analysis with emphasis on the untargeted nature of SWATH.

Determination of Isotopologue and Tandem Mass Isotopologue Ratios Using Gas Chromatography Chemical Ionization Time of Flight Mass Spectrometry - Methodology and Uncertainty of Measurement

The accurate and precise analysis of isotopologue and tandem mass isotopologue ratios in heavy stable isotope labeling experiments is a critical part of assessing absolute intracellular metabolic fluxes. Resulting from feeding the organism of interest with a specifically isotope-labeled substrate, the principal characteristics of these labeling experiments are the metabolites' non-naturally distributed isotopologue patterns. For the purpose of inferring metabolic rates by maximizing the fit between a priori simulated and experimentally obtained labeling patterns, ¹³C is the preferred stable isotope of use.The analysis of the obtained labeling patterns can be approached by different mass spectrometric approaches. Gas chromatography (GC) features broad metabolite coverage and excellent separation efficiency of biologically relevant isomers. These advantages compensate for laborious derivatization steps and the resulting need for interference correction for natural abundant isotopes.Here, we describe a workflow based on GC-high resolution mass spectrometry with chemical ionization for the analysis of carbon-isotopologue distributions and some positional labeling information of primary metabolites. To study the associated measurement uncertainty of the resulting ¹³C labeling patterns, guidance to uncertainty estimation according to the EURACHEM guidelines with Monte-Carlo simulation is provided.

High-resolution 13C metabolic flux analysis

Precise quantification of metabolic pathway fluxes in biological systems is of major importance in guiding efforts in metabolic engineering, biotechnology, microbiology, human health, and cell culture. ¹³C metabolic flux analysis (¹³C-MFA) is the predominant technique used for determining intracellular fluxes. Here, we present a protocol for ¹³C-MFA that incorporates recent advances in parallel labeling experiments, isotopic labeling measurements, and statistical analysis, as well as best practices developed through decades of experience. Experimental design to ensure that fluxes are estimated with the highest precision is an integral part of the protocol. The protocol is based on growing microbes in two (or more) parallel cultures with ¹³C-labeled glucose tracers, followed by gas chromatography-mass spectrometry (GC-MS) measurements of isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose. Fluxes are then estimated using software for ¹³C-MFA, such as Metran, followed by comprehensive statistical analysis to determine the goodness of fit and calculate confidence intervals of fluxes. The presented protocol can be completed in 4 d and quantifies metabolic fluxes with a standard deviation of ≤2%, a substantial improvement over previous implementations. The presented protocol is exemplified using an Escherichia coli ΔtpiA case study with full supporting data, providing a hands-on opportunity to step through a complex troubleshooting scenario. Although applications to prokaryotic microbial systems are emphasized, this protocol can be easily adjusted for application to eukaryotic organisms.

Regulatory mechanisms underlying coordination of amino acid and glucose catabolism in Escherichia coli

How microbes dynamically coordinate uptake and simultaneous utilization of nutrients in complex nutritional ecosystems is still an open question. Here, we develop a constraint-based modeling approach that exploits non-targeted exo-metabolomics data to unravel adaptive decision-making processes in dynamic nutritional environments. We thereby investigate metabolic adaptation of Escherichia coli to continuously changing conditions during batch growth in complex medium. Unexpectedly, model-based analysis of time resolved exo-metabolome data revealed that fastest growth coincides with preferred catabolism of amino acids, which, in turn, reduces glucose uptake and increases acetate overflow. We show that high intracellular levels of the amino acid degradation metabolites pyruvate and oxaloacetate can directly inhibit the phosphotransferase system (PTS), and reveal their functional role in mediating regulatory decisions for uptake and catabolism of alternative carbon sources. Overall, the proposed methodology expands the spectrum of possible applications of flux balance analysis to decipher metabolic adaptation mechanisms in naturally occurring habitats and diverse organisms.

The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) is the method of choice when a detailed inference of intracellular metabolic fluxes in living organisms under metabolic quasi-steady state conditions is desired. Being continuously developed since two decades, the technology made major contributions to the quantitative characterization of organisms in all fields of biotechnology and health-related research. ¹³C MFA, however, stands out from other "-omics sciences," in that it requires not only experimental-analytical data, but also mathematical models and a computational toolset to infer the quantities of interest, i.e., the metabolic fluxes. At present, these models cannot be conveniently exchanged between different labs. Here, we present the implementation-independent model description language FluxML for specifying ¹³C MFA models. The core of FluxML captures the metabolic reaction network together with atom mappings, constraints on the model parameters, and the wealth of data configurations. In particular, we describe the governing design processes that shaped the FluxML language. We demonstrate the utility of FluxML to represent many contemporary experimental-analytical requirements in the field of ¹³C MFA. The major aim of FluxML is to offer a sound, open, and future-proof language to unambiguously express and conserve all the necessary information for model re-use, exchange, and comparison. Along with FluxML, several powerful computational tools are supplied for easy handling, but also to maintain a maximum of flexibility. Altogether, the FluxML collection is an "all-around carefree package" for ¹³C MFA modelers. We believe that FluxML improves scientific productivity as well as transparency and therewith contributes to the efficiency and reproducibility of computational modeling efforts in the field of ¹³C MFA.

GC-QTOFMS with a low-energy electron ionization source for advancing isotopologue analysis in 13C-based metabolic flux analysis

For the study of different levels of (intra)cellular regulation and condition-dependent insight into metabolic activities, fluxomics experiments based on stable isotope tracer experiments using ¹³C have become a well-established approach. The experimentally obtained non-naturally distributed ¹³C labeling patterns of metabolite pools can be measured by mass spectrometric detection with front-end separation and can be consequently incorporated into biochemical network models. Here, despite a tedious derivatization step, gas chromatographic separation of polar metabolites is favorable because of the wide coverage range and high isomer separation efficiency. However, the typically employed electron ionization energy of 70 eV leads to significant fragmentation and consequently only low-abundant ions with an intact carbon backbone. Since these ions are considered a prerequisite for the analysis of the non-naturally distributed labeling patterns and further integration into modeling strategies, a softer ionization technique is needed. In the present work, a novel low energy electron ionization source is optimized for the analysis of primary metabolites and compared with a chemical ionization approach in terms of trueness, precision, and sensitivity.

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism

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 function of metabolic genes in the context of regulatory networks is not well understood. Here, the authors investigate the adaptive responses of E. coli after knockout of metabolic genes and highlight the influence of metabolite levels in the evolution of regulatory function.

Multiple Optimal Phenotypes Overcome Redox and Glycolytic Intermediate Metabolite Imbalances in Escherichia coli pgi Knockout Evolutions

A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.

Adaptation to the coupling of glycolysis to toxic methylglyoxal production in tpiA deletion strains of Escherichia coli requires synchronized and counterintuitive genetic changes

Methylglyoxal is a highly toxic metabolite that can be produced in all living organisms. Methylglyoxal was artificially elevated by removal of the tpiA gene from a growth optimized Escherichia coli strain. The initial response to elevated methylglyoxal and its toxicity was characterized, and detoxification mechanisms were studied using adaptive laboratory evolution. We found that: 1) Multi-omics analysis revealed biological consequences of methylglyoxal toxicity, which included attack on macromolecules including DNA and RNA and perturbation of nucleotide levels; 2) Counter-intuitive cross-talk between carbon starvation and inorganic phosphate signalling was revealed in the tpiA deletion strain that required mutations in inorganic phosphate signalling mechanisms to alleviate; and 3) The split flux through lower glycolysis depleted glycolytic intermediates requiring a host of synchronized and coordinated mutations in non-intuitive network locations in order to re-adjust the metabolic flux map to achieve optimal growth. Such mutations included a systematic inactivation of the Phosphotransferase System (PTS) and alterations in cell wall biosynthesis enzyme activity. This study demonstrated that deletion of major metabolic genes followed by ALE was a productive approach to gain novel insight into the systems biology underlying optimal phenotypic states.

SWATH Tandem Mass Spectrometry Workflow for Quantification of Mass Isotopologue Distribution of Intracellular Metabolites and Fragments Labeled with Isotopic 13C Carbon

Accurate quantification of mass isotopologue distribution (MID) of metabolites is a prerequisite for ¹³C-metabolic flux analysis. Currently used mass spectrometric (MS) techniques based on multiple reaction monitoring (MRM) place limitations on the number of MIDs that can be analyzed in a single run. Moreover, the deconvolution step results in amplification of error. Here, we demonstrate that SWATH MS/MS, a data independent acquisition (DIA) technique allows quantification of a large number of precursor and product MIDs in a single run. SWATH sequentially fragments all precursor ions in stacked mass isolation windows. Co-fragmentation of all precursor isotopologues in a single SWATH window yields higher sensitivity enabling quantification of MIDs of fragments with low abundance and lower systematic and random errors. We quantify the MIDs of 53 precursor and product ions corresponding to 19 intracellular metabolites from a dynamic ¹³C-labeling of a model cyanobacterium, Synechococcus sp. PCC 7002. The use of product MIDs resulted in an improved precision of many measured fluxes compared to when only precursor MIDs were used for flux analysis. The approach is truly untargeted and allows additional metabolites to be quantified from the same data.

RapidRIP quantifies the intracellular metabolome of 7 industrial strains of E. coli

Fast metabolite quantification methods are required for high throughput screening of microbial strains obtained by combinatorial or evolutionary engineering approaches. In this study, a rapid RIP-LC-MS/MS (RapidRIP) method for high-throughput quantitative metabolomics was developed and validated that was capable of quantifying 102 metabolites from central, amino acid, energy, nucleotide, and cofactor metabolism in less than 5 minutes. The method was shown to have comparable sensitivity and resolving capability as compared to a full length RIP-LC-MS/MS method (FullRIP). The RapidRIP method was used to quantify the metabolome of seven industrial strains of E. coli revealing significant differences in glycolytic, pentose phosphate, TCA cycle, amino acid, and energy and cofactor metabolites were found. These differences translated to statistically and biologically significant differences in thermodynamics of biochemical reactions between strains that could have implications when choosing a host for bioprocessing.

To be certain about the uncertainty: Bayesian statistics for 13 C metabolic flux analysis

¹³ C Metabolic Fluxes Analysis (¹³ C MFA) remains to be the most powerful approach to determine intracellular metabolic reaction rates. Decisions on strain engineering and experimentation heavily rely upon the certainty with which these fluxes are estimated. For uncertainty quantification, the vast majority of ¹³ C MFA studies relies on confidence intervals from the paradigm of Frequentist statistics. However, it is well known that the confidence intervals for a given experimental outcome are not uniquely defined. As a result, confidence intervals produced by different methods can be different, but nevertheless equally valid. This is of high relevance to ¹³ C MFA, since practitioners regularly use three different approximate approaches for calculating confidence intervals. By means of a computational study with a realistic model of the central carbon metabolism of E. coli, we provide strong evidence that confidence intervals used in the field depend strongly on the technique with which they were calculated and, thus, their use leads to misinterpretation of the flux uncertainty. In order to provide a better alternative to confidence intervals in ¹³ C MFA, we demonstrate that credible intervals from the paradigm of Bayesian statistics give more reliable flux uncertainty quantifications which can be readily computed with high accuracy using Markov chain Monte Carlo. In addition, the widely applied chi-square test, as a means of testing whether the model reproduces the data, is examined closer.

Review of emerging metabolomic tools and resources: 2015-2016

Data processing and analysis are major bottlenecks in high-throughput metabolomic experiments. Recent advancements in data acquisition platforms are driving trends toward increasing data size (e.g., petabyte scale) and complexity (multiple omic platforms). Improvements in data analysis software and in silico methods are similarly required to effectively utilize these advancements and link the acquired data with biological interpretations. Herein, we provide an overview of recently developed and freely available metabolomic tools, algorithms, databases, and data analysis frameworks. This overview of popular tools for MS and NMR-based metabolomics is organized into the following sections: data processing, annotation, analysis, and visualization. The following overview of newly developed tools helps to better inform researchers to support the emergence of metabolomics as an integral tool for the study of biochemistry, systems biology, environmental analysis, health, and personalized medicine.

Implementation of data-dependent isotopologue fragmentation in 13C-based metabolic flux analysis

A novel analytical approach based on liquid chromatography coupled to quadrupole time of flight mass spectrometry, employing data-dependent triggering for analysis of isotopologue and tandem mass isotopomer fractions of metabolites of the primary carbon metabolism was developed. The implemented QTOFMS method employs automated MS/MS triggering of higher abundant, biologically relevant isotopologues for generating positional information of the respective metabolite. Using this advanced isotopologue selective fragmentation approach enables the generation of significant tandem mass isotopomer data within a short cycle time without compromising sensitivity. Due to a lack of suitable reference material certified for isotopologue ratios, a Pichia pastoris cell extract with a defined 13C distribution as well as a cell extract from a 13C-based metabolic flux experiment were employed for proof of concept. Moreover, a method inter-comparison with an already established GC-CI-(Q)TOFMS approach was conducted. Both methods showed good agreement on isotopologue and tandem mass isotopomer distributions for the two different cell extracts. Graphical abstract Schematic overview of data-dependent isotopologue fragmentation for acquisition of isotopologue and tandem mass isotopomer fractions.

Comprehensive and accurate tracking of carbon origin of LC-tandem mass spectrometry collisional fragments for 13C-MFA

In recent years the benefit of measuring positionally resolved 13C-labeling enrichment from tandem mass spectrometry (MS/MS) collisional fragments for improved precision of 13C-Metabolic Flux Analysis (13C-MFA) has become evident. However, the usage of positional labeling information for 13C-MFA faces two challenges: (1) The mass spectrometric acquisition of a large number of potentially interfering mass transitions may hamper accuracy and sensitivity. (2) The positional identity of carbon atoms of product ions needs to be known. The present contribution addresses the latter challenge by deducing the maximal positional labeling information contained in LC-ESI-MS/MS spectra of product anions of central metabolism as well as product cations of amino acids. For this purpose, we draw on accurate mass spectrometry, selectively labeled standards, and published fragmentation pathways to structurally annotate all dominant mass peaks of a large collection of metabolites, some of which with a complete fragmentation pathway. Compiling all available information, we arrive at the most detailed map of carbon atom fate of LC-ESI-MS/MS collisional fragments yet, comprising 170 intense and structurally annotated product ions with unique carbon origin from 76 precursor ions of 72 metabolites. Our 13C-data proof that heuristic fragmentation rules often fail to yield correct fragment structures and we expose common pitfalls in the structural annotation of product ions. We show that the positionally resolved 13C-label information contained in the product ions that we structurally annotated allows to infer the entire isotopomer distribution of several central metabolism intermediates, which is experimentally demonstrated for malate using quadrupole-time-of-flight MS technology. Finally, the inclusion of the label information from a subset of these fragments improves flux precision in a Corynebacterium glutamicum model of the central carbon metabolism.

Tracing metabolism from lignocellulosic biomass and gaseous substrates to products with stable-isotopes

Engineered microbes offer a practical and sustainable alternative to traditional industrial approaches. To increase the economic feasibility of biological processes, microbial isolates are engineered to take up inexpensive feedstocks (including lignocellulosic biomass, syngas, methane, and carbon dioxide), and convert them into substrates of central metabolism and further into value-added products. To trace the metabolism of these feedstocks into products, isotopic tracers are applied together with isotopomer analysis techniques such as ¹³C-metabolic flux analysis to provide a detailed picture of pathway utilization. Flux data is then integrated with kinetic models and constraint-based approaches to identify metabolic bottlenecks, propose novel metabolic engineering strategies, and improve process performance.

Comprehensive metabolic modeling of multiple 13C-isotopomer data sets to study metabolism in perfused working hearts

In many forms of cardiomyopathy, alterations in energy substrate metabolism play a key role in disease pathogenesis. Stable isotope tracing in rodent heart perfusion systems can be used to determine cardiac metabolic fluxes, namely those relative fluxes that contribute to pyruvate, the acetyl-CoA pool, and pyruvate anaplerosis, which are critical to cardiac homeostasis. Methods have previously been developed to interrogate these relative fluxes using isotopomer enrichments of measured metabolites and algebraic equations to determine a predefined metabolic flux model. However, this approach is exquisitely sensitive to measurement error, thus precluding accurate relative flux parameter determination. In this study, we applied a novel mathematical approach to determine relative cardiac metabolic fluxes using ¹³C-metabolic flux analysis (¹³C-MFA) aided by multiple tracer experiments and integrated data analysis. Using ¹³C-MFA, we validated a metabolic network model to explain myocardial energy substrate metabolism. Four different ¹³C-labeled substrates were queried (i.e., glucose, lactate, pyruvate, and oleate) based on a previously published study. We integrated the analysis of the complete set of isotopomer data gathered from these mouse heart perfusion experiments into a single comprehensive network model that delineates substrate contributions to both pyruvate and acetyl-CoA pools at a greater resolution than that offered by traditional methods using algebraic equations. To our knowledge, this is the first rigorous application of ¹³C-MFA to interrogate data from multiple tracer experiments in the perfused heart. We anticipate that this approach can be used widely to study energy substrate metabolism in this and other similar biological systems.

Current state and challenges for dynamic metabolic modeling

While the stoichiometry of metabolism is probably the best studied cellular level, the dynamics in metabolism can still not be well described, predicted and, thus, engineered. Unknowns in the metabolic flux behavior arise from kinetic interactions, especially allosteric control mechanisms. While the stoichiometry of enzymes is preserved in vitro, their activity and kinetic behavior differs from the in vivo situation. Next to this challenge, it is infeasible to test the interaction of each enzyme with each intracellular metabolite in vitro exhaustively. As a consequence, the whole interacting metabolome has to be studied in vivo to identify the relevant enzymes properties. In this review we discuss current approaches for in vivo perturbation experiments, that is, stimulus response experiments using different setups and quantitative analytical approaches, including dynamic carbon tracing. Next to reliable and informative data, advanced modeling approaches and computational tools are required to identify kinetic mechanisms and their parameters.

GC-MS/MS survey of collision-induced dissociation of tert-butyldimethylsilyl-derivatized amino acids and its application to (13)C-metabolic flux analysis of Escherichia coli central metabolism

Stable isotope labeling experiments using mass spectrometry have been employed to investigate carbon flow levels (metabolic flux) in mammalian, plant, and microbial cells. To achieve a more precise (13)C-metabolic flux analysis ((13)C-MFA), novel fragmentations of tert-butyldimethylsilyl (TBDMS)-amino acids were investigated by gas chromatography-tandem mass spectrometry (GC-MS/MS). The product ion scan analyses of 15 TBDMS-amino acids revealed 24 novel fragment ions. The amino acid-derived carbons included in the five fragment ions were identified by the analyses of (13)C-labeled authentic standards. The identification of the fragment ion at m/z 170 indicated that the isotopic abundance of S-methyl carbon in methionine could be determined from the cleavage of C5 in the precursor of [M-159](+) (m/z 218). It was also confirmed that the precision of (13)C-MFA in Escherichia coli central carbon metabolism could be improved by introducing (13)C-labeling data derived from novel fragmentations. Graphical Abstract Novel collision-induced dissociation fragmentations of tert-butyldimethylsilyl amino acids were investigated and identified by GC-MS/MS.

Optimal tracers for parallel labeling experiments and 13C metabolic flux analysis: A new precision and synergy scoring system

¹³C-Metabolic flux analysis (¹³C-MFA) is a widely used approach in metabolic engineering for quantifying intracellular metabolic fluxes. The precision of fluxes determined by ¹³C-MFA depends largely on the choice of isotopic tracers and the specific set of labeling measurements. A recent advance in the field is the use of parallel labeling experiments for improved flux precision and accuracy. However, as of today, no systemic methods exist for identifying optimal tracers for parallel labeling experiments. In this contribution, we have addressed this problem by introducing a new scoring system and evaluating thousands of different isotopic tracer schemes. Based on this extensive analysis we have identified optimal tracers for ¹³C-MFA. The best single tracers were doubly ¹³C-labeled glucose tracers, including [1,6-¹³C]glucose, [5,6-¹³C]glucose and [1,2-¹³C]glucose, which consistently produced the highest flux precision independent of the metabolic flux map (here, 100 random flux maps were evaluated). Moreover, we demonstrate that pure glucose tracers perform better overall than mixtures of glucose tracers. For parallel labeling experiments the optimal isotopic tracers were [1,6-¹³C]glucose and [1,2-¹³C]glucose. Combined analysis of [1,6-¹³C]glucose and [1,2-¹³C]glucose labeling data improved the flux precision score by nearly 20-fold compared to widely use tracer mixture 80% [1-¹³C]glucose +20% [U-¹³C]glucose.