Structural insights into the subtype-selective antagonist binding to the M2 muscarinic receptor

2023 Julius Axelrod Symposium: Plant-Derived Molecules Acting on G Protein-Coupled Receptors

Plant extracts have played a significant role in traditional medicine for centuries, contributing to improved health and the treatment of various human illnesses. G protein-coupled receptors (GPCRs) are crucial in numerous physiologic functions, and there is growing evidence suggesting their involvement in the therapeutic effects of many plant extracts. In recent years, scientists have identified an expanding number of isolated molecules responsible for the biologic activity of these extracts, with many believed to act on GPCRs. This article critically reviews the evidence supporting the modulation of GPCR function by these plant-derived molecules through direct binding. Structural information is now available for some of these molecules, allowing for a comparison of their binding mode with that of endogenous GPCR ligands. The final section explores future trends and challenges, focusing on the identification of new plant-derived molecules with both orthosteric and allosteric binding modes, as well as innovative strategies for designing GPCR ligands inspired by these plant-derived compounds. In conclusion, plant-derived molecules are anticipated to play an increasingly vital role as therapeutic drugs and serve as templates for drug design. SIGNIFICANCE STATEMENT: This minireview summarizes the most pertinent publications on isolated plant-derived molecules interacting with G protein-coupled receptors (GPCRs) and comments on available structural information on GPCR/plant-derived ligand pairs. Future challenges and trends for the isolation and characterization of plant-derived molecules and drug design are discussed.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.

Xanomeline displays concomitant orthosteric and allosteric binding modes at the M4 mAChR

The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.

Role of Conserved Tyrosine Lid Residues in the Activation of the M2 Muscarinic Acetylcholine Receptor

The development of subtype selective small molecule drugs for the muscarinic acetylcholine receptor (mAChR) family has been challenging. The design of more selective ligands can be improved by understanding the structure and function of key amino acid residues that line ligand binding sites. Here we study the role of three conserved key tyrosine residues [Y1043.33, Y4036.51, and Y4267.39 (Ballesteros and Weinstein numbers in superscript)] at the human M2 mAChR, located at the interface between the orthosteric and allosteric binding sites of the receptor. We specifically focused on the role of the three tyrosine hydroxyl groups in the transition between the inactive and active conformations of the receptor by making phenylalanine point mutants. Single-point mutation at either of the three positions was sufficient to reduce the affinity of agonists by ∼100-fold for the M2 mAChR, whereas the affinity of antagonists remained largely unaffected. In contrast, neither of the mutations affected the efficacy of orthosteric agonists. When mutations were combined into double and triple M2 mAChR mutants, the affinity of antagonists was reduced by more than 100-fold compared with the wild-type M2 receptor. In contrast, the affinity of allosteric modulators, either negative or positive, was retained at all single and multiple mutations, but the degree of allosteric effect exerted on the endogenous ligand acetylcholine was affected at all mutants containing Y4267.39F. These findings will provide insights to consider when designing future mAChR ligands. SIGNIFICANCE STATEMENT: Structural studies demonstrated that three tyrosine residues between the orthosteric and allosteric sites of the M2 muscarinic acetylcholine receptor (mAChR) had different hydrogen bonding networks in the inactive and active conformations. The role of hydroxyl groups of the tyrosine residues on orthosteric and allosteric ligand pharmacology was unknown. We found that hydroxyl groups of the tyrosine residues differentially affected the molecular pharmacology of orthosteric and allosteric ligands. These results provide insights to consider when designing future mAChR ligands.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

Cymbopogon citratus Essential Oil: Extraction, GC-MS, Phytochemical Analysis, Antioxidant Activity, and In Silico Molecular Docking for Protein Targets Related to CNS

This study analyzed the chemical composition of Cymbopogon citratus essential oil from Puebla, México, assessed its antioxidant activity, and evaluated in silico protein-compound interactions related to central nervous system (CNS) physiology. GC-MS analysis identified myrcene (8.76%), Z-geranial (27.58%), and E-geranial (38.62%) as the main components, with 45 other compounds present, which depends on the region and growing conditions. DPPH and Folin-Ciocalteu assays using the leaves extract show a promising antioxidant effect (EC₅₀ = 48.5 µL EO/mL), reducing reactive oxygen species. The bioinformatic tool SwissTargetPrediction (STP) shows 10 proteins as potential targets associated with CNS physiology. Moreover, protein-protein interaction diagrams suggest that muscarinic and dopamine receptors are related to each other through a third party. Molecular docking reveals that Z-geranial has higher binding energy than M1 commercial blocker and blocks M2, but not M4 muscarinic acetylcholine receptors, whereas β-pinene and myrcene block M1, M2, and M4 receptors. These actions may positively affect cardiovascular activity, memory, Alzheimer's disease, and schizophrenia. This study highlights the significance of understanding natural product interactions with physiological systems to uncover potential therapeutic agents and advanced knowledge on their benefits for human health.

Density Functional Theory, Molecular Dynamics and AlteQ Studies Approaches of Baimantuoluoamide A and Baimantuoluoamide B to Identify Potential Inhibitors of Mpro Proteins: a Novel Target for the Treatment of SARS COVID-19

COVID-19 has resulted in epidemi conditions over the world. Despite efforts by scientists from all over the world to develop an effective va ine against this virus, there is presently no recognized cure for COVID-19. The most succeed treatments for various ailments come from natural components found in medicinal plants, which are also rucial for the development of new medications. This study intends to understand the role of the baimantuoluoamide A and baimantuoluoamide B molecules in the treatment of Covid19. Initially, density functional theory (DFT) used to explore their electronic potentials along with the Becke3-Lee-Yang-Parr (B3LYP) 6-311 + G(d, p) basis set. A number of characteristics, including the energy gap, hardness, local softness, electronegativity, and electrophilicity, have also been calculated to discuss the reactivity of mole ules. Using natural bond orbital, the title compound's bioactive nature and stability were investigated. Further, both compounds potential inhibitors with main protease (Mpro) proteins, molecular dynamics simulations and AlteQ investigations also studied.

Supplementary Information

The online version contains supplementary material available at 10.1134/S0021364023600039.

Structure-based pharmacophore modeling 2. Developing a novel framework for structure-based pharmacophore model generation and selection

Pharmacophore models are three-dimensional arrangements of molecular features required for biological activity that are used in ligand identification efforts for many biological targets, including G protein-coupled receptors (GPCR). Though GPCR are integral membrane proteins of considerable interest as targets for drug development, many of these receptors lack known ligands or experimentally determined structures necessary for ligand- or structure-based pharmacophore model generation, respectively. Thus, we here present a structure-based pharmacophore modeling approach that uses fragments placed with Multiple Copy Simultaneous Search (MCSS) to generate high-performing pharmacophore models in the context of experimentally determined, as well as modeled GPCR structures. Moreover, we have addressed the oft-neglected topic of pharmacophore model selection via development of a cluster-then-predict machine learning workflow. Herein score-based pharmacophore models were generated in experimentally determined and modeled structures of 13 class A GPCR and resulted in pharmacophore models exhibiting high enrichment factors when used to search a database containing 569 class A GPCR ligands. In addition, classification of pharmacophore models with the best performing cluster-then-predict logistic regression classifier resulted in positive predictive values (PPV) of 0.88 and 0.76 for selecting high enrichment pharmacophore models from among those generated in experimentally determined and modeled structures, respectively.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~ 1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

Apo2ph4: A Versatile Workflow for the Generation of Receptor-based Pharmacophore Models for Virtual Screening

Pharmacophore models are widely used as efficient virtual screening (VS) filters for the target-directed enrichment of large compound libraries. However, the generation of pharmacophore models that have the power to discriminate between active and inactive molecules traditionally requires structural information about ligand-target complexes or at the very least knowledge of one active ligand. The fact that the discovery of the first known active ligand of a newly investigated target represents a major hurdle at the beginning of every drug discovery project underscores the need for methods that are able to derive high-quality pharmacophore models even without the prior knowledge of any active ligand structures. In this work, we introduce a novel workflow, called apo2ph4, that enables the rapid derivation of pharmacophore models solely from the three-dimensional structure of the target receptor. The utility of this workflow is demonstrated retrospectively for the generation of a pharmacophore model for the M2 muscarinic acetylcholine receptor. Furthermore, in order to show the general applicability of apo2ph4, the workflow was employed for all 15 targets of the recently published LIT-PCBA dataset. Pharmacophore-based VS runs using the apo2ph4-derived models achieved a significant enrichment of actives for 13 targets. In the last presented example, a pharmacophore model derived from the etomidate site of the α1β2γ2 GABA_A receptor was used in VS campaigns. Subsequent in vitro testing of selected hits revealed that 19 out of 20 (95%) tested compounds were able to significantly enhance GABA currents, which impressively demonstrates the applicability of apo2ph4 for real-world drug design projects.

Fluorescence-Detection Size-Exclusion Chromatography-Based Thermostability Assay for Membrane Proteins

Green fluorescent proteins (GFPs) have lightened up almost every aspect of biological research including protein sciences. In the field of membrane protein structural biology, GFPs have been used widely to monitor membrane protein localization, expression level, the purification process and yield, and the stability inside the cells and in the test tube. Of particular interest is the fluorescence-detector size-exclusion chromatography-based thermostability assay (FSEC-TS). By simple heating and FSEC, the generally applicable method allows rapid assessment of the thermostability of GFP-fused membrane proteins without purification. Here we describe the experimental details and some typical results for the FSEC-TS method.

Molecular basis for selective activation of DREADD-based chemogenetics

Designer receptors exclusively activated by designer drugs (DREADDs) represent a powerful chemogenetic technology for the remote control of neuronal activity and cellular signalling1-4. The muscarinic receptor-based DREADDs are the most widely used chemogenetic tools in neuroscience research. The Gq-coupled DREADD (hM3Dq) is used to enhance neuronal activity, whereas the Gi/o-coupled DREADD (hM4Di) is utilized to inhibit neuronal activity5. Here we report four DREADD-related cryogenic electron microscopy high-resolution structures: a hM3Dq-miniGq complex and a hM4Di-miniGo complex bound to deschloroclozapine; a hM3Dq-miniGq complex bound to clozapine-N-oxide; and a hM3R-miniGq complex bound to iperoxo. Complemented with mutagenesis, functional and computational simulation data, our structures reveal key details of the recognition of DREADD chemogenetic actuators and the molecular basis for activation. These findings should accelerate the structure-guided discovery of next-generation chemogenetic tools.

A methodology for creating mutants of G-protein coupled receptors stabilized in active state by combining statistical thermodynamics and evolutionary molecular engineering

We challenged the stabilization of a G-protein coupled receptor (GPCR) in the active state solely by multiple amino-acid mutations without the agonist binding. For many GPCRs, the free energy of the active state is higher than that of the inactive state. When the inactive state is stabilized through the lowering of its free energy, the apparent midpoint temperature of thermal denaturation Tm exhibits a significant increase. However, this is not always the case for the stabilization of the active state. We constructed a modified version of our methodology combining statistical thermodynamics and evolutionary molecular engineering, which was recently developed for the inactive state. First, several residues to be mutated are determined using our statistical-thermodynamics theory. Second, a gene (mutant) library is constructed using Escherichia coli cells to efficiently explore most of the mutational space. Third, for the mutant screening, the mutants prepared in accordance with the library are expressed in engineered Saccharomyces cerevisiae YB14 cells which can grow only when a GPCR mutant stabilized in the active state has signaling function. For the adenosine A2A receptor tested, the methodology enabled us to sort out two triple mutants and a double mutant. It was experimentally corroborated that all the mutants exhibit much higher binding affinity for G protein than the wild type. Analyses indicated that the mutations make the structural characteristics shift toward those of the active state. However, only slight increases in Tm resulted from the mutations, suggesting the unsuitability of Tm to the stability measure for the active state.

Multitargeting nature of muscarinic orthosteric agonists and antagonists

Muscarinic receptors (mAChRs) are typical members of the G protein-coupled receptor (GPCR) family and exist in five subtypes from M₁ to M₅. Muscarinic receptor subtypes do not sufficiently differ in affinity to orthosteric antagonists or agonists; therefore, the analysis of receptor subtypes is complicated, and misinterpretations can occur. Usually, when researchers mainly specialized in CNS and peripheral functions aim to study mAChR involvement in behavior, learning, spinal locomotor networks, biological rhythms, cardiovascular physiology, bronchoconstriction, gastrointestinal tract functions, schizophrenia, and Parkinson’s disease, they use orthosteric ligands and they do not use allosteric ligands. Moreover, they usually rely on manufacturers’ claims that could be misleading. This review aimed to call the attention of researchers not deeply focused on mAChR pharmacology to this fact. Importantly, limited selective binding is not only a property of mAChRs but is a general attribute of most neurotransmitter receptors. In this review, we want to give an overview of the most common off-targets for established mAChR ligands. In this context, an important point is a mention the tremendous knowledge gap on off-targets for novel compounds compared to very well-established ligands. Therefore, we will summarize reported affinities and give an outline of strategies to investigate the subtype’s function, thereby avoiding ambiguous results. Despite that, the multitargeting nature of drugs acting also on mAChR could be an advantage when treating such diseases as schizophrenia. Antipsychotics are a perfect example of a multitargeting advantage in treatment. A promising strategy is the use of allosteric ligands, although some of these ligands have also been shown to exhibit limited selectivity. Another new direction in the development of muscarinic selective ligands is functionally selective and biased agonists. The possible selective ligands, usually allosteric, will also be listed. To overcome the limited selectivity of orthosteric ligands, the recommended process is to carefully examine the presence of respective subtypes in specific tissues via knockout studies, carefully apply “specific” agonists/antagonists at appropriate concentrations and then calculate the probability of a specific subtype involvement in specific functions. This could help interested researchers aiming to study the central nervous system functions mediated by the muscarinic receptor.

A Liquid Chromatography-Mass Spectrometry Method to Study the Interaction between Membrane Proteins and Low-Molecular-Weight Compound Mixtures

Molecular interaction analysis is an essential technique for the study of biomolecular functions and the development of new drugs. Most current methods generally require manipulation to immobilize or label molecules, and require advance identification of at least one of the two molecules in the reaction. In this study, we succeeded in detecting the interaction of low-molecular-weight (LMW) compounds with a membrane protein mixture derived from cultured cells expressing target membrane proteins by using the size exclusion chromatography-mass spectrometry (SEC-MS) method under the condition of 0.001% lauryl maltose neopentyl glycol as detergent and atmospheric pressure chemical ionization. This method allowed us to analyze the interaction of a mixture of medicinal herbal ingredients with a mixture of membrane proteins to identify the two interacting ingredients. As it does not require specialized equipment (e.g., a two-dimensional liquid chromatography system), this SEC-MS method enables the analysis of interactions between LMW compounds and relatively high-expressed membrane proteins without immobilization or derivatization of the molecules.

Evolutionary association of receptor-wide amino acids with G protein-coupling selectivity in aminergic GPCRs

Evolutionary analyses of aminergic G protein–coupled receptors reveal receptor-wide potential determinants of G protein–coupling selectivity.

G protein-coupled receptors (GPCRs) induce signal transduction pathways through coupling to four main subtypes of G proteins (G_s, G_i, G_q, and G_12/13), selectively. However, G protein selective activation mechanisms and residual determinants in GPCRs have remained obscure. Herein, we performed extensive phylogenetic analysis and identified specifically conserved residues for the aminergic receptors having similar coupling profiles. By integrating our methodology of differential evolutionary conservation of G protein–specific amino acids with structural analyses, we identified specific activation networks for G_s, G_i1, G_o, and G_q. To validate that these networks could determine coupling selectivity we further analyzed G_s-specific activation network and its association with G_s selectivity. Through molecular dynamics simulations, we showed that previously uncharacterized Glycine at position 7x41 plays an important role in receptor activation and it may determine G_s coupling selectivity by facilitating a larger TM6 movement. Finally, we gathered our results into a comprehensive model of G protein selectivity called “sequential switches of activation” describing three main molecular switches controlling GPCR activation: ligand binding, G protein selective activation mechanisms, and G protein contact.

Activation and allosteric regulation of the orphan GPR88-Gi1 signaling complex

GPR88 is an orphan class A G-protein-coupled receptor that is highly expressed in the striatum and regulates diverse brain and behavioral functions. Here we present cryo-EM structures of the human GPR88-Gi1 signaling complex with or without a synthetic agonist (1R, 2R)-2-PCCA. We show that (1R, 2R)-2-PCCA is an allosteric modulator binding to a herein identified pocket formed by the cytoplasmic ends of transmembrane segments 5, 6, and the extreme C terminus of the α5 helix of Gi1. We also identify an electron density in the extracellular orthosteric site that may represent a putative endogenous ligand of GPR88. These structures, together with mutagenesis studies and an inactive state model obtained from metadynamics simulations, reveal a unique activation mechanism for GPR88 with a set of distinctive structure features and a water-mediated polar network. Overall, our results provide a structural framework for understanding the ligand binding, activation and signaling mechanism of GPR88, and will facilitate the innovative drug discovery for neuropsychiatric disorders and for deorphanization of this receptor.

MD Simulations Revealing Special Activation Mechanism of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor (GPCR) that is gaining much interest for its regulating role in the central nervous system and its value as a drug target. Structures of CB1 in inactive and active states have revealed conformational change details that are not common in other GPCRs. Here, we performed molecular dynamics simulations of CB1 in different ligand binding states and with mutations to reveal its activation mechanism. The conformational change of the “twin toggle switch” residues F200^3.36 and W356^6.48 that correlates with ligand efficacy is identified as a key barrier step in CB1 activation. Similar conformational change of residues 3.36/6.48 is also observed in melanocortin receptor 4, showing this “twin toggle switch” residue pair is crucial for the activation of multiple GPCR members.

Development of an Outward Proton Pumping Rhodopsin with a New Record in Thermostability by Means of Amino Acid Mutations

We have developed a methodology for identifying further thermostabilizing mutations for an intrinsically thermostable membrane protein. The methodology comprises the following steps: (1) identifying thermostabilizing single mutations (TSSMs) for residues in the transmembrane region using our physics-based method; (2) identifying TSSMs for residues in the extracellular and intracellular regions, which are in aqueous environment, using an empirical force field FoldX; and (3) combining the TSSMs identified in steps (1) and (2) to construct multiple mutations. The methodology is illustrated for thermophilic rhodopsin whose apparent midpoint temperature of thermal denaturation T_m is ∼91.8 °C. The TSSMs previously identified in step (1) were F90K, F90R, and Y91I with ΔT_m ∼5.6, ∼5.5, and ∼2.9 °C, respectively, and those in step (2) were V79K, T114D, A115P, and A116E with ΔT_m ∼2.7, ∼4.2, ∼2.6, and ∼2.3 °C, respectively (ΔT_m denotes the increase in T_m). In this study, we construct triple and quadruple mutants, F90K+Y91I+T114D and F90K+Y91I+V79K+T114D. The values of ΔT_m for these multiple mutants are ∼11.4 and ∼13.5 °C, respectively. T_m of the quadruple mutant (∼105.3 °C) establishes a new record in a class of outward proton pumping rhodopsins. It is higher than T_m of Rubrobacter xylanophilus rhodopsin (∼100.8 °C) that was the most thermostable in the class before this study.

Ligand-based G Protein Coupled Receptor pharmacophore modeling: Assessing the role of ligand function in model development

Integral membrane proteins in the G Protein-Coupled Receptor (GPCR) class are attractive drug development targets. However, computational methods applicable to ligand discovery for many GPCR targets are restricted by limited numbers of known ligands. Pharmacophore models can be developed using variously sized training sets and applied in database mining to prioritize candidate ligands for subsequent validation. This in silico study assessed the impact of key pharmacophore modeling decisions that arise when known ligand numbers for a target of interest are low. GPCR included in this study are the adrenergic alpha-1A, 1D and 2A, adrenergic beta 2 and 3, kappa, delta and mu opioid, serotonin 1A and 2A, and the muscarinic 1 and 2 receptors, all of which have rich ligand data sets suitable to assess the performance of protocols intended for application to GPCR with limited ligand data availability. Impact of ligand function, potency and structural diversity in training set selection was assessed to define when pharmacophore modeling targeting GPCR with limited known ligands becomes viable. Pharmacophore elements and pharmacophore model selection criteria were also assessed. Pharmacophore model assessment was based on percent pharmacophore model generation failure, as well as Güner-Henry enrichment and goodness-of-hit scores. Three of seven pharmacophore element schemes evaluated in MOE 2018.0101, Unified, PCHD, and CHD, showed substantially lower failure rates and higher enrichment scores than the others. Enrichment and GH scores were used to compare construction protocol for pharmacophore models of varying purposes- such as function specific versus nonspecific ligand identification. Notably, pharmacophore models constructed from ligands of mixed functions (agonists and antagonists) were capable of enriching hitlists with active compounds, and therefore can be used when available sets of known ligands are limited in number.

Vibrational spectroscopy analysis of ligand efficacy in human M2 muscarinic acetylcholine receptor (M2R)

The intrinsic efficacy of ligand binding to G protein-coupled receptors (GPCRs) reflects the ability of the ligand to differentially activate its receptor to cause a physiological effect. Here we use attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy to examine the ligand-dependent conformational changes in the human M2 muscarinic acetylcholine receptor (M2R). We show that different ligands affect conformational alteration appearing at the C=O stretch of amide-I band in M2R. Notably, ATR-FTIR signals strongly correlated with G-protein activation levels in cells. Together, we propose that amide-I band serves as an infrared probe to distinguish the ligand efficacy in M2R and paves the path to rationally design ligands with varied efficacy towards the target GPCR.

Guanidine Derivatives: How Simple Structural Modification of Histamine H3R Antagonists Has Led to the Discovery of Potent Muscarinic M2R/M4R Antagonists

This article describes the discovery of novel potent muscarinic receptor antagonists identified during a search for more active histamine H₃ receptor (H₃R) ligands. The idea was to replace the flexible seven methylene linker with a semirigid 1,4-cyclohexylene or p-phenylene substituted group of the previously described histamine H₃R antagonists ADS1017 and ADS1020. These simple structural modifications of the histamine H₃R antagonist led to the emergence of additional pharmacological effects, some of which unexpectedly showed strong antagonist potency at muscarinic receptors. This paper reports the routes of synthesis and pharmacological characterization of guanidine derivatives, a novel chemotype of muscarinic receptor antagonists binding to the human muscarinic M₂ and M₄ receptors (hM₂R and hM₄R, respectively) in nanomolar concentration ranges. The affinities of the newly synthesized ADS10227 (1-{4-{4-{[4-(phenoxymethyl)cyclohexyl]methyl}piperazin-1-yl}but-1-yl}-1-(benzyl)guanidine) at hM₂R and hM₄R were 2.8 nM and 5.1 nM, respectively.

Molecular Dynamics Simulations for the Determination of the Characteristic Structural Differences between Inactive and Active States of Wild Type and Mutants of the Orexin2 Receptor

The orexin2 receptor (OX2R), which is classified as a class A G protein-coupled receptor (GPCR), is the target of our study. We performed over 20 several-microsecond-scale molecular dynamics simulations of the wild type and mutants of OX2R to extract the characteristics of the structural changes taking place in the active state. We introduced mutations that exhibited the stable inactive state and the constitutively active state in class A GPCRs. In these simulations, significant characteristic structural changes were observed in the V309^6.40Y mutant, which corresponded to a constitutively active mutant. These conformational changes include the outward movement of the transmembrane helix 6 (TM6) and the inward movement of TM7, which are common structural changes in the activation of GPCRs. In addition, we extracted a suitable index for the quantitative evaluation of the active and inactive states of GPCRs, namely, the inter-atomic distance of Cα atoms between x(3.46) and Y(7.53). The structures of the inactive and active states solved by X-ray crystallography and cryo-electron microscopy can be classified using the inter-atomic distance. Furthermore, we clarified that the inward movement of TM7 requires the swapping of M305^6.36 on TM6 and L367^7.56 on TM7. Finally, we discussed the structural advantages of TM7 inward movement for GPCR activation.

Structure-Activity Relationship Studies on Oxazolo[3,4-a]pyrazine Derivatives Leading to the Discovery of a Novel Neuropeptide S Receptor Antagonist with Potent In Vivo Activity

Neuropeptide S modulates important neurobiological functions including locomotion, anxiety, and drug abuse through interaction with its G protein-coupled receptor known as neuropeptide S receptor (NPSR). NPSR antagonists are potentially useful for the treatment of substance abuse disorders against which there is an urgent need for new effective therapeutic approaches. Potent NPSR antagonists in vitro have been discovered which, however, require further optimization of their in vivo pharmacological profile. This work describes a new series of NPSR antagonists of the oxazolo[3,4-a]pyrazine class. The guanidine derivative 16 exhibited nanomolar activity in vitro and 5-fold improved potency in vivo compared to SHA-68, a reference pharmacological tool in this field. Compound 16 can be considered a new tool for research studies on the translational potential of the NPSergic system. An in-depth molecular modeling investigation was also performed to gain new insights into the observed structure-activity relationships and provide an updated model of ligand/NPSR interactions.

Allosteric communication regulates ligand-specific GPCR activity

G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.

Dibenzodiazepinone-type muscarinic receptor antagonists conjugated to basic peptides: Impact of the linker moiety and unnatural amino acids on M2R selectivity

The family of human muscarinic acetylcholine receptors (MRs) is characterized by a high sequence homology among the five subtypes (M₁R-M₅R), being the reason for a lack of subtype selective MR ligands. In continuation of our work on dualsteric dibenzodiazepinone-type M₂R antagonists, a series of M₂R ligands containing a dibenzodiazepinone pharmacophore linked to small basic peptides was synthesized (64 compounds). The linker moiety was varied with respect to length, number of basic nitrogens (0-2) and flexibility. Besides proteinogenic basic amino acids (Lys, Arg), shorter homologues of Lys and Arg, containing three and two methylene groups, respectively, as well as D-configured amino acids were incorporated. The type of linker had a marked impact on M₂R affinity and also effected M₂R selectivity. In contrast, the structure of the basic peptide rather determined M₂R selectivity than M₂R affinity. For example, the most M₂R selective compound (UR-CG188, 89) with picomolar M₂R affinity (pK_i 9.60), exhibited a higher M₂R selectivity (ratio of K_i M₁R/M₂R/M₃R/M₄R/M₅R: 110:1:5200:55:2300) compared to the vast majority of reported M₂R preferring MR ligands. For selected ligands, M₂R antagonism was confirmed in a M₂R miniG protein recruitment assay.

Vibrational analysis of acetylcholine binding to the M2 receptor

The M₂ muscarinic acetylcholine receptor (M₂R) is a prototypical G protein-coupled receptor (GPCR) that responds to acetylcholine (ACh) and mediates various cellular responses in the nervous system. We recently established Attenuated Total Reflection-Fourier Transform Infrared (ATR-FTIR) spectroscopy for ligand binding to M₂R reconstituted in lipid membranes, paving the way to understand the mechanism in atomic detail. However, the obtained difference FTIR spectra upon ligand binding contained ligand, protein, lipid, and water signals, so a vibrational assignment was needed for a thorough understanding. In the present study, we compared difference FTIR spectra between unlabeled and 2-¹³C labeled ACh, and assigned the bands at 1741 and 1246 cm⁻¹ as the C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O and C–O stretches of ACh, respectively. The CO stretch of ACh in M₂R is close to that in aqueous solution (1736 cm⁻¹), and much lower in frequency than the free CO stretch (1778–1794 cm⁻¹), indicating a strong hydrogen bond, which probably formed with N404^6.52. We propose that a water molecule bridges ACh and N404^6.52. The other ACh terminal is positively charged, and it interacts with negatively charged D103^3.32. The present study revealed that D103^3.32 is deprotonated (negatively charged) in both ACh-bound and free states, a suggested mechanism to stabilize the negative charge of D103^3.32 in the free M₂R.

We recently reported difference FTIR spectra upon binding of Ach to M₂R.

Exploration of stilbenoid trimers as potential inhibitors of sirtuin1 enzyme using a molecular docking and molecular dynamics simulation approach

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein. SIRT1 expression regulates cellular stress responses that lead to the development of cancer. Redocking showed a good native ligand pose with an RMSD value of 1.40 Å at the receptor active site's coordinates. The molecular docking score uses a grid score functional (kcal mol⁻¹), which shows results of 1NS: 79.56, TS1: −26.83, TS2: −87.77, and TS3: −83.67. The TS2 and TS3 candidates were chosen for further analysis because they had a lower grid score than the native ligand (1NS). Furthermore, prediction of binding free energy (kcal mol⁻¹) using the Quantum Mechanics/generalized Born Surface Area (QM/MM-GBSA) method shows the results of 1NS: −31.52 ± 0.39, TS2: −58.99 ± 0.34, and TS3: −43.38 ± 0.35. These results indicate that the TS2 and TS3 compounds have good potential as inhibitors of the SIRT1 enzyme. Additionally, the amino acid residues were responsible for the inhibition mechanism through hydrogen bond interactions at the molecular level, including ASP22, PHE91, PRO11, ILE165, ASP166, and VAL230. The observations made in this study provide theoretical information for exploring the stilbenoid trimers as anticancer agents by targeting the SIRT1 enzyme.

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

The optimization of a novel selective antagonist for human M2 muscarinic acetylcholine receptor

Muscarinic acetylcholine receptors (mAChRs) comprise five distinct subtypes denoted M₁ to M₅. The antagonism of M₂ subtype could increase the release of acetylcholine from vesicles into the synaptic cleft and improve postsynaptic functions in the hippocampus via M₁ receptor activation, displaying therapeutic potentials for Alzheimer's disease. However, drug development for M₂ antagonists is still challenged among different receptor subtypes. In this study, by optimizing a scaffold from virtual screening, we synthesized two focused libraries and generated up to 50 derivatives. By measuring potency and binding selectivity, we discovered a novel M₂ antagonist, ligand 47, featuring submicromolar IC₅₀, high M₂/M₄ selectivity (~30-fold) and suitable lipophilicity (cLogP = 4.55). Further study with these compounds also illustrates the structure-activity relationship of this novel scaffold. Our study could not only provide novel lead structure, which was easy to synthesize, but also offer valuable information for further development of selective M₂ ligands.

Structure of the human gonadotropin-releasing hormone receptor GnRH1R reveals an unusual ligand binding mode

Gonadotrophin-releasing hormone (GnRH), also known as luteinizing hormone-releasing hormone, is the main regulator of the reproductive system, acting on gonadotropic cells by binding to the GnRH1 receptor (GnRH1R). The GnRH-GnRH1R system is a promising therapeutic target for maintaining reproductive function; to date, a number of ligands targeting GnRH1R for disease treatment are available on the market. Here, we report the crystal structure of GnRH1R bound to the small-molecule drug elagolix at 2.8 Å resolution. The structure reveals an interesting N-terminus that could co-occupy the enlarged orthosteric binding site together with elagolix. The unusual ligand binding mode was further investigated by structural analyses, functional assays and molecular docking studies. On the other hand, because of the unique characteristic of lacking a cytoplasmic C-terminal helix, GnRH1R exhibits different microswitch structural features from other class A GPCRs. In summary, this study provides insight into the ligand binding mode of GnRH1R and offers an atomic framework for rational drug design.

Sodium ions allosterically modulate the M2 muscarinic receptor

G protein coupled receptors (GPCRs) play a key role in the vast majority of cellular signal transduction processes. Previous experimental evidence has shown that sodium ion (Na⁺) allosterically modulate several class A GPCRs and theoretical studies suggested that the same also holds true for muscarinic receptors. Here we examined, using Xenopus oocytes as an expression system, the effect of Na⁺ on a prototypical GPCR, the M2 muscarinic receptor (M2R). We found that removal of extracellular Na⁺ resulted in a decrease in the potency of ACh toward the M2R and that a conserved aspartate in transmembrane domain 2 is crucial for this effect. We further show that this allosteric effect of Na⁺ does not underlie the voltage-dependence of this receptor.

Theoretical identification of thermostabilizing amino acid mutations for G-protein-coupled receptors

Thermostabilization of a membrane proteins, especially G-protein-coupled receptors (GPCRs), is often necessary for biochemical applications and pharmaceutical studies involving structure-based drug design. Here we review our theoretical, physics-based method for identifying thermostabilizing amino acid mutations. Its novel aspects are the following: The entropic effect originating from the translational displacement of hydrocarbon groups within the lipid bilayer is treated as a pivotal factor; a reliable measure of thermostability is introduced and a mutation which enlarges the measure to a significant extent is chosen; and all the possible mutations can be examined with moderate computational effort. It was shown that mutating the residue at a position of N_BW = 3.39 (N_BW is the Ballesteros-Weinstein number) to Arg or Lys leads to the stabilization of significantly many different GPCRs of class A in the inactive state. Up to now, we have been successful in stabilizing several GPCRs and newly solving three-dimensional structures for the muscarinic acetylcholine receptor 2 (M2R), prostaglandin E receptor 4 (EP4), and serotonin 2A receptor (5-HT_2AR) using X-ray crystallography. The subjects to be pursued in future studies are also discussed.

Methodology for Further Thermostabilization of an Intrinsically Thermostable Membrane Protein Using Amino Acid Mutations with Its Original Function Being Retained

We develop a new methodology best suited to the identification of thermostabilizing mutations for an intrinsically stable membrane protein. The recently discovered thermophilic rhodopsin, whose apparent midpoint temperature of thermal denaturation T_m is measured to be ∼91.8 °C, is chosen as a paradigmatic target. In the methodology, we first regard the residues whose side chains are missing in the crystal structure of the wild type (WT) as the "residues with disordered side chains," which make no significant contributions to the stability, unlike the other essential residues. We then undertake mutating each of the residues with disordered side chains to another residue except Ala and Pro, and the resultant mutant structure is constructed by modifying only the local structure around the mutated residue. This construction is based on the postulation that the structure formed by the other essential residues, which is nearly optimized in such a highly stable protein, should not be modified. The stability changes arising from the mutations are then evaluated using our physics-based free-energy function (FEF). We choose the mutations for which the FEF is much lower than for the WT and test them by experiments. We successfully find three mutants that are significantly more stable than the WT. A double mutant whose T_m reaches ∼100 °C is also discovered.

Agonist-Specific Conformations of the M2 Muscarinic Acetylcholine Receptor Assessed by Molecular Dynamics

Binding of muscarinic ligands, both antagonists and agonists, and their effects on the conformation of the M₂ acetylcholine receptor were modeled in silico and compared to experimental data. After docking of antagonists to the M₂ receptor in an inactive conformation (3UON, 5ZK3, 5ZKB, or 5ZKB) and agonists in an active conformation (4MQS), 100 ns of conventional molecular dynamics (MD) followed by 500 ns of accelerated MD was run. Conventional MD revealed ligand-specific interactions with the receptor. Antagonists stabilized the receptor in an inactive conformation during accelerated MD. The receptor in complex with various agonists attained different conformations specific to individual agonists. The magnitude of the TM6 movement correlated with agonist efficacy at the non-preferential G_s pathway. The shape of the intracellular opening where the receptor interacts with a G-protein was different for the classical agonist carbachol, super-agonist iperoxo, and G_i/o-biased partial agonists JR-6 and JR-7, being compatible with experimentally observed agonist bias at the G-protein level. Moreover, a wash-resistant binding of the unique agonist xanomeline associated with interactions with membrane lipids was formed during accelerated MD. Thus, accelerated MD is suitable for modeling of ligand-specific receptor binding and receptor conformations that is essential for the design of experiments aimed at identification of the secondary binding sites and understanding molecular mechanisms underlying receptor activation.

The structural study of mutation-induced inactivation of human muscarinic receptor M4

Human muscarinic receptor M4 belongs to the class A subfamily of the G-protein-coupled receptors (GPCRs). M4 has emerged as an attractive drug target for the treatment of Alzheimer's disease and schizophrenia. Recent results showed that M4-mediated cholinergic transmission is related to motor symptoms in Parkinson's disease. Selective ligand design for the five muscarinic acetylcholine receptor (mAchR) subtypes currently remains challenging owing to the high sequence and structural similarity of their orthosteric binding pockets. In order to obtain M4-selective antagonists, a new approach was tried to lock M4 into an inactive form by rationally designing an N4497.49R mutation, which mimics the allosteric sodium binding in the conserved sodium site usually found in class A GPCRs. In addition, the crystal structure of the mutation-induced inactive M4 was determined. By comparative analysis with other mAchR structures, followed by functional assays, the N4497.49R mutation was shown to stabilize M4 into an inactive state. Virtual screening of a focused ligand library using the crystal structure showed that the inactive M4 prefers antagonists much more than agonists. This study provides a powerful mutation strategy to stabilize GPCRs in inactive states and facilitate their structure determination.

G Protein-Coupled Receptors in Asthma Therapy: Pharmacology and Drug Action

Asthma is a heterogeneous inflammatory disease of the airways that is associated with airway hyperresponsiveness and airflow limitation. Although asthma was once simply categorized as atopic or nonatopic, emerging analyses over the last few decades have revealed a variety of asthma endotypes that are attributed to numerous pathophysiological mechanisms. The classification of asthma by endotype is primarily routed in different profiles of airway inflammation that contribute to bronchoconstriction. Many asthma therapeutics target G protein-coupled receptors (GPCRs), which either enhance bronchodilation or prevent bronchoconstriction. Short-acting and long-acting β 2-agonists are widely used bronchodilators that signal through the activation of the β 2-adrenergic receptor. Short-acting and long-acting antagonists of muscarinic acetylcholine receptors are used to reduce bronchoconstriction by blocking the action of acetylcholine. Leukotriene antagonists that block the signaling of cysteinyl leukotriene receptor 1 are used as an add-on therapy to reduce bronchoconstriction and inflammation induced by cysteinyl leukotrienes. A number of GPCR-targeting asthma drug candidates are also in different stages of development. Among them, antagonists of prostaglandin D2 receptor 2 have advanced into phase III clinical trials. Others, including antagonists of the adenosine A2B receptor and the histamine H4 receptor, are in early stages of clinical investigation. In the past decade, significant research advancements in pharmacology, cell biology, structural biology, and molecular physiology have greatly deepened our understanding of the therapeutic roles of GPCRs in asthma and drug action on these GPCRs. This review summarizes our current understanding of GPCR signaling and pharmacology in the context of asthma treatment. SIGNIFICANCE STATEMENT: Although current treatment methods for asthma are effective for a majority of asthma patients, there are still a large number of patients with poorly controlled asthma who may experience asthma exacerbations. This review summarizes current asthma treatment methods and our understanding of signaling and pharmacology of G protein-coupled receptors (GPCRs) in asthma therapy, and discusses controversies regarding the use of GPCR drugs and new opportunities in developing GPCR-targeting therapeutics for the treatment of asthma.

Crystal structure of the M5 muscarinic acetylcholine receptor

The 5 subtypes of the muscarinic acetylcholine receptors (mAChRs) are expressed throughout the central and peripheral nervous system where they play a vital role in physiology and pathologies. Recently, the M₅ mAChR subtype has emerged as an exciting drug target for the treatment of drug addiction. We have determined the atomic structure of the M₅ mAChR bound to the clinically used inverse agonist tiotropium. The M₅ mAChR structure now allows for a full comparison of all 5 mAChR subtypes and reveals that small differences in the extracellular loop regions can mediate orthosteric and allosteric ligand selectivity. Together, these findings open the door for future structure-based design of selective drugs that target this therapeutically important class of receptors.

The human M₅ muscarinic acetylcholine receptor (mAChR) has recently emerged as an exciting therapeutic target for treating a range of disorders, including drug addiction. However, a lack of structural information for this receptor subtype has limited further drug development and validation. Here we report a high-resolution crystal structure of the human M₅ mAChR bound to the clinically used inverse agonist, tiotropium. This structure allowed for a comparison across all 5 mAChR family members that revealed important differences in both orthosteric and allosteric sites that could inform the rational design of selective ligands. These structural studies, together with chimeric swaps between the extracellular regions of the M₂ and M₅ mAChRs, provided structural insight into kinetic selectivity, where ligands show differential residency times between related family members. Collectively, our study provides important insights into the nature of orthosteric and allosteric ligand interaction across the mAChR family that could be exploited for the design of selective drugs.

Ligand Binding-Induced Structural Changes in the M2 Muscarinic Acetylcholine Receptor Revealed by Vibrational Spectroscopy

M₂ muscarinic acetylcholine receptor (M₂R) is a prototypical G protein-coupled receptor (GPCR) that responds to acetylcholine and mediates various cellular responses in the nervous system. Here, we used attenuated total reflection-Fourier transform infrared spectroscopy analyses on M₂R reconstituted in a lipid membrane to understand the molecular mechanism behind the ligand binding-induced conformational changes. Upon agonist binding, M₂R shows large spectral change of the amide-I band corresponding to backbone C═O stretch, which likely connects with the receptor activation in the lipid environment. These results pave the way to probe effects of different ligand binding on GPCRs using vibrational spectroscopy.

Structure of an allosteric modulator bound to the CB1 cannabinoid receptor

The CB1 receptor mediates the central nervous system response to cannabinoids, and is a drug target for pain, anxiety and seizures. CB1 also responds to allosteric modulators, which influence cannabinoid binding and efficacy. To understand the mechanism of these compounds, we solved the crystal structure of CB1 with the negative allosteric modulator (NAM) ORG27569 and the agonist CP55940. The structure reveals that the NAM binds to an extrahelical site within the inner leaflet of the membrane, which overlaps with a conserved site of cholesterol interaction in many G protein-coupled receptors (GPCRs). The ternary structure with ORG27569 and CP55940 captures an intermediate state of the receptor, in which aromatic residues at the base of the agonist-binding pocket adopt an inactive conformation despite the large contraction of the orthosteric pocket. The structure illustrates a potential strategy for drug modulation of CB1 and other class A GPCRs.

Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors

Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.

Harnessing Ion-Binding Sites for GPCR Pharmacology

Endogenous ions play important roles in the function and pharmacology of G-protein coupled receptors (GPCRs). Historically the evidence for ionic modulation of GPCR function dates to 1973 with studies of opioid receptors, where it was demonstrated that physiologic concentrations of sodium allosterically attenuated agonist binding. This Na+-selective effect was distinct from effects of other monovalent and divalent cations, with the latter usually counteracting sodium's negative allosteric modulation of binding. Since then, numerous studies documenting the effects of mono- and divalent ions on GPCR function have been published. While ions can act selectively and nonselectively at many sites in different receptors, the discovery of the conserved sodium ion site in class A GPCR structures in 2012 revealed the unique nature of Na+ site, which has emerged as a near-universal site for allosteric modulation of class A GPCR structure and function. In this review, we synthesize and highlight recent advances in the functional, biophysical, and structural characterization of ions bound to GPCRs. Taken together, these findings provide a molecular understanding of the unique roles of Na+ and other ions as GPCR allosteric modulators. We will also discuss how this knowledge can be applied to the redesign of receptors and ligand probes for desired functional and pharmacological profiles. SIGNIFICANCE STATEMENT: The function and pharmacology of GPCRs strongly depend on the presence of mono and divalent ions in experimental assays and in living organisms. Recent insights into the molecular mechanism of this ion-dependent allosterism from structural, biophysical, biochemical, and computational studies provide quantitative understandings of the pharmacological effects of drugs in vitro and in vivo and open new avenues for the rational design of chemical probes and drug candidates with improved properties.

Allostery in G protein-coupled receptors investigated by molecular dynamics simulations

G-protein-coupled receptors (GPCRs) are allosteric signaling machines that trigger distinct functional responses depending on the particular conformational state they adopt upon binding. This so-called GPCR functional selectivity is prompted by ligands of different efficacy binding at orthosteric or allosteric sites on the receptor, as well as by interactions with intracellular protein partners or other receptor types. Molecular dynamics (MD) simulations can provide important mechanistic, thermodynamic, and kinetic insights into these interactions at a level of molecular detail that is necessary to rightly inform modern drug discovery. Here, we review the most recent MD contributions to understanding GPCR allostery, with an emphasis on their strengths and limitations.