help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Biophys J. BioFAST on April 15, 2005.
doi:10.1529/biophysj.105.060194
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.060194v1
89/1/76    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Abeliovich, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Abeliovich, H.
Biophysical Journal 89:76-79 (2005)
© 2005 The Biophysical Society

An Empirical Extremum Principle for the Hill Coefficient in Ligand-Protein Interactions Showing Negative Cooperativity

Hagai Abeliovich

Department of Biochemistry and Food Science, Hebrew University of Jerusalem, Rehovot, Israel

Correspondence: Address reprint requests to H. Abeliovich, Dept. of Biochemistry and Food Science, Hebrew University of Jerusalem, Rehovot, Israel 76100. Tel.: 972-8-948-9060; E-mail: ahagai{at}agri.huji.ac.il.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The Hill coefficient (nH) is a central parameter in the study of ligand-protein interactions, which measures the degree of cooperativity between subunits that bind the ligand in multisubunit proteins. The most common usage of nH is as an estimate of the minimal number of interacting binding sites in positively cooperating systems. In the present study, a statistical interpretation of nH for a generalized system of multiple identical binding sites is developed. This interpretation is then applied to the derivation of an empirical extremum principle for nH in negatively cooperating systems of identical binding sites, which can be used for the estimation of the minimal number of interacting sites in such systems.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Dating back to the studies by Hill (1Go) on oxygen binding to hemoglobin, the phenomenon of cooperative binding of small molecule ligands to multimeric proteins has held a special role in the understanding of physiological regulation of protein and enzyme function. For systems exhibiting positive cooperativity, it has been demonstrated that the Hill coefficient, nH, places a lower limit for the number of interacting sites on the protein (2Go,3Go). Comprehensive theoretical studies of various models of cooperative and allosteric systems have been published, as well as more generalized treatises. In particular, J. Wyman and various colleagues have contributed much to the analysis of the binding curves of cooperative systems (4Go–6Go). Such analyses are based on the general partition function dictated by the Adair equation, and show that nH is related to the variance of the binding number (the average occupancy of the protein taken over the whole ensemble). These studies proved mathematically that nH is a lower limit on the number of interacting sites, in positively cooperating systems (3Go). This extremum principle is highly useful in that it allows experimentalists to estimate the minimum number of binding sites on the protein and to compare this result with direct physical measurements. Negative cooperativity, a phenomenon in which binding at one site interferes with binding to additional, identical sites on the protein, is less well understood (7Go–11Go). Specifically an extremum principle is unknown for such systems, and there is no method to estimate the minimal number of interacting binding sites. Indeed, under the assumptions of the Koshland-Nemethy-Filmer model, the minimal theoretical value of nH is 0. However, nH is inherently an empirical entity, and the experimental conditions impose additional constraints on its values. The identification of an empirical lower limit for nH is an important issue since one of the main stumbling blocks in the identification of negatively cooperating systems is the fact that a mixture of nonidentical binding sites will also yield an nH value <1, and unless comprehensive additional studies are undertaken to prove that the binding sites are a priori identical, it is not possible to identify a system as a negatively cooperating system using only the fact that nH < 1 (12Go). Here, I have extended previous analyses by demonstrating that nH can be described wholly in statistical terms as the ratio of the variance of the binding number to the variance of the binding number in a reference system consisting of an identical number of noninteracting binding sites. This additional generalization allows the derivation of an analogous extremum principle for systems with negative cooperative interactions between binding sites, showing that the minimum empirical value of nH in a negatively cooperating system of N identical binding sites is 1/N, assuming that one uses the empirical saturation point as the absolute saturation point, without reference to the theoretical saturation capacity of the system. Since such a restriction does not hold for a system of N independent nonidentical binding sites, this finding provides an additional tool for distinguishing bona fide negative cooperativity from a mixture of binding sites.


    THEORY
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Consider a macromolecule P, which contains N identical sites for ligand L. We define the general (macroscopic) formation constants Kn for the nth ligated species of the macromolecule as

(1)

The generalized Adair equation becomes

(2)
where {theta} denotes the fractional occupancy and <n> is the average binding number. As shown by others (2Go), the slope of the Hill plot for this scheme is

(3)
which simplifies to (2Go,4Go)

(4)
where <n2><n>2 is the variance of the binding number (2Go). I now wish to point out that the term N{theta}(1 – {theta}) is simply the variance of the binomial distribution case where the probability of a site to be bound by the ligand is {theta} and the probability for a site to be unoccupied is 1 – {theta}. This is the case for a noncooperative protein with the same fractional saturation at the same ligand concentration. Therefore the slope of the Hill plot for a generalized Adair equation can be written as

(5)
where {sigma} is the standard deviation of the binding number. This interpretation of the Hill plot holds for all values of {theta}. Since the Hill coefficient is defined as the slope of the Hill plot at , we can write for nH that

(6)

One can readily verify that this equation will reduce to N in the extremum case of infinite positive cooperativity; in this case, the only observable forms of the protein are the unligated and the fully ligated species. Thus, at a saturation of 0.5 the standard deviation of the occupancy will be equal toN/2, and we can write

(7)

Derivation of an empirical extremum principle for a negatively cooperating system of N identical binding sites
In contrast with the extremum case for positive cooperativity, in which no intermediately ligated species can be observed, I wish to define an empirical extremum scenario for negative cooperativity as a system with N identical binding sites where the first ligand binds with a finite dissociation constant, but subsequent bindings show an infinitely large dissociation constant. The experimental binding isotherm will reach an asymptote in which one ligand is bound. Subsequently, the observed (as opposed to the theoretical) maximal binding number is 1 and the standard deviation at the apparent fractional occupancy of 0.5, will be 0.5. By substituting these values into Eq. 6, we find for the slope of the Hill plot at {theta}app = 0.5 that

(8)

The minimal value of is 4: this can be readily derived by setting the derivative with respect to {theta} to 0 and identifying the value of {theta} at the extremum as 0.5. Hence the actual value of N{theta}(1 – {theta}) (Eq. 4) will not exceed N/4, and the slope of the Hill plot at the experimentally determined half-saturation point cannot be <1/N, and we can write

(9)

It is important to stress that this analysis is predicated on the fact that N (or more precisely, bmax, which equals N[P], where [P] is the protein concentration) is not explicitly known to the experimentalist; we assume that the experimentalist chooses the half-saturation point empirically according to standard curve-fitting procedures, at the concentration that yields one-half the extrapolated saturation value. If one uses the true value of N to define the half-saturation point (e.g., 2 in a tetramer) instead of the empirical saturation point, one is not able to derive this relationship, and the theoretical limit of nH becomes 0, as predicted by simple Koshland-Nemethy-Filmer models with infinite interaction between sites (13Go). A more general approach to this issue is to point out that a minimal variance would require the protein to be in one of two binding states; inclusion of additional binding states would increase the variance. Since the half-saturation point would be between the two binding states, and neglecting other binding states, one arrives again at a binomial two-state distribution. Since the minimal {Delta}n between any two states is 1, we arrive at the same 0.25 value for the minimal variance of a binding curve for a homogeneous population of sites at half-saturation.

Distinguishing negative cooperativity from a mixture of binding sites
One of the problematic aspects of negative cooperativity is that a Hill coefficient less than unity can be observed with a heterogeneous mixture of binding sites, without assuming interactions between the sites. If we can demonstrate that a mixture of sites can empirically reach arbitrarily small values of nH, this will add a tool for distinguishing the two cases.

One may derive an expression for the average fractional saturation of a mixture of N monomer classes with different dissociation constants Ki, as

(10)
where [PTi] denotes the total concentration of protein species i, Ki is the dissociation constant for the ith protein species with the ligand L, is the fractional saturation of the ith species (or {theta}i), and the summation runs over N protein species. This expression may be simplified by the introduction of the statistical weight {alpha}i, such that

(11)
and therefore

(12)

We can write for the slope of the Hill plot, from Eq. 3,

(13)
and differentiating Eq. 12 with respect to ln[L], we obtain

(14)

Note that in this case N is the number of species, and is not linked to properties of the individual protein molecules. Therefore, the slope cannot have a lower limit that is related to the number of subunits. An extremum situation for a mixture, analogous to that derived above for an N-mer of equivalent binding sites, occurs when one of the species (species 1 for the sake of discussion) has a much smaller dissociation constant than the rest such that there exists a range of [L] in the experimental isotherm where K1 << [L] << Ki holds, for all i != 1. This implies that species 1 is saturated under conditions where binding of the ligand to other species is negligible. As before, the experimentalist observes only the K1 isotherm, and higher terms in the summation are negligible at the half-saturation point of the K1 species. In this case, Eq. 14 reduces to

(15)
since all {theta}i {cong} 0 for i > 1. As the experimentalist has chosen the halfway point of the K1 isotherm as the halfway point of the whole titration, one may also approximate [L] {approx} K1, and the expression for the slope of the Hill plot at the experimental half-saturation point is then further simplified as

(16)

This treatment yields an arbitrarily small value of nH for a mixture, since there are no a priori minimal value constraints on {alpha}i that depend on molecular properties of the binding species.

Examples of nH values for known proteins and enzymes with negative cooperativity
To verify the relevance of these derivations, one may survey recent literature and compare nH values of characterized enzymes (Table 1). Of the 11 enzymes listed in Table 1 with negative cooperativity, and which have a known subunit composition, none have an nH value of <1/N. In addition, it should be noted that higher oligomerization states are able to reach lower nH values, as predicted from the derivation described above.


View this table:
[in this window]
[in a new window]
 
TABLE 1  N and nH values for documented protein-ligand and enzyme-substrate interactions showing negative cooperativity

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The extremum principle derived in this presentation is an empirically useful relationship that complements the well-known mathematical limit to the Hill coefficient at positive cooperativity. The experimentally derived values for known cases of negative cooperativity in which quaternary structures have been determined are in striking agreement to the central prediction, which limits nH to values >1/N in cases of true negative cooperativity. In Fig. 1, we observe a series of parallel vertical arrays of data points, each defining a line that represents a given N. In the dataset of Fig. 1, I have represented N = 2, 4, and 6, as described in Table 1. One can observe that there are six dimers in the data set, three tetramers, and two hexamers. Despite the paucity of higher oligomerization states, it is evident that as N increases, the data points reach lower values of nH, and in addition reach values closer to the theoretical limit of 1/N. In particular, the experimental value for V-ATPase (0.16 ± 0.03) is within experimental error of the limit (0.16666...). These results are in agreement with the predictions of the 1/N rule. This relationship can therefore be used as a rule of thumb for gauging the degree of coupling between sites. A particularly instructive case, which has been omitted from Table 1, is that of dopamine ß-mono-oxygenase. This is a tetrameric enzyme, which shows negative cooperativity with respect to ascorbate. The measured nH was found to be 0.15–0.30, and showed batch-dependent variability (14Go). Subsequently, it was shown that the protein exists as a mixture of dimers and tetramers, which have different kinetic parameters (15Go). Although both forms of the enzyme show negative cooperativity, the difference in kinetic parameters between the two forms creates a mixture of catalytic sites which leads, in some preparations, to a lower than expected value of nH. This demonstrates a practical application of the 1/N rule—that the observation of nH values <0.25 for a tetramer indeed correlates with the existence of a mixture of binding sites.



View larger version (10K):
[in this window]
[in a new window]
 
FIGURE 1  Graphic representation of nH versus 1/N for the interactions documented in Table 1. The diagonal line indicates the 1/N limiting value for nH.

 

    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
I thank Drs. Amnon Horovitz and Charles S. Weaver for helpful discussions and for critical reading of the manuscript. I also thank Prof. A. Ben-Shaul for his elegant lectures on statistical thermodynamics.

H.A. is supported by grant No. 496/03 from the Israel Science Foundation.

Submitted on January 26, 2005; accepted for publication April 4, 2005.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Hill, A. V. 1910. The possible effects of the aggregation of the molecules of haemoglobin on its oxygen dissociation curve. J. Physiol. (London). 40:4–7.

2. Wyman, J. J. 1964. Linked functions and reciprocal effects in hemoglobin: a second look. Adv. Protein Chem. 19:223–286.[Medline]

3. Weber, G., and S. J. Anderson. 1965. Multiplicity of binding. Range and practical test of Adair's equation. Biochemistry. 10:1942–1947.

4. Dahlquist, F. W. 1978. The meaning of Scatchard and Hill plots. Methods Enzymol. 48:270–299.[Medline]

5. Gill, S. J., H. T. Gaud, J. Wyman, and B. G. Barisas. 1978. Analysis of ligand binding curves in terms of species fractions. Biophys. Chem. 8:53–59.[CrossRef][Medline]

6. Goldbeter, A. 1976. Kinetic cooperativity in the concerted model for allosteric enzymes. Biophys. Chem. 4:159–169.[CrossRef][Medline]

7. Acerenza, L., and E. Mizraji. 1997. Cooperativity: a unified view. Biochim. Biophys. Acta. 1339:155–166.[CrossRef][Medline]

8. Bardsley, W. G., and J. Wyman. 1978. Concerning the thermodynamic definition and graphical manifestations of positive and negative cooperativity. J. Theor. Biol. 72:373–376.[CrossRef][Medline]

9. Conway, A., and D. E. Koshland, Jr. 1968. Negative cooperativity in enzyme action. The binding of diphosphopyridine nucleotide to glyceraldehyde 3-phosphate dehydrogenase. Biochemistry. 7:4011–4023.[CrossRef][Medline]

10. Levitzki, A., and D. E. Koshland, Jr. 1969. Negative cooperativity in regulatory enzymes. Proc. Natl. Acad. Sci. USA. 62:1121–1128.[Abstract/Free Full Text]

11. Levitzki, A., and D. E. Koshland, Jr. 1976. The role of negative cooperativity and half-of-the-sites reactivity in enzyme regulation. Curr. Topics Cell Regul. 10:1–40.[Medline]

12. Fersht, A. 1999. Structure and Mechanism in Protein Science. W.H. Freeman and Company, NY.

13. Levitzki, A. 1978. Quantitative Aspects of Allosteric Mechanisms. Springer-Verlag, Berlin, Germany.

14. Stewart, L. C., and J. P. Klinman. 1991. Cooperativity in the dopamine ß-mono-oxygenase reaction. Evidence for ascorbate regulation of enzyme activity. J. Biol. Chem. 266:11537–11543.[Abstract/Free Full Text]

15. Stewart, L., and J. P. Klinman. 1999. Kinetic parameters for dimeric and tetrameric forms of bovine dopamine ß-mono-oxygenase and their relationship to non-Michaelis-Menten behavior. FEBS Lett. 454:229–232.[CrossRef][Medline]

16. Lipps, G., and G. Krauss. 1999. Adenylosuccinate synthase from Saccharomyces cerevisiae: homologous overexpression, purification and characterization of the recombinant protein. Biochem. J. 341:537–543.[CrossRef][Medline]

17. Njalsson, R., S. Norgren, A. Larsson, C. S. Huang, M. E. Anderson, and J. L. Luo. 2001. Cooperative binding of {gamma}-glutamyl substrate to human glutathione synthetase. Biochem. Biophys. Res. Commun. 289:80–84.[CrossRef][Medline]

18. Murata, T., Y. Kakinuma, and I. Yamato. 2001. ATP-dependent affinity change of Na+-binding sites of V-ATPase. J. Biol. Chem. 276:48337–48440.[Abstract/Free Full Text]

19. Garnier, A., A. Berredjem, and B. Botton. 1997. Purification and characterization of the NAD-dependent glutamate dehydrogenase in the Ectomycorrhizal fungus laccaria bicolor (Maire) Orton. Fungal Genet. Biol. 22:168–176.[CrossRef][Medline]

20. Hinderlich, S., R. Stasche, R. Zeitler, and W. Reutter. 1997. A bifunctional enzyme catalyzes the first two steps in n-acetylneuraminic acid biosynthesis of rat liver. Purification and characterization of UDP-n-acetylglucosamine 2-epimerase/n-acetylmannosamine kinase. J. Biol. Chem. 272:24313–24318.[Abstract/Free Full Text]

21. Henis, Y. I., and A. Levitzki. 1980. Mechanism of negative cooperativity in glyceraldehyde-3-phosphate dehydrogenase deduced from ligand competition experiments. Proc. Natl. Acad. Sci. USA. 77:5055–5059.[Abstract/Free Full Text]

22. Caccuri, A. M., G. Antonini, P. Ascenzi, M. Nicotra, M. Nuccetelli, A. P. Mazzetti, G. Federici, M. Lo Bello, and G. Ricci. 1999. Temperature adaptation of glutathione S-transferase P1–1. A case for homotropic regulation of substrate binding. J. Biol. Chem. 274: 19276–19280.

23. Suzuki, Y., E. Moriyoshi, D. Tsuchiya, and H. Jingami. 2004. Negative cooperativity of glutamate binding in the dimeric metabotropic glutamate receptor subtype 1. J. Biol. Chem. 279:35526–35534.[Abstract/Free Full Text]

24. Kolodziej, A. F., T. Tan, and D. E. Koshland, Jr. 1996. Producing positive, negative, and no cooperativity by mutations at a single residue located at the subunit interface in the aspartate receptor of Salmonella typhimurium. Biochemistry. 35:14782–14792.[CrossRef][Medline]

25. Usova, E. V., and S. Eriksson. 1997. The effects of high salt concentrations on the regulation of the substrate specificity of human recombinant deoxycytidine kinase. Eur. J. Biochem. 248:762–766.[Medline]

26. Kanao, T., T. Fukui, H. Atomi, and T. Imanaka. 2002. Kinetic and biochemical analyses on the reaction mechanism of a bacterial ATP-citrate lyase. Eur. J. Biochem. 269:3409–3416.[Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.060194v1
89/1/76    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Abeliovich, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Abeliovich, H.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2005 by the Biophysical Society.