Metabolic control analysis

{{Short description|Mathematical model of biochemical pathways}}

File:MCAlogo.png

In biochemistry, 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.Fell D., (1997) Understanding the Control of Metabolism, Portland Press.Heinrich R. and Schuster S. (1996) The Regulation of Cellular Systems, Chapman and Hall.

{{Cite journal

| last1 = Salter | first1 = M.

| last2 = Knowles | first2 = R. G.

| last3 = Pogson | first3 = C. I.

| title = Metabolic control

| journal = Essays in Biochemistry

| volume = 28

| pages = 1–12

| year = 1994

| pmid = 7925313

}}

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|date=November 2010}}.

More recent workIngalls, B. P. (2004) A Frequency Domain Approach to Sensitivity Analysis of Biochemical Systems, Journal of Physical Chemistry B, 108, 1143-1152. has shown that MCA can be mapped directly on to classical control theory and are as such equivalent.

Biochemical systems theorySavageau M.A (1976) Biochemical systems analysis: a study of function and design in molecular biology, Reading, MA, Addison–Wesley. (BST) is a similar formalism, though with rather different objectives. Both are evolutions of an earlier theoretical analysis by Joseph Higgins.{{Cite journal

| doi = 10.1111/j.1749-6632.1963.tb13382.x

| last1 = Higgins | first1 = J.

| title = Analysis of sequential reactions

| journal = Annals of the New York Academy of Sciences

| volume = 108

| pages = 305–321

| year = 1963

| issue = 1 | pmid = 13954410

| bibcode = 1963NYASA.108..305H | s2cid = 30821044 }}

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.{{Cite journal |last=Feinberg |first=Martin |date=1987 |title=Chemical reaction network structure and the stability of complex isothermal reactors—I. The deficiency zero and deficiency one theorems |url=https://linkinghub.elsevier.com/retrieve/pii/0009250987800994 |journal=Chemical Engineering Science |language=en |volume=42 |issue=10 |pages=2229–2268 |doi=10.1016/0009-2509(87)80099-4}} Its emphasis is primarily on dynamic stability criteria{{Citation |last=Clarke |first=Bruce L. |title=Stability of Complex Reaction Networks |date=January 1980 |url=https://onlinelibrary.wiley.com/doi/10.1002/9780470142622.ch1 |work=Advances in Chemical Physics |volume=43 |pages=1–215 |editor-last=Prigogine |editor-first=I. |access-date=2023-12-06 |edition=1 |publisher=Wiley |language=en |doi=10.1002/9780470142622.ch1 |isbn=978-0-471-05741-3 |editor2-last=Rice |editor2-first=Stuart A.}} and related theorems associated with mass-action networks. In more recent years the field has also developed {{Cite journal |last=Shinar |first=Guy |last2=Alon |first2=Uri |last3=Feinberg |first3=Martin |date=2009 |title=Sensitivity and Robustness in Chemical Reaction Networks |url=https://www.jstor.org/stable/40233656 |journal=SIAM Journal on Applied Mathematics |volume=69 |issue=4 |pages=977–998 |issn=0036-1399}} a sensitivity analysis which is similar if not identical to MCA and BST.

Control coefficients

{{main article|Control coefficient (biochemistry)}}

A control coefficient{{Cite journal

| last1 = Kacser | first1 = H.

| last2 = Burns | first2 = J. A.

| title = The control of flux

| journal = Symposia of the Society for Experimental Biology

| volume = 27

| pages = 65–104

| year = 1973

| pmid = 4148886

}}

{{Cite journal

| doi = 10.1111/j.1432-1033.1974.tb03318.x

| last1 = Heinrich | first1 = R.

| last2 = Rapoport | first2 = T. A.

| title = A linear steady-state treatment of enzymatic chains. General properties, control and effector strength

| journal = European Journal of Biochemistry

| volume = 42

| issue = 1

| pages = 89–95

| year = 1974

| pmid = 4830198

| doi-access = free

}}{{cite journal | last1 = Burns | first1 = J.A. | last2 = Cornish-Bowden | first2 = A. | last3 = Groen | first3 = A.K. | last4 = Heinrich | first4 = R. | last5 = Kacser | first5 = H. | last6 = Porteous | first6 = J.W. | last7 = Rapoport | first7 = S.M. | last8 = Rapoport | first8 = T.A. | last9 = Stucki | first9 = J.W. | last10 = Tager | first10 = J.M. | last11 = Wanders | first11 = R.J.A. | last12 = Westerhoff | first12 = H.V. | year = 1985 | title = Control analysis of metabolic systems | journal = Trends Biochem. Sci. | volume = 10 | page = 16 | doi=10.1016/0968-0004(85)90008-8}} 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 ( v_i ) of step i. The two main control coefficients are the flux and concentration control coefficients. Flux control coefficients are defined by

: C^J_{v_i} = \left( \frac{dJ}{dp} \frac{p}{J} \right) \bigg/ \left( \frac{\partial v_i}{\partial p}\frac{p}{v_i} \right) = \frac{d\ln J}{d\ln v_i}

and concentration control coefficients by

: C^S_{v_i} = \left( \frac{dS}{dp} \frac{p}{S} \right) \bigg/ \left( \frac{\partial v_i}{\partial p} \frac{p}{v_i} \right) = \frac{d\ln S}{d\ln v_i}

File:SystemControl.png.

= Summation theorems =

{{main article|Summation theorems (biochemistry)}}

The flux control summation theorem was discovered independently by the Kacser/Burns group and the Heinrich/Rapoport group 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.

: \sum_i C^J_{v_i} = 1

: \sum_i C^s_{v_i} = 0

= Elasticity coefficients =

{{main article|Elasticity coefficient}}

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.

File:LocalControl.png.

= Connectivity theorems =

{{main article|Connectivity theorems (biochemistry)}}

The connectivity theorems 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 S_n is different from the local species S_m .

: \sum_i C^J_i \varepsilon^i_s = 0

: \sum_i C^{s_n}_i \varepsilon^i_{s_m} = 0 \quad n \neq m

: \sum_i C^{s_n}_i \varepsilon^i_{s_m} = -1 \quad n = m

Response Coefficient

{{main article|Response coefficient (biochemistry)}}

Kacser and Burns 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.{{cite journal |last1=Moreno-Sánchez |first1=Rafael |last2=Saavedra |first2=Emma |last3=Rodríguez-Enríquez |first3=Sara |last4=Olín-Sandoval |first4=Viridiana |title=Metabolic Control Analysis: A Tool for Designing Strategies to Manipulate Metabolic Pathways |journal=Journal of Biomedicine and Biotechnology |date=2008 |volume=2008 |pages=1–30 |doi=10.1155/2008/597913 |pmid=18629230 |pmc=2447884 |doi-access=free }}

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

R^X_m = C^X_{i} \varepsilon^{i}_m

where X is a chosen observable such as a flux or metabolite concentration, i is the step that the external factor targets, C^X_{i} is the control coefficient of the target steps, and \varepsilon^{i}_m is the elasticity of the target step with respect to the external factor m .

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 \varepsilon^{i}_m 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 C^X_{i} .

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, m , targets multiple sites of action, for example n sites, then the overall response in a phenotypic factor X , is the sum of the individual responses:

R_m^X=\sum_{i=1}^n C_{i}^X \varepsilon_m^{i}

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:

: X_o \rightarrow S \rightarrow X_1

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

: C^J_{v_1} + C^J_{v_2} = 1

: C^J_{v_1} \varepsilon^{v_1}_s + C^J_{v_2} \varepsilon^{v_2}_s = 0

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

: C^J_{v_1} = \frac{\varepsilon^{2}_s}{\varepsilon^{2}_s - \varepsilon^{1}_s}

: C^J_{v_2} = \frac{-\varepsilon^{1}_s}{\varepsilon^{2}_s - \varepsilon^{1}_s}

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 \varepsilon^{v_1}_s = 0 . In this case, the control coefficients reduce to

: C^J_{v_1} = 1

: C^J_{v_2} = 0

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:

: C^s_{v_1} = \frac{1}{\varepsilon^{2}_s - \varepsilon^1_s}

: C^s_{v_2} = \frac{-1}{\varepsilon^{2}_s - \varepsilon^1_s}

Three step pathway

Consider the simple three step pathway:

: X_o \rightarrow S_1 \rightarrow S_2 \rightarrow X_1

where X_o and X_1 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.

:

C^J_{e_1} = \varepsilon^{2}_1 \varepsilon^{3}_2 / D

:

C^J_{e_2} = -\varepsilon^{1}_1 \varepsilon^{3}_2 / D

:

C^J_{e_3} = \varepsilon^{1}_1 \varepsilon^{2}_2 / D

where D the denominator is given by

:

D = \varepsilon^{2}_1 \varepsilon^{3}_2 -\varepsilon^{1}_1 \varepsilon^{3}_2 + \varepsilon^{1}_1 \varepsilon^{2}_2

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 S_1

:

C^{S_1}_{e_1} = (\varepsilon^{3}_2 - \varepsilon^{2}_2) / D

:

C^{S_1}_{e_2} = - \varepsilon^{3}_2 / D

:

C^{S_1}_{e_3} = \varepsilon^{2}_2 / D

And for S_2

:

C^{S_2}_{e_1} = \varepsilon^{2}_1 / D

:

C^{S_2}_{e_2} = -\varepsilon^{1}_1 / D

:

C^{S_2}_{e_3} = (\varepsilon^{1}_1 - \varepsilon^{2}_1) / D

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 v_1 and v_2 are determined by two enzymes e_1 and e_2 respectively. Changing either enzyme will result in a change to the steady state level of x and the steady state reaction rates v. Consider a small change in e_1 of magnitude \delta e_1. This will have a number of effects, it will increase v_1 which in turn will increase x which in turn will increase v_2. Eventually the system will settle to a new steady state. We can describe these changes by focusing on the change in v_1 and v_2. The change in v_2, which we designate \delta v_2, came about as a result of the change \delta x. Because we are only considering small changes we can express the change \delta v_2 in terms of \delta x using the relation

:

\delta v_2 = \frac{\partial v_2}{\partial x} \delta x

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

:

\delta v_1 = \frac{\partial v_1}{\partial e_1} \delta e_1 + \frac{\partial v_1}{\partial x} \delta x

We have two equations, one describing the change in v_1 and the other in v_2. 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 \delta v_1 = \delta v_2. With this in mind we equate the two equations and write

: \frac{\partial v_2}{\partial x} \delta x = \frac{\partial v_1}{\partial e_1} \delta e_1 + \frac{\partial v_1}{\partial x} \delta x

Solving for the ratio \delta x/\delta e_1 we obtain:

: \frac{\delta x}{\delta e_1} = \dfrac{-\dfrac{\partial v_1}{\partial e_1}}{\dfrac{\partial v_2}{\partial x} - \dfrac{\partial v_1}{\partial x}}

In the limit, as we make the change \delta e_1 smaller and smaller, the left-hand side converges to the derivative dx/de_1:

: \lim_{\delta e_1 \rightarrow 0} \frac{\delta x}{\delta e_1} = \frac{dx}{de_1} = \dfrac{-\dfrac{\partial v_1}{\partial e_1}}{\dfrac{\partial v_2}{\partial x} - \dfrac{\partial v_1}{\partial x}}

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

: \frac{dx}{de_1} \frac{e_1}{x}= \frac{-\dfrac{\partial v_1}{\partial e_1}\dfrac{e_1}{v_1}}

{\dfrac{\partial v_2}{\partial x} \dfrac{x}{v_2} - \dfrac{\partial v_1}{\partial x} \dfrac{x}{v_1}}

The scaled derivatives on the right-hand side are the elasticities, \varepsilon^v_x and the scaled left-hand term is the scaled sensitivity coefficient or concentration control coefficient, C^x_{e}

: C^x_{e_1} = \frac{\varepsilon^1_{e_1}}{\varepsilon^2_x - \varepsilon^1_x}

We can simplify this expression further. The reaction rate v_1 is usually a linear function of e_1. For example, in the Briggs–Haldane equation, the reaction rate is given by v= e_1 k_{cat} x/(K_m + x). Differentiating this rate law with respect to e_1 and scaling yields \varepsilon^{v_{1}}_{e_1} = 1.

Using this result gives:

:

C^x_{e_1} = \frac{1}{\varepsilon^2_x - \varepsilon^1_x}

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

:

C^x_{e_2} = -\frac{1}{\varepsilon^2_x - \varepsilon^1_x}

The above expressions measure how much enzymes e_1 and e_2 control the steady state concentration of intermediate x. We can also consider how the steady state reaction rates v_1 and v_2 are affected by perturbations in e_1 and e_2. 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 J_1 and J_2. 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 J. Expressing the change in flux as a result of a perturbation in e_1 and taking the limit as before we obtain

:

C^J_{e_1} = \frac{\varepsilon^1_x}{\varepsilon^2_x - \varepsilon^1_x}, \quad C^J_{e_2} = \frac{-\varepsilon^1_x}{\varepsilon^2_x - \varepsilon^1_x}

The above expressions tell us how much enzymes e_1 and e_2 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:

\dfrac{{\bf dx}}{dt} = {\bf N} {\bf v} ({\bf x} (p), p)

where {\bf N} is the stoichiometry matrix, {\bf x} is a vector of chemical species, and {\bf p} 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{{cite journal |last1=Heinrich |first1=R. |last2=Rapoport |first2=S. M. |last3=Rapoport |first3=T. A. |title=Metabolic regulation and mathematical models |journal=Progress in Biophysics and Molecular Biology |date=1 January 1978 |volume=32 |issue=1 |pages=1–82 |doi=10.1016/0079-6107(78)90017-2 |pmid=343173 |doi-access=free }} and Reder and Mazat.{{cite journal |last1=Reder |first1=Christine |title=Metabolic control theory: A structural approach |journal=Journal of Theoretical Biology |date=November 1988 |volume=135 |issue=2 |pages=175–201 |doi=10.1016/S0022-5193(88)80073-0|pmid=3267767 |bibcode=1988JThBi.135..175R }} A detailed discussion of this approach can be found in Heinrich & Schuster{{cite journal |last1=Heinrich |first1=Reinhart |last2=Schuster |first2=Stefan |title=The Regulation of Cellular Systems |journal=SpringerLink |date=1996 |doi=10.1007/978-1-4613-1161-4 |isbn=978-1-4612-8492-5 |s2cid=10252429 |url=https://link.springer.com/book/10.1007/978-1-4613-1161-4 |language=en}} and Hofmeyr.{{cite journal |last1=Hofmeyr |first1=Jan-Hendrik |title=Metabolic control analysis in a nutshell |journal=Proceedings of the 2nd International Conference on Systems Biology|url=https://www.academia.edu/4610732}}

Properties of a linear pathway

{{main article|Linear biochemical pathway}}

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

File:Linear3Step.png

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:{{cite journal |last1=Heinrich |first1=Reinhart |last2=Rapoport |first2=Tom A. |title=A Linear Steady-State Treatment of Enzymatic Chains. General Properties, Control and Effector Strength |journal=European Journal of Biochemistry |date=February 1974 |volume=42 |issue=1 |pages=89–95 |doi=10.1111/j.1432-1033.1974.tb03318.x|doi-access=free }}{{cite book |last1=Savageau |first1=Michael |title=Biochemical systems analysis. A study of function and design in molecular biology. |date=1976 |publisher=Addison-Wesley}}{{cite journal |last1=Sauro |first1=Herbert |title=A brief note on the properties of linear pathways |journal=Biochemical Society Transactions |date=28 August 2020 |volume=48 |issue=4 |pages=1379–1395 |doi=10.1042/BST20190842}}

  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:

  • COPASI (GUI)
  • PySCeS{{cite journal |last1=Olivier |first1=B. G. |last2=Rohwer |first2=J. M. |last3=Hofmeyr |first3=J.-H. S. |title=Modelling cellular systems with PySCeS |journal=Bioinformatics |date=15 February 2005 |volume=21 |issue=4 |pages=560–561 |doi=10.1093/bioinformatics/bti046|pmid=15454409 |doi-access=free }} (Python)
  • SBW{{cite book |last1=Bergmann |first1=Frank T. |last2=Sauro |first2=Herbert M. |title=Proceedings of the 2006 Winter Simulation Conference |chapter=SBW - A Modular Framework for Systems Biology |date=December 2006 |pages=1637–1645 |doi=10.1109/WSC.2006.322938|isbn=1-4244-0501-7 }} (GUI)
  • libroadrunner{{cite journal |last1=Choi |first1=Kiri |last2=Medley |first2=J. Kyle |last3=König |first3=Matthias |last4=Stocking |first4=Kaylene |last5=Smith |first5=Lucian |last6=Gu |first6=Stanley |last7=Sauro |first7=Herbert M. |title=Tellurium: An extensible python-based modeling environment for systems and synthetic biology |journal=Biosystems |date=September 2018 |volume=171 |pages=74–79 |doi=10.1016/j.biosystems.2018.07.006|pmid=30053414 |pmc=6108935 }} (Python)
  • VCell

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{{cite journal |last1=Drengstig |first1=Tormod |last2=Kjosmoen |first2=Thomas |last3=Ruoff |first3=Peter |title=On the Relationship between Sensitivity Coefficients and Transfer Functions of Reaction Kinetic Networks |journal=The Journal of Physical Chemistry B |date=19 May 2011 |volume=115 |issue=19 |pages=6272–6278 |doi=10.1021/jp200578e|pmid=21520979 }} 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{{cite book |last1=Nise |first1=Norman S. |title=Control systems engineering |date=2019 |location=Hoboken, NJ |isbn=978-1119592921 |edition=Eighth, Wiley abridged print companion}} 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

References

{{reflist}}