Volume 1 Supplement 1

BioSysBio 2007: Systems Biology, Bioinformatics, Synthetic Biology

Open Access

Kinetic constraints on the sensitivity of large metabolic responses

BMC Systems Biology20071(Suppl 1):P46

DOI: 10.1186/1752-0509-1-S1-P46

Published: 8 May 2007


The quantitative study of metabolic responses in intact cells is essential for predicting the phenotypic consequences of genetic manipulations. Metabolic responses have been described within the framework of metabolic control analysis (MCA) [1]. One of the central goals of MCA is to determine how the responses of system variables, quantified by control coefficients, depend on the kinetic properties of the component reactions or groups of reactions (modules). Attempts to extend infinitesimal control analysis to large changes in the variables have been reviewed in previous publications [2, 3].

In the present contribution, we analyse how some kinetic features of the modules constraint the sensitivity of large metabolic responses. In particular, we explore the different kinds of flux control profile that can be achieved using modules that follow two different types of rate equations: Michaelis-Menten (M-M) and Hill.

Materials and methods

The mean control coefficient quantifies the sensitivity of a steady-state metabolic variable (W) to changes in a parameter (p). It is defined as C p W ¯ = ( ( W x W o ) / W o ) / ( ( p x p o ) / p o ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadGaG31aaaeacaENama4DdoeadnacaExhaaWcbGaG3jadaE3GWbaCaeacaENama4DdEfaxbaaaaGccGaG3HPaRladaENH9aqpcGaG3HPaRpacaE3cgaqaia4DdGaG3fWaaeacaE3aia4DlyaabGaG3nacaExadaqaia4DcWaG3n4vaC1aia4DBaaaleacaENama4DdIha4bqaja4DaOGaia4DykW6cWaG3zOeI0Iaia4DykW6cWaG3n4vaC1aia4DBaaaleacaENama4Dd+gaVbqaja4DaaGccGaG3zjkaiacaENLPaaaaeacaENama4DdEfaxnacaE3gaaWcbGaG3jadaE3GVbWBaeqcaEhaaaaakiacaENLOaGaia4DwMcaaaqaia4DdGaG3fWaaeacaE3aia4DlyaabGaG3nacaExadaqaia4DcWaG3niCaa3aia4DBaaaleacaENama4DdIha4bqaja4DaOGaia4DykW6cWaG3zOeI0Iaia4DykW6cWaG3niCaa3aia4DBaaaleacaENama4Dd+gaVbqaja4DaaGccGaG3zjkaiacaENLPaaaaeacaENama4DdchaWnacaE3gaaWcbGaG3jadaE3GVbWBaeqcaEhaaaaakiacaENLOaGaia4DwMcaaaaaaaa@9C7D@ . It represents the relative change in W divided by the relative change in p when the system goes from one initial (o) to one final steady state (x) [2].


In this study, we consider that the metabolic network can be split into a supply module (v1), that produces the intermediate S, and a demand module (v2), that consumes S (Figure 1) [4]. This top-down or modular MCA approach has been widely used to analyse the sensitivity properties of metabolic networks [5]. Firstly, the rates and vs. S were chosen to obey a reversible M-M and Hill equation respectively (see Figure 2 for details). The maximal activities of module 1 and 2 are modified by the factor r1 and r2, respectively (one at a time). The system starts at a reference steady-state o, where r1 = r2 = 1. After the parameter is changed (r1 or r2) the variables of the system (S and J) freely adjust to the final steady state. In Figure 2 we plot the mean flux and metabolite control coefficient for both modules when subject to large changes in the enzyme concentrations or activities. The infinitesimal control coefficients (r1 = r2 = r = 1) are C v 1 J = 0.58 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacWaGWo4qam0aiaiSDaaaleacacRamaiSdAha2jadaclIXaqmaeacacRamaiSdQeakbaakiacaclMcSUamaiSg2da9iacaclMcSUamaiSicdaWiadacRGUaGlcWaGWIynauJamaiSiIda4aaa@4A60@ , C v 2 J = 0.42 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacWaGyo4qam0aiaiMDaaaleacaIPamaiMdAha2jadaIjIYaGmaeacaIPamaiMdQeakbaakiacaIjMcSUamaiMg2da9iacaIjMcSUamaiMicdaWiadaIPGUaGlcWaGyIinaqJamaiMikdaYaaa@4904@ , and C v 1 S = | C v 2 S | = 1.05 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacWaGyo4qam0aiaiMDaaaleacaIPamaiMdAha2jadaIjIXaqmaeacaIPamaiMdofatbaakiacaIjMcSUamaiMg2da9iacaIjMcS+aiaiMemaabGaGykadaI5GdbWqdGaGy2baaSqaiaiMcWaGyoODayNamaiMikdaYaqaiaiMcWaGyo4uamfaaaGccGaGyA5bSlacaIPLiWoacGaGyIPaRladaIPH9aqpcGaGyIPaRladaIjIXaqmcWaGykOla4IamaiMicdaWiadaIjI1aqnaaa@64E6@ . But for large changes, which module shows the greatest control depends on the value of r. In the interval (0.56–1.8) the greatest flux control is in module 1, meanwhile outside this interval the flux is mainly controlled by the second module. Regarding the control of S, for r > 1a change in module 1 produces a greater effect on S than module 2, while for r < 1 module 2 has the greatest effect. On the other hand, if both modules fulfil M-M rate equation it can be demonstrated that for all values of r one of the modules has the greatest flux control (results not shown).
Figure 1

Metabolic system constituted by a supply module (v1) and a demand module (v2) kinked by one intermediate.

Figure 2

Mean control coefficients vs module activity, r. a) flux mean control coefficients and b) intermediate mean control coefficients. M-M equation parameters: X0 = 10, KS = 1, KP = 0.15, 1/Keq = 0 and Vm1 = 1. Hill equation parameters: KS = 0.15, h = 4 and Vm2 = 0.63.


We have analysed two different models to see if certain patterns of control for large changes could be, in principle, obtained using modules whose rates are governed by usual rate equations. We have shown that, if the rates of both modules obey hyperbolic kinetics, the flux control distribution is constrained, one module having the greatest control over all enzyme activity range. On the contrary, if the rate of the demand module follows a Hill kinetics, which module has the greatest control on the flux depends on the enzymatic activity change.

Authors’ Affiliations

School of Biosciences, University of Birmingham
Laboratorio de Biología de Sistemas, Universidad de la República


  1. Heinrich R, Schuster S: The regulation of cellular systems. 1996, Chapman & Hall, NY, 1View ArticleGoogle Scholar
  2. Ortega F, Acerenza L: Elasticity analysis and design for large metabolic responses produced by changes in enzyme activities. Biochem J. 2002, 367: 41-48. 10.1042/BJ20020520PubMedPubMed CentralView ArticleGoogle Scholar
  3. Acerenza L, Ortega F: Metabolic control analysis for large changes: extension to variable elasticity coefficients. IEE Proc Syst Biol. 2006, 153: 323-326. 10.1049/ip-syb:20060004.View ArticleGoogle Scholar
  4. Hofmeyr JS, Cornish-Bowden A: Regulating the cellular economy of supply and demand. FEBS Lett. 2000, 476: 47-51. 10.1016/S0014-5793(00)01668-9PubMedView ArticleGoogle Scholar
  5. Quant PA: Experimental application of top-down control analysis to metabolic systems. Trends Biochem Sci. 1993, 18: 26-30. 10.1016/0968-0004(93)90084-ZPubMedView ArticleGoogle Scholar


© Ortega and Acerenza; licensee BioMed Central Ltd. 2007

This article is published under license to BioMed Central Ltd.