Design, Synthesis and Biological Evaluation of 1,3,5-Triazine Derivatives Targeting hA1 and hA3 Adenosine Receptor.

Adenosine mediates various physiological activities in the body. Adenosine receptors (ARs) are widely expressed in tumors and the tumor microenvironment (TME), and they induce tumor proliferation and suppress immune cell function. There are four types of human adenosine receptor (hARs): hA1, hA2A, hA2B, and hA3. Both hA1 and hA3 AR play an important role in tumor proliferation. We designed and synthesized novel 1,3,5-triazine derivatives through amination and Suzuki coupling, and evaluated them for binding affinities to each hAR subtype. Compounds 9a and 11b showed good binding affinity to both hA1 and hA3 AR, while 9c showed the highest binding affinity to hA1 AR. In this study, we discovered that 9c inhibits cell viability, leading to cell death in lung cancer cell lines. Flow cytometry analysis revealed that 9c caused an increase in intracellular reactive oxygen species (ROS) and a depolarization of the mitochondrial membrane potential. The binding mode of 1,3,5-triazine derivatives to hA1 and hA3 AR were predicted by a molecular docking study.

Development of Bicyclo[3.1.0]hexane-Based A3 Receptor Ligands: Closing the Gaps in the Structure-Affinity Relationships.

The adenosine A3 receptor is a promising target for treating and diagnosing inflammation and cancer. In this paper, a series of bicyclo[3.1.0]hexane-based nucleosides was synthesized and evaluated for their P1 receptor affinities in radioligand binding studies. The study focused on modifications at 1-, 2-, and 6-positions of the purine ring and variations of the 5'-position at the bicyclo[3.1.0]hexane moiety, closing existing gaps in the structure-affinity relationships. The most potent derivative 30 displayed moderate A3AR affinity (Ki of 0.38 µM) and high A3R selectivity. A subset of compounds varied at 5'-position was further evaluated in functional P2Y1R assays, displaying no off-target activity.

Subtle Chemical Changes Cross the Boundary between Agonist and Antagonist: New A3 Adenosine Receptor Homology Models and Structural Network Analysis Can Predict This Boundary.

Distinguishing compounds' agonistic or antagonistic behavior would be of great utility for the rational discovery of selective modulators. We synthesized truncated nucleoside derivatives and discovered 6c (Ki = 2.40 nM) as a potent human A3 adenosine receptor (hA3AR) agonist, and subtle chemical modification induced a shift from antagonist to agonist. We elucidated this shift by developing new hA3AR homology models that consider the pharmacological profiles of the ligands. Taken together with molecular dynamics (MD) simulation and three-dimensional (3D) structural network analysis of the receptor-ligand complex, the results indicated that the hydrogen bonding with Thr943.36 and His2727.43 could make a stable interaction between the 3'-amino group with TM3 and TM7, and the corresponding induced-fit effects may play important roles in rendering the agonistic effect. Our results provide a more precise understanding of the compounds' actions at the atomic level and a rationale for the design of new drugs with specific pharmacological profiles.

Development of Covalent, Clickable Probes for Adenosine A1 and A3 Receptors.

Adenosine receptors are attractive therapeutic targets for multiple conditions, including ischemia-reperfusion injury and neuropathic pain. Adenosine receptor drug discovery efforts would be facilitated by the development of appropriate tools to assist in target validation and direct receptor visualization in different native environments. We report the development of the first bifunctional (chemoreactive and clickable) ligands for the adenosine A1 receptor (A1R) and adenosine A3 receptor (A3R) based on an orthosteric antagonist xanthine-based scaffold and on an existing structure-activity relationship. Bifunctional ligands were functional antagonists with nanomolar affinity and irreversible binding at the A1R and A3R. In-depth pharmacological profiling of these bifunctional ligands showed moderate selectivity over A2A and A2B adenosine receptors. Once bound to the receptor, ligands were successfully "clicked" with a cyanine-5 fluorophore containing the complementary "click" partner, enabling receptor detection. These bifunctional ligands are expected to aid in the understanding of A1R and A3R localization and trafficking in native cells and living systems.

Pharmacological characterisation of novel adenosine A3 receptor antagonists.

The adenosine A3 receptor (A3R) belongs to a family of four adenosine receptor (AR) subtypes which all play distinct roles throughout the body. A3R antagonists have been described as potential treatments for numerous diseases including asthma. Given the similarity between (adenosine receptors) orthosteric binding sites, obtaining highly selective antagonists is a challenging but critical task. Here we screen 39 potential A3R, antagonists using agonist-induced inhibition of cAMP. Positive hits were assessed for AR subtype selectivity through cAMP accumulation assays. The antagonist affinity was determined using Schild analysis (pA2 values) and fluorescent ligand binding. Structure-activity relationship investigations revealed that loss of the 3-(dichlorophenyl)-isoxazolyl moiety or the aromatic nitrogen heterocycle with nitrogen at α-position to the carbon of carboximidamide group significantly attenuated K18 antagonistic potency. Mutagenic studies supported by molecular dynamic simulations combined with Molecular Mechanics-Poisson Boltzmann Surface Area calculations identified the residues important for binding in the A3R orthosteric site. We demonstrate that K18, which contains a 3-(dichlorophenyl)-isoxazole group connected through carbonyloxycarboximidamide fragment with a 1,3-thiazole ring, is a specific A3R (< 1 µM) competitive antagonist. Finally, we introduce a model that enables estimates of the equilibrium binding affinity for rapidly disassociating compounds from real-time fluorescent ligand-binding studies. These results demonstrate the pharmacological characterisation of a selective competitive A3R antagonist and the description of its orthosteric binding mode. Our findings may provide new insights for drug discovery.

Adenosine Receptor Ligands: Coumarin-Chalcone Hybrids as Modulating Agents on the Activity of hARs.

Adenosine receptors (ARs) play an important role in neurological and psychiatric disorders such as Alzheimer's disease, Parkinson's disease, epilepsy and schizophrenia. The different subtypes of ARs and the knowledge on their densities and status are important for understanding the mechanisms underlying the pathogenesis of diseases and for developing new therapeutics. Looking for new scaffolds for selective AR ligands, coumarin-chalcone hybrids were synthesized (compounds 1-8) and screened in radioligand binding (hA1, hA2A and hA3) and adenylyl cyclase (hA2B) assays in order to evaluate their affinity for the four human AR subtypes (hARs). Coumarin-chalcone hybrid has been established as a new scaffold suitable for the development of potent and selective ligands for hA1 or hA3 subtypes. In general, hydroxy-substituted hybrids showed some affinity for the hA1, while the methoxy counterparts were selective for the hA3. The most potent hA1 ligand was compound 7 (Ki = 17.7 µM), whereas compound 4 was the most potent ligand for hA3 (Ki = 2.49 µM). In addition, docking studies with hA1 and hA3 homology models were established to analyze the structure-function relationships. Results showed that the different residues located on the protein binding pocket could play an important role in ligand selectivity.

Functional Expression of Adenosine A3 Receptor in Yeast Utilizing a Chimera with the A2AR C-Terminus.

The adenosine A3 receptor (A3R) is the only adenosine receptor subtype to be overexpressed in inflammatory and cancer cells and therefore is considered a novel and promising therapeutic target for inflammatory diseases and cancer. Heterologous expression of A3R at levels to allow biophysical characterization is a major bottleneck in structure-guided drug discovery efforts. Here, we apply protein engineering using chimeric receptors to improve expression and activity in yeast. Previously we had reported improved expression and trafficking of the chimeric A1R variant using a similar approach. In this report, we constructed chimeric A3/A2AR comprising the N-terminus and transmembrane domains from A3R (residues 1-284) and the cytoplasmic C-terminus of the A2AR (residues 291-412). The chimeric receptor showed approximately 2-fold improved expression with a 2-fold decreased unfolded protein response when compared to wild type A3R. Moreover, by varying culture conditions such as initial cell density and induction temperature a further 1.7-fold increase in total receptor yields was obtained. We observed native-like coupling of the chimeric receptor to Gai-Gpa1 in engineered yeast strains, activating the downstream, modified MAPK pathway. This strategy of utilizing chimeric receptor variants in yeast thus provides an exciting opportunity to improve expression and activity of "difficult-to-express" receptors, expanding the opportunity for utilizing yeast in drug discovery.

Truncated (N)-Methanocarba Nucleosides as Partial Agonists at Mouse and Human A3 Adenosine Receptors: Affinity Enhancement by N6-(2-Phenylethyl) Substitution.

Dopamine-derived N6-substituents, compared to N6-(2-phenylethyl), in truncated (N)-methanocarba (bicyclo[3.1.0]hexyl) adenosines favored high A3 adenosine receptor (AR) affinity/selectivity, e.g., C2-phenylethynyl analogue 15 (MRS7591, Ki = 10.9/17.8 nM, at human/mouse A3AR). 15 was a partial agonist in vitro (hA3AR, cAMP inhibition, 31% Emax; mA3AR, [35S]GTP-γ-S binding, 16% Emax) and in vivo and also antagonized hA3AR in vitro. Distal H-bonding substitutions of the N6-(2-phenylethyl) moiety particularly enhanced mA3AR affinity by polar interactions with the extracellular loops, predicted using docking and molecular dynamics simulation with newly constructed mA3AR and hA3AR homology models. These hybrid models were based on an inactive antagonist-bound hA1AR structure for the upper part of TM2 and an agonist-bound hA2AAR structure for the remaining TM portions. These species-independent A3AR-selective nucleosides are low efficacy partial agonists and novel, nuanced modulators of the A3AR, a drug target of growing interest.

A Taxicab geometry quantification system to evaluate the performance of in silico methods: a case study on adenosine receptors ligands.

Among still comparatively few G protein-coupled receptors, the adenosine A2A receptor has been co-crystallized with several ligands, agonists as well as antagonists. It can thus serve as a template with a well-described orthosteric ligand binding region for adenosine receptors. As not all subtypes have been crystallized yet, and in order to investigate the usability of homology models in this context, multiple adenosine A1 receptor (A1AR) homology models had been previously obtained and a library of lead-like compounds had been docked. As a result, a number of potent and one selective ligand toward the intended target have been identified. However, in in vitro experimental verification studies, many ligands also bound to the A2AAR and the A3AR subtypes. In this work we asked the question whether a classification of the ligands according to their selectivity was possible based on docking scores. Therefore, we built an A3AR homology model and docked all previously found ligands to all three receptor subtypes. As a metric, we employed an in vitro/in silico selectivity ranking system based on taxicab geometry and obtained a classification model with reasonable separation. In the next step, the method was validated with an external library of, selective ligands with similarly good performance. This classification system might also be useful in further screens.

Structure-Based Optimization of Coumarin hA3 Adenosine Receptor Antagonists.

Adenosine receptors participate in many physiological functions. Molecules that may selectively interact with one of the receptors are favorable multifunctional chemical entities to treat or decelerate the evolution of different diseases. 3-Arylcoumarins have already been studied as neuroprotective agents by our group. Here, differently 8-substituted 3-arylcoumarins are complementarily studied as ligands of adenosine receptors, performing radioligand binding assays. Among the synthesized compounds, selective A3 receptor antagonists were found. 3-(4-Bromophenyl)-8-hydroxycoumarin (compound 4) displayed the highest potency and selectivity as A3 receptor antagonist (Ki = 258 nM). An analysis of its X-ray diffraction provided detailed information on its structure. Further evaluation of a selected series of compounds indicated that it is the nature and position of the substituents that determine their activity and selectivity. Theoretical modeling calculations corroborate and explain the experimental data, suggesting this novel scaffold can be involved in the generation of candidates as multitarget drugs.

Application of Fluorescent Purinoceptor Antagonists for Bioluminescence Resonance Energy Transfer Assays and Fluorescent Microscopy.

Fluorescent antagonists offer the ability to interrogate G protein-coupled receptor pharmacology. With resonance energy transfer techniques, fluorescent antagonists can be implemented to monitor receptor-ligand interactions using assays originally designed for radiolabeled probes. The fluorescent nature of these antagonists also enables the localization and distribution of the receptors to be visualized in living cells. Here, we describe the generation of modified purinergic receptors with the NanoLuc luciferase or SNAP-tag, using the P1 adenosine A3 receptor as an example. We also describe the procedure of characterizing a novel fluorescent purinergic antagonist using ligand-mediated bioluminescence resonance energy transfer assays and confocal microscopy.

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

We investigated the Gi-coupled A3 adenosine receptor (A3AR) 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) A3ARs 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 W6.48, and Na+ ions and water molecules present. Active-coupled state signatures similar to the X-ray structures of ß2 adrenergic receptor-Gs protein and A2AAR-mini-Gs and the recently solved cryo-EM A1AR-Gi complexes were found. Our MD analysis suggests that constitutive activation might arise from the D1073.49-R1083.50 ionic lock destabilization in R and the D1073.49-R1113.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.

Insights to the Binding of a Selective Adenosine A3 Receptor Antagonist Using Molecular Dynamic Simulations, MM-PBSA and MM-GBSA Free Energy Calculations, and Mutagenesis.

Adenosine A3 receptor (A3R) is a promising drug target cancer and for a number of other conditions like inflammatory diseases, including asthma and rheumatoid arthritis, glaucoma, chronic obstructive pulmonary disease, and ischemic injury. Currently, there is no experimentally determined structure of A3R. We explored the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18), which is a new specific and competitive antagonist at the orthosteric binding site of A3R. MD simulations and MM-GBSA calculations of the WT A3R in complex with K18 combined with in vitro mutagenic studies show that the most plausible binding conformation for the dichlorophenyl group of K18 is oriented toward trans-membrane helices (TM) 5, 6 and reveal important residues for binding. Further, MM-GBSA calculations distinguish mutations that reduce or maintain or increase antagonistic activity. Our studies show that selectivity of K18 toward A3R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V1695.30 is considered to be a selectivity filter for A3R binders, when it was mutated to glutamic acid, K18 maintained antagonistic potency, in agreement with our previous results obtained for agonists binding profile investigation. Mutation of the direct interacting residue L903.32 in the low region and the remote L2647.35 in the middle/upper region to alanine increases antagonistic potency, suggesting an empty space in the orthosteric area available for increasing antagonist potency. These results approve the computational model for the description of K18 binding at A3R, which we previously performed for agonists binding to A3R, and the design of more effective antagonists based on K18.

Structural Characterization of Agonist Binding to an A3 Adenosine Receptor through Biomolecular Simulations and Mutagenesis Experiments.

The adenosine A3 receptor (A3R) binds adenosine and is a drug target against cancer cell proliferation. Currently, there is no experimental structure of A3R. Here, we have generated a molecular model of A3R in complex with two agonists, the nonselective 1-(6-amino-9H-purin-9-yl)-1-deoxy-N-ethyl-ß-d-ribofuranuronamide (NECA) and the selective 1-deoxy-1-[6-[[(3-iodophenyl)methyl]amino]-9H-purin-9-yl]-N-methyl-ß-d-ribofuranuronamide (IB-MECA). Molecular dynamics simulations of the wild-type A3R in complex with both agonists, combined with in vitro mutagenic studies revealed important residues for binding. Further, molecular mechanics-generalized Born surface area calculations were able to distinguish mutations that reduce or negate agonistic activity from those that maintained or increased the activity. Our studies reveal that selectivity of IB-MECA toward A3R requires not only direct interactions with residues within the orthosteric binding area but also with remote residues. Although V1695.30 is considered to be a selectivity filter for A3R binders, when it was mutated to glutamic acid or alanine, the activity of IB-MECA increased by making new van der Waals contacts with TM5. This result may have implications in the design of new A3R agonists.

Free-Energy Calculations for Bioisosteric Modifications of A3 Adenosine Receptor Antagonists.

Adenosine receptors are a family of G protein-coupled receptors with increased attention as drug targets on different indications. We investigate the thermodynamics of ligand binding to the A3 adenosine receptor subtype, focusing on a recently reported series of diarylacetamidopyridine inhibitors via molecular dynamics simulations. With a combined approach of thermodynamic integration and one-step perturbation, we characterize the impact of the charge distribution in a central heteroaromatic ring on the binding affinity prediction. Standard charge distributions according to the GROMOS force field yield values in good agreement with the experimental data and previous free energy calculations. Subsequently, we examine the thermodynamics of inhibitor binding in terms of the energetic and entropic contributions. The highest entropy penalties are found for inhibitors with methoxy substituents in meta position of the aryl groups. This bulky group restricts rotation of aromatic rings attached to the pyrimidine core which leads to two distinct poses of the ligand. Our predictions support the previously proposed binding pose for the o-methoxy ligand, yielding in this case a very good correlation with the experimentally measured affinities with deviations below 4 kJ/mol.

7-Amino-2-aryl/hetero-aryl-5-oxo-5,8-dihydro[1,2,4]triazolo[1,5-a]pyridine-6-carbonitriles: Synthesis and adenosine receptor binding studies.

A series of novel 7-amino-5-oxo-2-substituted-aryl/hetero-aryl-5,8-dihydro[1,2,4]triazolo[1,5-a]pyridine-6-carbonitriles (4a-4t) was synthesized, characterized and evaluated for their binding affinity and selectivity towards hA1 , hA2A , hA2B and hA3 adenosine receptors (ARs). Compound 4a with a phenyl ring at 2-position of the triazolo moiety of the scaffold showed high affinity and selectivity for hA1 AR (Ki hA1  = 0.076 µM, hA2A  = 25.6 µM and hA3  > 100 µM). Introduction of various electron donating and withdrawing groups at different positions of the phenyl ring resulted in drastic reduction in affinity and selectivity towards all the ARs, except compound 4b with a 4-hydroxyphenyl group at 2-position. Interestingly, the replacement of the phenyl ring with a smaller heterocyclic thiophene ring (π excessive system) resulted in further improvement of affinity for hA1 AR of compound 4t (Ki hA1  = 0.051 µM, hA2A  = 9.01 µM and hA3  > 13.9 µM) while retaining the significant selectivity against all other AR subtypes similar to compound 4a. The encouraging results for compounds 4a and 4t indicate that substitution at 2-position of the scaffold with π-excessive systems other than thiophene may lead to even more potent and selective hA1 AR antagonists.

Comparative Study of Carborane- and Phenyl-Modified Adenosine Derivatives as Ligands for the A2A and A3 Adenosine Receptors Based on a Rigid in Silico Docking and Radioligand Replacement Assay.

Adenosine receptors are involved in many physiological processes and pathological conditions and are therefore attractive therapeutic targets. To identify new types of effective ligands for these receptors, a library of adenosine derivatives bearing a boron cluster or phenyl group in the same position was designed. The ligands were screened in silico to determine their calculated affinities for the A2A and A3 adenosine receptors. An virtual screening protocol based on the PatchDock web server was developed. In the first screening phase, the effects of the functional group (organic or inorganic modulator) on the adenosine ligand affinity for the receptors were determined. Then, the lead compounds were identified for each receptor in the second virtual screening phase. Two pairs of the most promising ligands, compounds 3 and 4, and two ligands with lower affinity scores (compounds 11 and 12, one with a boron cluster and one with a phenyl group) were synthesized and tested in a radioligand replacement assay for affinity to the A2A and A3 receptors. A reasonable correlation of in silico and biological assay results was observed. In addition, the effects of a phenyl group and boron cluster, which is new adenosine modifiers, on the adenosine ligand binding were compared.

Effect of Nitrogen Atom Substitution in A3 Adenosine Receptor Binding: N-(4,6-Diarylpyridin-2-yl)acetamides as Potent and Selective Antagonists.

We report the first family of 2-acetamidopyridines as potent and selective A3 adenosine receptor (AR) antagonists. The computer-assisted design was focused on the bioisosteric replacement of the N1 atom by a CH group in a previous series of diarylpyrimidines. Some of the generated 2-acetamidopyridines elicit an antagonistic effect with excellent affinity (Ki < 10 nM) and outstanding selectivity profiles, providing an alternative and simpler chemical scaffold to the parent series of diarylpyrimidines. In addition, using molecular dynamics and free energy perturbation simulations, we elucidate the effect of the second nitrogen of the parent diarylpyrimidines, which is revealed as a stabilizer of a water network in the binding site. The discovery of 2,6-diaryl-2-acetamidopyridines represents a step forward in the search of chemically simple, potent, and selective antagonists for the hA3AR, and exemplifies the benefits of a joint theoretical-experimental approach to identify novel hA3AR antagonists through succinct and efficient synthetic methodologies.

Alpha-bisabolol Promotes Glioma Cell Death by Modulating the Adenosinergic System.

BACKGROUND/AIM: Glioblastoma multiforme is the most malignant type of glioma. Alpha-bisabolol is an essential oil reported as a potent cell death agent. In the present work, we evaluated the effect of alpha-bisabolol on ecto-5'-nucleotidase/CD73, the most well-characterized enzymatic source of adenosine, present in lipid rafts. MATERIALS AND METHODS: Glioma cells were treated with alpha-bisabolol and, in some experiments, pre-treated with an A3 antagonist. MTT assay (viability), malachite green method (ecto-5'-nucleotidase/CD73 activity) and quantitative polymerase chain reaction (qPCR) (A3 mRNA) were carried out. RESULTS: Alpha-bisabolol led to a decrease in C6 and U138-MG glioma cells viability, accompanied by an increase in ecto-5'-NT/CD73 activity. Pre-treatment with an A3 antagonist reverted the effect of α-bisabolol with an increase of mRNA expression of this receptor. CONCLUSION: Our data indicated the participation of ecto-5'-nucleotidase/CD73 and A3 receptor in the anti-proliferative effect of α-bisabolol on glioma cells.

Structural Probing and Molecular Modeling of the A3 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 A1, A2A, A2B and A3, which belong to the G protein-coupled receptor (GPCR) superfamily. The human A3AR (hA3AR) subtype is implicated in several cytoprotective functions. Therefore, hA3AR 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 A3AR 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 A3AR and discuss the challenges associated with an accurate prediction of the receptor extracellular vestibule through homology modeling from the available X-ray templates.