Elementary modes

Consider a simple branched pathway with all three steps irreversible. Such a pathway will admit two elementary modes which are indicated in thicked (or red) reaction lines.

A simple branch pathway with all reactions irreversible containing just two elementary modes

Because both ${\displaystyle v_{2}}$ and ${\displaystyle v_{3}}$ are irreversible, and elementary mode lying on both these reactions is not possible since it would mean one reactions going against its thermodynamic direction. Each mode in this system satisfies the three conditions described above. The first condition is steady state, that is for each mode ${\displaystyle \mathbf {v} _{i}}$, it has to be true that ${\displaystyle \mathbf {N} \mathbf {v} _{i}=0}$.

Algebraically the two modes are given by:

${\displaystyle \mathbf {v} _{1}={\begin{bmatrix}1\\1\\0\end{bmatrix}}\quad \mathbf {v} _{2}={\begin{bmatrix}1\\0\\1\end{bmatrix}}}$

By substituting each of these vectors into ${\displaystyle \mathbf {N} \mathbf {v} _{i}=0}$, it is easy to show that condition one is satisfied. For condition two we must ensure that all reactions that are irreversible have positive entries in the corresponding elements of the elementary modes. Since all three reactions in the branch are irreversible and all entries in the elementary modes are positive, condition two is satisfied.

Finally, to satisfy condition three, we must ask whether we can decompose the two elementary modes into other paths that can sustain a steady state while using the same non-zero entries in the elementary mode. In this example, it is impossible to decompose the elementary modes any further without disrupting the ability to sustain a steady state. Therefore, with all three conditions satisfied, we can conclude that the two vectors shown above are elementary modes.

All possible flows through a network can be constructed from linear combinations of the elementary modes, that is:

${\displaystyle {\begin{aligned}\mathbf {v} &=\sum \lambda _{i}\mathbf {v} _{i}\\{\mbox{where}}&\\\lambda &\geq 0\\\end{aligned}}}$

such that the entire space of flows through a network can be described. ${\displaystyle \lambda _{i}}$ must be greater than or equal to zero to ensure that irreversible steps aren't inadvertently made to go in the reverse direction. For example, the following is a possible steady-state flow in the branched pathway.

${\displaystyle v=2.5{\begin{bmatrix}1\\1\\0\end{bmatrix}}+0.5{\begin{bmatrix}1\\0\\1\end{bmatrix}}={\begin{bmatrix}3.0\\2.5\\0.5\end{bmatrix}}}$

If one of the outflow steps in the simple branched pathway is made reversible, an additional elementary mode becomes available, representing the flow between the two outflow branches. An additional mode emerges because, with only the first two modes, it is impossible to represent a flow between the two branches because the scaling factor, ${\displaystyle \lambda _{i}}$, cannot be negative (which would be required to reverse the flow).

The Wikipedia page Metabolic pathway defines a pathway as "a metabolic pathway is a linked series of chemical reactions occurring within a cell". This means that any sequence of reactions can be labeled a metabolic pathway. However, as metabolism was being uncovered, groups of reactions were assigned specific labels, such as glycolysis, Krebs Cycle, or Serine biosynthesis. Often the categorization was based on common chemistry or identification of an input and output. For example, serine biosynthesis starts at 3-phosphoglycerate and ends at serine. This is a somewhat ad hoc means for defining pathways, particularly when pathways are dynamic structures, changing as environmental result in changes in gene expression. For example, the Kreb Cycle is often not cyclic as depicted in textbooks. In E. coli and other bacteria, it is only cyclic during aerobic growth on acetate or fatty acids.[3] Instead, under anaerobiosis, its enzymes function as two distinct biosynthetic pathways producing succinyl-CoA and α-ketoglutarate.

It has therefore been proposed[4] to define a pathway as either a single elementary mode or some combination of elementary modes. The added advantage is that the set of elementary modes is unique and non-decomposable to simpler pathways. A single elementary mode can therefore be thought of as an elementary pathway. Note that the set of elementary modes will change as the set of expressed enzymes change during transitions from one cell state to another.

Elementary modes, therefore, provide an unambiguous definition of a pathway.

Condition three relates to the non-decomposability of an elementary mode and is partly what makes elementary modes interesting. The two other important features as indicated before are pathway uniqueness and thermodynamic plausibility. Decomposition implies that it is possible to represent a mode as a combination of two or more other modes. For example, a mode ${\displaystyle \mathbf {e} _{1}}$ might be composed from two other modes, ${\displaystyle \mathbf {e} _{2}}$ and ${\displaystyle \mathbf {e} _{3}}$:

${\displaystyle \mathbf {e} _{1}=\lambda _{1}\mathbf {e} _{2}+\lambda _{2}\mathbf {e} _{3}}$

If a mode can be decomposed, does it mean that the mode is not an elementary mode? Condition three provides a rule to determine whether a decomposition means that a given mode is an elementary mode or not. If it is only possible to decompose a given mode by introducing enzymes that are not used in the mode, then the mode is elementary. That is, is there more than one way to generate a pathway (i.e., something that can sustain a steady state) with the enzymes currently used in the mode? If so, then the mode is not elementary. To illustrate this subtle condition, consider the pathway shown in below.

A stylized pathway map of glycolysis.

This pathway represents a stylized rendition of glycolysis. Step three and six are reversible and correspond to triose phosphate isomerase and glycerol 3-phosphate dehydrogenase, respectively.

The network has four elementary flux modes, which are shown in the figure below.

A stylized pathway map of glycolysis.

The elementary flux mode vectors are shown below:

${\displaystyle {\begin{array}{l}\qquad \ \mathbf {e} _{1}\quad \ \ \mathbf {e} _{2}\quad \mathbf {e} _{3}\quad \mathbf {e} _{4}\\{\begin{bmatrix}&{\phantom {-}}0&{\phantom {-}}1&{\phantom {-}}1&{\phantom {-}}1\\&{\phantom {-}}0&{\phantom {-}}1&{\phantom {-}}1&{\phantom {-}}1\\&{\phantom {-}}1&-1&{\phantom {-}}1&{\phantom {-}}0\\&{\phantom {-}}1&{\phantom {-}}0&{\phantom {-}}2&{\phantom {-}}1\\&{\phantom {-}}1&{\phantom {-}}0&{\phantom {-}}2&{\phantom {-}}1\\&-1&{\phantom {-}}2&{\phantom {-}}0&{\phantom {-}}1\\\end{bmatrix}}\\\end{array}}}$

Note that it is possible to have negative entries in the set of elementary modes because they will correspond to the reversible steps. Of interest is the observation that the fourth vector, ${\displaystyle e_{4}=[1\ 1\ 0\ 1\ 1\ 1\ ]^{T}}$ (where ${\displaystyle T}$ represents the transpose) can be formed from the sum of the first and second vectors. This suggests that the fourth vector is not an elementary mode.

${\displaystyle \mathbf {e} _{4}{\begin{bmatrix}{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}0{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\\end{bmatrix}}=\mathbf {e} _{1}{\begin{bmatrix}{\phantom {-}}0{\phantom {-}}\\{\phantom {-}}0{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\-1{\phantom {-}}\\\end{bmatrix}}+\mathbf {e} _{2}{\begin{bmatrix}{\phantom {-}}1{\phantom {-}}\\{\phantom {-}}1{\phantom {-}}\\-1{\phantom {-}}\\{\phantom {-}}0{\phantom {-}}\\{\phantom {-}}0{\phantom {-}}\\{\phantom {-}}2{\phantom {-}}\\\end{bmatrix}}}$

However, this decomposition only works because we have introduced a new enzyme, ${\displaystyle E_{4}}$ (triose phosphate isomerase) which is not used in ${\displaystyle \mathbf {e} _{4}}$. It is, in fact impossible to decompose ${\displaystyle \mathbf {e} _{4}}$ into pathways that can sustain a steady state with only the five steps, ${\displaystyle e_{1},e_{2},e_{4},e_{5}\ {\mbox{and}}\ e_{6}}$, used in the elementary mode. We conclude therefore that ${\displaystyle \mathbf {e} _{4}}$ is an elementary mode.

• Computational systems biology
• Metabolic control analysis
• Metabolic network
• Metabolic pathway
• Metagenomics