Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

The effect of net charge on the solubility, activity, and stability of ribonuclease Sa

The net charge and isoelectric pH (pI) of a protein depend on the content of ionizable groups and their pK values. Ribonuclease Sa (RNase Sa) is an acidic protein with a pI = 3.5 that contains no Lys residues. By replacing Asp and Glu residues on the surface of RNase Sa with Lys residues, we have created a 3K variant (D1K, D17K, E41K) with a pI = 6.4 and a 5K variant (3K + D25K, E74K) with a pI = 10.2. We show that pI values estimated using pK values based on model compound data can be in error by >1 pH unit, and suggest how the estimation can be improved. For RNase Sa and the 3K and 5K variants, the solubility, activity, and stability have been measured as a function of pH. We find that the pH of minimum solubility varies with the pI of the protein, but that the pH of maximum activity and the pH of maximum stability do not.     

pK values of histidine residues in ribonuclease Sa: effect of salt and net charge

The primary goal of this study was to gain a better understanding of the effect of environment and ionic strength on the pK values of histidine residues in proteins. The salt-dependence of pK values for two histidine residues in ribonuclease Sa (RNase Sa) (pI=3.5) and a variant in which five acidic amino acids have been changed to lysine (5K) (pI=10.2) was measured and compared to pK values of model histidine-containing peptides. The pK of His53 is elevated by two pH units (pK=8.61) in RNase Sa and by nearly one pH unit (pK=7.39) in 5K at low salt relative to the pK of histidine in the model peptides (pK=6.6). The pK for His53 remains elevated in 1.5M NaCl (pK=7.89). The elevated pK for His53 is a result of screenable electrostatic interactions, particularly with Glu74, and a non-screenable hydrogen bond interaction with water. The pK of His85 in RNase Sa and 5K is slightly below the model pK at low salt and merges with this value at 1.5M NaCl. The pK of His85 reflects mainly effects of long-range Coulombic interactions that are screenable by salt. The tautomeric states of the neutral histidine residues are changed by charge reversal. The histidine pK values in RNase Sa are always higher than the pK values in the 5K variant. These results emphasize that the net charge of the protein influences the pK values of the histidine residues. Structure-based pK calculations capture the salt-dependence relatively well but are unable to predict absolute histidine pK values.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

pKa predictions with a coupled finite difference Poisson-Boltzmann and Debye-Hückel method

Modeling charge interactions is important for understanding many aspects of biological structure and function, and continuum methods such as Finite Difference Poisson-Boltzmann (FDPB) are commonly employed. Calculations of pH-dependence have identified separate populations; surface groups that can be modeled with a simple Debye-Hückel (DH) model, and buried groups, with stronger resultant interactions that are dependent on detailed conformation. This observation led to the development of a combined FDPB and DH method for pK(a) prediction (termed FD/DH). This study reports application of this method to ionizable groups, including engineered buried charges, in staphylococcal nuclease. The data had been made available to interested research groups before publication of mutant structures and/or pK(a) values. Overall, FD/DH calculations perform as intended with low ΔpK(a) values for surface groups (RMSD between predicted and experimental pK(a) values of 0.74), and much larger ΔpK(a) values for the engineered internal groups, with RMSD = 1.64, where mutant structures were known and RMSD = 1.80, where they were modeled. The weaker resultant interactions of the surface groups are determined mostly by charge-charge interactions. For the buried groups, R(2) for correlation between predicted and measured ΔpK(a) values is 0.74, despite the high RMSDs. Charge-charge interactions are much less important, with the R(2) value for buried group ΔpK(a) values increasing to 0.80 when the term describing charge desolvation alone is used. Engineered charge burial delivers a relatively uniform, nonspecific effect, in terms of pK(a) . How the protein environment adapts in atomic detail to deliver this resultant effect is still an open question.     

Protein stabilization by salt bridges: concepts, experimental approaches and clarification of some misunderstandings

Salt bridges in proteins are bonds between oppositely charged residues that are sufficiently close to each other to experience electrostatic attraction. They contribute to protein structure and to the specificity of interaction of proteins with other biomolecules, but in doing so they need not necessarily increase a protein's free energy of unfolding. The net electrostatic free energy of a salt bridge can be partitioned into three components: charge-charge interactions, interactions of charges with permanent dipoles, and desolvation of charges. Energetically favorable Coulombic charge-charge interaction is opposed by often unfavorable desolvation of interacting charges. As a consequence, salt bridges may destabilize the structure of the folded protein. There are two ways to estimate the free energy contribution of salt bridges by experiment: the pK(a) approach and the mutation approach. In the pK(a) approach, the contribution of charges to the free energy of unfolding of a protein is obtained from the change of pK(a) of ionizable groups caused by altered electrostatic interactions upon folding of the protein. The pK(a) approach provides the relative free energy gained or lost when ionizable groups are being charged. In the mutation approach, the coupling free energy between interacting charges is obtained from a double mutant cycle. The coupling free energy is an indirect and approximate measure of the free energy of charge-charge interaction. Neither the pK(a) approach nor the mutation approach can provide the net free energy of a salt bridge. Currently, this is obtained only by computational methods which, however, are often prone to large uncertainties due to simplifying assumptions and insufficient structural information on which calculations are based. This state of affairs makes the precise thermodynamic quantification of salt bridge energies very difficult. This review is focused on concepts and on the assessment of experimental methods and does not cover the vast literature.     

Changing the net charge from negative to positive makes ribonuclease Sa cytotoxic

Ribonuclease Sa (pI = 3.5) from Streptomyces aureofaciens and its 3K (D1K, D17K, E41K) (pI = 6.4) and 5K (3K + D25K, E74K) (pI = 10.2) mutants were tested for cytotoxicity. The 5K mutant was cytotoxic to normal and v-ras-transformed NIH3T3 mouse fibroblasts, but RNase Sa and 3K were not. The structure, stability, and activity of the three proteins are comparable, but the net charge at pH 7 increases from -7 for RNase Sa to -1 for 3K and to +3 for 5K. These results suggest that a net positive charge is a key determinant of ribonuclease cytotoxicity. The cytotoxic 5K mutant preferentially attacks v-ras-NIH3T3 fibroblasts, suggesting that mammalian cells expressing the ras-oncogene are potential targets for ribonuclease-based drugs.     

Empirical relationships between protein structure and carboxyl pKa values in proteins

Relationships between protein structure and ionization of carboxyl groups were investigated in 24 proteins of known structure and for which 115 aspartate and 97 glutamate pK(a) values are known. Mean pK(a) values for aspartates and glutamates are < or = 3.4 (+/-1.0) and 4.1 (+/-0.8), respectively. For aspartates, mean pK(a) values are 3.9 (+/-1.0) and 3.1 (+/-0.9) in acidic (pI < 5) and basic (pI > 8) proteins, respectively, while mean pK(a) values for glutamates are approximately 4.2 for acidic and basic proteins. Burial of carboxyl groups leads to dispersion in pK(a) values: pK(a) values for solvent-exposed groups show narrow distributions while values for buried groups range from < 2 to 6.7. Calculated electrostatic potentials at the carboxyl groups show modest correlations with experimental pK(a) values and these correlations are not improved by including simple surface-area-based terms to account for the effects of desolvation. Mean aspartate pK(a) values decrease with increasing numbers of hydrogen bonds but this is not observed at glutamates. Only 10 pK(a) values are > 5.5 and most are found in active sites or ligand-binding sites. These carboxyl groups are buried and usually accept no more than one hydrogen bond. Aspartates and glutamates at the N-termini of helices have mean pK(a) values of 2.8 (+/-0.5) and 3.4 (+/-0.6), respectively, about 0.6 units less than the overall mean values.     

A buried lysine that titrates with a normal pKa: role of conformational flexibility at the protein-water interface as a determinant of pKa values

Previously we reported that Lys, Asp, and Glu residues at positions 66 and 92 in staphylococcal nuclease (SNase) titrate with pK(a) values shifted by up to 5 pK(a) units in the direction that promotes the neutral state. In contrast, the internal Lys-38 in SNase titrates with a normal pK(a). The crystal structure of the L38K variant shows that the side chain of Lys-38 is buried. The ionizable moiety is approximately 7 A from solvent and ion paired with Glu-122. This suggests that the pK(a) value of Lys-38 is normal because the energetic penalty for dehydration is offset by a favorable Coulomb interaction. However, the pK(a) of Lys-38 was also normal when Glu-122 was replaced with Gln or with Ala. Continuum electrostatics calculations were unable to reproduce the pK(a) of Lys-38 unless the protein was treated with an artificially high dielectric constant, consistent with structural reorganization being responsible for the normal pK(a) value of Lys-38. This reorganization must be local because circular dichroism and NMR spectroscopy indicate that the L38K protein is native-like under all conditions studied. In molecular dynamics simulations, the ion pair between Lys-38 and Glu-122 is unstable. The simulations show that a minor rearrangement of a loop is sufficient to allow penetration of water to the amino moiety of Lys-38. This illustrates both the important roles of local flexibility and water penetration as determinants of pK(a) values of ionizable groups buried near the protein-water interface, and the challenges faced by structure-based pK(a) calculations in reproducing these effects.     

Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsin.

The effects of solvation and charge-charge interactions on the pKa of ionizable groups in bacteriorhodopsin have been studied using a macroscopic dielectric model with atom-level detail. The calculations are based on the atomic model for bacteriorhodopsin recently proposed by Henderson et al. Even if the structural data are not resolved at the atomic level, such calculations can indicate the quality of the model, outline some general aspects of electrostatic interactions in membrane proteins, and predict some features. The effects of structural uncertainties on the calculations have been investigated by conformational sampling. The results are in reasonable agreement with experimental measurements of several unusually large pKa shifts (e.g. the experimental findings that Asp96 and Asp115 are protonated in the ground state over a wide pH range). In general, we find that the large unfavorable desolvation energies of forming charges in the protein interior must be compensated by strong favorable charge-charge interactions, with the result that the titrations of many ionizable groups are strongly coupled to each other. We find several instances of complex titration behavior due to strong electrostatic interactions between titrating sites, and suggest that such behavior may be common in proton transfer systems. We also propose that they can help to resolve structural ambiguities in the currently available density map. In particular, we find better agreement between theory and experiment when a structural ambiguity in the position of the Arg82 side-chain is resolved in favor of a position near the Schiff base.     

Multi-modality of pI distribution in whole proteome

Multi-modality of pI distribution is a common feature in different whole proteomes. Some researchers considered it relate to the proteins with different subcellular locations, indicating the result of natural selection. We explored the pI distribution of predicted proteomes (including animals, plants, bacterium, archaeans) and random proteome [random protein sequences constructed according to the special amino acid composition and molecular weight (MW) distribution of human predicted proteome]. Our results suggest that the multi-modality is the result of discrete pK(R) values for different amino acids. Amino acid composition and MW distribution of a proteome also contributes to the specific pI distribution. Although protein subcellular location was related to pI value, our analyses revealed that comparing with the random proteome, neither the multi-modality phenomenon nor the distribution bias of pI values is caused by subcellular location. It seems that the multi-modality distribution is just a mathematical fun. The blank region near the neutral pI was caused by the absence of amino acids with neutral pK(R), and suggests that the selection of amino acids with ionizable side chain might be restricted by the requirement for a special pH environment during the origin of life. From this point of view, the special distribution was the result of natural selection.     

Rearrangement of charge-charge interactions in variant ubiquitins as detected by double-mutant cycles and NMR

Previous studies of ubiquitin disclosed numerous charge-charge interactions on the protein's surface. To investigate how neighboring residues influence the strength of these interactions, double-mutant cycles are combined with pK(a) determinations by 2D NMR. More specifically, the environment around the Asp21-Lys29 ion pair has been altered through mutations at position 25, which is an asparagine in mammalian ubiquitin and a positively-charged residue in many other ubiquitin-like proteins. The pK(a) value of Asp21 decreases by 0.4 to 0.7 pH unit when Asn25 is substituted with a positively charged residue, suggesting a new and favorable ion pair interaction between positions 21 and 25. However, analysis of double mutants reveals that the favorable interaction between Asp21 and Lys29 is weakened when position 25 is a positively charged residue. Interestingly, while the pK(a) value of His25 in the N25H variant agrees with model compound values, additional mutants reveal that this agreement is fortuitous, resulting from a balance of favorable and unfavorable interactions; similar results were observed previously for Glu34 in ubiquitin and His8 in staphylococcal nuclease. Ionizable groups may thus have pK(a) values similar to model compound values and yet still be involved in significant interactions with other protein groups. One surprising result of introducing positively charged residues at position 25 is a new interaction between Lys29 and Glu18, an interaction not present in wild-type ubiquitin. This unanticipated result illustrates a key advantage of using NMR to determine pK(a) values for many residues simultaneously in the variant proteins. Overall, the strength of an interaction between two residues at the surface of ubiquitin is sensitive to the identity of neighboring residues. The results also demonstrate that relatively conservative and common point mutations such as substitutions of polar with charged residues and vice versa can have effects on interactions beyond the site of mutation per se.     

Charge-charge interactions in the denatured state influence the folding kinetics of ribonuclease Sa

Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

pK values of the ionizable groups of proteins

We have used potentiometric titrations to measure the pK values of the ionizable groups of proteins in alanine pentapeptides with appropriately blocked termini. These pentapeptides provide an improved model for the pK values of the ionizable groups in proteins. Our pK values determined in 0.1 M KCl at 25 degrees C are: 3.67+/-0.03 (alpha-carboxyl), 3.67+/-0.04 (Asp), 4.25+/-0.05 (Glu), 6.54+/-0.04 (His), 8.00+/-0.03 (alpha-amino), 8.55+/-0.03 (Cys), 9.84+/-0.11 (Tyr), and 10.40+/-0.08 (Lys). The pK values of some groups differ from the Nozaki and Tanford (N & T) pK values often used in the literature: Asp (3.67 this work vs. 4.0 N & T); His (6.54 this work vs. 6.3 N & T); alpha-amino (8.00 this work vs. 7.5 N & T); Cys (8.55 this work vs. 9.5 N & T); and Tyr (9.84 this work vs. 9.6 N & T). Our pK values will be useful to those who study pK perturbations in folded and unfolded proteins, and to those who use theory to gain a better understanding of the factors that determine the pK values of the ionizable groups of proteins.     

Thermal versus guanidine-induced unfolding of ubiquitin. An analysis in terms of the contributions from charge-charge interactions to protein stability.

We have characterized the guanidine-induced unfolding of both yeast and bovine ubiquitin at 25 degrees C and in the acidic pH range on the basis of fluorescence and circular dichroism measurements. Unfolding Gibbs energy changes calculated by linear extrapolation from high guanidine unfolding data are found to depend very weakly on pH. A simple explanation for this result involves the two following assumptions: (1) charged atoms of ionizable groups are exposed to the solvent in native ubiquitin (as supported by accessible surface area calculations), and Gibbs energy contributions associated with charge desolvation upon folding (a source of pK shifts) are small; (2) charge-charge interactions (another source of pK shifts upon folding) are screened out in concentrated guanidinium chloride solutions. We have also characterized the thermal unfolding of both proteins using differential scanning calorimetry. Unfolding Gibbs energy changes calculated from the calorimetric data do depend strongly on pH, a result that we attribute to the pH dependence of charge-charge interactions (not eliminated in the absence of guanidine). In fact, we find good agreement between the difference between the two series of experimental unfolding Gibbs energy changes (determined from high guanidine unfolding data by linear extrapolation and from thermal denaturation data in the absence of guanidine) and the theoretical estimates of the contribution from charge-charge interactions to the Gibbs energy change for ubiquitin unfolding obtained by using the solvent-accessibility-corrected Tanford-Kirkwood model, together with the Bashford-Karplus (reduced-set-of-sites) approximation. This contribution is found to be stabilizing at neutral pH, because most charged groups on the native protein interact mainly with groups of the opposite charge, a fact that, together with the absence of large charge-desolvation contributions, may explain the high stability of ubiquitin at neutral pH. In general, our analysis suggests the possibility of enhancing protein thermal stability by adequately redesigning the distribution of solvent-exposed, charged residues on the native protein surface.     

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.     

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.     

MCCE analysis of the pKas of introduced buried acids and bases in staphylococcal nuclease

The pK(a)s of 96 acids and bases introduced into buried sites in the staphylococcal nuclease protein (SNase) were calculated using the multiconformation continuum electrostatics (MCCE) program and the results compared with experimental values. The pK(a)s are obtained by Monte Carlo sampling of coupled side chain protonation and position as a function of pH. The dependence of the results on the protein dielectric constant (ε(prot)) in the continuum electrostatics analysis and on the Lennard-Jones non-electrostatics parameters was evaluated. The pK(a)s of the introduced residues have a clear dependence on ε(prot,) whereas native ionizable residues do not. The native residues have electrostatic interactions with other residues in the protein favoring ionization, which are larger than the desolvation penalty favoring the neutral state. Increasing ε(prot) scales both terms, which for these residues leads to small changes in pK(a). The introduced residues have a larger desolvation penalty and negligible interactions with residues in the protein. For these residues, changing ε(prot) has a large influence on the calculated pK(a). An ε(prot) of 8-10 and a Lennard-Jones scaling of 0.25 is best here. The X-ray crystal structures of the mutated proteins are found to provide somewhat better results than calculations carried out on mutations made in silico. Initial relaxation of the in silico mutations by Gromacs and extensive side chain rotamer sampling within MCCE can significantly improve the match with experiment.