CB1 Cannabinoid Receptor Signaling and Biased Signaling

The CB1 cannabinoid receptor is a G-protein coupled receptor highly expressed throughout the central nervous system that is a promising target for the treatment of various disorders, including anxiety, pain, and neurodegeneration. Despite the wide therapeutic potential of CB1, the development of drug candidates is hindered by adverse effects, rapid tolerance development, and abuse potential. Ligands that produce biased signaling—the preferential activation of a signaling transducer in detriment of another—have been proposed as a strategy to dissociate therapeutic and adverse effects for a variety of G-protein coupled receptors. However, biased signaling at the CB1 receptor is poorly understood due to a lack of strongly biased agonists. Here, we review studies that have investigated the biased signaling profile of classical cannabinoid agonists and allosteric ligands, searching for a potential therapeutic advantage of CB1 biased signaling in different pathological states. Agonist and antagonist bound structures of CB1 and proposed mechanisms of action of biased allosteric modulators are used to discuss a putative molecular mechanism for CB1 receptor activation and biased signaling. Current studies suggest that allosteric binding sites on CB1 can be explored to yield biased ligands that favor or hinder conformational changes important for biased signaling.     

Identification of a Conformational Equilibrium That Determines the Efficacy and Functional Selectivity of the μ‐Opioid Receptor

G‐protein‐coupled receptor (GPCR) ligands impart differing degrees of signaling in the G‐protein and arrestin pathways, in phenomena called “biased signaling”. However, the mechanism underlying the biased signaling of GPCRs is still unclear, although crystal structures of GPCRs bound to the G protein or arrestin are available. In this study, we observed the NMR signals from methionine residues of the μ‐opioid receptor (μOR) in the balanced‐ and biased‐ligand‐bound states. We found that the intracellular cavity of μOR exists in an equilibrium between closed and multiple open conformations with coupled conformational changes on the transmembrane helices 3, 5, 6, and 7, and that the population of each open conformation determines the G‐protein‐ and arrestin‐mediated signaling levels in each ligand‐bound state. These findings provide insight into the biased signaling of GPCRs and will be helpful for development of analgesics that stimulate μOR with reduced tolerance and dependence.     

Activation of G-protein-coupled receptors in cell-derived plasma membranes supported on porous beads

G-protein-coupled receptors (GPCRs) are ubiquitous mediators of signal transduction across cell membranes and constitute a very important class of therapeutic targets. In order to study the complex biochemical signaling network coupling to the intracellular side of GPCRs, it is necessary to engineer and control the downstream signaling components, which is difficult to realize in living cells. We have developed a bioanalytical platform enabling the study of GPCRs in their native membrane transferred inside-out from live cells to lectin-coated beads, with both membrane sides of the receptor being accessible for molecular interactions. Using heterologously expressed adenosine A(2A) receptor carrying a yellow fluorescent protein, we showed that the tethered membranes comprised fully functional receptors in terms of ligand and G protein binding. The interactions between the different signaling partners during the formation and subsequent dissociation of the ternary signaling complex on single beads could be observed in real time using multicolor fluorescence microscopy. This approach of tethering inside-out native membranes accessible from both sides is straightforward and readily applied to other transmembrane proteins. It represents a generic platform suitable for ensemble as well as single-molecule measurements to investigate signaling processes at plasma membranes.     

Pharmacological Characterization of µ-Opioid Receptor Agonists with Biased G Protein or β-Arrestin Signaling,and Computational Study of Conformational Changes during Receptor Activation

In recent years, G protein vs. β-arrestin biased agonism at opioid receptors has been proposed as an opportunity to produce antinociception with reduced adverse effects. However, at present this approach is highly debated, a reason why more information about biased ligands is required. While the practical relevance of bias in the case of µ-opioid receptors (MOP) still needs to be validated, it remains important to understand the basis of this bias of MOP (and other GPCRs). Recently, we reported two cyclopeptides with high affinity for MOP, the G protein biased Dmt-c[d-Lys-Phe-pCF₃-Phe-Asp]NH₂ (F-81), and the β-arrestin 2 biased Dmt-c[d-Lys-Phe-Asp]NH₂ (C-33), as determined by calcium mobilization assay and bioluminescence resonance energy transfer-based assay. The biased character of F-81 and C-33 has been further analyzed in the [³⁵S]GTPγS binding assay in human MOP-expressing cells, and the PathHunter enzyme complementation assay, used to measure β-arrestin 2 recruitment. To investigate the structural features of peptide-MOP complexes, we performed conformational analysis by NMR spectroscopy, molecular docking, and molecular dynamics simulation. These studies predicted that the two ligands form alternative complexes with MOP, engaging specific ligand–receptor contacts. This would induce different displays of the cytosolic side of the seven-helices bundle, in particular by stabilizing different angulations of helix 6, that could favor intracellular coupling to either G protein or β-arrestin.     

Differences between G-Protein-Stabilized Agonist–GPCR Complexes and their Nanobody-Stabilized Equivalents

Protein nanobodies have been used successfully as surrogates for unstable G-proteins in order to crystallize G-protein-coupled receptors (GPCRs) in their active states. We used molecular dynamics (MD) simulations, including metadynamics enhanced sampling, to investigate the similarities and differences between GPCR–agonist ternary complexes with the α-subunits of the appropriate G-proteins and those with the protein nanobodies (intracellular binding partners, IBPs) used for crystallization. In two of the three receptors considered, the agonist-binding mode differs significantly between the two alternative ternary complexes. The ternary-complex model of GPCR activation entails enhancement of ligand binding by bound IBPs: Our results show that IBP-specific changes can alter the agonist binding modes and thus also the criteria for designing GPCR agonists.     

A Photoswitchable Dualsteric Ligand Controlling Receptor Efficacy

The investigation of the mode and time course of the activation of G-protein-coupled receptors (GPCRs), in particular muscarinic acetylcholine (mACh or M) receptors, is still in its infancy despite the tremendous therapeutic relevance of M receptors and GPCRs in general. We herein made use of a dualsteric ligand that can concomitantly interact with the orthosteric, that is, the neurotransmitter, binding site and an allosteric one. We synthetically incorporated a photoswitchable (photochromic) azobenzene moiety. We characterized the photophysical properties of this ligand called BQCAAI and investigated its applicability as a pharmacological tool compound with a set of FRET techniques at the M₁ receptor. BQCAAI proved to be an unprecedented molecular tool; it is the first photoswitchable dualsteric ligand, and its activity can be regulated by light. We also applied BQCCAI to investigate the time course of several receptor activation processes.     

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ～75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.     

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.     

Mechanistic Studies on the Stereoselectivity of the Serotonin 5‐HT1A Receptor

G‐protein‐coupled receptors (GPCRs) are involved in a wide range of physiological processes, and they have attracted considerable attention as important targets for developing new medicines. A central and largely unresolved question in drug discovery, which is especially relevant to GPCRs, concerns ligand selectivity: Why do certain molecules act as activators (agonists) whereas others, with nearly identical structures, act as blockers (antagonists) of GPCRs? To address this question, we employed all‐atom, long‐timescale molecular dynamics simulations to investigate how two diastereomers (epimers) of dihydrofuroaporphine bind to the serotonin 5‐HT_1A receptor and exert opposite effects. By using molecular interaction fingerprints, we discovered that the agonist could mobilize nearby amino acid residues to act as molecular switches for the formation of a continuous water channel. In contrast, the antagonist epimer remained firmly stabilized in the binding pocket.     

Molecular dynamics simulations and structure-based network analysis reveal structural and functional aspects of G-protein coupled receptor dimer interactions

A significant amount of experimental evidence suggests that G-protein coupled receptors (GPCRs) do not act exclusively as monomers but also form biologically relevant dimers and oligomers. However, the structural determinants, stoichiometry and functional importance of GPCR oligomerization remain topics of intense speculation. In this study we attempted to evaluate the nature and dynamics of GPCR oligomeric interactions. A representative set of GPCR homodimers were studied through Coarse-Grained Molecular Dynamics simulations, combined with interface analysis and concepts from network theory for the construction and analysis of dynamic structural networks. Our results highlight important structural determinants that seem to govern receptor dimer interactions. A conserved dynamic behavior was observed among different GPCRs, including receptors belonging in different GPCR classes. Specific GPCR regions were highlighted as the core of the interfaces. Finally, correlations of motion were observed between parts of the dimer interface and GPCR segments participating in ligand binding and receptor activation, suggesting the existence of mechanisms through which dimer formation may affect GPCR function. The results of this study can be used to drive experiments aimed at exploring GPCR oligomerization, as well as in the study of transmembrane protein–protein interactions in general.     

Unwinding of the C‐Terminal Residues of Neuropeptide Y is critical for Y2 Receptor Binding and Activation

Despite recent breakthroughs in the structural characterization of G‐protein‐coupled receptors (GPCRs), there is only sparse data on how GPCRs recognize larger peptide ligands. NMR spectroscopy, molecular modeling, and double‐cycle mutagenesis studies were integrated to obtain a structural model of the peptide hormone neuropeptide Y (NPY) bound to its human G‐protein‐coupled Y₂ receptor (Y₂R). Solid‐state NMR measurements of specific isotope‐labeled NPY in complex with in vitro folded Y₂R reconstituted into phospholipid bicelles provided the bioactive structure of the peptide. Guided by solution NMR experiments, it could be shown that the ligand is tethered to the second extracellular loop by hydrophobic contacts. The C‐terminal α‐helix of NPY, which is formed in a membrane environment in the absence of the receptor, is unwound starting at T³² to provide optimal contacts in a deep binding pocket within the transmembrane bundle of the Y₂R.     

Influencing receptor-ligand binding mechanisms with multivalent ligand architecture

Multivalent ligands can function as inhibitors or effectors of biological processes. Potent inhibitory activity can arise from the high functional affinities of multivalent ligand-receptor interactions. Effector functions, however, are influenced not only by apparent affinities but also by alternate factors, including the ability of a ligand to cluster receptors. Little is known about the molecular features of a multivalent ligand that determine whether it will function as an inhibitor or effector. We envisioned that, by altering multivalent ligand architecture, ligands with preferences for different binding mechanisms would be generated. To this end, a series of 28 ligands possessing structural diversity was synthesized. This series provides the means to explore the effects of ligand architecture on the inhibition and clustering of a model protein, the lectin concanavalin A (Con A). The structural parameters that were varied include scaffold shape, size, valency, and density of binding elements. We found that ligands with certain architectures are effective inhibitors, but others mediate receptor clustering. Specifically, high molecular weight, polydisperse polyvalent ligands are effective inhibitors of Con A binding, whereas linear oligomeric ligands generated by the ring-opening metathesis polymerization have structural properties that favor clustering. The shape of a multivalent ligand also influences specific aspects of receptor clustering. These include the rate at which the receptor is clustered, the number of receptors in the clusters, and the average interreceptor distance. Our results indicate that the architecture of a multivalent ligand is a key parameter in determining its activity as an inhibitor or effector. Diversity-oriented syntheses of multivalent ligands coupled with effective assays that can be used to compare the contributions of different binding parameters may afford ligands that function by specific mechanisms.     

Molecular Alliance—From Orthosteric and Allosteric Ligands to Dualsteric/Bitopic Agonists at G Protein Coupled Receptors

Cell‐membrane‐spanning G protein coupled receptors (GPCRs) belong to the most important therapeutic target structures. Endogenous transmitters bind from the outer side of the membrane to the “orthosteric” binding site either deep in the binding pocket or at the extracellular N‐terminal end of the receptor protein. Exogenous modulators that utilize a different, “allosteric”, binding site unveil a pathway to receptor subtype‐selectivity. However, receptor activation through the orthosteric area is often more powerful. Recently there has been evidence that orthosteric/allosteric, in other words “dualsteric”, hybrid compounds unite subtype selectivity and receptor activation. These “bitopic” modulators channelreceptor activation and subsequent intracellular signaling into a subset of possible routes. This concept offers access to GPCR modulators with an unprecedented receptor‐subtype and signaling selectivity profile and, as a consequence, to drugs with fewer side effects.     

Design of Drug Efficacy Guided by Free Energy Simulations of the β2-Adrenoceptor

G-protein-coupled receptors (GPCRs) play important roles in physiological processes and are modulated by drugs that either activate or block signaling. Rational design of the pharmacological efficacy profiles of GPCR ligands could enable the development of more efficient drugs, but is challenging even if high-resolution receptor structures are available. We performed molecular dynamics simulations of the β₂ adrenergic receptor in active and inactive conformations to assess if binding free energy calculations can predict differences in ligand efficacy for closely related compounds. Previously identified ligands were successfully classified into groups with comparable efficacy profiles based on the calculated shift in ligand affinity upon activation. A series of ligands were then predicted and synthesized, leading to the discovery of partial agonists with nanomolar potencies and novel scaffolds. Our results demonstrate that free energy simulations enable design of ligand efficacy and the same approach can be applied to other GPCR drug targets.     

Differences between G‐Protein‐Stabilized Agonist–GPCR Complexes and their Nanobody‐Stabilized Equivalents

Protein nanobodies have been used successfully as surrogates for unstable G‐proteins in order to crystallize G‐protein‐coupled receptors (GPCRs) in their active states. We used molecular dynamics (MD) simulations, including metadynamics enhanced sampling, to investigate the similarities and differences between GPCR–agonist ternary complexes with the α‐subunits of the appropriate G‐proteins and those with the protein nanobodies (intracellular binding partners, IBPs) used for crystallization. In two of the three receptors considered, the agonist‐binding mode differs significantly between the two alternative ternary complexes. The ternary‐complex model of GPCR activation entails enhancement of ligand binding by bound IBPs: Our results show that IBP‐specific changes can alter the agonist binding modes and thus also the criteria for designing GPCR agonists.     

Monitoring the Reversibility of GPCR Signaling by Combining Photochromic Ligands with Label-free Impedance Analysis

G protein-coupled cell surface receptors (GPCR) trigger complex intracellular signaling cascades upon agonist binding. Classic pharmacological assays provide information about binding affinities, activation or blockade at different stages of the signaling cascade, but real time dynamics and reversibility of these processes remain often disguised. We show that combining photochromic NPY receptor ligands, which can be toggled in their receptor activation ability by irradiation with light of different wavelengths, with whole cell label-free impedance assays allows observing the cell response to receptor activation and its reversibility over time. The concept demonstrated on NPY receptors may be well applicable to many other GPCRs providing a deeper insight into the time course of intracellular signaling processes.     

Plasmon Waveguide Resonance: Principles,Applications and Historical Perspectives on Instrument Development

Plasmon waveguide resonance (PWR) is a variant of surface plasmon resonance (SPR) that was invented about two decades ago at the University of Arizona. In addition to the characterization of the kinetics and affinity of molecular interactions, PWR possesses several advantages relative to SPR, namely, the ability to monitor both mass and structural changes. PWR allows anisotropy information to be obtained and is ideal for the investigation of molecular interactions occurring in anisotropic-oriented thin films. In this review, we will revisit main PWR applications, aiming at characterizing molecular interactions occurring (1) at lipid membranes deposited in the sensor and (2) in chemically modified sensors. Among the most widely used applications is the investigation of G-protein coupled receptor (GPCR) ligand activation and the study of the lipid environment’s impact on this process. Pioneering PWR studies on GPCRs were carried out thanks to the strong and effective collaboration between two laboratories in the University of Arizona leaded by Dr. Gordon Tollin and Dr. Victor J. Hruby. This review provides an overview of the main applications of PWR and provides a historical perspective on the development of instruments since the first prototype and continuous technological improvements to ongoing and future developments, aiming at broadening the information obtained and expanding the application portfolio.     

Could the presence of sodium ion influence the accuracy and precision of the ligand-posing in the human A<Subscript>2A</Subscript> adenosine receptor orthosteric binding site using a molecular docking approach? Insights from <Emphasis Type="Italic">Dockbench</Emphasis>

The allosteric modulation of G protein-coupled receptors (GPCRs) by sodium ions has received considerable attention as crystal structures of several receptors, in their inactive conformation, show a Na⁺ ion bound to specific residues which, in the human A_2A adenosine receptor (hA_2A AR), are Ser91^3.39, Trp246^6.48, Asn280^7.45, and Asn284^7.49. A cluster of water molecules completes the coordination of the sodium ion in the putative allosteric site. It is absolutely consolidated that the progress made in the field of GPCRs structural determination has increased the adoption of docking-driven approaches for the identification or the optimization of novel potent and selective ligands. Despite the extensive use of docking protocols in virtual screening approaches, to date, almost any of these studies have been carried out without taking into account the presence of the sodium cation and its first solvation shell in the putative allosteric binding site. In this study, we have focused our attention on determining how the presence of sodium ion binding and additionally its first hydration sphere, in hA_2AAR could influence the ligand positioning accuracy during molecular docking simulations for most of the available resting and activated hA_2A AR crystal structures, using DockBench as a comparative benchmarking tool and implementing a new correlation coefficient (EM). This work provides indications on the evidence that the posing performance (accuracy and/or precision) of the docking protocols in reproducing the crystallographic poses of different hA_2A AR antagonists is generally increased in the presence of the sodium cation and its first solvation shell, in agreement with experimental observations. Consequently, the inclusion of sodium ion and its first solvation shell should be considered in order to facilitate the selection of new potential ligands in all molecular docking-based virtual screening protocols that aim to find novel GPCRs antagonists and inverse agonists.     

G protein betagamma subunits as targets for small molecule therapeutic development

G proteins mediate the action of G protein coupled receptors (GPCRs), a major target of current pharmaceuticals and a major target of interest in future drug development. Most pharmaceutical interest has been in the development of selective GPCR agonists and antagonists that activate or inhibit specific GPCRs. Some recent thinking has focused on the idea that some pathologies are the result of the actions of an array of GPCRs suggesting that targeting single receptors may have limited efficacy. Thus, targeting pathways common to multiple GPCRs that control critical pathways involved in disease has potential therapeutic relevance. G protein betagamma subunits released from some GPCRs upon receptor activation regulate a variety of downstream pathways to control various aspects of mammalian physiology. There is evidence from cell- based and animal models that excess Gbetagamma signaling can be detrimental and blocking Gbetagamma signaling has salutary effects in a number of pathological models. Gbetagamma regulates downstream pathways through modulation of enzymes that produce cellular second messengers or through regulation of ion channels by direct protein-protein interactions. Thus, blocking Gbetagamma functions requires development of small molecule agents that disrupt Gbetagamma protein interactions with downstream partners. Here we discuss evidence that small molecule targeting Gbetagamma could be of therapeutic value. The concept of disruption of protein-protein interactions by targeting a "hot spot" on Gbetagamma is delineated and the biochemical and virtual screening strategies for identification of small molecules that selectively target Gbetagamma functions are outlined. Evaluation of the effectiveness of virtual screening indicates that computational screening enhanced identification of true Gbetagamma binding molecules. However, further refinement of the approach could significantly improve the yield of Gbetagamma binding molecules from this screen that could result in multiple candidate leads for future drug development.