Evaluation of molecular modeling of agonist binding in light of the crystallographic structure of an agonist-bound A₂A adenosine receptor

Binding Characterization of Agonists and Antagonists by MCCS: A Case Study from Adenosine A2A Receptor

Characterizing the structural basis of ligand recognition of adenosine A_2A receptor (AA_2AR) will facilitate its rational design and development of small molecules with high affinity and selectivity, as well as optimal therapeutic effects for pain, cancers, drug abuse disorders, etc. In the present work, we applied our reported algorithm, molecular complex characterizing system (MCCS), to characterize the binding features of AA_2AR based on its reported 3D structures of protein-ligand complexes. First, we compared the binding score to the reported experimental binding affinities of each compound. Then, we calculated an output example of residue energy contribution using MCCS and compared the results with data obtained from MM/GBSA. The consistency in results indicated that MCCS is a powerful, fast, and accurate method. Sequentially, using a receptor-ligand data set of 57 crystallized structures of AA_2ARs, we characterized the binding features of the binding pockets in AA_2AR, summarized the key residues that distinguish antagonist from agonist, produced heatmaps of residue energy contribution for clustering various statuses of AA_2ARs, explored the selectivity between AA_2AR and AA_1AR, etc. All the information provided new insights into the protein features of AA_2AR and will facilitate its rational drug design.

Identification of Novel Chemical Entities for Adenosine Receptor Type 2A Using Molecular Modeling Approaches

Adenosine Receptor Type 2A (A_2AAR) plays a role in important processes, such as anti-inflammatory ones. In this way, the present work aimed to search for compounds by pharmacophore-based virtual screening. The pharmacokinetic/toxicological profiles of the compounds, as well as a robust QSAR, predicted the binding modes via molecular docking. Finally, we used molecular dynamics to investigate the stability of interactions from ligand-A_2AAR. For the search for A_2AAR agonists, the UK-432097 and a set of 20 compounds available in the BindingDB database were studied. These compounds were used to generate pharmacophore models. Molecular properties were used for construction of the QSAR model by multiple linear regression for the prediction of biological activity. The best pharmacophore model was used by searching for commercial compounds in databases and the resulting compounds from the pharmacophore-based virtual screening were applied to the QSAR. Two compounds had promising activity due to their satisfactory pharmacokinetic/toxicological profiles and predictions via QSAR (Diverset 10002403 pEC₅₀ = 7.54407; ZINC04257548 pEC₅₀ = 7.38310). Moreover, they had satisfactory docking and molecular dynamics results compared to those obtained for Regadenoson (Lexiscan^®), used as the positive control. These compounds can be used in biological assays (in vitro and in vivo) in order to confirm the potential activity agonist to A_2AAR.

Accelerating the Throughput of Affinity Mass Spectrometry-Based Ligand Screening toward a G Protein-Coupled Receptor

Affinity mass spectrometry (MS) enables rapid screening of compound mixtures for ligands bound to a specific protein target, yet its current throughput is limited to individually assay pools of 400-2000 compounds. Typical affinity MS screens implemented in pharmaceutical industry laboratories identify putative ligands based on qualitative analysis of compound binding to the target whereas no quantitative information is acquired to discriminate high- and low-affinity ligands in the screening phase. Furthermore, these screens require purification of a stabilized form of the protein target, which poses a great challenge for membrane receptor targets. Here, we describe a new, potentially general affinity MS strategy that allows screening of 20,000 compounds in one pool for highly efficient ligand discovery toward a G protein-coupled receptor (GPCR) target. Quantitative measurement of compound binding to the receptor enables high-affinity ligand selection using both the purified receptor and receptor-embedded cell membranes. This high-throughput, label-free and quantitative affinity MS screen resulted in discovery of three new antagonists of the A2A adenosine receptor.

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.

Structure-Activity Relationships of the Sustained Effects of Adenosine A2A Receptor Agonists Driven by Slow Dissociation Kinetics

The duration of action of adenosine A_2A receptor (A2A) agonists is critical for their clinical efficacy, and we sought to better understand how this can be optimized. The in vitro temporal response profiles of a panel of A2A agonists were studied using cAMP assays in recombinantly (CHO) and endogenously (SH-SY5Y) expressing cells. Some agonists (e.g., 3cd; UK-432,097) but not others (e.g., 3ac; CGS-21680) demonstrated sustained wash-resistant agonism, where residual receptor activation continued after washout. The ability of an antagonist to reverse pre-established agonist responses was used as a surrogate read-out for agonist dissociation kinetics, and together with radioligand binding studies suggested a role for slow off-rate in driving sustained effects. One compound, 3ch, showed particularly marked sustained effects, with a reversal t_1/2 > 6 hours and close to maximal effects that remained for at least 5 hours after washing. Based on the structure-activity relationship of these compounds, we suggest that lipophilic N6 and bulky C2 substituents can promote stable and long-lived binding events leading to sustained agonist responses, although a high compound logD is not necessary. This provides new insight into the binding interactions of these ligands and we anticipate that this information could facilitate the rational design of novel long-acting A2A agonists with improved clinical efficacy.

Structurally Enabled Discovery of Adenosine A2A Receptor Antagonists

Over the past decade there has been a revolution in the field of G protein-coupled receptor (GPCR) structural biology. Many years of innovative research from different areas have come together to fuel this significant change in the fortunes of this field, which for many years was characterized by the paucity of high-resolution structures. The determination to succeed has been in part due to the recognized importance of these proteins as drug targets, and although the pharmaceutical industry has been focusing on these receptors, it can be justifiably argued and demonstrated that many of the approved and commercially successful GPCR drugs can be significantly improved to increase efficacy and/or reduce undesired side effects. In addition, many validated targets in this class remain to be drugged. It is widely recognized that application of structure-based drug design approaches can help medicinal chemists a long way toward discovering better drugs. The achievement of structural biologists in providing high-resolution insight is beginning to transform drug discovery efforts, and there are a number of GPCR drugs that have been discovered by use of structural information that are in clinical development. This review aims to highlight the key developments that have brought success to GPCR structure resolution efforts and exemplify the practical application of structural information for the discovery of adenosine A_2A receptor antagonists that have potential to treat multiple conditions.

The discovery of a selective and potent A2a agonist with extended lung retention

Although the anti-inflammatory role of the A2a receptor is well established, controversy remains with regard to the therapeutic value for A2a agonists in treatment of inflammatory lung diseases, also as a result of unwanted A2a-mediated cardiovascular effects. In this paper, we describe the discovery and characterization of a new, potent and selective A2a agonist (compound 2) with prolonged lung retention and limited systemic exposure following local administration. To support the lead optimization chemistry program with compound selection and profiling, multiple in vitro and in vivo assays were used, characterizing compound properties, pharmacodynamics (PD), and drug concentrations. Particularly, pharmacokinetic-PD modeling was applied to quantify the effects on the cardiovascular system, and an investigative toxicology study in rats was performed to explore potential myocardial toxicities. Compound 2, in comparison to a reference A2a agonist, UK-432,097, demonstrated higher solubility, lower lipophilicity, lower plasma protein binding, high rat lung retention (28% remaining after 24 h), and was efficacious in a lung inflammatory rat model following intratracheal dosing. Despite these properties, compound 2 did not provide a sufficient therapeutic index, that is, separation of local anti-inflammatory efficacy in the lung from systemic side effects in the cardiovascular system. The plasma concentration that resulted in induction of hypotension (half maximal effective concentration; EC50 0.5 nmol/L) correlated to the in vitro A2a potency (rIC50 0.6 nmol/L). Histopathological lesions in the heart were observed at a dose level which is threefold above the efficacious dose level in the inflammatory rat lung model. In conclusion, compound 2 is a highly potent and selective A2a agonist with significant lung retention after intratracheal administration. Despite its local anti-inflammatory efficacy in rat lung, small margins to the cardiovascular effects suggested limited therapeutic value of this compound for treatment of inflammatory lung disease by the inhaled route.

Design, synthesis, and biological evaluation of novel 2-((2-(4-(substituted)phenylpiperazin-1-yl)ethyl)amino)-5'-N-ethylcarboxamidoadenosines as potent and selective agonists of the A2A adenosine receptor

Stimulation of A2A adenosine receptors (AR) promotes anti-inflammatory responses in animal models of allergic rhinitis, asthma, chronic obstructive pulmonary disease, and rheumatic diseases. Herein we describe the results of a research program aimed at identifying potent and selective agonists of the A2AAR as potential anti-inflammatory agents. The recent crystallographic analysis of A2AAR agonists and antagonists in complex with the receptor provided key information on the structural determinants leading to receptor activation or blocking. In light of this, we designed a new series of 2-((4-aryl(alkyl)piperazin-1-yl)alkylamino)-5'-N-ethylcarboxamidoadenosines with high A2AAR affinity, activation potency and selectivity obtained by merging distinctive structural elements of known agonists and antagonists of the investigated target. Docking-based SAR optimization allowed us to identify compound 42 as one of the most potent and selective A2A agonist discovered so far (Ki hA2AAR = 4.8 nM, EC50 hA2AAR = 4.9 nM, Ki hA1AR > 10 000 nM, Ki hA3AR = 1487 nM, EC50 hA2BAR > 10 000 nM).

Molecular docking screening using agonist-bound GPCR structures: probing the A2A adenosine receptor

Crystal structures of G protein-coupled receptors (GPCRs) have recently revealed the molecular basis of ligand binding and activation, which has provided exciting opportunities for structure-based drug design. The A2A adenosine receptor (A2AAR) is a promising therapeutic target for cardiovascular diseases, but progress in this area is limited by the lack of novel agonist scaffolds. We carried out docking screens of 6.7 million commercially available molecules against active-like conformations of the A2AAR to investigate whether these structures could guide the discovery of agonists. Nine out of the 20 predicted agonists were confirmed to be A2AAR ligands, but none of these activated the ARs. The difficulties in discovering AR agonists using structure-based methods originated from limited atomic-level understanding of the activation mechanism and a chemical bias toward antagonists in the screened library. In particular, the composition of the screened library was found to strongly reduce the likelihood of identifying AR agonists, which reflected the high ligand complexity required for receptor activation. Extension of this analysis to other pharmaceutically relevant GPCRs suggested that library screening may not be suitable for targets requiring a complex receptor-ligand interaction network. Our results provide specific directions for the future development of novel A2AAR agonists and general strategies for structure-based drug discovery.

Conserved disulfide bond is not essential for the adenosine A2A receptor: Extracellular cysteines influence receptor distribution within the cell and ligand-binding recognition

G protein-coupled receptors (GPCRs) are integral membrane proteins involved in cellular signaling and constitute major drug targets. Despite their importance, the relationship between structure and function of these receptors is not well understood. In this study, the role of extracellular disulfide bonds on the trafficking and ligand-binding activity of the human A2A adenosine receptor was examined. To this end, cysteine-to-alanine mutations were conducted to replace individual and both cysteines in three disulfide bonds present in the first two extracellular loops. Although none of the disulfide bonds were essential for the formation of plasma membrane-localized active GPCR, loss of the disulfide bonds led to changes in the distribution of the receptor within the cell and changes in the ligand-binding affinity. These results indicate that in contrast to many class A GPCRs, the extracellular disulfide bonds of the A2A receptor are not essential, but can modulate the ligand-binding activity, by either changing the conformation of the extracellular loops or perturbing the interactions of the transmembrane domains.

New selective A2A agonists and A3 antagonists for human adenosine receptors: synthesis, biological activity and molecular docking studies

Synthesis and pharmacological characterization of a new series of adenosine derivatives on the four human adenosine receptors are reported.

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.

Structure-Based Design of Reactive Nucleosides for Site-Specific Modification of the A2A Adenosine Receptor

Adenosine receptors (ARs) are members of the G protein-coupled receptor (GPCR) superfamily and have shown much promise as therapeutic targets. We have used an agonist-bound A2AAR X-ray crystallographic structure to design a chemically reactive agonist for site-specific chemical modification of the receptor. To further explore and chemically engineer its binding cavity, a 2-nitrophenyl active ester was attached through an elongated chain at adenine C2 position. This general structure was designed for irreversible transfer of a terminal acyl group to a nucleophilic amino group on the A2AAR. Preincubation with several O-acyl derivatives prevented radioligand binding that was not regenerated upon extensive washing. In silico receptor docking suggested two lysine residues (second extracellular loop) as potential target sites for an O-acetyl derivative (MRS5854, 3a), and site-directed mutagenesis indicated that K153 but not K150 is essential. Similarly, a butyl azide for click reaction was incorporated in the active ester moiety (3b). These promising results indicate a stable, covalent modification of the receptor by several reactive adenosine derivatives, which could be chemical tools for future imaging, structural probing, and drug discovery. Thus, structure-based ligand design has guided the site-specific modification of a GPCR.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.

Crystal structures of the A2A adenosine receptor and their use in medicinal chemistry

New insights into drug design are derived from the X-ray crystallographic structures of G protein-coupled receptors (GPCRs), and the adenosine receptors (ARs) are at the forefront of this effort. The 3D knowledge of receptor binding and activation promises to enable drug discovery for GPCRs in general, and specifically for the ARs. The predictability of modeling based on the X-ray structures of the A_2AAR has been well demonstrated in the identification, design and modification of both known and novel AR agonists and antagonists. It is expected that structure-based design of drugs acting through ARs will provide new avenues to clinically useful agents.

Graph analysis of β2 adrenergic receptor structures: a "social network" of GPCR residues

Purpose

G protein-coupled receptors (GPCRs) are a superfamily of membrane proteins of vast pharmaceutical interest. Here, we describe a graph theory-based analysis of the structure of the β₂ adrenergic receptor (β₂ AR), a prototypical GPCR. In particular, we illustrate the network of direct and indirect interactions that link each amino acid residue to any other residue of the receptor.

Methods

Networks of interconnected amino acid residues in proteins are analogous to social networks of interconnected people. Hence, they can be studied through the same analysis tools typically employed to analyze social networks – or networks in general – to reveal patterns of connectivity, influential members, and dynamicity. We focused on the analysis of closeness-centrality, which is a measure of the overall connectivity distance of the member of a network to all other members.

Results

The residues endowed with the highest closeness-centrality are located in the middle of the seven transmembrane domains (TMs). In particular, they are mostly located in the middle of TM2, TM3, TM6 or TM7, while fewer of them are located in the middle of TM1, TM4 or TM5. At the cytosolic end of TM6, the centrality detected for the active structure is markedly lower than that detected for the corresponding residues in the inactive structures. Moreover, several residues acquire centrality when the structures are analyzed in the presence of ligands. Strikingly, there is little overlap between the residues that acquire centrality in the presence of the ligand in the blocker-bound structures and the agonist-bound structures.

Conclusions

Our results reflect the fact that the receptor resembles a bow tie, with a rather tight knot of closely interconnected residues and two ends that fan out in two opposite directions: one toward the extracellular space, which hosts the ligand binding cavity, and one toward the cytosol, which hosts the G protein binding cavity. Moreover, they underscore how interaction network is by the conformational rearrangements concomitant with the activation of the receptor and by the presence of agonists or blockers.

Electronic supplementary material

The online version of this article (doi:10.1186/2193-9616-1-16) contains supplementary material, which is available to authorized users.

Structure-based approaches to ligands for G-protein-coupled adenosine and P2Y receptors, from small molecules to nanoconjugates

Adenosine receptor (ARs) and P2Y receptors (P2YRs) that respond to extracellular nucleosides/nucleotides are associated with new directions for therapeutics. The X-ray structures of the A2AAR complexes with agonists and antagonists are examined in relationship to the G-protein-coupled receptor (GPCR) superfamily and applied to drug discovery. Much of the data on AR ligand structure from early SAR studies now are explainable from the A2AAR X-ray crystallography. The ligand-receptor interactions in related GPCR complexes can be identified by means of modeling approaches, e.g., molecular docking. Thus, molecular recognition in binding and activation processes has been studied effectively using homology modeling and applied to ligand design. Virtual screening has yielded new nonnucleoside AR antagonists, and existing ligands have been improved with knowledge of the receptor interactions. New agonists are being explored for central nervous system and peripheral therapeutics based on in vivo activity, such as chronic neuropathic pain. Ligands for receptors more distantly related to the X-ray template, i.e., P2YRs, have been introduced and are mainly used as pharmacological tools for elucidating the physiological role of extracellular nucleotides. Other ligand tools for drug discovery include fluorescent probes, radioactive probes, multivalent probes, and functionalized nanoparticles.

Conformation guides molecular efficacy in docking screens of activated β-2 adrenergic G protein coupled receptor

A prospective, large library virtual screen against an activated β2-adrenergic receptor (β2AR) structure returned potent agonists to the exclusion of inverse-agonists, providing the first complement to the previous virtual screening campaigns against inverse-agonist-bound G protein coupled receptor (GPCR) structures, which predicted only inverse-agonists. In addition, two hits recapitulated the signaling profile of the co-crystal ligand with respect to the G protein and arrestin mediated signaling. This functional fidelity has important implications in drug design, as the ability to predict ligands with predefined signaling properties is highly desirable. However, the agonist-bound state provides an uncertain template for modeling the activated conformation of other GPCRs, as a dopamine D2 receptor (DRD2) activated model templated on the activated β2AR structure returned few hits of only marginal potency.

Sphingosine-1-phosphate and its receptors: structure, signaling, and influence

The sphingosine-1-phosphate (S1P) receptor signaling system has biological and medical importance and is the first lipid G protein-coupled receptor (GPCR) structure to be solved to 2.8-Å resolution. S1P binds to five high-affinity GPCRs generating multiple downstream signals that play essential roles in vascular development and endothelial integrity, control of cardiac rhythm, and routine oral treatment of multiple sclerosis. Genetics, chemistry, and now structural biology have advanced this integrated biochemical system. The S1P receptors have a novel N-terminal fold that occludes access to the binding pocket from the extracellular environment as well as orthosteric and bitopic ligands with very different physicochemical properties. S1P receptors and metabolizing enzymes have been deleted, inducibly deleted, and knocked in as tagged or altered receptors in mice. An array of genetic models allows analysis of integrated receptor function in vivo. We can now directly understand causal relationships among protein expression, signal, and control points in physiology and pathology.

Probing GPCR structure: adenosine and P2Y nucleotide receptors

The adenosine receptors (ARs) provide an example of how to accurately predict ligand recognition, even prior to the availability of a crystallographic structure. Homology modeling has been used to gain structural insight, in conjunction with site-directed mutagenesis, and structure-activity relationships of small molecular ligands. Recent X-ray structures greatly improved the accuracy of knowledge of AR ligand recognition and furthermore characterized conformational changes induced by receptor activation. Now, homology modeling extends these structural insights to related GPCRs and suggests new ligand structures. This strategy is also being applied to the eight subtypes of P2Y receptors for extracellular nucleotides, which lack X-ray structures and are best modeled by homology to the CXCR4 (peptide) receptor. Neoceptors, as studied for three of the four AR subtypes, create a molecular complementarity between a mutant receptor and a chemically tailored agonist ligand to selectively enhance affinity, implying direct physical contact and thus validating docking hypotheses.

Stabilised G protein-coupled receptors in structure-based drug design: a case study with adenosine A2A receptor

The first example of structure-based drug design with stabilised GPCRs has enabled the identification of a preclinical candidate for the treatment of Parkinson's disease.

The GPCR Network: a large-scale collaboration to determine human GPCR structure and function

G protein-coupled receptors (GPCRs) are targeted by ∼30-40% of marketed drugs, and their key roles in normal physiology and in disease demonstrate that an understanding of their structure and function is valuable to researchers in both basic science and drug discovery. However, until recently, detailed structural information on this protein family was limited by challenges in X-ray crystallographic analysis of such membrane proteins. The GPCR Network was created in 2010 with the goal of structurally characterizing 15-25 representative human GPCRs within 5 years, based on an active outreach programme addressing an interdisciplinary community of scientists interested in GPCR structure, chemistry and biology. Here, we provide an overview of how this collaborative effort has enabled the structural determination and characterization of eight human GPCRs so far, and discuss some of the challenges that remain in gaining more detailed insights into structure-function relationships in this receptor superfamily.

Structural basis for allosteric regulation of GPCRs by sodium ions

Pharmacological responses of G protein-coupled receptors (GPCRs) can be fine-tuned by allosteric modulators. Structural studies of such effects have been limited due to the medium resolution of GPCR structures. We reengineered the human A(2A) adenosine receptor by replacing its third intracellular loop with apocytochrome b(562)RIL and solved the structure at 1.8 angstrom resolution. The high-resolution structure allowed us to identify 57 ordered water molecules inside the receptor comprising three major clusters. The central cluster harbors a putative sodium ion bound to the highly conserved aspartate residue Asp(2.50). Additionally, two cholesterols stabilize the conformation of helix VI, and one of 23 ordered lipids intercalates inside the ligand-binding pocket. These high-resolution details shed light on the potential role of structured water molecules, sodium ions, and lipids/cholesterol in GPCR stabilization and function.

New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors

Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A(2A) adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Structure-guided design of A(3) adenosine receptor-selective nucleosides: combination of 2-arylethynyl and bicyclo[3.1.0]hexane substitutions

(N)-Methanocarba adenosine 5'-methyluronamides containing known A(3) AR (adenosine receptor)-enhancing modifications, i.e., 2-(arylethynyl)adenine and N(6)-methyl or N(6)-(3-substituted-benzyl), were nanomolar full agonists of human (h) A(3)AR and highly selective (K(i) ∼0.6 nM, N(6)-methyl 2-(halophenylethynyl) analogues 13 and 14). Combined 2-arylethynyl-N(6)-3-chlorobenzyl substitutions preserved A(3)AR affinity/selectivity in the (N)-methanocarba series (e.g., 3,4-difluoro full agonist MRS5698 31, K(i) 3 nM, human and mouse A(3)) better than that for ribosides. Polyaromatic 2-ethynyl N(6)-3-chlorobenzyl analogues, such as potent linearly extended 2-p-biphenylethynyl MRS5679 34 (K(i) hA(3) 3.1 nM; A(1), A(2A), inactive) and fluorescent 1-pyrene adduct MRS5704 35 (K(i) hA(3) 68.3 nM), were conformationally rigid; receptor docking identified a large, mainly hydrophobic binding region. The vicinity of receptor-bound C2 groups was probed by homology modeling based on recent X-ray structure of an agonist-bound A(2A)AR, with a predicted helical rearrangement requiring an agonist-specific outward displacement of TM2 resembling opsin. Thus, the X-ray structure of related A(2A)AR is useful in guiding the design of new A(3)AR agonists.

Optimization of adenosine 5'-carboxamide derivatives as adenosine receptor agonists using structure-based ligand design and fragment screening

Structures of G protein-coupled receptors (GPCRs) have a proven utility in the discovery of new antagonists and inverse agonists modulating signaling of this important family of clinical targets. Applicability of active-state GPCR structures to virtual screening and rational optimization of agonists, however, remains to be assessed. In this study of adenosine 5' derivatives, we evaluated the performance of an agonist-bound A(2A) adenosine receptor (AR) structure in retrieval of known agonists and then employed the structure to screen for new fragments optimally fitting the corresponding subpocket. Biochemical and functional assays demonstrate high affinity of new derivatives that include polar heterocycles. The binding models also explain modest selectivity gain for some substituents toward the closely related A(1)AR subtype and the modified agonist efficacy of some of these ligands. The study suggests further applicability of in silico fragment screening to rational lead optimization in GPCRs.