Extracellular Loop 2 of the Adenosine A1 Receptor Has a Key Role in Orthosteric Ligand Affinity and Agonist Efficacy

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.

Dynamic allosteric networks drive adenosine A1 receptor activation and G-protein coupling

G-protein coupled receptors (GPCRs) present specific activation pathways and signaling among receptor subtypes. Hence, an extensive knowledge of the structural dynamics of the receptor is critical for the development of therapeutics. Here, we target the adenosine A1 receptor (A1R), for which a negligible number of drugs have been approved. We combine molecular dynamics simulations, enhanced sampling techniques, network theory and pocket detection to decipher the activation pathway of A1R, decode the allosteric networks and identify transient pockets. The A1R activation pathway reveal hidden intermediate and pre-active states together with the inactive and fully-active states observed experimentally. The protein energy networks computed throughout these conformational states successfully unravel the extra and intracellular allosteric centers and the communication pathways that couples them. We observe that the allosteric networks are dynamic, being increased along activation and fine-tuned in presence of the trimeric G-proteins. Overlap of transient pockets and energy networks uncover how the allosteric coupling between pockets and distinct functional regions of the receptor is altered along activation. By an in-depth analysis of the bridge between activation pathway, energy networks and transient pockets, we provide a further understanding of A1R. This information can be useful to ease the design of allosteric modulators for A1R.

Structure-based virtual screening discovers potent and selective adenosine A1 receptor antagonists

Development of subtype-selective leads is essential in drug discovery campaigns targeting G protein-coupled receptors (GPCRs). Herein, a structure-based virtual screening approach to rationally design subtype-selective ligands was applied to the A₁ and A_2A adenosine receptors (A₁R and A_2AR). Crystal structures of these closely related subtypes revealed a non-conserved subpocket in the binding sites that could be exploited to identify A₁R selective ligands. A library of 4.6 million compounds was screened computationally against both receptors using molecular docking and 20 A₁R selective ligands were predicted. Of these, seven antagonized the A₁R with micromolar activities and several compounds displayed slight selectivity for this subtype. Twenty-seven analogs of two discovered scaffolds were designed, resulting in antagonists with nanomolar potency and up to 76-fold A₁R-selectivity. Our results show the potential of structure-based virtual screening to guide discovery and optimization of subtype-selective ligands, which could facilitate the development of safer drugs.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~ 1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

The full activation mechanism of the adenosine A1 receptor revealed by GaMD and Su-GaMD simulations

The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A₁ receptor (A₁R), we successfully revealed the full activation mechanism of A₁R, including adenosine (Ado)-A₁R recognition, preactivation of A₁R, and A₁R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A₁R initiates conformational changes and the preactivation of A₁R. In turn, the binding of G_i2 to the intracellular side of A₁R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A₁R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A₁R in the future.

Cancer-Related Somatic Mutations in Transmembrane Helices Alter Adenosine A1 Receptor Pharmacology

Overexpression of the adenosine A₁ receptor (A₁AR) has been detected in various cancer cell lines. However, the role of A₁AR in tumor development is still unclear. Thirteen A₁AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. We have investigated the pharmacology of the mutations located at the 7-transmembrane domain using a yeast system. Concentration-growth curves were obtained with the full agonist CPA and compared to the wild type hA₁AR. H78L^3.23 and S246T^6.47 showed increased constitutive activity, while only the constitutive activity of S246T^6.47 could be reduced to wild type levels by the inverse agonist DPCPX. Decreased constitutive activity was observed on five mutant receptors, among which A52V^2.47 and W188C^5.46 showed a diminished potency for CPA. Lastly, a complete loss of activation was observed in five mutant receptors. A selection of mutations was also investigated in a mammalian system, showing comparable effects on receptor activation as in the yeast system, except for residues pointing toward the membrane. Taken together, this study will enrich the view of the receptor structure and function of A₁AR, enlightening the consequences of these mutations in cancer. Ultimately, this may provide an opportunity for precision medicine for cancer patients with pathological phenotypes involving these mutations.

Cancer-related somatic mutations alter adenosine A1 receptor pharmacology-A focus on mutations in the loops and C-terminus

G protein-coupled receptors (GPCRs) are known to be involved in tumor progression and metastasis. The adenosine A₁ receptor (A₁ AR) has been detected to be over-expressed in various cancer cell lines. However, the role of A₁ AR in tumor development is not yet well characterized. A series of A₁ AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. In this study, we have investigated the pharmacology of mutations located outside of the 7-transmembrane domain by using a "single-GPCR-one-G protein" yeast system. Concentration-growth curves were obtained with the full agonist CPA for 12 mutant receptors and compared to the wild-type hA₁ AR. Most mutations located at the extracellular loops (EL) reduced the levels of constitutive activity of the receptor and agonist potency. For mutants at the intracellular loops (ILs) of the receptor, an increased constitutive activity was found for mutant receptor L211R^5.69 , while a decreased constitutive activity and agonist response were found for mutant receptor L113F^34.51 . Lastly, mutations identified on the C-terminus did not significantly influence the pharmacological function of the receptor. A selection of mutations was also investigated in a mammalian system. Overall, similar effects on receptor activation compared to the yeast system were found with mutations located at the EL, but some contradictory effects were observed for mutations located at the IL. Taken together, this study will enrich the insight of A₁ AR structure and function, enlightening the consequences of these mutations in cancer. Ultimately, this may provide potential precision medicine in cancer treatment.

GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling

We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A₁ receptor (A₁AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A₁AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.

Positive allosteric mechanisms of adenosine A1 receptor-mediated analgesia

The adenosine A₁ receptor (A₁R) is a promising therapeutic target for non-opioid analgesic agents to treat neuropathic pain^1,2. However, development of analgesic orthosteric A₁R agonists has failed because of a lack of sufficient on-target selectivity as well as off-tissue adverse effects³. Here we show that [2-amino-4-(3,5-bis(trifluoromethyl)phenyl)thiophen-3-yl)(4-chlorophenyl)methanone] (MIPS521), a positive allosteric modulator of the A₁R, exhibits analgesic efficacy in rats in vivo through modulation of the increased levels of endogenous adenosine that occur in the spinal cord of rats with neuropathic pain. We also report the structure of the A₁R co-bound to adenosine, MIPS521 and a G_i2 heterotrimer, revealing an extrahelical lipid-detergent-facing allosteric binding pocket that involves transmembrane helixes 1, 6 and 7. Molecular dynamics simulations and ligand kinetic binding experiments support a mechanism whereby MIPS521 stabilizes the adenosine-receptor-G protein complex. This study provides proof of concept for structure-based allosteric drug design of non-opioid analgesic agents that are specific to disease contexts.

Multisite Model of Allosterism for the Adenosine A1 Receptor

Despite being a target for about one-third of approved drugs, G protein-coupled receptors (GPCRs) still represent a tremendous reservoir for therapeutic strategies against diseases. For example, several cardiovascular and central nervous system conditions could benefit from clinical agents that activate the adenosine 1 receptor (A₁R); however, the pursuit of A₁R agonists for clinical use is usually impeded by both on- and off-target side effects. One of the possible strategies to overcome this issue is the development of positive allosteric modulators (PAMs) capable of selectively enhancing the effect of a specific receptor subtype and triggering functional selectivity (a phenomenon also referred to as bias). Intriguingly, besides enforcing the effect of agonists upon binding to an allosteric site, most of the A₁R PAMs display intrinsic partial agonism and orthosteric competition with antagonists. To rationalize this behavior, we simulated the binding of the prototypical PAMs PD81723 and VCP171, the full-agonist NECA, the antagonist 13B, and the bitopic agonist VCP746. We propose that a single PAM can bind several A₁R sites rather than a unique allosteric pocket, reconciling the structure-activity relationship and the mutagenesis results.

Development and Application of Subtype-Selective Fluorescent Antagonists for the Study of the Human Adenosine A1 Receptor in Living Cells

The adenosine A₁ receptor (A₁AR) is a G-protein-coupled receptor (GPCR) that provides important therapeutic opportunities for a number of conditions including congestive heart failure, tachycardia, and neuropathic pain. The development of A₁AR-selective fluorescent ligands will enhance our understanding of the subcellular mechanisms underlying A₁AR pharmacology facilitating the development of more efficacious and selective therapies. Herein, we report the design, synthesis, and application of a novel series of A₁AR-selective fluorescent probes based on 8-functionalized bicyclo[2.2.2]octylxanthine and 3-functionalized 8-(adamant-1-yl) xanthine scaffolds. These fluorescent conjugates allowed quantification of kinetic and equilibrium ligand binding parameters using NanoBRET and visualization of specific receptor distribution patterns in living cells by confocal imaging and total internal reflection fluorescence (TIRF) microscopy. As such, the novel A₁AR-selective fluorescent antagonists described herein can be applied in conjunction with a series of fluorescence-based techniques to foster understanding of A₁AR molecular pharmacology and signaling in living cells.

Deciphering the Agonist Binding Mechanism to the Adenosine A1 Receptor

Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A1R) was determined in the active, Gi protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A1R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N6-cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A1R agonists.

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.

New Insights into Key Determinants for Adenosine 1 Receptor Antagonists Selectivity Using Supervised Molecular Dynamics Simulations

Adenosine receptors (ARs), like many otherGprotein-coupledreceptors (GPCRs), are targets of primary interest indrug design. However, one of the main limits for the development of drugs for this class of GPCRs is the complex selectivity profile usually displayed by ligands. Numerous efforts have been madefor clarifying the selectivity of ARs, leading to the development of many ligand-based models. The structure of the AR subtype A₁ (A₁AR) has been recently solved, providing important structural insights. In the present work, we rationalized the selectivity profile of two selective A₁AR and A_2AAR antagonists, investigating their recognition trajectories obtained by Supervised Molecular Dynamics from an unbound state and monitoring the role of the water molecules in the binding site.

Medicinal Chemistry and Therapeutic Potential of Agonists, Antagonists and Allosteric Modulators of A1 Adenosine Receptor: Current Status and Perspectives

Adenosine is a purine nucleoside, responsible for the regulation of a wide range of physiological and pathophysiological conditions by binding with four G-protein-coupled receptors (GPCRs), namely A1, A2A, A2B and A3 adenosine receptors (ARs). In particular, A1 AR is ubiquitously present, mediating a variety of physiological processes throughout the body, thus represents a promising drug target for the management of various pathological conditions. Agonists of A1 AR are found to be useful for the treatment of atrial arrhythmia, angina, type-2 diabetes, glaucoma, neuropathic pain, epilepsy, depression and Huntington's disease, whereas antagonists are being investigated for the treatment of diuresis, congestive heart failure, asthma, COPD, anxiety and dementia. However, treatment with full A1 AR agonists has been associated with numerous challenges like cardiovascular side effects, off-target activation as well as desensitization of A1 AR leading to tachyphylaxis. In this regard, partial agonists of A1 AR have been found to be beneficial in enhancing insulin sensitivity and subsequently reducing blood glucose level, while avoiding severe CVS side effects and tachyphylaxis. Allosteric enhancer of A1 AR is found to be potent for the treatment of neuropathic pain, culminating the side effects related to off-target tissue activation of A1 AR. This review provides an overview of the medicinal chemistry and therapeutic potential of various agonists/partial agonists, antagonists and allosteric modulators of A1 AR, with a particular emphasis on their current status and future perspectives in clinical settings.

Retrospective ensemble docking of allosteric modulators in an adenosine G-protein-coupled receptor

BACKGROUND

Ensemble docking has proven useful in drug discovery and development. It increases the hit rate by incorporating receptor flexibility into molecular docking as demonstrated on important drug targets including G-protein-coupled receptors (GPCRs). Adenosine A₁ receptor (A₁AR) is a key GPCR that has been targeted for treating cardiac ischemia-reperfusion injuries, neuropathic pain and renal diseases. Development of allosteric modulators, compounds binding to distinct and less conserved GPCR target sites compared with agonists and antagonists, has attracted increasing interest for designing selective drugs of the A₁AR. Despite significant advances, more effective approaches are needed to discover potent and selective allosteric modulators of the A₁AR.

METHODS

Ensemble docking that integrates Gaussian accelerated molecular dynamic (GaMD) simulations and molecular docking using Autodock has been implemented for retrospective docking of known positive allosteric modulators (PAMs) in the A₁AR.

RESULTS

Ensemble docking outperforms docking of the receptor cryo-EM structure. The calculated docking enrichment factors (EFs) and the area under the receiver operating characteristic curves (AUC) are significantly increased.

CONCLUSIONS

Receptor ensembles generated from GaMD simulations are able to increase the success rate of discovering PAMs of A₁AR. It is important to account for receptor flexibility through GaMD simulations and flexible docking.

GENERAL SIGNIFICANCE

Ensemble docking is a promising approach for drug discovery targeting flexible receptors.

Structure of the agonist 12-HHT in its BLT2 receptor-bound state

G Protein-Coupled receptors represent the main communicating pathway for signals from the outside to the inside of most of eukaryotic cells. They define the largest family of integral membrane receptors at the surface of the cells and constitute the main target of the current drugs on the market. The low affinity leukotriene receptor BLT2 is a receptor involved in pro- and anti-inflammatory pathways and can be activated by various unsaturated fatty acid compounds. We present here the NMR structure of the agonist 12-HHT in its BLT2-bound state and a model of interaction of the ligand with the receptor based on a conformational homology modeling associated with docking simulations. Put into perspective with the data obtained with leukotriene B4, our results illuminate the ligand selectivity of BLT2 and may help define new molecules to modulate the activity of this receptor.

Allosteric Modulation of Cannabinoid Receptor 1-Current Challenges and Future Opportunities

The cannabinoid receptor type 1 (CB1R), a G protein-coupled receptor (GPCR), plays an essential role in the control of many physiological processes such as hunger, memory loss, gastrointestinal activity, catalepsy, fear, depression, and chronic pain. Therefore, it is an attractive target for drug discovery to manage pain, neurodegenerative disorders, obesity, and substance abuse. However, the psychoactive adverse effects, generated by CB1R activation in the brain, limit the use of the orthosteric CB1R ligands as drugs. The discovery of CB1R allosteric modulators during the last decade provided new tools to target the CB1R. Moreover, application of the site-directed mutagenesis in combination with advanced physical methods, especially X-ray crystallography and computational modeling, has opened new horizons for understanding the complexity of the structure, function, and activity of cannabinoid receptors. In this paper, we present the latest advances in research on the CB1R, its allosteric modulation and allosteric ligands, and their translational potential. We focused on structural essentials of the cannabinoid 1 receptor- ligand (drug) interactions, as well as modes of CB1R signaling regulation.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

Mechanistic Insights into Specific G Protein Interactions with Adenosine Receptors

Coupling between G-protein-coupled receptors (GPCRs) and the G proteins is a key step in cellular signaling. Despite extensive experimental and computational studies, the mechanism of specific GPCR-G protein coupling remains poorly understood. This has greatly hindered effective drug design of GPCRs that are primary targets of ∼¹/₃ of currently marketed drugs. Here, we have employed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method to decipher the mechanism of the GPCR-G protein interactions. Adenosine receptors (ARs) were used as model systems based on very recently determined cryo-EM structures of the A₁AR and A_2AAR coupled with the G_i and G_s proteins, respectively. Changing the G_i protein to the G_s led to increased fluctuations in the A₁AR and agonist adenosine (ADO), while agonist 5'-N-ethylcarboxamidoadenosine (NECA) binding in the A_2AAR could be still stabilized upon changing the G_s protein to the G_i. Free energy calculations identified one stable low-energy conformation for each of the A₁AR-G_i and A_2AAR-G_s complexes as in the cryo-EM structures, similarly for the A_2AAR-G_i complex. In contrast, the ADO agonist and G_s protein sampled multiple conformations in the A₁AR-G_s system. GaMD simulations thus indicated that the A₁AR preferred to couple with the G_i protein to the G_s, while the A_2AAR could couple with both the G_s and G_i proteins, being highly consistent with experimental findings of the ARs. More importantly, detailed analysis of the atomic simulations showed that the specific AR-G protein coupling resulted from remarkably complementary residue interactions at the protein interface, involving mainly the receptor transmembrane 6 helix and the G_α α5 helix and α4-β6 loop. In summary, the GaMD simulations have provided unprecedented insights into the dynamic mechanism of specific GPCR-G protein interactions at an atomistic level.

Molecular Dynamics Simulations of Adenosine Receptors: Advances, Applications and Trends

:

Adenosine receptors (ARs) are transmembrane proteins that belong to the G protein-coupled receptors (GPCRs) superfamily and mediate the biological functions of adenosine. To date, four AR subtypes are known, namely A1, A2A, A2B and A3 that exhibit different signaling pathways, tissue localization, and mechanisms of activation. Moreover, the widespread ARs and their implication in numerous physiological and pathophysiological conditions had made them pivotal therapeutic targets for developing clinically effective agents.

:

The crystallographic success in identifying the 3D crystal structures of A2A and A1 ARs has dramatically enriched our understanding of their structural and functional properties such as ligand binding and signal transduction. This, in turn, has provided a structural basis for a larger contribution of computational methods, particularly molecular dynamics (MD) simulations, toward further investigation of their molecular properties and designing bioactive ligands with therapeutic potential. MD simulation has been proved to be an invaluable tool in investigating ARs and providing answers to some critical questions. For example, MD has been applied in studying ARs in terms of ligand-receptor interactions, molecular recognition, allosteric modulations, dimerization, and mechanisms of activation, collectively aiding in the design of subtype selective ligands.

:

In this review, we focused on the advances and different applications of MD simulations utilized to study the structural and functional aspects of ARs that can foster the structure-based design of drug candidates. In addition, relevant literature was briefly discussed which establishes a starting point for future advances in the field of drug discovery to this pivotal group of drug targets.

Probe dependence of allosteric enhancers on the binding affinity of adenosine A1 -receptor agonists at rat and human A1 -receptors measured using NanoBRET

BACKGROUND AND PURPOSE

Adenosine is a local mediator that regulates a number of physiological and pathological processes via activation of adenosine A₁ -receptors. The activity of adenosine can be regulated at the level of its target receptor via drugs that bind to an allosteric site on the A₁ -receptor. Here, we have investigated the species and probe dependence of two allosteric modulators on the binding characteristics of fluorescent and nonfluorescent A₁ -receptor agonists.

EXPERIMENTAL APPROACH

A Nano-luciferase (Nluc) BRET (NanoBRET) methodology was used. This used N-terminal Nluc-tagged A₁ -receptors expressed in HEK293T cells in conjunction with both fluorescent A₁ -receptor agonists (adenosine and NECA analogues) and a fluorescent antagonist CA200645.

KEY RESULTS

PD 81,723 and VCP171 elicited positive allosteric effects on the binding affinity of orthosteric agonists at both the rat and human A₁ -receptors that showed clear probe dependence. Thus, the allosteric effect on the highly selective partial agonist capadenoson was much less marked than for the full agonists NECA, adenosine, and CCPA in both species. VCP171 and, to a lesser extent, PD 81,723, also increased the specific binding of three fluorescent A₁ -receptor agonists in a species-dependent manner that involved increases in B_max and pK_D .

CONCLUSIONS AND IMPLICATIONS

These results demonstrate the power of the NanoBRET ligand-binding approach to study the effect of allosteric ligands on the binding of fluorescent agonists to the adenosine A₁ -receptor in intact living cells. Furthermore, our studies suggest that VCP171 and PD 81,723 may switch a proportion of A₁ -receptors to an active agonist conformation (R*).

Utility of an "Allosteric Site-Impaired" M2 Muscarinic Acetylcholine Receptor as a Novel Construct for Validating Mechanisms of Action of Synthetic and Putative Endogenous Allosteric Modulators

Muscarinic acetylcholine receptors (mAChRs) are exemplar models for understanding G protein-coupled receptor (GPCR) allostery, possessing a "common" allosteric site in an extracellular vestibule (ECV) for synthetic modulators including gallamine, strychnine, and brucine. In addition, there is intriguing evidence of endogenous peptides/proteins that may target this region at the M₂ mAChR. A common feature of synthetic and endogenous M₂ mAChR negative allosteric modulators (NAMs) is their cationic nature. Using a structure-based approach, we previously designed a mutant M₂ mAChR (N410K+T423K) to specifically abrogate binding of ECV cationic modulators (Dror et al., 2013). Herein, we used this "allosteric site-impaired" receptor to investigate allosteric interactions of synthetic modulators as well as basic peptides (poly-l-arginine, endogenously produced protamine, and major basic protein). Using [³H]N-methylscopolamine equilibrium and kinetic binding and functional assays of guanosine 5'-O-[γ-thio]triphosphate [³⁵S] binding and extracellular signal-regulated kinases 1 and 2 phosphorylation, we found modest effects of the mutations on potencies of orthosteric antagonists and an increase in the affinity of the cognate agonist, acetylcholine, likely reflecting the effect of the mutations on the access/egress of these ligands into the orthosteric pocket. More importantly, we noted a significant abrogation in affinity for all synthetic or peptidic modulators at the mutant mAChR, validating their allosteric nature. Collectively, these findings provide evidence for a hitherto-unappreciated role of endogenous cationic peptides interacting allosterically at the M₂ mAChR and identify the allosteric site-impaired GPCR as a tool for validating NAM activity as well as a potential candidate for future chemogenetic strategies to understand the physiology of endogenous allosteric substances.

New paradigms in adenosine receptor pharmacology: allostery, oligomerization and biased agonism

Adenosine receptors are a family of GPCRs containing four subtypes (A₁ , A_2A , A_2B and A₃ receptors), all of which bind the ubiquitous nucleoside adenosine. These receptors play an important role in physiology and pathophysiology and therefore represent attractive drug targets for a range of conditions. The theoretical framework surrounding drug action at adenosine receptors now extends beyond the notion of prototypical agonism and antagonism to encompass more complex pharmacological concepts. New paradigms include allostery, in which ligands bind a topographically distinct receptor site from that of the endogenous agonist, homomeric or heteromeric interactions across receptor oligomers and biased agonism, that is, ligand-dependent differential intracellular signalling. This review provides a concise overview of allostery, oligomerization and biased agonism at adenosine receptors and outlines how these paradigms may enhance future drug discovery endeavours focussed on the development of novel therapeutic agents acting at adenosine receptors.

LINKED ARTICLES

This article is part of a themed section on Molecular Pharmacology of GPCRs. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.21/issuetoc.

Identification of avian vasotocin receptor subtype-specific antagonists involved in the stress response of the chicken, Gallus gallus

Vasotocin 1a and 1b receptors (V1aR and V1bR) have been shown to play important roles in the neuroendocrine regulation of stress responses via the anterior pituitary (AP) of birds. To identify effective subtype-specific antagonists for the chicken V1aR (cV1aR) and cV1bR, potential antagonists to the mammalian V1R were screened against the cV1aR and cV1bR 3D structural models by molecular docking analysis with determination of binding pocket/amino acid residues involved in the interaction. The antagonistic effects of the selected ligands were examined by measuring pro-opiomelanocortin (POMC) heteronuclear RNA (hnPOMC) levels following the in vitro stress administration to primary chicken AP cells. Results of in silico analysis showed that the Manning compound and several other antagonists were bound to cV1bR with higher affinity than the natural agonist, arginine vasotocin (AVT). Similarities and differences in the antagonist-receptor binding interface with receptors were characterized for each ligand. Non-peptide mammalian V1bR antagonists, SSR-149415 and L-368899, were shown to be effective and had an additive effect in blocking POMC hnRNA expression in pituitary cell culture studies. SR-49059 antagonized the effect(s) of AVT/CRH on the downregulation of the cV1aR and the upregulation of the cCRH-R2 expression but not the cV1bR and cCRH-R1. The Manning compound antagonized the downregulation of cV1aR, cV1bR and cCRH-R1 and the upregulation of cCRH-R2 expression. The specificity of antagonists apparently resulted from unique differences in the interacting residues and their binding affinities. Collectively, these results provide valuable leads for future development of novel compounds capable of blocking or attenuating the AP stress response of avian species and perhaps other non-mammalian vertebrates as well.

Computer-aided design of multi-target ligands at A1R, A2AR and PDE10A, key proteins in neurodegenerative diseases

Compounds designed to display polypharmacology may have utility in treating complex diseases, where activity at multiple targets is required to produce a clinical effect. In particular, suitable compounds may be useful in treating neurodegenerative diseases by promoting neuronal survival in a synergistic manner via their multi-target activity at the adenosine A1 and A2A receptors (A1R and A2AR) and phosphodiesterase 10A (PDE10A), which modulate intracellular cAMP levels. Hence, in this work we describe a computational method for the design of synthetically feasible ligands that bind to A1 and A2A receptors and inhibit phosphodiesterase 10A (PDE10A), involving a retrosynthetic approach employing in silico target prediction and docking, which may be generally applicable to multi-target compound design at several target classes. This approach has identified 2-aminopyridine-3-carbonitriles as the first multi-target ligands at A1R, A2AR and PDE10A, by showing agreement between the ligand and structure based predictions at these targets. The series were synthesized via an efficient one-pot scheme and validated pharmacologically as A1R/A2AR-PDE10A ligands, with IC50 values of 2.4-10.0 μM at PDE10A and Ki values of 34-294 nM at A1R and/or A2AR. Furthermore, selectivity profiling of the synthesized 2-amino-pyridin-3-carbonitriles against other subtypes of both protein families showed that the multi-target ligand 8 exhibited a minimum of twofold selectivity over all tested off-targets. In addition, both compounds 8 and 16 exhibited the desired multi-target profile, which could be considered for further functional efficacy assessment, analog modification for the improvement of selectivity towards A1R, A2AR and PDE10A collectively, and evaluation of their potential synergy in modulating cAMP levels.

Synergistic Use of GPCR Modeling and SDM Experiments to Understand Ligand Binding

There is a substantial amount of historical ligand binding data available from site-directed mutagenesis (SDM) studies of many different GPCR subtypes. This information was generated prior to the wave of GPCR crystal structure, in an effort to understand ligand binding with a view to drug discovery. Concerted efforts to determine the atomic structure of GPCRs have proven extremely successful and there are now more than 80 GPCR crystal structure in the PDB database, many of which have been obtained in the presence of receptor ligands and associated G proteins. These structural data enable the generation of computational model structures for all GPCRs, including those for which crystal structures do not yet exist. The power of these models in designing novel ligands, especially those with improved residence times, and for better understanding receptor function can be enhanced tremendously by combining them synergistically with historic SDM ligand binding data. Here, we describe a protocol by which historic SDM binding data and receptor models may be used together to identify novel key residues for mutagenesis studies.

Impact of protein-ligand solvation and desolvation on transition state thermodynamic properties of adenosine A2A ligand binding kinetics

Ligand-protein binding kinetic rates are growing in importance as parameters to consider in drug discovery and lead optimization. In this study we analysed using surface plasmon resonance (SPR) the transition state (TS) properties of a set of six adenosine A2A receptor inhibitors, belonging to both the xanthine and the triazolo-triazine scaffolds. SPR highlighted interesting differences among the ligands in the enthalpic and entropic components of the TS energy barriers for the binding and unbinding events. To better understand at a molecular level these differences, we developed suMetaD, a novel molecular dynamics (MD)-based approach combining supervised MD and metadynamics. This method allows simulation of the ligand unbinding and binding events. It also provides the system conformation corresponding to the highest energy barrier the ligand is required to overcome to reach the final state. For the six ligands evaluated in this study their TS thermodynamic properties were linked in particular to the role of water molecules in solvating/desolvating the pocket and the small molecules. suMetaD identified kinetic bottleneck conformations near the bound state position or in the vestibule area. In the first case the barrier is mainly enthalpic, requiring the breaking of strong interactions with the protein. In the vestibule TS location the kinetic bottleneck is instead mainly of entropic nature, linked to the solvent behaviour.

Structure-Based Design of Potent and Selective Ligands at the Four Adenosine Receptors

The four receptors that signal for adenosine, A₁, A_2A, A_2B and A₃ ARs, belong to the superfamily of G protein-coupled receptors (GPCRs). They mediate a number of (patho)physiological functions and have attracted the interest of the biopharmaceutical sector for decades as potential drug targets. The many crystal structures of the A_2A, and lately the A₁ ARs, allow for the use of advanced computational, structure-based ligand design methodologies. Over the last decade, we have assessed the efficient synthesis of novel ligands specifically addressed to each of the four ARs. We herein review and update the results of this program with particular focus on molecular dynamics (MD) and free energy perturbation (FEP) protocols. The first in silico mutagenesis on the A₁AR here reported allows understanding the specificity and high affinity of the xanthine-antagonist 8-Cyclopentyl-1,3-dipropylxanthine (DPCPX). On the A_2AAR, we demonstrate how FEP simulations can distinguish the conformational selectivity of a recent series of partial agonists. These novel results are complemented with the revision of the first series of enantiospecific antagonists on the A_2BAR, and the use of FEP as a tool for bioisosteric design on the A₃AR.

A monoclonal antibody raised against a thermo-stabilised β1-adrenoceptor interacts with extracellular loop 2 and acts as a negative allosteric modulator of a sub-set of β1-adrenoceptors expressed in stable cell lines

Recent interest has focused on antibodies that can discriminate between different receptor conformations. Here we have characterised the effect of a monoclonal antibody (mAb3), raised against a purified thermo-stabilised turkey β₁-adrenoceptor (β₁AR-m23 StaR), on β₁-ARs expressed in CHO-K1 or HEK 293 cells. Immunohistochemical and radioligand-binding studies demonstrated that mAb3 was able to bind to ECL2 of the tβ₁-AR, but not its human homologue. Specific binding of mAb3 to tβ₁-AR was inhibited by a peptide based on the turkey, but not the human, ECL2 sequence. Studies with [³H]-CGP 12177 demonstrated that mAb3 prevented the binding of orthosteric ligands to a subset (circa 40%) of turkey β₁-receptors expressed in both CHO K1 and HEK 293 cells. MAb3 significantly reduced the maximum specific binding capacity of [³H]-CGP-12177 without influencing its binding affinity. Substitution of ECL2 of tβ₁-AR with its human equivalent, or mutation of residues D186S, P187D, Q188E prevented the inhibition of [³H]-CGP 12177 binding by mAb3. MAb3 also elicited a negative allosteric effect on agonist-stimulated cAMP responses. The identity of the subset of turkey β₁-adrenoceptors influenced by mAb3 remains to be established but mAb3 should become an important tool to investigate the nature of β₁-AR conformational states and oligomeric complexes.

Molecular Recognition of Agonists and Antagonists by the Nucleotide-Activated G Protein-Coupled P2Y2 Receptor

A homology model of the nucleotide-activated P2Y₂R was created based on the X-ray structures of the P2Y₁ receptor. Docking studies were performed, and receptor mutants were created to probe the identified binding interactions. Mutation of residues predicted to interact with the ribose (Arg110) and the phosphates of the nucleotide agonists (Arg265, Arg292) or that contribute indirectly to binding (Tyr288) abolished activity. The Y114F, R194A, and F261A mutations led to inactivity of diadenosine tetraphosphate and to a reduced response of UTP. Significant reduction in agonist potency was observed for all other receptor mutants (Phe111, His184, Ser193, Phe261, Tyr268, Tyr269) predicted to be involved in agonist recognition. An ionic lock between Asp185 and Arg292 that is probably involved in receptor activation interacts with the phosphate groups. The antagonist AR-C118925 and anthraquinones likely bind to the orthosteric site. The updated homology models will be useful for virtual screening and drug design.

Serotonin 2A receptor disulfide bridge integrity is crucial for ligand binding to different signalling states but not for its homodimerization

The serotonin 2_A (5-HT_2A) receptor is a G-protein coupled receptor (GPCR) with a conserved disulfide bridge formed by Cys¹⁴⁸ (transmembrane helix 3, TM3) and Cys²²⁷ (extracellular loop 2, ECL-2). We hypothesized that disulfide bridges may determine serotonin 5-HT_2A receptor functions such as receptor activation, functional selectivity and ligand recognition. We used the reducing agent dithiothreitol (DTT) to determine how the reduction of disulfide bridges affects radioligand binding, second messenger mobilization and receptor dimerization. A DTT-induced decrease in the number of binding sites (1190 ± 63.55 fmol/mg protein for control cells compared with 921.2 ± 60.84 fmol/mg protein for DTT-treated cells) as well as in the efficacy of both signalling pathways characterized was observed, although the affinity and potency were unchanged. Bioluminiscence resonance energy transfer (BRET) assays revealed the DTT treatment did not modify the homodimeric nature of serotonin 5-HT_2A receptors. In molecular dynamic simulations, the ECL-2 of the receptor with a broken cysteine bond adopts a wider variety of conformations, some of which protrude deeper into the receptor orthosteric binding pocket leading to collapse of the pocket. A shrunken binding pocket would be incapable of accommodating lysergic acid diethylamide (LSD). Our findings suggest that the decrease of efficacy may be due to disruption of disulfide bridge between TM3 and ECL-2. This reveals the integrity of the ECL-2 epitope, which should be explored in the development of novel ligands acting as allosteric modulators of serotonin 5-HT_2A receptors.

Priming GPCR signaling through the synergistic effect of two G proteins

Although individual G-protein-coupled receptors (GPCRs) are known to activate one or more G proteins, the GPCR-G-protein interaction is viewed as a bimolecular event involving the formation of a ternary ligand-GPCR-G-protein complex. Here, we present evidence that individual GPCR-G-protein interactions can reinforce each other to enhance signaling through canonical downstream second messengers, a phenomenon we term "GPCR priming." Specifically, we find that the presence of noncognate Gq protein enhances cAMP stimulated by two Gs-coupled receptors, β2-adrenergic receptor (β2-AR) and D₁ dopamine receptor (D₁-R). Reciprocally, Gs enhances IP₁ through vasopressin receptor (V_1A-R) but not α1 adrenergic receptor (α1-AR), suggesting that GPCR priming is a receptor-specific phenomenon. The C terminus of either the Gαs or Gαq subunit is sufficient to enhance Gα subunit activation and cAMP levels. Interaction of Gαs or Gαq C termini with the GPCR increases signaling potency, suggesting an altered GPCR conformation as the underlying basis for GPCR priming. We propose three parallel mechanisms involving (i) sequential G-protein interactions at the cognate site, (ii) G-protein interactions at distinct allosteric and cognate sites on the GPCR, and (iii) asymmetric GPCR dimers. GPCR priming suggests another layer of regulation in the classic GPCR ternary-complex model, with broad implications for the multiplicity inherent in signaling networks.

Role of the Second Extracellular Loop of the Adenosine A1 Receptor on Allosteric Modulator Binding, Signaling, and Cooperativity

Allosteric modulation of adenosine A1 receptors (A1ARs) offers a novel therapeutic approach for the treatment of numerous central and peripheral disorders; however, despite decades of research, there is a relative paucity of structural information regarding the A1AR allosteric site and mechanisms governing cooperativity with orthosteric ligands. We combined alanine-scanning mutagenesis of the A1AR second extracellular loop (ECL2) with radioligand binding and functional interaction assays to quantify effects on allosteric ligand affinity, cooperativity, and efficacy. Docking and molecular dynamics (MD) simulations were performed using an A1AR homology model based on an agonist-bound A2AAR structure. Substitution of E172ECL2 for alanine reduced the affinity of the allosteric modulators PD81723 and VCP171 for the unoccupied A1AR. Residues involved in cooperativity with the orthosteric agonist NECA were different in PD81723 and VCP171; positive cooperativity between PD81723 and NECA was reduced on alanine substitution of a number of ECL2 residues, including E170ECL2 and K173ECL2, whereas mutation of W146ECL2 and W156ECL2 decreased VCP171 cooperativity with NECA. Molecular modeling localized a likely allosteric pocket for both modulators to an extracellular vestibule that overlaps with a region used by orthosteric ligands as they transit into the canonical A1AR orthosteric site. MD simulations confirmed a key interaction between E172ECL2 and both modulators. Bound PD81723 is flanked by another residue, E170ECL2, which forms hydrogen bonds with adjacent K168ECL2 and K173ECL2. Collectively, our data suggest E172ECL2 is a key allosteric ligand-binding determinant, whereas hydrogen-bonding networks within the extracellular vestibule may facilitate the transmission of cooperativity between orthosteric and allosteric sites.