Demystifying P2Y1 Receptor Ligand Recognition through Docking and Molecular Dynamics Analyses

We performed a molecular modeling analysis of 100 nucleotide-like bisphosphates and 46 non-nucleotide arylurea derivatives previously reported as P2Y1R binders using the recently solved hP2Y1R structures. We initially docked the compounds at the X-ray structures and identified the binding modes of representative compounds highlighting key patterns in the structure–activity relationship (SAR). We subsequently subjected receptor complexes with selected key agonists (2MeSADP and MRS2268) and antagonists (MRS2500 and BPTU) to membrane molecular dynamics (MD) simulations (at least 200 ns run in triplicate, simulation time 0.6–1.6 μs per ligand system) while considering alternative protonation states of nucleotides. Comparing the temporal evolution of the ligand–protein interaction patterns with available site-directed mutagenesis (SDM) data and P2Y1R apo state simulation provided further SAR insights and suggested reasonable explanations for loss/gain of binding affinity as well as the most relevant charged species for nucleotide ligands. The MD analysis also predicted local conformational changes required for the receptor inactive state to accommodate nucleotide agonists.