The role of hydrophobic microenvironments in modulating pKa shifts in proteins

The determinants of carboxyl pKa values in turkey ovomucoid third domain

A computational methodology for protein pK(a) predictions, based on ab initio quantum mechanical treatment of part of the protein and linear Poisson-Boltzmann equation treatment of the bulk solvent, is presented. The method is used to predict and interpret the pK(a) values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in the serine protease inhibitor turkey ovomucoid third domain. All the predicted pK(a) values are within 0.5 pH units of experiment, with a root-mean-square deviation of 0.31 pH units. We show that the decreased pK(a) values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrogen bonds involving amide protons are shown to be particularly important, and the effect of hydrogen bonding is shown to be nonadditive. Hydrophobic effects are also shown to be important in raising the pK(a). Interactions with charged residues are shown to have relatively little effect on the carboxyl pK(a) values in this protein, in general agreement with experiment.

A mixed mechanistic-electrostatic model to explain pH dependence of glycosyl hydrolase enzyme activity

Glycosyl hydrolases are a vast group of enzymes that share a common topology at their active site with two acidic residues that are responsible for activity. In spite of their similarity, they exhibit a wide range of pH optima that must depend on other factors. Using structural and mechanistic knowledge about glycosyl hydrolases from families 7, 10, and 16, we have formulated a new mathematical model that can include not only the ionization behavior of the catalytic residues but also as many ionizable residues as desired in the active site. In addition, the model can incorporate electrostatic influences via acid dissociation equilibrium constants and chemical relationships such as hydrogen bonds. The results of the simulations indicate a clear shift in the pH dependence of activity for the enzymes only when a close interrelation (hydrogen bond) between the catalytic and auxiliary residues in the active site is taken into account. This explains the observations from mutagenesis studies that show this type of shift and cannot be explained by a purely electrostatic interaction theory. Moreover, the presence of the kind of chemical interaction proposed could provide stabilization of the activity in the presence of environmental, structural, pH and electrostatic variations. These findings and the implications for the design of new mutagenesis strategies are discussed. The results suggest a way to modify, via site-directed mutagenesis, the acid dissociation of one of the catalytic residues in the active site independently of the other, which could have clear advantages over the purely electrostatic modifications that usually affect both residues simultaneously.

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.

Role of the protein side-chain fluctuations on the strength of pair-wise electrostatic interactions: comparing experimental with computed pK(a)s

The effect of the protein side-chain fluctuations on the strength of electrostatic interactions was studied. The effect was modeled on 7 different crystal structures on the same enzyme as well as on 20 molecular dynamics snapshot structures. It was shown that the side-chain flexibility affects predominantly the magnitude of the strong pair-wise interactions, that is, the pair-wise interaction among ion pairs, and practically does not affect the interactions with the rest of the protein. This was used to suggest a correction function that should be applied to the original pair-wise electrostatic interaction to mimic the effects of the fluctuations. The procedure is applied on three ion pairs identified in lysozyme. It was shown that sampling different side-chain rotamers and modifying the strength of the pair-wise interaction energies makes calculated pK(a)s less sensitive to the fluctuations of the structure and improves the prediction accuracy.

Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration

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.