Energetics of charge-charge interactions between residues adjacent in sequence.

Statistical analysis of the residue separation between a pair of ionizable side chains within 4 Å of each other was performed on a set of 1560 non-homologous PDB structures. We found that the frequency of pairs of like charges (i.e., pairs consisting of acidic residues Asp and Glu or pairs consisting of basic residues Arg and Lys) is two orders of magnitude lower than the pairs of oppositely charged residues (salt-bridges). We also found that for pairs of like charges the distribution is skewed dramatically towards short residue separation (<3). On the basis of these observations, we hypothesize that at short residue separation the repulsion between charges does not contribute much to the protein stability and the effects are largely dominated by the long range charge-charge interactions with other ionizable groups in the protein molecule. To test this hypothesis, we incorporated various pairs of charged residues at position 63 and 64 of ubiquitin and compared the stabilities of these variants. We also performed calculations of the expected changes in the charge-charge interactions. A very good correlation between experimental changes in the stability of ubiquitin variants, and changes in the energy of charge-charge interactions provides support for the hypothesis that a pair of ionizable residues next to each other in sequence modulates protein stability via long range charge-charge interactions with the rest of the protein.

Protein stability and surface electrostatics: a charged relationship.

Engineering proteins to withstand a broad range of conditions continues to be a coveted objective, holding the potential to advance biomedicine, industry, and our understanding of disease. One way of achieving this goal lies in elucidating the underlying interactions that define protein stability. It has been shown that the hydrophobic effect, hydrogen bonding, and packing interactions between residues in the protein interior are dominant factors that define protein stability. The role of surface residues in protein stability has received much less attention. It has been believed that surface residues are not important for protein stability particularly because their interactions with the solvent should be similar in the native and unfolded states. In the case of surface charged residues, it was sometimes argued that solvent exposure meant that the high dielectric of the solvent will further decrease the strength of the charge-charge interactions. In this paper, we challenge the notion that the surface charged residues are not important for protein stability. We computationally redesigned sequences of five different proteins to optimize the surface charge-charge interactions. All redesigned proteins exhibited a significant increase in stability relative to their parent proteins, as experimentally determined by circular dichroism spectroscopy and differential scanning calorimetry. These results suggest that surface charge-charge interactions are important for protein stability and that rational optimization of charge-charge interactions on the protein surface can be a viable strategy for enhancing protein stability.

Both helical propensity and side-chain hydrophobicity at a partially exposed site in alpha-helix contribute to the thermodynamic stability of ubiquitin.

Improving helical propensity of residues was proposed as one of the approaches to increase protein stability. Here the contribution of the helix propensity and hydrophobicity of residues at partially buried positions of alpha-helix to the stability of a model protein-ubiquitin- is explored. Thermodynamic stabilities of 13 ubiquitin variants with substitutions at a partially buried helical residue were measured by differential scanning calorimetry. It was found that the dynamic range of stabilities for different amino acid residues at this partially buried position is 3 times larger than that expected based on helical propensity alone. Correlation analysis shows that both helical propensity and hydrophobicity are important in defining the relative stabilities of the studied ubiquitin variants. These results provide experimental evidence that partially buried positions are potentially useful sites for engineering proteins with enhanced thermostability.

Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions.

Using computational and sequence analysis of bacterial cold shock proteins, we designed a protein (CspB-TB) that has the core residues of mesophilic protein from Bacillus subtilis(CspB-Bs) and altered distribution of surface charged residues. This designed protein was characterized by circular dichroism spectroscopy, and found to have secondary and tertiary structure similar to that of CspB-Bs. The activity of the CspB-TB protein as measured by the affinity to a single-stranded DNA (ssDNA) template at 25 degrees C is somewhat higher than that of CspB-Bs. Furthermore, the decrease in the apparent binding constant to ssDNA upon increase in temperature is much more pronounced for CspB-Bs than for CspB-TB. Temperature-induced unfolding (as monitored by differential scanning calorimetry and circular dichroism spectroscopy) and urea-induced unfolding experiments were used to compare the stabilities of CspB-Bs and CspB-TB. It was found that CspB-TB is approximately 20 degrees C more thermostable than CspB-Bs. The thermostabilization of CspB-TB relative to CspB-Bs is achieved by decrease in the enthalpy and entropy of unfolding without affecting their temperature dependencies, i.e. these proteins have similar heat capacity changes upon unfolding. These changes in the thermodynamic parameters result in the global stability function, i.e. Gibbs energy, deltaG(T), that is shifted to higher temperatures with only small changes in the maximum stability. Such a mechanism of thermostabilization, although predicted from the basic thermodynamic considerations, has never been identified experimentally.

Contribution of surface salt bridges to protein stability: guidelines for protein engineering.

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior.

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.

Removal of surface charge-charge interactions from ubiquitin leaves the protein folded and very stable.

The contribution of solvent-exposed charged residues to protein stability was evaluated using ubiquitin as a model protein. We combined site-directed mutagenesis and specific chemical modifications to first replace all Arg residues with Lys, followed by carbomylation of Lys-amino groups. Under the conditions in which all carboxylic groups are protonated (at pH 2), the chemically modified protein is folded and very stable (DeltaG = 18 kJ/mol). These results indicate that surface charge-charge interactions are not an essential fundamental force for protein folding and stability.