Search

Recent advances in structural biology revealed that water molecules play a crucial structural role in the protein architecture and ligand binding of G protein-coupled receptors. In this work, we present an alternative approach to monitor the time-dependent organization of water molecules during the final stage of the ligand-receptor recognition process by means of membrane molecular dynamics simulations. We inspect the variation of fluid dynamics properties of water molecules upon ligand binding with the aim to correlate the results with the binding affinities. The outcomes of this analysis are transferred into a bidimensional graph called water fluid dynamics maps, that allow a fast graphical identification of protein "hot-spots" characterized by peculiar shape and electrostatic properties that can play a critical role in ligand binding. We hopefully believe that the proposed approach might represent a valuable tool for structure-based drug discovery that can be extended to cases where crystal structures are not yet available, or have not been solved at high resolution.

G protein-coupled receptors (GPCRs) represent the largest family of cell-surface receptors and about one-third of the actual targets of clinically used drugs. Following the progress made in the field of GPCRs structural determination, docking-based screening for novel potent and selective ligands is becoming an increasingly adopted strategy in the drug discovery process. However, this methodology is not yet able to anticipate the "bioactive" binding mode and discern it among other conformations. In the present work, we present a novel approach consisting in the integration of molecular docking and membrane MD simulations with the aim to merge the rapid sampling of ligand poses into in the binding site, typical of docking algorithms, with the thermodynamic accuracy of MD simulations in describing, at the molecular level, the stability a GPCR-ligand complex embedded into explicit lipid-water environment. To validate our approach, we have chosen as a key study the human A(2A) adenosine receptor (hA(2A) AR) and selected four receptor-antagonist complexes and one receptor-agonist complex that have been recently crystallized. In light of the obtained results, we believe that our novel strategy can be extended to other GPCRs and might represent a valuable tool to anticipate the "bioactive" conformation of high-affinity ligands.