Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.     

Coupled molecular dynamics and continuum electrostatic method to compute the ionization pKa’s of proteins as a function of pH. Test on a large set of proteins

A computational method, to predict the pKa values of the ionizable residues Asp, Glu, His, Tyr, and Lys of proteins, is presented here. Calculation of the electrostatic free-energy of the proteins is based on an efficient version of a continuum dielectric electrostatic model. The conformational flexibility of the protein is taken into account by carrying out molecular dynamics simulations of 10 ns in implicit water. The accuracy of the proposed method of calculation of pKa values is estimated from a test set of experimental pKa data for 297 ionizable residues from 34 proteins. The pKa-prediction test shows that, on average, 57, 86, and 95% of all predictions have an error lower than 0.5, 1.0, and 1.5 pKa units, respectively. This work contributes to our general understanding of the importance of protein flexibility for an accurate computation of pKa, providing critical insight about the significance of the multiple neutral states of acid and histidine residues for pKa-prediction, and may spur significant progress in our effort to develop a fast and accurate electrostatic-based method for pKa-predictions of proteins as a function of pH.     

Electrostatic forces in two lysozymes: calculations and measurements of histidine pKa values.

In order to examine the electrostatic forces in globular proteins, pKa values and their ionic strength dependence of His residues of hen egg white lysozyme (HEWL) and human lysozyme (HUML) were measured, and they were compared with those calculated numerically. pKa values of His residues in HEWL, HUML, and short oligopeptides were determined from chemical shift changes of His side chains by 1H-nmr measurements. The associated changes in pKa values in HEWL and HUML were calculated by solving the Poisson-Boltzmann equations numerically for macroscopic dielectric models. The calculated pKa changes and their ionic strength dependence agreed fairly well with the observed ones. The contribution from each residue of each alpha-helix dipole to the pKa values and their ionic strength dependence was analyzed using Green's reciprocity theorem. The results indicate that (1) the pKa of His residues are largely affected by surrounding ionized and polar groups; (2) the ionic strength dependence of the pKa values is determined by the overall charge distributions and their accessibilities to solvent; and (3) alpha-helix dipoles make a significant contribution to the pKa, when the His residue is close to the helix terminus and not fully exposed to the solvent.     

pH dependence of stability of the 10th human fibronectin type III domain: a computational study

We present detailed computational studies based on electrostatic calculations to evaluate the origins of pKa values and the pH dependence of stability for the 10th type III domain of human fibronectin (FNfn10). One of our goals is to validate the calculation protocols by comparison to experimental data (Koide, A.; Jordan, M. R.; Horner, S.; Batori, V.; Koide, S. Biochemistry 2001, 40, 10326-10333). Another goal is to evaluate the sensitivity of the calculated ionization free energies and apparent pKa values on local structural fluctuations, which do not alter the structural convergence to a particular architecture, by using a complete ensemble of solution NMR structures and the NMR average minimized structure of FNfn10 (Main, A. L.; Harvey, T. S.; Baron, M.; Boyd, J.; Campbell, I. D. Cell 1992, 71, 671-678). Our calculations demonstrate that, at high ionic strength, FNfn10 is more stable at low pH compared to neutral pH, in overall agreement with experimental data. This behavior is attributed to contributions from unfavorable Coulombic interactions in a surface patch for the pairs Asp7-Glu9 and Asp7-Asp23. The unfavorable interactions are decreased at low pH, where the acidic residues become neutral, and are further decreased at high ionic strength because of increased screening by salt ions. Elimination of the unfavorable interactions in the theoretical mutants Asp7Asn (D7N) and Asp7Lys (D7K) produce higher calculated stabilities at neutral pH and any ionic strength compared to the wild-type, in agreement with the experimental data. We also discuss subtleties in the calculated apparent pKa values and ionization free energies, which are not in agreement with the experimental data. This work demonstrates that comparative electrostatic calculations can provide rapid predictions of pH-dependent properties of proteins and can be significant aids in guiding the design of proteins with tailored properties.     

Nanometric design of extraordinary hydrophobic-induced pKa shifts for aspartic acid: Relevance to protein mechanisms

Commonly a key element enabling proteins to function is an amino acid residue or residues with functional side chains having shifted pKa values. This article reports the results on a set of protein-based polymers (model proteins) that exhibit hydrophobic folding and assembly transitions, and that have been designed for the purpose of achieving large hydrophobic-induced pKa shifts by selectively replacing Val residues by Phe residues. The high molecular weight polypentapeptides, actually poly (tricosapeptides) with six varied pentamers in fixed sequence, were designed and synthesized to have the same amino acid compositions but different proximities between a single aspartic acid residue and 5 Phe residues per 30 residues. With the 5 Phe residues distal from the Asp residue, the observed pKa shift was 2.9 when compared to the Val-containing reference. With the 5 Phe residues within 1 nm of the Asp residue, the pKa shift was 6.2. This represents a free energy of interaction of 8 kcal/moles. To our knowledge, this is the largest pKa shift documented for an Asp residue in a polypeptide– or protein–water system. Data are reviewed that do not support the usual electrostatic arguments for pKa shifts of charge–charge repulsion and/or unfavorable ion self-energies arising from displacement of water by hydrophobic moieties, but rather the data are interpreted to indicate the presence of an apolar–polar repulsive free energy of hydration, which results from a potentially highly cooperative competition between apolar and polar species for hydration. © 1994 John Wiley & Sons, Inc.     

Catalytic activity of Ntau-carboxymethylhistidine-12 ribonuclease: pH dependence.

The pH-dependence of RNAase A and of Ntau-carboxymethylhistidine-12-RNAase (ribonucleate 3'-pyrimidino-oligonucleotidohydrolase) catalysis was studied. Apparent acid dissociation constants were obtained by least squares analysis of the kinetics data. These dissociation constants were compared with pKa values of model imidazole compounds, and with pKa values of histidine residues 12 and 119 on the protein. The shapes of the kcat versus pH profiles for RNAase A and its carboxymethyl derivative are very similar, from which it is concluded that the mechanism of catalysis is closely similar in the two proteins. Apparent pKa values obtained from the kinetic data are higher for the carboxymethylated protein than for RNAase A, as are the pKa values of residues 12 and 119. The similar shifts are consistent with the conclusions that both these residues are functionally significant in native and modified enzyme, and that an unblocked tau-nitrogen on histidine-12 is not essential for activity. From the enzyme's catalytic dependence on pH, and the NMR determined pKa values we propose that histidine 12 and 119 function catalytically in their basic and acidic forms respectively.     

Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability

Computational techniques based on continuum electrostatics treatments have been successful in predicting and interpreting the pKa values of ionizable amino acids in folded proteins. Despite this progress, efforts to reproduce the pH-dependence of protein stability have met with only limited success: agreement with experimental results has been only qualitative. It has been argued previously that the most likely reason for discrepancies is the presence of residual electrostatic interactions in the unfolded state, which cause pKa values to be shifted from their model compound values. Here we show that by constructing atomistic models of the unfolded state with a simple molecular mechanics protocol that uses the native state as a starting point, much improved reproduction of pH effects on protein stability can be obtained. In contrast, when a fully extended model of the unfolded state is used, no such improvement is obtained, a result that suggests that local interactions with residues nearby in the sequence are not sufficient to properly account for the pKa shifts in the unfolded state. In comparison to model compound values, the pKa values of acidic residues in "native-like" unfolded states are typically found to be shifted downwards by approximately 0.3 pH unit, in good agreement with the average downward shift deduced from experimental measurements. Given its success in the present situation, the protocol employed here for developing simple models of the unfolded state may prove useful in other computer simulation applications.     

A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins.

An improved approach is presented for calculating pH-dependent electrostatic effects in proteins using sigmoidally screened Coulomb potentials (SCP). It is hypothesized that a key determinant of seemingly aberrant behavior in pKa shifts is due to the properties of the unique microenvironment around each residue. To help demonstrate this proposal, an approach is developed to characterize the microenvironments using the local hydrophobicity/hydrophilicity around each residue of the protein. The quantitative characterization of the microenvironments shows that the protein is a complex mosaic of differing dielectric regions that provides a physical basis for modifying the dielectric screening functions: in more hydrophobic microenvironments the screening decreases whereas the converse applies to more hydrophilic regions. The approach was applied to seven proteins providing more than 100 measured pKa values and yielded a root mean square deviation of 0.5 between calculated and experimental values. The incorporation of the local hydrophobicity characteristics into the algorithm allowed the resolution of some of the more intractable problems in the calculation of pKa. Thus, the divergent shifts of the pKa of Glu-35 and Asp-66 in hen egg white lysozyme, which are both about 90% buried, was correctly predicted. Mechanistically, the divergence occurs because Glu-35 is in a hydrophobic microenvironment, while Asp-66 is in a hydrophilic microenvironment. Furthermore, because the calculation of the microenvironmental effects takes very little CPU time, the computational speed of the SCP formulation is conserved. Finally, results from different crystal structures of a given protein were compared, and it is shown that the reliability of the calculated pKa values is sufficient to allow identification of conformations that may be more relevant for the solution structure.     

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.     

Electrostatics in the active site of an alpha-amylase.

The importance of electrostatics in catalysis has been emphasized in the literature for a large number of enzymes. We examined this hypothesis for the Bacillus licheniformis alpha-amylase by constructing site-directed mutants that were predicted to change the pKa values of the catalytic residues and thus change the pH-activity profile of the enzyme. To change the pKa of the catalytic residues in the active site, we constructed mutations that altered the hydrogen bonding network, mutations that changed the solvent accessibility, and mutations that altered the net charge of the molecule. The results show that changing the hydrogen bonding network near an active site residue or changing the solvent accessibility of an active site residue will very likely result in an enzyme with drastically reduced activity. The differences in the pH-activity profiles for these mutants were modest. pH-activity profiles of mutants which change the net charge on the molecule were significantly different from the wild-type pH-activity profile. The differences were, however, difficult to correlate with the electrostatic field changes calculated. In several cases we observed that pH-activity profiles shifted in the opposite direction compared to the shift predicted from electrostatic calculations. This strongly suggests that electrostatic effects cannot be solely responsible for the pH-activity profile of the B. licheniformis alpha-amylase.     

Free energy of transition for the individual alkaline conformers of yeast iso-1-cytochrome c

Direct protein electrochemistry was used to obtain the thermodynamic parameters of transition from the native (state III) to the alkaline (state IV) conformer for untrimethylated Saccharomyces cerevisiae iso-1-cytochrome c expressed in E. coli and its single and multiple lysine-depleted variants. In these variants, one or more of the lysine residues involved in axial Met substitution (Lys72, Lys73, and Lys79) was mutated to alanine. The aim of this work is to determine the thermodynamic affinity of each of the substituting lysines for the heme iron and evaluate the interplay of enthalpic and entropic factors. The equilibrium constants for the deprotonation reaction of Lys72, 73, and 79 were computed for the minimized MD average structures of the wild-type and mutated proteins, applying a modified Tanford-Kirkwood calculation. Solvent accessibility calculations for the substituting lysines in all variants were also performed. The transition enthalpy and entropy values within the protein series show a compensatory behavior, typical of a process involving extensive solvent reorganization effects. The experimental and theoretical data indicate that Lys72 most readily deprotonates and replaces M80 as the axial heme iron ligand, whereas Lys73 and Lys79 show comparably higher pKa values and larger transition free energies. A good correlation is found within the series between the lowest calculated Lys pKa value and the corresponding experimental pKa value, which can be interpreted as indicative of the deprotonating lysine itself acting as the triggering group for the conformational transition. The triple Lys to Ala mutant, in which no lysine residues are available for heme iron binding, features transition thermodynamics consistent with a hydroxide ion replacing the axial methionine ligand.     

Calculating pKa values in enzyme active sites

The ionization properties of the active-site residues in enzymes are of considerable interest in the study of the catalytic mechanisms of enzymes. Knowledge of these ionization constants (pKa values) often allows the researcher to identify the proton donor and the catalytic nucleophile in the reaction mechanism of the enzyme. Estimates of protein residue pKa values can be obtained by applying pKa calculation algorithms to protein X-ray structures. We show that pKa values accurate enough for identifying the proton donor in an enzyme active site can be calculated by considering in detail only the active-site residues and their immediate electrostatic interaction partners, thus allowing for a large decrease in calculation time. More specifically we omit the calculation of site-site interaction energies, and the calculation of desolvation and background interaction energies for a large number of pairs of titratable groups. The method presented here is well suited to be applied on a genomic scale, and can be implemented in most pKa calculation algorithms to give significant reductions in calculation time with little or no impact on the accuracy of the results. The work presented here has implications for the understanding of enzymes in general and for the design of novel biocatalysts.     

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.     

Structure-based pKa prediction provides a thermodynamic basis for the role of histidines in pH-induced conformational transitions in dengue virus

pH-induced conformational changes in dengue virus (DENV) are critical to its ability to infect host cells. The envelope protein heterodimers that make up the viral envelope shift from a dimer to a trimer conformation at low-pH during membrane fusion. Previous studies have suggested that the ionization of histidine residues at low-pH is central to this pH-induced conformational change. We sought out to use molecular modeling with structure-based pKa prediction to provide a quantitative basis for the role of histidines in pH-induced conformational changes and identify which histidine residues were primarily responsible for this transition. We combined existing crystallographic and cryo-electron microscopy data to construct templates of the dimer and trimer conformations for the mature and immature virus. We then generated homology models for the four DENV serotypes and carried out structure-based pKa prediction using Rosetta. Our results showed that the pKa values of a subset of conserved histidines in DENV successfully capture the thermodynamics necessary to drive pH-induced conformational changes during fusion. Here, we identified the structural determinants underlying these pKa values and compare our findings with previous experimental results.     

Apparent pKa shifts of titratable residues at high denaturant concentration and the impact on protein stability

Urea and guanidine-hydrochloride (GdnHCl) are frequently used for protein denaturation in order to determine the Gibbs free energy of folding and kinetic folding/unfolding parameters. Constant pH value is applied in the folding/unfolding experiments at different denaturant concentrations and steady protonation state of titratable groups is assumed in the folded and unfolded protein, respectively. The apparent side-chain pKa values of Asp, Glu, His and Lys in the absence and presence of 6 M urea and GdnHCl, respectively, have been determined by 1H-NMR. pKa values of all four residues are up-shifted by 0.3-0.5 pH units in presence of 6 M urea by comparison with pKa values of the residues dissolved in water. In the presence of 6 M GdnHCl, pKa values are down-shifted by 0.2-0.3 pH units in the case of acidic and up-shifted by 0.3-0.5 pH units in the case of basic residues. Shifted pKa values in the presence of denaturant may have a pronounced effect on the outcome of the protein stability obtained from denaturant unfolding experiments.     

The evolutionary origins and catalytic importance of conserved electrostatic networks within TIM-barrel proteins

Conservation of function is the basic tenet of protein evolution. Conservation of key electrostatic properties is a frequently employed mechanism that leads to conserved function. In a previous report, we identified several conserved electrostatic properties in four protein families and one functionally diverse enzyme superfamily. In this report, we demonstrate the evolutionary and catalytic importance of electrostatic networks in three ubiquitous metabolic enzymes: triosephosphate isomerase, enolase, and transaldolase. Evolutionary importance is demonstrated using phylogenetic motifs (sequence fragments that parallel the overall familial phylogeny). Phylogenetic motifs frequently correspond to both catalytic residues and conserved interactions that fine-tune catalytic residue pKa values. Further, in the case of triosephosphate isomerase, quantitative differences in the catalytic Glu169 pKa values parallel subfamily differentiation. Finally, phylogenetic motifs are shown to structurally cluster around the active sites of eight different TIM-barrel families. Depending upon the mechanistic requisites of each reaction catalyzed, interruptions to the canonical fold may or may not be identified as phylogenetic motifs.     

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.     

Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.     

Nuclear magnetic resonance titration curves of histidine ring protons. Human metmyoglobin and the effects of azide on human, horse, and sperm whale metmyoglobins.

Four titrating histidine ring C2 and C4 proton resonances are observed in 220 MHz proton NMR spectra of human metmyoglobin as a function of pH. Values of ionization constants determined from the NMR titration data using an equation describing a simple proton association-dissociation equilibrium are curves (1) 6.6, (2) 7.0, (3) 5.8, and (4) 7.4. Four histidine residues have also been found to be solvent-accessible in human metmyoglobin by carboxymethylation studies (Harris, C.M., and Hill, R.L. (1969) J. Biol. Chem. 244, 2195-2203). Two of the titration curves (3 and 4) deviate significantly from the chemical shift values normally observed for histidine C2 proton resonances. Curve 3, with a low pKa, is shifted downfield at high values of pH and also exhibits a second minor inflection with a pKa value of 8.8. On the other hand, the high pKa curve, 4, is shifted upfield at all values of pH. The characteristics of the NMR titration curves with the lowest and highest pKa values (3 and4) are very similar to curves observed previously with sperm whale and horse metmyoglobins (Cohen, J.S., Hagenmaier, H., Pollard, H., and Schechter, A.N. (1972) J. Mol. Biol. 71, 513-519). These results indicate that the histidine residues from which these curves are derived have unusual and characteristic environments in this series of homologous proteins. The NMR spectra of all three metmyoglobins are changed extensively as a result of azide ion binding, indicating conformational changes affecting the environments of several imidazole side chains. The presence of azide ion causes a selective downfield chemical shift for the low pKa curve and a selective upfield chemical shift for the high pKa curve in all three proteins. Azide also abolishes the second inflection seen in the low pKa curve at high pH. In addition to these effects, the presence of azide ion permits the observation of two additional titrating proton resonances for all three metmyoglobins. Increasing the azide to protein ratio at several fixed values of pH yields results which show that a slow exchange process is occurring with each of the metmyoglobins. In the azide titration studies the maximum changes in the NMR spectra occurred at approximately equimolar concentrations. The NMR results for these proteins in the absence and presence of azide ion are related to x-ray crystallographic studies of sperm whale metmyoglobin and the known alkylation properties of the histidine residues. Tentative assignments of the titrating resonances observed are suggested.     

Stabilization of a protein by removal of unfavorable abnormal pKa: substitution of undissociable residue for glutamic acid-35 in chicken lysozyme.

Glu35 in chicken lysozyme has an abnormally high pKa (6.1) partly due to the hydrophobic environment provided by Trp108. The relationship between protein stability and abnormal pKa was investigated in detail by using mutant lysozymes in which Glu35 was replaced by undissociable residues and an oppositely ionizable residue. It was found that lysozyme was stabilized at alkaline pH range by the replacement of Glu35 with an undissociable residue, Gln (E35Q lysozyme) or Al (E35A lysozyme). On the other hand, when Glu35 was replaced by His (E35H lysozyme), which could have an opposite charge to Glu by ionization, the introduced His35 was found to have an abnormally low pKa (3.6), leading to the destabilization of lysozyme at acidic pH. These observations are completely consistent with the situation that the environment around Glu35 is highly hydrophobic and therefore the placement of either a positive or negative charge in such an environment leads to destabilization of lysozyme. These observations also indicate that the replacement of an acidic residue having abnormally high pKa or a basic residue having abnormally low pKa by an undissociable residue is a very efficient and general method for stabilization of a protein.