Branched pathways

Common pattern in metabolism

Branched pathways, also known as branch points (not to be confused with the mathematical branch point), 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.

Branched pathways are present in numerous metabolic reactions, including glycolysis, the synthesis of lysine, glutamine, and penicillin,[1] and in the production of the aromatic amino acids.[2]

Simple Branch Pathway.  v 1 , v 2 {\displaystyle v_{1},v_{2}} and  v 3 {\displaystyle v_{3}} are the reaction rates for each arm of the branch.

In general, a single branch may have  b {\displaystyle b} producing branches and  d {\displaystyle d} consuming branches. If the intermediate at the branch point is given by  s i {\displaystyle s_{i}} , then the rate of change of  s i {\displaystyle s_{i}} is given by:

∑ i = 1 b v i − ∑ j = 1 d v j = d s i d t {\displaystyle \sum _{i=1}^{b}v_{i}-\sum _{j=1}^{d}v_{j}={\frac {ds_{i}}{dt}}}

At steady-state when  d s i / d t = 0 {\displaystyle ds_{i}/dt=0} the consumption and production rates must be equal:

∑ i = 1 b v i = ∑ j = 1 d v j {\displaystyle \sum _{i=1}^{b}v_{i}=\sum _{j=1}^{d}v_{j}}

Biochemical pathways can be investigated by computer simulation or by looking at the sensitivities, i.e. control coefficients for flux and species concentrations using metabolic control analysis.

## Elementary properties

[edit]

A simple branched pathway has one key property related to the conservation of mass. In general, the rate of change of the branch species based on the above figure is given by:

d s 1 d t = v 1 − ( v 2 + v 3 ) {\displaystyle {\frac {ds_{1}}{dt}}=v_{1}-(v_{2}+v_{3})}

At steady-state the rate of change of  S 1 {\displaystyle S_{1}} is zero. This gives rise to a steady-state constraint among the branch reaction rates:

v 1 = v 2 + v 3 {\displaystyle v_{1}=v_{2}+v_{3}}

Such constraints are key to computational methods such as flux balance analysis.

## Control properties of a branch pathway

[edit]

Branched pathways have unique control properties compared to simple linear chain or cyclic pathways. These properties can be investigated using metabolic control analysis. The fluxes can be controlled by enzyme concentrations  e 1 {\displaystyle e_{1}} ,  e 2 {\displaystyle e_{2}} , and  e 3 {\displaystyle e_{3}} respectively, described by the corresponding flux control coefficients. To do this the flux control coefficients with respect to one of the branch fluxes can be derived. The derivation is shown in a subsequent section. The flux control coefficient with respect to the upper branch flux,  J 2 {\displaystyle J_{2}} are given by:

C e 1 J 2 = ε 2 ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{1}}^{J_{2}}={\frac {\varepsilon _{2}}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 2 J 2 = ε 3 ( 1 − α ) − ε 1 ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{2}}^{J_{2}}={\frac {\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 3 J 2 = − ε 2 ( 1 − α ) ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{3}}^{J_{2}}={\frac {-\varepsilon _{2}(1-\alpha )}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

where  α {\displaystyle \alpha } is the fraction of flux going through the upper arm,  J 2 {\displaystyle J_{2}} , and  1 − α {\displaystyle 1-\alpha } the fraction going through the lower arm,  J 3 {\displaystyle J_{3}} .  ε 1 , ε 2 , {\displaystyle \varepsilon _{1},\varepsilon _{2},} and  ε 3 {\displaystyle \varepsilon _{3}} are the elasticities for  s 1 {\displaystyle s_{1}} with respect to  v 1 , v 2 , {\displaystyle v_{1},v_{2},} and  v 3 {\displaystyle v_{3}} respectively.

For the following analysis, the flux  J 2 {\displaystyle J_{2}} will be the observed variable in response to changes in enzyme concentrations.

There are two possible extremes to consider, either most of the flux goes through the upper branch  J 2 {\displaystyle J_{2}} or most of the flux goes through the lower branch,  J 3 {\displaystyle J_{3}} . The former, depicted in panel a), is the least interesting as it converts the branch in to a simple linear pathway. Of more interest is when most of the flux goes through  J 3 {\displaystyle J_{3}}

If most of the flux goes through  J 3 {\displaystyle J_{3}} , then  α → 0 {\displaystyle \alpha \rightarrow 0} and  1 − α → 1 {\displaystyle 1-\alpha \rightarrow 1} (condition (b) in the figure), the flux control coefficients for  J 2 {\displaystyle J_{2}} with respect to  e 2 {\displaystyle e_{2}} and  e 3 {\displaystyle e_{3}} can be written:

C e 2 J 2 → 1 {\displaystyle C_{e_{2}}^{J_{2}}\rightarrow 1}

C e 3 J 2 → ε 2 ε 1 − ε 3 {\displaystyle C_{e_{3}}^{J_{2}}\rightarrow {\frac {\varepsilon _{2}}{\varepsilon _{1}-\varepsilon _{3}}}}
Changes in flux control depending on whether the flux goes through the upper or lower branches. The system output is J2. If most of the flux goes through J2(a) then the pathway behaves like a simple linear change, where flux control on J3 is negligible and control is shared between J1 and J2. The other extreme is when most of the flux goes through J3 (b). This makes J2 highly sensitive to changes in J1 and J3 resulting is very high flux control, often greater than 1.0. Under these conditions the flux control at J3 is also negative since J3 can siphon off flux from J2.

That is,  e 2 {\displaystyle e_{2}} acquires proportional influence over its own flux,  J 2 {\displaystyle J_{2}} . Since  J 2 {\displaystyle J_{2}} only carries a very small amount of flux, any changes in  e 2 {\displaystyle e_{2}} will have little effect on  S {\displaystyle S} . Hence the flux through  e 2 {\displaystyle e_{2}} is almost entirely governed by the activity of  e 2 {\displaystyle e_{2}} . Because of the flux summation theorem and the fact that  C e 2 J 2 = 1 {\displaystyle C_{e_{2}}^{J_{2}}=1} , it means that the remaining two coefficients must be equal and opposite in value. Since  C e 1 J 2 {\displaystyle C_{e_{1}}^{J_{2}}} is positive,  C e 3 J 2 {\displaystyle C_{e_{3}}^{J_{2}}} must be negative. This also means that in this situation, there can be more than one Rate-limiting step (biochemistry) in a pathway.

Unlike a linear pathway, values for  C e 3 J 2 {\displaystyle C_{e_{3}}^{J_{2}}} and  C e 1 J 2 {\displaystyle C_{e_{1}}^{J_{2}}} are not bounded between zero and one. Depending on the values of the elasticities, it is possible for the control coefficients in a branched system to greatly exceed one.[3] This has been termed the branchpoint effect by some in the literature.[4]

### Example

[edit]

Flux control coefficients in a branched pathway where most flux goes through  J 3 {\displaystyle J_{3}} . Note that step 2 has almost proportional control over  J 2 {\displaystyle J_{2}} while steps 1 and 3 show greater than proportional control over  J 2 {\displaystyle J_{2}} .

The following branch pathway model (in antimony format) illustrates the case  J 1 {\displaystyle J_{1}} and  J 3 {\displaystyle J_{3}} have very high flux control and step J2 has proportional control.

J1: $Xo -> S1; e1*k1*Xo
J2:  S1 ->; e2*k3*S1/(Km1 + S1)
J3:  S1 ->; e3*k4*S1/(Km2 + S1)

k1 = 2.5;
k3 = 5.9; k4 = 20.75
Km1 = 4; Km2 = 0.02
Xo =5;
e1 = 1; e2 = 1; e3 = 1

A simulation of this model yields the following values for the flux control coefficients with respect to flux  J 2 {\displaystyle J_{2}}

## Branch point theorems

[edit]

In a linear pathway, only two sets of theorems exist, the summation and connectivity theorems. Branched pathways have an additional set of branch centric summation theorems. When combined with the connectivity theorems and the summation theorem, it is possible to derive the control equations shown in the previous section. The deviation of the branch point theorems is as follows.

1. Define the fractional flux through  J 2 {\displaystyle J_{2}} and  J 3 {\displaystyle J_{3}} as  α = J 2 / J 1 {\displaystyle \alpha =J_{2}/J_{1}} and  1 − α = J 3 / J 1 {\displaystyle 1-\alpha =J_{3}/J_{1}} respectively.
2. Increase  e 2 {\displaystyle e_{2}} by  δ e 2 {\displaystyle \delta e_{2}} . This will decrease  S 1 {\displaystyle S_{1}} and increase  J 1 {\displaystyle J_{1}} through relief of product inhibition.[5]
3. Make a compensatory change in  J 1 {\displaystyle J_{1}} by decreasing  e 1 {\displaystyle e_{1}} such that  S 1 {\displaystyle S_{1}} is restored to its original concentration (hence  δ S 1 = 0 {\displaystyle \delta S_{1}=0} ).
4. Since  e 1 {\displaystyle e_{1}} and  S 1 {\displaystyle S_{1}} have not changed,  δ J 1 = 0 {\displaystyle \delta J_{1}=0} .

Following these assumptions two sets of equations can be derived: the flux branch point theorems and the concentration branch point theorems.[6]

### Derivation

[edit]

From these assumptions, the following system equation can be produced:

C e 2 J 1 δ e 2 e 2 + C e 3 J 1 δ e 3 e 3 = δ J 1 J 1 = 0 {\displaystyle C_{e_{2}}^{J_{1}}{\frac {\delta e_{2}}{e_{2}}}+C_{e_{3}}^{J_{1}}{\frac {\delta e_{3}}{e_{3}}}={\frac {\delta J_{1}}{J_{1}}}=0}

Because  δ S 1 = 0 {\displaystyle \delta S_{1}=0} and, assuming that the flux rates are directly related to the enzyme concentration thus, the elasticities,  ε e i v {\displaystyle \varepsilon _{e_{i}}^{v}} , equal one, the local equations are:

δ v 2 v 2 = δ e 2 e 2 {\displaystyle {\frac {\delta v_{2}}{v_{2}}}={\frac {\delta e_{2}}{e_{2}}}}

δ v 3 v 3 = δ e 3 e 3 {\displaystyle {\frac {\delta v_{3}}{v_{3}}}={\frac {\delta e_{3}}{e_{3}}}}

Substituting  δ v i v i {\displaystyle {\frac {\delta v_{i}}{v_{i}}}} for  δ e i e i {\displaystyle {\frac {\delta e_{i}}{e_{i}}}} in the system equation results in:

C e 2 J 1 δ v 2 v 2 + C e 3 J 1 δ v 3 v 3 = 0 {\displaystyle C_{e_{2}}^{J_{1}}{\frac {\delta v_{2}}{v_{2}}}+C_{e_{3}}^{J_{1}}{\frac {\delta v_{3}}{v_{3}}}=0}

Conservation of mass dictates  δ J 1 = δ J 2 + δ J 3 {\displaystyle \delta J_{1}=\delta J_{2}+\delta J_{3}} since  δ J 1 = 0 {\displaystyle \delta J_{1}=0} then  δ v 2 = − δ v 3 {\displaystyle \delta v_{2}=-\delta v_{3}} . Substitution eliminates the  δ v 3 {\displaystyle \delta v_{3}} term from the system equation:

C e 2 J 1 δ v 2 v 2 − C e 3 J 1 δ v 2 v 3 = 0 {\displaystyle C_{e_{2}}^{J_{1}}{\frac {\delta v_{2}}{v_{2}}}-C_{e_{3}}^{J_{1}}{\frac {\delta v_{2}}{v_{3}}}=0}

Dividing out  δ v 2 v 2 {\displaystyle {\frac {\delta v_{2}}{v_{2}}}} results in:

C e 2 J 1 − C e 3 J 1 v 2 v 3 = 0 {\displaystyle C_{e_{2}}^{J_{1}}-C_{e_{3}}^{J_{1}}{\frac {v_{2}}{v_{3}}}=0}

v 2 {\displaystyle v_{2}} and  v 3 {\displaystyle v_{3}} can be substituted by the fractional rates giving:

C e 2 J 1 − C e 3 J 1 α 1 − α = 0 {\displaystyle C_{e_{2}}^{J_{1}}-C_{e_{3}}^{J_{1}}{\frac {\alpha }{1-\alpha }}=0}

Rearrangement yields the final form of the first flux branch point theorem:[6]

C e 2 J 1 ( 1 − α ) − C e 3 J 1 α = 0 {\displaystyle C_{e_{2}}^{J_{1}}(1-\alpha )-C_{e_{3}}^{J_{1}}{\alpha }=0}

Similar derivations result in two more flux branch point theorems and the three concentration branch point theorems.

### Flux branch point theorems

[edit]

C e 2 J 1 ( 1 − α ) − C e 3 J 1 ( α ) = 0 {\displaystyle C_{e_{2}}^{J_{1}}(1-\alpha )-C_{e_{3}}^{J_{1}}(\alpha )=0}

C e 1 J 2 ( 1 − α ) + C e 3 J 2 ( α ) = 0 {\displaystyle C_{e_{1}}^{J_{2}}(1-\alpha )+C_{e_{3}}^{J_{2}}(\alpha )=0}

C e 1 J 3 ( α ) + C e 2 J 3 = 0 {\displaystyle C_{e_{1}}^{J_{3}}(\alpha )+C_{e_{2}}^{J_{3}}=0}

### Concentration branch point theorems

[edit]

C e 2 S 1 ( 1 − α ) + C e 3 S 1 ( α ) = 0 {\displaystyle C_{e_{2}}^{S_{1}}(1-\alpha )+C_{e_{3}}^{S_{1}}(\alpha )=0}

C e 1 S 1 ( 1 − α ) + C e 3 S 1 = 0 {\displaystyle C_{e_{1}}^{S_{1}}(1-\alpha )+C_{e_{3}}^{S_{1}}=0}

C e 1 S 1 ( α ) + C e 2 S 1 = 0 {\displaystyle C_{e_{1}}^{S_{1}}(\alpha )+C_{e_{2}}^{S_{1}}=0}

Following the flux summation theorem[7] and the connectivity theorem[8] the following system of equations can be produced for the simple pathway.[9]

C e 1 J 1 + C e 2 J 1 + C e 3 J 1 = 1 {\displaystyle C_{e_{1}}^{J_{1}}+C_{e_{2}}^{J_{1}}+C_{e_{3}}^{J_{1}}=1}

C e 1 J 2 + C e 2 J 2 + C e 3 J 2 = 1 {\displaystyle C_{e_{1}}^{J_{2}}+C_{e_{2}}^{J_{2}}+C_{e_{3}}^{J_{2}}=1}

C e 1 J 3 + C e 2 J 3 + C e 3 J 3 = 1 {\displaystyle C_{e_{1}}^{J_{3}}+C_{e_{2}}^{J_{3}}+C_{e_{3}}^{J_{3}}=1}

C e 1 J 1 ε s v 1 + C e 2 J 1 ε s v 2 + C e 3 J 1 ε s v 3 = 0 {\displaystyle C_{e_{1}}^{J_{1}}\varepsilon _{s}^{v_{1}}+C_{e_{2}}^{J_{1}}\varepsilon _{s}^{v_{2}}+C_{e_{3}}^{J_{1}}\varepsilon _{s}^{v_{3}}=0}

C e 1 J 2 ε s v 1 + C e 2 J 2 ε s v 2 + C e 3 J 2 ε s v 3 = 0 {\displaystyle C_{e_{1}}^{J_{2}}\varepsilon _{s}^{v_{1}}+C_{e_{2}}^{J_{2}}\varepsilon _{s}^{v_{2}}+C_{e_{3}}^{J_{2}}\varepsilon _{s}^{v_{3}}=0}

C e 1 J 3 ε s v 1 + C e 2 J 3 ε s v 2 + C e 3 J 3 ε s v 3 = 0 {\displaystyle C_{e_{1}}^{J_{3}}\varepsilon _{s}^{v_{1}}+C_{e_{2}}^{J_{3}}\varepsilon _{s}^{v_{2}}+C_{e_{3}}^{J_{3}}\varepsilon _{s}^{v_{3}}=0}

Using these theorems plus flux summation and connectivity theorems values for the concentration and flux control coefficients can be determined using linear algebra.[6]

C e 1 J 2 = ε 2 ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{1}}^{J_{2}}={\frac {\varepsilon _{2}}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 2 J 2 = ε 3 ( 1 − α ) − ε 1 ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{2}}^{J_{2}}={\frac {\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 3 J 2 = − ε 2 ( 1 − α ) ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{3}}^{J_{2}}={\frac {-\varepsilon _{2}(1-\alpha )}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 1 S 1 = 1 ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{1}}^{S_{1}}={\frac {1}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 1 S 1 = − α ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{1}}^{S_{1}}={\frac {-\alpha }{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

C e 3 S 1 = − ( 1 − α ) ε 2 α + ε 3 ( 1 − α ) − ε 1 {\displaystyle C_{e_{3}}^{S_{1}}={\frac {-(1-\alpha )}{\varepsilon _{2}\alpha +\varepsilon _{3}(1-\alpha )-\varepsilon _{1}}}}

## See also

[edit]

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

## References

[edit]

1. ^ Heijnen, J. J.; van Gulik, W. M.; Shimizu, H.; Stephanopoulos, G. (2004-10-01). "Metabolic flux control analysis of branch points: an improved approach to obtain flux control coefficients from large perturbation data". Metabolic Engineering. 6 (4): 391–400\. doi:10.1016/j.ymben.2004.07.002. ISSN 1096-7176. PMID 15491867.
2. ^ "W_2022_Bis2a_Igo_Reading_15". Biology LibreTexts. 2021-12-13. Retrieved 2022-12-15.
3. ^ Kacser, H. (1 January 1983). "The control of enzyme systems in vivo : Elasticity analysis of the steady state". Biochemical Society Transactions. 11 (1): 35–40\. doi:10.1042/bst0110035. PMID 6825913.
4. ^ LaPorte, D C; Walsh, K; Koshland, D E (November 1984). "The branch point effect. Ultrasensitivity and subsensitivity to metabolic control". Journal of Biological Chemistry. 259 (22): 14068–14075\. doi:10.1016/S0021-9258(18)89857-X.
5. ^ Liu, Yan; Zhang, Fan; Jiang, Ling; Perry, J. Jefferson P.; Zhao, Zhihe; Liao, Jiayu (2021-12-15). "Product inhibition kinetics determinations - Substrate interaction affinity and enzymatic kinetics using one quantitative FRET assay". International Journal of Biological Macromolecules. 193 (Pt B): 1481–1487\. doi:10.1016/j.ijbiomac.2021.10.211. ISSN 0141-8130. PMID 34780893. S2CID 244107621.
6. ^ a b c Sauro, Herbert (2018). Systems Biology: An Introduction to Metabolic Control Analysis (1st ed.). Ambrosius Publishing. pp. 115–122\. ISBN 978-0-9824773-6-6.
7. ^ Agutter, Paul S. (2008-10-21). "The flux-summation theorem and the 'evolution of dominance'". Journal of Theoretical Biology. 254 (4): 821–825\. doi:10.1016/j.jtbi.2008.07.027. ISSN 0022-5193. PMID 18706429.
8. ^ Kacser, H.; Burns, J. A. (1973). "The control of flux". Symposia of the Society for Experimental Biology. 27: 65–104\. ISSN 0081-1386. PMID 4148886.
9. ^ Fell, David A.; Sauro, Herbert M. (1985). "Metabolic control and its analysis. Additional relationships between elasticities and control coefficients". European Journal of Biochemistry. 148 (3) (published May 1985): 555–561\. doi:10.1111/j.1432-1033.1985.tb08876.x. ISSN 0014-2956. PMID 3996393.