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 pK(a) 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 pK(a) 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 pK(a). Thus, the divergent shifts of the pK(a) 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 pK(a) values is sufficient to allow identification of conformations that may be more relevant for the solution structure.