Metabolic flux analysis in a nonstationary system: fed-batch fermentation of a high yielding strain of E. coli producing 1,3-propanediol

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.

Modulation of polyamine metabolic flux in adipose tissue alters the accumulation of body fat by affecting glucose homeostasis

The continued rise in obesity despite public education, awareness and policies indicates the need for mechanism-based therapeutic approaches to help control the disease. Our data, in conjunction with other studies, suggest an unexpected role for the polyamine catabolic enzyme spermidine/spermine-N1-acetyltransferase (SSAT) in fat homeostasis. Our previous studies showed that deletion of SSAT greatly exaggerates weight gain and that the transgenic overexpression suppresses weight gain in mice on a high-fat diet. This discovery is substantial but the underlying molecular linkages are only vaguely understood. Here, we used a comprehensive systems biology approach, on white adipose tissue (WAT), to discover that the partition of acetyl-CoA towards polyamine catabolism alters glucose homeostasis and hence, fat accumulation. Comparative proteomics and antibody-based expression studies of WAT in SSAT knockout, wild type and transgenic mice identified nine proteins with an increasing gradient across the genotypes, all of which correlate with acetyl-CoA consumption in polyamine acetylation. Adipose-specific SSAT knockout mice and global SSAT knockout mice on a high-fat diet exhibited similar growth curves and proteomic patterns in their WAT, confirming that attenuated consumption of acetyl-CoA in acetylation of polyamines in adipose tissue drives the obese phenotype of these mice. Analysis of protein expression indicated that the identified changes in the levels of proteins regulating acetyl-CoA consumption occur via the AMP-activated protein kinase pathway. Together, our data suggest that differential expression of SSAT markedly alters acetyl-CoA levels, which in turn trigger a global shift in glucose metabolism in adipose tissue, thus affecting the accumulation of body fat.

Towards better understanding of an industrial cell factory: investigating the feasibility of real-time metabolic flux analysis in Pichia pastoris

Background

Novel analytical tools, which shorten the long and costly development cycles of biopharmaceuticals are essential. Metabolic flux analysis (MFA) shows great promise in improving our understanding of the metabolism of cell factories in bioreactors, but currently only provides information post-process using conventional off-line methods. MFA combined with real time multianalyte process monitoring techniques provides a valuable platform technology allowing real time insights into metabolic responses of cell factories in bioreactors. This could have a major impact in the bioprocessing industry, ultimately improving product consistency, productivity and shortening development cycles.

Results

This is the first investigation using Near Infrared Spectroscopy (NIRS) in situ combined with metabolic flux modelling which is both a significant challenge and considerable extension of these techniques. We investigated the feasibility of our approach using the industrial workhorse Pichia pastoris in a simplified model system. A parental P. pastoris strain (i.e. which does not synthesize recombinant protein) was used to allow definition of distinct metabolic states focusing solely upon the prediction of intracellular fluxes in central carbon metabolism. Extracellular fluxes were determined using off-line conventional reference methods and on-line NIR predictions (calculated by multivariate analysis using the partial least squares algorithm, PLS). The results showed that the PLS-NIRS models for biomass and glycerol were accurate: correlation coefficients, R², above 0.90 and the root mean square error of prediction, RMSEP, of 1.17 and 2.90 g/L, respectively. The analytical quality of the NIR models was demonstrated by direct comparison with the standard error of the laboratory (SEL), which showed that performance of the NIR models was suitable for quantifying biomass and glycerol for calculating extracellular metabolite rates and used as independent inputs for the MFA (RMSEP lower than 1.5 × SEL). Furthermore, the results for the MFA from both datasets passed consistency tests performed for each steady state, showing that the precision of on-line NIRS is equivalent to that obtained by the off-line measurements.

Conclusions

The findings of this study show for the first time the potential of NIRS as an input generating for MFA models, contributing to the optimization of cell factory metabolism in real-time.

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and (13)C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.

Mapping global effects of the anti-sigma factor MucA in Pseudomonas fluorescens SBW25 through genome-scale metabolic modeling

Background

Alginate is an industrially important polysaccharide, currently produced commercially by harvesting of marine brown sea-weeds. The polymer is also synthesized as an exo-polysaccharide by bacteria belonging to the genera Pseudomonas and Azotobacter, and these organisms may represent an alternative alginate source in the future. The current work describes an attempt to rationally develop a biological system tuned for very high levels of alginate production, based on a fundamental understanding of the system through metabolic modeling supported by transcriptomics studies and carefully controlled fermentations.

Results

Alginate biosynthesis in Pseudomonas fluorescens was studied in a genomics perspective, using an alginate over-producing strain carrying a mutation in the anti-sigma factor gene mucA. Cells were cultivated in chemostats under nitrogen limitation on fructose or glycerol as carbon sources, and cell mass, growth rate, sugar uptake, alginate and CO₂ production were monitored. In addition a genome scale metabolic model was constructed and samples were collected for transcriptome analyses. The analyses show that polymer production operates in a close to optimal way with respect to stoichiometric utilization of the carbon source and that the cells increase the uptake of carbon source to compensate for the additional needs following from alginate synthesis. The transcriptome studies show that in the presence of the mucA mutation, the alg operon is upregulated together with genes involved in energy generation, genes on both sides of the succinate node of the TCA cycle and genes encoding ribosomal and other translation-related proteins. Strains expressing a functional MucA protein (no alginate production) synthesize cellular biomass in an inefficient way, apparently due to a cycle that involves oxidation of NADPH without ATP production. The results of this study indicate that the most efficient way of using a mucA mutant as a cell factory for alginate production would be to use non-growing conditions and nitrogen deprivation.

Conclusions

The insights gained in this study should be very useful for a future efficient production of microbial alginates.

Adaptive bi-level programming for optimal gene knockouts for targeted overproduction under phenotypic constraints

BACKGROUND

Optimization procedures to identify gene knockouts for targeted biochemical overproduction have been widely in use in modern metabolic engineering. Flux balance analysis (FBA) framework has provided conceptual simplifications for genome-scale dynamic analysis at steady states. Based on FBA, many current optimization methods for targeted bio-productions have been developed under the maximum cell growth assumption. The optimization problem to derive gene knockout strategies recently has been formulated as a bi-level programming problem in OptKnock for maximum targeted bio-productions with maximum growth rates. However, it has been shown that knockout mutants in fact reach the steady states with the minimization of metabolic adjustment (MOMA) from the corresponding wild-type strains instead of having maximal growth rates after genetic or metabolic intervention. In this work, we propose a new bi-level computational framework--MOMAKnock--which can derive robust knockout strategies under the MOMA flux distribution approximation.

METHODS

In this new bi-level optimization framework, we aim to maximize the production of targeted chemicals by identifying candidate knockout genes or reactions under phenotypic constraints approximated by the MOMA assumption. Hence, the targeted chemical production is the primary objective of MOMAKnock while the MOMA assumption is formulated as the inner problem of constraining the knockout metabolic flux to be as close as possible to the steady-state phenotypes of wide-type strains. As this new inner problem becomes a quadratic programming problem, a novel adaptive piecewise linearization algorithm is developed in this paper to obtain the exact optimal solution to this new bi-level integer quadratic programming problem for MOMAKnock.

RESULTS

Our new MOMAKnock model and the adaptive piecewise linearization solution algorithm are tested with a small E. coli core metabolic network and a large-scale iAF1260 E. coli metabolic network. The derived knockout strategies are compared with those from OptKnock. Our preliminary experimental results show that MOMAKnock can provide improved targeted productions with more robust knockout strategies.

Using multiple tracers for 13C metabolic flux analysis

(13)C-Metabolic flux analysis ((13)C-MFA) is a powerful technique for quantifying intracellular metabolic fluxes in living cells. These in vivo fluxes provide important information on the physiology of cells in culture, which can be used for metabolic engineering purposes and serve as inputs for systems biology modeling. The (13)C-MFA technique consists of several steps: (1) selecting appropriate tracers for a given system of interest, (2) performing isotopic labeling experiments, (3) measuring isotopic labeling distributions in metabolic products, (4) estimating metabolic fluxes using least-squares regression, and (5) evaluating the goodness of fit and computing confidence intervals for estimated fluxes. In this chapter, we provide guidelines for performing (13)C-MFA studies using multiple isotopic tracers, a technique that is especially useful for elucidating fluxes in complex biological systems where multiple carbon sources are present. Here, as an example, we describe key steps and decision points for designing (13)C-MFA studies for microbes grown on mixtures of glucose and xylose. The general concepts described in this chapter are applicable to many other biological systems. For example, the same procedures can be applied to design (13)C-MFA studies in mammalian cells, which are generally grown in complex media containing multiple substrates such as glucose and amino acids.

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.

Isotopically nonstationary 13C metabolic flux analysis

(13)C metabolic flux analysis (MFA) is a powerful approach for quantifying cell physiology based upon a combination of extracellular flux measurements and intracellular isotope labeling measurements. In this chapter, we present the method of isotopically nonstationary (13)C MFA (INST-MFA), which is applicable to systems that are at metabolic steady state, but are sampled during the transient period prior to achieving isotopic steady state following the introduction of a (13)C tracer. We describe protocols for performing the necessary isotope labeling experiments, for quenching and extraction of intracellular metabolites, for mass spectrometry (MS) analysis of metabolite labeling, and for computational flux estimation using INST-MFA. By combining several recently developed experimental and computational techniques, INST-MFA provides an important new platform for mapping carbon fluxes that is especially applicable to animal cell cultures, autotrophic organisms, industrial bioprocesses, high-throughput experiments, and other systems that are not amenable to steady-state (13)C MFA experiments.

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.

Isotopically nonstationary 13C flux analysis of Myc-induced metabolic reprogramming in B-cells

We assessed several methods of (13)C metabolic flux analysis (MFA) and found that isotopically nonstationary MFA achieved maximum flux resolution in cultured P493-6 B-cells, which have been engineered to provide tunable expression of the Myc oncoprotein. Comparison of metabolic flux maps obtained under oncogenic (High) and endogenous (Low) Myc expression levels revealed network-wide reprogramming in response to ectopic Myc expression. High Myc cells relied more heavily on mitochondrial oxidative metabolism than Low Myc cells and globally upregulated their consumption of amino acids relative to glucose. TCA cycle and amphibolic mitochondrial pathways exhibited 2- to 4-fold flux increases in High Myc cells, in contrast to modest increases in glucose uptake and lactate excretion. Because our MFA approach relied exclusively upon isotopic measurements of protein-bound amino acids and RNA-bound ribose, it is readily applicable to more complex tumor models that are not amenable to direct extraction and isotopic analysis of free intracellular metabolites.

Rational design of ¹³C-labeling experiments for metabolic flux analysis in mammalian cells

Background

¹³C-Metabolic flux analysis (¹³C-MFA) is a standard technique to probe cellular metabolism and elucidate in vivo metabolic fluxes. ¹³C-Tracer selection is an important step in conducting ¹³C-MFA, however, current methods are restricted to trial-and-error approaches, which commonly focus on an arbitrary subset of the tracer design space. To systematically probe the complete tracer design space, especially for complex systems such as mammalian cells, there is a pressing need for new rational approaches to identify optimal tracers.

Results

Recently, we introduced a new framework for optimal ¹³C-tracer design based on elementary metabolite units (EMU) decomposition, in which a measured metabolite is decomposed into a linear combination of so-called EMU basis vectors. In this contribution, we applied the EMU method to a realistic network model of mammalian metabolism with lactate as the measured metabolite. The method was used to select optimal tracers for two free fluxes in the system, the oxidative pentose phosphate pathway (oxPPP) flux and anaplerosis by pyruvate carboxylase (PC). Our approach was based on sensitivity analysis of EMU basis vector coefficients with respect to free fluxes. Through efficient grouping of coefficient sensitivities, simple tracer selection rules were derived for high-resolution quantification of the fluxes in the mammalian network model. The approach resulted in a significant reduction of the number of possible tracers and the feasible tracers were evaluated using numerical simulations. Two optimal, novel tracers were identified that have not been previously considered for ¹³C-MFA of mammalian cells, specifically [2,3,4,5,6-¹³C]glucose for elucidating oxPPP flux and [3,4-¹³C]glucose for elucidating PC flux. We demonstrate that ¹³C-glutamine tracers perform poorly in this system in comparison to the optimal glucose tracers.

Conclusions

In this work, we have demonstrated that optimal tracer design does not need to be a pure simulation-based trial-and-error process; rather, rational insights into tracer design can be gained through the application of the EMU basis vector methodology. Using this approach, rational labeling rules can be established a priori to guide the selection of optimal ¹³C-tracers for high-resolution flux elucidation in complex metabolic network models.

Expanding the chemical palate of cells by combining systems biology and metabolic engineering

The field of Metabolic Engineering has recently undergone a transformation that has led to a rapid expansion of the chemical palate of cells. Now, it is conceivable to produce nearly any organic molecule of interest using a cellular host. Significant advances have been made in the production of biofuels, biopolymers and precursors, pharmaceuticals and nutraceuticals, and commodity and specialty chemicals. Much of this rapid expansion in the field has been, in part, due to synergies and advances in the area of systems biology. Specifically, the availability of functional genomics, metabolomics and transcriptomics data has resulted in the potential to produce a wealth of new products, both natural and non-natural, in cellular factories. The sheer amount and diversity of this data however, means that uncovering and unlocking novel chemistries and insights is a non-obvious exercise. To address this issue, a number of computational tools and experimental approaches have been developed to help expedite the design process to create new cellular factories. This review will highlight many of the systems biology enabling technologies that have reduced the design cycle for engineered hosts, highlight major advances in the expanded diversity of products that can be synthesized, and conclude with future prospects in the field of metabolic engineering.

LC-MS data processing with MAVEN: a metabolomic analysis and visualization engine

MAVEN is an open-source software program for interactive processing of LC-MS-based metabolomics data. MAVEN enables rapid and reliable metabolite quantitation from multiple reaction monitoring data or high-resolution full-scan mass spectrometry data. It automatically detects and reports peak intensities for isotope-labeled metabolites. Menu-driven, click-based navigation allows visualization of raw and analyzed data. Here we provide a User Guide for MAVEN. Step-by-step instructions are provided for data import, peak alignment across samples, identification of metabolites that differ strongly between biological conditions, quantitation and visualization of isotope-labeling patterns, and export of tables of metabolite-specific peak intensities. Together, these instructions describe a workflow that allows efficient processing of raw LC-MS data into a form ready for biological analysis.

Metabolic cartography: experimental quantification of metabolic fluxes from isotopic labelling studies

For the past decade, flux maps have provided researchers with an in-depth perspective on plant metabolism. As a rapidly developing field, significant headway has been made recently in computation, experimentation, and overall understanding of metabolic flux analysis. These advances are particularly applicable to the study of plant metabolism. New dynamic computational methods such as non-stationary metabolic flux analysis are finding their place in the toolbox of metabolic engineering, allowing more organisms to be studied and decreasing the time necessary for experimentation, thereby opening new avenues by which to explore the vast diversity of plant metabolism. Also, improved methods of metabolite detection and measurement have been developed, enabling increasingly greater resolution of flux measurements and the analysis of a greater number of the multitude of plant metabolic pathways. Methods to deconvolute organelle-specific metabolism are employed with increasing effectiveness, elucidating the compartmental specificity inherent in plant metabolism. Advances in metabolite measurements have also enabled new types of experiments, such as the calculation of metabolic fluxes based on (13)CO(2) dynamic labelling data, and will continue to direct plant metabolic engineering. Newly calculated metabolic flux maps reveal surprising and useful information about plant metabolism, guiding future genetic engineering of crops to higher yields. Due to the significant level of complexity in plants, these methods in combination with other systems biology measurements are necessary to guide plant metabolic engineering in the future.

Integrating flux balance analysis into kinetic models to decipher the dynamic metabolism of Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 sequentially utilizes lactate and its waste products (pyruvate and acetate) during batch culture. To decipher MR-1 metabolism, we integrated genome-scale flux balance analysis (FBA) into a multiple-substrate Monod model to perform the dynamic flux balance analysis (dFBA). The dFBA employed a static optimization approach (SOA) by dividing the batch time into small intervals (i.e., ∼400 mini-FBAs), then the Monod model provided time-dependent inflow/outflow fluxes to constrain the mini-FBAs to profile the pseudo-steady-state fluxes in each time interval. The mini-FBAs used a dual-objective function (a weighted combination of "maximizing growth rate" and "minimizing overall flux") to capture trade-offs between optimal growth and minimal enzyme usage. By fitting the experimental data, a bi-level optimization of dFBA revealed that the optimal weight in the dual-objective function was time-dependent: the objective function was constant in the early growth stage, while the functional weight of minimal enzyme usage increased significantly when lactate became scarce. The dFBA profiled biologically meaningful dynamic MR-1 metabolisms: 1. the oxidative TCA cycle fluxes increased initially and then decreased in the late growth stage; 2. fluxes in the pentose phosphate pathway and gluconeogenesis were stable in the exponential growth period; and 3. the glyoxylate shunt was up-regulated when acetate became the main carbon source for MR-1 growth.

Predicting outcomes of steady-state ¹³C isotope tracing experiments using Monte Carlo sampling

BACKGROUND

Carbon-13 (13C) analysis is a commonly used method for estimating reaction rates in biochemical networks. The choice of carbon labeling pattern is an important consideration when designing these experiments. We present a novel Monte Carlo algorithm for finding the optimal substrate input label for a particular experimental objective (flux or flux ratio). Unlike previous work, this method does not require assumption of the flux distribution beforehand.

RESULTS

Using a large E. coli isotopomer model, different commercially available substrate labeling patterns were tested computationally for their ability to determine reaction fluxes. The choice of optimal labeled substrate was found to be dependent upon the desired experimental objective. Many commercially available labels are predicted to be outperformed by complex labeling patterns. Based on Monte Carlo Sampling, the dimensionality of experimental data was found to be considerably less than anticipated, suggesting that effectiveness of 13C experiments for determining reaction fluxes across a large-scale metabolic network is less than previously believed.

CONCLUSIONS

While 13C analysis is a useful tool in systems biology, high redundancy in measurements limits the information that can be obtained from each experiment. It is however possible to compute potential limitations before an experiment is run and predict whether, and to what degree, the rate of each reaction can be resolved.

Selection of tracers for 13C-metabolic flux analysis using elementary metabolite units (EMU) basis vector methodology

Metabolic flux analysis (MFA) is a powerful technique for elucidating in vivo fluxes in microbial and mammalian systems. A key step in (13)C-MFA is the selection of an appropriate isotopic tracer to observe fluxes in a proposed network model. Despite the importance of MFA in metabolic engineering and beyond, current approaches for tracer experiment design are still largely based on trial-and-error. The lack of a rational methodology for selecting isotopic tracers prevents MFA from achieving its full potential. Here, we introduce a new technique for tracer experiment design based on the concept of elementary metabolite unit (EMU) basis vectors. We demonstrate that any metabolite in a network model can be expressed as a linear combination of so-called EMU basis vectors, where the corresponding coefficients indicate the fractional contribution of the EMU basis vector to the product metabolite. The strength of this approach is the decoupling of substrate labeling, i.e. the EMU basis vectors, from the dependence on free fluxes, i.e. the coefficients. In this work, we demonstrate that flux observability inherently depends on the number of independent EMU basis vectors and the sensitivities of coefficients with respect to free fluxes. Specifically, the number of independent EMU basis vectors places hard limits on how many free fluxes can be determined in a model. This constraint is used as a guide for selecting feasible substrate labeling. In three example models, we demonstrate that by maximizing the number of independent EMU basis vectors the observability of a system is improved. Inspection of sensitivities of coefficients with respect to free fluxes provides additional constraints for proper selection of tracers. The present contribution provides a fresh perspective on an important topic in metabolic engineering, and gives practical guidelines and design principles for a priori selection of isotopic tracers for (13)C-MFA studies.

Optimization of 13C isotopic tracers for metabolic flux analysis in mammalian cells

Mammalian cells consume and metabolize various substrates from their surroundings for energy generation and biomass synthesis. Glucose and glutamine, in particular, are the primary carbon sources for proliferating cancer cells. While this combination of substrates generates static labeling patterns for use in (13)C metabolic flux analysis (MFA), the inability of single tracers to effectively label all pathways poses an obstacle for comprehensive flux determination within a given experiment. To address this issue we applied a genetic algorithm to optimize mixtures of (13)C-labeled glucose and glutamine for use in MFA. We identified tracer combinations that minimized confidence intervals in an experimentally determined flux network describing central carbon metabolism in tumor cells. Additional simulations were used to determine the robustness of the [1,2-(13)C(2)]glucose/[U-(13)C(5)]glutamine tracer combination with respect to perturbations in the network. Finally, we experimentally validated the improved performance of this tracer set relative to glucose tracers alone in a cancer cell line. This versatile method allows researchers to determine the optimal tracer combination to use for a specific metabolic network, and our findings applied to cancer cells significantly enhance the ability of MFA experiments to precisely quantify fluxes in higher organisms.

Collisional fragmentation of central carbon metabolites in LC-MS/MS increases precision of ¹³C metabolic flux analysis

Experimental determination of fluxes by (13)C-tracers relies on detection of (13)C-patterns in metabolites or by-products. In the field of (13)C metabolic flux analysis, the most recent developments point toward recording labeling patterns by liquid chromatography (LC)-mass spectrometry (MS)/MS directly in intermediates in central carbon metabolism (CCM) to increase temporal resolution. Surprisingly, the flux studies published so far with LC-MS measurements were based on intact metabolic intermediates-thus neglected the potential benefits of using positional information to improve flux estimation. For the first time, we exploit collisional fragmentation to obtain more fine-grained (13)C-data on intermediates of CCM and investigate their impact in (13)C metabolic flux analysis. For the case study of Bacillus subtilis grown in mineral medium with (13)C-labeled glucose, we compare the flux estimates obtained by iterative isotopologue balancing of (13)C-data obtained either by LC-MS/MS for solely intact intermediates or LC-MS/MS for intact and fragmented intermediates of CCM. We show that with LC-MS/MS data, fragment information leads to more precise estimates of fluxes in pentose phosphate pathway, glycolysis, and to the tricarboxylic acid cycle. Additionally, we present an efficient analytical strategy to rapidly acquire large sets of (13)C-patterns by tandem MS, and an in-depth analysis of the collisional fragmentation of primary intermediates. In the future, this catalogue will enable comprehensive in silico calculability analyses to identify the most sensitive measurements and direct experimental design.

Dynamic metabolic flux analysis (DMFA): a framework for determining fluxes at metabolic non-steady state

Metabolic flux analysis (MFA) is a key tool for measuring in vivo metabolic fluxes in systems at metabolic steady state. Here, we present a new method for dynamic metabolic flux analysis (DMFA) of systems that are not at metabolic steady state. The advantages of our DMFA method are: (1) time-series of metabolite concentration data can be applied directly for estimating dynamic fluxes, making data smoothing and estimation of average extracellular rates unnecessary; (2) flux estimation is achieved without integration of ODEs, or iterations; (3) characteristic metabolic phases in the fermentation data are identified automatically by the algorithm, rather than selected manually/arbitrarily. We demonstrate the application of the new DMFA framework in three example systems. First, we evaluated the performance of DMFA in a simple three-reaction model in terms of accuracy, precision and flux observability. Next, we analyzed a commercial glucose-limited fed-batch process for 1,3-propanediol production. The DMFA method accurately captured the dynamic behavior of the fed-batch fermentation and identified characteristic metabolic phases. Lastly, we demonstrate that DMFA can be used without any assumed metabolic network model for data reconciliation and detection of gross measurement errors using carbon and electron balances as constraints.

Improving protein secretion of a transglutaminase-secreting Corynebacterium glutamicum recombinant strain on the basis of 13C metabolic flux analysis

Corynebacterium glutamicum is known as a host species for amino acid production. This microorganism was recently noticed as a host that produces secreted proteins. In this study, we performed (13)C metabolic flux analysis ((13)C-MFA) on a recombinant C. glutamicum strain that secretes a heterologous transglutaminase (TGase) to improve TGase secretion. For the (13)C-MFA of a TGase-secreting C. glutamicum strain in batch cultivation, a (13)C-labeling experiment and measurement of mass isotopomer distributions of proteinogenic amino acids were performed, and metabolic fluxes were determined considering the changes in fractional (13)C-labeling of proteinogenic amino acids with respect to culture time. The TGase yield increased at the stationary phase but decreased toward its end. The results of (13)C-MFA revealed that the flux from glycolysis to the TCA cycle gradually increased during TGase secretion. We speculate that the NADH/NAD(+) ratio in the cells increases and that as a result, the specific glucose uptake rate decreases in the stationary phase because of the increased flux of the TCA cycle. Since it is expected that a decrease in the NADH/NAD(+) ratio would improve the TGase yield, we tried to enhance lactate production in a TGase-secreting C. glutamicum strain to decrease cellular NADH levels by increasing the pH level in the culture. The TGase yield increased in 1.4-fold by increasing the pH from 6.7 to 7.2, indicating that the TGase yield was successfully improved on the basis of the (13)C-MFA.

Eukaryotic metabolism: measuring compartment fluxes

Metabolic compartmentation represents a major characteristic of eukaryotic cells. The analysis of compartmented metabolic networks is complicated by separation and parallelization of pathways, intracellular transport, and the need for regulatory systems to mediate communication between interdependent compartments. Metabolic flux analysis (MFA) has the potential to reveal compartmented metabolic events, although it is a challenging task requiring demanding experimental techniques and sophisticated modeling. At present no ready-made solution can be provided to cope with the complexity of compartmented metabolic networks, but new powerful tools are emerging. This review gives an overview of different strategies to approach this issue, focusing on different MFA methods and highlighting the additional information that should be included to improve the outcome of an experiment and associate estimation procedures.

A novel deconvolution method for modeling UDP-N-acetyl-D-glucosamine biosynthetic pathways based on (13)C mass isotopologue profiles under non-steady-state conditions

Background

Stable isotope tracing is a powerful technique for following the fate of individual atoms through metabolic pathways. Measuring isotopic enrichment in metabolites provides quantitative insights into the biosynthetic network and enables flux analysis as a function of external perturbations. NMR and mass spectrometry are the techniques of choice for global profiling of stable isotope labeling patterns in cellular metabolites. However, meaningful biochemical interpretation of the labeling data requires both quantitative analysis and complex modeling. Here, we demonstrate a novel approach that involved acquiring and modeling the timecourses of ¹³C isotopologue data for UDP-N-acetyl-D-glucosamine (UDP-GlcNAc) synthesized from [U-¹³C]-glucose in human prostate cancer LnCaP-LN3 cells. UDP-GlcNAc is an activated building block for protein glycosylation, which is an important regulatory mechanism in the development of many prominent human diseases including cancer and diabetes.

Results

We utilized a stable isotope resolved metabolomics (SIRM) approach to determine the timecourse of ¹³C incorporation from [U-¹³C]-glucose into UDP-GlcNAc in LnCaP-LN3 cells. ¹³C Positional isotopomers and isotopologues of UDP-GlcNAc were determined by high resolution NMR and Fourier transform-ion cyclotron resonance-mass spectrometry. A novel simulated annealing/genetic algorithm, called 'Genetic Algorithm for Isotopologues in Metabolic Systems' (GAIMS) was developed to find the optimal solutions to a set of simultaneous equations that represent the isotopologue compositions, which is a mixture of isotopomer species. The best model was selected based on information theory. The output comprises the timecourse of the individual labeled species, which was deconvoluted into labeled metabolic units, namely glucose, ribose, acetyl and uracil. The performance of the algorithm was demonstrated by validating the computed fractional ¹³C enrichment in these subunits against experimental data. The reproducibility and robustness of the deconvolution were verified by replicate experiments, extensive statistical analyses, and cross-validation against NMR data.

Conclusions

This computational approach revealed the relative fluxes through the different biosynthetic pathways of UDP-GlcNAc, which comprises simultaneous sequential and parallel reactions, providing new insight into the regulation of UDP-GlcNAc levels and O-linked protein glycosylation. This is the first such analysis of UDP-GlcNAc dynamics, and the approach is generally applicable to other complex metabolites comprising distinct metabolic subunits, where sufficient numbers of isotopologues can be unambiguously resolved and accurately measured.

Construction of an E. Coli genome-scale atom mapping model for MFA calculations

Metabolic flux analysis (MFA) has so far been restricted to lumped networks lacking many important pathways, partly due to the difficulty in automatically generating isotope mapping matrices for genome-scale metabolic networks. Here we introduce a procedure that uses a compound matching algorithm based on the graph theoretical concept of pattern recognition along with relevant reaction information to automatically generate genome-scale atom mappings which trace the path of atoms from reactants to products for every reaction. The procedure is applied to the iAF1260 metabolic reconstruction of Escherichia coli yielding the genome-scale isotope mapping model imPR90068. This model maps 90,068 non-hydrogen atoms that span all 2,077 reactions present in iAF1260 (previous largest mapping model included 238 reactions). The expanded scope of the isotope mapping model allows the complete tracking of labeled atoms through pathways such as cofactor and prosthetic group biosynthesis and histidine metabolism. An EMU representation of imPR90068 is also constructed and made available.

Bridging the gap between fluxomics and industrial biotechnology

Metabolic flux analysis is a vital tool used to determine the ultimate output of cellular metabolism and thus detect biotechnologically relevant bottlenecks in productivity. ¹³C-based metabolic flux analysis (¹³C-MFA) and flux balance analysis (FBA) have many potential applications in biotechnology. However, noteworthy hurdles in fluxomics study are still present. First, several technical difficulties in both ¹³C-MFA and FBA severely limit the scope of fluxomics findings and the applicability of obtained metabolic information. Second, the complexity of metabolic regulation poses a great challenge for precise prediction and analysis of metabolic networks, as there are gaps between fluxomics results and other omics studies. Third, despite identified metabolic bottlenecks or sources of host stress from product synthesis, it remains difficult to overcome inherent metabolic robustness or to efficiently import and express nonnative pathways. Fourth, product yields often decrease as the number of enzymatic steps increases. Such decrease in yield may not be caused by rate-limiting enzymes, but rather is accumulated through each enzymatic reaction. Fifth, a high-throughput fluxomics tool hasnot been developed for characterizing nonmodel microorganisms and maximizing their application in industrial biotechnology. Refining fluxomics tools and understanding these obstacles will improve our ability to engineer highlyefficient metabolic pathways in microbial hosts.

13C-metabolic flux analysis for batch culture of Escherichia coli and its Pyk and Pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites

Since most bio-production processes are conducted in a batch or fed-batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied (13)C-metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild-type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild-type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose-6-phosphate dehydrogenase reaction was the rate-limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild-type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation.

Synergy between (13)C-metabolic flux analysis and flux balance analysis for understanding metabolic adaptation to anaerobiosis in E. coli

Genome-based Flux Balance Analysis (FBA) and steady-state isotopic-labeling-based Metabolic Flux Analysis (MFA) are complimentary approaches to predicting and measuring the operation and regulation of metabolic networks. Here, genome-derived models of Escherichia coli (E. coli) metabolism were used for FBA and ¹³C-MFA analyses of aerobic and anaerobic growths of wild-type E. coli (K-12 MG1655) cells. Validated MFA flux maps reveal that the fraction of maintenance ATP consumption in total ATP production is about 14% higher under anaerobic (51.1%) than aerobic conditions (37.2%). FBA revealed that an increased ATP utilization is consumed by ATP synthase to secrete protons from fermentation. The TCA cycle is shown to be incomplete in aerobically growing cells and submaximal growth is due to limited oxidative phosphorylation. An FBA was successful in predicting product secretion rates in aerobic culture if both glucose and oxygen uptake measurement were constrained, but the most-frequently predicted values of internal fluxes yielded from sampling the feasible space differ substantially from MFA-derived fluxes.

Resolving the TCA cycle and pentose-phosphate pathway of Clostridium acetobutylicum ATCC 824: Isotopomer analysis, in vitro activities and expression analysis

Solventogenic clostridia are an important class of microorganisms that can produce various biofuels. One of the bottlenecks in engineering clostridia stems from the fact that central metabolic pathways remain poorly understood. Here, we utilized the power of (13) C-based isotopomer analysis to re-examine central metabolic pathways of Clostridium acetobutylicum ATCC 824. We demonstrate using [1,2-(13) C]glucose, MS analysis of intracellular metabolites, and enzymatic assays that C. acetobutylicum has a split TCA cycle where only Re-citrate synthase (CS) contributes to the production of α-ketoglutarate via citrate. Furthermore, we show that there is no carbon exchange between α-ketoglutarate and fumarate and that the oxidative pentose-phosphate pathway (oxPPP) is inactive. Dynamic gene expression analysis of the putative Re-CS gene (CAC0970), its operon, and all glycolysis, pentose-phosphate pathway, and TCA cycle genes identify genes and their degree of involvement in these core pathways that support the powerful primary metabolism of this industrial organism.

Metabolic fluxes and beyond-systems biology understanding and engineering of microbial metabolism

The recent years have seen tremendous progress towards the understanding of microbial metabolism on a higher level of the entire functional system. Hereby, huge achievements including the sequencing of complete genomes and efficient post-genomic approaches provide the basis for a new, fascinating era of research-analysis of metabolic and regulatory properties on a global scale. Metabolic flux (fluxome) analysis displays the first systems oriented approach to unravel the physiology of microorganisms since it combines experimental data with metabolic network models and allows determining absolute fluxes through larger networks of central carbon metabolism. Hereby, fluxes are of central importance for systems level understanding because they fundamentally represent the cellular phenotype as integrated output of the cellular components, i.e. genes, transcripts, proteins, and metabolites. A currently emerging and promising area of research in systems biology and systems metabolic engineering is therefore the integration of fluxome data in multi-omics studies to unravel the multiple layers of control that superimpose the flux network and enable its optimal operation under different environmental conditions.

Soft constraints-based multiobjective framework for flux balance analysis

The current state of the art for linear optimization in Flux Balance Analysis has been limited to single objective functions. Since mammalian systems perform various functions, a multiobjective approach is needed when seeking optimal flux distributions in these systems. In most of the available multiobjective optimization methods, there is a lack of understanding of when to use a particular objective, and how to combine and/or prioritize mutually competing objectives to achieve a truly optimal solution. To address these limitations we developed a soft constraints based linear physical programming-based flux balance analysis (LPPFBA) framework to obtain a multiobjective optimal solutions. The developed framework was first applied to compute a set of multiobjective optimal solutions for various pairs of objectives relevant to hepatocyte function (urea secretion, albumin, NADPH, and glutathione syntheses) in bioartificial liver systems. Next, simultaneous analysis of the optimal solutions for three objectives was carried out. Further, this framework was utilized to obtain true optimal conditions to improve the hepatic functions in a simulated bioartificial liver system. The combined quantitative and visualization framework of LPPFBA is applicable to any large-scale metabolic network system, including those derived by genomic analyses.

Metabolic flux analysis for recombinant protein production by Pichia pastoris using dual carbon sources: Effects of methanol feeding rate

The intracellular metabolic fluxes through the central carbon pathways in the bioprocess for recombinant human erythropoietin (rHuEPO) production by Pichia pastoris (Mut(+)) were calculated to investigate the metabolic effects of dual carbon sources (methanol/sorbitol) and the methanol feed rate, and to obtain a deeper understanding of the regulatory circuitry of P. pastoris, using the established stoichiometry-based model containing 102 metabolites and 141 reaction fluxes. Four fed-batch operations with (MS-) and without (M-) sorbitol were performed at three different constant specific growth rates (h(-1)), and denoted as M-0.03, MS-0.02, MS-0.03, and MS-0.04. Considering the methanol consumption pathway, the M-0.03 and MS-0.02 conditions produced similar effects and had >85% of formaldehyde flux towards the assimilatory pathway. In contrast, the use of the dual carbon source condition generated a shift in metabolism towards the dissimilatory pathway that corresponded to the shift in dilution rate from MS-0.03 to MS-0.04, indicating that the methanol feed exceeded the metabolic requirements at the higher micro(0). Comparing M-0.03 and MS-0.03 conditions, which had the same methanol feeding rates, sorbitol addition increased the rHuEPO synthetic flux 4.4-fold. The glycolysis, gluconeogenesis, and PPP pathways worked uninterruptedly only at MS-0.02 condition. PPP and TCA cycles worked with the highest disturbances at MS-0.04 condition, which shows the stress of increased feeding rates of methanol on cell metabolism.

Mass spectrometry-based metabolomics of yeast

Driven by the advent of metabolomics, recent years have seen renewed interest in the investigation of yeast metabolism. Here we provide a practical guide to metabolomic analysis of yeast using liquid chromatography-mass spectrometry (LC-MS). We begin with background on LC-MS and its utility in studying yeast metabolism. We then describe key issues involved at each step of a typical yeast metabolomics experiment: in experimental design, cell culture, metabolite extraction, LC-MS, and data processing and analysis. Throughout, we highlight interdependencies between the steps that are relevant to developing an integrated workflow which effectively leverages LC-MS to reveal yeast biology.

Industrial systems biology

The chemical industry is currently undergoing a dramatic change driven by demand for developing more sustainable processes for the production of fuels, chemicals, and materials. In biotechnological processes different microorganisms can be exploited, and the large diversity of metabolic reactions represents a rich repository for the design of chemical conversion processes that lead to efficient production of desirable products. However, often microorganisms that produce a desirable product, either naturally or because they have been engineered through insertion of heterologous pathways, have low yields and productivities, and in order to establish an economically viable process it is necessary to improve the performance of the microorganism. Here metabolic engineering is the enabling technology. Through metabolic engineering the metabolic landscape of the microorganism is engineered such that there is an efficient conversion of the raw material, typically glucose, to the product of interest. This process may involve both insertion of new enzymes activities, deletion of existing enzyme activities, but often also deregulation of existing regulatory structures operating in the cell. In order to rapidly identify the optimal metabolic engineering strategy the industry is to an increasing extent looking into the use of tools from systems biology. This involves both x-ome technologies such as transcriptome, proteome, metabolome, and fluxome analysis, and advanced mathematical modeling tools such as genome-scale metabolic modeling. Here we look into the history of these different techniques and review how they find application in industrial biotechnology, which will lead to what we here define as industrial systems biology.