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.     

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.     

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.     

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.     

Role of the charge-charge interactions in defining stability and halophilicity of the CspB proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.     

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.     

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.     

Effects of charge-to-alanine substitutions on the stability of ribosomal protein L30e from Thermococcus celer

The ability to rationally engineer a protein with altered stability depends upon the detailed understanding of the role of noncovalent interactions in defining thermodynamic properties of proteins. In this paper, we used T. celer L30e as a model to address the question of the role of charge-charge interactions in defining the stability of this protein. A total of 26 single-site charge-to-alanine variants of this protein were generated, and the stability of these proteins was determined using thermal- and denaturant-induced unfolding. It was found that, although L30e is isolated from a thermophilic organism and is highly thermostable, some of the substitutions lead to a further increase in the transition temperature. Analysis of the effects of high ionic strength on the stabilities of L30e variants shows that the long-range charge-charge interactions are as important as the short-range (salt bridge) interactions. The changes in stabilities of the T. celer L30e protein variants were compared with the changes in the energy of charge-charge interactions calculated using different computational models. It was found that there is a good qualitative agreement between experimental and calculated data: for 70-80% (19-21 of 26, confidence p < 0.003) of the variants, computational models predict correctly the sign of the stability changes. In particular, computational models identify correctly those charged amino acid residue substitutions of which led to enhancement in thermostability. Thus, optimization of the charge-charge interactions might be a useful approach for the rational increase in protein stability.     

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.     

Charged residues at protein interaction interfaces: unexpected conservation and orchestrated divergence

Protein-protein interactions play an essential role in the functioning of cell. The importance of charged residues and their diverse role in protein-protein interactions have been well studied using experimental and computational methods. Often, charged residues located in protein interaction interfaces are conserved across the families of homologous proteins and protein complexes. However, on a large scale, it has been recently shown that charged residues are significantly less conserved than other residue types in protein interaction interfaces. The goal of this work is to understand the role of charged residues in the protein interaction interfaces through their conservation patterns. Here, we propose a simple approach where the structural conservation of the charged residue pairs is analyzed among the pairs of homologous binary complexes. Specifically, we determine a large set of homologous interactions using an interaction interface similarity measure and catalog the basic types of conservation patterns among the charged residue pairs. We find an unexpected conservation pattern, which we call the correlated reappearance, occurring among the pairs of homologous interfaces more frequently than the fully conserved pairs of charged residues. Furthermore, the analysis of the conservation patterns across different superkingdoms as well as structural classes of proteins has revealed that the correlated reappearance of charged residues is by far the most prevalent conservation pattern, often occurring more frequently than the unconserved charged residues. We discuss a possible role that the new conservation pattern may play in the long-range electrostatic steering effect.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

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.     

Ionic Interactions Promote Transmembrane Helix-Helix Association Depending on Sequence Context

Folding and oligomerization of integral membrane proteins frequently depend on specific interactions of transmembrane helices. Interacting amino acids of helix-helix interfaces may form complex motifs and exert different types of molecular forces. Here, a set of strongly self-interacting transmembrane domains (TMDs), as isolated from a combinatorial library, was found to contain basic and acidic residues, in combination with polar nonionizable amino acids and C-terminal GxxxG motifs. Mutational analyses of selected sequences and reconstruction of high-affinity interfaces confirmed the cooperation of these residues in homotypic interactions. Probing heterotypic interaction indicated the presence of interhelical charge-charge interactions. Furthermore, simple motifs of an ionizable residue and GxxxG are significantly overrepresented in natural TMDs, and a specific combination of these motifs exhibits high-affinity heterotypic interaction. We conclude that intramembrane charge-charge interactions depend on sequence context. Moreover, they appear important for homotypic and heterotypic interactions of numerous natural TMDs.     

Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions

Computational design of surface charge-charge interactions has been demonstrated to be an effective way to increase both the thermostability and the stability of proteins. To test the robustness of this approach for proteins with predominantly beta-sheet secondary structure, the chicken isoform of the Fyn SH3 domain was used as a model system. Computational analysis of the optimal distribution of surface charges showed that the increase in favorable energy per substitution begins to level off at five substitutions; hence, the designed Fyn sequence contained four charge reversals at existing charged positions and one introduction of a new charge. Three additional variants were also constructed to explore stepwise contributions of these substitutions to Fyn stability. The thermodynamic stabilities of the variants were experimentally characterized using differential scanning calorimetry and far-UV circular dichroism spectroscopy and are in very good agreement with theoretical predictions from the model. The designed sequence was found to have increased the melting temperature, DeltaT (m) = 12.3 +/- 0.2 degrees C, and stability, DeltaDeltaG(25 degrees C) = 7.1 +/- 2.2 kJ/mol, relative to the wild-type protein. The experimental data suggest that a significant increase in stability can be achieved through a very small number of amino acid substitutions. Consistent with a number of recent studies, the presented results clearly argue for a seminal role of surface charge-charge interactions in determining protein stability and suggest that the optimization of surface interactions can be an attractive strategy to complement algorithms optimizing interactions in the protein core to further enhance protein stability.     

Fluctuations between stabilizing and destabilizing electrostatic contributions of ion pairs in conformers of the c-Myc-Max leucine zipper

In solution proteins often exhibit backbone and side-chain flexibility. Yet electrostatic interactions in proteins are sensitive to motions. Hence, here we study the contribution of ion pairs toward protein stability in a range of conformers which sample the conformational space in solution. Specifically, we focus on the electrostatic contributions of ion pairs to the stability of each of the conformers in the NMR ensemble of the c-Myc-Max leucine zipper and to their average energy minimized structure. We compute the electrostatic contributions of inter- and intra-helical ion pairs and of an ion pair network. We find that the electrostatic contributions vary considerably among the 40 NMR conformers. Each ion pair, and the network, fluctuates between being stabilizing and being destabilizing. This fluctation reflects the variability in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other and with respect to other charged groups in the rest of the protein. Ion pair interactions in the c-Myc-Max leucine zipper in solution depend on the protein conformer which is analyzed. Hence, the overall stabilizing (or destabilizing) contribution of an ion pair is conformer population-dependent. This study indicates that free energy calculations performed using the continuum electrostatics methodology are sensitive to protein conformational details.     

Statistical studies of flexible nonhomogeneous polypeptide chains

Unfolded proteins attract increasing attention nowadays because of the accumulation of experimental evidence that they play an important role in different biological processes. Therefore, studies of various statistical properties of flexible protein-like polypeptide chains are becoming increasingly important as well. This paper presents distributions (histograms) of distances between atoms of titratable residues for flexible polypeptide chains with various residue compositions and with the hard-spheres potential taken into consideration. The factors influencing the parameters of the obtained histograms have been identified and analyzed. It was found that the sensitivity of the distributions with respect to the internal structure of intermediate residues increases with the number of residues between the considered charged residues. It was shown that branching at C(beta) atoms of the side chains of the intermediate residues is among the most considerable factors influencing the shape of the distance distribution and the average distance between atoms in flexible chains. Despite the model simplicity, the results of the calculations can be applied for systems with other types of interactions presented, and this was demonstrated for the charge-charge interactions. In particular, it was shown that those interactions have a significant effect on distances between the unlike charges, while such an effect for the like charges is much less pronounced. The comparison of predictions made on the basis of the presented calculations to some experimental data is also given, and possible applications of the theoretical concept described in the paper are discussed.     

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).     

The effect of protein relaxation on charge-charge interactions and dielectric constants of proteins.

The effect of the reorganization of the protein polar groups on charge-charge interaction and the corresponding effective dielectric constant (epsilon(eff)) is examined by the semimicroscopic version of the Protein Dipole Langevin Dipoles (PDLD/S) method within the framework of the Linear Response Approximation (LRA). This is done by evaluating the interactions between ionized residues in the reaction center of Rhodobacter sphaeroides, while taking into account the protein reorganization energy. It is found that an explicit consideration of the protein relaxation leads to a significant increase in epsilon(eff) and that semimicroscopic models that do not take this relaxation into account force one to use a large value for the so-called "protein dielectric constant," epsilon(p), of the Poisson-Boltzmann model or for the corresponding epsilon(in) in the PDLD/S model. An additional increase in epsilon(eff) is expected from the reorganization of ionized residues and from changes in the degree of water penetration. This finding provides further support for the idea that epsilon(in) (or epsilon(p)) represents contributions that are not considered explicitly. The present study also provides a systematic illustration of the nature of epsilon(eff), supporting our previously reported view that charge-charge interactions correspond to a large value of this "dielectric constant," even in protein interiors. It is also pointed out that epsilon(eff) for the interaction between ionizable groups in proteins is very different from the effective dielectric constant, epsilon'(eff), that determines the free energy of ion pairs in proteins (epsilon'(eff) reflects the effect of preoriented protein dipoles). Finally, the problems associated with the search for a general epsilon(in) are discussed. It is clarified that the epsilon(in) that reproduces the effect of protein relaxation on charge-charge interaction is not equal to the epsilon(in) that reproduces the corresponding effect upon formation of individual charges. This reflects fundamental inconsistencies in attempts to cast microscopic concepts in a macroscopic model. Thus one should either use a large epsilon(in) for charge-charge interactions and a small epsilon(in) for charge-dipole interactions or consider the protein relaxation microscopically.