BioMed Research International

BioMed Research International / 2010 / Article
Special Issue

Computational Approaches in Metabolic Engineering

View this Special Issue

Editorial | Open Access

Volume 2010 |Article ID 207414 | 7 pages |

Computational Approaches in Metabolic Engineering

Received31 Dec 2010
Accepted31 Dec 2010
Published07 Apr 2011

1. Introduction

Metabolic engineering involves the adjustment of metabolic and regulatory processes to improve desired cellular behaviors, such as the production of proteins and chemicals. Since cellular metabolic and regulatory networks are often large and complex, the construction and analysis of computational models of these networks can be useful for identifying current network states and evaluating the effects of network perturbations on desired phenotypes. This special issue includes papers that illustrate how computational approaches can be used in metabolic engineering. Here, we provide a brief overview of several established computational approaches that can be used to aid in the engineering of metabolic networks, while describing some of the exciting recent advances in these fields.

2. Descriptive Approaches for Identifying Intracellular Metabolic States

A variety of experimental measurements can be used to quantify the state of metabolic and regulatory networks, including flux analysis, gene expression, protein expression, metabolite concentration, enzyme activity, uptake and secretion rates, and transcription factor-DNA binding assays. Computational models and approaches can be useful for integrating and analyzing such datasets to quantify metabolic fluxes and uncover their regulatory properties.

2.1. 13C Labeling Experiments and Metabolic Flux Analysis

The ability to quantitatively map intracellular fluxes using metabolic flux analysis (MFA) is critical for identifying pathway bottlenecks and elucidating network regulation in biological systems, especially in engineered cells with nonnative metabolic capacities [1, 2]. 13C-MFA experiments involve feeding isotopically labeled substrates to cells, tissues, or whole organisms and subsequently measuring the patterns of isotope incorporation that occur in intracellular metabolites or secreted products [3]. Both mass spectrometry (MS) and nuclear magnetic resonance (NMR) can be used to quantify the relative abundance of different “isotopomers” (i.e., isotope isomers) associated with each measured biomolecule. Because different metabolic pathways often give rise to distinct isotope labeling patterns, isotopomer abundances can be used to infer the relative fluxes through these pathways [4]. As a result, isotopomer measurements obtained from MS or NMR, in combination with direct measurements of extracellular uptake or secretion rates, can be computationally analyzed to reconstruct comprehensive flux maps describing intracellular metabolism, which is the essence of MFA. A complete flux map is thus the phenotypic equivalent of the transcriptional map obtained from DNA microarrays, and tracking flux changes in response to targeted perturbations can provide important information about the distribution of kinetic and regulatory controls in metabolism [5].

The standard approach for computing a flux map involves a nonlinear least-squares regression to minimize the lack-of-fit between (i) experimentally measured and (ii) computationally simulated data. The latter are derived by solving a “forward problem”, which involves solving the isotopomer and metabolite balance equations for a particular set of flux parameters in order to calculate the relative abundance of all observable metabolite labeling patterns. Because the forward problem must be solved up to hundreds or even thousands of times to achieve an optimal fit, a great deal of effort has been placed on developing improved strategies to simulate the isotopic labeling produced by a particular network state. The Elementary Metabolite Unit (EMU) approach was developed by Antoniewicz et al. [6] to address precisely this problem. Through a novel decomposition of the isotopomer network, the algorithm systematically identifies the minimal set of variables required to simulate the available labeling measurements. The EMU approach has achieved 10-fold reductions in the system sizes required to simulate 13C labeling in both medium- and large-scale networks [6, 7]. While these gains are impressive, the true power of the EMU approach is in its ability to unlock entirely new isotope labeling strategies that were previously outside the reach of computational tractability. For instance, the number of isotopomer variables required to simulate mixed 13C, 18O, and 2H labeling in even a small-scale network of gluconeogenesis (24 reactions, 14 intracellular metabolites) climbs into the millions, while the number of EMU balances is in the mere hundreds [6].

Typically, MFA relies on the assumption of both metabolic and isotopic steady state. Achieving this situation experimentally involves (i) equilibrating the system in a stable metabolic state, (ii) introducing an isotopically labeled substrate without perturbing the metabolic steady state, (iii) allowing the system to establish a new isotopic steady state that reflects the underlying metabolic fluxes, and (iv) measuring isotopic labeling in the fully equilibrated system. Depending on the relative speed of metabolic and isotopic dynamics, however, other experimental scenarios can be envisioned. If the isotopic labeling responds quickly to any metabolic changes in the system, quasi-stationary MFA can be applied to obtain a series of snapshots that describe the variation in network fluxes over time [8, 9]. Conversely, if labeling occurs slowly but metabolism is maintained in a fixed state, isotopically nonstationary MFA (INST-MFA) can be used to estimate fluxes [10]. Finally, when measurements are obtained under both metabolically and isotopically nonstationary conditions, a fully dynamic modeling approach is required to estimate fluxes [11]. In the current issue, Lequeux et al. (G. Lequeux et al. “Dynamic metabolic flux analysis demonstrated on cultures where the limiting substrate is changed from carbon to nitrogen and vice versa”) have taken a different dynamic MFA approach by obtaining transient measurements of 10 different extracellular metabolites during the shift of E. coli cells from nitrogen- to carbon-limitation, or vice versa. Numerically differentiating these extracellular measurements allowed the authors to estimate dynamically changing uptake and secretion fluxes, which provided the measurements necessary to estimate all intracellular fluxes in their model as a function of time. This has the advantage of avoiding the complications imposed by measuring and fitting transient isotope labeling data, but at the same time provides limited redundancy to validate assumptions on cofactor balancing and respiration efficiency that must be invoked in order to close the system of balance equations [12].

By applying the EMU approach, Young et al. [13] have recently developed computational routines that achieve more than 5,000-fold speedup relative to prior INST-MFA algorithms. This opens the door to several practical applications that were previously intractable due to the computational complexity of INST-MFA, where the isotopomer balances are described by differential rather than algebraic equations. INST-MFA is ideally suited to systems that label slowly due to the presence of large intermediate pools or pathway bottlenecks. This approach not only avoids the additional time and cost of feeding isotope tracers over extended periods [14], but may become absolutely necessary in cases where the system cannot be held in a fixed metabolic state long enough to allow isotopic labeling to fully equilibrate. As a result, INST-MFA is expected to become an indispensible tool for extending MFA approaches to studies of mammalian systems [1517], industrial bioprocesses [10, 18], and other scenarios where obtaining a strict isotopic steady state may be impractical. Another emerging application of INST-MFA is its application to photoautotrophic metabolism, which is the process by which plants, algae, and cyanobacteria use light energy to fix carbon dioxide into complex organic molecules. Because photoautotrophs assimilate carbon solely from CO2, feeding 13CO2 will produce a uniform steady-state labeling pattern that is insensitive to fluxes. Thus, conventional steady-state 13C-MFA is incapable of quantifying autotrophic metabolic fluxes [19]. However, transient measurements of isotope incorporation following a step change from unlabeled to labeled CO2 can be used to estimate fluxes by applying INST-MFA [20]. This approach enables comprehensive flux analysis of photoautotrophic metabolism, complementing previous 13C-MFA studies of plants [21] and cyanobacteria [22] that were limited to heterotrophic or mixotrophic culture conditions, with sugar as the major carbon source. Taken together, these advances illustrate how combined progress in both analytical capabilities and computational techniques are driving MFA applications toward larger, more dynamic, and more complex biochemical networks, which encompass a growing variety of plant, animal, and microbial systems.

2.2. Using Metabolite and Gene/Protein Expression Measurements

In addition to 13C-MFA experiments, other large-scale measurements can be made which capture cellular metabolic states. These include measurements of metabolite concentrations and gene or protein expression. Metabolite concentrations can now be quantified using mass spectroscopy for hundreds of metabolites in a single condition. For example, Bennett and colleagues used LC-MS/MS to quantify over 100 metabolite concentrations in E. coli cells grown under three different conditions [23]. These metabolite concentrations can be analyzed to identify potential metabolic bottlenecks by evaluating enzyme saturation and estimating Gibbs-free energy changes of reactions. Bennett and colleagues compared measured metabolite concentrations to reported Michaelis-Menten kinetic parameters ( ) to determine whether individual reaction rates are substrate (where ) or enzyme limited (where ). Their analysis found that most substrate concentrations were higher than the reported values (83%) indicating that for many reactions the rates are limited by enzyme levels [23]. Metabolite concentrations can also be used to estimate thermodynamic properties of metabolic reactions since the change in Gibbs-free energy for a reaction ( ) is dependent on substrate and product concentrations [24, 25]. The values can be estimated from measured metabolite concentrations and then used to distinguish between those reactions that are operating close to equilibrium ( ), and those reactions that are far from equilibrium ( ) whose rates may be limited by regulation via enzyme kinetics [25]. Thermodynamic analysis using metabolite concentrations has been preformed for both E. coli and S. cerevisiae to identify these types of reactions [23, 25, 26]. Together the analysis of metabolite concentrations using both kinetic (if available) and thermodynamic information can be useful to identify metabolic bottlenecks or rate limiting reactions.

In addition to metabolite concentrations, gene expression and protein expression measurements can also be used to help elucidate metabolic fluxes and their regulation. Recent modeling efforts have used these types of measurements to improve predictions of fluxes through metabolic reactions [5, 2729]. With these approaches, gene expression data is used to place restrictions on flux values or flux changes. The GIMME method uses an expression threshold and prevents flux through reactions associated with genes whose expression is below the threshold. In this case, an inconsistency score is minimized, where penalties depend on the magnitude of the flux and how far the expression is below the threshold [27]. Another method, proposed by Shlomi and colleagues, instead groups reactions into high, medium, and low sets based on expression levels of associated genes [29]. A flux distribution is then identified that has flux through as many reactions in the high set as possible, and no flux through as many reactions as possible in the low set. In the third method (E-flux), the relative expression levels of genes from a given condition are used to place constraints on the upper limits flux values can take [28]. Another approach, developed by Moxley and colleagues, instead uses the changes in expression level between two conditions to estimate the changes in fluxes between two conditions [5]. In this case, the predicted flux change depends on the expression change (ΔmRNA) and some additional model parameters. All of these methods have been applied to gene expression measurements, but they likely can also be used with quantitative proteomics measurements.

3. Predictive Approaches for Improving Cellular Phenotypes

The modeling approaches described in the previous sections can provide descriptions of the metabolic fluxes by analyzing different types of experimental measurements. Once these metabolic states are known, other modeling approaches can be used to identify which environmental and/or genetic perturbations would improve cellular phenotypes, such as the production of desired chemicals. These predictive approaches include pathway-based and optimization-based methods, that take into account the structure and stoichiometry of metabolic networks, as well as, kinetic modeling approaches that also account for enzyme kinetics.

3.1. Pathway-Based Approaches

Identification of the relevant pathways of a metabolic network is essential for finding effective metabolic engineering strategies. These pathways can also help derive minimal media requirements for an organism and assess the robustness and redundancy of key metabolic pathways. In this special issue, F. Llaneras and J. Picó review and compare four established methods used to identify relevant metabolic pathways: extreme currents, elementary modes, extreme pathways, and minimal generators (F. Llaneras and J. Picó, “Which metabolic pathways generate and characterize the flux space? A comparison among elementary modes, extreme pathways and minimal generators”). The authors recommend elementary modes for determining the metabolic impact of gene knockouts and the required pathways of a network for producing desired chemicals. The calculation of elementary modes requires knowledge of reaction stoichiometry and reversibility, making the recent modeling advances in biochemical reaction thermodynamics important to this analysis. Several algorithms and software packages have been developed for calculating the elementary modes of a metabolic network, including Metatool [30, 31], FluxAnalyzer [32], and functionalities built into OptFlux [33]. However, the problem of combinatorial explosion when calculating elementary modes for large complex metabolic networks has been documented [34]. This has paved the way for clustering algorithms such as the Agglomeration of Common Motifs (ACoM) [35] in an attempt to give biological meaning to elementary modes and define relatedness between reactions. Application of elementary modes to rational metabolic engineering approaches has now enabled strain design algorithms such as those used by OptFlux [33] and the genetic algorithm-based method compiled by Boghigian et al. [36]. Elementary flux modes have been successfully used to engineer strains with a variety of desired phenotypes, including sugar coutilization [37], ethanol [37, 38], and carotenoid production [39]. Z. Chen et al., in this issue, use elementary modes to find strategies for improving the conversion of glycerol into succinate by considering the effects of oxygen utilization and genetic alterations (Z. Chen et al., “Elementary mode analysis for the rational design of efficient succinate conversion from glycerol by Escherichia coli”). Random sampling of flux distributions provides another way of investigating possible flux distributions through metabolic networks [40, 41], and A. De Martino et al. in this issue use both structural analysis and sampling to explore the robustness of human red blood cell (RBC) metabolism (A. De Martino et al., “Optimal fluxes, reaction replaceability, and response to enzymopathies in the human red blood cell”).

3.2. Optimization-Based Approaches

Alternatives to pathway-based approaches include optimization-based methods, which can also identify mutations that would improve desired phenotypes (e.g., increased production yields). Here, the models use the same stoichiometric and thermodynamic constraints as those used in pathway-based approaches, but solutions are identified which maximize or minimize a stated objective. To predict how metabolic fluxes will change in response to a genetic perturbation, a number of different approaches have been developed, including flux balance analysis (FBA, reviewed in [42]), minimization of metabolic adjustment (MOMA) [43], regulatory on/off mechanism (ROOM) [44], regulated flux balance analysis (rFBA and SR-FBA) [45, 46], and probabilistic regulation of metabolism (PROM) [47]. Most of these approaches are used to predict the immediate behavior of knockout strains [43, 44], with the exception being FBA, whose predictions more closely resemble strain behavior after cells have undergone adaptive evolution [48]. MOMA has been successfully used to engineer strains with increased production of a variety of products including lycopene [49, 50], valine [51], threonine [52], and polylactic acid [53].

To identify those genetic strategies that are predicted to have the highest chemical production levels, a large number of possible strategies must be considered. Bilevel approaches can be used to identify which are the best production strategies given a maximum number of genetic alterations without having to generate and store predictions for all possible strategies. These bilevel methods can be solved using integer programming and/or genetic algorithms [5457]. A variety of such bilevel methods have been proposed which use FBA, SR-FBA, and MOMA to predict mutant strain behaviors, and these approaches can consider genetic changes involving gene deletions, altered gene expression, or altered transcriptional regulation [54, 56, 5861]. These methods have been used to design strains for a variety of chemicals [55, 6264]. A more recent bilevel approach (OptFORCE) uses optimization to identify how metabolic fluxes must change to improve metabolite production, which is independent of any assumptions about what functions are used to predict cellular behavior [60].

3.3. Kinetic Modeling Approaches

Both types of models and methods described in the last two sections do not take enzyme kinetics into account. As a result, they are unable to predict how changes in kinetic properties, enzyme, and metabolite concentrations would affect fluxes through metabolic pathways. Kinetic models are needed to make these types of predictions since they capture the dependence of fluxes on metabolite and enzyme concentrations. These types of models can be analyzed to identify which changes are needed for improving cellular phenotypes. The classical framework for elucidating parameters responsible for the control of metabolic fluxes is metabolic control analysis (MCA), developed in the early 1970s independently by Kacser and Burns [65] and Heinrich and Rapoport [66]. Recently, Visser et al. developed an alternative approach called linlog kinetics [67, 68]. Here, all rate equations are modeled with the same basic mathematical structure in which the relationship between rates and enzyme levels is linear, while for metabolite levels, a linear combination of logarithmic terms is used. Young et al. [69] have recently expanded the cybernetic modeling framework of Ramkrishna [70] to incorporate metabolic pathway concepts derived from elementary mode analysis. This led to models that could predict both local and global control properties of metabolic networks in response to either dynamic environmental shifts or stable genetic manipulations. These models were applied to predict phenotypes of several recombinant E. coli strains, and they were found to provide good agreement with experimental data. Furthermore, because of the dynamic nature of these models, they were capable of simulating responses that are not readily addressed by purely stoichiometric models (e.g., allosteric or kinetic effects of intermediate metabolites, enzyme overexpressions and partial knockdowns, and time-dependent culture conditions). In this special issue, A. Yachie-Kinoshita et al. review the history of kinetic models for human red blood cells (RBCs) and describe an RBC metabolic model implemented in the E-cell simulation environment (A. Yachie-Kinoshita et al., “A metabolic model of human erythrocytes: practical application of the E-Cell Simulation Environment”). They discuss how this E-cell RBC model can be applied to predict RBC responses to hypoxic environments and long-term cold storage and identify enzymes whose altered activity could improve storage conditions for RBCs.

4. Concluding Remarks

A growing number of computational tools that facilitate the evaluation and improvement of strains for metabolic engineering are currently being developed and expanded. These methods can account for a wide range of experimental measurements to provide an improved understanding of metabolic states and current limitations, and they can be used to identify new engineering strategies for improved chemical production. The collection of papers in this special issue highlight several recent advances and underscore the emerging applications of these computational tools.

Jennifer L. Reed
Ryan S. Senger
Maciek R. Antoniewicz
Jamey D. Young


  1. U. Sauer, “Metabolic networks in motion: 13C-based flux analysis,” Molecular Systems Biology, vol. 2, article 62, 2006. View at: Publisher Site | Google Scholar
  2. J. Nielsen, “It is all about metabolic fluxes,” Journal of Bacteriology, vol. 185, no. 24, pp. 7031–7035, 2003. View at: Publisher Site | Google Scholar
  3. W. Wiechert, “C metabolic flux analysis,” Metabolic Engineering, vol. 3, no. 3, pp. 195–206, 2001. View at: Publisher Site | Google Scholar
  4. W. Wiechert, M. Möllney, N. Isermann, M. Wurzel, and A. A. De Graaf, “Bidirectional reaction steps in metabolic networks—III. Explicit solution and analysis of isotopomer labeling systems,” Biotechnology and Bioengineering, vol. 66, no. 2, pp. 69–85, 1999. View at: Publisher Site | Google Scholar
  5. J. F. Moxley, M. C. Jewett, M. R. Antoniewicz et al., “Linking high-resolution metabolic flux phenotypes and transcriptional regulation in yeast modulated by the global regulator Gcn4p,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 16, pp. 6477–6482, 2009. View at: Publisher Site | Google Scholar
  6. M. R. Antoniewicz, J. K. Kelleher, and G. Stephanopoulos, “Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions,” Metabolic Engineering, vol. 9, no. 1, pp. 68–86, 2007. View at: Publisher Site | Google Scholar
  7. P. F. Suthers, Y. J. Chang, and C. D. Maranas, “Improved computational performance of MFA using elementary metabolite units and flux coupling,” Metabolic Engineering, vol. 12, no. 2, pp. 123–128, 2010. View at: Publisher Site | Google Scholar
  8. M. R. Antoniewicz, D. F. Kraynie, L. A. Laffend, J. González-Lergier, J. K. Kelleher, and G. Stephanopoulos, “Metabolic flux analysis in a nonstationary systemml: fed-batch fermentation of a high yielding strain of E. coli producing 1,3-propanediol,” Metabolic Engineering, vol. 9, no. 3, pp. 277–292, 2007. View at: Publisher Site | Google Scholar
  9. Y. Noguchi, J. D. Young, J. O. Aleman, M. E. Hansen, J. K. Kelleher, and G. Stephanopoulos, “Effect of anaplerotic fluxes and amino acid availability on hepatic lipoapoptosis,” Journal of Biological Chemistry, vol. 284, no. 48, pp. 33425–33436, 2009. View at: Publisher Site | Google Scholar
  10. K. Nöh, K. Grönke, B. Luo, R. Takors, M. Oldiges, and W. Wiechert, “Metabolic flux analysis at ultra short time scale: isotopically non-stationary 13C labeling experiments,” Journal of Biotechnology, vol. 129, no. 2, pp. 249–267, 2007. View at: Publisher Site | Google Scholar
  11. S. A. Wahl, K. Nöh, and W. Wiechert, “13C labeling experiments at metabolic nonstationary conditions: an exploratory study,” BMC Bioinformatics, vol. 9, article 152, 2008. View at: Publisher Site | Google Scholar
  12. K. Schmidt, A. Marx, A. A. De Graaf et al., “13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches,” Biotechnology and Bioengineering, vol. 58, no. 2-3, pp. 254–257, 1998. View at: Publisher Site | Google Scholar
  13. J. D. Young, J. L. Walther, M. R. Antoniewicz, H. Yoo, and G. Stephanopoulos, “An elementary metabolite unit (EMU) based method of isotopically nonstationary flux analysis,” Biotechnology and Bioengineering, vol. 99, no. 3, pp. 686–699, 2008. View at: Publisher Site | Google Scholar
  14. Z. Zhao, K. Kuijvenhoven, C. Ras et al., “Isotopic non-stationary 13C gluconate tracer method for accurate determination of the pentose phosphate pathway split-ratio in Penicillium chrysogenum,” Metabolic Engineering, vol. 10, no. 3-4, pp. 178–186, 2008. View at: Publisher Site | Google Scholar
  15. K. Maier, U. Hofmann, M. Reuss, and K. Mauch, “Identification of metabolic fluxes in hepatic cells from transient C-labeling experiments—part II. Flux estimation,” Biotechnology and Bioengineering, vol. 100, no. 2, pp. 355–370, 2008. View at: Publisher Site | Google Scholar
  16. J. Munger, B. D. Bennett, A. Parikh et al., “Systems-level metabolic flux profiling identifies fatty acid synthesis as a target for antiviral therapy,” Nature Biotechnology, vol. 26, no. 10, pp. 1179–1186, 2008. View at: Publisher Site | Google Scholar
  17. K. Maier, U. Hofmann, A. Bauer et al., “Quantification of statin effects on hepatic cholesterol synthesis by transient 13C-flux analysis,” Metabolic Engineering, vol. 11, no. 4-5, pp. 292–309, 2009. View at: Publisher Site | Google Scholar
  18. S. Iwatani, Y. Yamada, and Y. Usuda, “Metabolic flux analysis in biotechnology processes,” Biotechnology Letters, vol. 30, no. 5, pp. 791–799, 2008. View at: Publisher Site | Google Scholar
  19. A. A. Shastri and J. A. Morgan, “A transient isotopic labeling methodology for C metabolic flux analysis of photoautotrophic microorganisms,” Phytochemistry, vol. 68, no. 16-18, pp. 2302–2312, 2007. View at: Publisher Site | Google Scholar
  20. W. Wiechert and K. Nöh, “From stationary to instationary metabolic flux analysis,” Advances in Biochemical Engineering/Biotechnology, vol. 92, pp. 145–172, 2005. View at: Google Scholar
  21. J. Schwender, “Metabolic flux analysis as a tool in metabolic engineering of plants,” Current Opinion in Biotechnology, vol. 19, no. 2, pp. 131–137, 2008. View at: Publisher Site | Google Scholar
  22. C. Yang, Q. Hua, and K. Shimizu, “Metabolic flux analysis in Synechocystis using isotope distribution from C-labeled glucose,” Metabolic Engineering, vol. 4, no. 3, pp. 202–216, 2002. View at: Publisher Site | Google Scholar
  23. B. D. Bennett, E. H. Kimball, M. Gao, R. Osterhout, S. J. Van Dien, and J. D. Rabinowitz, “Absolute metabolite concentrations and implied enzyme active site occupancy in Escherichia coli,” Nature Chemical Biology, vol. 5, no. 8, pp. 593–599, 2009. View at: Publisher Site | Google Scholar
  24. C. S. Henry, M. D. Jankowski, L. J. Broadbelt, and V. Hatzimanikatis, “Genome-scale thermodynamic analysis of Escherichia coli metabolism,” Biophysical Journal, vol. 90, no. 4, pp. 1453–1461, 2006. View at: Publisher Site | Google Scholar
  25. A. Kümmel, S. Panke, and M. Heinemann, “Putative regulatory sites unraveled by network-embedded thermodynamic analysis of metabolome data,” Molecular Systems Biology, vol. 2, article 2006.0034, 2006. View at: Google Scholar
  26. M. Klimacek, S. Krahulec, U. Sauer, and B. Nidetzky, “Limitations in xylose-fermenting saccharomyces cerevisiae, made evident through comprehensive metabolite profiling and thermodynamic analysis,” Applied and Environmental Microbiology, vol. 76, no. 22, pp. 7566–7574, 2010. View at: Publisher Site | Google Scholar
  27. S. A. Becker and B. O. Palsson, “Context-specific metabolic networks are consistent with experiments,” PLoS Computational Biology, vol. 4, no. 5, Article ID e1000082, 2008. View at: Publisher Site | Google Scholar
  28. C. Colijn, A. Brandes, J. Zucker et al., “Interpreting expression data with metabolic flux models: predicting Mycobacterium tuberculosis mycolic acid production,” PLoS Computational Biology, vol. 5, no. 8, article e1000489, 2009. View at: Publisher Site | Google Scholar
  29. T. Shlomi, M. N. Cabili, M. J. Herrgård, B. Ø. Palsson, and E. Ruppin, “Network-based prediction of human tissue-specific metabolism,” Nature Biotechnology, vol. 26, no. 9, pp. 1003–1010, 2008. View at: Publisher Site | Google Scholar
  30. A. von Kamp and S. Schuster, “Metatool 5.0: fast and flexible elementary modes analysis,” Bioinformatics, vol. 22, no. 15, pp. 1930–1931, 2006. View at: Publisher Site | Google Scholar
  31. S. Klamt, J. Gagneur, and A. von Kamp, “Algorithmic approaches for computing elementary modes in large biochemical reaction networks,” IEE Proceedings Systems Biology, vol. 152, no. 4, pp. 249–255, 2005. View at: Publisher Site | Google Scholar
  32. S. Klamt, J. Stelling, M. Ginkel, and E. D. Gilles, “FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps,” Bioinformatics, vol. 19, no. 2, pp. 261–269, 2003. View at: Publisher Site | Google Scholar
  33. I. Rocha, P. Maia, P. Evangelista et al., “OptFlux: an open-source software platform for in silico metabolic engineering,” BMC Systems Biology, vol. 4, article 45, 2010. View at: Publisher Site | Google Scholar
  34. S. Klamt and J. Stelling, “Combinatorial complexity of pathway analysis in metabolic networks,” Molecular Biology Reports, vol. 29, no. 1-2, pp. 233–236, 2002. View at: Publisher Site | Google Scholar
  35. S. Pérès, F. Vallée, M. Beurton-Aimar, and J. P. Mazat, “AComml: a classification method for elementary flux modes based on motif finding,” BioSystems, vol. 103, no. 3, pp. 410–419, 2011. View at: Publisher Site | Google Scholar
  36. B. A. Boghigian, H. Shi, K. Lee, and B. A. Pfeifer, “Utilizing elementary mode analysis, pathway thermodynamics, and a genetic algorithm for metabolic flux determination and optimal metabolic network design,” BMC Systems Biology, vol. 4, article 49, 2010. View at: Publisher Site | Google Scholar
  37. C. T. Trinh, P. Unrean, and F. Srienc, “Minimal Eschenchia coli cell for the most efficient production of ethanol from hexoses and pentoses,” Applied and Environmental Microbiology, vol. 74, no. 12, pp. 3634–3643, 2008. View at: Publisher Site | Google Scholar
  38. C. T. Trinh and F. Srienc, “Metabolic engineering of Escherichia coli for efficient conversion of glycerol to ethanol,” Applied and Environmental Microbiology, vol. 75, no. 21, pp. 6696–6705, 2009. View at: Publisher Site | Google Scholar
  39. P. Unrean, C. T. Trinh, and F. Srienc, “Rational design and construction of an efficient E. coli for production of diapolycopendioic acid,” Metabolic Engineering, vol. 12, no. 2, pp. 112–122, 2010. View at: Publisher Site | Google Scholar
  40. E. Almaas, B. Kovács, T. Vicsek, Z. N. Oltvai, and A. L. Barabási, “Global organization of metabolic fluxes in the bacterium Escherichia coli,” Nature, vol. 427, no. 6977, pp. 839–843, 2004. View at: Publisher Site | Google Scholar
  41. I. Thiele, N. D. Price, T. D. Vo, and B. Ø. Palsson, “Candidate metabolic network states in human mitochondria. Impact of diabetes, ischemia, and diet,” Journal of Biological Chemistry, vol. 280, no. 12, pp. 11683–11695, 2005. View at: Publisher Site | Google Scholar
  42. N. D. Price, J. L. Reed, and B. Ø. Palsson, “Genome-scale models of microbial cells: evaluating the consequences of constraints,” Nature Reviews Microbiology, vol. 2, no. 11, pp. 886–897, 2004. View at: Publisher Site | Google Scholar
  43. D. Segrè, D. Vitkup, and G. M. Church, “Analysis of optimality in natural and perturbed metabolic networks,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 23, pp. 15112–15117, 2002. View at: Publisher Site | Google Scholar
  44. T. Shlomi, O. Berkman, and E. Ruppin, “Regulatory on/off minimization of metabolic flux changes after genetic perturbations,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 21, pp. 7695–7700, 2005. View at: Publisher Site | Google Scholar
  45. M. W. Covert and B. Ø. Palsson, “Transcriptional regulation in constraints-based metabolic models of Escherichia coli,” Journal of Biological Chemistry, vol. 277, no. 31, pp. 28058–28064, 2002. View at: Publisher Site | Google Scholar
  46. T. Shlomi, Y. Eisenberg, R. Sharan, and E. Ruppin, “A genome-scale computational study of the interplay between transcriptional regulation and metabolism,” Molecular Systems Biology, vol. 3, article 101, 2007. View at: Publisher Site | Google Scholar
  47. S. Chandrasekaran and N. D. Price, “Probabilistic integrative modeling of genome-scale metabolic and regulatory networks in Escherichia coli and Mycobacterium tuberculosis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 41, pp. 17845–17850, 2010. View at: Publisher Site | Google Scholar
  48. S. S. Fong and B. Ø. Palsson, “Metabolic gene-deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes,” Nature Genetics, vol. 36, no. 10, pp. 1056–1058, 2004. View at: Publisher Site | Google Scholar
  49. H. Alper, Y. S. Jin, J. F. Moxley, and G. Stephanopoulos, “Identifying gene targets for the metabolic engineering of lycopene biosynthesis in Escherichia coli,” Metabolic Engineering, vol. 7, no. 3, pp. 155–164, 2005. View at: Publisher Site | Google Scholar
  50. H. Alper, K. Miyaoku, and G. Stephanopoulos, “Construction of lycopene-overproducing E. coli strains by combining systematic and combinatorial gene knockout targets,” Nature Biotechnology, vol. 23, no. 5, pp. 612–616, 2005. View at: Publisher Site | Google Scholar
  51. J. H. Park, K. H. Lee, T. Y. Kim, and S. Y. Lee, “Metabolic engineering of Escherichia coli for the production of L-valine based on transcriptome analysis and in silico gene knockout simulation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 19, pp. 7797–7802, 2007. View at: Publisher Site | Google Scholar
  52. K. H. Lee, J. H. Park, T. Y. Kim, H. U. Kim, and S. Y. Lee, “Systems metabolic engineering of Escherichia coli for L-threonine production,” Molecular Systems Biology, vol. 3, article 149, 2007. View at: Publisher Site | Google Scholar
  53. Y. K. Jung, T. Y. Kim, S. J. Park, and S. Y. Lee, “Metabolic engineering of Escherichia coli for the production of polylactic acid and its copolymers,” Biotechnology and Bioengineering, vol. 105, no. 1, pp. 161–171, 2010. View at: Publisher Site | Google Scholar
  54. A. P. Burgard, P. Pharkya, and C. D. Maranas, “OptKnock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization,” Biotechnology and Bioengineering, vol. 84, no. 6, pp. 647–657, 2003. View at: Publisher Site | Google Scholar
  55. J. A. Chemler, Z. L. Fowler, K. P. McHugh, and M. A. G. Koffas, “Improving NADPH availability for natural product biosynthesis in Escherichia coli by metabolic engineering,” Metabolic Engineering, vol. 12, no. 2, pp. 96–104, 2010. View at: Publisher Site | Google Scholar
  56. K. R. Patil, I. Rocha, J. Förster, and J. Nielsen, “Evolutionary programming as a platform for in silico metabolic engineering,” BMC Bioinformatics, vol. 6, article 308, 2005. View at: Publisher Site | Google Scholar
  57. D. S. Lun, G. Rockwell, N. J. Guido et al., “Large-scale identification of genetic design strategies using local search,” Molecular Systems Biology, vol. 5, article 296, 2009. View at: Publisher Site | Google Scholar
  58. P. Pharkya, A. P. Burgard, and C. D. Maranas, “OptStrain: a computational framework for redesign of microbial production systems,” Genome Research, vol. 14, no. 11, pp. 2367–2376, 2004. View at: Publisher Site | Google Scholar
  59. P. Pharkya and C. D. Maranas, “An optimization framework for identifying reaction activation/inhibition or elimination candidates for overproduction in microbial systems,” Metabolic Engineering, vol. 8, no. 1, pp. 1–13, 2006. View at: Publisher Site | Google Scholar
  60. S. Ranganathan, P. F. Suthers, and C. D. Maranas, “OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions,” PLoS Computational Biology, vol. 6, no. 4, article e1000744, 2010. View at: Publisher Site | Google Scholar
  61. J. Kim and J. L. Reed, “OptORF: optimal metabolic and regulatory perturbations for metabolic engineering of microbial strains,” BMC Systems Biology, vol. 4, article 53, 2010. View at: Publisher Site | Google Scholar
  62. M. A. Asadollahi, J. Maury, K. R. Patil, M. Schalk, A. Clark, and J. Nielsen, “Enhancing sesquiterpene production in Saccharomyces cerevisiae through in silico driven metabolic engineering,” Metabolic Engineering, vol. 11, no. 6, pp. 328–334, 2009. View at: Publisher Site | Google Scholar
  63. S. S. Fong, A. P. Burgard, C. D. Herring et al., “In silico design and adaptive evolution of Escherichia coli for production of lactic acid,” Biotechnology and Bioengineering, vol. 91, no. 5, pp. 643–648, 2005. View at: Publisher Site | Google Scholar
  64. A. R. Brochado, C. Matos, B. L. Møller, J. Hansen, U. H. Mortensen, and K. R. Patil, “Improved vanillin production in baker's yeast through in silico design,” Microbial Cell Factories, vol. 9, article 84, 2010. View at: Publisher Site | Google Scholar
  65. H. Kacser and J. A. Burns, “The control of flux,” Symposia of the Society for Experimental Biology, vol. 27, pp. 65–104, 1973. View at: Google Scholar
  66. R. Heinrich and T. A. Rapoport, “A linear steady-state treatment of enzymatic chains. General properties, control and effector strength,” European Journal of Biochemistry, vol. 42, no. 1, pp. 89–95, 1974. View at: Google Scholar
  67. D. Visser and J. J. Heijnen, “Dynamic simulation and metabolic re-design of a branched pathway using linlog kinetics,” Metabolic Engineering, vol. 5, no. 3, pp. 164–176, 2003. View at: Google Scholar
  68. D. Visser, J. W. Schmid, K. Mauch, M. Reuss, and J. J. Heijnen, “Optimal re-design of primary metabolism in Escherichia coli using linlog kinetics,” Metabolic Engineering, vol. 6, no. 4, pp. 378–390, 2004. View at: Google Scholar
  69. J. D. Young, K. L. Henne, J. A. Morgan, A. E. Konopka, and D. Ramkrishna, “Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control,” Biotechnology and Bioengineering, vol. 100, no. 3, pp. 542–559, 2008. View at: Publisher Site | Google Scholar
  70. D. Ramkrishna, “A cybernetic perspective of microbial growth,” in Foundations of Biochemical Engineering: Kinetics and Thermodynamics in Biological Systems, pp. 161–178, American Chemical Society, Washington, DC, USA, 1982. View at: Google Scholar

Copyright © 2010 Jennifer L. Reed et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1594 Views | 884 Downloads | 15 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19.