Experimental measurement of the effective dielectric in the hydrophobic core of a protein.

The dielectric inside a protein is a key physical determinant of the magnitude of electrostatic interactions in proteins. We have measured this dielectric phenomenologically, in terms of the dielectric that needs to be used with the Born equation in order to reproduce the observed pKa shifts induced by burial of an ionizable group in the hydrophobic core of a protein. Mutants of staphylococcal nuclease with a buried lysine residue at position 66 were engineered for this purpose. The pKa values of buried lysines were measured by difference potentiometry. The extent of coupling between the pKa and the global stability of the protein was evaluated by measuring pKa values in hyperstable forms of nuclease engineered to be 3.3 or 6.5 kcal mol-1 more stable than the wild type. The crystallographic structure of one mutant was determined to describe the environment of the buried lysine. The dielectrics that were measured range from 10 to 12. Published pKa values of buried ionizable residues in other proteins were analyzed in a similar fashion and the dielectrics obtained from these values are consistent with the ones measured in nuclease. These results argue strongly against the prevalent use of dielectrics of 4 or lower to describe the dielectric effect inside a protein in structure-based calculations of electrostatic energies with continuum dielectric models.

Contributions of the ionizable amino acids to the stability of staphylococcal nuclease.

To quantitate the contributions of the ionizable amino acids to the stability of the native state of staphylococcal nuclease, each of the 23 lysines, 5 arginines, 4 histidines, 12 glutamic acids, and 8 aspartic acids was substituted with both alanine and glycine. This collection of 104 mutant proteins was analyzed by guanidine hydrochloride (GuHCl) denaturation, using intrinsic tryptophan fluorescence to quantitate the equilibrium between native and denatured states. From the analysis of these data, each mutant protein's stability in the absence of denaturant (delta GH2O) and sensitivity to changes in denaturant concentration [mGuHCl = d(delta G)/d[GuHCl]] were obtained. Several general trends in these values suggest that electrostatic interactions make only a minor contribution to the net stability of this protein. For the residue pairs that form ten salt bridges and ten charged hydrogen bonds between side chains, no correlation was observed between the stability losses (delta delta G) accompanying alanine substitution of each member of the pair. Little or no significant correlation was found between the magnitude of the loss in stability and the local electrostatic potential calculated from the three-dimensional structure by numerical and model dependent solutions of the linearized Poisson-Boltzmann equation. The structural parameters which correlated most strongly with stability loss are measures of the extent of burial of the residue in the native structure, as was previously observed for alanine and glycine substitutions of large hydrophobic residues [Shortle et al. (1990) Biochemistry 29, 8033] and of the polar, uncharged residues [Green et al. (1992) Biochemistry 31, 5717]. These results suggest that the ionizable amino acids contribute to stability predominantly through packing and bonding interactions that do not depend on their electrostatic charge.

A calorimetric characterization of the salt dependence of the stability of the GCN4 leucine zipper.

The effects of different salts (LiCl, NaCl, ChoCl, KF, KCl, and KBr) on the structural stability of a 33-residue peptide corresponding to the leucine zipper region of GCN4 have been studied by high-sensitivity differential scanning calorimetry. These experiments have allowed an estimation of the salt dependence of the thermodynamic parameters that define the stability of the coiled coil. Independent of the nature of the salt, a destabilization of the coiled coil is always observed upon increasing salt concentration up to a maximum of approximately 0.5 M, depending on the specific cation or anion. At higher salt concentrations, this effect is reversed and a stabilization of the leucine zipper is observed. The effect of salt concentration is primarily entropic, judging from the lack of a significant salt dependence of the transition enthalpy. The salt dependence of the stability of the peptide is complex, suggesting the presence of specific salt effects at high salt concentrations in addition to the nonspecific electrostatic effects that are prevalent at lower salt concentrations. The data is consistent with the existence of specific interactions between anions and peptide with an affinity that follows a reverse size order (F- > Cl- > Br-). Under all conditions studied, the coiled coil undergoes reversible thermal unfolding that can be well represented by a reaction of the form N2<==>2U, indicating that the unfolding is a two-state process in which the helices are only stable when they are in the coiled coil conformation.

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions.

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.

pH-dependent processes in proteins.

Recent improvements in the understanding of electrostatic interactions in proteins serve as a focus for the general topic of pH-dependent processes in proteins. The general importance of pH-dependent processes is first set out in terms of hydrogen ion equilibria, stability, ligand interactions, assembly, dynamics, and events in related molecular systems. The development of various theoretical treatments includes various formalisms in addition to the solvent interface model developed by Shire et al. as an extension of the Tanford-Kirkwood treatment. A number of detailed applications of the model are presented and future potentialities are sketched.

Contributions of individual amino acid residues to the structural stability of cetacean myoglobins.

The acid-denaturation behavior of eleven cetacean myoglobins has been studied at two ionic strengths, 0.01 and 0.10 M, at 25.0 degrees C. The myoglobulins studied fall into four phylogenetic suborders, representing the sperm whales, dolphins, baleen whales, and beaked whales. The differences in response to acid denaturation among these closely related myoglobins are small but statistically significant. In three cases, free-energy differences between myoglobins can be ascribed to one amino acid difference and in three others to two differences. The differences in response were analyzed in terms of the changes in noncovalent interactions occurring in the native structure. The effects of changes in electrostatic interactions over the whole charge array were calculated for each myoglobin species by using the modified Tanford-Kirkwood theory. The predicted changes in stability correlated well with the experimental observations in most cases. When differences in hydrogen-bonding capability were considered at a first approximation, substantial effects were predicted. When these effects were taken in conjunction with the electrostatic interactions, the correlation with experiment was improved. Additionally, restrictions in motional freedom and packing constraints appeared to be significant in the single-site analysis. The detectable differences in stability due to single amino acid substitutions along with the small differences in stability between the cetacean suborders indicate that compensatory interactions provide the mechanism for the conservation of stability among the myoglobins studied.