Revisiting a receptor-based pharmacophore hypothesis for human A(2A) adenosine receptor antagonists

Identification of new potent A1 adenosine receptor antagonists using a multistage virtual screening approach

The use of antagonists for each adenosine receptor (AR) subtype as potent clinical candidates is of growing interest due to their involvement in the treatment of various diseases. The recent resolution of several A₁ and A_2A ARs X-ray structures provides opportunities for structure-based drug design. In this study, we describe the discovery of novel A₁AR antagonists by applying a multistage virtual screening approach, which is based on random forest (RF), e-pharmacophore modeling and docking methods. A multistage virtual screening approach was applied to screen the ChemDiv library (1,492,362 compounds). Among the final hits, 22 compounds were selected for further radioligand binding assay analysis against human A₁AR, and 18 compounds (81.82% success) exhibited nanomolar or low micromolar binding potency (K_i). Then, we selected six compounds (pK_i > 6) to further evaluate their antagonist profile in a cAMP functional assay, and we found that they had low micromolar antagonistic activity (pIC₅₀ = 5.51-6.38) for the A₁AR. Particularly, four of six compounds (pK_i > 6) showed very good affinity (pK_i = 6.11-7.13) and selectively (>100-fold) for A₁AR over A_2AAR. Moreover, the novelty analysis suggested that four of six compounds (pK_i > 6) were dissimilar to existing A₁AR antagonists and hence represented novel A₁AR antagonists. Further molecular docking and molecular dynamics (MD) studies showed that the three selective compounds 15, 20 and 22 were stabilized (RMS_lig value ≤ 2 Å) inside the binding pocket of A₁AR with similar orientations to the docking pose in 100-ns MD simulations, whereas they escaped from the binding area of A_2AAR with larger values of RMS_lig (RMS_lig ≥ 2 Å). We hope that these findings provide new insights into the discovery of drugs targeting A₁AR and facilitate research on new drugs and treatments for A₁AR-related human pathologies.

Breakthrough in GPCR Crystallography and Its Impact on Computer-Aided Drug Design

Recent crystallographic structures of G protein-coupled receptors (GPCRs) have greatly advanced our understanding of the recognition of their diverse agonist and antagonist ligands. We illustrate here how this applies to A_2A adenosine receptors (ARs) and to P2Y₁ and P2Y₁₂ receptors (P2YRs) for ADP. These X-ray structures have impacted the medicinal chemistry aimed at discovering new ligands for these two receptor families, including receptors that have not yet been crystallized but are closely related to the known structures. In this Chapter, we discuss recent structure-based drug design projects that led to the discovery of: (a) novel A₃AR agonists based on a highly rigidified (N)-methanocarba scaffold for the treatment of chronic neuropathic pain and other conditions, (b) fluorescent probes of the ARs and P2Y₁₄R, as chemical tools for structural probing of these GPCRs and for improving assay capabilities, and (c) new more drug-like antagonists of the inflammation-related P2Y₁₄R. We also describe the computationally enabled molecular recognition of positive (for A₃AR) and negative (P2Y₁R) allosteric modulators that in some cases are shown to be consistent with structure-activity relationship (SAR) data. Thus, computational modeling has become an essential tool for the design of purine receptor ligands.

Multivariate statistical analysis methods in QSAR

The emphasis of this review is particularly on multivariate statistical methods currently used in quantitative structure–activity relationship (QSAR) studies.

John Daly Lecture: Structure-guided Drug Design for Adenosine and P2Y Receptors

We establish structure activity relationships of extracellular nucleosides and nucleotides at G protein-coupled receptors (GPCRs), e.g. adenosine receptors (ARs) and P2Y receptors (P2YRs), respectively. We synthesize selective agents for use as pharmacological probes and potential therapeutic agents (e.g. A₃AR agonists for neuropathic pain). Detailed structural information derived from the X-ray crystallographic structures within these families enables the design of novel ligands, guides modification of known agonists and antagonists, and helps predict polypharmacology. Structures were recently reported for the P2Y₁₂ receptor (P2Y₁₂R), an anti-thrombotic target. Comparison of agonist-bound and antagonist-bound P2Y₁₂R indicates unprecedented structural plasticity in the outer portions of the transmembrane (TM) domains and the extracellular loops. Nonphosphate-containing ligands of the P2YRs, such as the selective P2Y₁₄R antagonist PPTN, are desired for bioavailability and increased stability. Also, A_2AAR structures are effectively applied to homology modeling of closely related A₁AR and A₃AR, which are not yet crystallized. Conformational constraint of normally flexible ribose with bicyclic analogues increased the ligand selectivity. Comparison of rigid A₃AR agonist congeners allows the exploration of interaction of specific regions of the nucleoside analogues with the target and off-target GPCRs, such as biogenic amine receptors. Molecular modeling predicts plasticity of the A₃AR at TM2 to accommodate highly rigidified ligands. Novel fluorescent derivatives of high affinity GPCR ligands are useful tool compounds for characterization of receptors and their oligomeric assemblies. Fluorescent probes are useful for characterization of GPCRs in living cells by flow cytometry and other methods. Thus, 3D knowledge of receptor binding and activation facilitates drug discovery.

Structural and energetic effects of A2A adenosine receptor mutations on agonist and antagonist binding

To predict structural and energetic effects of point mutations on ligand binding is of considerable interest in biochemistry and pharmacology. This is not only useful in connection with site-directed mutagenesis experiments, but could also allow interpretation and prediction of individual responses to drug treatment. For G-protein coupled receptors systematic mutagenesis has provided the major part of functional data as structural information until recently has been very limited. For the pharmacologically important A(2A) adenosine receptor, extensive site-directed mutagenesis data on agonist and antagonist binding is available and crystal structures of both types of complexes have been determined. Here, we employ a computational strategy, based on molecular dynamics free energy simulations, to rationalize and interpret available alanine-scanning experiments for both agonist and antagonist binding to this receptor. These computer simulations show excellent agreement with the experimental data and, most importantly, reveal the molecular details behind the observed effects which are often not immediately evident from the crystal structures. The work further provides a distinct validation of the computational strategy used to assess effects of point-mutations on ligand binding. It also highlights the importance of considering not only protein-ligand interactions but also those mediated by solvent water molecules, in ligand design projects.