Structure-based discovery of conformationally selective inhibitors of the serotonin transporter.

The serotonin transporter (SERT) removes synaptic serotonin and is the target of anti-depressant drugs. SERT adopts three conformations: outward-open, occluded, and inward-open. All known inhibitors target the outward-open state except ibogaine, which has unusual anti-depressant and substance-withdrawal effects, and stabilizes the inward-open conformation. Unfortunately, ibogaine's promiscuity and cardiotoxicity limit the understanding of inward-open state ligands. We docked over 200 million small molecules against the inward-open state of the SERT. Thirty-six top-ranking compounds were synthesized, and thirteen inhibited; further structure-based optimization led to the selection of two potent (low nanomolar) inhibitors. These stabilized an outward-closed state of the SERT with little activity against common off-targets. A cryo-EM structure of one of these bound to the SERT confirmed the predicted geometry. In mouse behavioral assays, both compounds had anxiolytic- and anti-depressant-like activity, with potencies up to 200-fold better than fluoxetine (Prozac), and one substantially reversed morphine withdrawal effects.

How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling.

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Angiotensin Analogs with Divergent Bias Stabilize Distinct Receptor Conformations.

"Biased" G protein-coupled receptor (GPCR) agonists preferentially activate pathways mediated by G proteins or ß-arrestins. Here, we use double electron-electron resonance spectroscopy to probe the changes that ligands induce in the conformational distribution of the angiotensin II type I receptor. Monitoring distances between 10 pairs of nitroxide labels distributed across the intracellular regions enabled mapping of four underlying sets of conformations. Ligands from different functional classes have distinct, characteristic effects on the conformational heterogeneity of the receptor. Compared to angiotensin II, the endogenous agonist, agonists with enhanced Gq coupling more strongly stabilize an "open" conformation with an accessible transducer-binding site. ß-arrestin-biased agonists deficient in Gq coupling do not stabilize this open conformation but instead favor two more occluded conformations. These data suggest a structural mechanism for biased ligand action at the angiotensin receptor that can be exploited to rationally design GPCR-targeting drugs with greater specificity of action.

Functional Selectivity in Cytokine Signaling Revealed Through a Pathogenic EPO Mutation.

Cytokines are classically thought to stimulate downstream signaling pathways through monotonic activation of receptors. We describe a severe anemia resulting from a homozygous mutation (R150Q) in the cytokine erythropoietin (EPO). Surprisingly, the EPO R150Q mutant shows only a mild reduction in affinity for its receptor but has altered binding kinetics. The EPO mutant is less effective at stimulating erythroid cell proliferation and differentiation, even at maximally potent concentrations. While the EPO mutant can stimulate effectors such as STAT5 to a similar extent as the wild-type ligand, there is reduced JAK2-mediated phosphorylation of select downstream targets. This impairment in downstream signaling mechanistically arises from altered receptor dimerization dynamics due to extracellular binding changes. These results demonstrate how variation in a single cytokine can lead to biased downstream signaling and can thereby cause human disease. Moreover, we have defined a distinct treatable form of anemia through mutation identification and functional studies.

Bias Factor and Therapeutic Window Correlate to Predict Safer Opioid Analgesics.

Biased agonism has been proposed as a means to separate desirable and adverse drug responses downstream of G protein-coupled receptor (GPCR) targets. Herein, we describe structural features of a series of mu-opioid-receptor (MOR)-selective agonists that preferentially activate receptors to couple to G proteins or to recruit ßarrestin proteins. By comparing relative bias for MOR-mediated signaling in each pathway, we demonstrate a strong correlation between the respiratory suppression/antinociception therapeutic window in a series of compounds spanning a wide range of signaling bias. We find that ßarrestin-biased compounds, such as fentanyl, are more likely to induce respiratory suppression at weak analgesic doses, while G protein signaling bias broadens the therapeutic window, allowing for antinociception in the absence of respiratory suppression.

GPCR-IPL score: multilevel featurization of GPCR-ligand interaction patterns and prediction of ligand functions from selectivity to biased activation.

G-protein-coupled receptors (GPCRs) mediate diverse cell signaling cascades after recognizing extracellular ligands. Despite the successful history of known GPCR drugs, a lack of mechanistic insight into GPCR challenges both the deorphanization of some GPCRs and optimization of the structure-activity relationship of their ligands. Notably, replacing a small substituent on a GPCR ligand can significantly alter extracellular GPCR-ligand interaction patterns and motion of transmembrane helices in turn to occur post-binding events of the ligand. In this study, we designed 3D multilevel features to describe the extracellular interaction patterns. Subsequently, these 3D features were utilized to predict the post-binding events that result from conformational dynamics from the extracellular to intracellular areas. To understand the adaptability of GPCR ligands, we collected the conformational information of flexible residues during binding and performed molecular featurization on a broad range of GPCR-ligand complexes. As a result, we developed GPCR-ligand interaction patterns, binding pockets, and ligand features as score (GPCR-IPL score) for predicting the functional selectivity of GPCR ligands (agonism versus antagonism), using the multilevel features of (1) zoomed-out 'residue level' (for flexible transmembrane helices of GPCRs), (2) zoomed-in 'pocket level' (for sophisticated mode of action) and (3) 'atom level' (for the conformational adaptability of GPCR ligands). GPCR-IPL score demonstrated reliable performance, achieving area under the receiver operating characteristic of 0.938 and area under the precision-recall curve of 0.907 (available in gpcr-ipl-score.onrender.com). Furthermore, we used the molecular features to predict the biased activation of downstream signaling (Gi/o, Gq/11, Gs and ß-arrestin) as well as the functional selectivity. The resulting models are interpreted and applied to out-of-set validation with three scenarios including the identification of a new MRGPRX antagonist.

Discovery of a functionally selective ghrelin receptor (GHSR1a) ligand for modulating brain dopamine.

SignificanceThe modulation of growth hormone secretagogue receptor-1a (GHSR1a) signaling is a promising strategy for treating brain conditions of metabolism, aging, and addiction. GHSR1a activation results in pleiotropic physiological outcomes through distinct and pharmacologically separable G protein- and ß-arrestin (ßarr)-dependent signaling pathways. Thus, pathway-selective modulation can enable improved pharmacotherapeutics that can promote therapeutic efficacy while mitigating side effects. Here, we describe the discovery of a brain-penetrant small molecule, N8279 (NCATS-SM8864), that biases GHSR1a conformations toward Gαq activation and reduces aberrant dopaminergic behavior in mice. N8279 represents a promising chemical scaffold to advance the development of better treatments for GHSR1a-related brain disorders involving the pathological dysregulation of dopamine.

A structural mechanism of nuclear receptor biased agonism.

Efforts to decrease the adverse effects of nuclear receptor (NR) drugs have yielded experimental agonists that produce better outcomes in mice. Some of these agonists have been shown to cause different, not just less intense, on-target transcriptomic effects; however, a structural explanation for such agonist-specific effects remains unknown. Here, we show that partial agonists of the NR peroxisome proliferator-associated receptor γ (PPARγ), which induce better outcomes in mice compared to clinically utilized type II diabetes PPARγ-binding drugs thiazolidinediones (TZDs), also favor a different group of coactivator peptides than the TZDs. We find that PPARγ full agonists can also be biased relative to each other in terms of coactivator peptide binding. We find differences in coactivator-PPARγ bonding between the coactivator subgroups which allow agonists to favor one group of coactivator peptides over another, including differential bonding to a C-terminal residue of helix 4. Analysis of all available NR-coactivator structures indicates that such differential helix 4 bonding persists across other NR-coactivator complexes, providing a general structural mechanism of biased agonism for many NRs. Further work will be necessary to determine if such bias translates into altered coactivator occupancy and physiology in cells.

Efficient functional localization of language regions in the brain.

Important recent advances in the cognitive neuroscience of language have been made using functional localizers to demarcate language-selective regions in individual brains. Although single-subject localizers offer insights that are unavailable in classic group analyses, they require additional scan time that imposes costs on investigators and participants. In particular, the unique practical challenges of scanning children and other special populations has led to less adoption of localizers for neuroimaging research with these theoretically and clinically important groups. Here, we examined how measurements of the spatial extent and functional response profiles of language regions are affected by the duration of an auditory language localizer. We compared how parametrically smaller amounts of data collected from one scanning session affected (i) consistency of group-level whole-brain parcellations, (ii) functional selectivity of subject-level activation in individually defined functional regions of interest (fROIs), (iii) sensitivity and specificity of subject-level whole-brain and fROI activation, and (iv) test-retest reliability of subject-level whole-brain and fROI activation. For many of these metrics, the localizer duration could be reduced by 50-75% while preserving the stability and reliability of both the spatial extent and functional response profiles of language areas. These results indicate that, for most measures relevant to cognitive neuroimaging studies, the brain's language network can be localized just as effectively with 3.5 min of scan time as it can with 12 min. Minimizing the time required to reliably localize the brain's language network allows more effective localizer use in situations where each minute of scan time is particularly precious.

Systematic assessment of chemokine ligand bias at the human chemokine receptor CXCR2 indicates G protein bias over ß-arrestin recruitment and receptor internalization.

BACKGROUND: The human CXC chemokine receptor 2 (CXCR2) is a G protein-coupled receptor (GPCR) interacting with multiple chemokines (i.e., CXC chemokine ligands CXCL1-3 and CXCL5-8). It is involved in inflammatory diseases as well as cancer. Consequently, much effort is put into the identification of CXCR2 targeting drugs. Fundamental research regarding CXCR2 signaling is mainly focused on CXCL8 (IL-8), which is the first and best described high-affinity ligand for CXCR2. Much less is known about CXCR2 activation induced by other chemokines and it remains to be determined to what extent potential ligand bias exists within this signaling system. This insight might be important to unlock new opportunities in therapeutic targeting of CXCR2. METHODS: Ligand binding was determined in a competition binding assay using labeled CXCL8. Activation of the ELR + chemokine-induced CXCR2 signaling pathways, including G protein activation, ß-arrestin1/2 recruitment, and receptor internalization, were quantified using NanoBRET-based techniques. Ligand bias within and between these pathways was subsequently investigated by ligand bias calculations, with CXCL8 as the reference CXCR2 ligand. Statistical significance was tested through a one-way ANOVA followed by Dunnett's multiple comparisons test. RESULTS: All chemokines (CXCL1-3 and CXCL5-8) were able to displace CXCL8 from CXCR2 with high affinity and activated the same panel of G protein subtypes (Gαi1, Gαi2, Gαi3, GαoA, GαoB, and Gα15) without any statistically significant ligand bias towards any one type of G protein. Compared to CXCL8, all other chemokines were less potent in ß-arrestin1 and -2 recruitment and receptor internalization while equivalently activating G proteins, indicating a G protein activation bias for CXCL1,-2,-3,-5,-6 and CXCL7. Lastly, with CXCL8 used as reference ligand, CXCL2 and CXCL6 showed ligand bias towards ß-arrestin1/2 recruitment compared to receptor internalization. CONCLUSION: This study presents an in-depth analysis of signaling bias upon CXCR2 stimulation by its chemokine ligands. Using CXCL8 as a reference ligand for bias index calculations, no ligand bias was observed between chemokines with respect to activation of separate G proteins subtypes or recruitment of ß-arrestin1/2 subtypes, respectively. However, compared to ß-arrestin recruitment and receptor internalization, CXCL1-3 and CXCL5-7 were biased towards G protein activation when CXCL8 was used as reference ligand.

Phosphoproteomic Analysis of Dopamine D2 Receptor Signaling Reveals Interplay of G Protein- and ß-Arrestin-Mediated Effects.

Leveraging biased signaling of G protein-coupled receptors has been proposed as a promising strategy for the development of drugs with higher specificity. However, the consequences of selectively targeting G protein- or ß-arrestin-mediated signaling on cellular functions are not comprehensively understood. In this study, we utilized phosphoproteomics to gain a systematic overview of signaling induced by the four biased and balanced dopamine D2 receptor (D2R) ligands MS308, BM138, quinpirole, and sulpiride in an in vitro D2R transfection model. Quantification of 14,160 phosphosites revealed a low impact of the partial G protein agonist MS308 on cellular protein phosphorylation, as well as surprising similarities between the balanced agonist quinpirole and the inverse agonist sulpiride. Analysis of the temporal profiles of ligand-induced phosphorylation events showed a transient impact of the G protein-selective agonist MS308, whereas the ß-arrestin-preferring agonist BM138 elicited a delayed, but more pronounced response. Functional enrichment analysis of ligand-impacted phosphoproteins and treatment-linked kinases confirmed multiple known functions of D2R signaling while also revealing novel effects, for example of MS308 on sterol regulatory element-binding protein-related gene expression. All raw data were deposited in MassIVE (MSV000089457).

Computational insights into ligand-induced G protein and ß-arrestin signaling of the dopamine D1 receptor.

The dopamine D1 receptor (D1R), is a class A G protein coupled-receptor (GPCR) which has been a promising drug target for psychiatric and neurological disorders such as Parkinson's disease (PD). Previous studies have suggested that therapeutic effects can be realized by targeting the ß-arrestin signaling pathway of dopamine receptors, while overactivation of the G protein-dependent pathways leads to side effects, such as dyskinesias. Therefore, it is highly desirable to develop a D1R ligand that selectively regulates the ß-arrestin pathway. Currently, most D1R agonists are signaling-balanced and stimulate both G protein and ß-arrestin pathways, with a few reports of G protein biased ligands. However, identification and characterization of ß-arrestin biased D1R agonists has been a challenge thus far. In this study, we implemented Gaussian accelerated molecular dynamics (GaMD) simulations to provide valuable computational insights into the possible underlying molecular mechanism of the different signaling properties of two catechol and two non-catechol D1R agonists that are either G protein biased or signaling-balanced. Dynamic network analysis further identified critical residues in the allosteric signaling network of D1R for each ligand at different conformational or binding states. Some of these residues are crucial for G protein or arrestin signals of GPCRs based on previous studies. Finally, we provided a molecular design strategy which can be utilized by medicinal chemists to develop potential ß-arrestin biased D1R ligands. The proposed hypotheses are experimentally testable and can guide the development of safer and more effective medications for a variety of CNS disorders.

Modulating the affinity and signaling bias of cannabinoid receptor 1 antagonists.

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor and a therapeutic target for metabolic disorders. Numerous CB1 antagonists have been developed, but their functional selectivities and bias towards G protein or ß-arrestin signaling have not been systemically characterized. In this study, we analyzed the binding affinities and downstream signaling of two series of pyrazole derivatives bearing 1-aminopiperidine (Series I) or 4-aminothiomorpholine 1,1-dioxide (Series II) moieties, as well as the well-known CB1 antagonists rimonabant and taranabant. Analyses of the results for the Series I and II derivatives showed that minor structure modifications to their functional groups and especially the incorporation of 1-aminopiperidine or 4-aminothiomorpholine 1,1-dioxide motifs can profoundly affect their bias toward G protein or ß-arrestin signaling, and that their binding affinity and functional activity can be disassociated. Docking and molecular dynamics simulations revealed that the binding modes of Series I and II antagonists differed primarily in that Series I antagonists formed an additional hydrogen bond with the receptor, whereas those in Series II formed a water bridge.

Dynamics of the Second Extracellular Loop Control Transducer Coupling of Peptide-Activated GPCRs.

Many peptide-activated rhodopsin-like GPCRs share a ß-hairpin folding motif in the extracellular loop 2 (ECL2), which interacts with the peptide ligand while at the same time being connected to transmembrane helix 3 (TM3) via a highly conserved disulfide bond. Currently, it remains unknown whether the coupling of the specifically shaped ECL2 to TM3 influences the activation of peptide-activated GPCRs. We investigated this possibility in a selection of peptide GPCRs with known structures. Most of the receptors with cysteine to alanine mutations folded like the respective wild-type and resided in the cell membrane, challenging pure folding stabilization by the disulfide bridge. G-protein signaling of the disulfide mutants was retained to a greater extent in secretin-like GPCRs than in rhodopsin-like GPCRs, while recruitment of arrestin was completely abolished in both groups, which may be linked to alterations in ligand residence time. We found a correlation between receptor activity of the neuropeptide Y2 receptor and alterations in ECL2 dynamics using engineered disulfide bridges or site-directed spin labeling and EPR spectroscopy. These data highlight the functional importance of the TM3-ECL2 link for the activation of specific signaling pathways in peptide-activated GPCRs, which might have implications for future drug discovery.

The Mechanisms of GPR55 Receptor Functional Selectivity during Apoptosis and Proliferation Regulation in Cancer Cells.

GPR55 is a non-canonical cannabinoid receptor, important for cancer proliferation. Depending on the ligand, it induces either cell proliferation or death. The objective of the study was to establish the mechanisms of this multidirectional signaling. Using the CRISPR-Cas9 system, the GPR55, CB1, CB2, and GPR18 receptor knockouts of the MDA-MB-231 line were obtained. After the CB2 receptor knockout, the pro-apoptotic activity of the pro-apoptotic ligand docosahexaenoyl dopamine (DHA-DA) slightly increased, while the pro-proliferative activity of the most active synthetic ligand of the GPR55 receptor (ML-184) completely disappeared. On the original cell line, the stimulatory effect of ML-184 was removed by the CB2 receptor blocker and by GPR55 receptor knockout. Thus, it can be confidently assumed that when proliferation is stimulated with the participation of the GPR55 receptor, a signal is transmitted from the CB2 receptor to the GPR55 receptor due to the formation of a heterodimer. GPR18 was additionally involved in the implementation of the pro-apoptotic effect of DHA-DA, while the CB1 receptor is not involved. In the implementation of the pro-apoptotic action of DHA-DA, the elimination of Gα13 led to a decrease in cytotoxicity. The obtained data provide novel details to the mechanism of the pro-proliferative action of GPR55.

Structural optimizations and bioevaluation of N-H aporphine analogues as Gq-biased and selective serotonin 5-HT2C receptor agonists.

The concept of subtype selectivity and functional bias has recently reshaped current GPCR drug discovery for G protein-coupled receptors. A series of new N-H aporphines with A-ring modifications have been synthesized, and their efficacy on 5-HT2 receptor activation was evaluated. SAR analysis led to the discovery of several more potent and selective 5-HT2C receptor agonists (e.g., 11b and 11f) with low nanomolar activity. Molecular docking analysis of this series of aporphines was in accordance with our SAR results. The functional selectivity of specific compounds was tested via both Gq-mediated calcium flux and ß-arrestin recruitment assays, which revealed that these compounds exhibited no ß-arrestin recruitment activity. Further ADMET study combined with behavioral assessment using a methamphetamine-induced hyperactivity model identified compound 11b and 11f possessing promising drug-like profiles and having antipsychotic potential. These agonists with an exclusive bias toward Gq signaling may serve as valuable pharmacological probes to facilitate the elucidation of therapeutically relevant 5-HT2C signaling pathways and the development of alternative antipsychotic medications.

Biased Ligands at the Kappa Opioid Receptor: Fine-Tuning Receptor Pharmacology.

The kappa opioid receptor (KOR) is a G protein-coupled receptor (GPCR) that can signal through multiple signaling pathways. KOR agonists are known to relieve pain and itch, as well as induce dysphoria, sedation, hallucinations, and diuresis. As is the case with many other GPCRs, specific signaling pathways downstream of the KOR have been linked to certain physiological responses induced by the receptor. Those studies motivated the search and discovery of a number of KOR ligands that preferentially activate one signaling pathway over another. Such compounds are termed functionally selective or biased ligands, and may present a way of inducing desired receptor effects with reduced adverse reactions. In this chapter, I review the molecular intricacies of KOR signaling and discuss the studies that have used biased signaling through the KOR as a way to selectively modulate in vivo physiology.

Conformational plasticity of the intracellular cavity of GPCR-G-protein complexes leads to G-protein promiscuity and selectivity.

While the dynamics of the intracellular surface in agonist-stimulated GPCRs is well studied, the impact of GPCR dynamics on G-protein selectivity remains unclear. Here, we combine molecular dynamics simulations with live-cell FRET and secondary messenger measurements, for 21 GPCR-G-protein combinations, to advance a dynamic model of the GPCR-G-protein interface. Our data show C terminus peptides of Gαs, Gαi, and Gαq proteins assume a small ensemble of unique orientations when coupled to their cognate GPCRs, similar to the variations observed in 3D structures of GPCR-G-protein complexes. The noncognate G proteins interface with latent intracellular GPCR cavities but dissociate due to weak and unstable interactions. Three predicted mutations in ß2-adrenergic receptor stabilize binding of noncognate Gαq protein in its latent cavity, allowing promiscuous signaling through both Gαs and Gαq in a dose-dependent manner. This demonstrates that latent GPCR cavities can be evolved, by design or nature, to tune G-protein selectivity, giving insights to pluridimensional GPCR signaling.

In Silico Study of Allosteric Communication Networks in GPCR Signaling Bias.

Signaling bias is a promising characteristic of G protein-coupled receptors (GPCRs) as it provides the opportunity to develop more efficacious and safer drugs. This is because biased ligands can avoid the activation of pathways linked to side effects whilst still producing the desired therapeutic effect. In this respect, a deeper understanding of receptor dynamics and implicated allosteric communication networks in signaling bias can accelerate the research on novel biased drug candidates. In this review, we aim to provide an overview of computational methods and techniques for studying allosteric communication and signaling bias in GPCRs. This includes (i) the detection of allosteric communication networks and (ii) the application of network theory for extracting relevant information pipelines and highly communicated sites in GPCRs. We focus on the most recent research and highlight structural insights obtained based on the framework of allosteric communication networks and network theory for GPCR signaling bias.

Functional Selectivity of Coumarin Derivates Acting via GPR55 in Neuroinflammation.

Anti-neuroinflammatory treatment has gained importance in the search for pharmacological treatments of different neurological and psychiatric diseases, such as depression, schizophrenia, Parkinson's disease, and Alzheimer's disease. Clinical studies demonstrate a reduction of the mentioned diseases' symptoms after the administration of anti-inflammatory drugs. Novel coumarin derivates have been shown to elicit anti-neuroinflammatory effects via G-protein coupled receptor GPR55, with possibly reduced side-effects compared to the known anti-inflammatory drugs. In this study, we, therefore, evaluated the anti-inflammatory capacities of the two novel coumarin-based compounds, KIT C and KIT H, in human neuroblastoma cells and primary murine microglia. Both compounds reduced PGE2-concentrations likely via the inhibition of COX-2 synthesis in SK-N-SH cells but only KIT C decreased PGE2-levels in primary microglia. The examination of other pro- and anti-inflammatory parameters showed varying effects of both compounds. Therefore, the differences in the effects of KIT C and KIT H might be explained by functional selectivity as well as tissue- or cell-dependent expression and signal pathways coupled to GPR55. Understanding the role of chemical residues in functional selectivity and specific cell- and tissue-targeting might open new therapeutic options in pharmacological drug development and might improve the treatment of the mentioned diseases by intervening in an early step of their pathogenesis.