HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Park, C. Articles by Marqusee, S. Search for Related Content PubMed PubMed Citation Articles by Park, C. Articles by Marqusee, S. Social Bookmarking                   What's this?

FOR THE RECORD

Analysis of the stability of multimeric proteins by effective G and effective m-values

Department of Molecular and Cell Biology and QB3 Institute, University of California, Berkeley, Berkeley, California 94720-3206, USA

Reprint requests to: Susan Marqusee, 215A Hildebrand Hall, Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720-3206, USA; e-mail: marqusee{at}uclink.berkeley.edu; fax: (510) 643-9290.

(RECEIVED April 14, 2004; FINAL REVISION May 27, 2004; ACCEPTED May 28, 2004)

	Abstract

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
Analyzing the stability of a multimeric protein is challenging because of the intrinsic difficulty in handling the mathematical model for the folded multimer-unfolded monomer equilibrium. To circumvent this problem, we introduce the concept of effective stability, G_eff (= –RTlnK_eff), where K_eff is the equilibrium constant expressed in monomer units. Analysis of the denaturant effect on G_eff gives new insight into the stability of multimeric proteins. When a multimeric protein is mostly folded, the dependence of effective stability on denaturant concentration (effective m-value) is simply the m-value of its monomeric unit. However, when the protein is mostly unfolded, its stability depends on denaturant concentration with the m-value of its multimeric form. We also find that the effective m-value at the C_m is a good approximation of the apparent m-value determined by fitting the equilibrium unfolding data from multimeric proteins with a two-state monomer model. Moreover, when the m-value of a monomeric unit is estimated from its size, the effective stability of a multimeric protein can be determined simply from C_m and this estimated m-value. These simple and intuitive approaches will allow a facile analysis of the stability of multimeric proteins. These analyses are also applicable for high-throughput analysis of protein stability on a proteomic scale.

Keywords: protein stability; m-value; multimeric protein; oligomeric protein

Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.04811004.

	Introduction

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
Determining protein stability is the first step in investigating the thermodynamic properties of proteins. The only requirement is that the protein establishes a reversible equilibrium between the folded and unfolded states. Most proteins whose stabilities have been studied so far are monomeric. Despite the prevalence of multimeric proteins, stability determinations of multimeric proteins are quite rare, mostly because of the inherent difficulty in the treatment and analysis of their equilibrium behavior. Enhancing our knowledge of the stability of multimeric proteins requires the development of new and facile ways to analyze their equilibria.

The equilibrium between an n-mer protein and its unfolded monomer without any intermediate state¹ is shown in equation 1:

(1)

where the equilibrium constant of the reaction is

The fraction of unfolded protein, f_D, is defined as in equation 2.

(2)

where P_t is the total protein concentration in monomer units. K_unf and G° can be expressed as functions of f_D:

(3)

(4)

The challenge in studying multimeric proteins results from their reaction order. Even for multimeric proteins, G° can be directly calculated from an experimentally determined f_D by equation 4 (Silinski et al. 2001). However, calculation of f_D from a given G° is nontrivial, because equation 4 is n-th order. Solving the equation is feasible for monomeric and dimeric proteins. However, even for trimeric proteins, deriving a general expression of f_D as a function of K_unf is challenging (Backmann et al. 1998; Silinski et al. 2001; Güthe et al. 2004):

For n > 3, the equation needs to be solved numerically (Johnson et al. 1995; Boudker et al. 1997; Panse et al. 2000). This problem precludes the development of a general description for the stability of multimeric proteins. Moreover, since K_unf is not unitless for multimeric proteins, it is much less intuitive and less useful for comparative purposes than is K_unf for monomeric proteins.

	The effective stability

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
Is there a more intuitive approach similar to that used for monomeric proteins to understand the thermodynamic stability of multimeric proteins? To address this question, we have developed an approach using an effective G (G_eff), instead of G° for the standard condition:

(5)

where f_D is defined as in equation 2. K_eff is simply the ratio of unfolded protein concentration to folded protein concentration in monomer units. G_eff and K_eff are pseudo-thermodynamic parameters where the actual order of chemical equation is ignored. Unlike G°, G_eff and K_eff are dependent on the total protein concentration of the system. However, these pseudo-thermodynamic parameters are more intuitive in that they relate directly to the amount of folded protein in the system. For example, G_eff is always zero when half of protein is unfolded, that is, f_D = 0.5, while G° of a multimeric protein has a nonzero value when f_D = 0.5.

Next, we derive a general equation for G_eff(H₂O). As mentioned earlier, analytical derivation of G_eff(H₂O) from G°(H₂O) in an n-mer system is mathematically impossible. If one assumes that in water f_D will be very small (a reasonable assumption for a stable protein under aqueous condition), then G°(H₂O) can be directly related to G_eff(H₂O). From equations 4 and 5,

Therefore,

(6)

Here, G°/n is the G° per monomer unit (Boudker et al. 1997), which is concentration-independent and an intrinsic property of the protein. The second term, however, is not protein-specific; it depends only on the total protein concentration and the number of monomer units in the multimer. This second term is identical for any n-meric protein under identical total protein concentration. Therefore, G_eff can be dissected into two terms: a protein-specific term, and a concentration-dependent term.

The separate terms in equation 6 allow a facile comparison of protein stability across proteins with different oligomeric states, which is not feasible with G°. Table 1 showsG°/n extracted from experimental data available for several multimeric proteins. While G° values vary significantly, G°/n values are all between 5 kcal/mole and 12 kcal/mole, which is also the usual range for monomeric protein stability. The G_eff values at 100 µM monomer concentration were also calculated for the same proteins (Table 1). When two proteins in different oligomeric states have the same G°/n values, the second term of the equation 6 explains the difference in G_eff, which is directly related to the unfavorable oligomerization at a nonstandard concentration. For example, G°/n of Trp aporepressor and 4-oxalocrotonate tautomerase (4-OT) are quite similar: 11.7 kcal/mole and 11.3 kcal/mole, respectively (Table 1). However, because Trp aporepressor is a dimer and 4-OT is a hexamer, folding of 4-OT involves a bigger loss of entropy than does Trp aporepressor. When the two proteins are at the same concentration, the effective stabilities of 4-OT and Trp aporepressor differ mostly in the second term of equation 6 (2.1 kcal/mole at 100 µM in monomer units).

View this table: [in this window] [in a new window]   Table 1. Comparison of G°, G°/n and Geff of various multimeric proteins

which shows that the mutational effect on G°/n reflects directly the change in G_eff at any protein concentration. This relationship also supports the utility of G_eff and G°/n in stability analysis of multimeric proteins.

	The effect of denaturant on Geff

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
The simplest method for determining protein stability is by monitoring chemical-induced denaturation profiles, combined with a two-state approximation and linear extrapolation (Greene and Pace 1974).

(7)

where D is the denaturant concentration and m is the dependence of the stability on denaturant. Unfolding equilibrium constants, K_unf, can be measured directly in the transition zone, where detectable amounts of folded and unfolded proteins are in equilibrium. Then, G°(H₂O) can be determined by extrapolating G° values (=RTlnK_unf) from varying concentrations of denaturant.

To show how denaturant affects G_eff of multimeric proteins, we define m_eff as:

(8)

By using equation 5,

(9)

f_D/D can be derived by using equations 4 and 7. The two equations are combined to yield:

(10)

By differentiating equation 10 with respect to D,

(11)

From equations 9 and 11,

(12)

This simple equation clearly shows the denaturant dependence of multimeric proteins. For a monomeric protein (n =1), m_eff is equal to m and constant with respect to denaturant. For multimeric proteins (n > 1), m_eff is not a constant. When the protein is mostly folded (f_D ??0), m_eff is close to m/n. When the protein is mostly unfolded (f_D ??1), m_eff is close to m.

Interestingly, m/n is approximately the expected m-value for a monomeric unit of an n-mer. Based on a study of reported m-values, Myers et al. demonstrated that m-values show a strong correlation with the size or the number of amino acids in a protein (Myers et al. 1995). It follows that the expected m-value for a hypothetical folded monomeric unit of an n-mer would be 1/n of the value of the n-mer protein (m/n). Therefore, when the protein is mostly folded, m_eff is approximately the expected m-value of the monomeric unit.

This denaturant-dependent behavior of m_eff is shown clearly in Figure 1. In this simulation, G_eff of each multimeric protein was determined numerically using equations 4 and 5. All the proteins were given the same expected m-value for a monomeric unit of 2 kcal/(mole•M) (m/n =2). For ease of comparison, all proteins were assigned to have the same C_m value (the denaturant concentration where f_D =0.5) of 3 M. As shown in Figure 1, under native condition (D < C_m), the denaturant dependence of G_eff shows a slope of 2, regardless of the oligomeric state. As predicted by equation 12, this slope corresponds to the expected m-value of the monomeric unit. Under denaturing conditions (D > C_m); however, the slope of the denaturant dependence for G_eff corresponds to the m-value of each oligomeric protein.

View larger version (21K): [in this window] [in a new window]   Figure 1. The effect of denaturant on effective G of multimeric proteins. n represents the oligomeric state of each protein. The m/n values and Cm of all the proteins were set to be 2 kcal/(mole•M) and 3 M, respectively, for comparison. G°(H2O) was determined with equation 14 by using these parameters. The fraction of unfolded protein, fD, was determined from G° by solving equation 4 numerically. Geff was calculated with this fD by equation 5. Total protein concentration was set to be 2 µM in monomer units. However, the result is indifferent to total protein concentration, because Cm was fixed.

which suggests that f_D is close to 1 at C_m^SUPREX. Therefore, m_eff determined by SUPREX is m_eff at f_D ??1, which is close to m/n as shown in equation 12.

	Analysis of denaturation data from multimeric proteins

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
The simplest treatment of chemically induced denaturation data from multimeric proteins would be to use the simple two-state monomer model. What m-values would result from this curve fitting? Even though this approach clearly uses an incorrect model, the curve fitting would yield the correct C_m, because this parameter is model-free. Moreover, only a limited range of G_eff can be determined by experiment—those in the transition zone where D is near C_m. Therefore, experimental G_eff values are unlikely to show the biphasic patterns shown in Figure 1. Rather, experimentally determined plots of G_eff versus D would show a more or less linear tendency unless n is very large. The m-value from the simplified curve-fitting with a monomer model would be close to the m_eff at C_m. From equation 12, m_eff at Cm can be calculated by setting f_D =0.5:

(13)

Table 2 shows that indeed the m-values obtained by this incorrect, yet simple, monomer model correlate extremely well with the expected m_eff at the midpoint determined by the above equation. The m-values in Table 2 were determined by fitting the f_D values used to generate Figure 1 with a two-state monomer model.

	Determination of Geff from Cm

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
Can G_eff(H₂O) be determined simply by knowing the mid-point? As mentioned earlier, with the statistical data from the study by Myers et al. (1995), one can estimate m-values simply from the size of proteins. In monomeric proteins, G_unf can be reasonably approximated by multiplying the estimated m-value and C_m. Is this approach valid for multimeric proteins? The simulated data in Figure 1 suggest that G_eff in H₂O might be approximated by C_m and m/n, the m-value for monomeric unit, regardless of the oligomeric state of a protein.

To verify this observation, we derive an equation for calculating G_eff from the C_m and m-value. In multimeric proteins, G° is not zero at the C_m. Rather, according to equation 4, G° at the C_m is expressed as (Backmann et al. 1998; Ragone 2000):

By using equation 7,

(14)

By using equations 6 and 14, G_eff can be expressed in a simple equation:

(15)

R(n) is a residual function of n. This function is responsible for the slight differences in G_eff(H₂O) values between different oligomers in Figure 1. For monomeric proteins, R(n) is zero. In all cases, R(n) is smaller than RTln2 (?0.41 kcal/mole at 25°C), which is quite negligible in most practical applications, and as shown in Figure 1. Because m/n is simply the m-value of a monomeric unit of the n-mer, G_eff(H₂O) can be deduced from two parameters: the m-value of a monomeric unit calculated from the statistical data by Myers et al. and the experimentally determined C_m. Surprisingly, knowing the number of monomers in the multimer is not necessary to determine G_eff(H₂O), because the two parameters, the m-value of the monomeric unit (m/n) and the C_m, are model-free.

	Conclusions

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
We show here that the use of effective G and effective m-value provides valuable information about the stabilities of multimeric proteins. G_eff is a useful parameter that gives an intuitive appreciation of the stability of multimeric proteins; one can calculate the fraction of unfolded protein from G_eff of a multimeric protein as easily as doing so from the G° of a monomeric protein. To understand the effect of denaturant on the stability of multimeric protein, we define m_eff as the derivative of G_eff with respect to denaturant concentration. Analysis of m_eff shows that denaturant affects G_eff in a biphasic way. When the m-value of an n-mer protein is m, m_eff of this protein is m/n under folding conditions and m under denaturing conditions. Thus, although direct analysis of G° is not feasible for multimeric proteins, G_eff and m_eff provide a facile alternative way to compare the stability of these proteins. Additionally,G_eff can be determined without prior knowledge of the oligomeric state of a protein. This unique property of G_eff could facilitate high-throughput analysis of protein stability on a proteomic scale.

	Footnotes

 
¹ An alternative unfolding mechanism is also observed in multimeric protein unfolding in which dissociation into monomers and unfolding are not coupled:

The unfolding transition in this mechanism is nothing but an unfolding transition of a monomeric protein, which can be analyzed with conventional methods. Therefore, here we only focus on the mechanism shown in equation 1.

	Acknowledgments

 
This work was supported by an NIH grant. We also thank D.E. Wildes and E.A. Shank for helpful discussions and comments on the manuscript. We also thank reviewers for their productive comments.

The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.

	References

TOP Abstract Introduction The effective stability The effect of denaturant... Analysis of denaturation data... Determination of {Delta}Geff... Conclusions References

 
Backmann, J., Schafer, G., Wyns, L., and Bonisch, H. 1998. Thermodynamics and kinetics of unfolding of the thermostable trimeric adenylate kinase from the archaeon Sulfolobus acidocaldarius. J. Mol. Biol. 284: 817–833.[CrossRef][Medline]

Boudker, O., Todd, M.J., and Freire, E. 1997. The structural stability of the cochaperonin GroES. J. Mol. Biol. 272: 770–779.[CrossRef][Medline]

Bowie, J.U. and Sauer, R.T. 1989. Equilibrium dissociation and unfolding of the arc repressor dimer. Biochemistry 28: 7139–7143.[CrossRef][Medline]

Ghaemmaghami, S., Fitzgerald, M.C., and Oas, T.G. 2000. A quantitative, high-throughput screen for protein stability. Proc. Natl. Acad. Sci. 97: 8296–8301.[Abstract/Free Full Text]

Gittelman, M.S. and Matthews, C.R. 1990. Folding and stability of trp aporepressor from Escherichia coli. Biochemistry 29: 7011–7020.[CrossRef][Medline]

Grant, S.K., Deckman, I.C., Culp, J.S., Minnich, M.D., Brooks, I.S., Hensley, P., Debouck, C., and Meek, T.D. 1992. Use of protein unfolding studies to determine the conformational and dimeric stabilities of HIV-1 and SIV proteases. Biochemistry 31: 9491–9501.[CrossRef][Medline]

Greene, R.F. and Pace, C.N. 1974. Urea and guanidine-hydrochloride denaturation of ribonuclease, lysozyme, -chymotrypsin, and -lactoglobulin. J. Biol. Chem. 249: 5388–5393.[Abstract/Free Full Text]

Güthe, S., Kapinos, L., Möglish, A., Meier, S., Grzesiek, S., and Kiefhaber, T. 2004. Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin. J. Mol. Biol. 337: 905–915.[CrossRef][Medline]

Johnson, C.R., Morin, P.E., Arrowsmith, C.H., and Freire, E. 1995. Thermodynamic analysis of the structural stability of the tetrameric oligomerization domain of p53 tumor suppressor. Biochemistry 34: 5309–5316.[CrossRef][Medline]

Myers, J.K., Pace, C.N., and Scholtz, J.M. 1995. Denaturant m values and heat capacity changes: Relation to changes in accessible surface areas of protein unfolding. Protein Sci. 4: 2138–2148.[Abstract]

Panse, V.G., Swaminathan, C.P., Aloor, J.J., Surolia, A., and Varadarajan, R. 2000. Unfolding thermodynamics of the tetrameric chaperone, SecB. Biochemistry 39: 2362–2369.[CrossRef][Medline]

Powell, K.D., Wales, T.E., and Fitzgerald, M.C. 2002. Thermodynamic stability measurements on multimeric proteins using a new H/D exchange- and matrix-assisted laser desorption/ionization (MALDI) mass spectrometry-based method. Protein Sci. 11: 841–851.[Abstract/Free Full Text]

Ragone, R. 2000. How the protein concentration affects unfolding curves of oligomers. Biopolymers 53: 221–225.[CrossRef][Medline]

Silinski, P., Allingham, M.J., and Fitzgerald, M.C. 2001. Guanidine-induced equilibrium unfolding of a homo-hexameric enzyme 4-oxalocrotonate tau-tomerase (4-OT). Biochemistry 40: 4493–4502.[CrossRef][Medline]

CiteULike    Connotea    Del.icio.us    Digg    Reddit    Technorati    What's this?

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Park, C. Articles by Marqusee, S. Search for Related Content PubMed PubMed Citation Articles by Park, C. Articles by Marqusee, S. Social Bookmarking                   What's this?