| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306; and Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104
Correspondence: Address reprint requests to Huan-Xiang Zhou, Institute of Molecular Biophysics, Florida State University, Tallahassee, FL 32306. Tel.: 850-644-7052; Fax: 850-644-0098; E-mail: hxzhou{at}csit.fsu.edu.
| ABSTRACT |
|---|
|
|
|---|
G) and experimental results have a correlation coefficient of 0.98. Bc-Csp appears to use destabilization of the unfolded state by unfavorable charge-charge interactions as a mechanism for increasing stability. Accounting for the effects of ionic strength and temperature on the electrostatic free energies in both the folded and the unfolded states, explanations for two important experimental observations are presented. The disparate ionic strength dependences of
G for Bc-Csp and Bs-CspB were attributed to the difference in the total charges (-2e and -6e, respectively). A main contribution to the much higher unfolding entropy of Bs-CspB was found to come from the less favorable electrostatic interactions in the folded state. These results should provide insight for understanding the thermostability of other thermophilic proteins. | INTRODUCTION |
|---|
|
|
|---|
|
|
![]() | (1) |
= (8
Ie2/
skBT)1/2,
s is the dielectric constant of water, and Q is the total charge on the molecule. According to Eq. 1, when
R << 1 (where R is the size of the protein), regardless of the specific distribution of the charges, the electrostatic free energy decreases with ionic strength; the larger the magnitude of the total charge, the steeper the decrease. Equation 1 allows us to attribute the disparate ionic strength dependences of the unfolding free energies of Bc-Csp and Bs-CspB to the difference in their total charges, -2e and -6e.
Because of the strong temperature dependence of the dielectric constant of water, electrostatic interactions are expected to contribute to the temperature dependence of the unfolding free energy (Zhou, 2002d
). The favorable interactions of a charge with water molecules implies that the hydration water molecules will be organized and thus have a lower entropy than bulk water. Upon protein folding, the solute-solvent interactions will in general become less favorable (giving rise to the desolvation cost) and solvent organization will be less, resulting in an entropy gain (illustrated in Fig. 3 A). When two individually solvated charges (as in the unfolded state) are brought together in the folded state, the entropy changes will be very different depending on whether the charges have the same or opposite signs (illustrated in Fig. 3, B and C). Two like charges next to each other will reinforce each other's solvent-organizing effect, but two opposite charges next to each other will cancel each other's solvent-organizing effect. These results indicate that favorable electrostatic interactions in the folded state may contribute a loss of unfolding entropy whereas unfavorable electrostatic interactions may contribute a gain in unfolding entropy. We therefore suggest that the observed higher unfolding entropies of Bc-CspB and Bc-Csp mutants such as R3E/L66E are a result of the less favorable electrostatic interactions in the folded state (relative to Bc-Csp). In agreement with experiment, the calculated value of T
S is greater than the magnitude of 
G, resulting in a higher unfolding enthalpy for Bs-CspB.
|
| THEORETICAL METHODS |
|---|
|
|
|---|
![]() | (2) |
|

Gel0. To isolate the electrostatic effect of the mutation, structural changes in the folded state were kept to a minimum. That is, only the mutated side chain was allowed to optimize its conformation within the existing environment using the program InsightII (Molecular Simulations, San Diego, CA.). The AMBER force field (Weiner et al., 1984
The PB equation was solved by the UHBD program (Madura et al., 1995
). The electrostatic potential
was calculated first from a 100 x 100 x 100 grid with 1.5 Å spacing centered at the geometric center of the wild-type protein. This was followed by a 140 x 140 x 140 grid with 0.5 Å spacing at the same center. A final round of focusing at the CB atom of a mutated side chain was introduced on a 60 x 60 x 60 grid with 0.25 Å spacing. The electrostatic energies of the folded protein and the isolated mutation residue were calculated by
![]() | (3) |

Gel0 was obtained by taking the difference in Gel between the isolated mutation residue (modeling the unfolded state) and the folded protein and then the change in these differences by the mutation. AMBER charges and radii were used. The solvent dielectric constant was set to that of water at a particular temperature (Archer and Wang, 1990
p) was set to 4. Here we also investigated the effect of varying this parameter by carrying our calculations with
p = 2, 4, and 10.

Gel0 can be decomposed into a solvation term for the mutation residue, a solvation term for the rest of the folded protein, and a term representing the interactions between the mutation residue and the rest of the protein. Thus
![]() | (4) |
Gintf. By multiplying the electrostatic potential of the mutation residue alone with the charges of other residues,
Gintf can be further decomposed into contributions from interactions with individual residues.
Twenty-seven mutations of Bc-Csp and Bs-CspB were studied in three sets of solvent conditions: 1), I = 100 mM and T = 70°C; 2), I = 2.1 M and T = 70°C; and 3), I = 100 mM and T = 25°C. For calculating Gel,0u and Gelf, Asp, Glu, and the C-terminal were unprotonated, whereas the N-terminal, Lys, and Arg were protonated. For calculating Gintu, the total energy of charge-charge interactions in the unfolded state (to model the effects of ionic strength and temperature), the same set of protonation states were assumed. However, for modeling the effects of charge mutations on Gintu, sampling of protonation states were allowed (Zhou, 2002a
).
Ionic strength dependence
To model the effects of ionic strength, Gelf for the folded protein, Gel,0u for all the charged residues, and Gintu were calculated at I = 100, 200, 300, 400, and 500 mM.
Temperature dependence
The electrostatic contribution to the unfolding entropy was calculated as
![]() | (5) |
| RESULTS AND DISCUSSION |
|---|
|
|
|---|

Gel for 21 mutations on Bc-Csp and six mutations on Bs-CspB at I = 100 mM and T = 70°C are listed in Table 1. Those with
p = 4 are compared to experimental data of Schmid and co-workers (Perl and Schmid, 2001
p is decreased to 2 or increased to 10, with the root-mean deviation increasing from 0.46 kcal/mol to 0.53 and 0.48 kcal/mol, respectively. Overall, the correlation coefficient between the calculation results with
p = 4 and experimental data is 0.98. Unless otherwise noted, the results below are for
p = 4.
|
|

Gel rises sharply. The calculated results on 
Gel quantitatively reproduce the measurements of Schmid and co-workers (Perl et al., 2000
G between the R3L and R3A mutants (-0.9 vs. -1.9 kcal/mol) suggests a nonpolar contribution as well, but this cannot be accounted for within our electrostatic model.
It is of interest to note that as much as 0.5 kcal/mol of the 2.7 kcal/mol decrease in folding stability of the R3E mutant comes from stabilizing the unfolded state. In the unfolded state, charge-charge interactions are dominated by residues close along the sequence, since these will have higher probabilities of sampling short distances (Zhou, 2002a
). For unfolded Bc-Csp, R3 is in an unfavorable position because the charges immediately preceding and after it (the N-terminal and residues K5 and K7) are all positive. Mutation to E3 will change unfavorable interactions with these positive charges into favorable ones and stabilize the unfolded state.
Of the 2.2 kcal/mol of -
Gel0 for the R3E mutations, 0.4 kcal/mol comes from higher desolvation cost for E3 (than R3), -0.1 kcal/mol comes from the difference in solvating the rest of the protein, and 1.9 kcal/mol comes from less favorable interactions with the protein environment. In particular, interactions with E46, E21, and the C-terminal account for 1.6, 0.6, and 0.5 kcal/mol, respectively, of -
Gel0. These are compensated by interactions with the N-terminal and K5, which contribute 1.3 kcal/mol to 
Gel0. Weak interactions with numerous other residues make up the remaining portion of 
Gel0.
The above details of the various contributions to the 
Gel0 result for R3E also help understanding the effects of additional mutations on the R3E mutant. Specifically, the R3E/E21A and R3E/E46A mutants have smaller magnitudes of 
Gel0 (-1.5 and -0.8 kcal/mol, respectively) mainly because an unfavorable interaction (with either E21 or E46) is now absent. In contrast, in the R3E/L66E mutant, an additional unfavorable interaction is introduced, and consequently -
Gel0 increases to 4.5 kcal/mol. In the R3E/E46A/L66E triple mutant (which resembles Bs-CspB the most among all R3E mutants), an unfavorable interaction is removed but another is introduced, and the value of 
Gel0, -2.94 kcal/mol, is intermediate. These results generally are in accord with experimental observations. However, there is disagreement regarding the interaction between E3 and E21. The fact that the measured 
G does not change upon the additional mutation E21A on the R3E mutant suggests absence of interaction between these two residues, whereas our calculations found an interaction energy of 0.6 kcal/mol. The apparent discrepancy may indicate the limitation of our calculations in not allowing for the mutated and neighboring residues to sample different rotamers.
Our calculations show that the unfolded state is also significantly stabilized by the E21A and E46A mutations. In both cases, the immediately preceding and following charges are all negative (E19 and D25 for the former and E42, E43, and E50 for the latter). All together, residual charge-charge interactions for the R3E/E46A/L66E triple mutant account for 0.76 kcal/mol, or 20%, of the total loss of 3.7 kcal/mol in folding stability.
In agreement with experimental observations (Perl et al., 2000
; Perl and Schmid, 2001
), most of the other 13 mutations on Bc-Csp have marginal effects on the folding stability. The increase in stability by converse mutations on Bs-CspB have been discussed previously (Dong and Zhou, 2002
).
Electrostatic interactions can be weakened by an increase in ionic strength or an increase in the dielectric constant of water brought by a decrease in temperature. Comparisons of calculated and experimental results for the effects of the mutations at I = 2.1 M and T = 70°C and at I = 100 mM and T = 25°C are displayed in Fig. 5, B and C. Again, there is reasonable overall agreement between calculation and experiment. It is of interest to note that, even at an ionic strength of 2.1 M, there are still substantial electrostatic interactions in the folded state, which may be attributed to the compactness of the protein. This finding casts doubt on attributing residual effects of charge mutations at a high ionic strength (such as 2 M) solely to a nonelectrostatic origin (such as hydrophobic interactions). At room temperature, our calculations on the effects of the mutations again indicate significant electrostatic contributions to the folding stability of Bc-Csp, but with smaller magnitudes as a result of the increase in the dielectric constant of water.
Effect of ionic strength on overall stability
The role of the total charge in the effect of ionic strength can be understood from the general result given in Eq. 1 for low ionic strengths and the specific result for a spherical protein:
![]() | (6) |
RKl-&1by2;(
R)/Kl+&1by2;(
R), with Kl+&1by2;(x) denoting modified Bessel functions; ri is the radial distance of charge qi;
ij is the angle between the radial vectors of charges qi and qj; and Pl(x) are the Legendre polynomials. Both Eq. 1 and Eq. 6 show that a dominant contribution to the ionic strength dependence of Gel is a term proportional to Q2. The total charges on Bc-Csp and Bs-CspB are -2e and -6e, respectively, thus the electrostatic free energy of folded Bs-CspB is expected to decrease with ionic strength much more rapidly than that of Bc-Csp. Our calculation results conform to this expectation. From I = 100 to 500 mM, Gel decreased 1.04 kcal/mol for folded Bc-Csp but 2.53 kcal/mol for folded Bs-CspB at T = 70°C. The electrostatic free energy of the unfolded state also decreased with ionic strength. Specifically, the solvation energies of all the charged groups decreased by 2.07 kcal/mol for unfolded Bc-Csp and by 2.30 kcal/mol for unfolded Bs-CspB, and the residual charge-charge interaction energy decreased by 0.36 kcal/mol for unfolded Bc-Csp and by 0.26 kcal/mol for unfolded Bs-CspB. Taken together,
Gel decreased by 1.4 kcal/mol for Bc-Csp but was unchanged for Bs-CspB.
Experimentally, Perl and Schmid (2001)
observed a slight decrease in the unfolding free energy from I = 100 to 500 mM for Bc-Csp and a significant increase in
G for Bs-CspB. For both proteins,
G increases from I = 500 mM to 2.1 M. The increase in stability by high ionic strengths is likely due to the Hofmeister effect of salts (Perl and Schmid, 2001
; Baldwin, 1996
). Our implicit assumption is that all protein variants studied here are stabilized to the same extent by the Hofmeister effect. With this assumption, we would predict that the difference in stability between Bc-Csp and Bs-CspB will decrease by 1.4 kcal/mol from I = 100 to 500 mM. This decrease is in quantitative agreement with the measurement of Schmid and co-workers.
The main source of the disparate ionic strength dependences of the unfolding free energy for Bc-Csp and Bs-CspB is the difference in the total charges. For Bs-CspB, the larger magnitude of the total charge leads to a sharper decrease with ionic strength in the electrostatic free energy of the folded state. This sharper decrease is at nearly the same pace with which the electrostatic free energy of the unfolded state varies with ionic strength, leading to an almost constant
G. On the other hand, Bc-Csp has a small net charge and the ionic strength dependence in the folded state is much weaker than in the unfolded state; consequently
G is found to decrease with ionic strength. A recent continuum model study by Dominy et al. (2002)
yielded similar results and reached the same conclusion on the role of the total charge. Within the PB model, the total charge as a determinant of the ionic strength dependence of
G is inevitable.
Two additional pieces of evidence support the attribution of the disparate ionic strength dependences of the stabilities of Bc-Csp and Bs-CspB to the difference in the total charges. On the calculation side, the R3E/L66E mutant of Bc-Csp, with an increase of total charge from -2e to -5e is found to have the same ionic strength dependence as Bs-CspB, i.e.,
Gel remains constant from I = 100 to 500 mM. On the experimental side, the ionic strength dependences of Bs-CspB and Bc-Csp were found to bracket a majority of the mutants with intermediate total charges (Perl and Schmid, 2001
).
Contributions of electrostatic interactions to temperature dependence of
G
The contribution of electrostatic interactions to the unfolding entropy was calculated according to Eq. 5 from the temperature dependence of the electrostatic free energy of the folded protein and the solvation energies of the charged residues and the residual charge-charge interaction energy in the unfolded state. With
p = 4,
Sel was 8.7 cal/mol/K for Bc-Csp but 29.9 cal/mol/K for Bs-CspB at I = 100 mM and T = 70°C. The difference in
Sel is 21.2 cal/mol/K, which accounts for more than half of the measured difference of 40 cal/mol/K in the unfolding entropy (Perl and Schmid, 2001
). Similarly, the change in
Sel by the R3E/L66E mutation is 19.4 cal/mol/K, accounting for half of the measured change of 41 cal/mol/K in the unfolding entropy. The magnitude of
Sel increases at
p = 2 and decreases at
p = 10. However, the differences in
Sel of the R3E/L66E mutant and Bs-CspB with Bc-Csp are nearly constant at
20 cal/mol/K.
We are now in a position to decompose the -5.0 kcal/mol calculated change in unfolding free energy by the R3EL66E mutation into entropic and enthalpic components. The entropic component T
Sel is 343 x 19.4 cal/mol = 6.7 kcal/mol, and the enthalpic component is thus 
Hel = 
Gel + T
Sel = 1.7 kca/mol. Though the calculated 
Hel is much smaller than the measured change in unfolding enthalpy (which was 9.1 kcal/mol), there is qualitative agreement in that the unfolding enthalpy of the R3EL66E mutant (and by extension Bs-CspB) is higher than that of Bc-Csp. We now see that the higher unfolding enthalpy of the R3EL66E mutant is a consequence of the higher unfolding entropy, which in turn may be attributed to the less favorable electrostatic interactions in the folded protein.
Implications for other thermophilic proteins
The effects of charge mutations, ionic strength, and temperature on the folding stability report complementary information on the contributions of electrostatic interactions. The effects of charge mutations usually are dominated by the desolvation cost and interactions with nearby charged and polar groups. Both our study and the study of Dominy et al. (2002)
indicate that the total charge plays a dominant role in the dependence on ionic strength. Our study further suggests that overall more favorable environments for charged residues lead to a smaller unfolding entropy.
A key to our approach is the explicit treatment of the unfolded state, which often is relegated to the status of an implicit "reference state". Only after accounting for the temperature and ionic strength dependence of the electrostatic free energy of the unfolded state can one properly obtain information on the unfolding entropy and enthalpy and the ionic strength dependence of the unfolding free energy. Moreover, in several cases (R3, E21, and E46), the thermophilic Bc-Csp appears to use destabilization of the unfolded state as a mechanism for increasing thermostability.
It is of great interest to see whether this intriguing mechanism is used by other thermophilic proteins. In addition, the control of ionic strength dependence of protein stability by the total charge and the contribution of electrostatic interactions to the unfolding entropy suggested by the present study warrant further testing. Applications of the present approach to other thermophilic proteins, with the interrogation of experiments, should provide a more comprehensive understanding on the role of electrostatic interactions in thermostability.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported in part by National Institutes of Health grant GM58187.
Submitted on July 29, 2002; accepted for publication December 4, 2002.
| REFERENCES |
|---|
|
|
|---|
Baldwin, R. L. 1996. How Hofmeister ion interactions affect protein stability. Biophys. J. 71:20562063.
Boschitsch, A. H., M. O. Fenley, and H.-X. Zhou. 2002. Fast boundary element method for the linear Poisson-Boltzmann equation. J. Phys. Chem. B. 106:27412754.
Delbruck, H., U. Mueller, D. Perl, F. X. Schmid, and U. Heinemann. 2001. Crystal structures of mutant forms of the Bacillus caldolyticus cold shock protein differing in thermal stability. J. Mol. Biol. 313:359369.[Medline]
Dominy, B. N., D. Perl, F. X. Schmid, and C. L. Brooks 3rd. 2002. The effects of ionic strength on protein stability: the cold shock protein family. J. Mol. Biol. 319:541554.[Medline]
Dong, F., and H.-X. Zhou. 2002. Electrostatic contributions to T4 lysozyme stability: solvent-exposed charges versus semi-buried salt bridges. Biophys. J. 83:13411347.
Jaenicke, R., and G. Bohm. 1998. The stability of proteins in extreme environments. Curr. Opin. Struct. Biol. 8:738748.[Medline]
Madura, J. D., J. M. Briggs, R. C. Wade, M. E. Davis, B. A. Luty, A. Ilin, J. Antosiewicz, M. K. Gilson, B. Bagheri, L. R. Scott, and J. A. McCammon. 1995. Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian Dynamics program. Comput. Phys. Commun. 91:5795.
Mueller, U., D. Perl, F. X. Schmid, and U. Heinemann. 2000. Thermal stability and atomic-resolution crystal structure of the Bacillus caldolyticus cold shock protein. J. Mol. Biol. 297:975988.[Medline]
Pace, C. N. 2000. Single surface stabilizer. Nat. Struct. Biol. 7:345346.[Medline]
Pace, C. N., R. W. Alston, and K. L. Shaw. 2000. Charge-charge interactions influence the denatured state ensemble and contribute to protein stability. Protein Sci. 9:13951398.[Abstract]
Perl, D., U. Mueller, U. Heinemann, and F. X. Schmid. 2000. Two exposed amino acid residues confer thermostability on a cold shock protein. Nat. Struct. Biol. 7:380383.[Medline]
Perl, D., and F. X. Schmid. 2001. Electrostatic stabilization of a thermophilic cold shock protein. J. Mol. Biol. 213:343357.
Perutz, M. F. 1978. Electrostatic effects in proteins. Science. 201:11871191.
Perutz, M. F., and H. Raidt. 1975. Stereochemical basis of heat stability in bacterial ferredoxins and in haemoglobin A2. Nature. 255:256259.[Medline]
Petsko, G. A. 2001. Structural basis of thermostability in hyperthermophilic proteins, or "there's more than one way to skin a cat". Methods Enzymol. 334:469478.[Medline]
Sanchez-Ruiz, J. M., and G. I. Makhatadze. 2001. To charge or not to charge? Trends Biotechnol. 19:132135.[Medline]
Schindelin, H., M. A. Marahiel, and U. Heinemann. 1993. Universal nucleic acid-binding domain revealed by crystal structure of the B. subtilis major cold-shock protein. Nature. 364:164168.[Medline]
Szilagyi, A., and P. Zavodszky. 2000. Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: results of a comprehensive survey. Structure. 8:493504.[Medline]
Vijayakumar, M., and H.-X. Zhou. 2001. Salt bridges stabilize the folded structure of barnase. J. Phys. Chem. B. 105:73347340.
Vogt, G., and P. Argos. 1997. Protein thermal stability, hydrogen bonds, and ion pairs. Fold. Des. 2:S40S46.[Medline]
Weiner, S. J., P. A. Kollman, D. A. Case, U. C. Singh, C. Ghio, G. Alagona, S. Prefeta, and P. Weiner. 1984. A new force field for molecular mechanical simulation of nucleic acids and proteins. J. Am. Chem. Soc. 106:765784.
Zhou, H.-X. 2002a. A Gaussian-chain model for treating residual charge-charge interactions in the unfolded state of proteins. Proc. Natl. Acad. Sci. USA. 99:35693574.
Zhou, H.-X. 2002b. Residual electrostatic effects in the unfolded state of the N-terminal domain of L9 can be attributed to nonspecific nonlocal charge-charge interactions. Biochemistry. 41:65336538.[Medline]
Zhou, H.-X. 2002c. Residual charge interactions in unfolded staphylococcal nuclease can be explained by the Gaussian-chain model. Biophys. J. 83:29812986.
Zhou, H.-X. 2002d. Toward the physical basis of thermophilic proteins: linking of enriched polar interactions and reduced heat capacity of unfolding. Biophys. J. 83:31263133.
This article has been cited by other articles:
![]() |
X. Huang and H.-X. Zhou Similarity and Difference in the Unfolding of Thermophilic and Mesophilic Cold Shock Proteins Studied by Molecular Dynamics Simulations Biophys. J., October 1, 2006; 91(7): 2451 - 2463. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Wunderlich and F. X. Schmid The correlation between protein stability and dipole moment: a critical test Protein Eng. Des. Sel., August 1, 2006; 19(8): 355 - 358. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Perez-Jimenez, R. Godoy-Ruiz, B. Ibarra-Molero, and J. M. Sanchez-Ruiz The Efficiency of Different Salts to Screen Charge Interactions in Proteins: A Hofmeister Effect? Biophys. J., April 1, 2004; 86(4): 2414 - 2429. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Torrez, M. Schultehenrich, and D. R. Livesay Conferring Thermostability to Mesophilic Proteins through Optimized Electrostatic Surfaces Biophys. J., November 1, 2003; 85(5): 2845 - 2853. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Dong, M. Vijayakumar, and H.-X. Zhou Comparison of Calculation and Experiment Implicates Significant Electrostatic Contributions to the Binding Stability of Barnase and Barstar Biophys. J., July 1, 2003; 85(1): 49 - 60. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |