- Open Access
Metabolic control analysis: biological applications and insights
Genome Biology volume 1, Article number: reviews1031.1 (2000)
Metabolic control analysis provides a robust mathematical and theoretical framework for describing metabolic and signaling pathways and networks, and for quantifying the controls over these processes. Its application has already shed light on some of the principles underlying the regulation of metabolic pathways, and it is well suited to the analysis of the types of data emerging from genomic studies.
Modeling biological processes and systems
Mathematical modeling allows us to examine an event, process or system that we are unable to observe or understand directly because of its timing, magnitude, location or complexity. Models enable us to view a process or system at different organizational levels (for example, molecular or organismal) 'simultaneously', and to test responses of the system and its components to perturbations. Even incomplete or limited models can pinpoint missing or incorrect pathways or components and can help ascertain the relative importance of pathways and components in different scenarios. Models can also elucidate underlying biological design principles, sometimes challenging existing scientific dogma. They do this by extending and integrating the effects of assumptions made at one organizational level to others, and by allowing the 'visualization' of hypothetical scenarios. The term model is used broadly for the purposes of this article, and is defined as any mathematical or theoretical framework used to describe a component, process or system. This description can take many forms and make use of a variety of mathematical techniques (reviewed in ).
To describe biological systems, which are naturally complex and integrated, properties at different organizational levels must be related to each other in a meaningful way. Thus, the properties of a system ('systems properties') must reflect underlying molecular design principles, and equations detailing molecular components must take into account systems-level constraints and contexts. One such constraint is capacity, the maximum allowable flux. The system, as an entity, is not just an assemblage of its individual parts but has 'emergent' properties of the whole . For a simple example, imagine a rubber ball that is cut into many small pieces. Much can be learned about the properties of the ball based on the individual cut pieces (for example, its elasticity), but we wouldn't know that the ball could roll.
Various modeling approaches have been used to study cellular metabolism and signaling (for reviews, see [3,4,5,6,7]). For the analysis of metabolism, such approaches include kinetic simulation, metabolic control analysis , biochemical systems theory [9,10], metabolic pathway analysis [11,12], and network analysis [6,13]. Signal transduction pathways and networks have, for the most part, been described qualitatively by the sets of expressed genes associated with the activation of a specific pathway [14,15]. Signaling has been modeled, however, for certain well-defined systems using detailed kinetics and neural-network-type approaches [2,16]. In addition, Krauss and Brand  have recently applied metabolic control analysis to signal transduction pathways. A description of the specific focuses and features of the above methods is not possible here. However, they can generally be grouped according to the organizational level(s) they describe and the type of data utilized. This article focuses on metabolic control analysis, a method that integrates 'local' kinetic information with systems-level information to quantitate proportions of control exerted by different components of a given pathway or system. The aim of this review is to provide the reader with the basic framework for understanding how (and why) metabolic control analysis can be used to examine specific systems and to elucidate fundamental underlying biological design principles, which are independent of a particular system.
Metabolic control analysis
Metabolic control analysis provides a robust mathematical and theoretical framework for describing metabolic and signaling pathways and networks, and for quantifying the controls over these processes. It can deal with systems of any complexity or architecture and does not require all system components to be known a priori, making it a valuable post-genomic tool. It was developed in the 1970s by Kacser and Burns  and Heinrich and Rapoport . Since then, dedicated researchers have expanded and advanced metabolic control analysis theory and applications, carefully defined the associated terms, and developed analytical and educational tools [20,21,22,23,24,25].
Metabolic control analysis uses equations based on the kinetics of enzymes (known as elasticity coefficients) to parameterize control coefficients (which describe the degree of control exerted by any given component on a particular output) resulting in response coefficients. The response to a perturbation is quantified by the summation of the response coefficients affecting the output of interest, which depends on both the control coefficients (systems properties) and elasticity coefficients (local properties). The fundamental equations are presented in Box 1 (also see Figure 1). It should be noted, explicitly, that the use of metabolic control analysis is not limited to linear pathways, but is also applicable to branching and cyclic pathways, and enzyme cascades.
Application to a specific system
Metabolic control analysis can be used to quantitate control exerted by different components of a specific system and pinpoint areas requiring further experimentation. For example, Krauss and Brand  recently used metabolic control analysis to quantitate the contributions of known and unknown signal transduction pathways in the early response of thymocytes to mitogen (concanavalin A, ConA) stimulation. To apply metabolic control analysis to a complex system containing both metabolic and signaling components, they gathered large parts of metabolism or signal transduction into 'black-box' groups of reactions, coupled signal transduction events to cellular variables, and limited the time frame of observation. In particular, thymocyte response to Con A stimulation was quantified by measuring steady-state respiration rates, and signaling routes (such as protein kinase C (PKC)) were grouped based on their sensitivity to specific inhibitors. The effects of these known and unknown signal transduction pathways on the mitochondrial membrane potential, a key intermediate in respiration, were also quantified. The analysis of the model system, presented in Figure 2, resulted in a quantitative topology of signaling routes involved in the early phase of mitogen stimulation of thymocytes. Novel findings, such as the significant role played by calcineurin signal transduction pathways (30% of total), highlight areas for future experimental work.
Elucidation of underlying biological design principles
The application of metabolic control analysis has altered our basic understanding of metabolic control. In particular, the belief that control over a pathway is dictated by a 'rate-limiting step' is now obsolete and is being removed from the biochemistry textbooks. Instead, it is replaced by the concept of shared control, where many - or theoretically all - enzymes in a pathway have a role in controlling the flux through the pathway. Implicit in this is the idea that the regulation of cell metabolism requires coordinate change in the activities of many enzymes ('multisite modulation' ). The validity of this notion has been supported by bioengineers' lack of success in increasing a particular flux (product yield) by overexpressing the 'rate-limiting' enzyme and success by overexpressing a group of enzymes in a pathway. For example, Niederberger et al.  found that overexpression of four of the five enzymes in the yeast biosynthetic pathway leading from chorismate to tryptophan was required to significantly increase (more than eightfold) the production of tryptophan.
If the coordinated expression of enzymes in a pathway is required to significantly increase its flux, this should be an underlying design principle of organisms. In fact, the coordinated induction of enzymes to increase metabolic flux through a pathway has long been observed in vivo. One of many examples provided by Fell  is the urea synthetic pathway. The rate of urea synthesis in rats responds proportionately to the amount of protein in the diet. When rats are fed on diets that increase urea output fourfold, eight of the enzymes measured increased significantly, including all four of the urea-cycle enzymes . With DNA microarrays and complete genome information, global expression data detailing the coordinated induction of pathway enzymes may be coupled with structural information on the organization of genes for pathway enzymes in operons or in clusters with common cis-acting elements (for example, see ). By performing these types of analyses on organisms responding to a variety of external effectors (such as nutrient conditions, pathogens, and so on), and on diverse organisms, this underlying design principle may be further explored.
Metabolic control analysis can also be used to explain why most mutations in diploid organisms are 'fully' recessive. Most enzymes have low-flux control coefficients; thus, a 50% reduction in enzyme concentration resulting from a null mutation in one allele of a diploid pair has little effect on the pathway flux. In addition, because pathway flux is a systems property, the influence of an alteration at one locus is measured in the whole system, minimizing the impact from any one reduction. Kacser and Burns  therefore posited the phenomenon of genetic dominance as the "inevitable consequence of the kinetic structure of enzyme networks" and not a result of natural selection. This conclusion was supported by Orr , who found the same extent of recessive mutations in artificial diploids created from the haploid organism Chlamydomonas reinhardtii (where the possibility of selection in the diploid was eliminated). The existence of a limited number of 'partially' recessive mutants, in which the heterozygote has an intermediate phenotype, is also consistent with metabolic control analysis. Theoretically, these enzymes (with high control coefficients) would be more likely to be a part of a very small pathway or the first enzyme of a branching pathway. Despite these studies, the 'inevitability' of dominance is still debated . With the availability of complete genome information for a number of diploid organisms, genomic information on natural variants (for example, the Arabidopsis ecotypes Columbia and Landsberg ), and numerous collections of mutants, this type of question can now be addressed on a global scale.
Future directions and challenges
As illustrated above, metabolic control analysis is particularly useful for describing the theoretical aspects of regulation. This utility will continue to expand in the post-genomic era, particularly with advances in the in vivo imaging and quantitation of proteins and metabolites (for example, using tracer nuclear magnetic resonance). Future modeling efforts will require the integration or sampling of current mathematical approaches, including metabolic control analysis, as well as the development of new theoretical approaches and tools. As models become more complex and integrated to reflect the sheer volume of simultaneously occurring reactions in a cell, the incorporation of Monte Carlo methods (random sampling) and finite element analysis (approximations based on subdivision into smaller, more manageable elements) is likely to be necessary. In addition, the platforms and databases required to construct models of increasing complexity need to be developed in an organized and collaborative manner and to be widely accessible . Access to the requisite computational resources will also become an issue. Perhaps an institute similar to the National Center for Atmospheric Research , which facilitates international global climate change research, could help coordinate and support biological modeling efforts.
Gershenfeld N: The Nature of Mathematical Modeling. Cambridge University Press;. 1999
Bhalla US, Iyengar R: Emergent properties of networks of biological signaling pathways. Science. 1999, 283: 381-387. 10.1016/S0009-2614(97)01391-2.
Collado-Vides J, Magasanik B, Smith TF: Integrative Approaches to Molecular Biology. Cambridge, MA: MIT Press;. 1996
Giersch C: Mathematical modeling of metabolism. Curr Opin Plant Biol. 2000, 3: 249-253. 10.1016/S1369-5266(00)00072-8.
Palsson B: The challenges of in silico biology. Nat Biotechnol. 2000, 18: 1147-1150. 10.1038/81125.
Fell DA, Wagner A: The small world of metabolism. Nat Biotechnol. 2000, 18: 1121-1122. 10.1038/81025.
Weng G, Bhalla US, Iyengar R: Complexity in biological signaling systems. Science. 1999, 284: 92-96. 10.1016/S0925-8388(98)00740-3.
Poolman MG, Fell DA, Thomas S: Modelling photosynthesis and its control. J Exp Bot. 2000, 51: 319-328. 10.1093/jexbot/51.suppl_1.319.
Ni T-C, Savageau M: Application of biochemical systems theory to metabolism in human red blood cells. J Biol Chem. 1996, 271: 7927-7941. 10.1074/jbc.271.14.7927.
Savageau MA: Power-law formalism: a canonical nonlinear approach to modeling and analysis. In World Congress of Nonlinear Analysts, 92, Vol 4. Edited by Lakshmikantham V. Berlin: Walter de Gruyter;. 1996
Schilling CH, Schuster S, Palsson BO, Heinrich R: Metabolic pathway analysis: basic concepts and scientific applications in the post-genomic era. Biotechnol Prog. 1999, 15: 296-303. 10.1021/bp990048k.
Edwards JS, Palsson BO: Systems properties of the Haemophilus influenzae Rd metabolic genotype. J Biol Chem. 1999, 274: 17410-17416. 10.1074/jbc.274.25.17410.
Jeong H., Tombor B, Albert R, Oltavai N, Barabasi A-L: The large-scale organization of metabolic networks. Nature. 2000, 407: 651-654. 10.1038/35036627.
Fambrough D, McClure K, Kazlauskas A, Lander ES: Diverse signaling pathways activated by growth factor receptors induce broadly overlapping, rather than independent, sets of genes. Cell. 1999, 97: 727-741.
Pawson T, Saxton TM: Signaling networks - do all roads lead to the same genes?. Cell. 1999, 97: 675-678.
Kholodenko BN, Demin OV, Moehren G, Hoek JB: Quantification of short term signaling by the epidermal growth factor receptor. J Biol Chem. 1999, 274: 30169-30181. 10.1074/jbc.274.42.30169.
Krauss S, Brand MD: Quantitation of signal transduction. FASEB J. 2000, 14: 2581-2588. 10.1096/fj.00-0064com.
Kacser H, Burns JA: Control of enzyme flux. Symp Soc Exp Biol. 1973, 27: 65-104.
Heinrich R, Rapoport TA: A linear steady-state treatment of enzymatic chains. Eur J Biochem. 1974, 42: 89-95.
Fell DA: Metabolic control analysis - a survey of its theoretical and experimental development. Biochem J. 1992, 286: 313-330.
Fell D: Understanding the Control of Metabolism. London: Portland Press;. 1997
Heinrich R, Schuster S: The Regulation of Cellular Systems. New York: Chapman and Hall;. 1996
Cornish-Bowden A: Metabolic control analysis in theory and practice. Adv Mol Cell Biol. 1995, 11: 21-64.
MCA website. [http://gepasi.dbs.aber.ac.uk/metab/mca_home.htm]
Bionet metabolic regulation newsgroup. [http://www.ibiblio.org/usenet-i/groups-html/bionet.metabolic-reg.html]
Fell DA, Thomas S: Physiological control of metabolic flux: the requirement for multisite modulation. Biochem J. 1995, 311: 35-39.
Niederberger P, Prasad R, Miozzari G, Kacser H: A strategy for increasing an in vivo flux by genetic manipulations of the tryptophan system of yeast. Biochem J. 1992, 287: 473-479.
Schimke RT: Adaptive characteristics of urea cycle enzymes in the rat. J Biol Chem. 1962, 237: 459-468.
Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM: Systematic determination of genetic network architecture. Nat Genet. 1999, 22: 281-285. 10.1038/10343.
Kacser H, Burns JA: The molecular basis of dominance. Genetics. 1981, 97: 639-666.
Orr HA: A test of Fisher's theory of dominance. Proc Natl Acad Sci USA. 1991, 88: 11413-11415.
Grossniklaus U, Madhusudhan MS, Nanjundiah V: Nonlinear enzyme kinetics can lead to high metabolic flux control coefficients: implications for the evolution of dominance. J Theor Biol. 1996, 182: 299-302. 10.1006/jtbi.1996.0167.
The Arabidopsis Information Resource. [http://www.arabidopsis.org]
Alliance for Cellular Signaling. [http://afcs.swmed.edu/]
National Center for Atmospheric Research. [http://www.ncar.ucar.edu/ncar/index.html]
My sincere thanks to David Fell, Stefan Krauss, Fred Ausubel and Julia Dewdney for their comments on drafts of this manuscript.
About this article
Cite this article
Wildermuth, M.C. Metabolic control analysis: biological applications and insights. Genome Biol 1, reviews1031.1 (2000) doi:10.1186/gb-2000-1-6-reviews1031
- Response Coefficient
- Control Coefficient
- Pathway Flux
- Metabolic Control Analysis
- Diploid Organism