Crystal Structure and Subsequent Ligand Design of a Nonriboside Partial Agonist Bound to the Adenosine A2A Receptor

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Discovery of novel thiophene-3-carboxamide derivatives as potential VEGFR-2 inhibitors with anti-angiogenic properties

VEGFR-2 is an attractive target for the development of anti-tumor drugs and plays a crucial role in tumor angiogenesis. This study reports a series of novel thiophene-3-carboxamide derivatives based on PAN-90806 as VEGFR-2 inhibitors, among which compound 14d exhibits excellent anti-proliferative activity against HCT116, MCF7, PC3, and A549 cell lines, and has effective VEGFR-2 inhibitory activity with an IC50 value of 191.1 nM. Additionally, CETSA results indicated that VEGFR-2 was a relevant target of compound 14d in the cell lines, and compound 14d could also inhibit VEGFR-2 protein phosphorylation in A549 cell line. Furthermore, compound 14d inhibited colony formation, cell migration, and HUVECs tube formation in a dose-dependent manner. The mechanism by which 14d induced cancer cell death involves blocking the cell cycle, increasing ROS production, inducing apoptosis, and dose-dependently reducing the levels of phosphorylated ERK and MEK. Molecular docking and molecular dynamics simulations had shown that compound 14d could stably bind to the active site of VEGFR-2. These results confirmed that compound 14d might be a promising lead compound for anti-angiogenesis.

Structural Insights into Partial Activation of the Prototypic G Protein-Coupled Adenosine A2A Receptor

The adenosine A2A receptor (A2AAR) belongs to the rhodopsin-like G protein-coupled receptor (GPCR) family, which constitutes the largest class of GPCRs. Partial agonists show reduced efficacy as compared to physiological agonists and can even act as antagonists in the presence of a full agonist. Here, we determined an X-ray crystal structure of the partial A2AAR agonist 2-amino-6-[(1H-imidazol-2-ylmethyl)sulfanyl]-4-p-hydroxyphenyl-3,5-pyridinedicarbonitrile (LUF5834) in complex with the A2AAR construct A2A-PSB2-bRIL, stabilized in its inactive conformation and being devoid of any mutations in the ligand binding pocket. The determined high-resolution structure (2.43 Å) resolved water networks and crucial binding pocket interactions. A direct hydrogen bond of the p-hydroxy group of LUF5834 with T883.36 was observed, an amino acid that was mutated to alanine in the most frequently used A2AAR crystallization constructs thus preventing the discovery of its interactions in most of the previous A2AAR co-crystal structures. G protein dissociation studies confirmed partial agonistic activity of LUF5834 as compared to that of the full agonist N-ethylcarboxamidoadenosine (NECA). In contrast to NECA, the partial agonist was still able to bind to the receptor construct locked in its inactive conformation by an S913.39K mutation, although with an affinity lower than that at the native receptor. This could explain the compound's partial agonistic activity: while full A2AAR agonists bind exclusively to the active conformation, likely following conformational selection, partial agonists bind to active as well as inactive conformations, showing higher affinity for the active conformation. This might be a general mechanism of partial agonism also applicable to other GPCRs.

Membrane mimetic-dependence of GPCR energy landscapes

We leveraged variable-temperature 19F-NMR spectroscopy to compare the conformational equilibria of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), across a range of temperatures ranging from lower temperatures typically employed in 19F-NMR experiments to physiological temperature. A2AAR complexes with partial agonists and full agonists showed large increases in the population of a fully active conformation with increasing temperature. NMR data measured at physiological temperature were more in line with functional data. This was pronounced for complexes with partial agonists, where the population of active A2AAR was nearly undetectable at lower temperature but became evident at physiological temperature. Temperature-dependent behavior of complexes with either full or partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. Cellular signaling experiments correlated with the temperature-dependent conformational equilibria of A2AAR in lipid nanodiscs but not in some detergents, underscoring the importance of the membrane environment in studies of GPCR function.

PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.

Crystal structure reveals the binding mode and selectivity of a photoswitchable ligand for the adenosine A2A receptor

Rational synthetic expansion of photoresponsive ligands is important for photopharmacological studies. Adenosine A2A receptor (A2AR) is stimulated by adenosine and related in Parkinson's disease and other diseases. Here, we report the crystal structure of the A2AR in complex with the novel photoresponsive ligand photoNECA (blue) at 3.34 Å resolution. PhotoNECA (blue) was designed for this structural study and the cell-based assay showed a photoresponsive and receptor selective characteristics of photoNECA (blue) for A2AR. The crystal structure explains the binding mode, photoresponsive mechanism and receptor selectivity of photoNECA (blue). Our study would promote not only the rational design of photoresponsive ligands but also dynamic structural studies of A2AR.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

Functional GPCR Expression in Eukaryotic LEXSY System

G protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins in the human genome, and represent one of the most important classes of drug targets. Their structural studies facilitate rational drug discovery. However, atomic structures of only about 20% of human GPCRs have been solved to date. Recombinant production of GPCRs for structural studies at a large scale is challenging due to their low expression levels and stability. Therefore, in this study, we explored the efficacy of the eukaryotic system LEXSY (Leishmania tarentolae) for GPCR production. We selected the human A2A adenosine receptor (A2AAR), as a model protein, expressed it in LEXSY, purified it, and compared with the same receptor produced in insect cells, which is the most popular expression system for structural studies of GPCRs. The A2AAR purified from both expression systems showed similar purity, stability, ligand-induced conformational changes and structural dynamics, with a remarkably higher protein yield in the case of LEXSY expression. Overall, our results suggest that LEXSY is a promising platform for large-scale production of GPCRs for structural studies.

A2A adenosine receptor agonists, antagonists, inverse agonists and partial agonists

The Gs-coupled A2A adenosine receptor (A2AAR) has been explored extensively as a pharmaceutical target, which has led to numerous clinical trials. However, only one selective A2AAR agonist (regadenoson, Lexiscan) and one selective A2AAR antagonist (istradefylline, Nouriast) have been approved by the FDA, as a pharmacological agent for myocardial perfusion imaging (MPI) and as a cotherapy for Parkinson's disease (PD), respectively. Adenosine is widely used in MPI, as Adenoscan. Despite numerous unsuccessful clinical trials, medicinal chemical activity around A2AAR ligands has accelerated recently, particularly through structure-based drug design. New drug-like A2AAR antagonists for PD and cancer immunotherapy have been identified, and many clinical trials have ensued. For example, imaradenant (AZD4635), a compound that was designed computationally, based on A2AAR X-ray structures and biophysical mapping. Mixed A2AAR/A2BAR antagonists are also hopeful for cancer treatment. A2AAR antagonists may also have potential as neuroprotective agents for treatment of Alzheimer's disease.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

New paradigms in purinergic receptor ligand discovery

The discovery and clinical implementation of modulators of adenosine, P2Y and P2X receptors (comprising nineteen subtypes) have progressed dramatically in ∼50 years since Burnstock's definition of purinergic signaling. Although most clinical trials of selective ligands (agonists and antagonists) of certain purinergic receptors failed, there is a renewed impetus to redirect efforts to new disease conditions and the discovery of more selective or targeted compounds with potentially reduced side effects, such as biased GPCR agonists. The elucidation of new receptor and enzyme structures is steering rational design of potent and selective agonists, antagonists, allosteric modulators and inhibitors. A2A adenosine receptor (AR) antagonists are being applied to neurodegenerative conditions and cancer immunotherapy. A3AR agonists have potential for treating chronic inflammation (e.g. psoriasis), stroke and pain, as well as cancer. P2YR modulators are being considered for treating inflammation, metabolic disorders, acute kidney injury, cancer, pain and other conditions, often with an immune mechanism. ADP-activated P2Y12R antagonists are widely used as antithrombotic drugs, while their repurposing toward neuroinflammation is considered. P2X3 antagonists have been in clinical trials for chronic cough. P2X7 antagonists have been in clinical trials for inflammatory diseases and depression (compounds that penetrate the blood-brain barrier). Thus, purinergic signaling is now recognized as an immense regulatory system in the body for rebalancing tissues and organs under stress, which can be adjusted by drug intervention for therapeutic purposes. The lack of success of many previous clinical trials can be overcome given more advanced pharmacokinetic and pharmacodynamic approaches, including structure-based drug design, prodrugs and biased signaling. This article is part of the Special Issue on "Purinergic Signaling: 50 years".

Intermediate-state-trapped mutants pinpoint G protein-coupled receptor conformational allostery

Understanding the roles of intermediate states in signaling is pivotal to unraveling the activation processes of G protein-coupled receptors (GPCRs). However, the field is still struggling to define these conformational states with sufficient resolution to study their individual functions. Here, we demonstrate the feasibility of enriching the populations of discrete states via conformation-biased mutants. These mutants adopt distinct distributions among five states that lie along the activation pathway of adenosine A2A receptor (A2AR), a class A GPCR. Our study reveals a structurally conserved cation-π lock between transmembrane helix VI (TM6) and Helix8 that regulates cytoplasmic cavity opening as a "gatekeeper" for G protein penetration. A GPCR activation process based on the well-discerned conformational states is thus proposed, allosterically micro-modulated by the cation-π lock and a previously well-defined ionic interaction between TM3 and TM6. Intermediate-state-trapped mutants will also provide useful information in relation to receptor-G protein signal transduction.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

Systematic analyses of the sequence conservation and ligand interaction patterns of purinergic P1 and P2Y receptors provide a structural basis for receptor selectivity

Purinergic receptors are membrane proteins that regulate numerous cellular functions by catalyzing reactions involving purine nucleotides or nucleosides. Among the three receptor families, i.e., P1, P2X, and P2Y, the P1 and P2Y receptors share common structural features of class A GPCR. Comprehensive sequence and structural analysis revealed that the P1 and P2Y receptors belong to two distinct groups. They exhibit different ligand-binding site features that can distinguish between specific activators. These specific amino acid residues in the binding cavity may be involved in the selectivity and unique pharmacological behavior of each subtype. In this study, we conducted a structure-based analysis of purinergic P1 and P2Y receptors to identify their evolutionary signature and obtain structural insights into ligand recognition and selectivity. The structural features of the P1 and P2Y receptor classes were compared based on sequence conservation and ligand interaction patterns. Orthologous protein sequences were collected for the P1 and P2Y receptors, and sequence conservation was calculated based on Shannon entropy to identify highly conserved residues. To analyze the ligand interaction patterns, we performed docking studies on the P1 and P2Y receptors using known ligand information extracted from the ChEMBL database. We analyzed how the conserved residues are related to ligand-binding sites and how the key interacting residues differ between P1 and P2Y receptors, or between agonists and antagonists. We extracted new similarities and differences between the receptor subtypes, and the results can be used for designing new ligands by predicting hotspot residues that are important for functional selectivity.

A Chemical Biological Approach to Study G Protein-Coupled Receptors: Labeling the Adenosine A1 Receptor Using an Electrophilic Covalent Probe

G protein-coupled receptors (GPCRs) have been known for decades as attractive drug targets. This has led to the development and approval of many ligands targeting GPCRs. Although ligand binding effects have been studied thoroughly for many GPCRs, there are multiple aspects of GPCR signaling that remain poorly understood. The reasons for this are the difficulties that are encountered upon studying GPCRs, for example, a poor solubility and low expression levels. In this work, we have managed to overcome some of these issues by developing an affinity-based probe for a prototypic GPCR, the adenosine A₁ receptor (A₁AR). Here, we show the design, synthesis, and biological evaluation of this probe in various biochemical assays, such as SDS-PAGE, confocal microscopy, and chemical proteomics.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A₁R/A₃R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A₁R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L250^6.51 to alanine increased the binding affinity of A17 at A₁R. We explored the structure-activity relationships against A₁R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

19F NMR: A promising tool for dynamic conformational studies of G protein-coupled receptors

Advances in X-ray crystallography and cryoelectron microscopy enabled unprecedented insights into the activation processes of G protein-coupled receptors (GPCRs). However, these static receptor structures provide limited information about dynamics and conformational transitions that play pivotal roles in mediating signaling diversity through the multifaceted interactions between ligands, receptors, and transducers. Developing NMR approaches to probe the dynamics of conformational transitions will push the frontier of receptor science toward a more comprehensive understanding of these signaling processes. Although much progress has been made during the last decades, it remains challenging to delineate receptor conformational states and interrogate the functions of the individual states at a quantitative level. Here we cover the progress of ¹⁹F NMR applications in GPCR conformational and dynamic studies during the past 20 years. Current challenges and limitations of ¹⁹F NMR for studying GPCR dynamics are also discussed, along with experimental strategies that will drive this field forward.

Purinergic receptors modulators: An emerging pharmacological tool for disease management

Purinergic signaling is mediated through extracellular nucleotides (adenosine 5'-triphosphate, uridine-5'-triphosphate, adenosine diphosphate, uridine-5'-diphosphate, and adenosine) that serve as signaling molecules. In the early 1990s, purines and pyrimidine receptors were cloned and characterized drawing the attention of scientists toward this aspect of cellular signaling. This signaling pathway is comprised of four subtypes of adenosine receptors (P1), eight subtypes of G-coupled protein receptors (P2YRs), and seven subtypes of ligand-gated ionotropic receptors (P2XRs). In current studies, the pathophysiology and therapeutic potentials of these receptors have been focused on. Various ligands, modulating the functions of purinergic receptors, are in current clinical practices for the treatment of various neurodegenerative disorders and cardiovascular diseases. Moreover, several purinergic receptors ligands are in advanced phases of clinical trials as a remedy for depression, epilepsy, autism, osteoporosis, atherosclerosis, myocardial infarction, diabetes, irritable bowel syndrome, and cancers. In the present study, agonists and antagonists of purinergic receptors have been summarized that may serve as pharmacological tools for drug design and development.

The Pharmacological Potential of Adenosine A2A Receptor Antagonists for Treating Parkinson's Disease

The adenosine A2A receptor subtype is recognized as a non-dopaminergic pharmacological target for the treatment of neurodegenerative disorders, notably Parkinson's disease (PD). The selective A2A receptor antagonist istradefylline is approved in the US and Japan as an adjunctive treatment to levodopa/decarboxylase inhibitors in adults with PD experiencing OFF episodes or a wearing-off phenomenon; however, the full potential of this drug class remains to be explored. In this article, we review the pharmacology of adenosine A2A receptor antagonists from the perspective of the treatment of both motor and non-motor symptoms of PD and their potential for disease modification.

International Union of Basic and Clinical Pharmacology. CXII: Adenosine Receptors: A Further Update

Our previous International Union of Basic and Clinical Pharmacology report on the nomenclature and classification of adenosine receptors (2011) contained a number of emerging developments with respect to this G protein-coupled receptor subfamily, including protein structure, protein oligomerization, protein diversity, and allosteric modulation by small molecules. Since then, a wealth of new data and results has been added, allowing us to explore novel concepts such as target binding kinetics and biased signaling of adenosine receptors, to examine a multitude of receptor structures and novel ligands, to gauge new pharmacology, and to evaluate clinical trials with adenosine receptor ligands. This review should therefore be considered a further update of our previous reports from 2001 and 2011. SIGNIFICANCE STATEMENT: Adenosine receptors (ARs) are of continuing interest for future treatment of chronic and acute disease conditions, including inflammatory diseases, neurodegenerative afflictions, and cancer. The design of AR agonists ("biased" or not) and antagonists is largely structure based now, thanks to the tremendous progress in AR structural biology. The A2A- and A2BAR appear to modulate the immune response in tumor biology. Many clinical trials for this indication are ongoing, whereas an A2AAR antagonist (istradefylline) has been approved as an anti-Parkinson agent.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Ribose and Non-Ribose A2A Adenosine Receptor Agonists: Do They Share the Same Receptor Recognition Mechanism?

Adenosine receptors have been a promising class of targets for the development of new therapies for several diseases. In recent years, a renewed interest in this field has risen, thanks to the implementation of a novel class of agonists that lack the ribose moiety, once considered essential for the agonistic profile. Recently, an X-ray crystal structure of the A_2A adenosine receptor has been solved, providing insights about the receptor activation from this novel class of agonists. Starting from this structural information, we have performed supervised molecular dynamics (SuMD) simulations to investigate the binding pathway of a non-nucleoside adenosine receptor agonist as well as one of three classic agonists. Furthermore, we analyzed the possible role of water molecules in receptor activation.

Recent advances in the synthesis of cyclic compounds using α,α-dicyanoolefins as versatile vinylogous nucleophiles

This article provides a review of the applications of α,α-dicyanoolefins as versatile vinylogous nucleophiles in the synthesis of various cyclic compounds, covering the literature from the past 13 years.

Deciphering conformational selectivity in the A2A adenosine G protein-coupled receptor by free energy simulations

Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A_2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A_2AAR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications.

The design of new ligands as chemical modulators of G protein-coupled receptors (GPCRs) has benefited considerably during the last years of advances in both the structural and computational biology disciplines. Within the last area, the use of free energy calculation methods has arisen as a computational tool to predict ligand affinities to explain structure-affinity relationships and guide lead optimization campaigns. However, our comprehension of the structural determinants of ligands with different pharmacological profile is scarce, and knowledge of the chemical modifications associated with an agonistic or antagonistic profile would be extremely valuable. We herein report an original implementation of the thermodynamic cycles associated with free energy perturbation (FEP) simulations, to mimic the conformational equilibrium between active and inactive GPCRs, and establish a framework to describe pharmacological profiles as a function of the ligands selectivity for a given receptor conformation. The advantage of this method resides into its simplicity of use, and the only consideration of active and inactive conformations of the receptor, with no simulation of the transitions between them. This model can accurately predict the pharmacological profile of series of full and partial agonists as opposed to antagonists of the A_2A adenosine receptor, and moreover, how certain mutations associated with modulation of basal activity can influence this pharmacological profiles, which enables our understanding of such clinically relevant mutations.