Metabolic control analysis

Last updated
Plot of steady state flux versus enzyme activity with flux control coefficients at various points. MCAlogo.png
Plot of steady state flux versus enzyme activity with flux control coefficients at various points.

Metabolic control analysis (MCA) is a mathematical framework for describing metabolic, signaling, and genetic pathways. MCA quantifies how variables, such as fluxes and species concentrations, depend on network parameters. In particular, it is able to describe how network-dependent properties, called control coefficients, depend on local properties called elasticities or Elasticity Coefficients. [1] [2] [3]

Contents

MCA was originally developed to describe the control in metabolic pathways but was subsequently extended to describe signaling and genetic networks. MCA has sometimes also been referred to as Metabolic Control Theory, but this terminology was rather strongly opposed by Henrik Kacser, one of the founders[ citation needed ].

More recent work [4] has shown that MCA can be mapped directly on to classical control theory and are as such equivalent.

Biochemical systems theory [5] (BST) is a similar formalism, though with rather different objectives. Both are evolutions of an earlier theoretical analysis by Joseph Higgins. [6]

Chemical reaction network theory is another theoretical framework that has overlap with both MCA and BST but is considerably more mathematically formal in its approach. [7] It's emphasis is primarily on dynamic stability criteria [8] and related theorems associated with mass-action networks. In more recent years the field has also developed [9] a sensitivity analysis which is similar if not identical to MCA and BST.

Control coefficients

A control coefficient [10] [11] [12] measures the relative steady state change in a system variable, e.g. pathway flux (J) or metabolite concentration (S), in response to a relative change in a parameter, e.g. enzyme activity or the steady-state rate () of step . The two main control coefficients are the flux and concentration control coefficients. Flux control coefficients are defined by

and concentration control coefficients by

Control coefficients measure the effect of perturbations in an enzyme on a steady-state observable such as a flux or metabolite concentration. Note that the effect of a perturbation can be positive or negative depending on context. In the figure, a perturbation is assumed to be at step three. Figure modified and redrawn from SystemControl.png
Control coefficients measure the effect of perturbations in an enzyme on a steady-state observable such as a flux or metabolite concentration. Note that the effect of a perturbation can be positive or negative depending on context. In the figure, a perturbation is assumed to be at step three. Figure modified and redrawn from

.

Summation theorems

The flux control summation theorem was discovered independently by the Kacser/Burns group [10] and the Heinrich/Rapoport group [11] in the early 1970s and late 1960s. The flux control summation theorem implies that metabolic fluxes are systemic properties and that their control is shared by all reactions in the system. When a single reaction changes its control of the flux this is compensated by changes in the control of the same flux by all other reactions.

Elasticity coefficients

The elasticity coefficient measures the local response of an enzyme or other chemical reaction to changes in its environment. Such changes include factors such as substrates, products, or effector concentrations. For further information, please refer to the dedicated page at elasticity coefficients.

Elasticities are local quantities that measure the effect of substrates, products, and effectors on a given reaction. The vertical black arrows represent enzyme catalysis. Figure redrawn and modified from LocalControl.png
Elasticities are local quantities that measure the effect of substrates, products, and effectors on a given reaction. The vertical black arrows represent enzyme catalysis. Figure redrawn and modified from

.

Connectivity theorems

The connectivity theorems [10] [11] are specific relationships between elasticities and control coefficients. They are useful because they highlight the close relationship between the kinetic properties of individual reactions and the system properties of a pathway. Two basic sets of theorems exists, one for flux and another for concentrations. The concentration connectivity theorems are divided again depending on whether the system species is different from the local species .

Response Coefficient

Kacser and Burns [10] introduced an additional coefficient that described how a biochemical pathway would respond the external environment. They termed this coefficient the response coefficient and designated it using the symbol R. The response coefficient is an important metric because it can be used to assess how much a nutrient or perhaps more important, how a drug can influence a pathway. This coefficient is therefore highly relevant to the pharmaceutical industry. [15]

The response coefficient is related to the core of metabolic control analysis via the response coefficient theorem, which is stated as follows:

where is a chosen observable such as a flux or metabolite concentration, is the step that the external factor targets, is the control coefficient of the target steps, and is the elasticity of the target step with respect to the external factor .

The key observation of this theorem is that an external factor such as a therapeutic drug, acts on the organism's phenotype via two influences: 1) How well the drug can affect the target itself through effective binding of the drug to the target protein and its effect on the protein activity. This effectiveness is described by the elasticity and 2) How well do modifications of the target influence the phenotype by transmission of the perturbation to the rest of the network. This is indicated by the control coefficient .

A drug action, or any external factor, is most effective when both these factors are strong. For example, a drug might be very effective at changing the activity of its target protein, however if that perturbation in protein activity is unable to be transmitted to the final phenotype then the effectiveness of the drug is greatly diminished.

If a drug or external factor, , targets multiple sites of action, for example sites, then the overall response in a phenotypic factor , is the sum of the individual responses:

Control equations

It is possible to combine the summation with the connectivity theorems to obtain closed expressions that relate the control coefficients to the elasticity coefficients. For example, consider the simplest non-trivial pathway:

We assume that and are fixed boundary species so that the pathway can reach a steady state. Let the first step have a rate and the second step . Focusing on the flux control coefficients, we can write one summation and one connectivity theorem for this simple pathway:

Using these two equations we can solve for the flux control coefficients to yield

Using these equations we can look at some simple extreme behaviors. For example, let us assume that the first step is completely insensitive to its product (i.e. not reacting with it), S, then . In this case, the control coefficients reduce to

That is all the control (or sensitivity) is on the first step. This situation represents the classic rate-limiting step that is frequently mentioned in textbooks. The flux through the pathway is completely dependent on the first step. Under these conditions, no other step in the pathway can affect the flux. The effect is however dependent on the complete insensitivity of the first step to its product. Such a situation is likely to be rare in real pathways. In fact the classic rate limiting step has almost never been observed experimentally. Instead, a range of limitingness is observed, with some steps having more limitingness (control) than others.

We can also derive the concentration control coefficients for the simple two step pathway:

Three step pathway

Consider the simple three step pathway:

where and are fixed boundary species, the control equations for this pathway can be derived in a similar manner to the simple two step pathway although it is somewhat more tedious.

where D the denominator is given by

Note that every term in the numerator appears in the denominator, this ensures that the flux control coefficient summation theorem is satisfied.

Likewise the concentration control coefficients can also be derived, for

And for

Note that the denominators remain the same as before and behave as a normalizing factor.

Derivation using perturbations

Control equations can also be derived by considering the effect of perturbations on the system. Consider that reaction rates and are determined by two enzymes and respectively. Changing either enzyme will result in a change to the steady state level of and the steady state reaction rates . Consider a small change in of magnitude . This will have a number of effects, it will increase which in turn will increase which in turn will increase . Eventually the system will settle to a new steady state. We can describe these changes by focusing on the change in and . The change in , which we designate , came about as a result of the change . Because we are only considering small changes we can express the change in terms of using the relation

where the derivative measures how responsive is to changes in . The derivative can be computed if we know the rate law for . For example, if we assume that the rate law is then the derivative is . We can also use a similar strategy to compute the change in as a result of the change . This time the change in is a result of two changes, the change in itself and the change in . We can express these changes by summing the two individual contributions:

We have two equations, one describing the change in and the other in . Because we allowed the system to settle to a new steady state we can also state that the change in reaction rates must be the same (otherwise it wouldn't be at steady state). That is we can assert that . With this in mind we equate the two equations and write

Solving for the ratio we obtain:

In the limit, as we make the change smaller and smaller, the left-hand side converges to the derivative :

We can go one step further and scale the derivatives to eliminate units. Multiplying both sides by and dividing both sides by yields the scaled derivatives:

The scaled derivatives on the right-hand side are the elasticities, and the scaled left-hand term is the scaled sensitivity coefficient or concentration control coefficient,

We can simplify this expression further. The reaction rate is usually a linear function of . For example, in the Briggs–Haldane equation, the reaction rate is given by . Differentiating this rate law with respect to and scaling yields .

Using this result gives:

A similar analysis can be done where is perturbed. In this case we obtain the sensitivity of with respect to :

The above expressions measure how much enzymes and control the steady state concentration of intermediate . We can also consider how the steady state reaction rates and are affected by perturbations in and . This is often of importance to metabolic engineers who are interested in increasing rates of production. At steady state the reaction rates are often called the fluxes and abbreviated to and . For a linear pathway such as this example, both fluxes are equal at steady-state so that the flux through the pathway is simply referred to as . Expressing the change in flux as a result of a perturbation in and taking the limit as before we obtain

The above expressions tell us how much enzymes and control the steady state flux. The key point here is that changes in enzyme concentration, or equivalently the enzyme activity, must be brought about by an external action.

Derivation using the systems equation

The control equations can also be derived in a more rigorous fashion using the systems equation:

where is the stoichiometry matrix, is a vector of chemical species, and is a vector of parameters (or inputs) that can influence the system. In metabolic control analysis the key parameters are the enzyme concentrations. This approach was popularized by Heinrich, Rapoport, and Rapoport [16] and Reder and Mazat. [17] A detailed discussion of this approach can be found in Heinrich & Schuster [18] and Hofmeyr. [19]

Properties of a linear pathway

A linear biochemical pathway is a chain of enzyme-catalyzed reaction steps. The figure below shows a three step pathway, with intermediates, and . In order to sustain a steady-state, the boundary species and are fixed.

Linear chain of four reactions catalyzed by enzymes e1 to e4 Linear3Step.png
Linear chain of four reactions catalyzed by enzymes e1 to e4

At steady-state the rate of reaction is the same at each step. This means there is an overall flux from X_o to X_1.

Linear pathways possess some well-known properties: [20] [21] [22]

  1. Flux control is biased towards the first few steps of the pathway. Flux control shifts more to the first step as the equilibrium constants become large.
  2. Flux control is small at reactions close to equilibrium.
  3. Assuming reversibly, flux control at a given step is proportional to the product of the equilibrium constants. For example, flux control at the second step in a three step pathway is proportional to the product of the second and third equilibrium constants.

In all cases, a rationale for these behaviors is given in terms of how elasticities transmit changes through a pathway.

Metabolic control analysis software

There are a number of software tools that can directly compute elasticities and control coefficients:

Relationship to Classical Control Theory

Classical Control theory is a field of mathematics that deals with the control of dynamical systems in engineered processes and machines. In 2004 Brian Ingalls published a paper [26] that showed that classical control theory and metabolic control analysis were identical. The only difference was that metabolic control analysis was confined to zero frequency responses when cast in the frequency domain whereas classical control theory imposes no such restriction. The other significant difference is that classical control theory [27] has no notion of stoichiometry and conservation of mass which makes it more cumbersome to use but also means it fails to recognize the structural properties inherent in stoichiometric networks which provide useful biological insights.

See also

Related Research Articles

<span class="mw-page-title-main">Dirac delta function</span> Generalized function whose value is zero everywhere except at zero

In mathematical analysis, the Dirac delta distribution, also known as the unit impulse, is a generalized function or distribution over the real numbers, whose value is zero everywhere except at zero, and whose integral over the entire real line is equal to one.

<span class="mw-page-title-main">Noether's theorem</span> Statement relating differentiable symmetries to conserved quantities

Noether's theorem or Noether's first theorem states that every differentiable symmetry of the action of a physical system with conservative forces has a corresponding conservation law. The theorem was proven by mathematician Emmy Noether in 1915 and published in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries over physical space.

In continuum mechanics, the infinitesimal strain theory is a mathematical approach to the description of the deformation of a solid body in which the displacements of the material particles are assumed to be much smaller than any relevant dimension of the body; so that its geometry and the constitutive properties of the material at each point of space can be assumed to be unchanged by the deformation.

<span class="mw-page-title-main">Hooke's law</span> Physical law: force needed to deform a spring scales linearly with distance

In physics, Hooke's law is an empirical law which states that the force needed to extend or compress a spring by some distance scales linearly with respect to that distance—that is, Fs = kx, where k is a constant factor characteristic of the spring, and x is small compared to the total possible deformation of the spring. The law is named after 17th-century British physicist Robert Hooke. He first stated the law in 1676 as a Latin anagram. He published the solution of his anagram in 1678 as: ut tensio, sic vis. Hooke states in the 1678 work that he was aware of the law since 1660.

In the calculus of variations and classical mechanics, the Euler–Lagrange equations are a system of second-order ordinary differential equations whose solutions are stationary points of the given action functional. The equations were discovered in the 1750s by Swiss mathematician Leonhard Euler and Italian mathematician Joseph-Louis Lagrange.

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

<span class="mw-page-title-main">Vapnik–Chervonenkis theory</span> Branch of statistical computational learning theory

Vapnik–Chervonenkis theory was developed during 1960–1990 by Vladimir Vapnik and Alexey Chervonenkis. The theory is a form of computational learning theory, which attempts to explain the learning process from a statistical point of view.

A directional derivative is a concept in multivariable calculus that measures the rate at which a function changes in a particular direction at a given point.

<span class="mw-page-title-main">Granular material</span> Conglomeration of discrete solid, macroscopic particles

A granular material is a conglomeration of discrete solid, macroscopic particles characterized by a loss of energy whenever the particles interact. The constituents that compose granular material are large enough such that they are not subject to thermal motion fluctuations. Thus, the lower size limit for grains in granular material is about 1 μm. On the upper size limit, the physics of granular materials may be applied to ice floes where the individual grains are icebergs and to asteroid belts of the Solar System with individual grains being asteroids.

In mathematics, in particular in algebraic geometry and differential geometry, Dolbeault cohomology (named after Pierre Dolbeault) is an analog of de Rham cohomology for complex manifolds. Let M be a complex manifold. Then the Dolbeault cohomology groups depend on a pair of integers p and q and are realized as a subquotient of the space of complex differential forms of degree (p,q).

Elastic energy is the mechanical potential energy stored in the configuration of a material or physical system as it is subjected to elastic deformation by work performed upon it. Elastic energy occurs when objects are impermanently compressed, stretched or generally deformed in any manner. Elasticity theory primarily develops formalisms for the mechanics of solid bodies and materials. The elastic potential energy equation is used in calculations of positions of mechanical equilibrium. The energy is potential as it will be converted into other forms of energy, such as kinetic energy and sound energy, when the object is allowed to return to its original shape (reformation) by its elasticity.

The rate of a chemical reaction is influenced by many different factors, such as temperature, pH, reactant, and product concentrations and other effectors. The degree to which these factors change the reaction rate is described by the elasticity coefficient. This coefficient is defined as follows:

The narrow escape problem is a ubiquitous problem in biology, biophysics and cellular biology.

Computational anatomy (CA) is the study of shape and form in medical imaging. The study of deformable shapes in computational anatomy rely on high-dimensional diffeomorphism groups which generate orbits of the form . In CA, this orbit is in general considered a smooth Riemannian manifold since at every point of the manifold there is an inner product inducing the norm on the tangent space that varies smoothly from point to point in the manifold of shapes . This is generated by viewing the group of diffeomorphisms as a Riemannian manifold with , associated to the tangent space at . This induces the norm and metric on the orbit under the action from the group of diffeomorphisms.

Control coefficients are used to describe how much influence a given reaction step has on the steady-state flux or species concentration level. In practice, this can be accomplished by changing the expression level of a given enzyme and measuring the resulting changes in flux and metabolite levels. Control coefficients form a central component of metabolic control analysis.

<span class="mw-page-title-main">Branched pathways</span> Common pattern in metabolism

Branched pathways, also known as branch points, are a common pattern found in metabolism. This is where an intermediate species is chemically made or transformed by multiple enzymatic processes. linear pathways only have one enzymatic reaction producing a species and one enzymatic reaction consuming the species.

Control coefficients measure the response of a biochemical pathway to changes in enzyme activity. The response coefficient, as originally defined by Kacser and Burns, is a measure of how external factors such as inhibitors, pharmaceutical drugs, or boundary species affect the steady-state fluxes and species concentrations. The flux response coefficient is defined by:

In metabolic control analysis, a variety of theorems have been discovered and discussed in the literature. The most well known of these are flux and concentration control coefficient summation relationships. These theorems are the result of the stoichiometric structure and mass conservation properties of biochemical networks. Equivalent theorems have not been found, for example, in electrical or economic systems.

The stoichiometric structure and mass-conservation properties of biochemical pathways gives rise to a series of theorems or relationships between the control coefficients and the control coefficients and elasticities. There are a large number of such relationships depending on the pathway configuration which have been documented and discovered by various authors. The term theorem has been used to describe these relationships because they can be proved in terms of more elementary concepts. The operational proofs in particular are of this nature.

Biological cells consume nutrients to sustain life. These nutrients are broken down to smaller molecules which are then reassembled into more complex structures required for life. The breakdown and reassembly of nutrients is called metabolism. An individual cell will contain thousands of different kinds of small molecules, such as sugars, lipids, amino acids, etc. The interconversion of these molecules is carried out by catalysts called enzymes. For example, E. coli contains 2,338 metabolic enzymes. These enzymes form a complex web of reactions forming pathways by which nutrients are converted. These pathways come in various forms. For example, enzymes can form simple linear sequences of reactions forming a linear biochemical pathway. Pathways can also show other patterns such as branches or cycles. A famous cyclic pathway is Kreb's cycle or the Calvin's cycle.

References

  1. Fell D., (1997) Understanding the Control of Metabolism, Portland Press.
  2. Heinrich R. and Schuster S. (1996) The Regulation of Cellular Systems, Chapman and Hall.
  3. Salter, M.; Knowles, R. G.; Pogson, C. I. (1994). "Metabolic control". Essays in Biochemistry. 28: 1–12. PMID   7925313.
  4. Ingalls, B. P. (2004) A Frequency Domain Approach to Sensitivity Analysis of Biochemical Systems, Journal of Physical Chemistry B, 108, 1143-1152.
  5. Savageau M.A (1976) Biochemical systems analysis: a study of function and design in molecular biology, Reading, MA, Addison–Wesley.
  6. Higgins, J. (1963). "Analysis of sequential reactions". Annals of the New York Academy of Sciences. 108 (1): 305–321. Bibcode:1963NYASA.108..305H. doi:10.1111/j.1749-6632.1963.tb13382.x. PMID   13954410. S2CID   30821044.
  7. Feinberg, Martin (1987). "Chemical reaction network structure and the stability of complex isothermal reactors—I. The deficiency zero and deficiency one theorems". Chemical Engineering Science. 42 (10): 2229–2268. doi:10.1016/0009-2509(87)80099-4.
  8. Clarke, Bruce L. (January 1980), Prigogine, I.; Rice, Stuart A. (eds.), "Stability of Complex Reaction Networks", Advances in Chemical Physics (1 ed.), Wiley, vol. 43, pp. 1–215, doi:10.1002/9780470142622.ch1, ISBN   978-0-471-05741-3 , retrieved 2023-12-06
  9. Shinar, Guy; Alon, Uri; Feinberg, Martin (2009). "Sensitivity and Robustness in Chemical Reaction Networks". SIAM Journal on Applied Mathematics. 69 (4): 977–998. ISSN   0036-1399.
  10. 1 2 3 4 Kacser, H.; Burns, J. A. (1973). "The control of flux". Symposia of the Society for Experimental Biology. 27: 65–104. PMID   4148886.
  11. 1 2 3 Heinrich, R.; Rapoport, T. A. (1974). "A linear steady-state treatment of enzymatic chains. General properties, control and effector strength". European Journal of Biochemistry. 42 (1): 89–95. doi: 10.1111/j.1432-1033.1974.tb03318.x . PMID   4830198.
  12. Burns, J.A.; Cornish-Bowden, A.; Groen, A.K.; Heinrich, R.; Kacser, H.; Porteous, J.W.; Rapoport, S.M.; Rapoport, T.A.; Stucki, J.W.; Tager, J.M.; Wanders, R.J.A.; Westerhoff, H.V. (1985). "Control analysis of metabolic systems". Trends Biochem. Sci. 10: 16. doi:10.1016/0968-0004(85)90008-8.
  13. Liebermeister, Wolfram (11 May 2022). "Structural Thermokinetic Modelling". Metabolites. 12 (5): 434. doi: 10.3390/metabo12050434 . PMC   9144996 . PMID   35629936.
  14. Liebermeister, Wolfram (2016). "Optimal enzyme rhythms in cells". arXiv: 1602.05167 .{{cite journal}}: Cite journal requires |journal= (help)
  15. Moreno-Sánchez, Rafael; Saavedra, Emma; Rodríguez-Enríquez, Sara; Olín-Sandoval, Viridiana (2008). "Metabolic Control Analysis: A Tool for Designing Strategies to Manipulate Metabolic Pathways". Journal of Biomedicine and Biotechnology. 2008: 1–30. doi: 10.1155/2008/597913 . PMC   2447884 . PMID   18629230.
  16. Heinrich, R.; Rapoport, S. M.; Rapoport, T. A. (1 January 1978). "Metabolic regulation and mathematical models". Progress in Biophysics and Molecular Biology. 32 (1): 1–82. doi: 10.1016/0079-6107(78)90017-2 . PMID   343173.
  17. Reder, Christine (November 1988). "Metabolic control theory: A structural approach". Journal of Theoretical Biology. 135 (2): 175–201. Bibcode:1988JThBi.135..175R. doi:10.1016/S0022-5193(88)80073-0. PMID   3267767.
  18. Heinrich, Reinhart; Schuster, Stefan (1996). "The Regulation of Cellular Systems". SpringerLink. doi:10.1007/978-1-4613-1161-4. ISBN   978-1-4612-8492-5. S2CID   10252429.
  19. Hofmeyr, Jan-Hendrik. "Metabolic control analysis in a nutshell". Proceedings of the 2nd International Conference on Systems Biology.
  20. Heinrich, Reinhart; Rapoport, Tom A. (February 1974). "A Linear Steady-State Treatment of Enzymatic Chains. General Properties, Control and Effector Strength". European Journal of Biochemistry. 42 (1): 89–95. doi: 10.1111/j.1432-1033.1974.tb03318.x .
  21. Savageau, Michael (1976). Biochemical systems analysis. A study of function and design in molecular biology. Addison-Wesley.
  22. Sauro, Herbert (28 August 2020). "A brief note on the properties of linear pathways". Biochemical Society Transactions. 48 (4): 1379–1395. doi:10.1042/BST20190842.
  23. Olivier, B. G.; Rohwer, J. M.; Hofmeyr, J.-H. S. (15 February 2005). "Modelling cellular systems with PySCeS". Bioinformatics. 21 (4): 560–561. doi: 10.1093/bioinformatics/bti046 . PMID   15454409.
  24. Bergmann, Frank T.; Sauro, Herbert M. (December 2006). "SBW - A Modular Framework for Systems Biology". Proceedings of the 2006 Winter Simulation Conference. pp. 1637–1645. doi:10.1109/WSC.2006.322938. ISBN   1-4244-0501-7.
  25. Choi, Kiri; Medley, J. Kyle; König, Matthias; Stocking, Kaylene; Smith, Lucian; Gu, Stanley; Sauro, Herbert M. (September 2018). "Tellurium: An extensible python-based modeling environment for systems and synthetic biology". Biosystems. 171: 74–79. doi:10.1016/j.biosystems.2018.07.006. PMC   6108935 . PMID   30053414.
  26. Drengstig, Tormod; Kjosmoen, Thomas; Ruoff, Peter (19 May 2011). "On the Relationship between Sensitivity Coefficients and Transfer Functions of Reaction Kinetic Networks". The Journal of Physical Chemistry B. 115 (19): 6272–6278. doi:10.1021/jp200578e. PMID   21520979.
  27. Nise, Norman S. (2019). Control systems engineering (Eighth, Wiley abridged print companion ed.). Hoboken, NJ. ISBN   978-1119592921.{{cite book}}: CS1 maint: location missing publisher (link)