The GPR139 reference agonists 1a and 7c, and tryptophan and phenylalanine share a common binding site

Research progress on the role of orphan receptor GPR139 in neuropsychiatric behaviours

The study of orphan G protein-coupled receptors (GPCRs) holds much promise for increasing our understanding of neuropsychiatric diseases and for the development of new therapeutic strategies for these diseases. GPR139 is an orphan GPCR expressed in the central nervous system, especially in areas of the brain that control movement, motivation, and reward, and those that regulate neuropsychiatric behaviour. This review provides information about the discovery, tissue expression, signal transduction pathways, and physiological functions of GPR139, as well as how GPR139 interacts with other GPCRs, which form heteromeric complexes that affect their pharmacology and function. We also discuss the utility and therapeutic potential of ligands that target GPR139, including the pharmacological properties of reported agonists and antagonists. Finally, we highlight the pathologic role of GPR139 in neuropsychiatric behaviour and its potential as a therapeutic target in neuropsychiatric disorders.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Pharmacological characterization of novel small molecule agonists and antagonists for the orphan receptor GPR139

The orphan G protein-coupled receptor GPR139 is predominantly expressed in the central nervous system and has attracted considerable interest as a therapeutic target. However, the biological role of this receptor remains somewhat elusive, in part due to the lack of quality pharmacological tools to investigate GPR139 function. In an effort to understand GPR139 signaling and to identify improved compounds, in this study we performed virtual screening and analog searches, in combination with multiple pharmacological assays. We characterized GPR139-dependent signaling using previously published reference agonists in Ca²⁺ mobilization and inositol monophosphate accumulation assays, as well as a novel real-time GPR139 internalization assay. For the four reference agonists tested, the rank order of potency was conserved across signaling and internalization assays: JNJ-63533054 > Compound 1a » Takeda > AC4 > DL43, consistent with previously reported values. We noted an increased efficacy of JNJ-63533054-mediated inositol monophosphate signaling and internalization, relative to Compound 1a. We then performed virtual screening for GPR139 agonist and antagonist compounds that were screened and validated in GPR139 functional assays. We identified four GPR139 agonists that were active in all assays, with similar or reduced potency relative to known compounds. Likewise, compound analogs selected based on GPR139 agonist and antagonist substructure searches behaved similarly to their parent compounds. Thus, we have characterized GPR139 signaling for multiple new ligands using G protein-dependent assays and a new real-time internalization assay. These data add to the GPR139 tool compound repertoire, which could be optimized in future medical chemistry campaigns.

The orphan G protein-coupled receptor, GPR139, is expressed in the hypothalamus and is involved in the regulation of body mass, blood glucose, and insulin

GPR139 is an orphan G-protein-coupled receptor that is expressed in restricted areas of the nervous system, including the hypothalamus. In this study, we hypothesized that GPR139 could be involved in the regulation of energy balance and metabolism. In the first part of the study, we confirmed that GPR139 is expressed in the hypothalamus and particularly in proopiomelanocortin and agouti-related peptide neurons of the mediobasal hypothalamus. Using a lentivirus with a short-hairpin RNA, we inhibited the expression of GPR139 bilaterally in the mediobasal hypothalamus of mice. The intervention promoted a 40% reduction in the hypothalamic expression of GPR139, which was accompanied by an increase in body mass, a reduction in fasting blood glucose levels, and an increase in insulin levels. In the hypothalamus, inhibition of GPR139 was accompanied by a reduction in the expression of orexin. As previous studies using a pharmacological antagonist of orexin showed a beneficial impact on type 2 diabetes and glucose metabolism, we propose that the inhibition of hypothalamic GPR139 could be acting indirectly through the orexin system to control systemic glucose and insulin. In conclusion, this study advances the characterization of GPR139 in the hypothalamus, demonstrating its involvement in the regulation of body mass, blood insulin, and glycemia.

Aromatic patterns: Tryptophan aromaticity as a catalyst for the emergence of life and rise of consciousness

Sunlight held the key to the origin of life on Earth. The earliest life forms, cyanobacteria, captured the sunlight to generate energy through photosynthesis. Life on Earth evolved in accordance with the circadian rhythms tied to sensitivity to sunlight patterns. A unique feature of cyanobacterial photosynthetic proteins and circadian rhythms' molecules, and later of nearly all photon-sensing molecules throughout evolution, is that the aromatic amino acid tryptophan (Trp) resides at the center of light-harvesting active sites. In this perspective, I review the literature and integrate evidence from different scientific fields to explore the role Trp plays in photon-sensing capabilities of living organisms through its resonance delocalization of π-electrons. The observations presented here are the product of apparently unrelated phenomena throughout evolution, but nevertheless share consistent patterns of photon-sensing by Trp-containing and Trp-derived molecules. I posit the unique capacity to transfer electrons during photosynthesis in the earliest life forms is conferred to Trp due to its aromaticity. I propose this ability evolved to assume more complex functions, serving as a host for mechanisms underlying mental aptitudes - a concept which provides a theoretical basis for defining the neural correlates of consciousness. The argument made here is that Trp aromaticity may have allowed for the inception of the mechanistic building blocks used to fabricate complexity in higher forms of life. More specifically, Trp aromatic non-locality may have acted as a catalyst for the emergence of consciousness by instigating long-range synchronization and stabilizing the large-scale coherence of neural networks, which mediate functional brain activity. The concepts proposed in this perspective provide a conceptual foundation that invites further interdisciplinary dialogue aimed at examining and defining the role of aromaticity (beyond Trp) in the emergence of life and consciousness.

The role of orphan receptor GPR139 in neuropsychiatric behavior

Orphan G protein Coupled Receptors (GPCRs) present attractive targets both for understanding neuropsychiatric diseases and for development of novel therapeutics. GPR139 is an orphan GPCR expressed in select brain circuits involved in controlling movement, motivation and reward. It has been linked to the opioid and dopamine neuromodulatory systems; however, its role in animal behavior and neuropsychiatric processes is poorly understood. Here we present a comprehensive behavioral characterization of a mouse model with a GPR139 null mutation. We show that loss of GPR139 in mice results in delayed onset hyperactivity and prominent neuropsychiatric manifestations including elevated stereotypy, increased anxiety-related traits, delayed acquisition of operant responsiveness, disruption of cued fear conditioning and social interaction deficits. Furthermore, mice lacking GPR139 exhibited complete loss of pre-pulse inhibition and developed spontaneous 'hallucinogenic' head-twitches, altogether suggesting schizophrenia-like symptomatology. Remarkably, a number of these behavioral deficits could be rescued by the administration of μ-opioid and D2 dopamine receptor (D2R) antagonists: naltrexone and haloperidol, respectively, suggesting that loss of neuropsychiatric manifestations in mice lacking GPR139 are driven by opioidergic and dopaminergic hyper-functionality. The inhibitory influence of GPR139 on D2R signaling was confirmed in cell-based functional assays. These observations define the role of GPR139 in controlling behavior and implicate in vivo actions of this receptor in the neuropsychiatric process with schizophrenia-like pathology.

Metabolite G-Protein Coupled Receptors in Cardio-Metabolic Diseases

G protein-coupled receptors (GPCRs) have originally been described as a family of receptors activated by hormones, neurotransmitters, and other mediators. However, in recent years GPCRs have shown to bind endogenous metabolites, which serve functions other than as signaling mediators. These receptors respond to fatty acids, mono- and disaccharides, amino acids, or various intermediates and products of metabolism, including ketone bodies, lactate, succinate, or bile acids. Given that many of these metabolic processes are dysregulated under pathological conditions, including diabetes, dyslipidemia, and obesity, receptors of endogenous metabolites have also been recognized as potential drug targets to prevent and/or treat metabolic and cardiovascular diseases. This review describes G protein-coupled receptors activated by endogenous metabolites and summarizes their physiological, pathophysiological, and potential pharmacological roles.

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.

Discovery of TAK-041: a Potent and Selective GPR139 Agonist Explored for the Treatment of Negative Symptoms Associated with Schizophrenia

The orphan G-protein-coupled receptor GPR139 is highly expressed in the habenula, a small brain nucleus that has been linked to depression, schizophrenia (SCZ), and substance-use disorder. High-throughput screening and a medicinal chemistry structure-activity relationship strategy identified a novel series of potent and selective benzotriazinone-based GPR139 agonists. Herein, we describe the chemistry optimization that led to the discovery and validation of multiple potent and selective in vivo GPR139 agonist tool compounds, including our clinical candidate TAK-041, also known as NBI-1065846 (compound 56). The pharmacological characterization of these GPR139 agonists in vivo demonstrated GPR139-agonist-dependent modulation of habenula cell activity and revealed consistent in vivo efficacy to rescue social interaction deficits in the BALB/c mouse strain. The clinical GPR139 agonist TAK-041 is being explored as a novel drug to treat negative symptoms in SCZ.

Free Energy Calculations for Protein-Ligand Binding Prediction

Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.

X-Ray Crystallography and Free Energy Calculations Reveal the Binding Mechanism of A2A Adenosine Receptor Antagonists

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A2A adenosine receptor (AR). Eight A2A AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2A AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A2A AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2A AR, an emerging target in immuno-oncology.

Pharmacology and function of the orphan GPR139 G protein-coupled receptor

G protein-coupled receptors (GPCRs) constitute the largest family of receptors and membrane proteins in the human genome with ~800 members of which half are olfactory. GPCRs are activated by a very broad range of endogenous signalling molecules and are involved in a plethora of physiological functions. All GPCRs contain a transmembrane domain, consisting of a bundle of seven α-helices spanning the cell membrane, and forming the majority of the known ortho- or allosteric ligand binding sites. Due to their many physiological functions and the accessible and druggable transmembrane pocket, GPCRs constitute the largest family of drug targets mediating the actions of 34% of currently marketed drugs. GPCRs activate one or more of the four G protein families (G_q/11 , G_i/o , G_s and G_12/13 ) and/or ß-arrestin. About a third of the non-olfactory GPCRs are referred to as orphan receptors which means that their endogenous agonist(s) have not yet been found or firmly established. In this MiniReview, we focus on the orphan GPR139 receptor, for which the aromatic amino acids L-Trp and L-Phe as well as ACTH/α-MSH-related peptides have been proposed as endogenous agonists. GPR139 has been reported to activate several G protein pathways of which G_q/11 is the primary one. The receptor shows the highest expression in the striatum, thalamus, hypothalamus, pituitary and habenula of the human, rat and mouse CNS. We review the surrogate agonists and antagonists that have been published as well as the agonist pharmacophore and binding site. Finally, the putative physiological functions and therapeutic potential are outlined.

In vivo Characterization of a Selective, Orally Available, and Brain Penetrant Small Molecule GPR139 Agonist

Recently, our group along with another demonstrated that GPR139 can be activated by L-phenylalanine (L-Phe) and L-tryptophan (L-Trp) at physiologically relevant concentrations. GPR139 is discretely expressed in brain, with highest expression in medial habenula. Not only are the endogenous ligands catecholamine/serotonin precursors, but GPR139 expressing areas can directly/indirectly regulate the activity of catecholamine/serotonin neurons. Thus, GPR139 appears expressed in an interconnected circuit involved in mood, motivation, and anxiety. The aim of this study was to characterize a selective and brain penetrant GPR139 agonist (JNJ-63533054) in relevant in vivo models. JNJ-63533054 was tested for its effect on c-fos activation in the habenula and dorsal striatum. In vivo microdialysis experiments were performed in freely moving rats to measure basal levels of serotonin or dopamine (DA) in prefrontal cortex (mPFC) and nucleus accumbens (NAc). Finally, the compound was profiled in behavioral models of anxiety, despair, and anhedonia. The agonist (10–30 mg/kg, p.o.) did not alter c-fos expression in medial habenula or dorsal striatum nor neurotransmitter levels in mPFC or NAc. JNJ-63533054 (10 mg/kg p.o.) produced an anhedonic-like effect on urine sniffing, but had no significant effect in tail suspension, with no interaction with imipramine, no effect on naloxone place aversion, and no effect on learned helplessness. In the marble burying test, the agonist (10 mg/kg p.o.) produced a small anxiolytic-like effect, with no interaction with fluoxetine, and no effect in elevated plus maze (EPM). Despite GPR139 high expression in medial habenula, an area with connections to limbic and catecholaminergic/serotoninergic areas, the GPR139 agonist had no effect on c-fos in medial habenula. It did not alter catecholamine/serotonin levels and had a mostly silent signal in in vivo models commonly associated with these pathways. The physiological function of GPR139 remains elusive.

Identification of a novel scaffold for a small molecule GPR139 receptor agonist

GPR139 is an orphan G protein-coupled receptor (GPCR) that is primarily expressed in the brain in regions known to regulate motor control and metabolism. Here, we screened a diverse 4,000 compound library in order to identify GPR139 agonists. We identified 11 initial hits in a calcium mobilization screen, including one compound, AC4, which contains a different chemical scaffold to what has previously been described for GPR139 agonists. Our mutagenesis data shows that AC4 interacts with the same hotspots in the binding site of GPR139 as those reported to interact with the reference agonists 1a and 7c. We additionally tested and validated 160 analogs in a calcium mobilization assay and found 5 compounds with improved potency compared to AC4. In total, we identified 36 GPR139 agonists with potencies in the nanomolar range (90–990 nM). The most potent compounds were confirmed as GPR139 agonists using an orthogonal ERK phosphorylation assay where they displayed a similar rank order of potency. Accordingly, we herein introduce multiple novel GPR139 agonists, including one with a novel chemical scaffold, which can be used as tools for future pharmacological and medicinal chemistry exploration of GPR139.

Mutagenesis of GPR139 reveals ways to create gain or loss of function receptors

GPR139 is a Gq-coupled receptor activated by the essential amino acids L-tryptophan (L-Trp) and L-phenylalanine (L-Phe). We carried out mutagenesis studies of the human GPR139 receptor to identify the critical structural motifs required for GPR139 activation. We applied site-directed and high throughput random mutagenesis approaches using a double addition normalization strategy to identify novel GPR139 sequences coding receptors that have altered sensitivity to endogenous ligands. This approach resulted in GPR139 clones with gain-of-function, reduction-of-function or loss-of-function mutations. The agonist pharmacology of these mutant receptors was characterized and compared to wild-type receptor using calcium mobilization, radioligand binding, and protein expression assays. The structure-activity data were incorporated into a homology model which highlights that many of the gain-of-function mutations are either in or immediately adjacent to the purported orthosteric ligand binding site, whereas the loss-of-function mutations were largely in the intracellular G-protein binding area or were disrupters of the helix integrity. There were also some reduction-of-function mutations in the orthosteric ligand binding site. These findings may not only facilitate the rational design of novel agonists and antagonists of GPR139, but also may guide the design of transgenic animal models to study the physiological function of GPR139.

Re-evaluation of Adrenocorticotropic Hormone and Melanocyte Stimulating Hormone Activation of GPR139 in Vitro

It is now well established that GPR139, a G-protein coupled receptor exclusively expressed in the brain and pituitary, is activated by the essential amino acids L-tryptophan (L-Trp) and L-phenylalanine (L-Phe) via G_αq-coupling. The in vitro affinity and potency values of L-Trp and L-Phe are within the physiological concentration ranges of L-Trp and L-Phe. A recent paper suggests that adrenocorticotropic hormone (ACTH), α and β melanocyte stimulating hormones (α-MSH and β-MSH) and derivatives α-MSH_1-9/α-MSH_1-10 can also activate GPR139 in vitro. We tested this hypothesis using guanosine 5′-O-(3-[³⁵S]thio)-triphosphate binding (GTPγS), calcium mobilization and [³H]JNJ-63533054 radioligand binding assays. In the GTPγS binding assay, α-MSH, α-MSH_1-9/α-MSH_1-10, and β-MSH had no effect on [³⁵S]GTPγS incorporation in cell membranes expressing GPR139 up to 30 μM in contrast to the concentration dependent activation produced by L-Trp, JNJ-63533054, and TC-09311 (two small molecule GPR139 agonists). ACTH slightly decreased the basal level of [³⁵S]GTPγS incorporation at 30 μM. In the GPR139 radioligand binding assay, a moderate displacement of [³H]JNJ-63533054 binding by ACTH and β-MSH was observed at 30 μM (40 and 30%, respectively); α-MSH, α-MSH_1-9/α-MSH_1-10 did not displace any specific binding at 30 μM. In three different host cell lines stably expressing GPR139, α-MSH, and β-MSH did not stimulate calcium mobilization in contrast to L-Trp, JNJ-63533054, and TC-09311. ACTH, α-MSH_1-9/α-MSH_1-10 only weakly stimulated calcium mobilization at 30 μM (<50% of EC₁₀₀). We then co-transfected GPR139 with the three melanocortin (MC) receptors (MC3R, MC4R, and MC5R) to test the hypothesis that ACTH, α-MSH, and β-MSH might stimulate calcium mobilization through a MCR/GPR139 interaction. All three MC peptides stimulated calcium response in cells co-transfected with GPR139 and MC3R, MC4R, or MC5R. The MC peptides did not stimulate calcium response in cells expressing MC3R or MC5R alone consistent with the G_s signaling transduction pathway of these receptors. In agreement with the previously reported multiple signaling pathways of MC4R, including G_q transduction pathway, the MC peptides produced a calcium response in cells expressing MC4R alone. Together, our findings do not support that GPR139 is activated by ACTH, α-MSH, and β-MSH at physiologically relevant concentration but we did unravel an in vitro interaction between GPR139 and the MCRs.

The G protein-coupled receptors deorphanization landscape

G protein-coupled receptors (GPCRs) are usually highlighted as being both the largest family of membrane proteins and the most productive source of drug targets. However, most of the GPCRs are understudied and hence cannot be used immediately for innovative therapeutic strategies. Besides, there are still around 100 orphan receptors, with no described endogenous ligand and no clearly defined function. The race to discover new ligands for these elusive receptors seems to be less intense than before. Here, we present an update of the various strategies employed to assign a function to these receptors and to discover new ligands. We focus on the recent advances in the identification of endogenous ligands with a detailed description of newly deorphanized receptors. Replication being a key parameter in these endeavors, we also discuss the latest controversies about problematic ligand-receptor pairings. In this context, we propose several recommendations in order to strengthen the reporting of new ligand-receptor pairs.

Elucidation of the Binding Mode of the Carboxyterminal Region of Peptide YY to the Human Y2 Receptor

Understanding the agonist-receptor interactions in the neuropeptide Y (NPY)/peptide YY (PYY) signaling system is fundamental for the design of novel modulators of appetite regulation. We report here the results of a multidisciplinary approach to elucidate the binding mode of the native peptide agonist PYY to the human Y₂ receptor, based on computational modeling, peptide chemistry and in vitro pharmacological analyses. The preserved binding orientation proposed for full-length PYY and five analogs, truncated at the amino terminus, explains our pharmacological results where truncations of the N-terminal proline helix showed little effect on peptide affinity. This was followed by receptor mutagenesis to investigate the roles of several receptor positions suggested by the modeling. As a complement, PYY-(3-36) analogs were synthesized with modifications at different positions in the common PYY/NPY C-terminal fragment (³²TRQRY³⁶-amide). The results were assessed and interpreted by molecular dynamics and Free Energy Perturbation (FEP) simulations of selected mutants, providing a detailed map of the interactions of the PYY/NPY C-terminal fragment with the transmembrane cavity of the Y₂ receptor. The amidated C-terminus would be stabilized by polar interactions with Gln288^6.55 and Tyr219^5.39, while Gln130^3.32 contributes to interactions with Q³⁴ in the peptide and T³² is close to the tip of TM7 in the receptor. This leaves the core, α-helix of the peptide exposed to make potential interactions with the extracellular loops. This model agrees with most experimental data available for the Y₂ system and can be used as a basis for optimization of Y₂ receptor agonists.

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.