Modeling the adenosine receptors: comparison of the binding domains of A2A agonists and antagonists

Bell-Evans model and steered molecular dynamics in uncovering the dissociation kinetics of ligands targeting G-protein-coupled receptors

Recently, academic and industrial scientific communities involved in kinetics-based drug development have become immensely interested in predicting the drug target residence time. Screening drug candidates in terms of their computationally predicted residence times, which is a measure of drug efficacy in vivo, and simultaneously assessing computational binding affinities are becoming inevitable. Non-equilibrium molecular simulation approaches are proven to be useful in this purpose. Here, we have implemented an optimized approach of combining the data derived from steered molecular dynamics simulations and the Bell-Evans model to predict the absolute residence times of the antagonist ZMA241385 and agonist NECA that target the A2A adenosine receptor of the G-protein-coupled receptor (GPCR) protein family. We have predicted the absolute ligand residence times on the timescale of seconds. However, our predictions were many folds shorter than those determined experimentally. Additionally, we calculated the thermodynamics of ligand binding in terms of ligand binding energies and the per-residue contribution of the receptor. Subsequently, binding pocket hotspot residues that would be important for further computational mutagenesis studies were identified. In the experiment, similar sets of residues were found to be in significant contact with both ligands under study. Our results build a strong foundation for further improvement of our approach by rationalizing the kinetics of ligand unbinding with the thermodynamics of ligand binding.

SINAPs: A Software Tool for Analysis and Visualization of Interaction Networks of Molecular Dynamics Simulations

As long as the structural study of molecular mechanisms requires multiple molecular dynamics reflecting contrasted bioactive states, the subsequent analysis of molecular interaction networks remains a bottleneck to be fairly treated and requires a user-friendly 3D view of key interactions. Structural Interaction Network Analysis Protocols (SINAPs) is a proprietary python tool developed to (i) quickly solve key interactions able to distinguish two protein states, either from two sets of molecular dynamics simulations or from two crystallographic structures, and (ii) render a user-friendly 3D view of these key interactions through a plugin of UCSF Chimera, one of the most popular open-source viewing software for biomolecular systems. Through two case studies, glucose transporter-1 (GLUT-1) and A2A adenosine receptor (A2AR), SINAPs easily pinpointed key interactions observed experimentally and relevant for their bioactivities. This very effective tool was thus applied to identify the amino acids involved in the molecular enzymatic mechanisms ruling the activation of an immunomodulator drug candidate, P28 glutathione-S-transferase (P28GST). SINAPs is freely available at https://github.com/ParImmune/SINAPs.

Current Therapeutic Molecules and Targets in Neurodegenerative Diseases Based on in silico Drug Design

Abstract: Background

As the number of elderly persons increases, neurodegenerative diseases are becoming ubiquitous. There is currently a great need for knowledge concerning management of old-age neurodegenerative diseases; the most important of which are: Alzheimer’s disease, Parkinson’s disease, Amyotrophic Lateral Sclerosis, and Huntington’s disease.

Objective

To summarize the potential of computationally predicted molecules and targets against neurodegenerative diseases.

Method

Review of literature published since 1997 against neurodegenerative diseases, utilizing as keywords: in silico, Alzheimer’s disease, Parkinson’s disease, Amyotrophic Lateral Sclerosis ALS, and Huntington’s disease was conducted.

Results and Conclusion

Due to the costs associated with experimentation and current ethical law, performing experiments directly on living organisms has become much more difficult. In this scenario, in silico techniques have been successful and have become powerful tools in the search to cure disease. Researchers use the Computer Aided Drug Design pipeline which: 1) generates 3-dimensional structures of target proteins through homology modeling 2) achieves stabilization through molecular dynamics simulation, and 3) exploits molecular docking through large compound libraries. Next generation sequencing is continually producing enormous amounts of raw sequence data while neuroimaging is producing a multitude of raw image data. To solve such pressing problems, these new tools and algorithms are required. This review elaborates precise in silico tools and techniques for drug targets, active molecules, and molecular docking studies, together with future prospects and challenges concerning possible breakthroughs in Alzheimer’s, Parkinson’s, Amyotrophic Lateral Sclerosis, and Huntington’s disease.

Comparison of the Human A2A Adenosine Receptor Recognition by Adenosine and Inosine: New Insight from Supervised Molecular Dynamics Simulations

Adenosine deaminase converts adenosine into inosine. In contrast to adenosine, relatively little attention has been paid to the physiological roles of inosine. Nevertheless, recent studies have demonstrated that inosine has neuroprotective, cardioprotective, immunomodulatory, and antidepressive effects. Inosine was recently shown to be a less potent agonist than adenosine at the A_2A adenosine receptor. To better depict the differences in the mechanisms of receptor recognition between adenosine and inosine, we carried out supervised molecular dynamics (SuMD) simulations, and the results are analyzed herein.

Supporting the Identification of Novel Fragment-Based Positive Allosteric Modulators Using a Supervised Molecular Dynamics Approach: A Retrospective Analysis Considering the Human A2A Adenosine Receptor as a Key Example

Structure-driven fragment-based (SDFB) approaches have provided efficient methods for the identification of novel drug candidates. This strategy has been largely applied in discovering several pharmacological ligand classes, including enzyme inhibitors, receptor antagonists and, more recently, also allosteric (positive and negative) modulators. Recently, Siegal and collaborators reported an interesting study, performed on a detergent-solubilized StaR adenosine A_2A receptor, describing the existence of both fragment-like negative allosteric modulators (NAMs), and fragment-like positive allosteric modulators (PAMs). From this retrospective study, our results suggest that Supervised Molecular Dynamics (SuMD) simulations can support, on a reasonable time scale, the identification of fragment-like PAMs following their receptor recognition pathways and characterizing the possible allosteric binding sites.

Exploring the Role of N6-Substituents in Potent Dual Acting 5'-C-Ethyltetrazolyladenosine Derivatives: Synthesis, Binding, Functional Assays, and Antinociceptive Effects in Mice ∇

Structural determinants of affinity of N⁶-substituted-5'-C-(ethyltetrazol-2-yl)adenosine and 2-chloroadenosine derivatives at adenosine receptor (AR) subtypes were studied with binding and molecular modeling. Small N⁶-cycloalkyl and 3-halobenzyl groups furnished potent dual acting A₁AR agonists and A₃AR antagonists. 4 was the most potent dual acting human (h) A₁AR agonist (K_i = 0.45 nM) and A₃AR antagonist (K_i = 0.31 nM) and highly selective versus A_2A; 11 and 26 were most potent at both h and rat (r) A₃AR. All N⁶-substituted-5'-C-(ethyltetrazol-2-yl)adenosine derivatives proved to be antagonists at hA₃AR but agonists at the rA₃AR. Analgesia of 11, 22, and 26 was evaluated in the mouse formalin test (A₃AR antagonist blocked and A₃AR agonist strongly potentiated). N⁶-Methyl-5'-C-(ethyltetrazol-2-yl)adenosine (22) was most potent, inhibiting both phases, as observed combining A₁AR and A₃AR agonists. This study demonstrated for the first time the advantages of a single molecule activating two AR pathways both leading to benefit in this acute pain model.

Molecular-Docking-Based Drug Design and Discovery

Currently 30-50% of drug targets are G Protein-Coupled Receptors (GPCRs). However, the clinical useful drugs for targeting GPCR have been limited by the lack of subtype selectivity or efficacy, leading to undesirable side effects. To develop subtype-selective GPCR ligands with desired molecular properties, better understanding is needed of the pharmacophore elements and of the binding mechanism required for subtype selectivity. To illustrate these issues, we describe here three successful applications to understand the binding mechanism associated with subtype selectivity: 5-HT2B (5-Hydroxytryptamine, 5-HT) serotonin receptor (HT2BR), H3 histamine receptor (H3HR) and A3 adenosine receptor (A3AR). The understanding of structure-function relationships among individual types and subtypes of GPCRs gained from such computational predictions combined with experimental validation and testing is expected the development of new highly selective and effective ligands to address such diseases while minimizing side-effects.

N6-(2-Hydroxyethyl)-Adenosine Exhibits Insecticidal Activity against Plutella xylostella via Adenosine Receptors

The diamondback moth, Plutella xylostella, is one of the most important pests of cruciferous crops. We have earlier shown that N6-(2-hydroxyethyl)-adenosine (HEA) exhibits insecticidal activity against P. xylostella. In the present study we investigated the possible mechanism of insecticidal action of HEA on P. xylostella. HEA is a derivative of adenosine, therefore, we speculated whether it acts via P. xylostella adenosine receptor (PxAdoR). We used RNAi approach to silence PxAdoR gene and used antagonist of denosine receptor (AdoR) to study the insecticidal effect of HEA. We cloned the whole sequence of PxAdoR gene. A BLAST search using NCBI protein database showed a 61% identity with the Drosophila adenosine receptor (DmAdoR) and a 32-35% identity with human AdoR. Though the amino acids sequence of PxAdoR was different compared to other adenosine receptors, most of the amino acids that are known to be important for adenosine receptor ligand binding and signaling were present. However, only 30% binding sites key residues was similar between PxAdoR and A1R. HEA, at a dose of 1 mg/mL, was found to be lethal to the second-instar larvae of P. xylostella, and a significant reduction of mortality and growth inhibition ratio were obtained when HEA was administered to the larvae along with PxAdoR-dsRNA or antagonist of AdoR (SCH58261) for 36, 48, or 60 h. Especially at 48 h, the rate of growth inhibition of the PxAdoR knockdown group was 3.5-fold less than that of the HEA group, and the corrected mortality of SCH58261 group was reduced almost 2-fold compared with the HEA group. Our findings show that HEA may exert its insecticidal activity against P. xylostella larvae via acting on PxAdoR.

Similarities and differences in affinity and binding modes of tricyclic pyrimido- and pyrazinoxanthines at human and rat adenosine receptors

A new series of 32 pyrimido- and 5 tetrahydropyrazino[2,1-f]purinediones was obtained and evaluated for their adenosine receptors (ARs) affinities. The 1,3-dibutyl derivative of 9-(4-(2-(dimethylamino)ethoxy)phenyl)-6,7,8,9-tetrahydropyrimido[1,2-f]purine-2,4(1H,3H)-dione was found to be the most potent A1 AR antagonist of the present series, showing selectivity over the other AR subtypes. The structure-activity for the obtained purinediones was established. Docking experiments of the investigated library to homology models of the human and rat A1 and A2A ARs allowed to compare the expected binding modes for selected compounds. The detailed analysis of binding cavities within individual AR subtypes indicated small but significant structural variations that may underlie the observed differences in binding affinities of purinediones at particular subtypes and species.

Novel homologated-apio adenosine derivatives as A3 adenosine receptor agonists: design, synthesis and molecular docking studies

A series of homologated-apio adenosine derivatives including homologated-apio IB-MECA and Cl-IB-MECA have been designed and synthesized successfully. The molecular modelling and docking studies of the compounds have been explored as A₃AR agonists.

Mutation-Guided Unbiased Modeling of the Fat Sensor GPR119 for High-Yield Agonist Screening

Recent benchmark studies have demonstrated the difficulties in obtaining accurate predictions of ligand binding conformations to comparative models of G-protein-coupled receptors. We have developed a data-driven optimization protocol, which integrates mutational data and structural information from multiple X-ray receptor structures in combination with a fully flexible ligand docking protocol to determine the binding conformation of AR231453, a small-molecule agonist, in the GPR119 receptor. Resulting models converge to one conformation that explains the majority of data from mutation studies and is consistent with the structure-activity relationship for a large number of AR231453 analogs. Another key property of the refined models is their success in separating active ligands from decoys in a large-scale virtual screening. These results demonstrate that mutation-guided receptor modeling can provide predictions of practical value for describing receptor-ligand interactions and drug discovery.

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).

Characteristic molecular vibrations of adenosine receptor ligands

Although the regulation of membrane receptor activation is known to be crucial for molecular signal transduction, the molecular mechanism underlying receptor activation is not fully elucidated. Here we study the physicochemical nature of membrane receptor behavior by investigating the characteristic molecular vibrations of receptor ligands using computational chemistry and informatics methods. By using information gain, t-tests, and support vector machines, we have identified highly informative features of adenosine receptor (AdoR) ligand and corresponding functional amino acid residues such as Asn (6.55) of AdoR that has informative significance and is indispensable for ligand recognition of AdoRs. These findings may provide new perspectives and insights into the fundamental mechanism of class A G protein-coupled receptor activation.

Free energy calculations of A2Aadenosine receptor mutation effects on agonist binding

A general computational scheme to evaluate the effects of single point mutations on ligand binding is reported.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Perturbation of fluid dynamics properties of water molecules during G protein-coupled receptor-ligand recognition: the human A2A adenosine receptor as a key study

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.

Allosteric sodium in class A GPCR signaling

Despite their functional and structural diversity, G-protein-coupled receptors (GPCRs) share a common mechanism of signal transduction via conformational changes in the seven-transmembrane (7TM) helical domain. New major insights into this mechanism come from the recent crystallographic discoveries of a partially hydrated sodium ion that is specifically bound in the middle of the 7TM bundle of multiple class A GPCRs. This review discusses the remarkable structural conservation and distinct features of the Na(+) pocket in this most populous GPCR class, as well as the conformational collapse of the pocket upon receptor activation. New insights help to explain allosteric effects of sodium on GPCR agonist binding and activation, and sodium's role as a potential co-factor in class A GPCR function.

Synthesis and anti-renal fibrosis activity of conformationally locked truncated 2-hexynyl-N(6)-substituted-(N)-methanocarba-nucleosides as A3 adenosine receptor antagonists and partial agonists

Truncated N⁶-substituted-(N)-methanocarba-adenosine derivatives with 2-hexynyl substitution were synthesized to examine parallels with corresponding 4′-thioadenosines. Hydrophobic N⁶ and/or C2 substituents were tolerated in A₃AR binding, but only an unsubstituted 6-amino group with a C2-hexynyl group promoted high hA_2AAR affinity. A small hydrophobic alkyl (4b and 4c) or N⁶-cycloalkyl group (4d) showed excellent binding affinity at the hA₃AR and was better than an unsubstituted free amino group (4a). A₃AR affinities of 3-halobenzylamine derivatives 4f–4i did not differ significantly, with K_i values of 7.8–16.0 nM. N⁶-Methyl derivative 4b (K_i = 4.9 nM) was a highly selective, low efficacy partial A₃AR agonist. All compounds were screened for renoprotective effects in human TGF-β1-stimulated mProx tubular cells, a kidney fibrosis model. Most compounds strongly inhibited TGF-β1-induced collagen I upregulation, and their A₃AR binding affinities were proportional to antifibrotic effects; 4b was most potent (IC₅₀ = 0.83 μM), indicating its potential as a good therapeutic candidate for treating renal fibrosis.

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.

Mapping the intramolecular signal transduction of G-protein coupled receptors

G-protein coupled receptors (GPCRs), a major gatekeeper of extracellular signals on plasma membrane, are unarguably one of the most important therapeutic targets. Given the recent discoveries of allosteric modulations, an allosteric wiring diagram of intramolecular signal transductions would be of great use to glean the mechanism of receptor regulation. Here, by evaluating betweenness centrality (CB ) of each residue, we calculate maps of information flow in GPCRs and identify key residues for signal transductions and their pathways. Compared with preexisting approaches, the allosteric hotspots that our CB -based analysis detects for A2 A adenosine receptor (A2 A AR) and bovine rhodopsin are better correlated with biochemical data. In particular, our analysis outperforms other methods in locating the rotameric microswitches, which are generally deemed critical for mediating orthosteric signaling in class A GPCRs. For A2 A AR, the inter-residue cross-correlation map, calculated using equilibrium structural ensemble from molecular dynamics simulations, reveals that strong signals of long-range transmembrane communications exist only in the agonist-bound state. A seemingly subtle variation in structure, found in different GPCR subtypes or imparted by agonist bindings or a point mutation at an allosteric site, can lead to a drastic difference in the map of signaling pathways and protein activity. The signaling map of GPCRs provides valuable insights into allosteric modulations as well as reliable identifications of orthosteric signaling pathways.

Conformationally constrained ortho-anilino diaryl ureas: discovery of 1-(2-(1'-neopentylspiro[indoline-3,4'-piperidine]-1-yl)phenyl)-3-(4-(trifluoromethoxy)phenyl)urea, a potent, selective, and bioavailable P2Y1 antagonist

Preclinical antithrombotic efficacy and bleeding models have demonstrated that P2Y1 antagonists are efficacious as antiplatelet agents and may offer a safety advantage over P2Y12 antagonists in terms of reduced bleeding liabilities. In this article, we describe the structural modification of the tert-butyl phenoxy portion of lead compound 1 and the subsequent discovery of a novel series of conformationally constrained ortho-anilino diaryl ureas. In particular, spiropiperidine indoline-substituted diaryl ureas are described as potent, orally bioavailable small-molecule P2Y1 antagonists with improved activity in functional assays and improved oral bioavailability in rats. Homology modeling and rat PK/PD studies on benchmark compound 3l will also be presented. Compound 3l was our first P2Y1 antagonist to demonstrate a robust oral antithrombotic effect with mild bleeding liability in the rat thrombosis and hemostasis models.

Synthesis, biological activity and molecular modelling studies of tricyclic alkylimidazo-, pyrimido- and diazepinopurinediones

Syntheses and biological activities of imidazo-, pyrimido- and diazepino[2,1-f]purinediones containing N-alkyl substituents (with straight, branched or unsaturated chains) are described. Tricyclic derivatives were synthesized by the cyclization of 8-bromo-substituted 7-(2-bromoethyl)-, 7-(3-chloropropyl)- or 7-(4-bromobutyl)-theophylline with primary amines under various conditions. Compound 22 with an ethenyl substituent was synthesized by dehydrohalogenation of 9-(2-bromoethyl)-1,3-dimethyltetrahydropyrimido[2,1-f]purinedione. The obtained derivatives (5-35) were initially evaluated for their affinity at rat A1 and A2A adenosine receptors (AR), showing moderate affinity for both adenosine receptor subtypes. The best ligands were diazepinopurinedione 28 (K i = 0.28 μM) with fivefold A2A selectivity and the non-selective A1/A2A AR ligand pyrimidopurinedione 35 (K i A1 = 0.28 μM and K i A2A = 0.30 μM). The compounds were also evaluated for their affinity at human A1, A2A, A2B and A3 ARs. All of the obtained compounds were docked to the A2A AR X-ray structure in complex with the xanthine-based, potent adenosine receptor antagonist-XAC. The likely interactions of imidazo-, pyrimido- and diazepino[2,1-f]purinediones with the residues forming the A2A binding pocket were discussed. Furthermore, the new compounds were tested in vivo as anticonvulsants in maximal electroshock, subcutaneous pentylenetetrazole (ScMet) and TOX tests in mice (i.p.). Pyrimidopurinediones showed anticonvulsant activity mainly in the ScMet test. The best derivative was compound 11, showing 100 % protection at a dose of 100 mg/kg without symptoms of neurotoxicity. Compounds 6, 7, 8 and 14 with short substituents showed neurotoxicity and caused death. In rat tests (p.o.), 9 was characterized by a high protection index (>13.3). AR affinity did not apparently correlate with the antiepileptic potency of the compounds.

Anticonvulsant activity of B2, an adenosine analog, on chemical convulsant-induced seizures

Epilepsy is a chronic neurological disorder characterized by recurrent seizures. However, approximately one-third of epilepsy patients still suffer from uncontrolled seizures. Effective treatments for epilepsy are yet to be developed. N⁶-(3-methoxyl-4-hydroxybenzyl) adenine riboside (B2) is a N⁶-substitued adenosine analog. Here we describe an investigation of the effects and mechanisms of B2 on chemical convulsant-induced seizures. Seizures were induced in mice by administration of 4-aminopyridine (4-AP), pentylenetetrazol (PTZ), picrotoxin, kainite acid (KA), or strychnine. B2 has a dose-related anticonvulsant effect in these chemical-induced seizure models. The protective effects of B2 include increased latency of seizure onset, decreased seizure occurrence, shorter seizure duration and reduced mortality rate. Radioligand binding and cAMP accumulation assays indicated that B2 might be a functional ligand for both adenosine A₁ and A_2A receptors. Furthermore, DPCPX, a selective A₁ receptor antagonist, but not SCH58261, a selective A_2A receptor antagonist, blocked the anticonvulsant effect of B2 on PTZ-induced seizure. c-Fos is a cellular marker for neuronal activity. Immunohistochemical and western blot analyses indicated that B2 significantly reversed PTZ-induced c-Fos expression in the hippocampus. Together, these results indicate that B2 has significant anticonvulsant effects. The anticonvulsant effects of B2 may be attributed to adenosine A₁ receptor activation and reduced neuronal excitability in the hippocampus. These observations also support that the use of adenosine receptor agonist may be a promising approach for the treatment of epilepsy.

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

The application of both structure- and ligand-based design approaches represents to date one of the most useful strategies in the discovery of new drug candidates. In the present paper, we investigated how the application of docking-driven conformational analysis can improve the predictive ability of 3D-QSAR statistical models. With the use of the crystallographic structure in complex with the high affinity antagonist ZM 241385 (4-(2-[7-amino-2-(2-furyl)[1,2,4]-triazolo[2,3-a][1,3,5]triazin-5-ylamino]ethyl)phenol), we revisited a general pharmacophore hypothesis for the human A(2A) adenosine receptor of a set of 751 known antagonists, by applying an integrated ligand- and structure-based approach. Our novel pharmacophore hypothesis has been validated by using an external test set of 29 newly synthesized human adenosine receptor antagonists.

Insight into the binding mode and the structural features of the pyrimidine derivatives as human A2A adenosine receptor antagonists

The interaction of 278 monocyclic and bicyclic pyrimidine derivatives with human A2A adenosine receptor (AR) was investigated by employing molecular dynamics, thermodynamic analysis and three-dimensional quantitative structure-activity relationship (3D-QSAR) approaches. The binding analysis reveals that the pyrimidine derivatives are anchored in TM2, 3, 5, 6 and 7 of A2A AR by the aromatic stacking and hydrogen bonding interactions. The key residues involving Phe168, Glu169, and Asn253 stabilize the monocyclic and bicyclic cores of inhibitors. The thermodynamic analysis by molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) approach also confirms the reasonableness of the binding modes. In addition, the ligand-/receptor-based comparative molecular similarity indices analysis (CoMSIA) models of high statistical significance were generated and the resulting contour maps correlate well with the structural features of the antagonists essential for high A2A AR affinity. A minor/bulky group with negative charge at C2/C6 of pyrimidine ring respectively enhances the activity for all these pyrimidine derivatives. Particularly, the higher electron density of the ring in the bicyclic derivatives, the more potent the antagonists. The obatined results might be helpful in rational design of novel candidate of A2A adenosine receptor antagonist for treatment of Parkinson's disease.

Discovery of 2-(phenoxypyridine)-3-phenylureas as small molecule P2Y1 antagonists

Two distinct G protein-coupled purinergic receptors, P2Y1 and P2Y12, mediate ADP-driven platelet activation. The clinical effectiveness of P2Y12 blockade is well established. Recent preclinical data suggest that P2Y1 and P2Y12 inhibition provide equivalent antithrombotic efficacy, while targeting P2Y1 has the potential for reduced bleeding liability. In this account, the discovery of a 2-(phenoxypyridine)-3-phenylurea chemotype that inhibited ADP-mediated platelet aggregation in human blood samples is described. Optimization of this series led to the identification of compound 16, 1-(2-(2-tert-butylphenoxy)pyridin-3-yl)-3-4-(trifluoromethoxy)phenylurea, which demonstrated a 68 ± 7% thrombus weight reduction in an established rat arterial thrombosis model (10 mg/kg plus 10 mg/kg/h) while only prolonging cuticle and mesenteric bleeding times by 3.3- and 3.1-fold, respectively, in provoked rat bleeding time models. These results suggest that a P2Y1 antagonist could potentially provide a safe and efficacious antithrombotic profile.

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.

Functional efficacy of adenosine A₂A receptor agonists is positively correlated to their receptor residence time

BACKGROUND AND PURPOSE The adenosine A(2A) receptor belongs to the superfamily of GPCRs and is a promising therapeutic target. Traditionally, the discovery of novel agents for the A(2A) receptor has been guided by their affinity for the receptor. This parameter is determined under equilibrium conditions, largely ignoring the kinetic aspects of the ligand-receptor interaction. The aim of this study was to assess the binding kinetics of A(2A) receptor agonists and explore a possible relationship with their functional efficacy. EXPERIMENTAL APPROACH We set up, validated and optimized a kinetic radioligand binding assay (a so-called competition association assay) at the A(2A) receptor from which the binding kinetics of unlabelled ligands were determined. Subsequently, functional efficacies of A(2A) receptor agonists were determined in two different assays: a novel label-free impedance-based assay and a more traditional cAMP determination. KEY RESULTS A simplified competition association assay yielded an accurate determination of the association and dissociation rates of unlabelled A(2A) receptor ligands at their receptor. A correlation was observed between the receptor residence time of A(2A) receptor agonists and their intrinsic efficacies in both functional assays. The affinity of A(2A) receptor agonists was not correlated to their functional efficacy. CONCLUSIONS AND IMPLICATIONS This study indicates that the molecular basis of different agonist efficacies at the A(2A) receptor lies within their different residence times at this receptor.

NHBA isolated from Gastrodia elata exerts sedative and hypnotic effects in sodium pentobarbital-treated mice

The rhizomes of Gastrodia elata have been used for the treatment of insomnia in oriental countries. N⁶-(4-hydroxybenzyl) adenine riboside (NHBA) was originally isolated from G. elata. For the first time we report a detailed study on the effects and mechanisms of NHBA on its sedative and hypnotic activity. Adenosine, an endogenous sleep factor, regulates sleep-wake cycle via interacting with adenosine A₁/A(2A) receptors. Using radioligand binding studies and cAMP accumulation assays, our results show that NHBA may be a functional ligand for the adenosine A₁ and A(2A) receptors. NHBA significantly decreases spontaneous locomotor activity and potentiates the hypnotic effect of sodium pentobarbital in mice. Sleep architecture analyses reveal that NHBA significantly decreases wakefulness time and increases NREM sleep times. However, NHBA does not affect the amount of REM sleep. Pretreatment with the adenosine A₁ receptor antagonist DPCPX or the A(2A) receptor antagonist SCH 58261 significantly reverses the increase in sleeping time induced by NHBA in sodium pentobarbital treated mice. Immunohistochemical studies show that NHBA increases c-Fos expression in GABAergic neurons of the ventrolateral preoptic area (VLPO), which suggests that NHBA activates the sleep center in the anterior hypothalamus. Altogether, these results indicate that NHBA produces significant sedative and hypnotic effects. Such effects might be mediated by the activation of adenosine A₁/A(2A) receptors and stimulation of the sleep center VLPO.

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.

Characterization of adenosine receptor in its native environment: insights from molecular dynamics simulations of palmitoylated/glycosylated, membrane-integrated human A(2B) adenosine receptor

Selective A(2B) receptor antagonists and agonists may play a role in important pathologies such as gastrointestinal, neurological (i.e., Alzheimer disease and dementia) and hypersensitive disorders (i.e., asthma), diabetes, atherosclerosis, restenosis and cancer. Hence, it is regarded as a good target for the development of clinically useful agents. In this study, the effects of lipid bilayer, N-acetylglucosamine and S-palmitoyl on the dynamic behavior of A(2B)AR model is explored. Homology modeling, molecular docking and molecular dynamics simulations were performed to explore structural features of A(2B)AR in the presence of lipid bilayer. Twenty ns MD simulation was performed on the constructed model inserted in a hydrated lipid bilayer to examine stability of the best model. OSIP339391 as the most potent antagonist was docked in the active site of the model. Another MD simulation was performed on the ligand-protein complex to explore effects of the bilayer on this complex. A similar procedure was performed for the modified protein with N-acetylglucosamine and S-palmitoyl moieties in its structure. Phe173 and Glu174 located in EL2 were determined to be involved in ligand-receptor interactions through π-π stacking and hydrogen bonding. Asn254 was crucial to form hydrogen-bonding. The reliability of the model was assessed through docking using both commercial and synthetic antagonists and an r(2) of 0.70 was achieved. Our results show that molecular dynamics simulations of palmitoylated/glycosylated, membrane-integrated human A(2B)AR in its native environment is a possible approach and this model can be used for designing potent and selective A(2B)AR antagonists.

Characterization of the Drosophila adenosine receptor: the effect of adenosine analogs on cAMP signaling in Drosophila cells and their utility for in vivo experiments

Adenosine receptors (AR) belonging to the G protein-coupled receptor family influence a wide range of physiological processes. Recent elucidation of the structure of human A2AR revealed the conserved amino acids necessary for contact with the Ado moiety. However, the selectivity of Ado analogs for AR subtypes is still not well understood. We have shown previously that the Drosophila adenosine receptor (DmAdoR) evokes an increase in cAMP and calcium concentration in heterologous cells. In this study, we have characterized the second-messenger stimulation by endogenous DmAdoR in a Drosophila neuroblast cell line and examined a number of Ado analogs for their ability to interact with DmAdoR. We show that Ado can stimulate cAMP but not calcium levels in Drosophila cells. We found one full and four partial DmAdoR agonists, as well as four antagonists. The employment of the full agonist, 2-chloroadenosine, in flies mimicked in vivo the phenotype of DmAdoR over-expression, whereas the antagonist, SCH58261, rescued the flies from the lethality caused by DmAdoR over-expression. Differences in pharmacological effect of the tested analogs between DmAdoR and human A2AR can be partially explained by the dissimilarity of specific key amino acid residues disclosed by the alignment of these receptors.

Discovery of 1,2,4-triazine derivatives as adenosine A(2A) antagonists using structure based drug design

Potent, ligand efficient, selective, and orally efficacious 1,2,4-triazine derivatives have been identified using structure based drug design approaches as antagonists of the adenosine A_2A receptor. The X-ray crystal structures of compounds 4e and 4g bound to the GPCR illustrate that the molecules bind deeply inside the orthosteric binding cavity. In vivo pharmacokinetic and efficacy data for compound 4k are presented, demonstrating the potential of this series of compounds for the treatment of Parkinson’s disease.

Structure of the adenosine A(2A) receptor in complex with ZM241385 and the xanthines XAC and caffeine

Methylxanthines, including caffeine and theophylline, are among the most widely consumed stimulant drugs in the world. These effects are mediated primarily via blockade of adenosine receptors. Xanthine analogs with improved properties have been developed as potential treatments for diseases such as Parkinson's disease. Here we report the structures of a thermostabilized adenosine A(2A) receptor in complex with the xanthines xanthine amine congener and caffeine, as well as the A(2A) selective inverse agonist ZM241385. The receptor is crystallized in the inactive state conformation as defined by the presence of a salt bridge known as the ionic lock. The complete third intracellular loop, responsible for G protein coupling, is visible consisting of extended helices 5 and 6. The structures provide new insight into the features that define the ligand binding pocket of the adenosine receptor for ligands of diverse chemotypes as well as the cytoplasmic regions that interact with signal transduction proteins.

A novel nonribose agonist, LUF5834, engages residues that are distinct from those of adenosine-like ligands to activate the adenosine A(2a) receptor

The recent publication of both the antagonist- and agonist-bound structures of the adenosine A(2A) receptor have revealed much about how a ligand may bind to a receptor and cause the conformational changes associated with agonist-mediated activation. In particular, the agonist-bound structure revealed key interactions between the ribose group of adenosine-derived agonists and amino acids in the receptor binding pocket that lead to receptor activation. However, agonists without a ribose group also exist, and we wondered whether such compounds occupy the same agonist binding site. Therefore we used a mutagenesis approach in this study to investigate the mode of binding of 2-amino-4-(4-hydroxyphenyl)- 6-(1H-imidazol-2-ylmethylsulfanyl)pyridine-3,5-dicarbonitrile (LUF5834), a potent partial agonist without a ribose moiety, compared with the adenosine-derived reference agonist 2-[p-(2-carboxyethyl)phenyl-ethylamino]-5'-N-ethylcarboxamidoadenosine (CGS21680). Mutation of the orthosteric residue Phe168 to alanine abrogated the function of both agonists. However, mutation to alanine of residues Thr88 and Ser277 shown by the crystal structures to interact with the ribose group of adenosine-like ligands had no effect on the potency of LUF5834. Furthermore, alanine mutation of Asn253, which makes a hydrogen-bonding interaction with the exocyclic nitrogen of the adenine ring, had minimal effect on LUF5834 affinity but removed agonist activity of this ligand. Mutation of other residues, such as the highly conserved Trp246 or Glu13, had significant deleterious effects on the function of CGS21680 but little effect on LUF5834. In summary, our findings suggest that this class of agonist interacts with distinct residues to activate the receptor compared with classic adenosine derived agonists.

Status of GPCR modeling and docking as reflected by community-wide GPCR Dock 2010 assessment

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe.

The A3 adenosine receptor as multifaceted therapeutic target: pharmacology, medicinal chemistry, and in silico approaches

Adenosine is an ubiquitous local modulator that regulates various physiological and pathological functions by stimulating four membrane receptors, namely A(1), A(2A), A(2B), and A(3). Among these G protein-coupled receptors, the A(3) subtype is found mainly in the lung, liver, heart, eyes, and brain in our body. It has been associated with cerebroprotection and cardioprotection, as well as modulation of cellular growth upon its selective activation. On the other hand, its inhibition by selective antagonists has been reported to be potentially useful in the treatment of pathological conditions including glaucoma, inflammatory diseases, and cancer. In this review, we focused on the pharmacology and the therapeutic implications of the human (h)A(3) adenosine receptor (AR), together with an overview on the progress of hA(3) AR agonists, antagonists, allosteric modulators, and radioligands, as well as on the recent advances pertaining to the computational approaches (e.g., quantitative structure-activity relationships, homology modeling, molecular docking, and molecular dynamics simulations) applied to the modeling of hA(3) AR and drug design. 

ss-TEA: Entropy based identification of receptor specific ligand binding residues from a multiple sequence alignment of class A GPCRs

BACKGROUND

G-protein coupled receptors (GPCRs) are involved in many different physiological processes and their function can be modulated by small molecules which bind in the transmembrane (TM) domain. Because of their structural and sequence conservation, the TM domains are often used in bioinformatics approaches to first create a multiple sequence alignment (MSA) and subsequently identify ligand binding positions. So far methods have been developed to predict the common ligand binding residue positions for class A GPCRs.

RESULTS

Here we present 1) ss-TEA, a method to identify specific ligand binding residue positions for any receptor, predicated on high quality sequence information. 2) The largest MSA of class A non olfactory GPCRs in the public domain consisting of 13324 sequences covering most of the species homologues of the human set of GPCRs. A set of ligand binding residue positions extracted from literature of 10 different receptors shows that our method has the best ligand binding residue prediction for 9 of these 10 receptors compared to another state-of-the-art method.

CONCLUSIONS

The combination of the large multi species alignment and the newly introduced residue selection method ss-TEA can be used to rapidly identify subfamily specific ligand binding residues. This approach can aid the design of site directed mutagenesis experiments, explain receptor function and improve modelling. The method is also available online via GPCRDB at http://www.gpcr.org/7tm/.

Biophysical mapping of the adenosine A2A receptor

A new approach to generating information on ligand receptor interactions within the binding pocket of G protein-coupled receptors has been developed, called Biophysical Mapping (BPM). Starting from a stabilized receptor (StaR), minimally engineered for thermostability, additional single mutations are then added at positions that could be involved in small molecule interactions. The StaR and a panel of binding site mutants are captured onto Biacore chips to enable characterization of the binding of small molecule ligands using surface plasmon resonance (SPR) measurement. A matrix of binding data for a set of ligands versus each active site mutation is then generated, providing specific affinity and kinetic information (K(D), k(on), and k(off)) of receptor-ligand interactions. This data set, in combination with molecular modeling and docking, is used to map the small molecule binding site for each class of compounds. Taken together, the many constraints provided by these data identify key protein-ligand interactions and allow the shape of the site to be refined to produce a high quality three-dimensional picture of ligand binding, thereby facilitating structure based drug design. Results of biophysical mapping of the adenosine A(2A) receptor are presented.

Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation

Adenosine receptors and β-adrenoceptors are G-protein-coupled receptors (GPCRs) that activate intracellular G proteins on binding the agonists adenosine or noradrenaline, respectively. GPCRs have similar structures consisting of seven transmembrane helices that contain well-conserved sequence motifs, indicating that they are probably activated by a common mechanism. Recent structures of β-adrenoceptors highlight residues in transmembrane region 5 that initially bind specifically to agonists rather than to antagonists, indicating that these residues have an important role in agonist-induced activation of receptors. Here we present two crystal structures of the thermostabilized human adenosine A(2A) receptor (A(2A)R-GL31) bound to its endogenous agonist adenosine and the synthetic agonist NECA. The structures represent an intermediate conformation between the inactive and active states, because they share all the features of GPCRs that are thought to be in a fully activated state, except that the cytoplasmic end of transmembrane helix 6 partially occludes the G-protein-binding site. The adenine substituent of the agonists binds in a similar fashion to the chemically related region of the inverse agonist ZM241385 (ref. 8). Both agonists contain a ribose group, not found in ZM241385, which extends deep into the ligand-binding pocket where it makes polar interactions with conserved residues in H7 (Ser 277(7.42) and His 278(7.43); superscripts refer to Ballesteros-Weinstein numbering) and non-polar interactions with residues in H3. In contrast, the inverse agonist ZM241385 does not interact with any of these residues and comparison with the agonist-bound structures indicates that ZM241385 sterically prevents the conformational change in H5 and therefore it acts as an inverse agonist. Comparison of the agonist-bound structures of A(2A)R with the agonist-bound structures of β-adrenoceptors indicates that the contraction of the ligand-binding pocket caused by the inward motion of helices 3, 5 and 7 may be a common feature in the activation of all GPCRs.