Biology:Connectivity theorems

From HandWiki

The stoichiometric structure and mass-conservation properties of biochemical pathways gives rise to a series of theorems or relationships between the control coefficients and the control coefficients and elasticities. There are a large number of such relationships depending on the pathway configuration (e.g. linear, branched or cyclic) which have been documented and discovered by various authors. The term theorem has been used to describe these relationships because they can be proved in terms of more elementary concepts. The operational proofs[1] in particular are of this nature.

The most well known of these theorems are the summation theorems for the control coefficients and the connectivity theorems which relate control coefficients to the elasticities. The focus of this page are the connectivity theorems.

When deriving the summation theorems, a thought experiment was conducted that involved manipulating enzyme activities such that concentrations were unaffected but fluxes changed. The connectivity theorems use the opposite thought experiment, that is enzyme activities are changed such that concentrations change but fluxes are unchanged.[1] This is an important observation that highlights the orthogonal nature of these two sets of theorem.[2]

As with the summation theorems, the connectivity theorems can also be proved using more rigorous mathematical approaches involving calculus and linear algebra.[3][4][5] Here the more intuitive and operational proofs will be used to prove the connectivity theorems.

Statement of the connectivity theorems

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 [math]\displaystyle{ S_n }[/math] is different from the local species [math]\displaystyle{ S_m }[/math].

[math]\displaystyle{ \sum_i C^J_i \varepsilon^i_s = 0 }[/math]
[math]\displaystyle{ \sum_i C^{s_n}_i \varepsilon^i_{s_m} = 0 \quad n \neq m }[/math]
[math]\displaystyle{ \sum_i C^{s_n}_i \varepsilon^i_{s_m} = -1 \quad n = m }[/math]

Proof

The operational proof for the flux connectivity theorem[1] relies on making perturbations to enzyme levels such that the pathway flux is unchanged but a single metabolite level is changed. This can be illustrated with the following pathway:

[math]\displaystyle{ \stackrel{v_1}{\longrightarrow} S_1 \stackrel{v_2}{\longrightarrow} S_2 \stackrel{v_3}{\longrightarrow} S_3 \stackrel{v_4}{\longrightarrow} }[/math]

Let us make a change to the rate through [math]\displaystyle{ v_2 }[/math] by increasing the concentration of enzyme [math]\displaystyle{ e_2 }[/math]. Assume [math]\displaystyle{ e_2 }[/math] is increased by an amount, [math]\displaystyle{ \delta e_2 }[/math]. This will result in a change to the steady-state of the pathway. The concentrations of [math]\displaystyle{ s_2, s_3 }[/math], and the flux, [math]\displaystyle{ J }[/math] through the pathway will increase, and the concentration of [math]\displaystyle{ s_1 }[/math] will decrease because it is upstream of the disturbance.

Impose a second change to the pathway such that the flux, [math]\displaystyle{ J }[/math] is restored to what it was before the original change. Since the flux increased when [math]\displaystyle{ e_2 }[/math] was changed, the flux can be decreased by decreasing one of the other enzyme levels. If the concentration of [math]\displaystyle{ e_3 }[/math] is decreased, this will reduce the flux. Decreasing [math]\displaystyle{ e_3 }[/math] will also cause the concentration of [math]\displaystyle{ s_2 }[/math] to further increase. However, [math]\displaystyle{ s_1 }[/math] and [math]\displaystyle{ s_3 }[/math] will change in the opposite direction compared to when [math]\displaystyle{ e_2 }[/math] was increased.

When [math]\displaystyle{ e_3 }[/math] is sufficiently changed so that the flux is restored to its original value, the concentrations of [math]\displaystyle{ s_1 }[/math] and [math]\displaystyle{ s_3 }[/math] will also be restored to their original values. It is only [math]\displaystyle{ s_2 }[/math] that will differ. This is true because the flux through [math]\displaystyle{ v_1 }[/math] is now the same as it was originally (since we’ve restored the flux), and [math]\displaystyle{ e_1 }[/math] has not been manipulated in anyway. This means that the concentration of [math]\displaystyle{ s_1 }[/math] and all species upstream of [math]\displaystyle{ s_1 }[/math] must be the same as they were before the modulations occurred. The same arguments apply to [math]\displaystyle{ s_3 }[/math] and all species downstream of [math]\displaystyle{ v_4 }[/math].

The net result is that [math]\displaystyle{ e_2 }[/math] has been increased by [math]\displaystyle{ \delta e_2 }[/math] resuling a change in flux of [math]\displaystyle{ \delta J }[/math]. The concentration of [math]\displaystyle{ e_3 }[/math] was decreased such that the flux was restored to it original value, [math]\displaystyle{ \delta J = 0 }[/math]. In the process, [math]\displaystyle{ s_2 }[/math] changed by [math]\displaystyle{ \delta s_2 }[/math] but neither [math]\displaystyle{ s_1 }[/math] or [math]\displaystyle{ s_3 }[/math]. In fact no other species in the entire system has changed other than [math]\displaystyle{ s_2 }[/math].

Perturbations involved in the connectivity theorem. a) First perturbation (orange) at step 2; b) Second counter perturbation at step 3, resulting in no change in flux, J but a change to species S2. No other species changes

This thought experiment can be expressed mathematically as follows. The system equations in terms of the flux control coefficients can be written as:

[math]\displaystyle{ \frac{\delta J}{J} = 0 = C^J_2 \frac{\delta e_2}{e_2} + C^J_3 \frac{\delta e_3}{e_3} }[/math]

There are only two terms because only [math]\displaystyle{ e_2 }[/math] and [math]\displaystyle{ e_3 }[/math] were changed.

The local change at each step can be written for [math]\displaystyle{ v_2 }[/math] and [math]\displaystyle{ v_2 }[/math] in terms of elasticities:

[math]\displaystyle{ 0 = \frac{\delta v_2}{v_2} = \frac{\delta e_2}{e_2} + \varepsilon^2_2 \frac{\delta s_2}{s_2} }[/math]

[math]\displaystyle{ 0 = \frac{\delta v_3}{v_3} = \frac{\delta e_3}{e_3} + \varepsilon^3_2 \frac{\delta s_2}{s_2} }[/math]

Note that [math]\displaystyle{ \delta e_2/e_2 }[/math] won't necessarily equal [math]\displaystyle{ \delta e_2/e_3 }[/math] and by construction both rates, [math]\displaystyle{ v_2 }[/math] and [math]\displaystyle{ v_3 }[/math] showed no change. Also by construction only [math]\displaystyle{ s_2 }[/math] changed.

The local equation can be rearranged as:

[math]\displaystyle{ \frac{\delta e_2}{e_2}=-\varepsilon_2^2 \frac{\delta s_2}{s_2} }[/math]

[math]\displaystyle{ \frac{\delta e_3}{e_3}=-\varepsilon_2^3 \frac{\delta s_2}{s_2} }[/math]

The right-hand sides can be inserted into the system equation the change in flux:

[math]\displaystyle{ 0=\frac{\delta J}{J}=-\left(C_{e_2}^J \varepsilon_2^2 \frac{\delta s_2}{s_2}+C_{e_3}^J \varepsilon_2^3 \frac{\delta s_2}{s_2}\right) }[/math]

Therefore:

[math]\displaystyle{ 0=\frac{\delta s_2}{s_2}\left(C_{e_2}^J \varepsilon_2^2+C_{e_3}^J \varepsilon_2^3\right) }[/math]

However, by construction of the perturbations, [math]\displaystyle{ \delta s_2/s_2 }[/math] does not equal zero, hence we arrive at the connectivity theorem:

[math]\displaystyle{ 0=C_{e_2}^J \varepsilon_2^2+C_{e_3}^J \varepsilon_2^3 }[/math]

The operational method can also be used for systems where a given metabolite can influence multiple steps. This would apply to cases such as branched systems or systems with negative feedback loops.

The same approach can be used to derive the concentration connectivity theorems except one can consider either the case that focuses on a single species or a second case where the system equation is written to consider the effect on a distance species.

Interpretation

The flux control coefficient connectivity theorem is the easiest to understand. Starting with a simple two step pathway:

[math]\displaystyle{ X_o \stackrel{v_1}{\longrightarrow} S_1 \stackrel{v_2}{\longrightarrow} X_1 }[/math]

where [math]\displaystyle{ X_o }[/math] and [math]\displaystyle{ X_1 }[/math] are fixed species so that the pathway can reach a steady-state. [math]\displaystyle{ v_1 }[/math] and [math]\displaystyle{ v_2 }[/math] are the reaction rates for the first and second steps.

We can write the flux connectivity theorem for this simple system as follows:

[math]\displaystyle{ C^J_1 \varepsilon^1_1 + C^J_2 \varepsilon^2_1 = 0 }[/math]

where [math]\displaystyle{ \varepsilon^1_1 }[/math] is the elasticity of the first step [math]\displaystyle{ v_1 }[/math] with respect to the species [math]\displaystyle{ S_1 }[/math] and [math]\displaystyle{ \varepsilon^2_1 }[/math] is the elasticity of the second step [math]\displaystyle{ v_2 }[/math] with respect to the species [math]\displaystyle{ S_1 }[/math]. It is easier to interpret the equation with a slight rearrangement to the following form:

[math]\displaystyle{ \frac{C^J_1}{C^J_2} = -\frac{\varepsilon^2_1}{\varepsilon^1_1} }[/math]

The equation indicates that the ratio of the flux control coefficients is inversely proportional to the elasticities. That is, a high flux control coefficient on step one is associated with a low elasticity [math]\displaystyle{ \varepsilon^1_1 }[/math] and vice versa. Likewise a high value for the flux control coefficient on step two is associated with a low elasticity [math]\displaystyle{ \varepsilon^2_1 }[/math].

This can be explained as follows: If [math]\displaystyle{ \varepsilon^1_1 }[/math] is high (in absolute terms, since it is negative) then a change at [math]\displaystyle{ v_1 }[/math] will be resisted by the elasticity, hence the flux control coefficient on step one will be low.

See also

References

  1. 1.0 1.1 1.2 Kacser, H.; Burns, J. A. (1973). "The control of flux". Symposia of the Society for Experimental Biology 27: 65–104. PMID 4148886. 
  2. Mazat, J P; Letellier, T; Reder, C (1990). "Metabolic control theory: the geometry of the triangle". Biomed Biochim Acta 49 (8–9): 801–810. PMID 2082923. 
  3. Heinrich, Reinhart; Schuster, Stefan (31 August 1996) (in English). The Regulation of Cellular Systems (1996th ed.). Springer. ISBN 978-0-412-03261-5. 
  4. Hofmeyr, Jan-Hendrik. "Metabolic control analysis in a nutshell". Proceedings of the 2nd International Conference on Systems Biology. https://www.academia.edu/4610732. 
  5. Mazat, Jean-Pierre; Reder, Christine; Letellier, T. (1990). "The Structural Approach to Metabolic Control Analysis II: Geometrical Aspects" (in en). Control of Metabolic Processes. Springer US. pp. 129–138. doi:10.1007/978-1-4757-9856-2_9. ISBN 978-1-4757-9856-2. https://link.springer.com/chapter/10.1007/978-1-4757-9856-2_9.