Catalog of MCA Control Patterns

Jannie Hofmeyr published the first catalog of control patterns in metabolic control analysis (MCA). His doctoral research.^[1] concerned the use of graphical patterns to elucidate chains of interaction in metabolic regulation, later published in the European Journal of Biochemistry.^[2] In his thesis, he cataloged 25 patterns for various biochemical networks. In later work, his research group, together with Carl D Christensen and Johann Rohwer, developed a Python based tool called SymCA that was part of the PySCeSToolbox toolkit ^[3][4] that could generate patterns automatically and symbolically from a description of the network. This software was used to generate the patterns shown below.

The control equations, especially the numerators of the equations, can give information on the relative importance and routes by which perturbations travel through a biochemical network^[5]

Control patterns describe how a perturbation to a given parameter affects the steady-state level of a given variable. For example, a concentration control coefficient can describe how the overexpression of a specific enzyme can influence steady-state metabolite concentrations. Flux control coefficients are similar in that they describe how a perturbation in a given enzyme affects steady-state flux through a pathway. Such coefficients can be written in terms of elasticity coefficients.

Elasticity coefficients are local properties that describe how a single reaction is influenced by changes in the substrates and products that might influence the rate. For example, given a reaction such as:

${\displaystyle S{\stackrel {v}{\longrightarrow }}P}$

we will assume it has a rate of reaction of ${\displaystyle v}$. This reaction rate can be influenced by changes in the concentrations of substrate ${\displaystyle S}$ or product ${\displaystyle P}$. This influence is measured by an elasticity which is defined as:

${\displaystyle \varepsilon _{s}^{v}={\frac {\partial v}{\partial s}}{\frac {s}{v}}}$

To make the notation manageable, a specific numbering scheme is used in the following patterns. If a substrate has an index of ${\displaystyle i}$, then the reaction index will be ${\displaystyle v_{i+1}}$. The product elasticity will also have an index of ${\displaystyle i+1}$. This means that a product elasticity will have identical subscripts and superscripts making them easy to identify. The source boundary species is always labeled zero as well as the label for the first reaction.

For example, the following fragment of a network illustrates this labeling:

${\displaystyle X_{o}{\stackrel {v_{1}}{\longrightarrow }}S_{1}{\stackrel {v_{2}}{\longrightarrow }}S_{2}{\stackrel {v_{3}}{\longrightarrow }}}$

then

${\displaystyle \varepsilon _{1}^{2}={\frac {\partial v_{2}}{\partial s_{1}}}{\frac {s_{1}}{v_{2}}},\quad \varepsilon _{2}^{2}={\frac {\partial v_{2}}{\partial s_{2}}}{\frac {s_{2}}{v_{2}}},\quad \varepsilon _{2}^{3}={\frac {\partial v_{3}}{\partial s_{2}}}{\frac {s_{2}}{v_{3}}}}$

${\displaystyle X_{o}{\stackrel {v_{1}}{\longrightarrow }}S_{1}{\stackrel {v_{2}}{\longrightarrow }}X_{1}}$

${\displaystyle C_{e_{1}}^{J}=1\qquad C_{e_{2}}^{J}=0}$

${\displaystyle C_{e_{1}}^{s_{1}}={\frac {1}{\varepsilon _{1}^{2}}}\qquad C_{e_{2}}^{s_{1}}={\frac {-1}{\varepsilon _{1}^{2}}}}$

${\displaystyle C_{v_{1}}^{J}={\frac {\varepsilon _{1}^{2}}{\varepsilon _{1}^{2}-\varepsilon _{1}^{1}}}\qquad C_{v_{2}}^{J}={\frac {-\varepsilon _{1}^{1}}{\varepsilon _{1}^{2}-\varepsilon _{1}^{1}}}}$

${\displaystyle C_{v_{1}}^{s_{1}}={\frac {1}{\varepsilon _{1}^{2}-\varepsilon _{1}^{1}}}\qquad C_{v_{2}}^{s_{1}}={\frac {-1}{\varepsilon _{1}^{2}-\varepsilon _{1}^{1}}}}$

${\displaystyle X_{o}{\stackrel {v_{1}}{\longrightarrow }}S_{1}{\stackrel {v_{2}}{\longrightarrow }}S_{2}{\stackrel {v_{3}}{\longrightarrow }}X_{1}}$

Denominator:

${\displaystyle d=\varepsilon _{1}^{2}\varepsilon _{2}^{3}}$

Assume that each of the following expressions is divided by d

${\displaystyle {\begin{array}{lll}C_{e_{1}}^{J}=1&C_{e_{2}}^{J}=0&C_{e_{3}}^{J}=0\end{array}}}$

${\displaystyle {\begin{array}{ll}C_{e_{1}}^{s_{1}}=\varepsilon _{2}^{3}&C_{e_{1}}^{s_{2}}=\varepsilon _{1}^{2}\\[6pt]C_{e_{2}}^{s_{1}}=-\varepsilon _{2}^{3}&C_{e_{2}}^{s_{2}}=0\\[6pt]C_{e_{2}}^{s_{2}}=0&C_{e_{3}}^{s_{2}}=-\varepsilon _{1}^{2}\end{array}}}$

Denominator:

${\displaystyle d=\varepsilon _{1}^{2}\varepsilon _{2}^{3}-\varepsilon _{1}^{1}\varepsilon _{2}^{3}+\varepsilon _{1}^{1}\varepsilon _{2}^{2}}$

Assume that each of the following expressions is divided by ${\displaystyle d}$

${\displaystyle {\begin{array}{lll}C_{e_{1}}^{J}=\varepsilon _{1}^{2}\varepsilon _{2}^{3}&C_{e_{2}}^{J}=-\varepsilon _{1}^{1}\varepsilon _{2}^{3}&C_{e_{3}}^{J}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}\\[6pt]\end{array}}}$

${\displaystyle {\begin{array}{ll}C_{e_{1}}^{s_{1}}=\varepsilon _{2}^{3}-\varepsilon _{2}^{2}&C_{e_{1}}^{s_{2}}=\varepsilon _{1}^{2}\\[6pt]C_{e_{2}}^{s_{1}}=-\varepsilon _{2}^{3}&C_{e_{2}}^{s_{2}}=-\varepsilon _{1}^{1}\\[6pt]C_{e_{3}}^{s_{1}}=\varepsilon _{2}^{2}&C_{e_{3}}^{s_{2}}=\varepsilon _{1}^{1}-\varepsilon _{1}^{2}\end{array}}}$

${\displaystyle X_{o}{\stackrel {v_{1}}{\longrightarrow }}S_{1}{\stackrel {v_{2}}{\longrightarrow }}S_{2}{\stackrel {v_{3}}{\longrightarrow }}S_{3}{\stackrel {v_{4}}{\longrightarrow }}X_{1}}$

Denominator:

${\displaystyle d=\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{3}-\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{4}+\varepsilon _{1}^{1}\varepsilon _{2}^{3}\varepsilon _{3}^{4}-\varepsilon _{1}^{2}\varepsilon _{2}^{3}\varepsilon _{3}^{4}}$

Assume that each of the following expressions is divided by ${\displaystyle d}$.

${\displaystyle {\begin{array}{lll}C_{e_{1}}^{J}=-\varepsilon _{1}^{2}\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{e_{2}}^{J}=\varepsilon _{1}^{1}\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{e_{3}}^{J}=-\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{4}&C_{e_{4}}^{J}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{3}\\[4pt]C_{e_{1}}^{s_{1}}=-\varepsilon _{2}^{2}\varepsilon _{3}^{3}+\varepsilon _{2}^{2}\varepsilon _{3}^{4}-\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{e_{2}}^{s_{1}}=-\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{e_{3}}^{s_{1}}=-\varepsilon _{2}^{2}\varepsilon _{3}^{4}&C_{e_{4}}^{s_{1}}=\varepsilon _{2}^{2}\varepsilon _{3}^{3}\\[4pt]C_{e_{1}}^{s_{2}}=\varepsilon _{1}^{2}\varepsilon _{3}^{3}-\varepsilon _{1}^{2}\varepsilon _{3}^{4}&C_{e_{2}}^{s_{2}}=-\varepsilon _{1}^{1}\varepsilon _{3}^{3}+\varepsilon _{1}^{1}\varepsilon _{3}^{4}&C_{e_{3}}^{s_{2}}=-\varepsilon _{1}^{1}\varepsilon _{3}^{4}+\varepsilon _{1}^{2}\varepsilon _{3}^{4}&C_{e_{4}}^{s_{2}}=\varepsilon _{1}^{1}\varepsilon _{3}^{3}-\varepsilon _{1}^{2}\varepsilon _{3}^{3}\\[4pt]C_{e_{1}}^{s_{3}}=-\varepsilon _{1}^{2}\varepsilon _{2}^{3}&C_{e_{2}}^{s_{3}}=\varepsilon _{1}^{1}\varepsilon _{2}^{3}&C_{e_{3}}^{s_{2}}=-\varepsilon _{1}^{1}\varepsilon _{2}^{2}&C_{e_{4}}^{s_{2}}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}-\varepsilon _{1}^{1}\varepsilon _{2}^{3}+\varepsilon _{1}^{2}\varepsilon _{2}^{3}\end{array}}}$

Denominator:

${\displaystyle d=\varepsilon _{1}^{1}\varepsilon _{2}^{2}-\varepsilon _{1}^{1}\varepsilon _{2}^{3}+\varepsilon _{1}^{2}\varepsilon _{2}^{3}-\varepsilon _{2}^{1}\varepsilon _{1}^{2}}$

Assume that each of the following expressions is divided by ${\displaystyle d}$.

${\displaystyle {\begin{array}{lll}C_{e_{1}}^{J}=\varepsilon _{1}^{2}\varepsilon _{2}^{3}&C_{e_{2}}^{J}=-\varepsilon _{1}^{1}\varepsilon _{2}^{3}&C_{e_{3}}^{J}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}-\varepsilon _{2}^{1}\varepsilon _{1}^{2}\\[4pt]C_{e_{1}}^{s_{1}}=\varepsilon _{2}^{3}-\varepsilon _{2}^{2}&C_{e_{2}}^{s_{1}}=-\varepsilon _{2}^{3}-\varepsilon _{2}^{1}&C_{e_{3}}^{s_{1}}=\varepsilon _{2}^{2}-\varepsilon _{2}^{1}\\[4pt]C_{e_{1}}^{s_{2}}=\varepsilon _{1}^{2}&C_{e_{2}}^{s_{2}}=-\varepsilon _{1}^{1}&C_{e_{3}}^{s_{2}}=\varepsilon _{1}^{1}-\varepsilon _{1}^{2}\\[4pt]\end{array}}}$

Denominator:

${\displaystyle d=\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{4}-\varepsilon _{1}^{1}\varepsilon _{2}^{3}\varepsilon _{3}^{4}-\varepsilon _{3}^{1}\varepsilon _{1}^{2}\varepsilon _{2}^{3}+\varepsilon _{1}^{2}\varepsilon _{2}^{3}\varepsilon _{3}^{4}-\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{3}}$

Assume that each of the following expressions is divided by ${\displaystyle d}$.

${\displaystyle {\begin{array}{llll}C_{v_{1}}^{J}=\varepsilon _{1}^{2}\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{v_{2}}^{J}=-\varepsilon _{1}^{1}\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{v_{3}}^{J}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{4}&C_{v_{4}}^{J}=-\varepsilon _{1}^{1}\varepsilon _{2}^{2}\varepsilon _{3}^{3}-\varepsilon _{3}^{1}\varepsilon _{1}^{2}\varepsilon _{2}^{3}\\C_{v_{1}}^{S_{1}}=\varepsilon _{2}^{2}\varepsilon _{3}^{3}-\varepsilon _{2}^{2}\varepsilon _{3}^{4}+\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{v_{2}}^{S_{1}}=\varepsilon _{3}^{1}\varepsilon _{2}^{3}-\varepsilon _{2}^{3}\varepsilon _{3}^{4}&C_{v_{3}}^{S_{1}}=-\varepsilon _{3}^{1}\varepsilon _{2}^{2}+\varepsilon _{2}^{2}\varepsilon _{3}^{4}&C_{v_{4}}^{S_{1}}=\varepsilon _{3}^{1}\varepsilon _{2}^{2}-\varepsilon _{3}^{1}\varepsilon _{2}^{3}-\varepsilon _{2}^{2}\varepsilon _{3}^{3}\\C_{v_{1}}^{S_{2}}=-\varepsilon _{1}^{2}\varepsilon _{3}^{3}+\varepsilon _{1}^{2}\varepsilon _{3}^{4}&C_{v_{2}}^{S_{2}}=\varepsilon _{1}^{1}\varepsilon _{3}^{3}-\varepsilon _{1}^{1}\varepsilon _{3}^{4}&C_{v_{3}}^{S_{2}}=\varepsilon _{B}^{1}\varepsilon _{3}^{4}+\varepsilon _{3}^{1}\varepsilon _{1}^{2}-\varepsilon _{1}^{2}\varepsilon _{3}^{4}&C_{v_{4}}^{S_{2}}=-\varepsilon _{1}^{1}\varepsilon _{3}^{3}-\varepsilon _{3}^{1}\varepsilon _{1}^{2}+\varepsilon _{1}^{2}\varepsilon _{3}^{3}\\C_{v_{1}}^{S_{3}}=\varepsilon _{1}^{2}\varepsilon _{2}^{3}&C_{v_{2}}^{S_{3}}=-\varepsilon _{1}^{1}\varepsilon _{2}^{3}&C_{v_{3}}^{S_{3}}=\varepsilon _{1}^{1}\varepsilon _{2}^{2}&C_{v_{4}}^{S_{3}}=-\varepsilon _{1}^{1}\varepsilon _{2}^{2}+\varepsilon _{1}^{1}\varepsilon _{2}^{3}-\varepsilon _{1}^{2}\varepsilon _{2}^{3}\end{array}}}$

At steady-state ${\displaystyle v_{1}=v_{2}+v_{3}}$, therefore define the following two terms:

${\displaystyle \alpha ={\frac {v_{2}}{v_{1}}}\quad 1-\alpha ={\frac {v_{3}}{v_{1}}}}$

Denominator:

${\displaystyle d=\varepsilon _{s}^{2}\alpha +\varepsilon _{s}^{3}(1-\alpha )-\varepsilon _{s}^{1}}$

Assume that each of the following expressions is divided by ${\displaystyle d}$.

${\displaystyle {\begin{array}{lll}C_{e_{1}}^{J_{1}}=\varepsilon _{s}^{3}(1-\alpha )+\varepsilon _{s}^{2}\alpha \\C_{e_{1}}^{J_{1}}=-\varepsilon _{s}^{1}\alpha \\C_{e_{1}}^{J_{1}}=-\varepsilon _{s}^{1}(1-\alpha )+\varepsilon _{s}^{2}\alpha \end{array}}}$

• Metabolic control analysis
• Branched pathways
• Elasticity coefficient
• Biochemical systems theory
• Control coefficient (biochemistry)
• Flux (metabolism)