Dataset Information

Probing the role of interfacial waters in protein-DNA recognition using a hybrid implicit/explicit solvation model.

ABSTRACT: When proteins bind to their DNA target sites, ordered water molecules are often present at the protein-DNA interface bridging protein and DNA through hydrogen bonds. What is the role of these ordered interfacial waters? Are they important determinants of the specificity of DNA sequence recognition, or do they act in binding in a primarily nonspecific manner, by improving packing of the interface, shielding unfavorable electrostatic interactions, and solvating unsatisfied polar groups that are inaccessible to bulk solvent? When modeling details of structure and binding preferences, can fully implicit solvent models be fruitfully applied to protein-DNA interfaces, or must the individualistic properties of these interfacial waters be accounted for? To address these questions, we have developed a hybrid implicit/explicit solvation model that specifically accounts for the locations and orientations of small numbers of DNA-bound water molecules, while treating the majority of the solvent implicitly. Comparing the performance of this model with that of its fully implicit counterpart, we find that explicit treatment of interfacial waters results in a modest but significant improvement in protein side-chain placement and DNA sequence recovery. Base-by-base comparison of the performance of the two models highlights DNA sequence positions whose recognition may be dependent on interfacial water. Our study offers large-scale statistical evidence for the role of ordered water for protein-DNA recognition, together with detailed examination of several well-characterized systems. In addition, our approach provides a template for modeling explicit water molecules at interfaces that should be extensible to other systems.

SUBMITTER: Li S

PROVIDER: S-EPMC3904374 | biostudies-literature | 2013 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Probing the role of interfacial waters in protein-DNA recognition using a hybrid implicit/explicit solvation model.

Li Shen S Bradley Philip P

Proteins 20130422 8

When proteins bind to their DNA target sites, ordered water molecules are often present at the protein-DNA interface bridging protein and DNA through hydrogen bonds. What is the role of these ordered interfacial waters? Are they important determinants of the specificity of DNA sequence recognition, or do they act in binding in a primarily nonspecific manner, by improving packing of the interface, shielding unfavorable electrostatic interactions, and solvating unsatisfied polar groups that are in ...[more]

PMID: 23444044

Similar Datasets

Project description:Accurate yet efficient computational models of solvent environment are central for most calculations that rely on atomistic modeling, such as prediction of protein-ligand binding affinities. In this study, we evaluate the accuracy of a recently developed generalized Born implicit solvent model, GBNSR6 (Aguilar et al. J. Chem. Theory Comput. 2010, 6, 3613-3639), in estimating the electrostatic solvation free energies (ΔG(pol)) and binding free energies (ΔΔG(pol)) for small protein-ligand complexes. We also compare estimates based on three different explicit solvent models (TIP3P, TIP4PEw, and OPC). The two main findings are as follows. First, the deviation (RMSD = 7.04 kcal/mol) of GBNSR6 binding affinities from commonly used TIP3P reference values is comparable to the deviations between explicit models themselves, e.g. TIP4PEw vs TIP3P (RMSD = 5.30 kcal/mol). A simple uniform adjustment of the atomic radii by a single scaling factor reduces the RMS deviation of GBNSR6 from TIP3P to within the above "error margin" - differences between ΔΔG(pol) estimated by different common explicit solvent models. The simple radii scaling virtually eliminates the systematic deviation (ΔΔG(pol)) between GBNSR6 and two out of the three explicit water models and significantly reduces the deviation from the third explicit model. Second, the differences between electrostatic binding energy estimates from different explicit models is disturbingly large; for example, the deviation between TIP4PEw and TIP3P estimates of ΔΔG(pol) values can be up to ∼50% or ∼9 kcal/mol, which is significantly larger than the "chemical accuracy" goal of ∼1 kcal/mol. The absolute ΔG(pol) calculated with different explicit models could differ by tens of kcal/mol. These discrepancies point to unacceptably high sensitivity of binding affinity estimates to the choice of common explicit water models. The absence of a clear "gold standard" among these models strengthens the case for the use of accurate implicit solvation models for binding energetics, which may be orders of magnitude faster.

Project description:Our group has implemented a smooth Gaussian-based dielectric function in DelPhi (J. Chem. Theory Comput. 2013, 9 (4), 2126-2136) which models the solute as an object with inhomogeneous dielectric permittivity and provides a smooth transition of dielectric permittivity from surface-bound water to bulk solvent. Although it is well-understood that the protein hydrophobic core is less polarizable than the hydrophilic protein surface, less attention is paid to the polarizability of water molecules inside the solute and on its surface. Here, we apply explicit water simulations to study the behavior of water molecules buried inside a protein and on the surface of that protein and contrast it with the behavior of the bulk water. We selected a protein that is experimentally shown to have five cavities, most of which are occupied by water molecules. We demonstrate through molecular dynamics (MD) simulations that the behavior of water in the cavity is drastically different from that in the bulk. These observations were made by comparing the mean residence times, dipole orientation relaxation times, and average dipole moment fluctuations. We also show that the bulk region has a nonuniform distribution of these tempo-spatial properties. From the perspective of continuum electrostatics, we argue that the dielectric "constant" in water-filled cavities of proteins and the space close to the molecular surface should differ from that assigned to the bulk water. This provides support for the Gaussian-based smooth dielectric model for solving electrostatics in the Poisson-Boltzmann equation framework. Furthermore, we demonstrate that using a well-parametrized Gaussian-based model with a single energy-minimized configuration of a protein can also reproduce its ensemble-averaged polar solvation energy. Thus, we argue that the Gaussian-based smooth dielectric model not only captures accurate physics but also provides an efficient way of computing ensemble-averaged quantities.

Dataset Information

Probing the role of interfacial waters in protein-DNA recognition using a hybrid implicit/explicit solvation model.

Publications

Probing the role of interfacial waters in protein-DNA recognition using a hybrid implicit/explicit solvation model.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets