A new scale for side-chain contribution to protein stability based on the empirical stability analysis of mutant proteins

Kazufumi Takano and Katsuhide Yutani,1

1 Institute for Protein Research, Osaka University, Yamadaoka, Suita,Osaka 565-0871, Japan

E-mail: yutani{at}protein.osaka-u.ac.jp


    Abstract
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
The hydrophobicity scales for amino acid side chains based on the transfer Gibbs energy ({Delta}{Delta}Gtrans) of amino acids from non-aqueous phases to water have been widely used to estimate the contribution of buried side chains to the conformational stability of proteins. In this paper, we propose a new scale for the side-chain contribution to protein stability, which is derived from data on protein denaturation experiments using systematic and comprehensive mutant proteins. In the experiments, the contribution of some physical properties were quantitatively determined as parameters in a unique equation representing the stability change ({Delta}{Delta}{Delta}{Delta}G) of mutant proteins as a function of the structural changes due to the mutations. These parameters are able conveniently to provide a scale for the side-chain contribution to protein stability. This new scale also has the advantage over the previously reported hydrophobicity scales of residues with the contributions of hydrogen bonds or secondary structural propensity. It may find practical application in algorithms for the prediction of protein structures.

Keywords: hydrophobicity/mutant protein/protein denaturation/protein stability/stability scale/transfer Gibbs energy


    Introduction
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
The hydrophobic effect is believed to play a major role in organizing the self-assembly of globular proteins (Kauzmann, 1959Go; Dill, 1990Go). To quantify this effect, many hydrophobicity scales for the amino acid residues have been proposed (Nozaki and Tanford, 1971Go; Wolfenden et al., 1981Go; Fauchere and Pliska, 1983Go; Rose et al., 1985Go; Damodaran and Song, 1986Go; Kawashima et al., 1999Go). These scales are based on the transfer Gibbs energy ({Delta}Gtrans) of the amino acid residues from a non-aqueous phase to an aqueous phase. The octanol scale taken from Fauchere and Pliska (1983) is probably the most commonly used for estimating the contribution of buried non-polar side chains to the conformational stability of proteins.

A comparison of several scales, however, reveals that significant differences exist among them. Especially the values for polar or aromatic side chains are highly dependent on the non-aqueous phase (Karplus, 1997Go). There is then the question of which non-aqueous phase best characterizes the protein interior (Dill, 1990Go) and even some doubt as to whether they can be applied to protein denaturation (Vajda et al., 1995Go). This evidence suggests that we need a hydrophobicity scale of amino acid residues obtained from data on protein denaturation experiments. To obtain hydrophobicity parameters from protein experiments, we will introduce an alternative method based on the quantifiable changes in protein stability due to mutations (Sharp et al., 1991Go; Vajda et al., 1995Go).

Recently, it has been proposed that the stability change ({Delta}{Delta}G) of each mutant human lysozyme can be precisely represented in a unique equation by considering the conformational changes due to the mutations (Funahashi et al., 1999Go; Takano et al., 1999bGo). This {Delta}{Delta}G calculation provides each contribution of several physical properties, such as the hydrophobic effect and hydrogen bonds, to protein stability and it is compatible with other proteins (Funahashi et al., 2001).

In this paper, we propose a new scale for the side-chain contribution to protein stability, which is derived from the empirical {Delta}{Delta}G calculation on mutant proteins (Funahashi et al., 2001). This new scale presented here is more valuable than the previous hydrophobicity scales; it is based on more precise data from experimental studies on proteins and it is also able to show the contribution of residues with hydrogen bonds and residues in secondary structures.


    Calculation of {Delta}{Delta}G of mutant proteins
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
It has been proposed that the changes in stability ({Delta}{Delta}G, kJ/mol) of each mutant protein are represented by a unique equation, as expressed by Equation 1Go (Funahashi et al., 2001). In this equation, the contributions of the hydrophobic effect ({Delta}{Delta}GHP), side-chain conformational entropy ({Delta}{Delta}Gconf), hydrogen bonds ({Delta}{Delta}GHB), water molecules ({Delta}{Delta}GH2O), secondary structural propensity ({Delta}{Delta}Gpro: {Delta}{Delta}Gpro{alpha} and {Delta}{Delta}Gproß) and cavity volume ({Delta}{Delta}Gcav) to the stability are represented by each parameter in terms of the conformational change due to the mutation, as expressed by Equations 2–7GoGoGoGoGoGo. These parameters can be generally applied to globular proteins (Funahashi et al., 2001).

(1)

(2)

(3)

(4)

(5)

(6)

(7)
where {Delta}ASANP and {Delta}ASAP represent the differences in the ASA (accessible surface area, Å2) of the non-polar (C and S) and polar (N and O) atoms, respectively, of all residues upon denaturation; {Delta}Sconf is the difference in the side-chain conformational entropy upon denaturation defined by Doig and Sternberg (1995); rpp, rpw and rww are the length (Å) of the intramolecular, protein–water and wate–water hydrogen bonds, respectively, introduced/deleted due to mutation; NH2O is the number of water molecules introduced; P{alpha} and Pß are the {alpha}-helix and ß-sheet propensities, respectively, of the residue defined by Chou and Fasman, (1978) [revised by Koehl and Levitt (1999)]; and Vcav represents the cavity volume (Å3). The coefficients in Equations 2Go and 4–7GoGoGoGo have been estimated by a least-squares fit of each {Delta}{Delta}G (Equation 1Go) to the experimental {Delta}{Delta}G values using the stability–structure database of mutant proteins (Funahashi et al., 2001). The meaning of each parameter has been discussed in the previous papers by Funahashi et al. (1999, 2001).


    Side-chain contribution to protein stability ({Delta}{Delta}Gaa) based on {Delta}{Delta}G calculation
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
The parameters of various stabilization factors used in the {Delta}{Delta}G calculation described above are appropriate for a basic understanding of protein stability. We are able to extract the general contributions of each amino acid side chain for the conformational stability of globular proteins relative to Gly ({Delta}{Delta}Gaa) from the parameters of the hydrophobic effect and side chain entropy:

(8)
If an amino acid residue is buried in the interior of a protein, the {Delta}{Delta}Gaa value is calculated using the ASA value of the amino acid residue (Miller et al., 1987Go) and the entropic effect of the amino acid residue (Doig and Sternberg, 1995Go), as shown in Table IGo. As the {Delta}{Delta}Gaa values include the entropic effect, the {Delta}{Delta}Gaa values might not correspond to the {Delta}Gtrans values (hydrophobicity), so they represent the contributions of each amino acid side chain to the conformational stability of the globular proteins. The {Delta}{Delta}GHP values might correspond to the {Delta}Gtrans values. The {Delta}Gtrans values of a vapor (Wolfenden et al., 1981Go), cyclohexane (Radzicka and Wolfenden, 1988) and octanol (Fauchere and Pliska, 1983Go) to water are also listed in Table IGo. The {Delta}{Delta}Gtrans value of a vapor to water is different from those of the other scales, because there is little van der Waals interaction in the vapor phase.


View this table:
[in this window]
[in a new window]
 
Table I. Side-chain contribution to protein stability ({Delta}{Delta}Gaa)a and transfer Gibbs energy values ({Delta}Gtrans) (kJ/mol)
 
In the {Delta}{Delta}Gaa values, all residues (except for Gly) compensate the entropic loss of the side chain with the hydrophobic effect, resulting in a more positive contribution to the stability than the Gly residue. This means that even polar side chains contribute favorably to the protein stability through the hydrophobic effect from their non-polar atoms. This point is different from previous transfer values; usually polar residues would contribute negatively to the stability.


    Side-chain contribution of residues with hydrogen bonds to protein stability ({Delta}{Delta}GaaHB)
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
The side chains of buried polar residues usually form hydrogen bonds with the protein atoms or water molecules. The hydrogen bonds contribute to the protein stability (Myers and Pace, 1996Go; Yamagata et al., 1998Go; Takano et al., 1999aGo,bGo). We are also able to estimate the contribution of residues with hydrogen bonds ({Delta}{Delta}GaaHB) using the parameters of the hydrogen bonds:

(9)
If the length of a hydrogen bond is 3.0 Å, the contributions of one intramolecular and one protein-water hydrogen bonds are 7.4 and 3.0 kJ/mol, respectively. The {Delta}{Delta}GaaHB values of each amino acid residue with several hydrogen bonds are summarized in Table IIGo.


View this table:
[in this window]
[in a new window]
 
Table II. Side-chain contribution of residues with hydrogen bonds to protein stability ({Delta}{Delta}GaaHB)a (kJ/mol)
 
In Table IIGo, the polar residues with hydrogen bonds contribute to the stability of proteins, similarly to hydrophobic residues. This means that a protein structure is stabilized not only by buried hydrophobic residues but also by buried polar residues with hydrogen bonds. For example, for the Val and Thr residues, which are isosteric residues, the Thr residue with one intramolecular hydrogen bond has a contribution (14.3 kJ/mol) comparable to the Val residue (15.3 kJ/mol). Furthermore, when a Thr residue has two intramolecular hydrogen bonds, it contributes (21.7 kJ/mol) more than a Val residue. In fact, a Thr residue forms about two hydrogen bonds on average in protein structures (McDonald and Thornton, 1994Go). Recently, Pace (Pace, 2001) has pointed out in another way that burial of a polar residue contributes more to protein stability than burial of a non-polar residue.


    Side-chain contribution of residues on secondary structures to protein stability ({Delta}{Delta}Gaa{alpha} and {Delta}{Delta}Gaaß)
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
In secondary structures, the {alpha}-helix and ß-sheet, certain amino acids are found more frequently, but others are found less frequently (Chou and Fasman, 1978Go). This tendency is known as {alpha}-helix or ß-sheet propensity. These propensities affect the protein stability (Pace and Scholtz, 1998Go; Street and Mayo, 1999Go), suggesting that the stability scales of residues in the secondary structures are different to those in the other regions. We also calculated the contributions of residues in the {alpha}-helix ({Delta}{Delta}Gaa{alpha}) and in the ß-sheet ({Delta}{Delta}Gaaß) relative to Gly using the parameters of the secondary structural propensities:

(10)

(11)
where {Delta}P{alpha} and {Delta}Pß are the {alpha}-helix and ß-sheet propensities, respectively, of the residue relative to Gly defined by Chou and Fasman (Chou and Fasman, 1978Go) [revised by Koehl and Levitt (1999)]. The {Delta}{Delta}Gaa{alpha} and {Delta}{Delta}Gaaß values are summarized in Table IIIGo. If a residue on a secondary structure forms hydrogen bonds, the contribution of hydrogen bonds could definitely be added as described above.


View this table:
[in this window]
[in a new window]
 
Table III. Side-chain contribution of residues on secondary structures to protein stability ({Delta}{Delta}Gaa{alpha} and {Delta}{Delta}Gaaß)a (kJ/mol)
 
Conclusion

It is not easy to predict the tertiary structure of a protein from the amino acid sequence, although it is one of the most important subjects in life science research. One problem is that the computational approaches to predicting protein structures mostly use the hydrophobicity scales based on the transfer Gibbs energy as the contribution of the amino acid residues to protein stability. We now propose a new scale for the side chain contribution to stability based on an experimental study of proteins. This scale is different to the previous scales, especially for polar residues. We expect that this scale would be useful to theoretical researchers in improving structural predictions. Furthermore, the contributions of residues in various environments, such as residues with hydrogen bonds or residues in secondary structures, are also estimated. This scale will also contribute significantly to improving the protein stability.


    Notes
 
1 To whom correspondence should be addressed. E-mail: yutani{at}protein.osaka-u.ac.jp Back


    Acknowledgments
 
This work was supported in part by Fellowships from the Japan Society for the Promotion of Science for Young Scientists (K.T.) and by a Grant-in-Aid for Scientific Research on Priority Areas (C) `Genome Information Science' from the Ministry of Education, Science, Sports and Culture of Japan (K.Y.).


    References
 Top
 Abstract
 Introduction
 Calculation of {Delta}{Delta}G...
 Side-chain contribution to...
 Side-chain contribution of...
 Side-chain contribution of...
 References
 
Chou,P.Y. and Fasman,G.D. (1978) Adv. Enzymol., 47, 45–148.[Medline]

Damodaran,S. and Song,K.B. (1986) J. Biol. Chem., 261, 7220–7222.[Abstract/Free Full Text]

Dill,K.A. (1990) Biochemistry, 29, 7133–7155.[ISI][Medline]

Doig,A.J. and Sternberg,M.J.E. (1995) Protein Sci., 4, 2247–2251.[Abstract/Free Full Text]

Fauchere,J.-L. and Pliska,V. (1983) Eur. J. Med. Chem., 18, 369–375.[ISI]

Funahashi,J., Takano,K., Yamagata,Y. and Yutani,K. (1999) Protein Eng., 12, 841–850.[Abstract/Free Full Text]

Funahashi,J., Takano,K. and Yutani,K. (2001) Protein Eng., 14, 127–134.[Abstract/Free Full Text]

Karplus,P.A. (1997) Protein Sci., 6, 1302–1307.[Abstract/Free Full Text]

Kauzmann,W. (1959) Adv. Protein Chem., 14, 1–63.[ISI]

Kawashima,S., Ogata,H. and Kanehisa,M. (1999) Nucleic Acids Res., 27, 368–369.[Abstract/Free Full Text]

Koehl,P. and Levitt,M. (1999) Proc. Natl Acad. Sci. USA, 96, 12524–12529.[Abstract/Free Full Text]

McDonald,I.K. and Thornton,J.M. (1994) J. Mol. Biol., 238, 777–793.[ISI][Medline]

Miller,S., Janin,J., Lesk,A.M. and Chothia,C. (1987) J. Mol. Biol., 196, 641–656.[ISI][Medline]

Myers,J.K. and Pace,C.N. (1996) Biophys. J., 71, 2033–2039.[Abstract]

Nozaki,Y. and Tanford,C. (1971) J. Biol. Chem., 246, 2211–2217.[Abstract/Free Full Text]

Pace,C.N. (2001) Biochemistry, 40, 310–313.[ISI][Medline]

Pace,C.N. and Scholtz,J.M. (1998) Biophys. J., 75, 422–427.[Abstract/Free Full Text]

Radzika,A. and Wolfenden,R. (1988) Biochemistry, 27, 1644–1670.

Rose,G.D., Geselowitz,A.R., Lesser,G.J., Lee,R.H. and Zehfus,M.H. (1985) Science, 229, 834–838.[ISI][Medline]

Sharp,K.A., Nicholls,A., Friedman,R. and Honig,B. (1991) Biochemistry, 30, 9686–9697.[ISI][Medline]

Street,A.G. and Mayo,S.L. (1999) Proc. Natl Acad. Sci. USA, 96, 9074–9076.[Abstract/Free Full Text]

Takano,K., Yamagata,Y., Kubota,M., Funahashi,J., Fujii,S. and Yutani,K. (1999a) Biochemistry, 38, 6623–6629.[ISI][Medline]

Takano,K., Yamagata,Y., Funahashi,J., Hioki,Y., Kuramitsu,S. and Yutani,K. (1999b) Biochemistry, 38, 12698–12708.[ISI][Medline]

Vajda,S., Weng,Z. and DeLisi,C. (1995) Protein Eng., 11, 1081–1092.

Wolfenden,R., Andersson,L., Cullis,P.M. and Southgate,C.C.B. (1981) Biochemistry, 20, 849–855.[ISI][Medline]

Yamagata,Y., Kubota,M., Sumikawa,Y., Funahashi,J., Takano,K., Fujii,S. and Yutani,K. (1998) Biochemistry, 37, 9355–9362.[ISI][Medline]

Received February 9, 2001; revised May 25, 2001; accepted July 8, 2001.