An integrated computational approach for metabolic flux analysis coupled with inference of tandem-MS collisional fragments

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.

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.

Efficient Modeling of MS/MS Data for Metabolic Flux Analysis

Metabolic flux analysis (MFA) is a widely used method for quantifying intracellular metabolic fluxes. It works by feeding cells with isotopic labeled nutrients, measuring metabolite isotopic labeling, and computationally interpreting the measured labeling data to estimate flux. Tandem mass-spectrometry (MS/MS) has been shown to be useful for MFA, providing positional isotopic labeling data. Specifically, MS/MS enables the measurement of a metabolite tandem mass-isotopomer distribution, representing the abundance in which certain parent and product fragments of a metabolite have different number of labeled atoms. However, a major limitation in using MFA with MS/MS data is the lack of a computationally efficient method for simulating such isotopic labeling data. Here, we describe the tandemer approach for efficiently computing metabolite tandem mass-isotopomer distributions in a metabolic network, given an estimation of metabolic fluxes. This approach can be used by MFA to find optimal metabolic fluxes, whose induced metabolite labeling patterns match tandem mass-isotopomer distributions measured by MS/MS. The tandemer approach is applied to simulate MS/MS data in a small-scale metabolic network model of mammalian methionine metabolism and in a large-scale metabolic network model of E. coli. It is shown to significantly improve the running time by between two to three orders of magnitude compared to the state-of-the-art, cumomers approach. We expect the tandemer approach to promote broader usage of MS/MS technology in metabolic flux analysis. Implementation is freely available at www.cs.technion.ac.il/~tomersh/methods.html

(13)C metabolic flux analysis of recombinant expression hosts

Identifying host cell metabolic phenotypes that promote high recombinant protein titer is a major goal of the biotech industry. (13)C metabolic flux analysis (MFA) provides a rigorous approach to quantify these metabolic phenotypes by applying isotope tracers to map the flow of carbon through intracellular metabolic pathways. Recent advances in tracer theory and measurements are enabling more information to be extracted from (13)C labeling experiments. Sustained development of publicly available software tools and standardization of experimental workflows is simultaneously encouraging increased adoption of (13)C MFA within the biotech research community. A number of recent (13)C MFA studies have identified increased citric acid cycle and pentose phosphate pathway fluxes as consistent markers of high recombinant protein expression, both in mammalian and microbial hosts. Further work is needed to determine whether redirecting flux into these pathways can effectively enhance protein titers while maintaining acceptable glycan profiles.