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

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

Systematic evaluation of parameters for genome-scale metabolic models of cultured mammalian cells

Genome-scale metabolic models describe cellular metabolism with mechanistic detail. Given their high complexity, such models need to be parameterized correctly to yield accurate predictions and avoid overfitting. Effective parameterization has been well-studied for microbial models, but it remains unclear for higher eukaryotes, including mammalian cells. To address this, we enumerated model parameters that describe key features of cultured mammalian cells – including cellular composition, bioprocess performance metrics, mammalian-specific pathways, and biological assumptions behind model formulation approaches. We tested these parameters by building thousands of metabolic models and evaluating their ability to predict the growth rates of a panel of phenotypically diverse Chinese Hamster Ovary cell clones. We found the following considerations to be most critical for accurate parameterization: (1) cells limit metabolic activity to maintain homeostasis, (2) cell morphology and viability change dynamically during a growth curve, and (3) cellular biomass has a particular macromolecular composition. Depending on parameterization, models predicted different metabolic phenotypes, including contrasting mechanisms of nutrient utilization and energy generation, leading to varying accuracies of growth rate predictions. Notably, accurate parameter values broadly agreed with experimental measurements. These insights will guide future investigations of mammalian metabolism.     

Interrogating the effect of enzyme kinetics on metabolism using differentiable constraint-based models

Metabolic models are typically characterized by a large number of parameters. Traditionally, metabolic control analysis is applied to differential equation-based models to investigate the sensitivity of predictions to parameters. A corresponding theory for constraint-based models is lacking, due to their formulation as optimization problems. Here, we show that optimal solutions of optimization problems can be efficiently differentiated using constrained optimization duality and implicit differentiation. We use this to calculate the sensitivities of predicted reaction fluxes and enzyme concentrations to turnover numbers in an enzyme-constrained metabolic model of Escherichia coli. The sensitivities quantitatively identify rate limiting enzymes and are mathematically precise, unlike current finite difference based approaches used for sensitivity analysis. Further, efficient differentiation of constraint-based models unlocks the ability to use gradient information for parameter estimation. We demonstrate this by improving, genome-wide, the state-of-the-art turnover number estimates for E. coli. Finally, we show that this technique can be generalized to arbitrarily complex models. By differentiating the optimal solution of a model incorporating both thermodynamic and kinetic rate equations, the effect of metabolite concentrations on biomass growth can be elucidated. We benchmark these metabolite sensitivities against a large experimental gene knockdown study, and find good alignment between the predicted sensitivities and in vivo metabolome changes. In sum, we demonstrate several applications of differentiating optimal solutions of constraint-based metabolic models, and show how it connects to classic metabolic control analysis.     

Metabolomics assisted metabolic network modeling and network wide analysis of metabolites in microbiology

Metabolomics is the science of qualitatively and quantitatively analyzing low molecular weight metabolites occur in a given biological system. It provides valuable information to elucidate the functional roles and relations of different metabolites in a metabolic pathway. In recent years, a large amount of research on microbial metabolomics has been conducted. It has become a useful tool for achieving highly efficient synthesis of target metabolites. At the same time, many studies have been conducted over the years in order to integrate metabolomics data into metabolic network modeling, which has yielded many exciting results. Additionally, metabolomics also shows great advantages in analyzing the relationship of metabolites network wide. Integrating metabolomics data into metabolic network construction and applying it in network wide analysis of cell metabolism would further improve our ability to control cellular metabolism and optimize the design of cell factories for the overproduction of valuable biochemicals. This review will examine recent progress in the application of metabolomics approaches in metabolic network modeling and network wide analysis of microbial cell metabolism.     

Elevated pCO2 affects the lactate metabolic shift in CHO cell culture processes

The shift from lactate production to consumption in CHO cell metabolism is a key event during cell culture cultivations and is connected to increased culture longevity and final product titers. However, the mechanisms controlling this metabolic shift are not yet fully understood. Variations in lactate metabolism have been mainly reported to be induced by process pH and availability of substrates like glucose and glutamine. The aim of this study was to investigate the effects of elevated pCO₂ concentrations on the lactate metabolic shift phenomena in CHO cell culture processes. In this publication, we show that at elevated pCO₂ in batch and fed‐batch cultures, the lactate metabolic shift was absent in comparison to control cultures at lower pCO₂ values. Furthermore, through metabolic flux analysis we found a link between the lactate metabolic shift and the ratio of NADH producing and regenerating intracellular pathways. This ratio was mainly affected by a reduced oxidative capacity of cultures at elevated pCO₂. The presented results are especially interesting for large‐scale and perfusion processes where increased pCO₂ concentrations are likely to occur. Our results suggest, that so far unexplained metabolic changes may be connected to increased pCO₂ accumulation in larger scale fermentations. Finally, we propose several mechanisms through which increased pCO₂ might affect the cell metabolism and briefly discuss methods to enable the lactate metabolic shift during cell cultivations.     

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 ¹³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.     

Physiology of microbial cells and metabolic engineering

This review is devoted to the problems of the physiology and cell biology of microorganisms in relation to metabolic engineering. The latter is considered as a branch of fundamental and applied biotechnology aimed at controlling microbial metabolism by methods of genetic engineering and classical genetics and based on intimate knowledge of cell metabolism. Attention is also given to the problems associated with the metabolic limitation of microbial biosyntheses, analysis and control of metabolic fluxes, rigidity of metabolic pathways, the role of pleiotropic (global) regulatory systems in the control of metabolic fluxes, and prospects of physiological and evolutionary approaches in metabolic engineering.     

Lethality and synthetic lethality in the genome-wide metabolic network of Escherichia coli

Recent genomic analyses on the cellular metabolic network show that reaction flux across enzymes are diverse and exhibit power-law behavior in its distribution. While intuition might suggest that the reactions with larger fluxes are more likely to be lethal under the blockade of its catalysing gene products or gene knockouts, we find, by in silico flux analysis, that the lethality rarely has correlations with the flux level owing to the widespread backup pathways innate in the genome-wide metabolism of Escherichia coli. Lethal reactions, of which the deletion generates cascading failure of following reactions up to the biomass reaction, are identified in terms of the Boolean network scheme as well as the flux balance analysis. The avalanche size of a reaction, defined as the number of subsequently blocked reactions after its removal, turns out to be a useful measure of lethality. As a means to elucidate phenotypic robustness to a single deletion, we investigate synthetic lethality in reaction level, where simultaneous deletion of a pair of nonlethal reactions leads to the failure of the biomass reaction. Synthetic lethals identified via flux balance and Boolean scheme are consistently shown to act in parallel pathways, working in such a way that the backup machinery is compromised.     

Exploring local structural organization of metabolic networks using subgraph patterns

Metabolic networks of many cellular organisms share global statistical features. Their connectivity distributions follow the long-tailed power law and show the small-world property. In addition, their modular structures are organized in a hierarchical manner. Although the global topological organization of metabolic networks is well understood, their local structural organization is still not clear. Investigating local properties of metabolic networks is necessary to understand the nature of metabolism in living organisms. To identify the local structural organization of metabolic networks, we analysed the subgraphs of metabolic networks of 43 organisms from three domains of life. We first identified the network motifs of metabolic networks and identified the statistically significant subgraph patterns. We then compared metabolic networks from different domains and found that they have similar local structures and that the local structure of each metabolic network has its own taxonomical meaning. Organisms closer in taxonomy showed similar local structures. In addition, the common substrates of 43 metabolic networks were not randomly distributed, but were more likely to be constituents of cohesive subgraph patterns.     

Genetically constrained metabolic flux analysis

Significant progress has been made in using existing metabolic databases to estimate metabolic fluxes. Traditional metabolic flux analysis generally starts with a predetermined metabolic network. This approach has been employed successfully to analyze the behaviors of recombinant strains by manually adding or removing the corresponding pathway(s) in the metabolic map. The current work focuses on the development of a new framework that utilizes genomic and metabolic databases, including available genetic/regulatory network structures and gene chip expression data, to constrain metabolic flux analysis. The genetic network consisting of the sensing/regulatory circuits will activate or deactivate a specific set of genes in response to external stimulus. The activation and/or repression of this set of genes will result in different gene expression levels that will in turn change the structure of the metabolic map. Hence, the metabolic map will automatically "adapt" to the external stimulus as captured by the genetic network. This adaptation selects a subnetwork from the pool of feasible reactions and so performs what we term "environmentally driven dimensional reduction." The Escherichia coli oxygen and redox sensing/regulatory system, which controls the metabolic patterns connected to glycolysis and the TCA cycle, was used as a model system to illustrate the proposed approach.     

Structural analysis of metabolic networks based on flux centrality

Metabolic reactions are fundamental to living organisms, and a large number of reactions simultaneously occur at a given time in living cells transforming diverse metabolites into each other. There has been an ongoing debate on how to classify metabolites with respect to their importance for metabolic performance, usually based on the analysis of topological properties of genome scale metabolic networks. However, none of these studies have accounted quantitatively for flux in metabolic networks, thus lacking an important component of a cell’s biochemistry.We therefore analyzed a genome scale metabolic network of Escherichia coli by comparing growth under 19 different growth conditions, using flux balance analysis and weighted network centrality investigation. With this novel concept of flux centrality we generated metabolite rankings for each particular growth condition. In contrast to the results of conventional analysis of genome scale metabolic networks, different metabolites were top-ranking dependent on the growth condition. At the same time, several metabolites were consistently among the high ranking ones. Those are associated with pathways that have been described by biochemists as the most central part of metabolism, such as glycolysis, tricarboxylic acid cycle and pentose phosphate pathway. The values for the average path length of the analyzed metabolite networks were between 10.5 and 12.6, supporting recent findings that the metabolic network of E. coli is not a small-world network.     

Genomically and biochemically accurate metabolic reconstruction of Methanosarcina barkeri Fusaro,iMG746

Methanosarcina barkeri is an Archaeon that produces methane anaerobically as the primary byproduct of its metabolism. M. barkeri can utilize several substrates for ATP and biomass production including methanol, acetate, methyl amines, and a combination of hydrogen and carbon dioxide. In 2006, a metabolic reconstruction of M. barkeri, iAF692, was generated based on a draft genome annotation. The iAF692 reconstruction enabled the first genome-Scale simulations for Archaea. Since the publication of the first metabolic reconstruction of M. barkeri, additional genomic, biochemical, and phenotypic data have clarified several metabolic pathways. We have used this newly available data to improve the M. barkeri metabolic reconstruction. Modeling simulations using the updated model, iMG746, have led to increased accuracy in predicting gene knockout phenotypes and simulations of batch growth behavior. We used the model to examine knockout lethality data and make predictions about metabolic regulation under different growth conditions. Thus, the updated metabolic reconstruction of M. barkeri metabolism is a useful tool for predicting cellular behavior, studying the methanogenic lifestyle, guiding experimental studies, and making predictions relevant to metabolic engineering applications.     

Genome-scale metabolic modeling of P. thermoglucosidasius NCIMB 11955 reveals metabolic bottlenecks in anaerobic metabolism

Parageobacillus thermoglucosidasius represents a thermophilic, facultative anaerobic bacterial chassis, with several desirable traits for metabolic engineering and industrial production. To further optimize strain productivity, a systems level understanding of its metabolism is needed, which can be facilitated by a genome-scale metabolic model. Here, we present p-thermo, the most complete, curated and validated genome-scale model (to date) of Parageobacillus thermoglucosidasius NCIMB 11955. It spans a total of 890 metabolites, 1175 reactions and 917 metabolic genes, forming an extensive knowledge base for P. thermoglucosidasius NCIMB 11955 metabolism. The model accurately predicts aerobic utilization of 22 carbon sources, and the predictive quality of internal fluxes was validated with previously published ¹³C-fluxomics data. In an application case, p-thermo was used to facilitate more in-depth analysis of reported metabolic engineering efforts, giving additional insight into fermentative metabolism. Finally, p-thermo was used to resolve a previously uncharacterised bottleneck in anaerobic metabolism, by identifying the minimal required supplemented nutrients (thiamin, biotin and iron(III)) needed to sustain anaerobic growth. This highlights the usefulness of p-thermo for guiding the generation of experimental hypotheses and for facilitating data-driven metabolic engineering, expanding the use of P. thermoglucosidasius as a high yield production platform.     

Drug-efficacy depends on the inhibitor type and the target position in a metabolic network--a systematic study

Drug discovery usually focuses on candidate molecules that affect individual reactions with presumed essential functions in the cellular reaction network, especially in the development of diseases. Unfortunately, appropriately designed drugs often fail to show the expected biological effect, since the multitude of interactions in the biochemical reaction network buffers the individual changes or causes significant side effects. We address this problem through a computational approach, which considers the effect of drug application within a generalized biochemical pathway and by studying the effect of changes regarding the type and strength of inhibitors on the reduction of flux. This allows us to systematically search for the appropriate target and for type and concentration of the optimal inhibitor. We propose the flux selectivity as a measure for the discrimination of the effect on different pathways. Since the calculation of the flux selectivity is based on flux control coefficients that are calculated in the non-affected state, it is also a means for predicting the inhibitor efficacy. Furthermore, we will propose how to increase discriminative inhibition in the case of a parasitic disease by using multi-target drugs.This work is devoted to the memorial of our teacher Reinhart Heinrich, who made important contributions to the investigation of the regulation of metabolic networks, namely by introducing and applying the concept of metabolic control.     

Optimizing metabolic pathways by screening for feasible synthetic reactions

Background

Reconstruction of genome-scale metabolic networks has resulted in models capable of reproducing experimentally observed biomass yield/growth rates and predicting the effect of alterations in metabolism for biotechnological applications. The existing studies rely on modifying the metabolic network of an investigated organism by removing or inserting reactions taken either from evolutionary similar organisms or from databases of biochemical reactions (e.g., KEGG). A potential disadvantage of these knowledge-driven approaches is that the result is biased towards known reactions, as such approaches do not account for the possibility of including novel enzymes, together with the reactions they catalyze.

Results

Here, we explore the alternative of increasing biomass yield in three model organisms, namely Bacillus subtilis, Escherichia coli, and Hordeum vulgare, by applying small, chemically feasible network modifications. We use the predicted and experimentally confirmed growth rates of the wild-type networks as reference values and determine the effect of inserting mass-balanced, thermodynamically feasible reactions on predictions of growth rate by using flux balance analysis.

Conclusions

While many replacements of existing reactions naturally lead to a decrease or complete loss of biomass production ability, in all three investigated organisms we find feasible modifications which facilitate a significant increase in this biological function. We focus on modifications with feasible chemical properties and a significant increase in biomass yield. The results demonstrate that small modifications are sufficient to substantially alter biomass yield in the three organisms. The method can be used to predict the effect of targeted modifications on the yield of any set of metabolites (e.g., ethanol), thus providing a computational framework for synthetic metabolic engineering.     

基于酶约束的代谢网络模型研究进展及其应用

基因组尺度代谢网络模型已经成功地应用于指导代谢工程改造,但由于传统通量平衡分析法仅考虑化学计量学和反应方向约束,模拟得到的是理论最优结果,对一些现象如代谢溢流、底物层级利用等无法准确描述。近年来人们通过在代谢网络模型中引入新的蛋白量、热力学等约束发展了新的约束优化计算方法,可以更准确真实地模拟细胞在不同条件下的代谢行为。文中主要对近年来提出的多种酶约束模型进行评述,对酶约束引入的基本思路、酶约束的数学方程表示及优化目标设定、引入酶约束后对代谢通量计算结果的影响及酶约束模型在代谢工程菌种改造中的应用等进行了全面深入的介绍,并提出了已有各种方法存在的主要问题,展望了相关方法的未来发展方向。通过引入新的约束,代谢网络模型能够更精确模拟和预测细胞在环境和基因扰动下的代谢行为,为代谢工程菌种改造提供更准确可靠的指导。     

Modeling metabolic adjustment in <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> upon treatment with isoniazid

Complex biological systems exhibit a property of robustness at all levels of organization. Through different mechanisms, the system tries to sustain stress such as due to starvation or drug exposure. To explore whether reconfiguration of the metabolic networks is used as a means to achieve robustness, we have studied possible metabolic adjustments in Mtb upon exposure to isoniazid (INH), a front-line clinical drug. The redundancy in the genome of M. tuberculosis (Mtb) makes it an attractive system to explore if alternate routes of metabolism exist in the bacterium. While the mechanism of action of INH is well studied, its effect on the overall metabolism is not well characterized. Using flux balance analysis, inhibiting the fluxes flowing through the reactions catalyzed by Rv1484, the target of INH, significantly changes the overall flux profiles. At the pathway level, activation or inactivation of certain pathways distant from the target pathway, are seen. Metabolites such as NADPH are shown to reduce drastically, while fatty acids tend to accumulate. The overall biomass also decreases with increasing inhibition levels. Inhibition studies, pathway level clustering and comparison of the flux profiles with the gene expression data indicate the activation of folate metabolism, ubiquinone metabolism, and metabolism of certain amino acids. This analysis provides insights useful for target identification and designing strategies for combination therapy. Insights gained about the role of individual components of a system and their interactions will also provide a basis for reconstruction of whole systems through synthetic biology approaches.     

Analysis of phenetic trees based on metabolic capabilites across the three domains of life

Here, we used data of complete genomes to study comparatively the metabolism of different species. We built phenetic trees based on the enzymatic functions present in different parts of metabolism. Seven broad metabolic classes, comprising a total of 69 metabolic pathways, were comparatively analyzed for 27 fully sequenced organisms of the domains Eukarya, Bacteria and Archaea. Phylogenetic profiles based on the presence/absence of enzymatic functions for each metabolic class were determined and distance matrices for all the organisms were then derived from the profiles. Unrooted phenetic trees based upon the matrices revealed the distribution of the organisms according to their metabolic capabilities, reflecting the ecological pressures and adaptations that those species underwent during their evolution. We found that organisms that are closely related in phylogenetic terms could be distantly related metabolically and that the opposite is also true. For example, obligate bacterial pathogens were usually grouped together in our metabolic trees, demonstrating that obligate pathogens share common metabolic features regardless of their diverse phylogenetic origins. The branching order of proteobacteria often did not match their classical phylogenetic classification and Gram-positive bacteria showed diverse metabolic affinities. Archaea were found to be metabolically as distant from free-living bacteria as from eukaryotes, and sometimes were placed close to the metabolically highly specialized group of obligate bacterial pathogens. Metabolic trees represent an integrative approach for the comparison of the evolution of the metabolism and its correlation with the evolution of the genome, helping to find new relationships in the tree of life.     

13C metabolic flux analysis of three divergent extremely thermophilic bacteria: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252

Thermophilic organisms are being increasingly investigated and applied in metabolic engineering and biotechnology. The distinct metabolic and physiological characteristics of thermophiles, including broad substrate range and high uptake rates, coupled with recent advances in genetic tool development, present unique opportunities for strain engineering. However, poor understanding of the cellular physiology and metabolism of thermophiles has limited the application of systems biology and metabolic engineering tools to these organisms. To address this concern, we applied high resolution ¹³C metabolic flux analysis to quantify fluxes for three divergent extremely thermophilic bacteria from separate phyla: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252. We performed 18 parallel labeling experiments, using all singly labeled glucose tracers for each strain, reconstructed and validated metabolic network models, measured biomass composition, and quantified precise metabolic fluxes for each organism. In the process, we resolved many uncertainties regarding gaps in pathway reconstructions and elucidated how these organisms maintain redox balance and generate energy. Overall, we found that the metabolisms of the three thermophiles were highly distinct, suggesting that adaptation to growth at high temperatures did not favor any particular set of metabolic pathways. All three strains relied heavily on glycolysis and TCA cycle to generate key cellular precursors and cofactors. None of the investigated organisms utilized the Entner-Doudoroff pathway and only one strain had an active oxidative pentose phosphate pathway. Taken together, the results from this study provide a solid foundation for future model building and engineering efforts with these and related thermophiles.