Molecular basis for high-affinity agonist binding in GPCRs

GPCR Activation States Induced by Nanobodies and Mini-G Proteins Compared by NMR Spectroscopy

In this work, we examine methyl nuclear magnetic resonance (NMR) spectra of the methionine ε-[¹³CH₃] labelled thermostabilized β₁ adrenergic receptor from turkey in association with a variety of different effectors, including mini-G_s and nanobody 60 (Nb60), which have not been previously studied in complex with β₁ adrenergic receptor (β₁AR) by NMR. Complexes with pindolol and Nb60 induce highly similar inactive states of the receptor, closely resembling the resting state conformational ensemble. We show that, upon binding of mini-Gs or nanobody 80 (Nb80), large allosteric changes throughout the receptor take place. The conformation of tβ₁AR stabilized by the native-like mini-G_s protein is highly similar to the conformation induced by the currently used surrogate Nb80. Interestingly, in both cases residual dynamics are present, which were not observed in the resting states. Finally, we reproduce a pharmaceutically relevant situation, where an antagonist abolishes the interaction of the receptor with the mini-G protein in a competitive manner, validating the functional integrity of our preparation. The presented system is therefore well suited for reproducing the individual steps of the activation cycle of a G protein-coupled receptor (GPCR) in vitro and serves as a basis for functional and pharmacological characterizations of more native-like systems in the future.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

Structure of the class D GPCR Ste2 dimer coupled to two G proteins

G-protein-coupled receptors (GPCRs) are divided phylogenetically into six classes^1,2, denoted A to F. More than 370 structures of vertebrate GPCRs (belonging to classes A, B, C and F) have been determined, leading to a substantial understanding of their function³. By contrast, there are no structures of class D GPCRs, which are found exclusively in fungi where they regulate survival and reproduction. Here we determine the structure of a class D GPCR, the Saccharomyces cerevisiae pheromone receptor Ste2, in an active state coupled to the heterotrimeric G protein Gpa1-Ste4-Ste18. Ste2 was purified as a homodimer coupled to two G proteins. The dimer interface of Ste2 is formed by the N terminus, the transmembrane helices H1, H2 and H7, and the first extracellular loop ECL1. We establish a class D1 generic residue numbering system (CD1) to enable comparisons with orthologues and with other GPCR classes. The structure of Ste2 bears similarities in overall topology to class A GPCRs, but the transmembrane helix H4 is shifted by more than 20 Å and the G-protein-binding site is a shallow groove rather than a cleft. The structure provides a template for the design of novel drugs to target fungal GPCRs, which could be used to treat numerous intractable fungal diseases⁴.

Different conformational responses of the β2-adrenergic receptor-Gs complex upon binding of the partial agonist salbutamol or the full agonist isoprenaline

G protein-coupled receptors (GPCRs) are responsible for most cytoplasmic signaling in response to extracellular ligands with different efficacy profiles. Various spectroscopic techniques have identified that agonists exhibiting varying efficacies can selectively stabilize a specific conformation of the receptor. However, the structural basis for activation of the GPCR-G protein complex by ligands with different efficacies is incompletely understood. To better understand the structural basis underlying the mechanisms by which ligands with varying efficacies differentially regulate the conformations of receptors and G proteins, we determined the structures of β2AR-Gαs$\beta $γ bound with partial agonist salbutamol or bound with full agonist isoprenaline using single-particle cryo-electron microscopy at resolutions of 3.26 Å and 3.80 Å, respectively. Structural comparisons between the β2AR-Gs-salbutamol and β2AR-Gs-isoprenaline complexes demonstrated that the decreased binding affinity and efficacy of salbutamol compared with those of isoprenaline might be attributed to weakened hydrogen bonding interactions, attenuated hydrophobic interactions in the orthosteric binding pocket and different conformational changes in the rotamer toggle switch in TM6. Moreover, the observed stronger interactions between the intracellular loop 2 or 3 (ICL2 or ICL3) of β2AR and Gαs with binding of salbutamol versus isoprenaline might decrease phosphorylation in the salbutamol-activated β2AR-Gs complex. From the observed structural differences between these complexes of β2AR, a mechanism of β2AR activation by partial and full agonists is proposed to provide structural insights into β2AR desensitization.

Synthesis and dopamine receptor pharmacological evaluations on ring C ortho halogenated 1-phenylbenzazepines

A series of 1-phenylbenzazepines containing bromine or chlorine substituents at the ortho position of the appended phenyl ring (2'-monosubstituted or 2',6'- disubstituted patterns) were synthesized and evaluated for affinity towards dopamine D₁R, D₂R and D₅R. As is typical of the 1-phenylbenzazepine scaffold, the compounds displayed selectivity towards D₁R and D₅R; analogs generally lacked affinity for D₂R. Interestingly, 2',6'-dichloro substituted analogs showed modest D₅R versus D₁R selectivity whereas this selectivity was reversed in compounds with a 2'-halo substitution pattern. Compound 10a was identified as a D₁R antagonist (K_i = 14 nM; IC₅₀ = 9.4 nM).

Allosteric Mechanisms of Nonadditive Substituent Contributions to Protein-Ligand Binding

Quantifying chemical substituent contributions to ligand-binding free energies is challenging due to nonadditive effects. Protein allostery is a frequent cause of nonadditivity, but the underlying allosteric mechanisms often remain elusive. Here, we propose a general NMR-based approach to elucidate such mechanisms and we apply it to the HCN4 ion channel, whose cAMP-binding domain is an archetypal conformational switch. Using NMR, we show that nonadditivity arises not only from concerted conformational transitions, but also from conformer-specific effects, such as steric frustration. Our results explain how affinity-reducing functional groups may lead to affinity gains if combined. Surprisingly, our approach also reveals that nonadditivity depends markedly on the receptor conformation. It is negligible for the inhibited state but highly significant for the active state, opening new opportunities to tune potency and agonism of allosteric effectors.

Crystal structure of the human oxytocin receptor

The peptide hormone oxytocin modulates socioemotional behavior and sexual reproduction via the centrally expressed oxytocin receptor (OTR) across several species. Here, we report the crystal structure of human OTR in complex with retosiban, a nonpeptidic antagonist developed as an oral drug for the prevention of preterm labor. Our structure reveals insights into the detailed interactions between the G protein-coupled receptor (GPCR) and an OTR-selective antagonist. The observation of an extrahelical cholesterol molecule, binding in an unexpected location between helices IV and V, provides a structural rationale for its allosteric effect and critical influence on OTR function. Furthermore, our structure in combination with experimental data allows the identification of a conserved neurohypophyseal receptor-specific coordination site for Mg²⁺ that acts as potent, positive allosteric modulator for agonist binding. Together, these results further our molecular understanding of the oxytocin/vasopressin receptor family and will facilitate structure-guided development of new therapeutics.

Molecular basis of β-arrestin coupling to formoterol-bound β1-adrenoceptor

The β₁-adrenoceptor (β₁AR) is a G-protein-coupled receptor (GPCR) that couples¹ to the heterotrimeric G protein G_s. G-protein-mediated signalling is terminated by phosphorylation of the C terminus of the receptor by GPCR kinases (GRKs) and by coupling of β-arrestin 1 (βarr1, also known as arrestin 2), which displaces G_s and induces signalling through the MAP kinase pathway². The ability of synthetic agonists to induce signalling preferentially through either G proteins or arrestins-known as biased agonism³-is important in drug development, because the therapeutic effect may arise from only one signalling cascade, whereas the other pathway may mediate undesirable side effects⁴. To understand the molecular basis for arrestin coupling, here we determined the cryo-electron microscopy structure of the β₁AR-βarr1 complex in lipid nanodiscs bound to the biased agonist formoterol⁵, and the crystal structure of formoterol-bound β₁AR coupled to the G-protein-mimetic nanobody⁶ Nb80. βarr1 couples to β₁AR in a manner distinct to that⁷ of G_s coupling to β₂AR-the finger loop of βarr1 occupies a narrower cleft on the intracellular surface, and is closer to transmembrane helix H7 of the receptor when compared with the C-terminal α5 helix of G_s. The conformation of the finger loop in βarr1 is different from that adopted by the finger loop of visual arrestin when it couples to rhodopsin⁸. β₁AR coupled to βarr1 shows considerable differences in structure compared with β₁AR coupled to Nb80, including an inward movement of extracellular loop 3 and the cytoplasmic ends of H5 and H6. We observe weakened interactions between formoterol and two serine residues in H5 at the orthosteric binding site of β₁AR, and find that formoterol has a lower affinity for the β₁AR-βarr1 complex than for the β₁AR-G_s complex. The structural differences between these complexes of β₁AR provide a foundation for the design of small molecules that could bias signalling in the β-adrenoceptors.

Psi4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

A high-resolution description of β1-adrenergic receptor functional dynamics and allosteric coupling from backbone NMR

Signal transmission and regulation of G-protein-coupled receptors (GPCRs) by extra- and intracellular ligands occurs via modulation of complex conformational equilibria, but their exact kinetic details and underlying atomic mechanisms are unknown. Here we quantified these dynamic equilibria in the β₁-adrenergic receptor in its apo form and seven ligand complexes using ¹H/¹⁵N NMR spectroscopy. We observe three major exchanging conformations: an inactive conformation (C_i), a preactive conformation (C_p) and an active conformation (C_a), which becomes fully populated in a ternary complex with a G protein mimicking nanobody. The C_i ↔ C_p exchange occurs on the microsecond scale, the C_p ↔ C_a exchange is slower than ~5 ms and only occurs in the presence of two highly conserved tyrosines (Y^5.58, Y^7.53), which stabilize the active conformation of TM6. The C_p→C_a chemical shift changes indicate a pivoting motion of the entire TM6 that couples the effector site to the orthosteric ligand pocket.

Signal transmission and regulation of G-protein-coupled receptors (GPCRs) by ligands occurs via modulation of complex conformational equilibria. Here authors quantify these equilibria and their dynamics in the β1-adrenergic receptor in its apo form and seven ligand complexes using NMR.

Truncated (N)-Methanocarba Nucleosides as Partial Agonists at Mouse and Human A3 Adenosine Receptors: Affinity Enhancement by N6-(2-Phenylethyl) Substitution

Dopamine-derived N6-substituents, compared to N6-(2-phenylethyl), in truncated (N)-methanocarba (bicyclo[3.1.0]hexyl) adenosines favored high A3 adenosine receptor (AR) affinity/selectivity, e.g., C2-phenylethynyl analogue 15 (MRS7591, Ki = 10.9/17.8 nM, at human/mouse A3AR). 15 was a partial agonist in vitro (hA3AR, cAMP inhibition, 31% Emax; mA3AR, [35S]GTP-γ-S binding, 16% Emax) and in vivo and also antagonized hA3AR in vitro. Distal H-bonding substitutions of the N6-(2-phenylethyl) moiety particularly enhanced mA3AR affinity by polar interactions with the extracellular loops, predicted using docking and molecular dynamics simulation with newly constructed mA3AR and hA3AR homology models. These hybrid models were based on an inactive antagonist-bound hA1AR structure for the upper part of TM2 and an agonist-bound hA2AAR structure for the remaining TM portions. These species-independent A3AR-selective nucleosides are low efficacy partial agonists and novel, nuanced modulators of the A3AR, a drug target of growing interest.

The atomistic level structure for the activated human κ-opioid receptor bound to the full Gi protein and the MP1104 agonist

The kappa opioid receptor (κOR) is an important target for pain therapeutics to reduce depression and other harmful side effects of existing medications. The analgesic activity is mediated by κOR signaling through the adenylyl cyclase-inhibitory family of Gi protein. Here, we report the three-dimensional (3D) structure for the active state of human κOR complexed with both heterotrimeric Gi protein and MP1104 agonist. This structure resulted from long molecular dynamics (MD) and metadynamics (metaMD) simulations starting from the 3.1-Å X-ray structure of κOR-MP1104 after replacing the nanobody with the activated Gi protein and from the 3.5-Å cryo-EM structure of μOR-Gi complex after replacing the 168 missing residues. Using MD and metaMD we discovered interactions to the Gi protein with strong anchors to two intracellular loops and transmembrane helix 6 of the κOR. These anchors strengthen the binding, contributing to a contraction in the binding pocket but an expansion in the cytoplasmic region of κOR to accommodate G protein. These remarkable changes in κOR structure reveal that the anchors are essential for activation.

Haloperidol bound D2 dopamine receptor structure inspired the discovery of subtype selective ligands

The D₂ dopamine receptor (DRD2) is one of the most well-established therapeutic targets for neuropsychiatric and endocrine disorders. Most clinically approved and investigational drugs that target this receptor are known to be subfamily-selective for all three D₂-like receptors, rather than subtype-selective for only DRD2. Here, we report the crystal structure of DRD2 bound to the most commonly used antipsychotic drug, haloperidol. The structures suggest an extended binding pocket for DRD2 that distinguishes it from other D₂-like subtypes. A detailed analysis of the structures illuminates key structural determinants essential for DRD2 activation and subtype selectivity. A structure-based and mechanism-driven screening combined with a lead optimization approach yield DRD2 highly selective agonists, which could be used as chemical probes for studying the physiological and pathological functions of DRD2 as well as promising therapeutic leads devoid of promiscuity.

The D₂ dopamine receptor (DRD2) is one of the most well-established therapeutic targets for neuropsychiatric and endocrine disorders. Here, the authors report the crystal structure of the antipsychotic drug haloperidol bound to DRD2 via an extended binding pocket that distinguishes it from other D₂-like subtypes.

Conformational plasticity of ligand-bound and ternary GPCR complexes studied by 19F NMR of the β1-adrenergic receptor

G-protein-coupled receptors (GPCRs) are allosteric signaling proteins that transmit an extracellular stimulus across the cell membrane. Using ¹⁹F NMR and site-specific labelling, we investigate the response of the cytoplasmic region of transmembrane helices 6 and 7 of the β₁-adrenergic receptor to agonist stimulation and coupling to a G_s-protein-mimetic nanobody. Agonist binding shows the receptor in equilibrium between two inactive states and a pre-active form, increasingly populated with higher ligand efficacy. Nanobody coupling leads to a fully active ternary receptor complex present in amounts correlating directly with agonist efficacy, consistent with partial agonism. While for different agonists the helix 6 environment in the active-state ternary complexes resides in a well-defined conformation, showing little conformational mobility, the environment of the highly conserved NPxxY motif on helix 7 remains dynamic adopting diverse, agonist-specific conformations, implying a further role of this region in receptor function. An inactive nanobody-coupled ternary receptor form is also observed.

Structural insights into melatonin receptors

The long-anticipated high-resolution structures of the human melatonin G protein-coupled receptors MT₁ and MT₂ , involved in establishing and maintaining circadian rhythm, were obtained in complex with two melatonin analogs and two approved anti-insomnia and antidepression drugs using X-ray free-electron laser serial femtosecond crystallography. The structures shed light on the overall conformation and unusual structural features of melatonin receptors, as well as their ligand binding sites and the melatonergic pharmacophore, thereby providing insights into receptor subtype selectivity. The structures revealed an occluded orthosteric ligand binding site with a membrane-buried channel for ligand entry in both receptors, and an additional putative ligand entry path in MT₂ from the extracellular side. This unexpected ligand entry mode contributes to facilitating the high specificity with which melatonin receptors bind their cognate ligand and exclude structurally similar molecules such as serotonin, the biosynthetic precursor of melatonin. Finally, the MT₂ structure allowed accurate mapping of type 2 diabetes-related single-nucleotide polymorphisms, where a clustering of residues in helices I and II on the protein-membrane interface was observed which could potentially influence receptor oligomerization. The role of receptor oligomerization is further discussed in light of the differential interaction of MT₁ and MT₂ with GPR50, a regulatory melatonin coreceptor. The melatonin receptor structures will facilitate design of selective tool compounds to further dissect the specific physiological function of each receptor subtype as well as provide a structural basis for next-generation sleeping aids and other drugs targeting these receptors with higher specificity and fewer side effects.

Azi-medetomidine: Synthesis and Characterization of a Novel α2 Adrenergic Photoaffinity Ligand

Agonists at the α₂ adrenergic receptor produce sedation, increase focus, provide analgesia, and induce centrally mediated hypotension and bradycardia, yet neither their dynamic interactions with adrenergic receptors nor their modulation of neuronal circuit activity is completely understood. Photoaffinity ligands of α₂ adrenergic agonists have the potential both to capture discrete moments of ligand-receptor interactions and to prolong naturalistic drug effects in discrete regions of tissue in vivo. We present here the synthesis and characterization of a novel α₂ adrenergic agonist photolabel based on the imidazole medetomidine called azi-medetomidine. Azi-medetomidine shares protein association characteristics with its parent compound in experimental model systems and by molecular dynamics simulation of interactions with the α_2A adrenergic receptor. Azi-medetomidine acts as an agonist at α_2A adrenergic receptors, and produces hypnosis in Xenopus laevis tadpoles. Azi-medetomidine competes with the α₂ agonist clonidine at α_2A adrenergic receptors, which is potentiated by photolabeling, and azi-medetomidine labels moieties on the α_2A adrenergic receptor as determined by mass spectrometry in a manner consistent with a simulated model. This novel α₂ adrenergic agonist photolabel can serve as a powerful tool for in vitro and in vivo investigations of adrenergic signaling.

A Selectivity Study of FFAR4/FFAR1 Agonists by Molecular Modeling

FFAR4 has been considered as a potential target for metabolic diseases, including diabetes. Some compounds with biphenyl scaffold, represented by compound SR13 reported by our group, showed significant FFAR4 selectivity. However, the molecular basis for their selectivity has not been definitely disclosed. This study provided insights into the protein-ligand interactions between agonists and FFAR4/FFAR1 by molecular modeling. The important residues identified were consistent with those found in experimental studies. Moreover, the results proposed that the selectivity of SR13 between FFAR4 and FFAR1 depended on whether it can enter the ligand-binding site through the entrance region by adopting its preferential conformation. The big difference between the preferential conformation of SR13 and the narrow entrance region determined its poor agonist activity against FFAR1. These findings will facilitate the further development of selective FFAR4 agonists.

High Pressure Shifts the β1-Adrenergic Receptor to the Active Conformation in the Absence of G Protein

G protein-coupled receptors (GPCRs) are versatile chemical sensors, which transmit the signal of an extracellular binding event across the plasma membrane to the intracellular side. This function is achieved via the modulation of highly dynamical equilibria of various conformational receptor states. Here we have probed the effect of pressure on the conformational equilibria of a functional thermostabilized β₁-adrenergic GPCR (β₁AR) by solution NMR. High pressure induces a large shift in the conformational equilibrium (midpoint ∼600 bar) from the preactive conformation of agonist-bound β₁AR to the fully active conformation, which under normal pressure is only populated when a G protein or a G protein-mimicking nanobody (Nb) binds to the intracellular side of the β₁AR·agonist complex. No such large effects are observed for an antagonist-bound β₁AR or the ternary β₁AR·agonist·Nb80 complex. The detected structural changes of agonist-bound β₁AR around the orthosteric ligand binding pocket indicate that the fully active receptor occupies an ∼100 ų smaller volume than that of its preactive form. Most likely, this volume reduction is caused by the compression of empty (nonhydrated) cavities in the ligand binding pocket and the center of the receptor, which increases the ligand receptor interactions and explains the ∼100-fold affinity increase of agonists in the presence of G protein. The finding that isotropic pressure induces a directed motion from the preactive to the fully active GPCR conformation provides evidence of the high mechanical robustness of this important functional switch.