Structures of the M1 and M2 muscarinic acetylcholine receptor/G-protein complexes

Research progress in myocardial function and diseases related to muscarinic acetylcholine receptor (Review)

Muscarinic acetylcholine (ACh) receptors (also known as M receptors) are widely distributed in all organs and tissues of the body, mainly playing a role in cholinergic nerve conduction. There are five known subtypes of muscarinic ACh receptors, but their pharmacological mechanisms of action on myocardial function have remained to be clearly defined. Functional myocardial diseases and myocardial injuries, such as arrhythmia, myocardial ischemia, myocarditis and myocardial fibrosis, may be affected by muscarinic ACh receptors. This article reviews the research progress of the regulation of myocardial function by muscarinic ACh receptors and related diseases, with the aim of developing better strategies and providing references for further revealing and clarifying the signal transduction and mechanisms of muscarinic ACh receptors in cardiomyocytes, and finding potential myocardial protective drugs that act on muscarinic ACh receptors.

The many facets of biased signaling: Mechanisms and possible therapeutic implications

Receptor-mediated cell activation frequently results in a plethora of effects, and interestingly, not all agonists that act on a given receptor activate all of those actions to the same extent. Biased agonism refers to this fact, i.e., the possibility to activate only a part of the receptor's signaling capabilities. It is worth mentioning that Biased Signaling is an integral concept that includes the system (organisms, isolated tissues, or cells), the individual receptor studied, and the ligands. It should be remembered that the system's genetic expression profile defines the type, abundance, and cellular localization of proteins that participate in signaling. This short review will be focused on G protein receptors, but biased signaling occurs in many other receptor types. Biased signaling can be related to the G proteins and β-arrestins available. Similarly, enzymes that catalyze receptor posttranslational modifications, such as phosphorylation, acylation, or ubiquitination, can play a role. G protein-coupled receptor signaling occurs at the plasma membrane, but it is well-established that endosomal signaling is a functional reality. Therefore, paying attention to cellular elements that participate in receptor endosomal traffic and destination (recycling to the plasma membrane/ degradation) is pertinent. There is still much to be known about these bias mechanisms, which are essential for basic knowledge of receptor drug action and for treating many pathological entities.

Structural insight into the self-activation and G-protein coupling of P2Y2 receptor

Purinergic P2Y2 receptor (P2Y2R) represents a typically extracellular ATP and UTP sensor for mediating purinergic signaling. Despite its importance as a pharmacological target, the molecular mechanisms underlying ligand recognition and G-protein coupling have remained elusive due to lack of structural information. In this study, we determined the cryo-electron microscopy (cryo-EM) structures of the apo P2Y2R in complex with Gq, ATP-bound P2Y2R in complex with Gq or Go, and UTP-bound P2Y4R in complex with Gq. These structures reveal the similarities and distinctions of ligand recognition within the P2Y receptor family. Furthermore, a comprehensive analysis of G-protein coupling reveals that P2Y2R exhibits promiscuity in coupling with both Gq and Go proteins. Combining molecular dynamics simulations and signaling assays, we elucidate the molecular mechanisms by which P2Y2R differentiates pathway-specific Gq or Go coupling through distinct structural components on the intracellular side. Strikingly, we identify a helix-like segment within the N-terminus that occupies the orthosteric ligand-binding pocket of P2Y2R, accounting for its self-activation. Taken together, these findings provide a molecular framework for understanding the activation mechanism of P2Y2R, encompassing ligand recognition, G-protein coupling, and a novel N-terminus-mediated self-activation mechanism.

Cyclic peptide inhibitors function as molecular glues to stabilize Gq/11 heterotrimers

Heterotrimeric Gα:Gβγ G proteins function as molecular switches downstream of G protein-coupled receptors (GPCRs). They alternate between a heterotrimeric GDP-bound OFF-state and a GTP-bound ON-state in which GαGTP is separated from the Gβγ dimer. Consequently, pharmacological tools to securely prevent the OFF-ON transition are of utmost importance to investigate their molecular switch function, specific contribution to GPCR signal transduction, and potential as drug targets. FR900359 (FR) and YM-254890 (YM), two natural cyclic peptides and highly specific inhibitors of Gq/11 heterotrimers, are exactly such tools. To date, their efficient and long-lasting inhibition of Gq/11 signaling has been attributed solely to a wedge-like binding to Gα, thereby preventing separation of the GTPase and α-helical domains and thus GDP release. Here, we use X-ray crystallography, biochemical and signaling assays, and BRET-based biosensors to show that FR and YM also function as stabilizers of the Gα:Gβγ subunit interface. Our high-resolution structures reveal a network of residues in Gα and two highly conserved amino acids in Gβ that are targeted by FR and YM to glue the Gβγ complex to the inactive GαGDP subunit. Unlike all previously developed nucleotide-state specific inhibitors that sequester Gα in its OFF-state but compete with Gβγ, FR and YM actively promote the inhibitory occlusion of GαGDP by Gβγ. In doing so, they securely lock the entire heterotrimer, not just Gα, in its inactive state. Our results identify FR and YM as molecular glues for Gα and Gβγ that combine simultaneous binding to both subunits with inhibition of G protein signaling.

Developing muscarinic receptor M1 classification models utilizing transfer learning and generative AI techniques

Muscarinic receptor subtype 1 (M1) is a G protein-coupled receptor (GPCR) and a key pharmacological target for peripheral neuropathy, chronic obstructive pulmonary disease, nerve agent exposures, and cognitive disorders. Screening and identifying compounds with potential to interact with M1 will aid in rational drug design for these disorders. In this work, we developed machine learning-based M1 classification models utilizing publicly available bioactivity data. As inactive compounds are rarely reported in the literature, we encountered the problem of imbalanced datasets. We investigated two strategies to overcome this bottleneck: 1) transfer learning and 2) using generative models to oversample the inactive class. Our analysis shows that these approaches reduced misclassification of the inactive class not only for M1 but also for other GPCR targets. Overall, we have developed classification models for M1 receptor that will enable rapid screening of large chemical databases and advance drug discovery.

Signal termination of the chemokine receptor CCR9 is governed by an arrestin-independent phosphorylation mechanism

The C-C chemokine receptor type 9 (CCR9) coordinates immune cell migration from the thymus to the small intestine along gradients of the chemokine CCL25. Receptor dysregulation is associated with a variety of inflammatory bowel diseases such as Crohn's and ulcerative colitis, whereas aberrant CCR9 overexpression correlates with tumor metastasis. Despite being an attractive therapeutic target, attempts to clinically antagonize CCR9 have been unsuccessful. This highlights the need for a deeper understanding of its specific regulatory mechanisms and signaling pathways. CCR9 is a G protein-coupled receptor (GPCR) and activates Gi and Gq pathways. Unexpectedly, live-cell bioluminescence resonance energy transfer assays reveal only limited G protein activation, and signaling is rapidly terminated. Truncating the receptor C terminus significantly enhanced G protein coupling, highlighting a regulatory role of this domain. Signal suppression was not because of canonical arrestin-coordinated desensitization. Rather, removal of GPCR kinase phosphorylation led to sustained and robust G protein activation by CCR9. Using site-directed mutagenesis, we identified specific phosphorylation motifs that attenuate G protein coupling. Receptor internalization did not correlate with G protein activation capabilities. Instead, CCR9 phosphorylation disrupted the interaction of G protein heterotrimers with the receptor. This interference may lead to rapid loss of productive coupling and downstream signaling as phosphorylation would effectively render the receptor incapable of G protein coupling. An arrestin-independent, phosphorylation-driven deactivation mechanism could complement arrestin-dependent regulation of other GPCRs and have consequences for therapeutically targeting these receptors.

Cryoelectron microscopy as a tool for illuminating activation mechanisms of human class A orphan G protein-coupled receptors

G protein-coupled receptors (GPCRs) are critically important medicinal targets, and the cryogenic electron microscopy (cryo-EM) revolution is providing novel high-resolution GPCR structures at a rapid pace. Orphan G protein-coupled receptors (oGPCRs) are a group of approximately 100 nonolfactory GPCRs for which endogenous ligands are unknown or not validated. The absence of modulating ligands adds difficulties to understanding the physiologic significance of oGPCRs and in the determination of high-resolution structures of isolated receptors that could facilitate drug discovery. Despite the challenges, cryo-EM structures of oGPCR-G protein complexes are emerging. This is being facilitated by numerous developments to stabilize GPCR-G protein complexes such as the use of dominant-negative G proteins, mini-G proteins, complex-stabilizing nanobodies or antibody fragments, and protein tethering methods. Moreover, many oGPCRs are constitutively active, which can facilitate complex formation in the absence of a known activating ligand. Consequently, in addition to providing templates for drug discovery, active oGPCR structures shed light on constitutive GPCR activation mechanisms. These comprise self-activation, whereby mobile extracellular portions of the receptor act as tethered agonists by occupying a canonical orthosteric-binding site in the transmembrane core, constitutive activity due to alterations to conserved molecular switches that stabilize inactive states of GPCRs, as well as receptors activated by cryptic ligands that are copurified with the receptor. Cryo-EM structures of oGPCRs are now being determined at a rapid pace and are expected to be invaluable tools for oGPCR drug discovery. SIGNIFICANCE STATEMENT: Orphan G protein-coupled receptors (GPCRs) provide large untapped potential for development of new medicines. Many of these receptors display constitutive activity, enabling structure determination and insights into observed GPCR constitutive activity including (1) self-activation by mobile receptor extracellular portions that function as tethered agonists, (2) modification of conserved motifs canonically involved in receptor quiescence and/or activation, and (3) activation by cryptic lipid ligands. Collectively, these studies advance fundamental understanding of GPCR function and provide opportunities for novel drug discovery.

Structural basis for the activation of proteinase-activated receptors PAR1 and PAR2

Members of the proteinase-activated receptor (PAR) subfamily of G protein-coupled receptors (GPCRs) play critical roles in processes like hemostasis, thrombosis, development, wound healing, inflammation, and cancer progression. Comprising PAR1-PAR4, these receptors are specifically activated by protease cleavage at their extracellular amino terminus, revealing a 'tethered ligand' that self-activates the receptor. This triggers complex intracellular signaling via G proteins and beta-arrestins, linking external protease signals to cellular functions. To date, direct structural visualization of these ligand-receptor complexes has been limited. Here, we present structural snapshots of activated PAR1 and PAR2 bound to their endogenous tethered ligands, revealing a shallow and constricted orthosteric binding pocket. Comparisons with antagonist-bound structures show minimal conformational changes in the TM6 helix and larger movements of TM7 upon activation. These findings reveal a common activation mechanism for PAR1 and PAR2, highlighting critical residues involved in ligand recognition. Additionally, the structure of PAR2 bound to a pathway selective antagonist, GB88, demonstrates how potent orthosteric engagement can be achieved by a small molecule mimicking the endogenous tethered ligand's interactions.

The Prosperity and Adversity of M4 Muscarinic Acetylcholine Receptor Activators in the Treatment of Neuropsychiatric Disorders

Since the serendipitous discovery of chlorpromazine in the 1950s, almost all current anti-schizophrenia drugs utilize the same mode of action by blocking the dopamine receptors in the brain. Unfortunately, these only treat part of the symptoms and are ineffective in almost 30% of patients. The recent FDA approval of Cobenfy, a coformulation of xanomeline, a M1/M4 muscarinic acetylcholine receptor (mAChR) agonist, and a peripherally restricted pan-mAChR blocker, has propelled the M4R as a validated and novel antipsychotic target. With >25 years of history in developing xanomeline, significant challenges remain in developing M4R activators, either at the ACh orthosteric binding site or allosterically via secondary less-conserved binding sites. Herein, we summarize recent successes and failures of M4R agonists and positive allosteric modulators, along with the progress in structure-activity relationship studies on both orthosteric and allosteric scaffolds to offer pathways for future therapeutics to this novel biological target for neuropsychiatric disorders.

5-HT2A receptors: Pharmacology and functional selectivity

Serotonin 5-HT2A receptors were one of the first serotonin receptors to be pharmacologically characterized. In mammals, they are expressed throughout the body in nearly every cell and tissue type, with the highest density in cortical layer V of the brain. They are involved in several aspects of normal physiological processes and behaviors and have been implicated in the etiology of neuropsychiatric diseases such as schizophrenia. Atypical antipsychotics have targeted blockade of 5-HT2A receptors as part of their therapeutic mechanism. More recently, 5-HT2A receptors have come to prominence for their role as the primary target for psychedelic drugs, which activate this receptor subtype to produce their characteristic behavioral effects. 5-HT2A receptor agonists like psilocybin, dimethyltryptamine, and lysergic acid diethylamide have each demonstrated long-lasting therapeutic efficacy in clinical trials for psychiatric disorders such as major depression and substance use disorders. There is a significant effort in both academia and industry to develop new agonists of 5-HT2A receptors with therapeutic efficacy. There are 3 primary scaffolds for agonists: tryptamines, ergolines, and phenylalkylamines, each engaging different subsets of amino acid residues in the receptor binding pocket. Differences can lead to differential responses between ligands for functionally selective outcomes. Here, we provide a historical perspective on 5-HT2A receptors, their key structural features and motifs involved in ligand-receptor interactions, and how these interactions can affect signaling pathways downstream of the receptor. Understanding how ligands interact with the 5-HT2A receptor will fundamentally inform future drug discovery to optimize therapeutics for a variety of disorders. SIGNIFICANCE STATEMENT: Psychedelic drugs have demonstrated long-lasting therapeutic efficacy for several conditions in multiple clinical trials. Their target, serotonin 5-HT2A receptors, are GPCRs with complex pharmacology. Having knowledge of how ligands interact with 5-HT2A receptors in the orthosteric binding pocket at the structural level to induce specific signal transduction pathways will inform on efforts to design and develop functionally selective drugs to potentially treat a variety of diseases.

Structural basis of the cysteinyl leukotriene receptor type 2 activation by LTD4

The G protein-coupled cysteinyl leukotriene receptor CysLT2R plays intricate roles in the physiology and pathogenesis of inflammation-related processes. It has garnered increasing attention as a potential therapeutic target for atopic asthma, brain injury, central nervous system disorders, and various types of cancer. In this study, we present the cryo-electron microscopy structure of the cysteinyl leukotriene D4 (LTD4)-bound human CysLT2R in complex with a Gαq protein, adopting an active conformation at a resolution of 3.15 Å. The structure elucidates a spacious polar pocket designed to accommodate the two branched negative ends of LTD4 and reveals a lateral ligand access route into the orthosteric pocket located on transmembrane domain helix (TM) 4 and 5. Furthermore, our findings highlight the crucial role of transmembrane domain helix 3 in sensing agonist moieties, representing the pivotal mechanism of receptor activation for both CysLT1R and CysLT2R. Collectively, the insights derived from our structural investigation establish a foundation for comprehending CysLT2R activation by its endogenous ligand LTD4, offering a rational basis for the design of drugs targeting CysLT2R.

Discovering Key Activation Hotspots in the M2 Muscarinic Receptor

The M2 muscarinic receptor (M2R) is a prototypical G protein-coupled receptor (GPCR) that serves as a model system for understanding ligand recognition and GPCR activation. Here, using vibrational spectroscopy, we identify the mechanisms governing M2R activation by its native agonist, acetylcholine. Combined with mutagenesis, computational chemistry, and organic synthetic chemistry, our analyses found that the precise distance between acetylcholine and Asn404, one of the amino acids constituting the ligand-binding site, is important for M2R activation and that the N404Q mutant undergoes partial active state-like conformational changes. We discovered that a water molecule bridging acetylcholine and Asn404 forms a precise and flexible hydrogen bond network, triggering the outward movement of transmembrane helix 6 in M2R. Consistent with this observation, disruptions in this hydrogen bond network via chemical modification at the α- or β-position of acetylcholine failed to activate M2R. Collectively, our findings pinpoint Asn404 as a critical residue that both senses acetylcholine binding and induces M2R activation.

Influence of expression and purification protocols on Gα biochemical activity: kinetics of plant and mammalian G protein cycles

Heterotrimeric G proteins, composed of Gα, Gβ, and Gγ subunits, are a class of signal transduction complexes with broad roles in human health and agriculturally relevant plant physiological and developmental traits. In the classic paradigm, guanine nucleotide binding to the Gα subunit regulates the activation status of the complex. We sought to develop improved methods for heterologous expression and rapid purification of Gα subunits, initially targeting GPA1, the sole canonical Gα subunit of the model plant species, Arabidopsis thaliana. Compared to conventional methods, our expression methodology and rapid StrepII-tag mediated purification facilitates substantially higher yield, and isolation of protein with increased GTP binding and hydrolysis activities. Human GNAI1 purified using our approach displayed the expected binding and hydrolysis activities, indicating our protocol is applicable to mammalian Gα subunits, potentially including those for which purification of enzymatically active protein has been historically problematic. We subsequently utilized domain swaps of GPA1 and human GNAO1 to demonstrate that the inherent instability of GPA1 is a function of the interaction between the Ras and helical domains. Additionally, we found that GPA1-GNAO1 domain swaps partially uncouple the instability from the rapid nucleotide binding kinetics displayed by GPA1. In summary, our work provides insights into methods to optimally study heterotrimeric G proteins, and reveals roles of the helical domain in Gα kinetics and stability.

Functional dynamics of G protein-coupled receptors reveal new routes for drug discovery

G protein-coupled receptors (GPCRs) are the largest human membrane protein family that transduce extracellular signals into cellular responses. They are major pharmacological targets, with approximately 26% of marketed drugs targeting GPCRs, primarily at their orthosteric binding site. Despite their prominence, predicting the pharmacological effects of novel GPCR-targeting drugs remains challenging due to the complex functional dynamics of these receptors. Recent advances in X-ray crystallography, cryo-electron microscopy, spectroscopic techniques and molecular simulations have enhanced our understanding of receptor conformational dynamics and ligand interactions with GPCRs. These developments have revealed novel ligand-binding modes, mechanisms of action and druggable pockets. In this Review, we highlight such aspects for recently discovered small-molecule drugs and drug candidates targeting GPCRs, focusing on three categories: allosteric modulators, biased ligands, and bivalent and bitopic compounds. Although studies so far have largely been retrospective, integrating structural data on ligand-induced receptor functional dynamics into the drug discovery pipeline has the potential to guide the identification of drug candidates with specific abilities to modulate GPCR interactions with intracellular effector proteins such as G proteins and β-arrestins, enabling more tailored selectivity and efficacy profiles.

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein-coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

Structural basis for the constitutive activity of the melanocortin receptor family

The constitutive activity of melanocortin receptors (MCRs) is integral to several physiological processes. The unliganded cryo-electron microscopy structures of MC1R, MC2R, MC3R, MC4R, and MC5R in complex with Gs proteins determined at global resolutions of 2.98 Å, 3.01 Å, 2.75 Å, 3.12 Å, and 2.86 Å, respectively, revealed that their binding poses and interactions with Gs are similar to those of agonist-bound MCRs. The extracellular regions of the transmembrane helices (TMs) exhibit distinct conformational rearrangements, characterized by varying shifts of TM3 and outward displacements of TM4. These variations represent unique structural features of constitutively activated MCRs. Unassigned electron densities were observed within the orthosteric pockets where extensive interactions with cognate ligands occur. In addition, zinc ions, but not calcium, positively regulated MC4R activity in a dose-dependent manner. Our findings provide valuable insights into the molecular mechanisms underlying MCR basal activity and highlight the role of divalent ions in receptor activation.

Visualizing agonist-induced M2 receptor activation regulated by aromatic ring dynamics

Despite the growing number of G protein-coupled receptor (GPCR) structures being resolved, the dynamic process of how GPCRs transit from the inactive toward the active state remains unclear. In this study, comprehensive molecular dynamics simulations were performed to explore how ligand binding modulates the conformational dynamics of the M2 muscarinic acetylcholine receptor (M2R). We observed a sequential occurrence of structural changes in the inactive-to-active transition of M2R induced by a superagonist iperoxo, which includes the orthosteric binding site contraction, the TM6 opening into an intermediate conformation, and a further structural change toward full activation upon binding to G protein or a G protein mimetic nanobody. Two activation intermediates were identified, which show structural features different from those reported for apo-GPCRs. Moreover, our results suggest that stabilization of a specific W4006.48 conformation and enhanced F3966.44 dynamics are crucial for activation, whereas distinct side-chain rotamer equilibriums of Y2065.58 in the cytoplasmic cavity are correlated with agonist efficacies. Our work provides atomic-level structural insights into the agonist-induced M2R activation pathway and highlights a mechanism by which ligand efficacy can be encoded and transduced in the form of aromatic ring dynamics.

How Ligands Achieve Biased Signaling toward Arrestins

G protein-coupled receptors (GPCRs) mediate the effects of various endogenous and extracellular stimuli through multiple transducers, including heterotrimeric G proteins, GPCR kinases (GRKs), and arrestins. Biased signaling, which preferentially activates certain G protein or GRK/arrestin signaling pathways, provides great opportunities for developing drugs with enhanced therapeutic efficacy and minimized side effects. In this Review, we review studies addressing the structural dynamics of GPCRs bound to balanced and biased ligands and current consensus on how ligand-receptor interactions determine signaling outcomes. We also examine the conformational changes in GPCRs when in complex with G proteins, arrestins, and GRKs, highlighting a more profound impact of signal transducers on receptor rearrangements compared with biased ligands. This evidence supports the idea that biased signaling can be achieved through the promotion of multiple conformational states by biased agonists and the stabilization of specific active conformations by individual signal transducers.

Structural insights into the engagement of lysophosphatidic acid receptor 1 with different G proteins

Lysophosphatidic acid (LPA) and sphingosine-1-phosphate (S1P) are bioactive lysophospholipids derived from cell membranes that activate the endothelial differentiation gene family of G protein-coupled receptors. Activation of these receptors triggers multiple downstream signaling cascades through G proteins such as Gi/o, Gq/11, and G12/13. Therefore, LPA and S1P mediate several physiological processes, including cytoskeletal dynamics, neurite retraction, cell migration, cell proliferation, and intracellular ion fluxes. The basis for the G-protein coupling selectivity of EDG receptors, however, remains unknown. Here, we present cryo-electron microscopy structures of LPA-activated LPA1 in complexes with Gi, Gq, and G13 heterotrimers. Comparison of the three LPA1-G protein structures shows clearly different conformations of intracellular loop 2 (ICL2) and ICL3 that are likely induced by the different Gα protein interfaces. Interestingly, this G-protein interface interaction is a common feature of LPA and S1P receptors. Our findings provide clues to understanding the promiscuity of G-protein coupling in the endothelial differentiation gene family.

The structure and function of the ghrelin receptor coding for drug actions

Drugs targeting the ghrelin receptor hold therapeutic potential in anorexia, obesity and diabetes. However, developing effective drugs is challenging. To tackle this common issue across a broad drug target, this study aims to understand how anamorelin, the only approved drug targeting the ghrelin receptor, operates compared to other synthetic drugs. Our research elucidated the receptor's structure with anamorelin and miniGq, unveiling anamorelin's superagonistic activity. We demonstrated that ligands with distinct chemical structures uniquely bind to the receptor, resulting in diverse conformations and biasing signal transduction. Moreover, our study showcased the utility of structural information in effectively identifying natural genetic variations altering drug action and causing severe functional deficiencies, offering a basis for selecting the right medication on the basis of the individual's genomic sequence. Thus, by building on structural analysis, this study enhances the foundational framework for selecting therapeutic agents targeting the ghrelin receptor, by effectively leveraging signaling bias and genetic variations.

G protein-coupled receptors: A golden key to the treasure-trove of neurodegenerative diseases

G protein-coupled receptors (GPCRs) are a class of transmembrane proteins that distribute in various organs extensively. They can regulate physiological functions such as perception, neurotransmission and endocrinology through the synergies of signaling pathways. At present, Food and Drug Administration (FDA) have approved more than 500 drugs targeting GPCRs to treat a variety of conditions, including neurological diseases, gastrointestinal diseases and tumors. Conformational diversity and dynamic changes make GPCRs a star target for the treatment of neurodegenerative diseases. Moreover, GPCRs can also open biased signaling pathways for G protein and β-arrestin, which has unique functional selectivity and the possibility of overcoming side effects. Some studies believe that biased drugs will be the mainstream direction of drug innovation in the future. To disclose the essential role and research process of GPCRs in neurodegenerative diseases, we firstly reviewed several pivotal GPCRs and their mediated signaling pathways in Alzheimer's disease (AD), Parkinson's disease (PD) and Amyotrophic lateral sclerosis (ALS). Then we focused on the biased signaling pathway of GPCRs in these diseases. Finally, we updated the GPCR drugs under research for the treatment of neurodegenerative diseases in the clinical trials or approval. This review could provide valuable targets for precision therapy to cope with the dysfunction of neurodegenerative diseases in the future.

Cryo-EM Structures and AlphaFold3 Models of Histamine Receptors Reveal Diverse Ligand Binding and G Protein Bias

Background: The four subtypes of G protein-coupled receptors (GPCRs) regulated by histamine play critical roles in various physiological and pathological processes, such as allergy, gastric acid secretion, cognitive and sleep disorders, and inflammation. Previous experimental structures of histamine receptors (HRs) with agonists and antagonists exhibited multiple conformations for the ligands and G protein binding. However, the structural basis for HR regulation and signaling remains elusive. Methods: We determined the cryo-electron microscopy (cryo-EM) structure of the H4R-histamine-Gi complex at 2.9 Å resolution, and predicted the models for all four HRs in the ligand-free apo and G protein subtype binding states using AlphaFold3 (AF3). Results: By comparing our H4R structure with the experimental HR structures and the computational AF3 models, we elucidated the distinct histamine binding modes and G protein interfaces, and proposed the essential roles of Y6.51 and Q7.42 in receptor activation and the intracellular loop 2 (ICL2) in G protein bias. Conclusions: Our findings deciphered the molecular mechanisms underlying the regulation of different HRs, from the extracellular ligand-binding pockets and transmembrane motifs to the intracellular G protein coupling interfaces. These insights are expected to facilitate selective drug discovery targeting HRs for diverse therapeutic purposes.

BaNDyT: Bayesian Network Modeling of Molecular Dynamics Trajectories

Bayesian network modeling (BN modeling, or BNM) is an interpretable machine learning method for constructing probabilistic graphical models from the data. In recent years, it has been extensively applied to diverse types of biomedical data sets. Concurrently, our ability to perform long-time scale molecular dynamics (MD) simulations on proteins and other materials has increased exponentially. However, the analysis of MD simulation trajectories has not been data-driven but rather dependent on the user's prior knowledge of the systems, thus limiting the scope and utility of the MD simulations. Recently, we pioneered using BNM for analyzing the MD trajectories of protein complexes. The resulting BN models yield novel fully data-driven insights into the functional importance of the amino acid residues that modulate proteins' function. In this report, we describe the BaNDyT software package that implements the BNM specifically attuned to the MD simulation trajectories data. We believe that BaNDyT is the first software package to include specialized and advanced features for analyzing MD simulation trajectories using a probabilistic graphical network model. We describe here the software's uses, the methods associated with it, and a comprehensive Python interface to the underlying generalist BNM code. This provides a powerful and versatile mechanism for users to control the workflow. As an application example, we have utilized this methodology and associated software to study how membrane proteins, specifically the G protein-coupled receptors, selectively couple to G proteins. The software can be used for analyzing MD trajectories of any protein as well as polymeric materials.

Cryo-EM reveals a new allosteric binding site at the M5 mAChR

The M5 muscarinic acetylcholine receptor (M5 mAChR) represents a promising therapeutic target for neurological disorders. However, the high conservation of its orthosteric binding site has posed significant challenges for drug development. While selective positive allosteric modulators (PAMs) offer a potential solution, a structural understanding of the M5 mAChR and its allosteric binding sites has remained limited. Here, we present a 2.8 Å cryo-electron microscopy structure of the M5 mAChR complexed with heterotrimeric Gq protein and the agonist iperoxo, completing the active-state structural characterization of the mAChR family. To identify the binding site of M5-selective PAMs, we implemented an integrated approach combining mutagenesis, pharmacological assays, structural biology, and molecular dynamics simulations. Our mutagenesis studies revealed that selective M5 PAMs bind outside previously characterized M5 mAChR allosteric sites. Subsequently, we obtained a 2.1 Å structure of M5 mAChR co-bound with acetylcholine and the selective PAM VU6007678, revealing a novel allosteric pocket at the extrahelical interface between transmembrane domains 3 and 4 that was confirmed through mutagenesis and simulations. These findings demonstrate the diverse mechanisms of allosteric regulation in mAChRs and highlight the value of integrating pharmacological and structural approaches to identify novel allosteric binding sites.

G protein peptidomimetics reveal allosteric effects and stepwise interactions in ghrelin receptor-G protein coupling

G protein-coupled receptor (GPCR) signaling is a dynamic process involving various conformational intermediates in addition to those captured in static three-dimensional structures. Here, we used newly developed G protein peptidomimetics to characterize the interactions of the ghrelin receptor (GHSR) with G proteins. Coupling to the G protein peptidomimetic not only affected the conformational features of the cytoplasmic regions of the receptor where the G protein binds but also allosterically affected the extracellular ligand-binding pocket. These conformational and allosteric changes increased the affinity of G protein-coupled GHSR for the endogenous agonist ghrelin. In addition, our data identified different complexes along the G protein activation pathway that differed in the engagement of the Gαq C-terminal helix. Given that this helix is the main link between the activated receptor and the Gα nucleotide-binding pocket, these findings suggested a stepwise process involving distinct states in GPCR-catalyzed G protein activation. Collectively, our results provide evidence for the dynamic behavior of GPCR-G protein signaling complexes, with such dynamics most likely contributing to signaling selectivity and/or efficacy.

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a new Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

Balancing G protein selectivity and efficacy in the adenosine A2A receptor

The adenosine A2A receptor (A2AR) engages several G proteins, notably Go and its cognate Gs protein. This coupling promiscuity is facilitated by a dynamic ensemble, revealed by 19F nuclear magnetic resonance imaging of A2AR and G protein. Two transmembrane helix 6 (TM6) activation states, formerly associated with partial and full agonism, accommodate the differing volumes of Gs and Go. While nucleotide depletion biases TM7 toward a fully active state in A2AR-Gs, A2AR-Go is characterized by a dynamic inactive/intermediate fraction. Molecular dynamics simulations reveal that the NPxxY motif, a highly conserved switch, establishes a unique configuration in the A2AR-Go complex, failing to stabilize the helix-8 interface with Gs, and adoption of the active state. The resulting TM7 dynamics hamper G protein coupling, suggesting kinetic gating may be responsible for reduced efficacy in the noncognate G protein complex. Thus, dual TM6 activation states enable greater diversity of coupling partners while TM7 dynamics dictate coupling efficacy.

Flavors of GPCR signaling bias

GPCRs are inherently flexible molecules existing in an equilibrium of multiple conformations. Binding of GPCR agonists shifts this equilibrium. Certain agonists can increase the fraction of active-like conformations that predispose the receptor to coupling to a particular signal transducer or a select group of transducers. Such agonists are called biased, in contrast to balanced agonists that facilitate signaling via all transducers the receptor couples to. These biased agonists preferentially channel the signaling of a GPCR to particular G proteins, GRKs, or arrestins. Preferential activation of particular G protein or arrestin subtypes can be beneficial, as it would reduce unwanted on-target side effects, widening the therapeutic window. However, biasing GPCRs has two important limitations: a) complete bias is impossible due to inherent flexibility of GPCRs; b) receptor-independent functions of signal transducer proteins cannot be directly affected by GPCR ligands or differential receptor barcoding by GRK phosphorylation. This article is part of the Special Issue on "Ligand Bias".

A non-canonical mechanism of GPCR activation

The goal of designing safer, more effective drugs has led to tremendous interest in molecular mechanisms through which ligands can precisely manipulate the signaling of G-protein-coupled receptors (GPCRs), the largest class of drug targets. Decades of research have led to the widely accepted view that all agonists-ligands that trigger GPCR activation-function by causing rearrangement of the GPCR's transmembrane helices, opening an intracellular pocket for binding of transducer proteins. Here we demonstrate that certain agonists instead trigger activation of free fatty acid receptor 1 by directly rearranging an intracellular loop that interacts with transducers. We validate the predictions of our atomic-level simulations by targeted mutagenesis; specific mutations that disrupt interactions with the intracellular loop convert these agonists into inverse agonists. Further analysis suggests that allosteric ligands could regulate the signaling of many other GPCRs via a similar mechanism, offering rich possibilities for precise control of pharmaceutically important targets.

BaNDyT: Bayesian Network modeling of molecular Dynamics Trajectories

Bayesian network modeling (BN modeling, or BNM) is an interpretable machine learning method for constructing probabilistic graphical models from the data. In recent years, it has been extensively applied to diverse types of biomedical datasets. Concurrently, our ability to perform long-timescale molecular dynamics (MD) simulations on proteins and other materials has increased exponentially. However, the analysis of MD simulation trajectories has not been data-driven but rather dependent on the user's prior knowledge of the systems, thus limiting the scope and utility of the MD simulations. Recently, we pioneered using BNM for analyzing the MD trajectories of protein complexes. The resulting BN models yield novel fully data-driven insights into the functional importance of the amino acid residues that modulate proteins' function. In this report, we describe the BaNDyT software package that implements the BNM specifically attuned to the MD simulation trajectories data. We believe that BaNDyT is the first software package to include specialized and advanced features for analyzing MD simulation trajectories using a probabilistic graphical network model. We describe here the software's uses, the methods associated with it, and a comprehensive Python interface to the underlying generalist BNM code. This provides a powerful and versatile mechanism for users to control the workflow. As an application example, we have utilized this methodology and associated software to study how membrane proteins, specifically the G protein-coupled receptors, selectively couple to G proteins. The software can be used for analyzing MD trajectories of any protein as well as polymeric materials.

Ligand bias at the muscarinic acetylcholine receptor family: Opportunities and challenges

Muscarinic acetylcholine receptors (mAChRs) are G protein-coupled receptors (GPCRs) that are activated by the endogenous neurotransmitter, acetylcholine (ACh). Disruption of mAChR signalling has been associated with a variety of neurological disorders and non-neurological diseases. Consequently, the development of agonists and antagonists of the mAChRs has been a major avenue in drug discovery. Unfortunately, mAChR ligands are often associated with on-target side effects for two reasons. The first reason is due to the high sequence conservation at the orthosteric ACh binding site among all five receptor subtypes (M1-M5), making on-target subtype selectivity a major challenge. The second reason is due to on-target side effects of mAChR drugs that are associated with the pleiotropic nature of mAChR signalling at the level of a single mAChR subtype. Indeed, there is growing evidence that within the myriad of signalling events produced by mAChR ligands, some will have therapeutic benefits, whilst others may promote cholinergic side effects. This paradigm of drug action, known as ligand bias or biased agonism, is an attractive feature for next-generation mAChR drugs, as it holds the promise of developing drugs devoid of on-target adverse effects. Although relatively simple to detect and even quantify in vitro, ligand bias, as observed in recombinant systems, does not always translate to in vivo systems, which remains a major hurdle in GPCR drug discovery, including the mAChR family. Here we report recent studies that have attempted to detect and quantify ligand bias at the mAChR family, and briefly discuss the challenges associated with biased agonist drug development. This article is part of the Special Issue on "Ligand Bias".

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

Understanding the kinetics of LSD in receptors and subsequent induced signaling is crucial for comprehending both the psychoactive and therapeutic effects of LSD. Despite extensive research on LSD's interactions with serotonin 2A and 2B receptors, its behavior on other targets, including dopamine receptors, has remained elusive. Here, we present cryo-EM structures of LSD/PF6142-bound dopamine D1 receptor (DRD1)-legobody complexes, accompanied by a β-arrestin-mimicking nanobody, NBA3, shedding light on the determinants of G protein coupling versus β-arrestin coupling. Structural analysis unveils a distinctive binding mode of LSD in DRD1, particularly with the ergoline moiety oriented toward TM4. Kinetic investigations uncover an exceptionally rapid dissociation rate of LSD in DRD1, attributed to the flexibility of extracellular loop 2 (ECL2). Moreover, G protein can stabilize ECL2 conformation, leading to a significant slowdown in ligand's dissociation rate. These findings establish a solid foundation for further exploration of G protein-coupled receptor (GPCR) dynamics and their relevance to signal transduction.

Molecular mechanism of the endothelin receptor type B interactions with Gs, Gi, and Gq

The endothelin receptor type B (ETB) exhibits promiscuous coupling with various heterotrimeric G protein subtypes including Gs, Gi/o, Gq/11, and G12/13. Recent fluorescence and structural studies have raised questions regarding the coupling efficiencies and determinants of these G protein subtypes. Herein, by utilizing an integrative approach, combining hydrogen/deuterium exchange mass spectrometry and NanoLuc Binary Technology-based cellular systems, we investigated conformational changes of Gs, Gi, and Gq triggered by ETB activation. ETB coupled to Gi and Gq but not with Gs. We underscored the critical roles of specific regions, including the C terminus of Gα and intracellular loop 2 (ICL2) of ETB in ETB-Gi1 or ETB-Gq coupling. Although The C terminus of Gα is essential for ETB-Gi1 and ETB-Gq coupling, ETB ICL2 influences Gq-coupling but not Gi1-coupling. Our results suggest a differential coupling efficiency of ETB with Gs, Gi1, and Gq, accompanied by distinct conformational changes in G proteins upon ETB-induced activation.

Clinical and Molecular Profiling in GNAO1 Permits Phenotype-Genotype Correlation

BACKGROUND

Defects in GNAO1, the gene encoding the major neuronal G-protein Gαo, are related to neurodevelopmental disorders, epilepsy, and movement disorders. Nevertheless, there is a poor understanding of how molecular mechanisms explain the different phenotypes.

OBJECTIVES

We aimed to analyze the clinical phenotype and the molecular characterization of GNAO1-related disorders.

METHODS

Patients were recruited in collaboration with the Spanish GNAO1 Association. For patient phenotyping, direct clinical evaluation, analysis of homemade-videos, and an online questionnaire completed by families were analyzed. We studied Gαo cellular expression, the interactions of the partner proteins, and binding to guanosine triphosphate (GTP) and G-protein-coupled receptors (GPCRs).

RESULTS

Eighteen patients with GNAO1 genetic defects had a complex neurodevelopmental disorder, epilepsy, central hypotonia, and movement disorders. Eleven patients showed neurological deterioration, recurrent hyperkinetic crisis with partial recovery, and secondary complications leading to death in three cases. Deep brain stimulation improved hyperkinetic crisis, but had inconsistent benefits in dystonia. The molecular defects caused by pathogenic Gαo were aberrant GTP binding and hydrolysis activities, an inability to interact with cellular binding partners, and reduced coupling to GPCRs. Decreased localization of Gαo in the plasma membrane was correlated with the phenotype of "developmental and epileptic encephalopathy 17." We observed a genotype-phenotype correlation, pathogenic variants in position 203 were related to developmental and epileptic encephalopathy, whereas those in position 209 were related to neurodevelopmental disorder with involuntary movements. Milder phenotypes were associated with other molecular defects such as del.16q12.2q21 and I344del.

CONCLUSION

We highlight the complexity of the motor phenotype, which is characterized by fluctuations throughout the day, and hyperkinetic crisis with a distinct post-hyperkinetic crisis state. We confirm a molecular-based genotype-phenotype correlation for specific variants. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

An overview of sphingosine-1-phosphate receptor 2: Structure, biological function, and small-molecule modulators

Over the past decade, there has been a notable increase in research on sphingosine-1-phosphate receptor 2 (S1PR2), which is a type of G-protein-coupled receptor. Upon activation by S1P or other ligands, S1PR2 initiates downstream signaling pathways such as phosphoinositide 3-kinase (PI3K), Mitogen-activated protein kinase (MAPK), Rho/Rho-associated coiled-coil containing kinases (ROCK), and others, contributing to the diverse biological functions of S1PR2 and playing a pivotal role in various physiological processes and disease progressions, such as multiple sclerosis, fibrosis, inflammation, and tumors. Due to the extensive biological functions of S1PR2, many S1PR2 modulators, including agonists and antagonists, have been developed and discovered by pharmaceutical companies (e.g., Novartis and Galapagos NV) and academic medicinal chemists for disease diagnosis and treatment. However, few reviews have been published that comprehensively overview the functions and regulators of S1PR2. Herein, we provide an in-depth review of the advances in the function of S1PR2 and its modulators. We first summarize the structure and biological function of S1PR2 and its pathological role in human diseases. We then focus on the discovery approach, design strategy, development process, and biomedical application of S1PR2 modulators. Additionally, we outline the major challenges and future directions in this field. Our comprehensive review will aid in the discovery and development of more effective and clinically applicable S1PR2 modulators.

Structural insights into ligand recognition, selectivity, and activation of bombesin receptor subtype-3

Bombesin receptor subtype-3 (BRS3) is an important orphan G protein-coupled receptor that regulates energy homeostasis and insulin secretion. As a member of the bombesin receptor (BnR) family, the lack of known endogenous ligands and high-resolution structure has hindered the understanding of BRS3 signaling and function. We present two cryogenic electron microscopy (cryo-EM) structures of BRS3 in complex with the heterotrimeric Gq protein in its active states: one bound to the pan-BnR agonist BA1 and the other bound to the synthetic BRS3-specific agonist MK-5046. These structures reveal the architecture of the orthosteric ligand pocket underpinning molecular recognition and provide insights into the structural basis for BRS3's selectivity and low affinity for bombesin peptides. Examination of conserved micro-switches suggests a shared activation mechanism among BnRs. Our findings shed light on BRS3's ligand selectivity and signaling mechanisms, paving the way for exploring its therapeutic potential for diabetes, obesity, and related metabolic disorders.

The plasma membrane inner leaflet PI(4,5)P2 is essential for the activation of proton-activated chloride channels

Proton-activated chloride (PAC) channels, ubiquitously expressed in tissues, regulate intracellular Cl- levels and cell death following acidosis. However, molecular mechanisms and signaling pathways involved in PAC channel modulation are largely unknown. Herein, we determine that phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2] of the plasma membrane inner leaflet is essential for the proton activation of PAC channels. PI(4,5)P2 depletion by activating phosphatidylinositol 5-phosphatases or Gq protein-coupled muscarinic receptors substantially inhibits human PAC currents. In excised inside-out patches, PI(4,5)P2 application to the cytoplasmic side increases the currents. Structural simulation reveals that the putative PI(4,5)P2-binding site is localized within the cytosol in resting state but shifts to the cell membrane's inner surface in an activated state and interacts with inner leaflet PI(4,5)P2. Alanine neutralization of basic residues near the membrane-cytosol interface of the transmembrane helice 2 significantly attenuates PAC currents. Overall, our study uncovers a modulatory mechanism of PAC channel through inner membrane PI(4,5)P2.

Structural basis for hormone recognition and distinctive Gq protein coupling by the kisspeptin receptor

Kisspeptin signaling through its G protein-coupled receptor, KISS1R, plays an indispensable role in regulating reproduction via the hypothalamic-pituitary-gonadal axis. Dysregulation of this pathway underlies severe disorders like infertility and precocious puberty. Here, we present cryo-EM structures of KISS1R bound to the endogenous agonist kisspeptin-10 and a synthetic analog TAK-448. These structures reveal pivotal interactions between peptide ligands and KISS1R extracellular loops for receptor activation. Both peptides exhibit a conserved binding mode, unveiling their common activation mechanism. Intriguingly, KISS1R displays a distinct 40° angular deviation in its intracellular TM6 region compared to other Gq-coupled receptors, enabling distinct interactions with Gq. This study reveals the molecular intricacies governing ligand binding and activation of KISS1R, while highlighting its exceptional ability to couple with Gq. Our findings pave the way for structure-guided design of therapeutics targeting this physiologically indispensable receptor.

GOLEM: Automated and Robust Cryo-EM-Guided Ligand Docking with Explicit Water Molecules

A detailed understanding of ligand-protein interaction is essential for developing rational drug-design strategies. In recent years, technological advances in cryo-electron microscopy (cryo-EM) brought a new era to the structural determination of biological macromolecules and assemblies at high resolution, marking cryo-EM as a promising tool for studying ligand-protein interactions. However, even in high-resolution cryo-EM results, the densities for the bound small-molecule ligands are often of lower quality due to their relatively dynamic and flexible nature, frustrating their accurate coordinate assignment. To address the challenge of ligand modeling in cryo-EM maps, here we report the development of GOLEM (Genetic Optimization of Ligands in Experimental Maps), an automated and robust ligand docking method that predicts a ligand's pose and conformation in cryo-EM maps. GOLEM employs a Lamarckian genetic algorithm to perform a hybrid global/local search for exploring the ligand's conformational, orientational, and positional space. As an important feature, GOLEM explicitly considers water molecules and places them at optimal positions and orientations. GOLEM takes into account both molecular energetics and the correlation with the cryo-EM maps in its scoring function to optimally place the ligand. We have validated GOLEM against multiple cryo-EM structures with a wide range of map resolutions and ligand types, returning ligand poses in excellent agreement with the densities. As a VMD plugin, GOLEM is free of charge and accessible to the community. With these features, GOLEM will provide a valuable tool for ligand modeling in cryo-EM efforts toward drug discovery.

Control of G protein-coupled receptor function via membrane-interacting intrinsically disordered C-terminal domains

G protein-coupled receptors (GPCRs) control intracellular signaling cascades via agonist-dependent coupling to intracellular transducers including heterotrimeric G proteins, GPCR kinases (GRKs), and arrestins. In addition to their critical interactions with the transmembrane core of active GPCRs, all three classes of transducers have also been reported to interact with receptor C-terminal domains (CTDs). An underexplored aspect of GPCR CTDs is their possible role as lipid sensors given their proximity to the membrane. CTD-membrane interactions have the potential to control the accessibility of key regulatory CTD residues to downstream effectors and transducers. Here, we report that the CTDs of two closely related family C GPCRs, metabotropic glutamate receptor 2 (mGluR2) and mGluR3, bind to membranes and that this interaction can regulate receptor function. We first characterize CTD structure with NMR spectroscopy, revealing lipid composition-dependent modes of membrane binding. Using molecular dynamics simulations and structure-guided mutagenesis, we then identify key conserved residues and cancer-associated mutations that modulate CTD-membrane binding. Finally, we provide evidence that mGluR3 transducer coupling is controlled by CTD-membrane interactions in live cells, which may be subject to regulation by CTD phosphorylation and changes in membrane composition. This work reveals an additional mechanism of GPCR modulation, suggesting that CTD-membrane binding may be a general regulatory mode throughout the broad GPCR superfamily.

CASTpFold: Computed Atlas of Surface Topography of the universe of protein Folds

Geometric and topological properties of protein structures, including surface pockets, interior cavities and cross channels, are of fundamental importance for proteins to carry out their functions. Computed Atlas of Surface Topography of proteins (CASTp) is a widely used web server for locating, delineating, and measuring these geometric and topological properties of protein structures. Recent developments in AI-based protein structure prediction such as AlphaFold2 (AF2) have significantly expanded our knowledge on protein structures. Here we present CASTpFold, a continuation of CASTp that provides accurate and comprehensive identifications and quantifications of protein topography. It now provides (i) results on an expanded database of proteins, including the Protein Data Bank (PDB) and non-singleton representative structures of AlphaFold2 structures, covering 183 million AF2 structures; (ii) functional pockets prediction with corresponding Gene Ontology (GO) terms or Enzyme Commission (EC) numbers for AF2-predicted structures and (iii) pocket similarity search function for surface and protein-protein interface pockets. The CASTpFold web server is freely accessible at https://cfold.bme.uic.edu/castpfold/.

Cryo-EM advances in GPCR structure determination

G-protein-coupled receptors (GPCRs) constitute a prominent superfamily in humans and are categorized into six classes (A-F) that play indispensable roles in cellular communication and therapeutics. Nonetheless, their structural comprehension has been limited by challenges in high-resolution data acquisition. This review highlights the transformative impact of cryogenic electron microscopy (cryo-EM) on the structural determinations of GPCR-G-protein complexes. Specific technologies, such as nanobodies and mini-G-proteins, stabilize complexes and facilitate structural determination. We discuss the structural alterations upon receptor activation in different GPCR classes, revealing their diverse mechanisms. This review highlights the robust foundation for comprehending GPCR function and pave the way for future breakthroughs in drug discovery and therapeutic targeting.

Allosteric modulators of M1 muscarinic receptors enhance acetylcholine efficacy and decrease locomotor activity and turning behaviors in zebrafish

Allosteric modulation of muscarinic acetylcholine receptors (mAChR) has been identified as a potential strategy for regulating cholinergic signaling in the treatment of various neurological disorders. Most positive allosteric modulators (PAMs) of mAChR enhance agonist affinity and potency, while very few PAMs (e.g., amiodarone) selectively enhance G protein coupling efficacy. The key structural features of amiodarone responsible for enhancement of mAChR efficacy were examined in CHO cells expressing M1 receptors. Subsequent incorporation of these structural features into previously identified allosteric modulators of potency (i.e., n-benzyl isatins) generated ligands that demonstrated similar or better enhancement of mAChR efficacy, lower in vivo toxicity, and higher allosteric binding affinity relative to amiodarone. Notable ligands include 8a, c which respectively demonstrated the strongest binding affinity and the most robust enhancement of mAChR efficacy as calculated from an allosteric operational model. Amiodarone derivatives and hybrid ligands were additionally screened in wildtype zebrafish (Danio rerio) to provide preliminary in vivo toxicity data as well as to observe effects on locomotor and turning behaviors relative to other mAChR PAMs. Several compounds, including 8a, c, reduced locomotor activity and increased measures of turning behaviors in zebrafish, suggesting that allosteric modulation of muscarinic receptor efficacy might be useful in the treatment of repetitive behaviors associated with autism spectrum disorder (ASD) and other neuropsychiatric disorders.

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Structural basis for recognition of 26RFa by the pyroglutamylated RFamide peptide receptor

The neuropeptide 26RFa, a member of the RF-amide peptide family, activates the pyroglutamylated RF-amide peptide receptor (QRFPR), a class A GPCR. The 26RFa/QRFPR system plays critical roles in energy homeostasis, making QRFPR an attractive drug target for treating obesity, diabetes, and eating disorders. However, the lack of structural information has hindered our understanding of the peptide recognition and regulatory mechanism of QRFPR, impeding drug design efforts. In this study, we determined the cryo-EM structure of the Gq-coupled QRFPR bound to 26RFa. The structure reveals a unique assembly mode of the extracellular region of the receptor and the N-terminus of the peptide, and elucidates the recognition mechanism of the C-terminal heptapeptide of 26RFa by the transmembrane binding pocket of QRFPR. The study also clarifies the similarities and distinctions in the binding pattern of the RF-amide moiety in five RF-amide peptides and the RY-amide segment in neuropeptide Y. These findings deepen our understanding of the RF-amide peptide recognition, aiding in the rational design of drugs targeting QRFPR and other RF-amide peptide receptors.

Molecular Determinant Underlying Selective Coupling of Primary G-Protein by Class A GPCRs

G-protein-coupled receptors (GPCRs) transmit downstream signals predominantly via G-protein pathways. However, the conformational basis of selective coupling of primary G-protein remains elusive. Histamine receptors H2R and H3R couple with Gs- or Gi-proteins respectively. Here, three cryo-EM structures of H2R-Gs and H3R-Gi complexes are presented at a global resolution of 2.6-2.7 Å. These structures reveal the unique binding pose for endogenous histamine in H3R, wherein the amino group interacts with E2065.46 of H3R instead of the conserved D1143.32 of other aminergic receptors. Furthermore, comparative analysis of the H2R-Gs and H3R-Gi complexes reveals that the structural geometry of TM5/TM6 determines the primary G-protein selectivity in histamine receptors. Machine learning (ML)-based structuromic profiling and functional analysis of class A GPCR-G-protein complexes illustrate that TM5 length, TM5 tilt, and TM6 outward movement are key determinants of the Gs and Gi/o selectivity among the whole Class A family. Collectively, the findings uncover the common structural geometry within class A GPCRs that determines the primary Gs- and Gi/o-coupling selectivity.

Identification of Gα12-vs-Gα13-coupling determinants and development of a Gα12/13-coupled designer GPCR

G-protein-coupled receptors (GPCRs) transduce diverse signals into the cell by coupling to one or several Gα subtypes. Of the 16 Gα subtypes in human cells, Gα12 and Gα13 belong to the G12 subfamily and are reported to be functionally different. Notably, certain GPCRs display selective coupling to either Gα12 or Gα13, highlighting their significance in various cellular contexts. However, the structural basis underlying this selectivity remains unclear. Here, using a Gα12-coupled designer receptor exclusively activated by designer drugs (DREADD; G12D) as a model system, we identified residues in the α5 helix and the receptor that collaboratively determine Gα12-vs-Gα13 selectivity. Residue-swapping experiments showed that G12D distinguishes differences between Gα12 and Gα13 in the positions G.H5.09 and G.H5.23 in the α5 helix. Molecular dynamics simulations observed that I378G.H5.23 in Gα12 interacts with N1032.39, S1693.53 and Y17634.53 in G12D, while H364G.H5.09 in Gα12 interact with Q2645.71 in G12D. Screening of mutations at these positions in G12D identified G12D mutants that enhanced coupling with Gα12 and to an even greater extent with Gα13. Combined mutations, most notably the dual Y17634.53H and Q2645.71R mutant, further enhanced Gα12/13 coupling, thereby serving as a potential Gα12/13-DREADD. Such novel Gα12/13-DREADD may be useful in future efforts to develop drugs that target Gα12/13 signaling as well as to identify their therapeutic indications.

CASTpFold: Computed Atlas of Surface Topography of the universe of protein Folds

Geometric and topological properties of protein structures, including surface pockets, interior cavities, and cross channels, are of fundamental importance for proteins to carry out their functions. Computed Atlas of Surface Topography of proteins (CASTp) is a widely used web server for locating, delineating, and measuring these geometric and topological properties of protein structures. Recent developments in AI-based protein structure prediction such as AlphaFold2 (AF2) have significantly expanded our knowledge on protein structures. Here we present CASTpFold, a continuation of CASTp that provides accurate and comprehensive identifications and quantifications of protein topography. It now provides (i) results on an expanded database of proteins, including the Protein Data Bank (PDB) and non-singleton representative structures of AlphaFold2 structures, covering 183 million AF2 structures; (ii) functional pockets prediction with corresponding Gene Ontology (GO) terms or Enzyme Commission (EC) numbers for AF2-predicted structures; and (iii) pocket similarity search function for surface and protein-protein interface pockets. The CASTpFold web server is freely accessible at https://cfold.bme.uic.edu/castpfold/.

Control of G protein-coupled receptor function via membrane-interacting intrinsically disordered C-terminal domains

G protein-coupled receptors (GPCRs) control intracellular signaling cascades via agonist-dependent coupling to intracellular transducers including heterotrimeric G proteins, GPCR kinases (GRKs), and arrestins. In addition to their critical interactions with the transmembrane core of active GPCRs, all three classes of transducers have also been reported to interact with receptor C-terminal domains (CTDs). An underexplored aspect of GPCR CTDs is their possible role as lipid sensors given their proximity to the membrane. CTD-membrane interactions have the potential to control the accessibility of key regulatory CTD residues to downstream effectors and transducers. Here we report that the CTDs of two closely related family C GPCRs, metabotropic glutamate receptor 2 (mGluR2) and mGluR3, bind to membranes and that this interaction can regulate receptor function. We first characterize CTD structure with NMR spectroscopy, revealing lipid composition-dependent modes of membrane binding. Using molecular dynamics simulations and structure-guided mutagenesis, we then identify key conserved residues and cancer-associated mutations that modulate CTD-membrane binding. Finally, we provide evidence that mGluR3 transducer coupling is controlled by CTD-membrane interactions in live cells, which may be subject to regulation by CTD phosphorylation and changes in membrane composition. This work reveals a novel mechanism of GPCR modulation, suggesting that CTD-membrane binding may be a general regulatory mode throughout the broad GPCR superfamily.

Agonist-selective activation of individual G-proteins by muscarinic receptors

Selective activation of individual subtypes of muscarinic receptors is a promising way to safely alleviate a wide range of pathological conditions in the central nervous system and the periphery as well. The flexible G-protein interface of muscarinic receptors allows them to interact with several G-proteins with various efficacy, potency, and kinetics. Agonists biased to the particular G-protein mediated pathway may result in selectivity among muscarinic subtypes and, due to the non-uniform expression of individual G-protein alpha subunits, possibly achieve tissue specificity. Here, we demonstrate that novel tetrahydropyridine-based agonists exert specific signalling profiles in coupling with individual G-protein α subunits. These signalling profiles profoundly differ from the reference agonist carbachol. Moreover, coupling with individual Gα induced by these novel agonists varies among subtypes of muscarinic receptors which may lead to subtype selectivity. Thus, the novel tetrahydropyridine-based agonist can contribute to the elucidation of the mechanism of pathway-specific activation of muscarinic receptors and serve as a starting point for the development of desired selective muscarinic agonists.