Isotopomer analysis using GC-MS

Optimised collection of non-uniformly sampled 2D-HSQC NMR spectra for use in metabolic flux analysis

Nuclear magnetic resonance (NMR) spectroscopy is integral to metabolic studies; yet, it can suffer from the long acquisition times required to collect data of sufficient signal strength and resolution. The use of non-uniform sampling (NUS) allows faster collection of NMR spectra without loss of spectral integrity. When planning experimental methodologies to perform metabolic flux analysis (MFA) of cell metabolism, a variety of options are available for the acquisition of NUS NMR data. Before beginning data collection, decisions have to be made regarding selection of pulse sequence, number of transients and NUS specific parameters such as the sampling level and sampling schedule. Poor choices will impact data quality, which may have a negative effect on the subsequent analysis and biological interpretation. Herein, we describe factors that should be considered when setting up non-uniformly sampled 2D-¹ H,¹³ C HSQC NMR experiments for MFA and provide a standard protocol for users to follow.

Parallel isotope differential modeling for instationary 13C fluxomics at the genome scale

BACKGROUND

A precise map of the metabolic fluxome, the closest surrogate to the physiological phenotype, is becoming progressively more important in the metabolic engineering of photosynthetic organisms for biofuel and biomass production. For photosynthetic organisms, the state-of-the-art method for this purpose is instationary 13C fluxomics, which has arisen as a sibling of transcriptomics or proteomics. Instationary 13C data processing requires solving high-dimensional nonlinear differential equations and leads to large computational and time costs when its scope is expanded to a genome-scale metabolic network.

RESULT

Here, we present a parallelized method to model instationary 13C labeling data. The elementary metabolite unit (EMU) framework is reorganized to allow treating individual mass isotopomers and breaking up of their networks into strongly connected components (SCCs). A variable domain parallel algorithm is introduced to process ordinary differential equations in a parallel way. 15-fold acceleration is achieved for constant-step-size modeling and ~ fivefold acceleration for adaptive-step-size modeling.

CONCLUSION

This algorithm is universally applicable to isotope granules such as EMUs and cumomers and can substantially accelerate instationary 13C fluxomics modeling. It thus has great potential to be widely adopted in any instationary 13C fluxomics modeling.

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.

Advances in metabolic flux analysis toward genome-scale profiling of higher organisms

Methodological and technological advances have recently paved the way for metabolic flux profiling in higher organisms, like plants. However, in comparison with omics technologies, flux profiling has yet to provide comprehensive differential flux maps at a genome-scale and in different cell types, tissues, and organs. Here we highlight the recent advances in technologies to gather metabolic labeling patterns and flux profiling approaches. We provide an opinion of how recent local flux profiling approaches can be used in conjunction with the constraint-based modeling framework to arrive at genome-scale flux maps. In addition, we point at approaches which use metabolomics data without introduction of label to predict either non-steady state fluxes in a time-series experiment or flux changes in different experimental scenarios. The combination of these developments allows an experimentally feasible approach for flux-based large-scale systems biology studies.

The rate of lactate production from glucose in hearts is not altered by per-deuteration of glucose

This study was designed to determine whether perdeuterated glucose experiences a kinetic isotope effect (KIE) as glucose passes through glycolysis and is further oxidized in the tricarboxylic acid (TCA) cycle. Metabolism of deuterated glucose was investigated in two groups of perfused rat hearts. The control group was supplied with a 1:1 mixture of [U-¹³C₆]glucose and [1,6-¹³C₂]glucose, while the experimental group received [U-¹³C₆,U-²H₇]glucose and [1,6-¹³C₂]glucose. Tissue extracts were analyzed by ¹H, ²H and proton-decoupled ¹³C NMR spectroscopy. Extensive ²H-¹³C scalar coupling plus chemical shift isotope effects were observed in the proton-decoupled ¹³C NMR spectra of lactate, alanine and glutamate. A small but measureable (∼8%) difference in the rate of conversion of [U-¹³C₆]glucose vs. [1,6-¹³C₂]glucose to lactate, likely reflecting rates of CC bond breakage in the aldolase reaction, but conversion of [U-¹³C₆]glucose versus [U-¹³C₆,U-²H₇]glucose to lactate did not differ. This shows that the presence of deuterium in glucose does not alter glycolytic flux. However, there were two distinct effects of deuteration on metabolism of glucose to alanine and oxidation of glucose in the TCA. First, alanine undergoes extensive exchange of methyl deuterons with solvent protons in the alanine amino transferase reaction. Second, there is a substantial kinetic isotope effect in metabolism of [U-¹³C₆,U-²H₇]glucose to alanine and glutamate. In the presence of [U-¹³C₆,U-²H₇]glucose, alanine and lactate are not in rapid exchange with the same pool of pyruvate. These studies indicate that the appearance of hyperpolarized ¹³C-lactate from hyperpolarized [U-¹³C₆,U-²H₇]glucose is not substantially influenced by a deuterium kinetic isotope effect.

Quantitative 13 C-isotope labelling-based analysis to elucidate the influence of environmental parameters on the production of fermentative aromas during wine fermentation

Nitrogen and lipids are key nutrients of grape must that influence the production of fermentative aromas by wine yeast, and we have previously shown that a strong interaction exists between these two nutrients. However, more than 90% of the acids and higher alcohols (and their acetate ester derivatives) were derived from intermediates produced by the carbon central metabolism (CCM). The objective of this study was to determine how variations in nitrogen and lipid resources can modulate the contribution of nitrogen and carbon metabolisms for the production of fermentative aromas. A quantitative analysis of metabolism using ¹³C‐labelled leucine and valine showed that nitrogen availability affected the part of the catabolism of N‐containing compounds, the formation of α‐ketoacids from CCM and the redistribution of fluxes around these precursors, explaining the optimum production of higher alcohols occurring at an intermediate nitrogen content. Moreover, nitrogen content modulated the total production of acids and higher alcohols differently, through variations in the redox state of cells. We also demonstrated that the phytosterol content, modifying the intracellular availability of acetyl‐CoA, can influence the flux distribution, especially the formation of higher alcohols and the conversion of α‐ketoisovalerate to α‐ketoisocaproate.

Accurate and sensitive determination of molar fractions of 13C-Labeled intracellular metabolites in cell cultures grown in the presence of isotopically-labeled glucose

This work describes a methodology based on multiple linear regression and GC-MS for the determination of molar fractions of isotopically-labeled intracellular metabolites in cell cultures. Novel aspects of this work are: i) the calculation of theoretical isotopic distributions of the different isotopologues from an experimentally measured value of % 13C enrichment of the labeled precursor ii) the calculation of the contribution of lack of mass resolution of the mass spectrometer and different fragmentation mechanism such as the loss or gain of hydrogen atoms in the EI source to measure the purity of the selected cluster for each metabolite and iii) the validation of the methodology not only by the analysis of gravimetrically prepared mixtures of isotopologues but also by the comparison of the obtained molar fractions with experimental values obtained by GC-Combustion-IRMS based on ¹³C/¹²C isotope ratio measurements. The method is able to measure molar fractions for twenty-eight intracellular metabolites derived from glucose metabolism in cell cultures grown in the presence of ¹³C-labeled Glucose. The validation strategies demonstrate a satisfactory accuracy and precision of the proposed procedure. Also, our results show that the minimum value of ¹³C incorporation that can be accurately quantified is significantly influenced by the calculation of the spectral purity of the measured cluster and the number of ¹³C atoms of the labeled precursor. The proposed procedure was able to accurately quantify gravimetrically prepared mixtures of natural and labeled glucose molar fractions of 0.07% and mixtures of natural and labeled glycine at molar fractions down to 0.7%. The method was applied to initial studies of glucose metabolism of different prostate cancer cell lines.

The importance of accurately correcting for the natural abundance of stable isotopes

The use of isotopically labeled tracer substrates is an experimental approach for measuring in vivo and in vitro intracellular metabolic dynamics. Stable isotopes that alter the mass but not the chemical behavior of a molecule are commonly used in isotope tracer studies. Because stable isotopes of some atoms naturally occur at non-negligible abundances, it is important to account for the natural abundance of these isotopes when analyzing data from isotope labeling experiments. Specifically, a distinction must be made between isotopes introduced experimentally via an isotopically labeled tracer and the isotopes naturally present at the start of an experiment. In this tutorial review, we explain the underlying theory of natural abundance correction of stable isotopes, a concept not always understood by metabolic researchers. We also provide a comparison of distinct methods for performing this correction and discuss natural abundance correction in the context of steady state 13C metabolic flux, a method increasingly used to infer intracellular metabolic flux from isotope experiments.

Management of Multiple Nitrogen Sources during Wine Fermentation by Saccharomyces cerevisiae

During fermentative growth in natural and industrial environments, Saccharomyces cerevisiae must redistribute the available nitrogen from multiple exogenous sources to amino acids in order to suitably fulfill anabolic requirements. To exhaustively explore the management of this complex resource, we developed an advanced strategy based on the reconciliation of data from a set of stable isotope tracer experiments with labeled nitrogen sources. Thus, quantifying the partitioning of the N compounds through the metabolism network during fermentation, we demonstrated that, contrary to the generally accepted view, only a limited fraction of most of the consumed amino acids is directly incorporated into proteins. Moreover, substantial catabolism of these molecules allows for efficient redistribution of nitrogen, supporting the operative de novo synthesis of proteinogenic amino acids. In contrast, catabolism of consumed amino acids plays a minor role in the formation of volatile compounds. Another important feature is that the α-keto acid precursors required for the de novo syntheses originate mainly from the catabolism of sugars, with a limited contribution from the anabolism of consumed amino acids. This work provides a comprehensive view of the intracellular fate of consumed nitrogen sources and the metabolic origin of proteinogenic amino acids, highlighting a strategy of distribution of metabolic fluxes implemented by yeast as a means of adapting to environments with changing and scarce nitrogen resources.IMPORTANCE A current challenge for the wine industry, in view of the extensive competition in the worldwide market, is to meet consumer expectations regarding the sensory profile of the product while ensuring an efficient fermentation process. Understanding the intracellular fate of the nitrogen sources available in grape juice is essential to the achievement of these objectives, since nitrogen utilization affects both the fermentative activity of yeasts and the formation of flavor compounds. However, little is known about how the metabolism operates when nitrogen is provided as a composite mixture, as in grape must. Here we quantitatively describe the distribution through the yeast metabolic network of the N moieties and C backbones of these nitrogen sources. Knowledge about the management of a complex resource, which is devoted to improvement of the use of the scarce N nutrient for growth, will be useful for better control of the fermentation process and the sensory quality of wines.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.

Quantification of Microbial Phenotypes

Metabolite profiling technologies have improved to generate close to quantitative metabolomics data, which can be employed to quantitatively describe the metabolic phenotype of an organism. Here, we review the current technologies available for quantitative metabolomics, present their advantages and drawbacks, and the current challenges to generate fully quantitative metabolomics data. Metabolomics data can be integrated into metabolic networks using thermodynamic principles to constrain the directionality of reactions. Here we explain how to estimate Gibbs energy under physiological conditions, including examples of the estimations, and the different methods for thermodynamics-based network analysis. The fundamentals of the methods and how to perform the analyses are described. Finally, an example applying quantitative metabolomics to a yeast model by 13C fluxomics and thermodynamics-based network analysis is presented. The example shows that (1) these two methods are complementary to each other; and (2) there is a need to take into account Gibbs energy errors. Better estimations of metabolic phenotypes will be obtained when further constraints are included in the analysis.

A simplified calculation procedure for mass isotopomer distribution analysis (MIDA) based on multiple linear regression

We have developed a novel, rapid and easy calculation procedure for Mass Isotopomer Distribution Analysis based on multiple linear regression which allows the simultaneous calculation of the precursor pool enrichment and the fraction of newly synthesized labelled proteins (fractional synthesis) using linear algebra. To test this approach, we used the peptide RGGGLK as a model tryptic peptide containing three subunits of glycine. We selected glycine labelled in two ¹³ C atoms (¹³ C₂ -glycine) as labelled amino acid to demonstrate that spectral overlap is not a problem in the proposed methodology. The developed methodology was tested first in vitro by changing the precursor pool enrichment from 10 to 40% of ¹³ C₂ -glycine. Secondly, a simulated in vivo synthesis of proteins was designed by combining the natural abundance RGGGLK peptide and 10 or 20% ¹³ C₂ -glycine at 1 : 1, 1 : 3 and 3 : 1 ratios. Precursor pool enrichments and fractional synthesis values were calculated with satisfactory precision and accuracy using a simple spreadsheet. This novel approach can provide a relatively rapid and easy means to measure protein turnover based on stable isotope tracers. Copyright © 2016 John Wiley & Sons, Ltd.

Identifying Metabolic Subpopulations from Population Level Mass Spectrometry

Metabolism underlies many important cellular decisions, such as the decisions to proliferate and differentiate, and defects in metabolic signaling can lead to disease and aging. In addition, metabolic heterogeneity can have biological consequences, such as differences in outcomes and drug susceptibilities in cancer and antibiotic treatments. Many approaches exist for characterizing the metabolic state of a population of cells, but technologies for measuring metabolism at the single cell level are in the preliminary stages and are limited. Here, we describe novel analysis methodologies that can be applied to established experimental methods to measure metabolic variability within a population. We use mass spectrometry to analyze amino acid composition in cells grown in a mixture of (12)C- and (13)C-labeled sugars; these measurements allow us to quantify the variability in sugar usage and thereby infer information about the behavior of cells within the population. The methodologies described here can be applied to a large range of metabolites and macromolecules and therefore have the potential for broad applications.

13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production

Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many ¹³C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms

Natural isotope correction of MS/MS measurements for metabolomics and (13)C fluxomics

Fluxomics and metabolomics are crucial tools for metabolic engineering and biomedical analysis to determine the in vivo cellular state. Especially, the application of (13)C isotopes allows comprehensive insights into the functional operation of cellular metabolism. Compared to single MS, tandem mass spectrometry (MS/MS) provides more detailed and accurate measurements of the metabolite enrichment patterns (tandem mass isotopomers), increasing the accuracy of metabolite concentration measurements and metabolic flux estimation. MS-type data from isotope labeling experiments is biased by naturally occurring stable isotopes (C, H, N, O, etc.). In particular, GC-MS(/MS) requires derivatization for the usually non-volatile intracellular metabolites introducing additional natural isotopes leading to measurements that do not directly represent the carbon labeling distribution. To make full use of LC- and GC-MS/MS mass isotopomer measurements, the influence of natural isotopes has to be eliminated (corrected). Our correction approach is analyzed for the two most common applications; (13)C fluxomics and isotope dilution mass spectrometry (IDMS) based metabolomics. Natural isotopes can have an impact on the calculated flux distribution which strongly depends on the substrate labeling and the actual flux distribution. Second, we show that in IDMS based metabolomics natural isotopes lead to underestimated concentrations that can and should be corrected with a nonlinear calibration. Our simulations indicate that the correction for natural abundance in isotope based fluxomics and quantitative metabolomics is essential for correct data interpretation.

White biotechnology: State of the art strategies for the development of biocatalysts for biorefining

White biotechnology is a term that is now often used to describe the implementation of biotechnology in the industrial sphere. Biocatalysts (enzymes and microorganisms) are the key tools of white biotechnology, which is considered to be one of the key technological drivers for the growing bioeconomy. Biocatalysts are already present in sectors such as the chemical and agro-food industries, and are used to manufacture products as diverse as antibiotics, paper pulp, bread or advanced polymers. This review proposes an original and global overview of highly complementary fields of biotechnology at both enzyme and microorganism level. A certain number of state of the art approaches that are now being used to improve the industrial fitness of biocatalysts particularly focused on the biorefinery sector are presented. The first part deals with the technologies that underpin the development of industrial biocatalysts, notably the discovery of new enzymes and enzyme improvement using directed evolution techniques. The second part describes the toolbox available by the cell engineer to shape the metabolism of microorganisms. And finally the last part focuses on the 'omic' technologies that are vital for understanding and guide microbial engineering toward more efficient microbial biocatalysts. Altogether, these techniques and strategies will undoubtedly help to achieve the challenging task of developing consolidated bioprocessing (i.e. CBP) readily available for industrial purpose.

Non-targeted Tracer Fate Detection

Stable isotopes have been used to trace atoms through metabolism and quantify metabolic fluxes for several decades. Only recently non-targeted stable isotope labeling approaches have emerged as a powerful tool to gain insights into metabolism. However, the manual detection of isotopic enrichment for a non-targeted analysis is tedious and time consuming. To overcome this limitation, the non-targeted tracer fate detection (NTFD) algorithm for the automated metabolome-wide detection of isotopic enrichment has been developed. NTFD detects and quantifies isotopic enrichment in the form of mass isotopomer distributions (MIDs) in an automated manner, providing the means to trace functional groups, determine MIDs for metabolic flux analysis, or detect tracer-derived molecules in general. Here, we describe the algorithmic background of NTFD, discuss practical considerations for the freely available NTFD software package, and present potential applications of non-targeted stable isotope labeling analysis.

Protein enrichment of an Opuntia ficus-indica cladode hydrolysate by cultivation of Candida utilis and Kluyveromyces marxianus

BACKGROUND

The cladodes of Opuntia ficus-indica (prickly pear cactus) have a low protein content; for use as a balanced feed, supplementation with other protein sources is therefore desirable. We investigated protein enrichment by cultivation of the yeasts Candida utilis and Kluyveromyces marxianus in an enzymatic hydrolysate of the cladode biomass.

RESULTS

Dilute acid pretreatment and enzymatic hydrolysis of sun-dried cladodes resulted in a hydrolysate containing (per litre) 45.5 g glucose, 6.3 g xylose, 9.1 g galactose, 10.8 g arabinose and 9.6 g fructose. Even though K. marxianus had a much higher growth rate and utilized l-arabinose and d-galactose more completely than C. utilis, its biomass yield coefficient was lower due to ethanol and ethyl acetate production despite aerobic cultivation. Yeast cultivation more than doubled the protein content of the hydrolysate, with an essential amino acid profile superior to sorghum and millet grains.

CONCLUSIONS

This K. marxianus strain was weakly Crabtree positive. Despite its low biomass yield, its performance compared well with C. utilis. This is the first report showing that the protein content and quality of O. ficus-indica cladode biomass could substantially be improved by yeast cultivation, including a comparative evaluation of C. utilis and K. marxianus.

Tracking the metabolic pulse of plant lipid production with isotopic labeling and flux analyses: Past, present and future

Metabolism is comprised of networks of chemical transformations, organized into integrated biochemical pathways that are the basis of cellular operation, and function to sustain life. Metabolism, and thus life, is not static. The rate of metabolites transitioning through biochemical pathways (i.e., flux) determines cellular phenotypes, and is constantly changing in response to genetic or environmental perturbations. Each change evokes a response in metabolic pathway flow, and the quantification of fluxes under varied conditions helps to elucidate major and minor routes, and regulatory aspects of metabolism. To measure fluxes requires experimental methods that assess the movements and transformations of metabolites without creating artifacts. Isotopic labeling fills this role and is a long-standing experimental approach to identify pathways and quantify their metabolic relevance in different tissues or under different conditions. The application of labeling techniques to plant science is however far from reaching it potential. In light of advances in genetics and molecular biology that provide a means to alter metabolism, and given recent improvements in instrumentation, computational tools and available isotopes, the use of isotopic labeling to probe metabolism is becoming more and more powerful. We review the principal analytical methods for isotopic labeling with a focus on seminal studies of pathways and fluxes in lipid metabolism and carbon partitioning through central metabolism. Central carbon metabolic steps are directly linked to lipid production by serving to generate the precursors for fatty acid biosynthesis and lipid assembly. Additionally some of the ideas for labeling techniques that may be most applicable for lipid metabolism in the future were originally developed to investigate other aspects of central metabolism. We conclude by describing recent advances that will play an important future role in quantifying flux and metabolic operation in plant tissues.

CBFA: phenotype prediction integrating metabolic models with constraints derived from experimental data

Background

Flux analysis methods lie at the core of Metabolic Engineering (ME), providing methods for phenotype simulation that allow the determination of flux distributions under different conditions. Although many constraint-based modeling software tools have been developed and published, none provides a free user-friendly application that makes available the full portfolio of flux analysis methods.

Results

This work presents Constraint-based Flux Analysis (CBFA), an open-source software application for flux analysis in metabolic models that implements several methods for phenotype prediction, allowing users to define constraints associated with measured fluxes and/or flux ratios, together with environmental conditions (e.g. media) and reaction/gene knockouts. CBFA identifies the set of applicable methods based on the constraints defined from user inputs, encompassing algebraic and constraint-based simulation methods. The integration of CBFA within the OptFlux framework for ME enables the utilization of different model formats and standards and the integration with complementary methods for phenotype simulation and visualization of results.

Conclusions

A general-purpose and flexible application is proposed that is independent of the origin of the constraints defined for a given simulation. The aim is to provide a simple to use software tool focused on the application of several flux prediction methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-014-0123-1) contains supplementary material, which is available to authorized users.

Computational strategies for a system-level understanding of metabolism

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

OpenFLUX2: (13)C-MFA modeling software package adjusted for the comprehensive analysis of single and parallel labeling experiments

BACKGROUND

Steady-state (13)C-based metabolic flux analysis ((13)C-MFA) is the most powerful method available for the quantification of intracellular fluxes. These analyses include concertedly linked experimental and computational stages: (i) assuming the metabolic model and optimizing the experimental design; (ii) feeding the investigated organism using a chosen (13)C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the (13)C-patterns of intracellular metabolites; and (iv) computing flux parameters that minimize the differences between observed and simulated measurements, followed by evaluating flux statistics. In its early stages, (13)C-MFA was performed on the basis of data obtained in a single labeling experiment (SLE) followed by exploiting the developed high-performance computational software. Recently, the advantages of parallel labeling experiments (PLEs), where several LEs are conducted under the conditions differing only by the tracer(s) choice, were demonstrated, particularly with regard to improving flux precision due to the synergy of complementary information. The availability of an open-source software adjusted for PLE-based (13)C-MFA is an important factor for PLE implementation.

RESULTS

The open-source software OpenFLUX, initially developed for the analysis of SLEs, was extended for the computation of PLE data. Using the OpenFLUX2, in silico simulation confirmed that flux precision is improved when (13)C-MFA is implemented by fitting PLE data to the common model compared with SLE-based analysis. Efficient flux resolution could be achieved in the PLE-mediated analysis when the choice of tracer was based on an experimental design computed to minimize the flux variances from different parts of the metabolic network. The analysis provided by OpenFLUX2 mainly includes (i) the optimization of the experimental design, (ii) the computation of the flux parameters from LEs data, (iii) goodness-of-fit testing of the model's adequacy, (iv) drawing conclusions concerning the identifiability of fluxes and construction of a contribution matrix reflecting the relative contribution of the measurement variances to the flux variances, and (v) precise determination of flux confidence intervals using a fine-tunable and convergence-controlled Monte Carlo-based method.

CONCLUSIONS

The developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state (13)C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.

Isotope labeling pattern study of central carbon metabolites using GC/MS

Determination of fluxes by (13)C tracer experiments depends on monitoring the (13)C labeling pattern of metabolites during isotope experiments. In metabolome-based (13)C metabolic flux analysis, liquid chromatography combined with mass spectrometry or tandem mass spectrometry (LC/MS or LC/MS/MS, respectively) has been mainly used as an analytical platform for isotope pattern studies of central carbon metabolites. However, gas chromatography with mass spectrometry (GC/MS) has several advantages over LC/MS, such as high sensitivity, low cost, ease of operation, and availability of mass spectra databases for comparison. In this study, analysis of isotope pattern for central carbon metabolites using GC/MS was demonstrated. First, a proper set of mass ions for central carbon metabolites was selected based on carbon backbone information and structural isomers of mass fragment ions. A total of 34 mass fragment ions was selected and used for the quantification of 25 central carbon metabolites. Then, to quantify isotope fractions, a natural mass isotopomer library for selected mass fragment ions was constructed and subtracted from isotopomer mass spectra data. The results revealed a surprisingly high abundance of partially labeled (13)C intermediates, such as 56.4% of fructose 6-phosphate and 47.6% of dihydroxyacetone phosphate at isotopic steady state, which were generated in the pentose phosphate pathway. Finally, dynamic changes of isotope fragments of central metabolites were monitored with a U-(13)C glucose stimulus response experiment in Kluyveromyces marxianus. With a comprehensive study of isotope patterns of central carbon metabolites using GC/MS, 25 central carbon metabolites and their isotopic fractions were successfully quantified. Dynamic and precise acquisition of isotope pattern can then be used in combination with proper kinetic models to calculate metabolic fluxes.

isoMETLIN: a database for isotope-based metabolomics

The METLIN metabolite database has become one of the most widely used resources in metabolomics for making metabolite identifications. However, METLIN is not designed to identify metabolites that have been isotopically labeled. As a result, unbiasedly tracking the transformation of labeled metabolites with isotope-based metabolomics is a challenge. Here, we introduce a new database, called isoMETLIN (http://isometlin.scripps.edu/), that has been developed specifically to identify metabolites incorporating isotopic labels. isoMETLIN enables users to search all computed isotopologues derived from METLIN on the basis of mass-to-charge values and specified isotopes of interest, such as (13)C or (15)N. Additionally, isoMETLIN contains experimental MS/MS data on hundreds of isotopomers. These data assist in localizing the position of isotopic labels within a metabolite. From these experimental MS/MS isotopomer spectra, precursor atoms can be mapped to fragments. The MS/MS spectra of additional isotopomers can then be computationally generated and included within isoMETLIN. Given that isobaric isotopomers cannot be separated chromatographically or by mass but are likely to occur simultaneously in a biological system, we have also implemented a spectral-mixing function in isoMETLIN. This functionality allows users to combine MS/MS spectra from various isotopomers in different ratios to obtain a theoretical MS/MS spectrum that matches the MS/MS spectrum from a biological sample. Thus, by searching MS and MS/MS experimental data, isoMETLIN facilitates the identification of isotopologues as well as isotopomers from biological samples and provides a platform to drive the next generation of isotope-based metabolomic studies.

A peptide-based method for 13C Metabolic Flux Analysis in microbial communities

The study of intracellular metabolic fluxes and inter-species metabolite exchange for microbial communities is of crucial importance to understand and predict their behaviour. The most authoritative method of measuring intracellular fluxes, ¹³C Metabolic Flux Analysis (¹³C MFA), uses the labeling pattern obtained from metabolites (typically amino acids) during ¹³C labeling experiments to derive intracellular fluxes. However, these metabolite labeling patterns cannot easily be obtained for each of the members of the community. Here we propose a new type of ¹³C MFA that infers fluxes based on peptide labeling, instead of amino acid labeling. The advantage of this method resides in the fact that the peptide sequence can be used to identify the microbial species it originates from and, simultaneously, the peptide labeling can be used to infer intracellular metabolic fluxes. Peptide identity and labeling patterns can be obtained in a high-throughput manner from modern proteomics techniques. We show that, using this method, it is theoretically possible to recover intracellular metabolic fluxes in the same way as through the standard amino acid based ¹³C MFA, and quantify the amount of information lost as a consequence of using peptides instead of amino acids. We show that by using a relatively small number of peptides we can counter this information loss. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

Microbial communities underlie a variety of important biochemical processes ranging from underground cave formation to gold mining or the onset of obesity. Metabolic fluxes describe how carbon and energy flow through the microbial community and therefore provide insights that are rarely captured by other techniques, such as metatranscriptomics or metaproteomics. The most authoritative method to measure fluxes for pure cultures consists of feeding the cells a labeled carbon source and deriving the fluxes from the ensuing metabolite labeling pattern (typically amino acids). Since we cannot easily separate cells of metabolite for each species in a community, this approach is not generally applicable to microbial communities. Here we present a method to derive fluxes from the labeling of peptides, instead of amino acids. This approach has the advantage that peptides can be assigned to each species in a community in a high-throughput fashion through modern proteomic methods. We show that, by using this method, it is theoretically possible to recover the same amount of information as through the standard approach, if enough peptides are used. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

Decoding the dynamics of cellular metabolism and the action of 3-bromopyruvate and 2-deoxyglucose using pulsed stable isotope-resolved metabolomics

BACKGROUND

Cellular metabolism is highly dynamic and continuously adjusts to the physiological program of the cell. The regulation of metabolism appears at all biological levels: (post-) transcriptional, (post-) translational, and allosteric. This regulatory information is expressed in the metabolome, but in a complex manner. To decode such complex information, new methods are needed in order to facilitate dynamic metabolic characterization at high resolution.

RESULTS

Here, we describe pulsed stable isotope-resolved metabolomics (pSIRM) as a tool for the dynamic metabolic characterization of cellular metabolism. We have adapted gas chromatography-coupled mass spectrometric methods for metabolomic profiling and stable isotope-resolved metabolomics. In addition, we have improved robustness and reproducibility and implemented a strategy for the absolute quantification of metabolites.

CONCLUSIONS

By way of examples, we have applied this methodology to characterize central carbon metabolism of a panel of cancer cell lines and to determine the mode of metabolic inhibition of glycolytic inhibitors in times ranging from minutes to hours. Using pSIRM, we observed that 2-deoxyglucose is a metabolic inhibitor, but does not directly act on the glycolytic cascade.

Reliable Metabolic Flux Estimation in Escherichia coli Central Carbon Metabolism Using Intracellular Free Amino Acids

¹³C metabolic flux analysis (MFA) is a tool of metabolic engineering for investigation of in vivo flux distribution. A direct ¹³C enrichment analysis of intracellular free amino acids (FAAs) is expected to reduce time for labeling experiments of the MFA. Measurable FAAs should, however, vary among the MFA experiments since the pool sizes of intracellular free metabolites depend on cellular metabolic conditions. In this study, minimal ¹³C enrichment data of FAAs was investigated to perform the FAAs-based MFA. An examination of a continuous culture of Escherichia coli using ¹³C-labeled glucose showed that the time required to reach an isotopically steady state for FAAs is rather faster than that for conventional method using proteinogenic amino acids (PAAs). Considering 95% confidence intervals, it was found that the metabolic flux distribution estimated using FAAs has a similar reliability to that of the PAAs-based method. The comparative analysis identified glutamate, aspartate, alanine and phenylalanine as the common amino acids observed in E. coli under different culture conditions. The results of MFA also demonstrated that the ¹³C enrichment data of the four amino acids is required for a reliable analysis of the flux distribution.

Global metabolic response of Enterococcus faecalis to oxygen

Oxygen and oxidative stress have become relevant components in clarifying the mechanism that weakens bacterial cells in parallel to the mode of action of bactericidal antibiotics. Given the importance of oxidative stress in the overall defense mechanism of bacteria and their apparent role in the antimicrobial mode of action, it is important to understand how bacteria respond to this stress at a metabolic level. The aim of this study was to determine the impact of oxygen on the metabolism of the facultative anaerobe Enterococcus faecalis using continuous culture, metabolomics, and (13)C enrichment of metabolic intermediates. When E. faecalis was rapidly transitioned from anaerobic to aerobic growth, cellular metabolism was directed toward intracellular glutathione production and glycolysis was upregulated 2-fold, which increased the supply of critical metabolite precursors (e.g., glycine and glutamate) for sulfur metabolism and glutathione biosynthesis as well as reducing power for cellular respiration in the presence of hemin. The ultimate metabolic response of E. faecalis to an aerobic environment was the upregulation of fatty acid metabolism and benzoate degradation, which was linked to important changes in the bacterial membrane composition as evidenced by changes in membrane fatty acid composition and the reduction of membrane-associated demethylmenaquinone. These key metabolic pathways associated with the response of E. faecalis to oxygen may represent potential new targets to increase the susceptibility of this bacterium to bactericidal drugs.

GC-MS-based ¹³C metabolic flux analysis

The in vivo analysis of metabolic fluxes has become a valuable method for the investigation of microorganisms. It turned out especially useful in industrial biotechnology for the prediction of beneficial genetic targets for rational strain optimization. Here, we describe in detail the procedure for the state-of-the-art approach of (13)C metabolic flux analysis comprising steady-state cultivations in a mineral salt medium with (13)C-labeled substrates, GC-MS measurement for labeling analysis, as well as metabolic modeling using the open-source software OpenFlux.

Deregulation of feedback inhibition of phosphoenolpyruvate carboxylase for improved lysine production in Corynebacterium glutamicum

Allosteric regulation of phosphoenolpyruvate carboxylase (PEPC) controls the metabolic flux distribution of anaplerotic pathways. In this study, the feedback inhibition of Corynebacterium glutamicum PEPC was rationally deregulated, and its effect on metabolic flux redistribution was evaluated. Based on rational protein design, six PEPC mutants were designed, and all of them showed significantly reduced sensitivity toward aspartate and malate inhibition. Introducing one of the point mutations (N917G) into the ppc gene, encoding PEPC of the lysine-producing strain C. glutamicum LC298, resulted in ∼37% improved lysine production. In vitro enzyme assays and (13)C-based metabolic flux analysis showed ca. 20 and 30% increases in the PEPC activity and corresponding flux, respectively, in the mutant strain. Higher demand for NADPH in the mutant strain increased the flux toward pentose phosphate pathway, which increased the supply of NADPH for enhanced lysine production. The present study highlights the importance of allosteric regulation on the flux control of central metabolism. The strategy described here can also be implemented to improve other oxaloacetate-derived products.

Analysis of metabolic flux using dynamic labelling and metabolic modelling

Metabolic fluxes and the capacity to modulate them are a crucial component of the ability of the plant cell to react to environmental perturbations. Our ability to quantify them and to attain information concerning the regulatory mechanisms that control them is therefore essential to understand and influence metabolic networks. For all but the simplest of flux measurements labelling methods have proven to be the most informative. Both steady-state and dynamic labelling approaches have been adopted in the study of plant metabolism. Here the conceptual basis of these complementary approaches, as well as their historical application in microbial, mammalian and plant sciences, is reviewed, and an update on technical developments in label distribution analyses is provided. This is supported by illustrative cases studies involving the kinetic modelling of secondary metabolism. One issue that is particularly complex in the analysis of plant fluxes is the extensive compartmentation of the plant cell. This problem is discussed from both theoretical and experimental perspectives, and the current approaches used to address it are assessed. Finally, current limitations and future perspectives of kinetic modelling of plant metabolism are discussed.

Kinetic isotope effects significantly influence intracellular metabolite (13) C labeling patterns and flux determination

Rigorous mathematical modeling of carbon-labeling experiments allows estimation of fluxes through the pathways of central carbon metabolism, yielding powerful information for basic scientific studies as well as for a wide range of applications. However, the mathematical models that have been developed for flux determination from (13) C labeling data have commonly neglected the influence of kinetic isotope effects on the distribution of (13) C label in intracellular metabolites, as these effects have often been assumed to be inconsequential. We have used measurements of the (13) C isotope effects on the pyruvate dehydrogenase enzyme from the literature to model isotopic fractionation at the pyruvate node and quantify the modeling errors expected to result from the assumption that isotope effects are negligible. We show that under some conditions kinetic isotope effects have a significant impact on the (13) C labeling patterns of intracellular metabolites, and the errors associated with neglecting isotope effects in (13) C-metabolic flux analysis models can be comparable in size to measurement errors associated with GC-MS. Thus, kinetic isotope effects must be considered in any rigorous assessment of errors in (13) C labeling data, goodness-of-fit between model and data, confidence intervals of estimated metabolic fluxes, and statistical significance of differences between estimated metabolic flux distributions.

Software LS-MIDA for efficient mass isotopomer distribution analysis in metabolic modelling

BACKGROUND

The knowledge of metabolic pathways and fluxes is important to understand the adaptation of organisms to their biotic and abiotic environment. The specific distribution of stable isotope labelled precursors into metabolic products can be taken as fingerprints of the metabolic events and dynamics through the metabolic networks. An open-source software is required that easily and rapidly calculates from mass spectra of labelled metabolites, derivatives and their fragments global isotope excess and isotopomer distribution.

RESULTS

The open-source software "Least Square Mass Isotopomer Analyzer" (LS-MIDA) is presented that processes experimental mass spectrometry (MS) data on the basis of metabolite information such as the number of atoms in the compound, mass to charge ratio (m/e or m/z) values of the compounds and fragments under study, and the experimental relative MS intensities reflecting the enrichments of isotopomers in 13C- or 15 N-labelled compounds, in comparison to the natural abundances in the unlabelled molecules. The software uses Brauman's least square method of linear regression. As a result, global isotope enrichments of the metabolite or fragment under study and the molar abundances of each isotopomer are obtained and displayed.

CONCLUSIONS

The new software provides an open-source platform that easily and rapidly converts experimental MS patterns of labelled metabolites into isotopomer enrichments that are the basis for subsequent observation-driven analysis of pathways and fluxes, as well as for model-driven metabolic flux calculations.

Highlighting the tricarboxylic acid cycle: liquid and gas chromatography-mass spectrometry analyses of (13)C-labeled organic acids

The tricarboxylic acid (TCA) cycle is involved in the complete oxidation of organic acids to carbon dioxide in aerobic cells. It not only uses the acetyl-CoA derived from glycolysis but also uses breakdown products of proteins, fatty acids, and nucleic acids. Therefore, the TCA cycle involves numerous carbon fluxes through central metabolism to produce reductant power and transfer the generated electrons to the aerobic electron transport system where energy is formed by oxidative phosphorylation. Although the TCA cycle plays a crucial role in aerobic organisms and tissues, the lack of direct isotopic labeling information in its intermediates (organic acids) makes the quantification of its metabolic fluxes rather approximate. This is the major technical gap that this study intended to fill. In this work, we established and validated liquid and gas chromatography-mass spectrometry methods to determine (13)C labeling in organic acids involved in the TCA cycle using scheduled multiple reaction monitoring and single ion monitoring modes, respectively. Labeled samples were generated using maize embryos cultured with [(13)C]glucose or [(13)C]glutamine. Once steady-state labeling was reached, (13)C-labeled organic acids were extracted and purified. When applying our mass spectrometric methods to those extracts, mass isotopomer abundances of seven major organic acids were successfully determined.

13C-based metabolic flux analysis: fundamentals and practice

Isotope-based metabolic flux analysis is one of the emerging technologies applied to system level metabolic phenotype characterization in metabolic engineering. Among the developed approaches, (13)C-based metabolic flux analysis has been established as a standard tool and has been widely applied to quantitative pathway characterization of diverse biological systems. To implement (13)C-based metabolic flux analysis in practice, comprehending the underlying mathematical and computational modeling fundamentals is of importance along with carefully conducted experiments and analytical measurements. Such knowledge is also crucial when designing (13)C-labeling experiments and properly acquiring key data sets essential for in vivo flux analysis implementation. In this regard, the modeling fundamentals of (13)C-labeling systems and analytical data processing are the main topics we will deal with in this chapter. Along with this, the relevant numerical optimization techniques are addressed to help implementation of the entire computational procedures aiming at (13)C-based metabolic flux analysis in vivo.

Parallel labeling experiments and metabolic flux analysis: Past, present and future methodologies

Radioactive and stable isotopes have been applied for decades to elucidate metabolic pathways and quantify carbon flow in cellular systems using mass and isotope balancing approaches. Isotope-labeling experiments can be conducted as a single tracer experiment, or as parallel labeling experiments. In the latter case, several experiments are performed under identical conditions except for the choice of substrate labeling. In this review, we highlight robust approaches for probing metabolism and addressing metabolically related questions though parallel labeling experiments. In the first part, we provide a brief historical perspective on parallel labeling experiments, from the early metabolic studies when radioisotopes were predominant to present-day applications based on stable-isotopes. We also elaborate on important technical and theoretical advances that have facilitated the transition from radioisotopes to stable-isotopes. In the second part of the review, we focus on parallel labeling experiments for (13)C-metabolic flux analysis ((13)C-MFA). Parallel experiments offer several advantages that include: tailoring experiments to resolve specific fluxes with high precision; reducing the length of labeling experiments by introducing multiple entry-points of isotopes; validating biochemical network models; and improving the performance of (13)C-MFA in systems where the number of measurements is limited. We conclude by discussing some challenges facing the use of parallel labeling experiments for (13)C-MFA and highlight the need to address issues related to biological variability, data integration, and rational tracer selection.