Cooperativity

When a substrate binds to one enzymatic subunit, the rest of the subunits are stimulated and become active. Ligands can either have positive cooperativity, negative cooperativity, or non-cooperativity.

 

An example of positive cooperativity is the binding of oxygen to hemoglobin. One oxygen molecule can bind to the ferrous iron of a heme molecule in each of the four chains of a hemoglobin molecule. Deoxy-hemoglobin has a relatively low affinity for oxygen, but when one molecule binds to a single heme, the oxygen affinity increases, allowing the second molecule to bind more easily, and the third and fourth even more easily. The oxygen affinity of 3-oxy-hemoglobin is ~300 times greater than that of deoxy-hemoglobin. This behavior leads the affinity curve of hemoglobin to be sigmoidal, rather than hyperbolic as with the monomeric myoglobin. By the same process, the ability for hemoglobin to lose oxygen increases as fewer oxygen molecules are bound.^[1] See also Oxygen-hemoglobin dissociation curve.

Negative cooperativity means that the opposite will be true; as ligands bind to the protein, the protein's affinity for the ligand will decrease, i.e. it becomes less likely for the ligand to bind to the protein. An example of this occurring is the relationship between glyceraldehyde-3-phosphate and the enzyme glyceraldehyde-3-phosphate dehydrogenase.

Homotropic cooperativity refers to the fact that the molecule causing the cooperativity is the one that will be affected by it. Heterotropic cooperativity is where a third party substance causes the change in affinity. Homotropic or heterotropic cooperativity could be of both positives as well as negative types depend upon whether it support or oppose further binding of the ligand molecules to the enzymes.^[3]

Cooperativity is not only a phenomenon of ligand binding, but also applies anytime energetic interactions make it easier or more difficult for something to happen involving multiple units as opposed to with single units. (That is, easier or more difficult compared with what is expected when only accounting for the addition of multiple units). For example, unwinding of DNA involves cooperativity: Portions of DNA must unwind in order for DNA to carry out replication, transcription and recombination. Positive cooperativity among adjacent DNA nucleotides makes it easier to unwind a whole group of adjacent nucleotides than it is to unwind the same number of nucleotides spread out along the DNA chain. The cooperative unit size is the number of adjacent bases that tend to unwind as a single unit due to the effects of positive cooperativity. This phenomenon applies to other types of chain molecules as well, such as the folding and unfolding of proteins and in the "melting" of phospholipid chains that make up the membranes of cells. Subunit cooperativity is measured on the relative scale known as Hill's Constant.

A simple and widely used model for molecular interactions is the Hill equation, which provides a way to quantify cooperative binding by describing the fraction of saturated ligand binding sites as a function of the ligand concentration.

The Hill coefficient is a measure of ultrasensitivity (i.e. how steep is the response curve).

From an operational point of view the Hill coefficient can be estimated as:

${\displaystyle n_{H}={\frac {\log(81)}{\ce {\log(EC90/EC10)}}}}$ .

where ${\displaystyle {\ce {EC90}}}$  and ${\displaystyle {\ce {EC10}}}$  are the input values needed to produce the 10% and 90% of the maximal response, respectively.

Global sensitivity measures such as the Hill coefficient do not characterise the local behaviours of the s-shaped curves. Instead, these features are well captured by the response coefficient measure ^[4] defined as:

${\displaystyle R(x)={\frac {x}{y}}{\frac {dy}{dx}}}$ 

In systems biology, such responses are referred to as elasticities.

Altszyler et al. (2017) have shown that these ultrasensitivity measures can be linked by the following equation:^[5]

${\displaystyle n_{H}=2{\frac {\int _{\ce {\log(EC10)}}^{\ce {\log(EC90)}}R_{f}(I)d(\log I)}{\ce {\log(EC90)-\log(EC10)}}}=2\langle R_{f}\rangle _{\ce {EC10,EC90}}}$ 

where ${\displaystyle \langle X\rangle _{a,b}}$  denoted the mean value of the variable x over the range [a,b].

Consider two coupled ultrasensitive modules, disregarding effects of sequestration of molecular components between layers. In this case, the expression for the system's dose-response curve, F, results from the mathematical composition of the functions, ${\displaystyle f_{i}}$ , which describe the input/output relationship of isolated modules ${\displaystyle i=1,2}$ :

${\displaystyle F(I_{1})=f_{2}{\big (}f_{1}(I_{1}){\big )}}$ 

Brown et al. (1997) ^[6][5] have shown that the local ultrasensitivity of the different layers combines multiplicatively:

${\displaystyle R(x)=R_{2}(f_{1}(x)).R_{1}(x)}$ .

In connection with this result, Ferrell et al. (1997) ^[7] showed, for Hill-type modules, that the overall cascade global ultrasensitivity had to be less than or equal to the product of the global ultrasensitivity estimations of each cascade's layer,^[5]

${\displaystyle n\leq n_{1}.n_{2}}$ ,

where ${\displaystyle n_{1}}$  and ${\displaystyle n_{2}}$  are the Hill coefficient of modules 1 and 2 respectively.

Altszyler et al. (2017) ^[5] have shown that the cascade's global ultrasensitivity can be analytically calculated:

${\displaystyle n=2\overbrace {\langle R_{2}\rangle _{X10_{2},X90_{2}}} ^{\nu _{2}}\overbrace {\langle R_{1}\rangle _{X10_{1},X90_{1}}} ^{\nu _{1}}=2\,\nu _{2}\,\nu _{1}}$ 

where ${\displaystyle X10_{i}}$  and ${\displaystyle X90_{i}}$  delimited the Hill input's working range of the composite system, i.e. the input values for the i-layer so that the last layer (corresponding to ${\displaystyle i=2}$  in this case) reached the 10% and 90% of it maximal output level. It followed this equation that the system's Hill coefficient n could be written as the product of two factors, ${\displaystyle \nu _{1}}$  and ${\displaystyle \nu _{2}}$ , which characterized local average sensitivities over the relevant input region for each layer: ${\displaystyle [X10_{i},X90_{i}]}$ , with ${\displaystyle i=1,2}$  in this case.

For the more general case of a cascade of N modules, the Hill coefficient can be expressed as:

${\displaystyle n=\nu _{N}\,\nu _{N-1}...\nu _{1}}$ ,

Several authors have reported the existence of supramultiplicative behavior in signaling cascades ^[8][9](i.e. the ultrasensitivity of the combination of layers is higher than the product of individual ultrasensitivities), but in many cases the ultimate origin of supramultiplicativity remained elusive. Altszyler et al. (2017)^[5] framework naturally suggested a general scenario where supramultiplicative behavior could take place. This could occur when, for a given module, the corresponding Hill's input working range was located in an input region with local ultrasensitivities higher than the global ultrasensitivity of the respective dose-response curve.