A glutamic acid was buried in the hydrophobic core of staphylococcal nuclease by replacement of Val-66. Its pK(a) was measured with equilibrium thermodynamic methods. It was 4.3 units higher than the pK(a) of Glu in water. This increase was comparable to the DeltapK(a) of 4.9 units measured previously for a lysine buried at the same location. According to the Born formalism these DeltapK(a) are energetically equivalent to the transfer of a charged group from water to a medium of dielectric constant of 12. In contrast, the static dielectric constants of dry protein powders range from 2 to 4. In the crystallographic structure of the V66E mutant, a chain of water molecules was seen that hydrates the buried Glu-66 and links it with bulk solvent. The buried water molecules have never previously been detected in >20 structures of nuclease. The structure and the measured energetics constitute compelling and unprecedented experimental evidence that solvent penetration can contribute significantly to the high apparent polarizability inside proteins. To improve structure-based calculations of electrostatic effects with continuum methods, it will be necessary to learn to account quantitatively for the contributions by solvent penetration to dielectric effects in the protein interior.

Lys-66 and Glu-66, buried in the hydrophobic interior of staphylococcal nuclease by mutagenesis, titrate with pK(a) values of 5.7 and 8.8, respectively (Dwyer et al., Biophys. J. 79:1610-1620; García-Moreno E. et al., Biophys. Chem. 64:211-224). Continuum calculations with static structures reproduced the pK(a) values when the protein interior was treated with a dielectric constant (epsilon(in)) of 10. This high apparent polarizability can be rationalized in the case of Glu-66 in terms of internal water molecules, visible in crystallographic structures, hydrogen bonded to Glu-66. The water molecules are absent in structures with Lys-66; the high polarizability cannot be reconciled with the hydrophobic environment surrounding Lys-66. Equilibrium thermodynamic experiments showed that the Lys-66 mutant remained folded and native-like after ionization of the buried lysine. The high polarizability must therefore reflect water penetration, minor local structural rearrangement, or both. When in pK(a) calculations with continuum methods, the internal water molecules were treated explicitly, and allowed to relax in the field of the buried charged group, the pK(a) values of buried residues were reproduced with epsilon(in) in the range 4-5. The calculations show that internal waters can modulate pK(a) values of buried residues effectively, and they support the hypothesis that the buried Lys-66 is in contact with internal waters even though these are not seen crystallographically. When only the one or two innermost water molecules were treated explicitly, epsilon(in) of 5-7 reproduced the pK(a) values. These values of epsilon(in) > 4 imply that some conformational reorganization occurs concomitant with the ionization of the buried groups.

No abstract

The structure of a complex of staphylococcal nuclease with Ca2+ and deoxythymidine 3',5'-bisphosphate (pdTp) has been refined by stereochemically restrained least-squares minimization to a crystallographic R value of 0.161 at 1.65 A resolution. The estimated root-mean-square (rms) error in the coordinates is 0.16 A. The final model comprises 1082 protein atoms, one calcium ion, the pdTp molecule, and 82 solvent water molecules; it displays an rms deviation from ideality of 0.017 A for bond distances and 1.8 degrees for bond angles. The mean distance between corresponding alpha carbons in the refined and unrefined structures is 0.6 A; we observe small but significant differences between the refined and unrefined models in the turn between residues 27 and 30, the loop between residues 44 and 50, the first helix, and the extended strand between residues 112 and 117 which forms part of the active site binding pocket. The details of the calcium liganding and solvent structure in the active site are clearly shown in the final electron density map. The structure of the catalytic site is consistent with the mechanism that has been proposed for this enzyme. However, we note that two lysines from a symmetry-related molecule in the crystal lattice may play an important role in determining the geometry of inhibitor binding, and that only one of the two required calcium ions is observed in the crystal structure; thus, caution is advised in extrapolating from the structure of the complex of enzyme and inhibitor to that of enzyme and substrate.