Generic GPCR residue numbers - aligning topology maps while minding the gaps

AlloViz: A tool for the calculation and visualisation of protein allosteric communication networks

Allostery, the presence of functional interactions between distant parts of proteins, is a critical concept in the field of biochemistry and molecular biology, particularly in the context of protein function and regulation. Understanding the principles of allosteric regulation is essential for advancing our knowledge of biology and developing new therapeutic strategies. This paper presents AlloViz, an open-source Python package designed to quantitatively determine, analyse, and visually represent allosteric communication networks on the basis of molecular dynamics (MD) simulation data. The software integrates well-known techniques for understanding allosteric properties simplifying the process of accessing, rationalising, and representing protein allostery and communication routes. It overcomes the inefficiency of having multiple methods with heterogeneous implementations and showcases the advantages of using MD simulations and multiple replicas to obtain statistically sound information on protein dynamics; it also enables the calculation of "consensus-like" scores aggregating methods that consider multiple structural aspects of allosteric networks. We demonstrate the features of AlloViz on two proteins: β-arrestin 1, a key player for regulating G protein-coupled receptor (GPCR) signalling, and the protein tyrosine phosphatase 1B, an important pharmaceutical target for allosteric inhibitors. The software includes comprehensive documentation and examples, tutorials, and a user-friendly graphical interface.

5-HT6 receptor neutral antagonists protect astrocytes: A lesson from 2-phenylpyrrole derivatives

The serotonin type 6 receptor (5-HT6R) displays a strong constitutive activity, suggesting it participates largely in the physiological and pathological processes controlled by the receptor. The active states of 5-HT6R engage particular signal transduction pathways that lead to different biological responses. In this study, we present the development of 5-HT6R neutral antagonists at Gs signaling built upon the 2-phenylpyrrole scaffold. Using molecular dynamics simulations, we outline the relationship between the exposure of the basic center of the molecules and their ability to target the agonist-activated state of the receptor. Our study identifies compound 30 as a potent and selective neutral antagonist at 5-HT6R-operated Gs signaling. Furthermore, we demonstrate the cytoprotective effects of 30 and structurally diverse 5-HT6R neutral antagonists at Gs signaling in C8-D1A cells and human astrocytes exposed to rotenone. This effect is not observed for 5-HT6R agonists or inverse agonists. In light of these findings, we propose compound 30 as a valuable molecular probe to study the biological effects associated with the agonist-activated state of 5-HT6R and provide insight into the glioprotective properties of 5-HT6R neutral antagonists at Gs signaling.

Structure-guided discovery of orexin receptor-binding PET ligands

Molecular imaging using positron emission tomography (PET) can serve as a promising tool for visualizing biological targets in the brain. Insights into the expression pattern and the in vivo imaging of the G protein-coupled orexin receptors OX1R and OX2R will further our understanding of the orexin system and its role in various physiological and pathophysiological processes. Guided by crystal structures of our lead compound JH112 and the approved hypnotic drug suvorexant bound to OX1R and OX2R, respectively, we herein describe the design and synthesis of two novel radioligands, [18F]KD23 and [18F]KD10. Key to the success of our structural modifications was a bioisosteric replacement of the triazole moiety with a fluorophenyl group. The 19F-substituted analog KD23 showed high affinity for the OX1R and selectivity over OX2R, while the high affinity ligand KD10 displayed similar Ki values for both subtypes. Radiolabeling starting from the respective pinacol ester precursors resulted in excellent radiochemical yields of 93% and 88% for [18F]KD23 and [18F]KD10, respectively, within 20 min. The new compounds will be useful in PET studies aimed at subtype-selective imaging of orexin receptors in brain tissue.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

The versatile binding landscape of the TAAR1 pocket for LSD and other antipsychotic drug molecules

Increasing global concerns about psychoactive substance addiction and psychotic disorders highlight the need for comprehensive research into the structure-function relationship governing ligand recognition between these substances and their receptors in the brain. Recent studies indicate the significant involvement of trace amine-associated receptor 1 (TAAR1) in the signaling regulation of the hallucinogen lysergic acid diethylamide (LSD) and other antipsychotic drugs. This study presents structures of the TAAR1-Gs protein complex recognizing LSD, which exhibits a polypharmacological profile, and the partial agonist RO5263397, which is a drug candidate for schizophrenia and addiction. Moreover, we elucidate the cross-species recognition and partial activation mechanism for TAAR1, which holds promising implications from a drug discovery perspective. Through mutagenesis, functional studies, and molecular dynamics (MD) simulations, we provide a comprehensive understanding of a versatile TAAR1 pocket in recognizing various ligands as well as in the ligand-free state, underpinning the structural basis of its high adaptability. These findings offer valuable insights for the design of antipsychotic drugs.

Relevance of G protein-coupled receptor (GPCR) dynamics for receptor activation, signalling bias and allosteric modulation

G protein-coupled receptors (GPCRs) are one of the major drug targets. In recent years, computational drug design for GPCRs has mainly focused on static structures obtained through X-ray crystallography, cryogenic electron microscopy (cryo-EM) or in silico modelling as a starting point for virtual screening campaigns. However, GPCRs are highly flexible entities with the ability to adopt different conformational states that elicit different physiological responses. Including this knowledge in the drug discovery pipeline can help to tailor novel conformation-specific drugs with an improved therapeutic profile. In this review, we outline our current knowledge about GPCR dynamics that is relevant for receptor activation, signalling bias and allosteric modulation. Ultimately, we highlight new technological implementations such as time-resolved X-ray crystallography and cryo-EM as well as computational algorithms that can contribute to a more comprehensive understanding of receptor dynamics and its relevance for GPCR functionality.

Single nucleotide variations encoding missense mutations in G protein-coupled receptors may contribute to autism

Autism is a neurodevelopmental condition with a range of symptoms that vary in intensity and severity from person to person. Genetic sequencing has identified thousands of genes containing mutations in autistic individuals, which may contribute to the development of autistic symptoms. Several of these genes encode G protein-coupled receptors (GPCRs), which are cell surface expressed proteins that transduce extracellular messages to the intracellular space. Mutations in GPCRs can impact their function, resulting in aberrant signalling within cells and across neurotransmitter systems in the brain. This review summarises the current knowledge on autism-associated single nucleotide variations encoding missense mutations in GPCRs and the impact of these genetic mutations on GPCR function. For some autism-associated mutations, changes in GPCR expression levels, ligand affinity, potency and efficacy have been observed. However, for many the functional consequences remain unknown. Thus, further work to characterise the functional impacts of the genetically identified mutations is required. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Discovery and characterization of small-molecule TGR5 ligands with agonistic activity

The Takeda G protein-coupled receptor 5 (TGR5) is activated endogenously by primary and secondary bile acids. This receptor is considered a candidate target for addressing inflammatory and metabolic disorders. We have targeted TGR5 with structure-based methods for ligand finding using the recently solved experimental structures, as well as structures obtained from molecular dynamics simulations. Through addressing the orthosteric as well as a putative allosteric site, we identified agonists and positive allosteric modulators. While the predicted binding locations were not in line with their efficacy, our work contributes activating small-molecule ligands that we have thoroughly characterized in vitro.

AlphaFold2 structures guide prospective ligand discovery

AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.

New Mechanism for the Apoptosis of Human Neuroblastoma Cells by the Interaction between Fluorene-9-Bisphenol and the G Protein-Coupled Estrogen Receptor 1

Fluorene-9-bisphenol (BHPF) is an emerging contaminant. Presently, there is no report on its interaction with G protein-coupled estrogen receptor 1 (GPER). By using an integrated toxicity research scenario that combined theoretical study with experimental methods, BHPF was found to inhibit the GPER-mediated effect via direct receptor binding. Molecular dynamics simulations found that Trp2726.48 and Glu2756.51 be the key amino acids of BHPF binding with GPER. Moreover, the calculation indicated that BHPF was a suspected GPER inhibitor, which neither can activate GPER nor is able to form water channels of GPER. The role of two residues was successfully verified by following gene knockout and site-directed mutagenesis assays. Further in vitro assays showed that BHPF could attenuate the increase in intracellular concentration of free Ca2+ induced by G1-activated GPER. Besides, BHPF showed an enhanced cytotoxicity compared with G15, indicating that BHPF might be a more potent GPER inhibitor than G15. In addition, a statistically significant effect on the mRNA level of GPER was observed for BHPF. In brief, the present study proposes that BHPF be a GPER inhibitor, and its GPER molecular recognition mechanism has been revealed, which is of great significance for the health risk and assessment of BHPF.

Intracellular Allosteric Antagonist of the Olfactory Receptor OR51E2

Olfactory receptors are members of class A (rhodopsin-like) family of G protein-coupled receptors (GPCRs). Their expression and function have been increasingly studied in nonolfactory tissues, and many have been identified as potential therapeutic targets. In this manuscript, we focus on the discovery of novel ligands for the olfactory receptor family 51 subfamily E2 (OR51E2). We performed an artificial intelligence-based virtual drug screen of a ∼2.2 million small molecule library. Cell-based functional assay identified compound 80 (C80) as an antagonist and inverse agonist, and detailed pharmacological analysis revealed C80 acts as a negative allosteric modulator by significantly decreasing the agonist efficacy, while having a minimal effect on receptor affinity for agonist. C80 binds to an allosteric binding site formed by a network of nine residues localized in the intracellular parts of transmembrane domains 3, 5, 6, 7, and H8, which also partially overlaps with a G protein binding site. Mutational experiments of residues involved in C80 binding uncovered the significance of the C2406.37 position in blocking the activation-related conformational change and keeping the receptor in the inactive form. Our study provides a mechanistic understanding of the negative allosteric action of C80 on agonist-ctivated OR51E2. We believe the identification of the antagonist of OR51E2 will enable a multitude of studies aiming to determine the functional role of this receptor in specific biologic processes. SIGNIFICANCE STATEMENT: OR51E2 has been implicated in various biological processes, and its antagonists that can effectively modulate its activity have therapeutic potential. Here we report the discovery of a negative allosteric modulator of OR51E2 and provide a mechanistic understanding of its action. We demonstrate that this modulator has an inhibitory effect on the efficacy of the agonist for the receptor and reveal a network of nine residues that constitute its binding pocket, which also partially overlaps with the G protein binding site.

Benchmarking AlphaFold-Generated Structures of Chemokine-Chemokine Receptor Complexes

AlphaFold and AlphaFold-Multimer have become two essential tools for the modeling of unknown structures of proteins and protein complexes. In this work, we extensively benchmarked the quality of chemokine-chemokine receptor structures generated by AlphaFold-Multimer against experimentally determined structures. Our analysis considered both the global quality of the model, as well as key structural features for chemokine recognition. To study the effects of template and multiple sequence alignment parameters on the results, a new prediction pipeline called LIT-AlphaFold (https://github.com/LIT-CCM-lab/LIT-AlphaFold) was developed, allowing extensive input customization. AlphaFold-Multimer correctly predicted differences in chemokine binding orientation and accurately reproduced the unique binding orientation of the CXCL12-ACKR3 complex. Further, the predictions of the full receptor N-terminus provided insights into a putative chemokine recognition site 0.5. The accuracy of chemokine N-terminus binding mode prediction varied between complexes, but the confidence score permitted the distinguishing of residues that were very likely well positioned. Finally, we generated a high-confidence model of the unsolved CXCL12-CXCR4 complex, which agreed with experimental mutagenesis and cross-linking data.

Insights on the G protein-coupled receptor helix 8 solution structure and orientation using a neurotensin receptor 1 peptide

G-protein coupled receptors (GPCRs) are the largest class of membrane proteins encoded in the human genome with high pharmaceutical relevance and implications to human health. These receptors share a prevalent architecture of seven transmembrane helices followed by an intracellular, amphipathic helix 8 (H8) and a disordered C-terminal tail (Ctail). Technological advancements have led to over 1000 receptor structures in the last two decades, yet frequently H8 and the Ctail are conformationally heterogeneous or altogether absent. Here we synthesize a peptide comprising the neurotensin receptor 1 (NTS1) H8 and Ctail (H8-Ctail) to investigate its structural stability, conformational dynamics, and orientation in the presence of detergent and phospholipid micelles, which mimic the membrane. Circular dichroism (CD) and nuclear magnetic resonance (NMR) measurements confirm that zwitterionic 1,2-diheptanoyl-sn-glycero-3-phosphocholine is a potent stabilizer of H8 structure, whereas the commonly-used branched detergent lauryl maltose neopentyl glycol (LMNG) is unable to completely stabilize the helix - even at amounts four orders of magnitude greater than its critical micellar concentration. We then used NMR spectroscopy to assign the backbone chemical shifts. A series of temperature and lipid titrations were used to define the H8 boundaries as F376-R392 from chemical shift perturbations, changes in resonance intensity, and chemical-shift-derived phi/psi angles. Finally, the H8 azimuthal and tilt angles, defining the helix orientation relative of the membrane normal were measured using paramagnetic relaxation enhancement NMR. Taken together, our studies reveal the H8-Ctail region is sensitive to membrane physicochemical properties and is capable of more adaptive behavior than previously suggested by static structural techniques.

Structural basis for selectivity and antagonism in extracellular GPCR-nanobodies

G protein-coupled receptors (GPCRs) are pivotal therapeutic targets, but their complex structure poses challenges for effective drug design. Nanobodies, or single-domain antibodies, have emerged as a promising therapeutic strategy to target GPCRs, offering advantages over traditional small molecules and antibodies. However, an incomplete understanding of the structural features enabling GPCR-nanobody interactions has limited their development. In this study, we investigate VUN701, a nanobody antagonist targeting the atypical chemokine receptor 3 (ACKR3). We determine that an extended CDR3 loop is required for ACKR3 binding. Uncommon in most nanobodies, an extended CDR3 is prevalent in GPCR-targeting nanobodies. Combining experimental and computational approaches, we map an inhibitory ACKR3-VUN701 interface and define a distinct conformational mechanism for GPCR inactivation. Our results provide insights into class A GPCR-nanobody selectivity and suggest a strategy for the development of these new therapeutic tools.

Structural basis of Δ9-THC analog activity at the Cannabinoid 1 receptor

Δ9-tetrahydrocannabinol (THC) is the principal psychoactive compound derived from the cannabis plant Cannabis sativa and approved for emetic conditions, appetite stimulation and sleep apnea relief. THC's psychoactive actions are mediated primarily by the cannabinoid receptor CB1. Here, we determine the cryo-EM structure of HU210, a THC analog and widely used tool compound, bound to CB1 and its primary transducer, Gi1. We leverage this structure for docking and 1,000 ns molecular dynamics simulations of THC and 10 structural analogs delineating their spatiotemporal interactions at the molecular level. Furthermore, we pharmacologically profile their recruitment of Gi and β-arrestins and reversibility of binding from an active complex. By combining detailed CB1 structural information with molecular models and signaling data we uncover the differential spatiotemporal interactions these ligands make to receptors governing potency, efficacy, bias and kinetics. This may help explain the actions of abused substances, advance fundamental receptor activation studies and design better medicines.

Stepwise activation of a metabotropic glutamate receptor

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain that is linked via a cysteine-rich domain to their 7-transmembrane domain1. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the extracellular ligand-binding domain to the G protein-coupling 7-transmembrane domain2. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging, we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

Promiscuous G-protein activation by the calcium-sensing receptor

The human calcium-sensing receptor (CaSR) detects fluctuations in the extracellular Ca2+ concentration and maintains Ca2+ homeostasis1,2. It also mediates diverse cellular processes not associated with Ca2+ balance3-5. The functional pleiotropy of CaSR arises in part from its ability to signal through several G-protein subtypes6. We determined structures of CaSR in complex with G proteins from three different subfamilies: Gq, Gi and Gs. We found that the homodimeric CaSR of each complex couples to a single G protein through a common mode. This involves the C-terminal helix of each Gα subunit binding to a shallow pocket that is formed in one CaSR subunit by all three intracellular loops (ICL1-ICL3), an extended transmembrane helix 3 and an ordered C-terminal region. G-protein binding expands the transmembrane dimer interface, which is further stabilized by phospholipid. The restraint imposed by the receptor dimer, in combination with ICL2, enables G-protein activation by facilitating conformational transition of Gα. We identified a single Gα residue that determines Gq and Gs versus Gi selectivity. The length and flexibility of ICL2 allows CaSR to bind all three Gα subtypes, thereby conferring capacity for promiscuous G-protein coupling.

Venom-inspired somatostatin receptor 4 (SSTR4) agonists as new drug leads for peripheral pain conditions

Persistent pain affects one in five people worldwide, often with severely debilitating consequences. Current treatment options, which can be effective for mild or acute pain, are ill-suited for moderate-to-severe persistent pain, resulting in an urgent need for new therapeutics. In recent years, the somatostatin receptor 4 (SSTR 4 ), which is expressed in sensory neurons of the peripheral nervous system, has emerged as a promising target for pain relief. However, the presence of several closely related receptors with similar ligand-binding surfaces complicates the design of receptor-specific agonists. In this study, we report the discovery of a potent and selective SSTR 4 peptide, consomatin Fj1, derived from extensive venom gene datasets from marine cone snails. Consomatin Fj1 is a mimetic of the endogenous hormone somatostatin and contains a minimized binding motif that provides stability and drives peptide selectivity. Peripheral administration of synthetic consomatin Fj1 provided analgesia in mouse models of postoperative and neuropathic pain. Using structure-activity studies, we designed and functionally evaluated several Fj1 analogs, resulting in compounds with improved potency and selectivity. Our findings present a novel avenue for addressing persistent pain through the design of venom-inspired SSTR 4 -selective pain therapeutics.

One Sentence Summary

Venom peptides from predatory marine mollusks provide new leads for treating peripheral pain conditions through a non-opioid target.

Calcium-Driven In Silico Inactivation of a Human Olfactory Receptor

Conformational changes as well as molecular determinants related to the activation and inactivation of olfactory receptors are still poorly understood due to the intrinsic difficulties in the structural determination of this GPCR family. Here, we perform, for the first time, the in silico inactivation of human olfactory receptor OR51E2, highlighting the possible role of calcium in this receptor state transition. Using molecular dynamics simulations, we show that a divalent ion in the ion binding site, coordinated by two acidic residues at positions 2.50 and 3.39 conserved across most ORs, stabilizes the receptor in its inactive state. In contrast, protonation of the same two acidic residues is not sufficient to drive inactivation within the microsecond timescale of our simulations. Our findings suggest a novel molecular mechanism for OR inactivation, potentially guiding experimental validation and offering insights into the possible broader role of divalent ions in GPCR signaling.

Structural insights into the unexpected agonism of tetracyclic antidepressants at serotonin receptors 5-HT1eR and 5-HT1FR

Serotonin [5-hydroxytryptamine (5-HT)] acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT1eR have confirmed roles in native tissue and are validated drug targets. Despite 5-HT1eR's therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT1eR's pharmacology in relation to the highly homologous 5-HT1FR, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT1eR/5-HT1FR agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed antimigraine properties. Using cryo-EM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT1eR distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT1eR and 5-HT1FR contribute to the agonist activity of these antidepressants.

Cardiotoxic effects of common and emerging drugs: role of cannabinoid receptors

Drug-induced cardiotoxicity has become one of the most common and detrimental health concerns, which causes significant loss to public health and drug resources. Cannabinoid receptors (CBRs) have recently achieved great attention for their vital roles in the regulation of heart health and disease, with mounting evidence linking CBRs with the pathogenesis and progression of drug-induced cardiotoxicity. This review aims to summarize fundamental characteristics of two well-documented CBRs (CB1R and CB2R) from aspects of molecular structure, signaling and their functions in cardiovascular physiology and pathophysiology. Moreover, we describe the roles of CB1R and CB2R in the occurrence of cardiotoxicity induced by common drugs such as antipsychotics, anti-cancer drugs, marijuana, and some emerging synthetic cannabinoids. We highlight the 'yin-yang' relationship between CB1R and CB2R in drug-induced cardiotoxicity and propose future perspectives for CBR-based translational medicine toward cardiotoxicity curation and clinical monitoring.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

AlphaFold2 structures template ligand discovery

AlphaFold2 (AF2) and RosettaFold have greatly expanded the number of structures available for structure-based ligand discovery, even though retrospective studies have cast doubt on their direct usefulness for that goal. Here, we tested unrefined AF2 models prospectively, comparing experimental hit-rates and affinities from large library docking against AF2 models vs the same screens targeting experimental structures of the same receptors. In retrospective docking screens against the σ2 and the 5-HT2A receptors, the AF2 structures struggled to recapitulate ligands that we had previously found docking against the receptors' experimental structures, consistent with published results. Prospective large library docking against the AF2 models, however, yielded similar hit rates for both receptors versus docking against experimentally-derived structures; hundreds of molecules were prioritized and tested against each model and each structure of each receptor. The success of the AF2 models was achieved despite differences in orthosteric pocket residue conformations for both targets versus the experimental structures. Intriguingly, against the 5-HT2A receptor the most potent, subtype-selective agonists were discovered via docking against the AF2 model, not the experimental structure. To understand this from a molecular perspective, a cryoEM structure was determined for one of the more potent and selective ligands to emerge from docking against the AF2 model of the 5-HT2A receptor. Our findings suggest that AF2 models may sample conformations that are relevant for ligand discovery, much extending the domain of applicability of structure-based ligand discovery.

Identification of oleic acid as an endogenous ligand of GPR3

Although GPR3 plays pivotal roles in both the nervous system and metabolic processes, such as cold-induced thermogenesis, its endogenous ligand remains elusive. Here, by combining structural approach (including cryo-electron microscopy), mass spectrometry analysis, and functional studies, we identify oleic acid (OA) as an endogenous ligand of GPR3. Our study reveals a hydrophobic tunnel within GPR3 that connects the extracellular side of the receptor to the middle of plasma membrane, enabling fatty acids to readily engage the receptor. Functional studies demonstrate that OA triggers downstream Gs signaling, whereas lysophospholipids fail to activate the receptor. Moreover, our research reveals that cold stimulation induces the secretion of OA in mice, subsequently activating Gs/cAMP/PKA signaling in brown adipose tissue. Notably, brown adipose tissues from Gpr3 knockout mice do not respond to OA during cold stimulation, reinforcing the significance of GPR3 in this process. Finally, we propose a "born to be activated and cold to enhance" model for GPR3 activation. Our study provides a starting framework for the understanding of GPR3 signaling in cold-stimulated thermogenesis.

GPR101: Modeling a constitutively active receptor linked to X-linked acrogigantism

GPR101 is a G protein-coupled receptor (GPCR) implicated in a rare form of genetic gigantism known as X-linked acrogigantism, or X-LAG. In particular, X-LAG patients harbor microduplications in the long arm of the X-chromosome that invariably include the GPR101 gene. Duplications of the GPR101 gene lead to the formation of a new chromatin domain that causes over-expression of the receptor in the pituitary tumors of the patients. Notably, GPR101 is a constitutively active receptor, which stimulates cells to produce the second messenger cyclic AMP (cAMP) in the absence of ligands. Moreover, GPR101 was recently reported to constitutively activate not only the cAMP pathway via Gs, but also other G protein subunits (Gq/11 and G12/13). Hence, chemicals that block the constitutive activity of GPR101, known as inverse agonists, have the potential to be useful for the development of pharmacological tools for the treatment of X-LAG. In this study, we provide structural insights into the putative structure of GPR101 based on in-house built homology models, as well as third party models based on the machine learning methods AlphaFold and AlphaFold-Multistate. Moreover, we report a molecular dynamics study, meant to further probe the constitutive activity of GPR101. Finally, we provide a structural comparison with the closest GPCRs, which suggests that GPR101 does not share their natural ligands. While this manuscript was under review, cryo-electron microscopy structures of GPR101 were reported. These structures are expected to enable computer-aided ligand discovery efforts targeting GPR101.

Classifying tetraspanins: A universal system for numbering residues and a proposal for naming structural motifs and subfamilies

All tetraspanins have four transmembrane domains (TMs). The large extracellular loop (LEL) that connects the third and fourth TMs contains multiple secondary structures together with the family's signature Cys-Cys-Gly motif. These intriguing membrane proteins are involved in diverse and incompletely understood cellular processes including cell adhesion, tissue differentiation, immune cell maturation and host-parasite interactions. Here we present a classification system that accurately describes the position of each amino acid within its primary sequence based on both sequence and topological conservation of the TMs and LEL. This builds on the numbering systems that have been used in the G protein-coupled receptor (GPCR) field for nearly three decades and which have aided the understanding of GPCR structure/activity relationships and ligand interactions. The high-resolution structures of the tetraspanins CD81, CD9, CD53 and Tspan15 were used to validate the structural relevance of our new tetraspanin classification system. Modelling of all tetraspanin LELs highlighted flexibility in LEL disulfide bonding across the family and suggests that the structural arrangement of tetraspanin LELs is more complex than previously thought. We therefore propose a new subfamily naming system that addresses this added complexity and facilitates the systematic classification of human tetraspanins, shedding light on all structural motifs within the family. We anticipate that our universal tetraspanin classification system will enable progress in defining how sequence and structure inform function.

Nanobody-Mediated Dualsteric Engagement of the Angiotensin Receptor Broadens Biased Ligand Pharmacology

Dualsteric G protein-coupled receptor (GPCR) ligands are a class of bitopic ligands that consist of an orthosteric pharmacophore, which binds to the pocket occupied by the receptor's endogenous agonist, and an allosteric pharmacophore, which binds to a distinct site. These ligands have the potential to display characteristics of both orthosteric and allosteric ligands. To explore the signaling profiles that dualsteric ligands of the angiotensin II type 1 receptor (AT1R) can access, we ligated a 6e epitope tag-specific nanobody (single-domain antibody fragment) to angiotensin II (AngII) and analogs that show preferential allosteric coupling to Gq (TRV055, TRV056) or β-arrestin (TRV027). While the nanobody itself acts as a probe-specific neutral or negative allosteric ligand of N-terminally 6e-tagged AT1R, nanobody conjugation to orthosteric ligands had varying effects on Gq dissociation and β-arrestin plasma membrane recruitment. The potency of certain AngII analogs was enhanced up to 100-fold, and some conjugates behaved as partial agonists, with up to a 5-fold decrease in maximal efficacy. Nanobody conjugation also biased the signaling of TRV055 and TRV056 toward Gq, suggesting that Gq bias at AT1R can be modulated through molecular mechanisms distinct from those previously elucidated. Both competition radioligand binding experiments and functional assays demonstrated that orthosteric antagonists (angiotensin receptor blockers) act as non-competitive inhibitors of all these nanobody-peptide conjugates. This proof-of-principle study illustrates the array of pharmacological patterns that can be achieved by incorporating neutral or negative allosteric pharmacophores into dualsteric ligands. Nanobodies directed toward linear epitopes could provide a rich source of allosteric reagents for this purpose. SIGNIFICANCE STATEMENT: Here we engineer bitopic (dualsteric) ligands for epitope-tagged angiotensin II type 1 receptor by conjugating angiotensin II or its biased analogs to an epitope-specific nanobody (antibody fragment). Our data demonstrate that nanobody-mediated interactions with the receptor N-terminus endow angiotensin analogs with properties of allosteric modulators and provide a novel mechanism to increase the potency, modulate the maximal effect, or alter the bias of ligands.

Pharmacological Characterization of the Zebrafish (Danio Rerio) Histamine H1 Receptor Reveals the Involvement of the Second Extracellular Loop in the Binding of Histamine

The zebrafish (Danio rerio) histamine H1 receptor gene (zfH1R) was cloned in 2007 and reported to be involved in fish locomotion. Yet, no detailed characterization of its pharmacology and signaling properties have so far been reported. In this study, we pharmacologically characterized the zfH1R expressed in HEK-293T cells by means of [3H]-mepyramine binding and G protein-signaling assays. The zfH1R [dissociation constant (KD), 0.7 nM] displayed similar affinity for the antagonist [3H]-mepyramine as the human histamine H1 receptor (hH1R) (KD, 1.5 nM), whereas the affinity for histamine is 100-fold higher than for the human H1R. The zfH1R couples to Gαq/11 proteins and activates several reporter genes, i.e., NFAT, NFϰB, CRE, VEGF, COX-2, SRE, and AP-1, and zfH1R-mediated signaling is prevented by the Gαq/11 inhibitor YM-254890 and the antagonist mepyramine. Molecular modeling of the zfH1R and human H1R shows that the binding pockets are identical, implying that variations along the ligand binding pathway could underly the differences in histamine affinity instead. Targeting differentially charged residues in extracellular loop 2 (ECL2) using site-directed mutagenesis revealed that Arg21045x55 is most likely involved in the binding process of histamine in zfH1R. This study aids the understanding of the pharmacological differences between H1R orthologs and the role of ECL2 in histamine binding and provides fundamental information for the understanding of the histaminergic system in the zebrafish. SIGNIFICANCE STATEMENT: The use of the zebrafish as in vivo models in neuroscience is growing exponentially, which asks for detailed characterization of the aminergic neurotransmitter systems in this model. This study is the first to pharmacologically characterize the zebrafish histamine H1 receptor after expression in HEK-293T cells. The results show a high pharmacological and functional resemblance with the human ortholog but also reveal interesting structural differences and unveils an important role of the second extracellular loop in histamine binding.

The Pivotal Distinction between Antagonists' and Agonists' Binding into Dopamine D4 Receptor-MD and FMO/PIEDA Studies

The dopamine D4 receptor (D4R) is a promising therapeutic target in widespread diseases, and the search for novel agonists and antagonists appears to be clinically relevant. The mechanism of binding to the receptor (R) for antagonists and agonists varies. In the present study, we conducted an in-depth computational study, teasing out key similarities and differences in binding modes, complex dynamics, and binding energies for D4R agonists and antagonists. The dynamic network method was applied to investigate the communication paths between the ligand (L) and G-protein binding site (GBS) of human D4R. Finally, the fragment molecular orbitals with pair interaction energy decomposition analysis (FMO/PIEDA) scheme was used to estimate the binding energies of L-R complexes. We found that a strong salt bridge with D3.32 initiates the inhibition of the dopamine D4 receptor. This interaction also occurs in the binding of agonists, but the change in the receptor conformation to the active state starts with interaction with cysteine C3.36. Such a mechanism may arise in the case of agonists unable to form a hydrogen bond with the serine S5.46, considered, so far, to be crucial in the activation of GPCRs. The energy calculations using the FMO/PIEDA method indicate that antagonists show higher residue occupancy of the receptor binding site than agonists, suggesting they could form relatively more stable complexes. Additionally, antagonists were characterized by repulsive interactions with S5.46 distinguishing them from agonists.

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.

Deep learning-based design and screening of benzimidazole-pyrazine derivatives as adenosine A2B receptor antagonists

The Adenosine A2B receptor (A2BAR) is considered a novel potential target for the immunotherapy of cancer, and A2BAR antagonists have an inhibitory effect on tumor growth, proliferation, and metastasis. In our previous studies, we identified a class of benzimidazole-pyrazine scaffolds whose derivatives exhibited the antagonistic effect but lacked subtype selectivity towards A2BAR. In this work, we developed a scaffold-based protocol that incorporates a deep generative model and multilayer virtual screening to design benzimidazole-pyrazine derivatives as potential selective A2BAR antagonists. By utilizing a generative model with reported A2BAR antagonists as the training set, we built up a scaffold-focused library of benzimidazole-pyrazine derivatives and processed a virtual screening protocol to discover potential A2BAR antagonists. Finally, five molecules with different Bemis-Murcko scaffolds were identified and exhibited higher binding free energies than the reference molecule 12o. Further computational analysis revealed that the 3-benzyl derivative ABA-1266 presented high selectivity toward A2BAR and showed preferred draggability, providing future potent development of selective A2BAR antagonists.Communicated by Ramaswamy H. Sarma.

Molecular determinants of ligand efficacy and potency in GPCR signaling

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) bind to extracellular ligands and drugs and modulate intracellular responses through conformational changes. Despite their importance as drug targets, the molecular origins of pharmacological properties such as efficacy (maximum signaling response) and potency (the ligand concentration at half-maximal response) remain poorly understood for any ligand-receptor-signaling system. We used the prototypical adrenaline-β2 adrenergic receptor-G protein system to reveal how specific receptor residues decode and translate the information encoded in a ligand to mediate a signaling response. We present a data science framework to integrate pharmacological and structural data to uncover structural changes and allosteric networks relevant for ligand pharmacology. These methods can be tailored to study any ligand-receptor-signaling system, and the principles open possibilities for designing orthosteric and allosteric compounds with defined signaling properties.

Snapshot of the cannabinoid receptor 1-arrestin complex unravels the biased signaling mechanism

Cannabis activates the cannabinoid receptor 1 (CB1), which elicits analgesic and emotion regulation benefits, along with adverse effects, via Gi and β-arrestin signaling pathways. However, the lack of understanding of the mechanism of β-arrestin-1 (βarr1) coupling and signaling bias has hindered drug development targeting CB1. Here, we present the high-resolution cryo-electron microscopy structure of CB1-βarr1 complex bound to the synthetic cannabinoid MDMB-Fubinaca (FUB), revealing notable differences in the transducer pocket and ligand-binding site compared with the Gi protein complex. βarr1 occupies a wider transducer pocket promoting substantial outward movement of the TM6 and distinctive twin toggle switch rearrangements, whereas FUB adopts a different pose, inserting more deeply than the Gi-coupled state, suggesting the allosteric correlation between the orthosteric binding pocket and the partner protein site. Taken together, our findings unravel the molecular mechanism of signaling bias toward CB1, facilitating the development of CB1 agonists.

Understanding the molecular mechanisms of odorant binding and activation of the human OR52 family

Structural and mechanistic studies on human odorant receptors (ORs), key in olfactory signaling, are challenging because of their low surface expression in heterologous cells. The recent structure of OR51E2 bound to propionate provided molecular insight into odorant recognition, but the lack of an inactive OR structure limited understanding of the activation mechanism of ORs upon odorant binding. Here, we determined the cryo-electron microscopy structures of consensus OR52 (OR52cs), a representative of the OR52 family, in the ligand-free (apo) and octanoate-bound states. The apo structure of OR52cs reveals a large opening between transmembrane helices (TMs) 5 and 6. A comparison between the apo and active structures of OR52cs demonstrates the inward and outward movements of the extracellular and intracellular segments of TM6, respectively. These results, combined with molecular dynamics simulations and signaling assays, shed light on the molecular mechanisms of odorant binding and activation of the OR52 family.

Recognition of methamphetamine and other amines by trace amine receptor TAAR1

Trace amine-associated receptor 1 (TAAR1), the founding member of a nine-member family of trace amine receptors, is responsible for recognizing a range of biogenic amines in the brain, including the endogenous β-phenylethylamine (β-PEA)1 as well as methamphetamine2, an abused substance that has posed a severe threat to human health and society3. Given its unique physiological role in the brain, TAAR1 is also an emerging target for a range of neurological disorders including schizophrenia, depression and drug addiction2,4,5. Here we report structures of human TAAR1-G-protein complexes bound to methamphetamine and β-PEA as well as complexes bound to RO5256390, a TAAR1-selective agonist, and SEP-363856, a clinical-stage dual agonist for TAAR1 and serotonin receptor 5-HT1AR (refs. 6,7). Together with systematic mutagenesis and functional studies, the structures reveal the molecular basis of methamphetamine recognition and underlying mechanisms of ligand selectivity and polypharmacology between TAAR1 and other monoamine receptors. We identify a lid-like extracellular loop 2 helix/loop structure and a hydrogen-bonding network in the ligand-binding pockets, which may contribute to the ligand recognition in TAAR1. These findings shed light on the ligand recognition mode and activation mechanism for TAAR1 and should guide the development of next-generation therapeutics for drug addiction and various neurological disorders.

Ligand selectivity hotspots in serotonin GPCRs

Serotonin is a neurotransmitter regulating numerous physiological processes also modulated by drugs, for example, schizophrenia, depression, migraine, and obesity. However, these drugs typically have adverse effects caused by promiscuous binding across 12 serotonin and more than 20 homologous receptors. Recently, structures of the entire serotonin receptor family uncovered molecular ligand recognition. Here, we present a map of 19 'selectivity hotspots', that is, nonconserved binding site residues governing selectivity via favorable target interactions or repulsive 'off-target' contacts. Furthermore, we review functional rationale from observed ligand-binding affinities and mutagenesis effects. Unifying knowledge underlying specific probes and drugs is critical toward the functional characterization of different receptors and alleviation of adverse effects.

Modeling the structure of the transmembrane domain of T1R3, a subunit of the sweet taste receptor, with neohesperidin dihydrochalcone using molecular dynamics simulation

Neohesperidin dihydrochalcone (NHDC) is a sweetener, which interacts with the transmembrane domain (TMD) of the T1R3 subunit of the human sweet taste receptor. Although NHDC and a sweet taste inhibitor lactisole share similar structural motifs, they have opposite effects on the receptor. This study involved the creation of an NHDC-docked model of T1R3 TMD through mutational analyses followed by in silico simulations. When certain NHDC derivatives were docked to the model, His7345.44 was demonstrated to play a crucial role in activating T1R3 TMD. The NHDC-docked model was then compared with a lactisole-docked inactive form, several residues were characterized as important for the recognition of NHDC; however, most of them were distinct from those of lactisole. Residues such as His6413.33 and Gln7947.38 were found to be oriented differently. This study provides useful information that will facilitate the design of sweeteners and inhibitors that interact with T1R3 TMD.

Engineered odorant receptors illuminate structural principles of odor discrimination

A central challenge in olfaction is understanding how the olfactory system detects and distinguishes odorants with diverse physicochemical properties and molecular configurations. Vertebrate animals perceive odors via G protein-coupled odorant receptors (ORs). In humans, ~400 ORs enable the sense of smell. The OR family is composed of two major classes: Class I ORs are tuned to carboxylic acids while Class II ORs, representing the vast majority of the human repertoire, respond to a wide variety of odorants. How ORs recognize chemically diverse odorants remains poorly understood. A fundamental bottleneck is the inability to visualize odorant binding to ORs. Here, we uncover fundamental molecular properties of odorant-OR interactions by employing engineered ORs crafted using a consensus protein design strategy. Because such consensus ORs (consORs) are derived from the 17 major subfamilies of human ORs, they provide a template for modeling individual native ORs with high sequence and structural homology. The biochemical tractability of consORs enabled four cryoEM structures of distinct consORs with unique ligand recognition properties. The structure of a Class I consOR, consOR51, showed high structural similarity to the native human receptor OR51E2 and yielded a homology model of a related member of the human OR51 family with high predictive power. Structures of three Class II consORs revealed distinct modes of odorant-binding and activation mechanisms between Class I and Class II ORs. Thus, the structures of consORs lay the groundwork for understanding molecular recognition of odorants by the OR superfamily.

Structure-Based Discovery of Mouse Trace Amine-Associated Receptor 5 Antagonists

Trace amine-associated receptors (TAARs) were discovered in 2001 as new members of class A G protein-coupled receptors (GPCRs). With the only exception of TAAR1, TAAR members (TAAR2-9, also known as noncanonical olfactory receptors) were originally described exclusively in the olfactory epithelium and believed to mediate the innate perception of volatile amines. However, most noncanonical olfactory receptors are still orphan receptors. Given its recently discovered nonolfactory expression and therapeutic potential, TAAR5 has been the focus of deorphanization campaigns that led to the discovery of a few druglike antagonists. Here, we report four novel TAAR5 antagonists identified through high-throughput screening, which, along with the four ligands published in the literature, constituted our starting point to design a computational strategy for the identification of TAAR5 ligands. We developed a structure-based virtual screening protocol that allowed us to identify three new TAAR5 antagonists with a hit rate of 10%. Despite lacking an experimental structure, we accurately modeled the TAAR5 binding site by integrating comparative sequence- and structure-based analyses of serotonin receptors with homology modeling and side-chain optimization. In summary, we have identified seven new TAAR5 antagonists that could serve as lead candidates for the development of new treatments for depression, anxiety, and neurodegenerative diseases.

Development of Fluorescent AF64394 Analogues Enables Real-Time Binding Studies for the Orphan Class A GPCR GPR3

The orphan G protein-coupled receptor (oGPCR) GPR3 represents a potential drug target for the treatment of Alzheimer's disease and metabolic disorders. However, the limited toolbox of pharmacological assays hampers the development of advanced ligands. Here, we developed a signaling pathway-independent readout of compound-GPR3 interaction. Starting from computational binding pose predictions of the most potent GPR3 ligand, we designed a series of fluorescent AF64394 analogues and assessed their suitability for BRET-based binding studies. The most potent ligand, 45 (UR-MB-355), bound to GPR3 and closely related receptors, GPR6 and GPR12, with similar submicromolar affinities. Furthermore, we found that 45 engages GPR3 in a distinct mode compared to AF64394, and coincubation studies with the GPR3 agonist diphenyleneiodonium chloride revealed allosteric modulation of 45 binding. These insights provide new cues for the pharmacological manipulation of GPR3 activity. This novel binding assay will foster the development of future drugs acting through these pharmacologically attractive oGPCRs.

Side-chain dynamics of the α1B -adrenergic receptor determined by NMR via methyl relaxation

G protein-coupled receptors (GPCRs) are medically important membrane proteins that sample inactive, intermediate, and active conformational states characterized by relatively slow interconversions (~μs-ms). On a faster timescale (~ps-ns), the conformational landscape of GPCRs is governed by the rapid dynamics of amino acid side chains. Such dynamics are essential for protein functions such as ligand recognition and allostery. Unfortunately, technical challenges have almost entirely precluded the study of side-chain dynamics for GPCRs. Here, we investigate the rapid side-chain dynamics of a thermostabilized α1B -adrenergic receptor (α1B -AR) as probed by methyl relaxation. We determined order parameters for Ile, Leu, and Val methyl groups in the presence of inverse agonists that bind orthosterically (prazosin, tamsulosin) or allosterically (conopeptide ρ-TIA). Despite the differences in the ligands, the receptor's overall side-chain dynamics are very similar, including those of the apo form. However, ρ-TIA increases the flexibility of Ile1764×56 and possibly of Ile2145×49 , adjacent to Pro2155×50 of the highly conserved P5×50 I3×40 F6×44 motif crucial for receptor activation, suggesting differences in the mechanisms for orthosteric and allosteric receptor inactivation. Overall, increased Ile side-chain rigidity was found for residues closer to the center of the membrane bilayer, correlating with denser packing and lower protein surface exposure. In contrast to two microbial membrane proteins, in α1B -AR Leu exhibited higher flexibility than Ile side chains on average, correlating with the presence of Leu in less densely packed areas and with higher protein-surface exposure than Ile. Our findings demonstrate the feasibility of studying receptor-wide side-chain dynamics in GPCRs to gain functional insights.

Rules and mechanisms governing G protein coupling selectivity of GPCRs

G protein-coupled receptors (GPCRs) convert extracellular stimuli into intracellular signaling by coupling to heterotrimeric G proteins of four classes: Gi/o, Gq, Gs, and G12/13. However, our understanding of the G protein selectivity of GPCRs is incomplete. Here, we quantitatively measure the enzymatic activity of GPCRs in living cells and reveal the G protein selectivity of 124 GPCRs with the exact rank order of their G protein preference. Using this information, we establish a classification of GPCRs by functional selectivity, discover the existence of a G12/13-coupled receptor, G15-coupled receptors, and a variety of subclasses for Gi/o-, Gq-, and Gs-coupled receptors, culminating in development of the predictive algorithm of G protein selectivity. We further identify the structural determinants of G protein selectivity, allowing us to synthesize non-existent GPCRs with de novo G protein selectivity and efficiently identify putative pathogenic variants.

Reanalysis of a μ opioid receptor crystal structure reveals a covalent adduct with BU72

BACKGROUND

The first crystal structure of the active μ opioid receptor (μOR) exhibited several unexplained features. The ligand BU72 exhibited many extreme deviations from ideal geometry, along with unexplained electron density. I previously showed that inverting the benzylic configuration resolved these problems, establishing revised stereochemistry of BU72 and its analog BU74. However, another problem remains unresolved: additional unexplained electron density contacts both BU72 and a histidine residue in the N-terminus, revealing the presence of an as-yet unidentified atom.

RESULTS

These short contacts and uninterrupted density are inconsistent with non-covalent interactions. Therefore, BU72 and μOR form a covalent adduct, rather than representing two separate entities as in the original model. A subsequently proposed magnesium complex is inconsistent with multiple lines of evidence. However, oxygen fits the unexplained density well. While the structure I propose is tentative, similar adducts have been reported previously in the presence of reactive oxygen species. Moreover, known sources of reactive oxygen species were present: HEPES buffer, nickel ions, and a sequence motif that forms redox-active nickel complexes. This motif contacts the unexplained density. The adduct exhibits severe strain, and the tethered N-terminus forms contacts with adjacent residues. These forces, along with the nanobody used as a G protein substitute, would be expected to influence the receptor conformation. Consistent with this, the intracellular end of the structure differs markedly from subsequent structures of active μOR bound to Gi protein.

CONCLUSIONS

Later Gi-bound structures are likely to be more accurate templates for ligand docking and modelling of active G protein-bound μOR. The possibility of reactions like this should be considered in the choice of protein truncation sites and purification conditions, and in the interpretation of excess or unexplained density.

Structural Insights into the Unexpected Agonism of Tetracyclic Antidepressants at Serotonin Receptors 5-HT1eR and 5-HT1FR

Serotonin (5-hydroxytryptamine, 5-HT) acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT1eR have confirmed roles in native tissue and are validated drug targets. Despite 5-HT1eR's therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT1eR's pharmacology in relation to the highly homologous 5-HT1FR, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT1e/1FR agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed anti-migraine properties. Using cryoEM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT1eR distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT1eR and 5-HT1FR contribute to the agonist activity of these antidepressants.

Synthesis, computational simulations and biological evaluation of new dual 5HT1A/5HT7 receptor ligands based on purine-2,6-dione scaffold

The new dual 5HT1A/5HT7 receptor ligands were designed based on the purine-2,6-dione scaffold with the fluorine atom. Twenty-one new derivatives were synthesized, and their structure-activity relationship was summarized. Compound 11 (7-(2-(3-fluorophenyl)-2-oxoethyl)-8-((4-(4-(2-methoxyphenyl)piperazin-1-yl)butyl)amino)-1,3-dimethyl-3,7-dihydro-1H-purine-2,6-dione) showed the highest affinity to 5HT1AR and 5HT7R, and was the most potent antagonist of 5-HT1AR (Kb = 0.26 ± 0.1 nM) which activity can be to reference compound NAN-190 (Kb = 0.26 ± 0.1 nM). The experimentally established physicochemical parameters of compound 11 showed that compound, as slightly ionized in the blood, could penetrate the blood-brain barrier. A molecular docking study showed that the fluorine substitution introduces additional stabilization effects on binding to 5HT1A/5HT7Rs. In animal assays of depression and anxiety, compound 11 revealed activity in terms of dosage compared to marketed psychotropics such as fluoxetine, citalopram, and sertraline.

The development of diphenyleneiodonium analogs as GPR3 agonists

G protein-coupled receptor 3 (GPR3) is an orphan receptor potentially involved in many important physiological processes such as drug abuse, neuropathic pain, and anxiety and depression related disorders. Pharmacological studies of GPR3 have been limited due to the restricted number of known agonists and inverse agonists for this constitutively active receptor. In this medicinal chemistry study, we report the discovery of GPR3 agonists based off the diphenyleneiodonium (DPI) scaffold. The most potent full agonist was the 3-trifluoromethoxy analog (32) with an EC50 of 260 nM and 90% efficacy compared to DPI. Investigation of a homology model of GPR3 from multiple sequence alignment resulted in the finding of a binding site rich in potential π-π and π-cation interactions stabilizing DPI-scaffold agonists. MMGBSA free energy analysis showed a good correlation with trends in observed EC50s. DPI analogs retained the same high receptor selectivity for GPR3 over GPR6 and GPR12 as observed with DPI. Collectively, the DPI analog series shows that order of magnitude improvements in potency with the scaffold were attainable; however, attempts to replace the iodonium ion to make the scaffold more druggable failed.

The Lateral Metalation of Isoxazolo[3,4-d]pyridazinones towards Hit-to-Lead Development of Selective Positive Modulators of Metabotropic Glutamate Receptors

Isoxazolo[3,4-d] pyridazinones ([3,4-d]s) were previously shown to have selective positive modulation at the metabotropic glutamate receptor (mGluR) Subtypes 2 and 4, with no functional cross-reactivity at mGluR1a, mGluR5, or mGluR8. Additional analogs were prepared to access more of the allosteric pocket and achieve higher binding affinity, as suggested by homology modeling. Two different sets of analogs were generated. One uses the fully formed [3,4-d] with an N6-aryl with and without halogens. These underwent successful selective lateral metalation and electrophilic quenching (LM&EQ) at the C3 of the isoxazole. In a second set of analogs, a phenyl group was introduced at the C4 position of the [3,4-d] ring via a condensation of 4-phenylacetyl-3-ethoxcarbonyl-5-methyl isoxazole with the corresponding hydrazine to generate the 3,4-ds 2b and 2j to 2n.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

Step-wise activation of a Family C GPCR

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain (ECD) that is linked via a cysteine-rich domain (CRDs) to their 7-transmembrane (TM) domain. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the ECD to the G protein-coupling TM. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

De Novo Design of κ-Opioid Receptor Antagonists Using a Generative Deep-Learning Framework

Likely effective pharmacological interventions for the treatment of opioid addiction include attempts to attenuate brain reward deficits during periods of abstinence. Pharmacological blockade of the κ-opioid receptor (KOR) has been shown to abolish brain reward deficits in rodents during withdrawal, as well as to reduce the escalation of opioid use in rats with extended access to opioids. Although KOR antagonists represent promising candidates for the treatment of opioid addiction, very few potent selective KOR antagonists are known to date and most of them exhibit significant safety concerns. Here, we used a generative deep-learning framework for the de novo design of chemotypes with putative KOR antagonistic activity. Molecules generated by models trained with this framework were prioritized for chemical synthesis based on their predicted optimal interactions with the receptor. Our models and proposed training protocol were experimentally validated by binding and functional assays.