Misplaced Pages

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:

${\displaystyle R_{x}^{J}={\frac {dJ}{dx}}{\frac {x}{J}}}$

where ${\displaystyle J}$ is the steady-state pathway flux. Similarly, the concentration response coefficient is defined by the expression:

${\displaystyle R_{x}^{s}={\frac {ds}{dx}}{\frac {x}{s}}}$

where in both cases ${\displaystyle x}$ is the concentration of the external factor. The response coefficient measures how sensitive a pathway is to changes in external factors other than enzyme activities.

The flux response coefficient is related to control coefficients and elasticities through the following relationship:

${\displaystyle R_{x}^{J}=\sum _{i=1}^{n}C_{e_{i}}^{J}\varepsilon _{x}^{v_{i}}}$

Likewise, the concentration response coefficient is related by the following expression:

${\displaystyle R_{x}^{s}=\sum _{i=1}^{n}C_{e_{i}}^{s}\varepsilon _{x}^{v_{i}}}$

The summation in both cases accounts for cases where a given external factor, ${\displaystyle x}$, can act at multiple sites. For example, a given drug might act on multiple protein sites. The overall response is the sum of the individual responses.

These results show that the action of an external factor, such as a drug, has two components:

1. The elasticity indicates how potent the drug is at affecting the activity of the target site itself.
2. The control coefficient indicates how any perturbation at the target site will propagate to the rest of the system and thereby affect the phenotype.

When designing drugs for therapeutic action, both aspects must therefore be considered.

Proof of Response Theorem

There are various ways to prove the response theorems:

Proof by perturbation

The perturbation proof by Kacser and Burns is given as follows.

Given the simple linear pathway catalyzed by two enzymes ${\displaystyle e_{1}}$ and ${\displaystyle e_{2}}$:

${\displaystyle X{\stackrel {e_{1}}{\longrightarrow }}S{\stackrel {e_{2}}{\longrightarrow }}}$

where ${\displaystyle X}$ is the fixed boundary species. Let us increase the concentration of enzyme ${\displaystyle e_{1}}$ by an amount ${\displaystyle \delta e_{1}}$. This will cause the steady state flux and concentration of ${\displaystyle S}$, and all downstream species beyond ${\displaystyle e_{2}}$ to increase. The concentration of ${\displaystyle X}$ is now decreased such that the flux and steady-state concentration of ${\displaystyle S}$ is restored back to their original values. These changes allow one to write down the following local and systems equations for the changes that occurred:

${\displaystyle {\begin{array}{r}\left.{\dfrac {\delta v_{1}}{v_{1}}}=\varepsilon _{x}^{1}{\dfrac {\delta x}{x}}+\varepsilon _{e_{1}}^{1}{\dfrac {\delta e_{1}}{e_{1}}}=0\right\}{\text{ Local equation }}\\\left.{\dfrac {\delta J}{J}}=R_{x}^{J}{\dfrac {\delta x}{x}}+C_{e_{1}}^{J}{\dfrac {\delta e_{1}}{e_{1}}}=0\right\}{\text{ System equation }}\end{array}}}$

There is no ${\displaystyle s}$ term in either equation because the concentration of ${\displaystyle s}$ is unchanged. Both right-hand sides of the equations are guaranteed to be zero by construction. The term ${\displaystyle \delta e_{1}/e_{1}}$ can be eliminated by combining both equations. If we also assume that the reaction rate for an enzyme-catalyzed reaction is proportional to the enzyme concentration, then ${\displaystyle \varepsilon _{e_{1}}^{1}=1}$, therefore:

${\displaystyle 0=R_{x}^{J}{\frac {\delta x}{x}}-C_{e_{1}}^{J}\varepsilon _{x}^{1}{\frac {\delta x}{x}}}$

Since ${\displaystyle \delta e_{1}/e_{1}\neq 0}$

this yields:

${\displaystyle R_{x}^{J}=C_{e_{1}}^{J}\varepsilon _{x}^{1}}$.

This proof can be generalized to the case where ${\displaystyle X}$ may act at multiple sites.

Pure algebraic proof

The pure algebraic proof is more complex and requires consideration of the system equation:

${\displaystyle {\bf {N}}{\bf {v}}(s(p),p)=0}$

where ${\displaystyle {\bf {N}}}$ is the stoichiometry matrix and ${\displaystyle {\bf {v}}}$ the rate vector. In this derivation, we assume there are no conserved moieties in the network, but this doesn't invalidate the proof. Using the chain rule and differentiating with respect to ${\displaystyle p}$ yields, after rearrangement:

${\displaystyle {\dfrac {ds}{dp}}=\left^{-1}{\dfrac {\partial v}{\partial p}}}$

The inverted term is the unscaled control coefficient so that after scaling, it is possible to write:

${\displaystyle R_{p}^{s}=C_{v}^{s}\varepsilon _{p}^{v}}$

To derive the flux response coefficient theorem, we must use the additional equation:

${\displaystyle {\bf {v}}={\bf {v}}({\bf {s}}(p),p)}$

See also

• Control coefficient (biochemistry)
• Elasticity coefficient
• Metabolic control analysis

References

1. ^ Kacser, H; Burns, JA (1973). "The control of flux". Symposia of the Society for Experimental Biology. 27: 65–104. PMID 4148886.
2. Cascante, Marta; Boros, Laszlo G.; Comin-Anduix, Begoña; de Atauri, Pedro; Centelles, Josep J.; Lee, Paul W.-N. (March 2002). "Metabolic control analysis in drug discovery and disease". Nature Biotechnology. 20 (3): 243–249. doi:10.1038/nbt0302-243. PMID 11875424. S2CID 3937563.
3. 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.
4. Hofmeyr, Jan-hendrik S. (2001). "Metabolic control analysis in a nutshell". In Proceedings of the 2 Nd International Conference on Systems Biology: 291–300. CiteSeerX 10.1.1.324.922.

• Biochemistry methods
• Metabolism
• Mathematical and theoretical biology
• Systems biology