Systems-level analysis of genome-scalein silico metabolic models using MetaFluxNet

The systems-level analysis of microbes with myriad of heterologous data generated by omics technologies has been applied to improve our understanding of cellular function and physiology and consequently to enhance production of various bioproducts. At the heart of this revolution residesin silico genome-scale metabolic model. In order to fully exploit the power of genome-scale model, a systematic approach employing user-friendly software is required. Metabolic flux analysis of genome-scale metabolic network is becoming widely employed to quantify the flux distribution and validate model-driven hypotheses. Here we describe the development of an upgraded MetaFluxNet which allows (1) construction of metabolic models connected to metabolic databases, (2) calculation of fluxes by metabolic flux analysis, (3) comparative flux analysis with flux-profile visualization, (4) the use of metabolic flux analysis markup language to enable models to be exchanged efficiently, and (5) the exporting of data from constraints-based flux analysis into various formats. MetaFluxNet also allows cellular physiology to be predicted and strategies for strain improvement to be developed from genome-based information on flux distributions. This integrated software environment promises to enhance our understanding on metabolic network at a whole organism level and to establish novel strategies for improving the properties of organisms for various biotechnological applications.     

Strategies for investigating the plant metabolic network with steady-state metabolic flux analysis: lessons from an Arabidopsis cell culture and other systems

Steady-state (13)C metabolic flux analysis (MFA) is currently the experimental method of choice for generating flux maps of the compartmented network of primary metabolism in heterotrophic and mixotrophic plant tissues. While statistically robust protocols for the application of steady-state MFA to plant tissues have been developed by several research groups, the implementation of the method is still far from routine. The effort required to produce a flux map is more than justified by the information that it contains about the metabolic phenotype of the system, but it remains the case that steady-state MFA is both analytically and computationally demanding. This article provides an overview of principles that underpin the implementation of steady-state MFA, focusing on the definition of the metabolic network responsible for redistribution of the label, experimental considerations relating to data collection, the modelling process that allows a set of metabolic fluxes to be deduced from the labelling data, and the interpretation of flux maps. The article draws on published studies of Arabidopsis cell cultures and other systems, including developing oilseeds, with the aim of providing practical guidance and strategies for handling the issues that arise when applying steady-state MFA to the complex metabolic networks encountered in plants.     

Insights into plant metabolic networks from steady-state metabolic flux analysis

Steady-state metabolic flux analysis (MFA) is an experimental approach that allows the measurement of multiple fluxes in the core network of primary carbon metabolism. It is based on isotopic labelling experiments, and although well established in the analysis of micro-organisms, and some mammalian systems, the extension of the method to plant cells has been challenging because of the extensive subcellular compartmentation of the metabolic network. Despite this difficulty there has been substantial progress in developing robust protocols for the analysis of heterotrophic plant metabolism by steady-state MFA, and flux maps have now been published that reflect the metabolic phenotypes of excised root tips, developing embryos and cotyledons, hairy root cultures, and cell suspensions under a variety of physiological conditions. There has been a steady improvement in the quality, extent and statistical reliability of these analyses, and new information is emerging on the performance of the plant metabolic network and the contributions of specific pathways.     

Tandem mass spectrometry: a novel approach for metabolic flux analysis

The goal of metabolic flux analysis (MFA) is the accurate estimation of intracellular fluxes in metabolic networks. Here, we introduce a new method for MFA based on tandem mass spectrometry (MS) and stable-isotope tracer experiments. We demonstrate that tandem MS provides more labeling information than can be obtained from traditional full scan MS analysis and allows estimation of fluxes with better precision. We present a modeling framework that takes full advantage of the additional labeling information obtained from tandem MS for MFA. We show that tandem MS data can be computed for any network model, any compound and any tandem MS fragmentation using linear mapping of isotopomers. The inherent advantages of tandem MS were illustrated in two network models using simulated and literature data. Application of tandem MS increased the observability of the models and improved the precision of estimated fluxes by 2- to 5-fold compared to traditional MS analysis.     

A small metabolic flux model to identify transient metabolic regulations in Saccharomyces cerevisiae

The understanding of dynamic metabolic regulations is important for physiological studies and strain characterization tasks. The present study combined transient experiments with online metabolic flux analysis (MFA) in order to quantify metabolic regulations, namely carbon catabolite repression of respiration and transient acetic-acid production, in Saccharomyces cerevisiae during aerobic growth on glucose. The aim was to investigate which additional information can be gained from using a small metabolic flux model to study transient growth provoked by shift-up and shift-down experiments, compared to online monitoring alone. The MFA model allowed us to propose new correlations between pathways of the central metabolism. A linear correlation between glycolytic flux and respiratory capacity holds for shift-down and shift-up experiments. This confirmed that respiratory functions were subjected to carbon catabolite repression and suggested that respiratory capacity is controlled by the glycolytic flux rather than the glucose influx. Furthermore, the model showed that control of repression of respiration by the glycolytic flux was a dynamic phenomenon. Co-factor balancing within the MFA model showed that transient acetic-acid production indicated a transient limitation in another part of the central metabolism but not in oxidative phosphorylation. However, at super-critical growth rates and when coupling of anabolism and catabolism is resumed, the limitation shifts to oxidative phosphorylation, with the consequence that ethanol is formed. The online application of small metabolic flux models to transient experiments enhanced the physiological insight into transient growth and opens up the use of transient experiments as an efficient tool to understand dynamic metabolic regulations.     

Hybrid dynamic/static method for large-scale simulation of metabolism

Background  

Many computer studies have employed either dynamic simulation or metabolic flux analysis (MFA) to predict the behaviour of biochemical pathways. Dynamic simulation determines the time evolution of pathway properties in response to environmental changes, whereas MFA provides only a snapshot of pathway properties within a particular set of environmental conditions. However, owing to the large amount of kinetic data required for dynamic simulation, MFA, which requires less information, has been used to manipulate large-scale pathways to determine metabolic outcomes.     

The thermodynamic meaning of metabolic exchange fluxes

Metabolic flux analysis (MFA) deals with the experimental determination of steady-state fluxes in metabolic networks. An important feature of the ¹³C MFA method is its capability to generate information on both directions of bidirectional reaction steps given by exchange fluxes. The biological interpretation of these exchange fluxes and their relation to thermodynamic properties of the respective reaction steps has never been systematically investigated. As a central result, it is shown here that for a general class of enzyme reaction mechanisms the quotients of net and exchange fluxes measured by ¹³C MFA are coupled to Gibbs energies of the reaction steps. To establish this relation the concept of apparent flux ratios of enzymatic isotope-labeling networks is introduced and some computing rules for these flux ratios are given. Application of these rules reveals a conceptional pitfall of ¹³C MFA, which is the inherent dependency of measured exchange fluxes on the chosen tracer atom. However, it is shown that this effect can be neglected for typical biochemical reaction steps under physiological conditions. In this situation, the central result can be formulated as a two-sided inequality relating fluxes, pool sizes, and standard Gibbs energies. This relation has far-reaching consequences for metabolic flux analysis, quantitative metabolomics, and network thermodynamics.     

Towards high throughput metabolic flux analysis in plants

Research on plant metabolism is currently experiencing the common use of various omics methods creating valuable information on the concentrations of the cell's constituents. However, little is known about in vivo reaction rates, which can be determined by Metabolic Flux Analysis (MFA), a combination of isotope labeling experiments and computer modeling of the metabolic network. Large-scale applications of this method so far have been hampered by tedious procedures of tissue culture, analytics, modeling and simulation. By streamlining the workflow of MFA, the throughput of the method could be significantly increased. We propose strategies for these improvements on various sub-steps which will move flux analysis to the medium-throughput range and closer to established methods such as metabolite profiling. Furthermore, this may enable novel applications of MFA, for example screening plant populations for traits related to the flux phenotype.     

Parallel labeling experiments with [U-(13)C]glucose validate E. coli metabolic network model for (13)C metabolic flux analysis

(13)C-metabolic flux analysis (MFA) is a widely used method for measuring intracellular metabolic fluxes in living cells. (13)C MFA relies on several key assumptions: (1) the assumed metabolic network model is complete, in that it accounts for all significant enzymatic and transport reactions; (2) (13)C-labeling measurements are accurate and precise; and (3) enzymes and transporters do not discriminate between (12)C- and (13)C-labeled metabolites. In this study, we tested these inherent assumptions of (13)C MFA for wild-type E. coli by parallel labeling experiments with [U-(13)C]glucose as tracer. Cells were grown in six parallel cultures in custom-constructed mini-bioreactors, starting from the same inoculum, on medium containing different mixtures of natural glucose and fully labeled [U-(13)C]glucose, ranging from 0% to 100% [U-(13)C]glucose. Macroscopic growth characteristics of E. coli showed no observable kinetic isotope effect. The cells grew equally well on natural glucose, 100% [U-(13)C]glucose, and mixtures thereof. (13)C MFA was then used to determine intracellular metabolic fluxes for several metabolic network models: an initial network model from literature; and extended network models that accounted for potential dilution effects of isotopic labeling. The initial network model did not give statistically acceptable fits and produced inconsistent flux results for the parallel labeling experiments. In contrast, an extended network model that accounted for dilution of intracellular CO(2) by exchange with extracellular CO(2) produced statistically acceptable fits, and the estimated metabolic fluxes were consistent for the parallel cultures. This study illustrates the importance of model validation for (13)C MFA. We show that an incomplete network model can produce statistically unacceptable fits, as determined by a chi-square test for goodness-of-fit, and return biased metabolic fluxes. The validated metabolic network model for E. coli from this study can be used in future investigations for unbiased metabolic flux measurements.     

The topology of metabolic isotope labeling networks

Background  

Metabolic Flux Analysis (MFA) based on isotope labeling experiments (ILEs) is a widely established tool for determining fluxes in metabolic pathways. Isotope labeling networks (ILNs) contain all essential information required to describe the flow of labeled material in an ILE. Whereas recent experimental progress paves the way for high-throughput MFA, large network investigations and exact statistical methods, these developments are still limited by the poor performance of computational routines used for the evaluation and design of ILEs. In this context, the global analysis of ILN topology turns out to be a clue for realizing large speedup factors in all required computational procedures.     

Thermodynamics-based metabolic flux analysis

A new form of metabolic flux analysis (MFA) called thermodynamics-based metabolic flux analysis (TMFA) is introduced with the capability of generating thermodynamically feasible flux and metabolite activity profiles on a genome scale. TMFA involves the use of a set of linear thermodynamic constraints in addition to the mass balance constraints typically used in MFA. TMFA produces flux distributions that do not contain any thermodynamically infeasible reactions or pathways, and it provides information about the free energy change of reactions and the range of metabolite activities in addition to reaction fluxes. TMFA is applied to study the thermodynamically feasible ranges for the fluxes and the Gibbs free energy change, Delta(r)G', of the reactions and the activities of the metabolites in the genome-scale metabolic model of Escherichia coli developed by Palsson and co-workers. In the TMFA of the genome scale model, the metabolite activities and reaction Delta(r)G' are able to achieve a wide range of values at optimal growth. The reaction dihydroorotase is identified as a possible thermodynamic bottleneck in E. coli metabolism with a Delta(r)G' constrained close to zero while numerous reactions are identified throughout metabolism for which Delta(r)G' is always highly negative regardless of metabolite concentrations. As it has been proposed previously, these reactions with exclusively negative Delta(r)G' might be candidates for cell regulation, and we find that a significant number of these reactions appear to be the first steps in the linear portion of numerous biosynthesis pathways. The thermodynamically feasible ranges for the concentration ratios ATP/ADP, NAD(P)/NAD(P)H, and H(extracellular)(+)/H(intracellular)(+) are also determined and found to encompass the values observed experimentally in every case. Further, we find that the NAD/NADH and NADP/NADPH ratios maintained in the cell are close to the minimum feasible ratio and maximum feasible ratio, respectively.     

A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.     

Combining metabolic flux analysis tools and 13C NMR to estimate intracellular fluxes of cultured astrocytes

In this work, brain cell metabolism was investigated by (13)C NMR spectroscopy and metabolic flux analysis (MFA). Monotypic cultures of astrocytes were incubated with labeled glucose for 38 h, and the distribution of the label was analyzed by (13)C NMR spectroscopy. The analysis of the spectra reveals two distinct physiological states characterized by different ratios of pyruvate carboxylase to pyruvate dehydrogenase activities (PC/PDH). Intracellular flux distributions for both metabolic states were estimated by MFA using the isotopic information and extracellular rate measurements as constraints. The model was subsequently checked with the consistency index method. From a biological point of view, the occurrence of the two physiological states appears to be correlated with the presence or absence of extracellular glutamate. Concerning the model, it can be stated that the metabolic network and the set of constraints adopted provide a consistent and robust characterization of the astrocytic metabolism, allowing for the calculation of central intracellular fluxes such as pyruvate recycling, the anaplerotic flux mediated by pyruvate carboxylase, and the glutamine formation through glutamine synthetase.     

Insights into metabolic efficiency from flux analysis

The efficiency of carbon and energy flows throughout metabolism defines the potential for growth and reproductive success of plants. Understanding the basis for metabolic efficiency requires relevant definitions of efficiency as well as measurements of biochemical functions through metabolism. Here insights into the basis of efficiency provided by (13)C-based metabolic flux analysis (MFA) as well as the uses and limitations of efficiency in predictive flux balance analysis (FBA) are highlighted. (13)C-MFA studies have revealed unusual features of central metabolism in developing green seeds for the efficient use of light to conserve carbon and identified metabolic inefficiencies in plant metabolism due to dissipation of ATP by substrate cycling. Constraints-based FBA has used efficiency to guide the prediction of the growth and actual internal flux distribution of plant systems. Comparisons in a few cases have been made between flux maps measured by (13)C-based MFA and those predicted by FBA assuming one or more maximal efficiency parameters. These studies suggest that developing plant seeds and photoautotrophic microorganisms may indeed have patterns of metabolic flux that maximize efficiency. MFA and FBA are synergistic toolsets for uncovering and explaining the metabolic basis of efficiencies and inefficiencies in plant systems.     

Computational tools for isotopically instationary 13C labeling experiments under metabolic steady state conditions

(13)C metabolic flux analysis (MFA) has become an important and powerful tool for the quantitative analysis of metabolic networks in the framework of metabolic engineering. Isotopically instationary (13)C MFA under metabolic stationary conditions is a promising refinement of classical stationary MFA. It accounts for the experimental requirements of non-steady-state cultures as well as for the shortening of the experimental duration. This contribution extends all computational methods developed for classical stationary (13)C MFA to the instationary situation by using high-performance computing methods. The developed tools allow for the simulation of instationary carbon labeling experiments (CLEs), sensitivity calculation with respect to unknown parameters, fitting of the model to the measured data, statistical identifiability analysis and an optimal experimental design facility. To explore the potential of the new approach all these tools are applied to the central metabolism of Escherichia coli. The achieved results are compared to the outcome of the stationary counterpart, especially focusing on statistical properties. This demonstrates the specific strengths of the instationary method. A new ranking method is proposed making both an a priori and an a posteriori design of the sampling times available. It will be shown that although still not all fluxes are identifiable, the quality of flux estimates can be strongly improved in the instationary case. Moreover, statements about the size of some immeasurable pool sizes can be made.     

Validation-based model selection for 13C metabolic flux analysis with uncertain measurement errors

Accurate measurements of metabolic fluxes in living cells are central to metabolism research and metabolic engineering. The gold standard method is model-based metabolic flux analysis (MFA), where fluxes are estimated indirectly from mass isotopomer data with the use of a mathematical model of the metabolic network. A critical step in MFA is model selection: choosing what compartments, metabolites, and reactions to include in the metabolic network model. Model selection is often done informally during the modelling process, based on the same data that is used for model fitting (estimation data). This can lead to either overly complex models (overfitting) or too simple ones (underfitting), in both cases resulting in poor flux estimates. Here, we propose a method for model selection based on independent validation data. We demonstrate in simulation studies that this method consistently chooses the correct model in a way that is independent on errors in measurement uncertainty. This independence is beneficial, since estimating the true magnitude of these errors can be difficult. In contrast, commonly used model selection methods based on the χ²-test choose different model structures depending on the believed measurement uncertainty; this can lead to errors in flux estimates, especially when the magnitude of the error is substantially off. We present a new approach for quantification of prediction uncertainty of mass isotopomer distributions in other labelling experiments, to check for problems with too much or too little novelty in the validation data. Finally, in an isotope tracing study on human mammary epithelial cells, the validation-based model selection method identified pyruvate carboxylase as a key model component. Our results argue that validation-based model selection should be an integral part of MFA model development.     

Investigation of metabolic objectives in cultured hepatocytes

Using optimization based methods to predict fluxes in metabolic flux balance models has been a successful approach for some microorganisms, enabling construction of in silico models and even inference of some regulatory motifs. However, this success has not been translated to mammalian cells. The lack of knowledge about metabolic objectives in mammalian cells is a major obstacle that prevents utilization of various metabolic engineering tools and methods for tissue engineering and biomedical purposes. In this work, we investigate and identify possible metabolic objectives for hepatocytes cultured in vitro. To achieve this goal, we present a special data-mining procedure for identifying metabolic objective functions in mammalian cells. This multi-level optimization based algorithm enables identifying the major fluxes in the metabolic objective from MFA data in the absence of information about critical active constraints of the system. Further, once the objective is determined, active flux constraints can also be identified and analyzed. This information can be potentially used in a predictive manner to improve cell culture results or clinical metabolic outcomes. As a result of the application of this method, it was found that in vitro cultured hepatocytes maximize oxygen uptake, coupling of urea and TCA cycles, and synthesis of serine and urea. Selection of these fluxes as the metabolic objective enables accurate prediction of the flux distribution in the system given a limited amount of flux data; thus presenting a workable in silico model for cultured hepatocytes. It is observed that an overall homeostasis picture is also emergent in the findings.     

Interpretation of metabolic flux maps by limitation potentials and constrained limitation sensitivities

Two new concepts, "Limitation Potential" and "Constraint Limitation Sensitivity" are introduced that use definitions derived from metabolic flux analysis (MFA) and metabolic network analysis (MNA). They are applied to interpret a measured flux distribution in the context of all possible flux distributions and thus combine MFA with MNA. The proposed measures are used to quantify and compare the influence of intracellular fluxes on the production yield. The methods are purely based on the stoichiometry of the network and constraints that are given from irreversible fluxes. In contrast to metabolic control analysis (MCA), within this approach no information about the kinetic mechanisms are needed. A limitation potential (LP) is defined as the reduction of the reachable (theoretical) maximum by a measured flux. Measured fluxes that strongly narrow the reachable maximum are assumed to be limiting as the network has no ability to counterbalance the restriction due to the observed flux. In a second step, the sensitivity of the reduced maximum is regarded. This measure provides information about the necessitated changes to reach higher yields. The methods are applied to interpret the capabilities of a network based on measured fluxes for a L-phenylalanine producer. The strain was examined by a series of experiments and three flux maps of the production phase are analyzed. It can be shown that the reachable yield is drastically reduced by the measured efflux into the TCA cycle, while the oxidative pentose-phosphate pathway only plays a secondary role on the reachable maximum.     

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.     

Metabolic flux analysis as a tool in metabolic engineering of plants

Methods of metabolic flux analysis (MFA) provide insights into the theoretical capabilities of metabolic networks and allow probing the in vivo performance of cellular metabolism. In recent years, an increasing awareness has developed that network analysis methods within the systems biology toolbox are serving to improve our understanding and ability to manipulate metabolism. In this minireview the potential of MFA to increase the chances of success in metabolic engineering of plants is presented, recent progress related to engineering and flux analysis in central metabolism of plants is discussed, and some recent advances in flux analysis methodology are highlighted.