Antagonist binding and induced conformational dynamics of GPCR A2A adenosine receptor

Selective Deuteration Improves the Affinity of Adenosine A2A Receptor Ligands: A Computational Case Study with Istradefylline and Caffeine

We used a range of computational techniques to assess the effect of selective C-H deuteration on the antagonist istradefylline affinity for the adenosine A_2A receptor, which was discussed relative to its structural analogue caffeine, a well-known and likely the most widely used stimulant. The obtained results revealed that smaller caffeine shows high receptor flexibility and exchanges between two distinct poses, which agrees with crystallographic data. In contrast, the additional C8-trans-styryl fragment in istradefylline locks the ligand within a uniform binding pose, while contributing to the affinity through the C-H···π and π···π contacts with surface residues, which, together with its much lower hydration prior to binding, enhances the affinity over caffeine. In addition, the aromatic C8-unit shows a higher deuteration sensitivity over the xanthine part, so when both of its methoxy groups are d₆-deuterated, the affinity improvement is -0.4 kcal mol^-1, which surpasses the overall affinity gain of -0.3 kcal mol^-1 in the perdeuterated d₉-caffeine. Yet, the latter predicts around 1.7-fold potency increase, being relevant for its pharmaceutical implementations, and also those within the coffee and energy drink production industries. Still, the full potential of our strategy is achieved in polydeuterated d₁₉-istradefylline, whose A_2A affinity improves by -0.6 kcal mol^-1, signifying a 2.8-fold potency increase that strongly promotes it as a potential synthetic target. This knowledge supports deuterium application in drug design, and while the literature already reports about over 20 deuterated drugs currently in the clinical development, it is easily foreseen that more examples will hit the market in the years to come. With this in mind, we propose that the devised computational methodology, involving the ONIOM division of the QM region for the ligand and the MM region for its environment, with an implicit quantization of nuclear motions relevant for the H/D exchange, allows fast and efficient estimates of the binding isotope effects in any biological system.

Single Stabilizing Point Mutation Enables High-Resolution Co-Crystal Structures of the Adenosine A2A Receptor with Preladenant Conjugates

The G protein-coupled adenosine A_2A receptor (A_2A AR) is an important new (potential) drug target in immuno-oncology, and for neurodegenerative diseases. Preladenant and its derivatives belong to the most potent A_2A AR antagonists displaying exceptional selectivity. While crystal structures of the human A_2A AR have been solved, mostly using the A_2A -StaR2 protein that bears 9 point mutations, co-crystallization with Preladenant derivatives has so far been elusive. We developed a new A_2A AR construct harboring a single point mutation (S91^3.39 K) which renders it extremely thermostable. This allowed the co-crystallization of two novel Preladenant derivatives, the polyethylene glycol-conjugated (PEGylated) PSB-2113, and the fluorophore-labeled PSB-2115. The obtained crystal structures (2.25 Å and 2.6 Å resolution) provide explanations for the high potency and selectivity of Preladenant derivatives. They represent the first crystal structures of a GPCR in complex with PEG- and fluorophore-conjugated ligands. The applied strategy is predicted to be applicable to further class A GPCRs.

Revisiting the Allosteric Regulation of Sodium Cation on the Binding of Adenosine at the Human A2A Adenosine Receptor: Insights from Supervised Molecular Dynamics (SuMD) Simulations

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A2A receptor (A2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Role of Extracellular Loops and Membrane Lipids for Ligand Recognition in the Neuronal Adenosine Receptor Type 2A: An Enhanced Sampling Simulation Study

Human G-protein coupled receptors (GPCRs) are important targets for pharmaceutical intervention against neurological diseases. Here, we use molecular simulation to investigate the key step in ligand recognition governed by the extracellular domains in the neuronal adenosine receptor type 2A (hA_2AR), a target for neuroprotective compounds. The ligand is the high-affinity antagonist (4-(2-(7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-ylamino)ethyl)phenol), embedded in a neuronal membrane mimic environment. Free energy calculations, based on well-tempered metadynamics, reproduce the experimentally measured binding affinity. The results are consistent with the available mutagenesis studies. The calculations identify a vestibular binding site, where lipids molecules can actively participate to stabilize ligand binding. Bioinformatic analyses suggest that such vestibular binding site and, in particular, the second extracellular loop, might drive the ligand toward the orthosteric binding pocket, possibly by allosteric modulation. Taken together, these findings point to a fundamental role of the interaction between extracellular loops and membrane lipids for ligands’ molecular recognition and ligand design in hA_2AR.

Exploring the mechanism of F282L mutation-caused constitutive activity of GPCR by a computational study

G-protein-coupled receptors (GPCRs) are important drug targets and generally activated by ligands. However, some experiments found that GPCRs also give rise to constitutive activity through some mutations (viz., CAM), which are usually associated with different kinds of diseases. However, the mechanisms of CAMs and their roles in interactions with drug-ligands are unclear in experiments. Herein, we used microsecond molecular dynamics simulations to study the effect of one important F282L mutation on β2AR in order to address the questions above. With the aid of principle component and correlation analysis, our results revealed that the F282L mutation could increase the instability of the overall structure, increase the dramatic fluctuations of NPxxY and extracellular loops, and decrease restraint of the helices through weakening interhelical H-bonding and correlations between residues, which could partly contribute to the constitutive activity reported by the experiments. The observations from the protein structure network (PSN) analysis indicate that the mutant exhibits less information flow than the wild β2AR and weakens the role of TM5 and TM6 in the signal transmission, but it enhances the impact of TM3 on the orthosteric pathway and TM4 on the allosteric one. In addition, the results from the virtual screening reveal that the mutant prefers to select agonists rather than antagonists, similar to the active state but opposite of the inactive state, further confirming that the F282L mutation advances the activation of β2AR. Our observations provide valuable information for understanding the mechanism of the mutation-caused constitutive activity of GPCR and related drug-design.

Investigation of the conformational dynamics of the apo A2A adenosine receptor

The activation/deactivation processes for G-protein coupled receptors (GPCRs) have been computationally studied for several different classes, including rhodopsin, the β2 adrenergic receptor, and the M2 muscarinic receptor. Despite determined cocrystal structures of the adenosine A2A receptor (A2A AR) in complex with antagonists, agonists and an antibody, the deactivation process of this GPCR is not completely understood. In this study, we investigate the convergence of two apo simulations, one starting with an agonist-bound conformation (PDB: 3QAK)(14) and the other starting with an antagonist-bound conformation (PDB: 3EML)(11) . Despite the two simulations not completely converging, we were able to identify distinct intermediate steps of the deactivation process characterized by the movement of Y288(7.53) in the NPxxY motif. We find that Y288(7.53) contributes to the process by forming hydrogen bonds to residues in transmembrane helices 2 and 7 and losing these interactions upon full deactivation. Y197(5.58) also plays a role in the process by forming a hydrogen bond only once the side chain moves from the lipid interface to the middle of the helical bundle.

Binding of the Antagonist Caffeine to the Human Adenosine Receptor hA2AR in Nearly Physiological Conditions

Lipid composition may significantly affect membrane proteins function, yet its impact on the protein structural determinants is not well understood. Here we present a comparative molecular dynamics (MD) study of the human adenosine receptor type 2A (hA(2A)R) in complex with caffeine--a system of high neuro-pharmacological relevance--within different membrane types. These are POPC, mixed POPC/POPE and cholesterol-rich membranes. 0.8-μs MD simulations unambiguously show that the helical folding of the amphipathic helix 8 depends on membrane contents. Most importantly, the distinct cholesterol binding into the cleft between helix 1 and 2 stabilizes a specific caffeine-binding pose against others visited during the simulation. Hence, cholesterol presence (~33%-50% in synaptic membrane in central nervous system), often neglected in X-ray determination of membrane proteins, affects the population of the ligand binding poses. We conclude that including a correct description of neuronal membranes may be very important for computer-aided design of ligands targeting hA(2A)R and possibly other GPCRs.

Communication over the network of binary switches regulates the activation of A2A adenosine receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.

As the key signal transmitters of a number of physiological processes, G-protein coupled receptors (GPCRs) are arguably one of the most important therapeutic targets. Orchestration of the intra-molecular signaling across transmembrane domain is key for the function of GPCRs. To investigate the microscopic underpinnings of intramolecular signaling that regulates the activation of GPCRs, we performed molecular dynamics simulations of the receptor in three distinct ligand-bound states using A_2A adenosine receptor as a model system of GPCRs. Statistical analyses on the dynamics of and correlation among the 10 “binary switches” reveal that the three receptor states retain distinct dynamic properties. The antagonist- and agonist-bound forms of the receptors explore vastly different conformational space, and the apo form lies between them, yet located closer to the antagonist-bound form. In regard to the agonist-binding triggered activation mechanism, the correlation map among the 10 binary switches unequivocally shows that direct sensing of agonist ligand by the indole ring of W246 actuates the rest of intramolecular signaling.

On the role of A₂A and D₂ receptors in control of cocaine and food-seeking behaviors in rats

Recent studies indicate that adenosine may influence dopamine neurotransmission via A2A receptors which antagonistically interact with D2 receptor-mediated signaling in the brain. We examined the effects of selective A2A receptor ligands such as the agonist CGS 21680 and the antagonists KW 6002 or SCH 58261 as well as of the D2-like receptor antagonist raclopride on reinstatement of cocaine seeking induced by cocaine, the cocaine-conditioned cue, or the D2-like receptor agonist quinpirole in rats. For comparison, effects of the A2A receptor ligands on reinstatement of food seeking were also studied. CGS 21680 significantly attenuated the reinstatement of cocaine (ip) seeking, and even more potently it reduced quinpirole (ip) or the cue-induced relapse of cocaine seeking as well as cue-induced food seeking. A potent reduction toward the cocaine-, quinpirole-, or cue-induced reinstatement of cocaine seeking was seen with raclopride. Pretreatment with KW 6002 or SCH 58261 reinstated cocaine seeking, and such increases were blocked by raclopride. In the higher doses, KW 6002 or SCH 58261 evoked food-seeking. In combination with the subthreshold dose of cocaine (2.5 mg/kg) or with the cue, low doses of KW 6002 but not SCH 58261 reinstated cocaine-seeking behavior, while none of the A2A receptor antagonists affected the cue-induced food-seeking behavior. The results indicate that A2A activation and D2-like receptor blockade counteract cocaine and food relapse. It is proposed that A2A receptor- and D2 receptor-mediated adenosine and dopamine signaling antagonistically interact in the striato-pallidal GABA neurons to regulate cocaine and food-seeking behavior.

An insight into antagonist binding and induced conformational dynamics of class B GPCR corticotropin-releasing factor receptor 1

The binding of small-molecule antagonists, CP-376395 and MTIP, would induce conformational dynamics behaviors of CRF₁R.