Search

The orphan G protein-coupled receptor (GPCR) GPR161 is enriched in primary cilia, where it plays a central role in suppressing Hedgehog signaling1. GPR161 mutations lead to developmental defects and cancers2,3,4. The fundamental basis of how GPR161 is activated, including potential endogenous activators and pathway-relevant signal transducers, remains unclear. To elucidate GPR161 function, we determined a cryogenic-electron microscopy structure of active GPR161 bound to the heterotrimeric G protein complex Gs. This structure revealed an extracellular loop 2 that occupies the canonical GPCR orthosteric ligand pocket. Furthermore, we identify a sterol that binds to a conserved extrahelical site adjacent to transmembrane helices 6 and 7 and stabilizes a GPR161 conformation required for Gscoupling. Mutations that prevent sterol binding to GPR161 suppress cAMP pathway activation. Surprisingly, these mutants retain the ability to suppress GLI2 transcription factor accumulation in cilia, a key function of ciliary GPR161 in Hedgehog pathway suppression. By contrast, a protein kinase A-binding site in the GPR161 C-terminus is critical in suppressing GLI2 ciliary accumulation. Our work highlights how unique structural features of GPR161 interface with the Hedgehog pathway and sets a foundation to understand the broader role of GPR161 function in other signaling pathways.

The neuropeptide substance P (SP) is important in pain and inflammation. SP activates the neurokinin-1 receptor (NK1R) to signal via Gq and Gs proteins. Neurokinin A also activates NK1R, but leads to selective Gq signaling. How two stimuli yield distinct G protein signaling at the same G protein-coupled receptor remains unclear. We determined cryogenic-electron microscopy structures of active NK1R bound to SP or the Gq-biased peptide SP6-11. Peptide interactions deep within NK1R are critical for receptor activation. Conversely, interactions between SP and NK1R extracellular loops are required for potent Gs signaling but not Gq signaling. Molecular dynamics simulations showed that these superficial contacts restrict SP flexibility. SP6-11, which lacks these interactions, is dynamic while bound to NK1R. Structural dynamics of NK1R agonists therefore depend on interactions with the receptor extracellular loops and regulate G protein signaling selectivity. Similar interactions between other neuropeptides and their cognate receptors may tune intracellular signaling.

Thyroid hormones are vital in metabolism, growth and development(1). Thyroid hormone synthesis is controlled by thyrotropin (TSH), which acts at the thyrotropin receptor (TSHR)(2). In patients with Graves’ disease, autoantibodies that activate the TSHR pathologically increase thyroid hormone activity(3). How autoantibodies mimic thyrotropin function remains unclear. Here we determined cryo-electron microscopy structures of active and inactive TSHR. In inactive TSHR, the extracellular domain lies close to the membrane bilayer. Thyrotropin selects an upright orientation of the extracellular domain owing to steric clashes between a conserved hormone glycan and the membrane bilayer. An activating autoantibody from a patient with Graves’ disease selects a similar upright orientation of the extracellular domain. Reorientation of the extracellular domain transduces a conformational change in the seven-transmembrane-segment domain via a conserved hinge domain, a tethered peptide agonist and a phospholipid that binds within the seven-transmembrane-segment domain. Rotation of the TSHR extracellular domain relative to the membrane bilayer is sufficient for receptor activation, revealing a shared mechanism for other glycoprotein hormone receptors that may also extend to other G-protein-coupled receptors with large extracellular domains.

Choosing the drug to fit the protein Many approved drugs bind to G protein–coupled receptors (GPCRs). A challenge in targeting GPCRs is that different ligands preferentially activate different signaling pathways. Two papers show how biased signaling arises for the angiotensin II type 1 receptor that couples to two signaling partners (G proteins and arrestins). Suomivuori et al. used large-scale atomistic simulations to show that coupling to the two pathways is through two distinct GPCR conformations and that extracellular ligands favor one or the other conformation. Wingler et al. present crystal structures of the same receptor bound to ligands with different bias profiles. These structures show conformational changes in and around the binding pocket that match those observed in simulations. This work could provide a framework for the rational design of drugs that are more effective and have fewer side effects. Science , this issue p. 881 , p. 888

Human cytomegalovirus (HCMV) encodes G protein–coupled receptors (GPCRs) US28 and US27, which facilitate viral pathogenesis through engagement of host G proteins. Here we report cryo–electron microscopy structures of US28 and US27 forming nonproductive and productive complexes with Gi and Gq, respectively, exhibiting unusual features with functional implications. The “orphan” GPCR US27 lacks a ligand-binding pocket and has captured a guanosine diphosphate–bound inactive Gi through a tenuous interaction. The docking modes of CX3CL1-US28 and US27 to Gi favor localization to endosome-like curved membranes, where US28 and US27 can function as nonproductive Gi sinks to attenuate host chemokine-dependent Gi signaling. The CX3CL1-US28-Gq/11 complex likely represents a trapped intermediate during productive signaling, providing a view of a transition state in GPCR–G protein coupling for signaling. Our collective results shed new insight into unique G protein–mediated HCMV GPCR structural mechanisms, compared to mammalian GPCR counterparts, for subversion of host immunity., Description, The HCMV GPCR-human G protein structures reveal unusual host-pathogen interactions likely critical for HCMV’s pathogenesis.

Drugs targeting the μ-opioid receptor (μOR) are the most effective analgesics available but are also associated with fatal respiratory depression through a pathway that remains unclear. Here we investigated the mechanistic basis of action of lofentanil (LFT) and mitragynine pseudoindoxyl (MP), two μOR agonists with different safety profiles. LFT, one of the most lethal opioids, and MP, a kratom plant derivative with reduced respiratory depression in animal studies, exhibited markedly different efficacy profiles for G protein subtype activation and β-arrestin recruitment. Cryo-EM structures of μOR-Gi1 complex with MP (2.5 Å) and LFT (3.2 Å) revealed that the two ligands engage distinct subpockets, and molecular dynamics simulations showed additional differences in the binding site that promote distinct active-state conformations on the intracellular side of the receptor where G proteins and β-arrestins bind. These observations highlight how drugs engaging different parts of the μOR orthosteric pocket can lead to distinct signaling outcomes.

Neuronal spiking activity is routinely recorded using genetically encoded calcium indicators (GECIs), but calcium imaging is limited in temporal resolution and does not report subthreshold voltage changes. Genetically encoded voltage indicators (GEVIs) offer better temporal resolution and subthreshold sensitivity, but spike detection with fast GEVIs has required specialized imaging equipment. Here, we report the ASAP4 subfamily of genetically encoded voltage indicators (GEVIs) that brighten in response to membrane depolarization, inverting the fluorescence-voltage relationship of previous ASAP GEVIs. Two variants, ASAP4b and ASAP4e, feature 128% and 178% fluorescence increases over 100-mV of depolarization, respectively, facilitating spike detection in single trials in vivo with standard 1 and 2-photon imaging systems. Simultaneous voltage and calcium imaging confirms improved temporal resolution and spike discernment by ASAP4 GEVIs. Thus, positively tuned ASAP4 voltage indicators enable recording of neuronal spiking activity using similar equipment as calcium imaging, while providing higher temporal resolution.One Sentence SummaryUpward ASAPs increase detection capability of GEVIs in vivo.

Serotonin (5-hydroxytryptamine [5-HT]) 5-HT2-family receptors represent essential targets for lysergic acid diethylamide (LSD) and all other psychedelic drugs. Although the primary psychedelic drug effects are mediated by the 5-HT

The neuropeptide Substance P (SP) is important in pain and inflammation. SP activates the neurokinin-1 receptor (NK1R) to signal via Gqand Gsproteins. Neurokinin A also activates NK1R, but leads to selective Gqsignaling. How two stimuli yield distinct G-protein signaling at the same G-protein-coupled-receptor remains unclear. We determined cryo-EM structures of active NK1R bound to SP or the Gq-biased peptide SP6-11. Peptide interactions deep within NK1R are critical for receptor activation. Conversely, interactions between SP and NK1R extracellular loops are required for potent Gssignaling but not Gqsignaling. Molecular dynamics simulations showed that these superficial contacts restrict SP flexibility deep in the NK1R pocket. SP6-11, which lacks these interactions, is dynamic while bound to NK1R. Structural dynamics of NK1R agonists therefore depend on interactions with the receptor extracellular loops and regulate G-protein signaling selectivity. Similar interactions between other neuropeptides and their cognate receptors may tune intracellular signaling.

Choosing the drug to fit the protein Many approved drugs bind to G protein–coupled receptors (GPCRs). A challenge in targeting GPCRs is that different ligands preferentially activate different signaling pathways. Two papers show how biased signaling arises for the angiotensin II type 1 receptor that couples to two signaling partners (G proteins and arrestins). Suomivuori et al. used large-scale atomistic simulations to show that coupling to the two pathways is through two distinct GPCR conformations and that extracellular ligands favor one or the other conformation. Wingler et al. present crystal structures of the same receptor bound to ligands with different bias profiles. These structures show conformational changes in and around the binding pocket that match those observed in simulations. This work could provide a framework for the rational design of drugs that are more effective and have fewer side effects. Science , this issue p. 881 , p. 888

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Rational design of light-capturing properties requires understanding the molecular and electronic structure of chromophores in their native chemical or biological environment. We employ here large-scale quantum chemical calculations to study the light-capturing properties of retinal in recently designed human cellular retinol binding protein II (hCRBPII) variants (Wang et al. Science, 2012, 338, 1340-1343). Our calculations show that these proteins absorb across a large part of the visible spectrum by combined polarization and electrostatic effects. These effects stabilize the ground or excited state energy levels of the retinal by perturbing the Schiff-base or β-ionone moieties of the chromophore, which in turn modulates the amount of charge transfer within the molecule. Based on the predicted tuning principles, we design putative in silico mutations that further shift the absorption properties of retinal in hCRBPII towards the ultraviolet and infrared regions of the spectrum.

Biased signaling, in which different ligands that bind to the same G protein-coupled receptor preferentially trigger distinct signaling pathways, holds great promise for the design of safer and more effective drugs. Its structural mechanism remains unclear, however, hampering efforts to design drugs with desired signaling profiles. Here, we use extensive atomic-level molecular dynamics simulations to determine how arrestin bias and G protein bias arise at the angiotensin II type 1 receptor. The receptor adopts two major signaling conformations, one of which couples almost exclusively to arrestin, whereas the other also couples effectively to a G protein. A long-range allosteric network allows ligands in the extracellular binding pocket to favor either of the two intracellular conformations. Guided by this computationally determined mechanism, we designed ligands with desired signaling profiles.

The light-harvesting chlorophyll (Chl) molecules of photosynthetic systems form the basis for light-driven energy conversion. In biological environments, the Chl chromophores occur in two distinct diastereotopic configurations, where the alpha and beta configurations have a magnesium-ligating histidine residue and a 17-propionic acid moiety on the opposite side or on the same side of the Chl ring, respectively. Although beta-ligated Chl dimers occupy conserved positions around the reaction center of photosystem I (PSI), the functional relevance of the alpha/beta configuration of the ligation is poorly understood. We employ here correlated ab initio calculations using the algebraic-diagrammatic construction through second order (ADC(2)) and the approximate second-order coupled cluster (CC2) methods in combination with the reduced virtual space (RVS) approach in studies of the intrinsic excited-state properties of alpha-ligated and beta-ligated Chl dimers of PSI. Our ab initio calculations suggest that the absorption of the alpha-ligated reaction-center Chl dimer of PSI is redshifted by 0.13-0.14 eV in comparison to the beta-ligated dimers due to combined excitonic coupling and strain effects. We also show that time-dependent density functional theory (TDDFT) calculations using range-separated density functionals underestimate the absorption shift between the alpha- and beta-ligated dimers. Our findings may provide a molecular starting point for understanding the energy flow in natural photosynthetic systems, as well as a blueprint for developing new molecules that convert sunlight into other forms of energy.

Accurate predictions of photoexcitation properties are a major challenge for modern methods of theoretical chemistry. We show here how approximate coupled-cluster singles and doubles (CC2) calculations in combination with the reduced virtual space (RVS) approach can be employed in studies of excited states of large biomolecular systems. The RVS-CC2 approach is used for accurately predicting optical properties of the p-hydroxybenzylidene-dihydroimidazolinone (p-HBDI) chromophore embedded in green fluorescent protein (GFP) models using quantum mechanical calculations in combination with large basis sets. We study the lowest excited states for the isolated and protein-embedded chromophore in two different protonation states, and show how omitting high-lying virtual orbitals in the RVS calculation of excitation energies renders large-scale CC2 studies computationally feasible. We also discuss how the error introduced by the RVS approach can be systematically estimated and controlled. The obtained CC2 excitation energies of 3.13-3.27 and 2.69-2.77 eV for the two protonation states of different protein models are in excellent agreement with the maxima of the experimental absorption spectra of 3.12-3.14 and 2.61-2.64 eV, respectively. Thus, the calculated energy splitting between the excited states of the two protonation states is 0.44-0.52 eV, which agrees very well with the experimental value of 0.48-0.51 eV. The calculations at the RVS-CC2 level on the protein models show the importance of using large QM regions in studies of biochromophores embedded in proteins.

The conversion of light energy into ion gradients across biological membranes is one of the most fundamental reactions in primary biological energy transduction. Recently, the structure of the first light-activated Na+ pump, Krokinobacter eikastus rhodopsin 2 (KR2), was resolved at atomic resolution [Kato HE, et al. (2015) Nature 521: 48-53]. To elucidate its molecular mechanism for Na+ pumping, we perform here extensive classical and quantum molecular dynamics (MD) simulations of transient photocycle states. Our simulations show how the dynamics of key residues regulate water and ion access between the bulk and the buried light-triggered retinal site. We identify putative Na+ binding sites and show how protonation and conformational changes gate the ion through these sites toward the extracellular side. We further show by correlated ab initio quantum chemical calculations that the obtained putative photocycle intermediates are in close agreement with experimental transient optical spectroscopic data. The combined results of the ion translocation and gating mechanisms in KR2 may provide a basis for the rational design of novel light-driven ion pumps with optogenetic applications.

Chlorophylls are light-capturing units found in photosynthetic proteins. We study here the ground and excited state properties of monomeric, dimeric, and tetrameric models of the special chlorophyll/bacteriochlorophyll (Chl/BChl) pigment (P) centers P700 and P680/P870 of type I and type II photosystems, respectively. In the excited state calculations, we study the performance of the algebraic diagrammatic construction through second-order (ADC(2)) method in combination with the reduced virtual space (RVS) approach and the recently developed Laplace-transformed scaled-opposite-spin (LT-SOS) algorithm, which allows us, for the first time, to address multimeric effects at correlated ab initio levels using large basis sets. At the LT-SOS-RVS-ADC(2)/def2-TZVP level, we obtain vertical excitation energies (VEEs) of 2.00-2.07 and 1.52-1.62 eV for the P680/P700 and the P870 pigment models, respectively, which agree well with the experimental absorption maxima of 1.82, 1.77, and 1.43 eV for P680, P700, and P870, respectively. In the P680/P870 models, we find that the photoexcitation leads to a π → π* transition in which the exciton is delocalized between the adjacent Chl/BChl molecules of the central pair, whereas the exciton is localized to a single chlorophyll molecule in the P700 model. Consistent with experiments, the calculated excitonic splittings between the central pairs of P680, P700, and P870 models are 80, 200, and 400 cm(-1), respectively. The calculations show that the electron affinity of the radical cation of the P680 model is 0.4 V larger than for the P870 model and 0.2 V larger than for P700. The chromophore stacking interaction is found to strongly influence the electron localization properties of the light-absorbing pigments, which may help to elucidate mechanistic details of the charge separation process in type I and type II photosystems.

Structural changes in the metal site of the copper-zinc superoxide dismutase (SOD1) are involved in the various mechanisms proposed for the pathogenesis of the SOD1-linked familial form of amyotrophic lateral sclerosis (ALS). Elucidating how the metal site of SOD1 can be disrupted by ALS-linked mutations is important for a better understanding of the pathogenesis of the disease and for developing more efficient treatments. Residue D124, a second-sphere ligand of the copper and zinc ions, is known from experimental studies to be essential for the integrity of the metal-site structure. In this work, we used density functional theory calculations and molecular dynamics simulations to elucidate which factors keep D124 attached to the metal site and how structural changes may disrupt the binding between D124 and the metal first-sphere ligands. The calculations show that D124 is kept attached to the metal site in a kinetic trap. The exclusion of solvent molecules by the electrostatic loop of the protein is found to create the binding of D124 to the metal site. The calculations also indicate that changes in the structure of the electrostatic loop of the protein can weaken the D124-metal site interaction, lowering the affinity of the zinc site for the metal. Destabilization of the electrostatic loop of SOD1 has been previously shown to be a common property of ALS-linked variants of the protein, but its role in the pathogenesis of SOD1-linked ALS has not been elucidated.