Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites

Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity

Asymmetric structural motions of the homomeric alpha 7 nicotinic receptor ligand binding domain revealed by molecular dynamics simulation

AbstractA homology model of the ligand binding domain of the α7 nicotinic receptor is constructed based on the acetylcholine-binding protein crystal structure. This structure is refined in a 10ns molecular dynamics simulation. The modeled structure proves fairly resilient, with no significant changes at the secondary or tertiary structural levels. The hypothesis that the acetylcholine-binding protein template is in the activated or desensitized state, and the absence of a bound agonist in the simulation suggests that the structure may also be relaxing from this state to the activatable state. Candidate motions that take place involve not only the side chains of residues lining the binding sites, but also the subunit positions that determine the overall shape of the receptor. In particular, two nonadjacent subunits move outward, whereas their partners counterclockwise to them move inward, leading to a marginally wider interface between themselves and an overall asymmetric structure. This in turn affects the binding sites, producing two that are more open and characterized by distinct side-chain conformations of W54 and L118, although motions of the side chains of all residues in every binding site still contribute to a reduction in binding site size, especially the outward motion of W148, which hinders acetylcholine binding. The Cys loop at the membrane interface also displays some flexibility. Although the short simulation timescale is unlikely to sample adequately all the conformational states, the pattern of observed motions suggests how ligand binding may correlate with larger-scale subunit motions that would connect with the transmembrane region that controls the passage of ions. Furthermore, the shape of the asymmetry with binding sites of differing affinity for acetylcholine, characteristic of other nicotinic receptors, may be a natural property of the relaxed, activatable state of α7

Protein Dynamics of the HIF-2α PAS-B Domain upon Heterodimerization and Ligand Binding

Hypoxia-Inducible Factor (HIF) transcription factors are heterodimeric proteins involved in the regulation of oxygen homeostatis. Their upregulation has been related to several tumors with a remarkably poor clinical outcome. The recent discovery of a druggable cavity in the HIF-2α PAS-B domain has opened an unprecedented opportunity for targeting the HIF-2α transcription factor in view of pharmaceutical strategies. Coincidentally, a novel compound able to selectively disrupt the HIF heterodimerization with a submicromolar activity has been reported. In this work, we investigated the molecular mechanisms responsible for the inhibition by comparing the dynamical features of the HIF-2α PAS-B monomer and the HIF-2α PAS-B/HIF-1β PAS-B complex, in the ligand-bound and -unbound states. Plain and biased Molecular Dynamics were used to characterize the differential conformational changes both structurally and energetically