Docking studies of agonists and antagonists suggest an activation pathway of the A3 adenosine receptor

Interaction of Purine and its Derivatives with A1, A2-Adenosine Receptors and Vascular Endothelial Growth Factor Receptor-1 (Vegf-R1) as a Therapeutic Alternative to Treat Cancer

BACKGROUND

There are several studies that indicate that cancer development may be conditioned by the activation of some biological systems that involve the interaction of different biomolecules, such as adenosine and vascular endothelial growth factor. These biomolecules have been targeted of some drugs for treat of cancer; however, there is little information on the interaction of purine derivatives with adenosine and vascular endothelial growth factor receptor (VEGF-R1).

OBJECTIVE

The aim of this research was to determine the possible interaction of purine (1: ) and their derivatives (2-31: ) with A1, A2-adenosine receptors, and VEGF-R1.

METHODS

Theoretical interaction of purine and their derivatives with A1, A2-adenosine receptors and VEGF-R1 was carried out using the 5uen, 5mzj and 3hng proteins as theoretical tools. Besides, adenosine, cgs-15943, rolofylline, cvt-124, wrc-0571, luf-5834, cvt-6883, AZD-4635, cabozantinib, pazopanib, regorafenib, and sorafenib drugs were used as controls.

RESULTS

The results showed differences in the number of aminoacid residues involved in the interaction of purine and their derivatives with 5uen, 5mzj and 3hng proteins compared with the controls. Besides, the inhibition constants (Ki) values for purine and their derivatives 5: , 9: , 10: , 14: , 15: , 16: , and 20: were lower compared with the controls CONCLUSIONS: Theoretical data suggest that purine and their derivatives 5: , 9: , 10: , 14: , 15: , 16: , and 20: could produce changes in cancer cell growth through inhibition of A1, A2-adenosine receptors and VEGFR-1 inhibition. These data indicate that these purine derivatives could be a therapeutic alternative to treat some types of cancer.

Development of Pharmacophore Models for the Important Off-Target 5-HT2B Receptor

Toxicity is a major cause of attrition in the development of pharmaceuticals, and the off-target effects are a frequent contributor. The 5-HT_2B receptor agonism is known to be responsible for a variety of safety concerns including valvular heart disease which was the cause for the withdrawal of several compounds from the market. An early detection of potential binding to this receptor is thus desirable. Herein, we present the identification of key amino acid residues in the active site of 5-HT_2B by molecular dynamics simulations, the development of pharmacophore models and their performance on in-house data, and a structurally highly diverse subset of Enamine REAL labeled for 5-HT_2B activity by a machine learning model. These models may be used as filters employed on screening compound sets for the early filtration of compounds with potential 5-HT_2B off-target liabilities.

Efficiency of Homology Modeling Assisted Molecular Docking in G-protein Coupled Receptors

BACKGROUND

Molecular docking is in regular practice to assess ligand affinity on a target protein crystal structure. In the absence of protein crystal structure, the homology modeling or comparative modeling is the best alternative to elucidate the relationship details between a ligand and protein at the molecular level. The development of accurate homology modeling (HM) and its integration with molecular docking (MD) is essential for successful, rational drug discovery.

OBJECTIVE

The G-protein coupled receptors (GPCRs) are attractive therapeutic targets due to their immense role in human pharmacology. The GPCRs are membrane-bound proteins with the complex constitution, and the understanding of their activation and inactivation mechanisms is quite challenging. Over the past decade, there has been a rapid expansion in the number of solved G-protein-coupled receptor (GPCR) crystal structures; however, the majority of the GPCR structures remain unsolved. In this context, HM guided MD has been widely used for structure-based drug design (SBDD) of GPCRs.

METHODS

The focus of this review is on the recent (i) developments on HM supported GPCR drug discovery in the absence of GPCR crystal structures and (ii) application of HM in understanding the ligand interactions at the binding site, virtual screening, determining receptor subtype selectivity and receptor behaviour in comparison with GPCR crystal structures.

RESULTS

The HM in GPCRs has been extremely challenging due to the scarcity in template structures. In such a scenario, it is difficult to get accurate HM that can facilitate understanding of the ligand-receptor interactions. This problem has been alleviated to some extent by developing refined HM based on incorporating active /inactive ligand information and inducing protein flexibility. In some cases, HM proteins were found to outscore crystal structures.

CONCLUSION

The developments in HM have been highly operative to gain insights about the ligand interaction at the binding site and receptor functioning at the molecular level. Thus, HM guided molecular docking may be useful for rational drug discovery for the GPCRs mediated diseases.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.

A3 adenosine receptor activation mechanisms: molecular dynamics analysis of inactive, active, and fully active states

We investigated the Gi-coupled A₃ adenosine receptor (A₃AR) activation mechanism by running 7.2 µs of molecular dynamics (MD) simulations. Based on homology to G protein-coupled receptor (GPCR) structures, three constitutively active mutant (CAM) and the wild-type (WT) A₃ARs in the apo form were modeled. Conformational signatures associated with three different receptor states (inactive R, active R*, and bound to Gi protein mimic) were predicted by analyzing and comparing the CAMs with WT receptor and by considering site-directed mutagenesis data. Detected signatures that were correlated with receptor state included: Persistent salt-bridges involving key charged residues for activation (including a novel, putative ionic lock), rotameric state of conserved W^6.48, and Na⁺ ions and water molecules present. Active-coupled state signatures similar to the X-ray structures of β₂ adrenergic receptor-Gs protein and A_2AAR-mini-Gs and the recently solved cryo-EM A₁AR-Gi complexes were found. Our MD analysis suggests that constitutive activation might arise from the D107^3.49-R108^3.50 ionic lock destabilization in R and the D107^3.49-R111^3.53 ionic lock stabilization in R* that presumably lowers the energy barrier associated with an R to R* transition. This study provides new opportunities to understand the underlying interactions of different receptor states of other Gi protein-coupled GPCRs.

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.

Steered molecular dynamics simulation of the binding of the bovine auxilin J domain to the Hsc70 nucleotide-binding domain

The Hsp70 and Hsp40 chaperone machine plays critical roles in protein folding, membrane translocation, and protein degradation by binding and releasing protein substrates in a process that utilizes ATP. The activities of the Hsp70 family of chaperones are recruited and stimulated by the J domains of Hsp40 chaperones. However, structural information on the Hsp40-Hsp70 complex is lacking, and the molecular details of this interaction are yet to be elucidated. Here we used steered molecular dynamics (SMD) simulations to investigate the molecular interactions that occur during the dissociation of the auxilin J domain from the Hsc70 nucleotide-binding domain (NBD). The changes in energy observed during the SMD simulation suggest that electrostatic interactions are the dominant type of interaction. Additionally, we found that Hsp70 mainly interacts with auxilin through the surface residues Tyr866, Arg867, and Lys868 of helix II, His874, Asp876, Lys877, Thr879, and Gln881 of the HPD loop, and Phe891, Asn895, Asp896, and Asn903 of helix III. The conservative residues Tyr866, Arg867, Lys868, His874, Asp876, Lys877, and Phe891 were also found in a previous study to be indispensable to the catalytic activity of the DnaJ J domain and the binding of it with the NBD of DnaK. The in silico identification of the importance of auxilin residues Asn895, Asp896, and Asn903 agrees with previous mutagenesis and NMR data suggesting that helix III of the J domain of the T antigen interacts with Hsp70. Furthermore, our data indicate that Thr879 and Gln881 from the HPD loop are also important as they mediate the interaction between the bovine auxilin J domain and Hsc70.

N6-Substituted 5'-N-Methylcarbamoyl-4'-selenoadenosines as Potent and Selective A3 Adenosine Receptor Agonists with Unusual Sugar Puckering and Nucleobase Orientation

Potent and selective A₃ adenosine receptor (AR) agonists were identified by the replacement of 4'-oxo- or 4'-thionucleosides with bioisosteric selenium. Unlike previous agonists, 4'-seleno analogues preferred a glycosidic syn conformation and South sugar puckering, as shown in the X-ray crystal structure of 5'-N-methylcarbamoyl derivative 3p. Among the compounds tested, N⁶-3-iodobenzyl analogue 3d was found to be the most potent A₃AR full agonist (K_i = 0.57 nM), which was ≥800- and 1900-fold selective for A₁AR and A_2AAR, respectively. In the N⁶-cycloalkyl series, 2-Cl analogues generally exhibited better hA₃AR affinity than 2-H analogues, whereas 2-H > 2-Cl in the N⁶-3-halobenzyl series. N⁷ isomers 3t and 3u were much weaker in binding than corresponding N⁹ isomers, but compound 3t lacked A₃AR activation, appearing to be a weak antagonist. 2-Cl-N⁶-3-iodobenzyl analogue 3p inhibited chemoattractant-induced migration of microglia/monocytes without inducing cell death at ≤50 μM. This suggests the potential for the development of 4'-selenonucleoside A₃AR agonists as novel antistroke agents.

Structural Probing and Molecular Modeling of the A₃ Adenosine Receptor: A Focus on Agonist Binding

Adenosine is an endogenous modulator exerting its functions through the activation of four adenosine receptor (AR) subtypes, termed A₁, A_2A, A_2B and A₃, which belong to the G protein-coupled receptor (GPCR) superfamily. The human A₃AR (hA₃AR) subtype is implicated in several cytoprotective functions. Therefore, hA₃AR modulators, and in particular agonists, are sought for their potential application as anti-inflammatory, anticancer, and cardioprotective agents. Structure-based molecular modeling techniques have been applied over the years to rationalize the structure–activity relationships (SARs) of newly emerged A₃AR ligands, guide the subsequent lead optimization, and interpret site-directed mutagenesis (SDM) data from a molecular perspective. In this review, we showcase selected modeling-based and guided strategies that were applied to elucidate the binding of agonists to the A₃AR and discuss the challenges associated with an accurate prediction of the receptor extracellular vestibule through homology modeling from the available X-ray templates.

Different Types of Molecular Docking Based on Variations of Interacting Molecules

Molecular docking plays an important role in drug discovery research by facilitating target identification, target validation, virtual screening for lead identification and lead optimization. Depending upon the nature of the disease of interest, targets can be either protein or DNA while drugs are mostly organic small molecules. Different types of molecular docking techniques like protein-protein or protein-DNA or protein-small molecule or DNA-small molecule are employed for achieving the above mentioned objectives. This chapter provides a clear idea of the position of molecular docking in drug discovery with detailed discussion on different types of molecular docking based on the varieties of interacting partners. Subsequently the authors provide a detailed list of tools that can be used for docking in drug discovery and discus some examples of molecular docking in drug discovery before concluding with a remark on future areas of improvement in molecular docking related to drug discovery.

Different Types of Molecular Docking Based on Variations of Interacting Molecules

Molecular docking plays an important role in drug discovery research by facilitating target identification, target validation, virtual screening for lead identification and lead optimization. Depending upon the nature of the disease of interest, targets can be either protein or DNA while drugs are mostly organic small molecules. Different types of molecular docking techniques like protein-protein or protein-DNA or protein-small molecule or DNA-small molecule are employed for achieving the above mentioned objectives. This chapter provides a clear idea of the position of molecular docking in drug discovery with detailed discussion on different types of molecular docking based on the varieties of interacting partners. Subsequently the authors provide a detailed list of tools that can be used for docking in drug discovery and discus some examples of molecular docking in drug discovery before concluding with a remark on future areas of improvement in molecular docking related to drug discovery.

Molecular Determinants of CGS21680 Binding to the Human Adenosine A2A Receptor

The adenosine A2A receptor (A(2A)R) plays a key role in transmembrane signaling mediated by the endogenous agonist adenosine. Here, we describe the crystal structure of human A2AR thermostabilized in an active-like conformation bound to the selective agonist 2-[p-(2-carboxyethyl)phenylethyl-amino]-5'-N-ethylcarboxamido adenosine (CGS21680) at a resolution of 2.6 Å. Comparison of A(2A)R structures bound to either CGS21680, 5'-N-ethylcarboxamido adenosine (NECA), UK432097 [6-(2,2-diphenylethylamino)-9-[(2R,3R,4S,5S)-5-(ethylcarbamoyl)-3,4-dihydroxy-tetrahydrofuran-2-yl]-N-[2-[[1-(2-pyridyl)-4-piperidyl]carbamoylamino]ethyl]purine-2-carboxamide], or adenosine shows that the adenosine moiety of the ligands binds to the receptor in an identical fashion. However, an extension in CGS21680 compared with adenosine, the (2-carboxyethyl)phenylethylamino group, binds in an extended vestibule formed from transmembrane regions 2 and 7 (TM2 and TM7) and extracellular loops 2 and 3 (EL2 and EL3). The (2-carboxyethyl)phenylethylamino group makes van der Waals contacts with side chains of amino acid residues Glu169(EL2), His264(EL3), Leu267(7.32), and Ile274(7.39), and the amine group forms a hydrogen bond with the side chain of Ser67(2.65). Of these residues, only Ile274(7.39) is absolutely conserved across the human adenosine receptor subfamily. The major difference between the structures of A(2A)R bound to either adenosine or CGS21680 is that the binding pocket narrows at the extracellular surface when CGS21680 is bound, due to an inward tilt of TM2 in that region. This conformation is stabilized by hydrogen bonds formed by the side chain of Ser67(2.65) to CGS21680, either directly or via an ordered water molecule. Mutation of amino acid residues Ser67(2.65), Glu169(EL2), and His264(EL3), and analysis of receptor activation either in the presence or absence of ligands implicates this region in modulating the level of basal activity of A(2A)R.

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. 

Predicted structures of agonist and antagonist bound complexes of adenosine A3 receptor

We used the GEnSeMBLE Monte Carlo method to predict ensemble of the 20 best packings (helix rotations and tilts) based on the neutral total energy (E) from a vast number (10 trillion) of potential packings for each of the four subtypes of the adenosine G protein-coupled receptors (GPCRs), which are involved in many cytoprotective functions. We then used the DarwinDock Monte Carlo methods to predict the binding pose for the human A(3) adenosine receptor (hAA(3)R) for subtype selective agonists and antagonists. We found that all four A(3) agonists stabilize the 15th lowest conformation of apo-hAA(3)R while also binding strongly to the 1st and 3rd. In contrast the four A(3) antagonists stabilize the 2nd or 3rd lowest conformation. These results show that different ligands can stabilize different GPCR conformations, which will likely affect function, complicating the design of functionally unique ligands. Interestingly all agonists lead to a trans χ1 angle for W6.48 that experiments on other GPCRs associate with G-protein activation while all 20 apo-AA(3)R conformations have a W6.48 gauche+ χ1 angle associated experimentally with inactive GPCRs for other systems. Thus docking calculations have identified critical ligand-GPCR structures involved with activation. We found that the predicted binding site for selective agonist Cl-IB-MECA to the predicted structure of hAA(3)R shows favorable interactions to three subtype variable residues, I253(6.58), V169(EL2), and Q167(EL2), while the predicted structure for hAA(2A)R shows weakened to the corresponding amino acids: T256(6.58), E169(EL2), and L167(EL2), explaining the observed subtype selectivity.

Predicted structures and dynamics for agonists and antagonists bound to serotonin 5-HT2B and 5-HT2C receptors

Subtype 2 serotonin (5-hydroxytryptamine, 5-HT) receptors are major drug targets for schizophrenia, feeding disorders, perception, depression, migraines, hypertension, anxiety, hallucinogens, and gastrointestinal dysfunctions. (1) We report here the predicted structure of 5-HT2B and 5-HT2C receptors bound to highly potent and selective 5-HT2B antagonist PRX-08066 3, (pKi: 30 nM), including the key binding residues [V103 (2.53), L132 (3.29), V190 (4.60), and L347 (6.58)] determining the selectivity of binding to 5-HT2B over 5-HT2A. We also report structures of the endogenous agonist (5-HT) and a HT2B selective antagonist 2 (1-methyl-1-1,6,7,8-tetrahydro-pyrrolo[2,3-g]quinoline-5-carboxylic acid pyridine-3-ylamide). We examine the dynamics for the agonist- and antagonist-bound HT2B receptors in explicit membrane and water finding dramatically different patterns of water migration into the NPxxY motif and the binding site that correlates with the stability of ionic locks in the D(E)RY region.

Toward failure analyses in systems biology

Parallels between designed and biological systems with respect to formal failure analyses have been presented. Failure analysis in designed systems depends on an identified, limited set of parameters or operation variables with high predictive value. In contrast, the biological systems pose problems in identification of operation variables and the identified variables may not be accurate predictors of failure. The difficulty in parameter identification is because of large numbers of components and the inability to envelope variables at each compartment or contour level. Contour level maps for biological systems are currently non-existent, and most failure models are based on very limited, unilateral operation variables (a mutant gene). Operation variable identification within each contour level will enhance failure analyses of complex biological systems.

Design, synthesis, and binding of homologated truncated 4'-thioadenosine derivatives at the human A3 adenosine receptors

We synthesized homologated truncated 4'-thioadenosine analogues 3 in which a methylene (CH(2)) group was inserted in place of the glycosidic bond of a potent and selective A(3) adenosine receptor antagonist 2. The analogues were designed to induce maximum binding interaction in the binding site of the A(3) adenosine receptor. However, all homologated nucleosides were devoid of binding affinity at all subtypes of adenosine receptors, indicating that free rotation through the single bond allowed the compound to adopt an indefinite number of conformations, disrupting the favorable binding interaction essential for receptor recognition.

Medicinal chemistry of the A3 adenosine receptor: agonists, antagonists, and receptor engineering

A(3) adenosine receptor (A(3)AR) ligands have been modified to optimize their interaction with the A(3)AR. Most of these modifications have been made to the N(6) and C2 positions of adenine as well as the ribose moiety, and using a combination of these substitutions leads to the most efficacious, selective, and potent ligands. A(3)AR agonists such as IB-MECA and Cl-IB-MECA are now advancing into Phase II clinical trials for treatments targeting diseases such as cancer, arthritis, and psoriasis. Also, a wide number of compounds exerting high potency and selectivity in antagonizing the human (h)A(3)AR have been discovered. These molecules are generally characterized by a notable structural diversity, taking into account that aromatic nitrogen-containing monocyclic (thiazoles and thiadiazoles), bicyclic (isoquinoline, quinozalines, (aza)adenines), tricyclic systems (pyrazoloquinolines, triazoloquinoxalines, pyrazolotriazolopyrimidines, triazolopurines, tricyclic xanthines) and nucleoside derivatives have been identified as potent and selective A(3)AR antagonists. Probably due to the "enigmatic" physiological role of A(3)AR, whose activation may produce opposite effects (for example, concerning tissue protection in inflammatory and cancer cells) and may produce effects that are species dependent, only a few molecules have reached preclinical investigation. Indeed, the most advanced A(3)AR antagonists remain in preclinical testing. Among the antagonists described above, compound OT-7999 is expected to enter clinical trials for the treatment of glaucoma, while several thiazole derivatives are in development as antiallergic, antiasthmatic and/or antiinflammatory drugs.

Structure-activity relationships of truncated D- and l-4'-thioadenosine derivatives as species-independent A3 adenosine receptor antagonists

Novel D- and l-4'-thioadenosine derivatives lacking the 4'-hydroxymethyl moiety were synthesized, starting from d-mannose and d-gulonic gamma-lactone, respectively, as potent and selective species-independent A 3 adenosine receptor (AR) antagonists. Among the novel 4'-truncated 2-H nucleosides tested, a N(6)-(3-chlorobenzyl) derivative 7c was the most potent at the human A 3 AR (K i = 1.5 nM), but a N(6)-(3-bromobenzyl) derivative 7d showed the optimal species-independent binding affinity.

Selective A(3) adenosine receptor antagonists derived from nucleosides containing a bicyclo[3.1.0]hexane ring system

We have prepared 50-modified derivatives of adenosine and a corresponding (N)-methanocarba nucleoside series containing a bicyclo[3.1.0]hexane ring system in place of the ribose moiety. The compounds were examined in binding assays at three subtypes of adenosine receptors (ARs) and in functional assays at the A3 AR. The H-bonding ability of a group of 9-riboside derivatives containing a 50-uronamide moiety was reduced by modification of the NH; however these derivatives did not display the desired activity as selective A3 AR antagonists, as occurs with 50-N,N-dimethyluronamides. However, truncated (N)-methanocarba analogues lacking a 40-hydroxymethyl group were highly potent and selective antagonists of the human A3 AR. The compounds were synthesized from D-ribose using a reductive free radical decarboxylation of a 50-carboxy intermediate. A less efficient synthetic approach began with L-ribose, which was similar to the published synthesis of (N)-methanocarba A3AR agonists. Compounds 33b-39b (N6-3-halobenzyl and related arylalkyl derivatives) were potent A3AR antagonists with binding Ki values of 0.7-1.4 nM. In a functional assay of [35S]GTPcS binding, 33b (3-iodobenzyl) completely inhibited stimulation by NECA with a KB of 8.9 nM. Thus, a highly potent and selective series of A3AR antagonists has been described.

Selective human adenosine A3 antagonists based on pyrido[2,1-f]purine-2,4-diones: novel features of hA3 antagonist binding

Based on our previous results on the potent antagonist effect of 1H,3H-pyrido[2,1-f]purine-2,4-diones at the human A(3) adenosine receptor, new series of this family of compounds have been synthesized and evaluated in radioligand binding studies against the human A(1), A(2A), A(2B), and A(3) receptors. A remarkable improvement in potency, and most noticeable, in selectivity has been achieved, as exemplified by the 3-cyclopropylmethyl-8-methoxy-1-(4-methylbenzyl)-1H,3H-pyrido[2,1-f]purine-2,4-dione (10) that combines a very high affinity at hA(3) (K(i)=2.24 nM), with lack of affinity for the A(1), A(2A), and A(2B) receptors. On the basis of the published hA(3) receptor model (PDB 1OEA), molecular modeling studies, including molecular dynamics (MD) simulations, have been performed to depict the binding mode of the 1 H,3H-pyrido[2,1-f]purine-2,4-diones and to justify the selectivity against the other adenosine receptors. These studies have led to novel features of the cavity where our antagonists are bound so that the cavity is lined by the hydrogen-bonded Gln 167-Asn 250 pair and by the highly conserved Phe 168.

Molecular modeling of adenosine receptors: new results and trends

Adenosine is a ubiquitous neuromodulator, which carries out its biological task by stimulating four cell surface receptors (A(1), A(2A), A(2B), and A(3)). Adenosine receptors (ARs) are members of the superfamily of G protein-coupled receptors (GPCRs). Their discovery opened up new avenues for potential drug treatment of a variety of conditions such as asthma, neurodegenerative disorders, chronic inflammatory diseases, and many other physiopathological states that are believed to be associated with changes in adenosine levels. Knowledge of the 3D structure of ARs could be of great help in the task of understanding their function and in the rational design of specific ligands. However, since GPCRs are membrane-bound proteins, high-resolution structural characterization is still an extremely difficult task. For this reason, great importance has been placed on molecular modeling studies and, particularly in the last few years, on homology modeling (HM) techniques. The publication of the first high-resolution crystal structure for bovine rhodopsin (bRh), a GPCR superfamily member, provides the option of utilizing HM to generate 3D models based on detailed structural information. In this review we report, analyze, and compare the main experimental data, computational HM procedures and validation methods used for ARs, describing in detail the most successful results.

Molecular modeling of the second extracellular loop of G-protein coupled receptors and its implication on structure-based virtual screening

The current study describes the validation of high-throughput modeling procedures for the construction of the second extracellular loop (ecl2) of all nonolfactory human G Protein-coupled receptors. Our modeling flowchart is based on the alignment of essential residues determining the particular ecl2 fold observed in the bovine rhodopsin (bRho) crystal structure. For a set of GPCR targets, the dopamine D2 receptor (DRD2), adenosine A3 receptor (AA3R), and the thromboxane A2 receptor (TA2R), the implications of including ecl2 atomic coordinates is evaluated in terms of structure-based virtual screening accuracy: the suitability of the 3D models to distinguish between known antagonists and randomly chosen decoys using automated docking approaches. The virtual screening results of different models describing increasingly exhaustive receptor representations (seven helices only, seven helices and ecl2 loop, full model) have been compared. Explicit modeling of the ecl2 loop was found to be important in only one of three test cases whereas a loopless model was shown to be accurate enough in the two other receptors. An exhaustive comparison of ecl2 loops of 365 receptors to that of bRho suggests that explicit ecl2 loop modeling should be reserved to receptors where loop building can be guided by experimental restraints.

Probing distal regions of the A2B adenosine receptor by quantitative structure-activity relationship modeling of known and novel agonists

The binding modes at the A 2B adenosine receptor (AR) of 72 derivatives of adenosine and its 5'- N-methyluronamide with diverse substitutions at the 2 and N (6) positions were studied using a molecular modeling approach. The compounds in their receptor-docked conformations were used to build CoMFA and CoMSIA quantitative structure-activity relationship models. Various parameters, including different types of atomic charges, were examined. The best statistical parameters were obtained with a joint CoMFA and CoMSIA model: R (2) = 0.960, Q (2) = 0.676, SEE = 0.175, F = 158, and R (2) test = 0.782 for an independent test set containing 18 compounds. On the basis of the modeling results, four novel adenosine analogues, having elongated or bulky substitutions at N (6) position and/or 2 position, were synthesized and evaluated biologically. All of the proposed compounds were potent, full agonists at the A 2B AR in adenylate cyclase studies. Thus, in support of the modeling, bulky substitutions at both positions did not prevent A 2B AR activation, which predicts separate regions for docking of these moieties.

Power-law models of signal transduction pathways

The mathematical modelling of signal transduction pathways has become a valuable aid to understanding the complex interactions involved in intracellular signalling mechanisms. An important aspect of the mathematical modelling process is the selection of the model type and structure. Until recently, the convention has been to use a standard kinetic model, often with the Michaelis-Menten steady state assumption. However this model form, although valuable, is only one of a number of choices, and the aim of this article is to consider the mathematical structure and essential features of an alternative model form--the power-law model. Specifically, we analyse how power-law models can be applied to increase our understanding of signal transduction pathways when there may be limited prior information. We distinguish between two kinds of power law models: a) Detailed power-law models, as a tool for investigating pathways when the structure of protein-protein interactions is completely known, and; b) Simplified power-law models, for the analysis of systems with incomplete structural information or insufficient quantitative data for generating detailed models. If sufficient data of high quality are available, the advantage of detailed power-law models is that they are more realistic representations of non-homogenous or crowded cellular environments. The advantages of the simplified power-law model formulation are illustrated using some case studies in cell signalling. In particular, the investigation on the effects of signal inhibition and feedback loops and the validation of structural hypotheses are discussed.

Neoceptors: reengineering GPCRs to recognize tailored ligands

Efforts to model and reengineer the putative binding sites of G-protein-coupled receptors (GPCRs) have led to an approach that combines small-molecule 'classical' medicinal chemistry and gene therapy. In this approach, complementary structural changes (e.g. based on novel ionic or H-bonds) are made in the receptor and ligand for the selective enhancement of affinity. Thus, a modified receptor (neoceptor) is designed for activation by tailor-made agonists that do not interact with the native receptor. The neoceptor is no longer activated by the native agonist, but rather functions as a scaffold for the docking of novel small molecules (neoligands). In theory, the approach could verify the accuracy of GPCR molecular modeling, the investigation of signaling, the design of small molecules to rescue disease-related mutations, and small-molecule-directed gene therapy. The neoceptor-neoligand pairing could offer spatial specificity by delivering the neoceptor to a target site, and temporal specificity by administering neoligand when needed.

Action of nucleosides and nucleotides at 7 transmembrane-spanning receptors

Ribose ring-constrained nucleosides and nucleotides to act at cell-surface purine recesptors have been designed and synthesized. At the P2Y1 nucleotide receptor and the A3 adenosine receptor (AR) the North envelope conformation of ribose is highly preferred. We have applied mutagenesis and rhodopsin-based homology modeling to the study of purine receptors and used the structural insights gained to assist in the design of novel ligands. Two subgroups of P2Y receptors have been defined, containing different sets of cationic residues for coordinating the phosphate groups. Modeling/mutagenesis of adenosine receptors has focused on determinants of intrinsic efficacy in adenosine derivatives and on a conserved Trp residue (6.48) which is involved in the activation process. The clinical use of adenosine agonists as cytoprotective agents has been limited by the widespread occurrence of ARs, thus, leading to undesirable side effects of exogenously administered adenosine derivatives. In order to overcome the inherent nonselectivity of activating the native receptors, we have introduced the concept of neoceptors. By this strategy, intended for eventual use in gene therapy, the putative ligand binding site of a G protein-coupled receptor is reengineered for activation by synthetic agonists (neoligands) built to have a structural complementarity. Using a rational design process we have identified neoceptor-neoligand pairs which are pharmacologically orthogonal with respect to the native species.

Amino-pyrrolidine tricarboxylic acids give new insight into group III metabotropic glutamate receptor activation mechanism

Like most class C G-protein-coupled receptors, metabotropic glutamate (mGlu) receptors possess a large extracellular domain where orthosteric ligands bind. Crystal structures revealed that this domain, called Venus FlyTrap (VFT), adopts a closed or open conformation upon agonist or antagonist binding, respectively. We have described amino-pyrrolidine tricarboxylic acids (APTCs), including (2S,4S)-4-amino-1-[(E)-3-carboxyacryloyl]pyrrolidine-2,4-dicarboxylic acid (FP0429), as new selective group III mGlu agonists. Whereas FP0429 is an almost full mGlu4 agonist, it is a weak and partial agonist of the closely related mGlu8 subtype. To get more insight into the activation mechanism of mGlu receptors, we aimed to elucidate why FP0429 behaves differently at these two highly homologous receptors by focusing on two residues within the binding site that differ between mGlu4 and mGlu8. Site-directed mutagenesis of Ser157 and Gly158 of mGlu4 into their mGlu8 homologs (Ala) turned FP0429 into a weak partial agonist. Conversely, introduction of Ser and Gly residues into mGlu8 increased FP0429 efficacy. Docking of FP0429 in mGlu4 VFT 3D model helped us characterize the role of each residue. Indeed, mGlu4 Ser157 seems to have an important role in FP0429 binding, whereas Gly158 may allow a deeper positioning of this agonist in the cavity of lobe I, thereby ensuring optimal interactions with lobe II residues in the fully closed state of the VFT. In contrast, the presence of a methyl group in mGlu8 (Ala instead of Gly) weakens the interactions with the lobe II residues. This probably results in a less stable or a partially closed form of the mGlu8 VFT, leading to partial receptor activation.

Molecular dynamics simulation of the human adenosine A3 receptor: agonist induced conformational changes of Trp243

The adenosine A(3) receptor together with rhodopsin belongs to Class A of the G-protein coupled receptors. As the crystal structure of bovine rhodopsin represents the dark (inactive) state of the receptor, the details of GPCR activation are still unknown. In this molecular dynamics study we investigate how the homology model of the human adenosine A(3) receptor responds to ligand exposure. To this end we placed the homology model in a POPC membrane model. After equilibrating for 13 ns an agonist (Cl-IB-MECA) and an inverse agonist (PSB-10) were placed inside the putative binding pocket. In the following 10 ns molecular dynamics simulation we observed a different behaviour of the side-chain torsions of Trp243(6.48), depending on the presence or absence of the agonist or inverse agonist. This conformational change of Trp243 correlates with the assumed influence of ligands on receptor activation. Other predicted conformational changes of the receptor could not be observed yet. So Trp243 may represent the first switch in receptor activation.