Conserved waters mediate structural and functional activation of family A (rhodopsin-like) G protein-coupled receptors

Functional Water Channels Within the TSH Receptor: A New Paradigm for TSH Action With Disease Implications

The thyroid-stimulating hormone receptor (TSHR) transmembrane domain (TMD) is found in the plasma membrane and consists of lipids and water molecules. To understand the role of TSHR-associated water molecules, we used molecular dynamic simulations of the TMD and identified a network of putative receptor-associated transmembrane water channels. This result was confirmed with extended simulations of the full-length TSHR with and without TSH ligand binding. While the transport time observed in the simulations via the TSHR protein was slower than via the lipid bilayer itself, we found that significantly more water traversed via the TSHR than via the lipid bilayer, which more than doubled with the binding of TSH. Using rat thyroid cells (FRTL-5) and a calcein fluorescence technique, we measured cell volumes after blockade of aquaporins 1 and 4, the major thyroid cell water transporters. TSH showed a dose-dependent ability to influence water transport, and similar effects were observed with stimulating TSHR autoantibodies. Small molecule TSHR agonists, which are allosteric activators of the TMD, also enhanced water transport, illustrating the role of the TMD in this phenomenon. Furthermore, the water channel pathway was also mapped across 2 activation motifs within the TSHR TMD, suggesting how water movement may influence activation of the receptor. In pathophysiological conditions such as hypothyroidism and hyperthyroidism where TSH concentrations are highly variable, this action of TSH may greatly influence water movement in thyroid cells and many other extrathyroidal sites where the TSHR is expressed, thus affecting normal cellular function.

A short story on how chromophore is hydrolyzed from rhodopsin for recycling

The photocycle of visual opsins is essential to maintain the light sensitivity of the retina. The early physical observations of the rhodopsin photocycle by Böll and Kühne in the 1870s inspired over a century's worth of investigations on rhodopsin biochemistry. A single photon isomerizes the Schiff-base linked 11-cis-retinylidene chromophore of rhodopsin, converting it to the all-trans agonist to elicit phototransduction through photoactivated rhodopsin (Rho*). Schiff base hydrolysis of the agonist is a key step in the photocycle, not only diminishing ongoing phototransduction but also allowing for entry and binding of fresh 11-cis chromophore to regenerate the rhodopsin pigment and maintain light sensitivity. Many challenges have been encountered in measuring the rate of this hydrolysis, but recent advancements have facilitated studies of the hydrolysis within the native membrane environment of rhodopsin. These techniques can now be applied to study hydrolysis of agonist in other opsin proteins that mediate phototransduction or chromophore turnover. In this review, we discuss the progress that has been made in characterizing the rhodopsin photocycle and the journey to characterize the hydrolysis of its all-trans-retinylidene agonist.

Structural basis for the allosteric modulation of rhodopsin by nanobody binding to its extracellular domain

Rhodopsin is a prototypical G protein-coupled receptor (GPCR) critical for vertebrate vision. Research on GPCR signaling states has been facilitated using llama-derived nanobodies (Nbs), some of which bind to the intracellular surface to allosterically modulate the receptor. Extracellularly binding allosteric nanobodies have also been investigated, but the structural basis for their activity has not been resolved to date. Here, we report a library of Nbs that bind to the extracellular surface of rhodopsin and allosterically modulate the thermodynamics of its activation process. Crystal structures of Nb2 in complex with native rhodopsin reveal a mechanism of allosteric modulation involving extracellular loop 2 and native glycans. Nb2 binding suppresses Schiff base deprotonation and hydrolysis and prevents intracellular outward movement of helices five and six - a universal activation event for GPCRs. Nb2 also mitigates protein misfolding in a disease-associated mutant rhodopsin. Our data show the power of nanobodies to modulate the photoactivation of rhodopsin and potentially serve as therapeutic agents for disease-associated rhodopsin misfolding.

Structural dynamics of chemokine receptors

Membrane proteins such as G protein-coupled receptors (GPCRs) are involved in awide range of physiological and pathological cellular processes. Binding of extracellular signals to GPCRs, including hormones, neurotransmitters, peptides and proteins, can activate intracellular signaling cascades via G protein interaction. Chemokine receptors are key GPCRs implicated in cancers, immune responses, cell migration and inflammation. Specifically, the CCR5 and CXCR4 chemokine receptors serve as important therapeutic targets against Human Immunodeficiency virus (HIV) entry into human cells. Maraviroc and Vicriviroc, two clinically used HIV entry inhibitors, are antagonists of the CCR5 receptor. These drugs block HIV entry, but ultimately resistance develops, due to emergence of viruses that can utilize the CXCR4 co-receptor. Unfortunately, development of chemokine receptor antagonists as selective drugs of HIV infection has been greatly hindered as their target orthosteric site is conserved among different receptor subtypes. Accordingly, it is important to understand the structural dynamics of these receptors to develop more effective therapeutics. In this chapter, we describe the latest advances in studies of these two key chemokine receptors with respect to their structures, dynamics and function.

Counterion at an atypical position: Investigating the mechanism of photoisomerization in jellyfish rhodopsin

The position of the counterion in animal rhodopsins plays a crucial role in maintaining visible light sensitivity and facilitating the photoisomerization of their retinal chromophore. The counterion displacement is thought to be closely related to the evolution of rhodopsins, with different positions found in invertebrates and vertebrates. Interestingly, box jellyfish rhodopsin (JelRh) acquired the counterion in transmembrane 2 (TM2) independently. This is a unique feature, as in most animal rhodopsins, the counterion is found in a different location. In this study, we used Fourier Transform Infrared spectroscopy to examine the structural changes that occur in the early photointermediate state of JelRh. We aimed to determine whether the photochemistry of JelRh is similar to that of other animal rhodopsins by comparing its spectra to those of vertebrate bovine rhodopsin (BovRh) and invertebrate squid rhodopsin (SquRh). We observed that the N-D stretching band of the retinal Schiff base was similar to that of BovRh, indicating the interaction between the Schiff base and the counterion is similar in both rhodopsins, despite their different counterion positions. Furthermore, we found that the chemical structure of the retinal in JelRh is similar to that in BovRh, including the changes in the hydrogen-out-of-plane band that indicates a retinal distortion. Overall, the protein conformational changes induced by the photoisomerization of JelRh yielded spectra that resemble an intermediate between BovRh and SquRh, suggesting a unique spectral property of JelRh, and making it the only animal rhodopsin with a counterion in TM2 and an ability to activate G_s protein.

Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor

The United States is experiencing the most profound and devastating opioid crisis in history, with the number of deaths involving opioids, including prescription and illegal opioids, continuing to climb over the past two decades. This severe public health issue is difficult to combat as opioids remain a crucial treatment for pain, and at the same time, they are also highly addictive. Opioids act on the opioid receptor, which in turn activates its downstream signaling pathway that eventually leads to an analgesic effect. Among the four types of opioid receptors, the µ subtype is primarily responsible for the analgesic cascade. This review describes available 3D structures of the µ opioid receptor in the protein data bank and provides structural insights for the binding of agonists and antagonists to the receptor. Comparative analysis on the atomic details of the binding site in these structures was conducted and distinct binding interactions for agonists, partial agonists, and antagonists were observed. The findings in this article deepen our understanding of the ligand binding activity and shed some light on the development of novel opioid analgesics which may improve the risk benefit balance of existing opioids.

Molecular basis of ligand selectivity for melatonin receptors

Sleep disorders in adults are related to adverse health effects such as reduced quality of life and increased mortality. About 30-40% of adults are suffering from different sleep disorders. The human melatonin receptors (MT1 and MT2) are family A G protein-coupled receptors that respond to the neurohormone melatonin MEL which regulates circadian rhythm and sleep. Many efforts have been made to develop drugs targeting melatonin receptors to treat insomnia, circadian rhythm disorders, and even cancer. However, designing subtype-selective melatonergic drugs remains challenging due to their high similarities in both sequences and structures. MEL (a function-selective compound with a bulky β-naphthyl group) behaves as an MT2-selective antagonist, whereas it is an agonist of MT1. Here, molecular dynamics simulations were used to investigate the ligand selectivity of MT receptors at the atomic level. We found that the binding conformation of MEL differs in different melatonin receptors. In MT1, the naphthalene ring of MEL forms a structure perpendicular to the membrane surface. In contrast, there is a 130° angle between the naphthalene ring of MEL and the membrane surface in MT2. Because of this conformational difference, the MEL leads to a constant water channel in MT1 which activates the receptor. However, MEL hinders the formation of continuous water channels, resulting in an inactive state of MT2. Furthermore, we found that A117^3.29 in MT2 is a crucial amino acid capable of hindering the conformational flip of the MEL molecule. These results, coupled with previous functional data, reveal that although MT1 and MT2 share highly similar orthosteric ligand-binding pockets, they also display distinctive features that could be used to design selective compounds. Our findings provide new insights into functionally selective melatonergic drug development for sleep disorders.

Key Residues in δ Opioid Receptor Allostery Explored by the Elastic Network Model and the Complex Network Model Combined with the Perturbation Method

Opioid receptors, a kind of G protein-coupled receptors (GPCRs), mainly mediate an analgesic response via allosterically transducing the signal of endogenous ligand binding in the extracellular domain to couple to effector proteins in the intracellular domain. The δ opioid receptor (DOP) is associated with emotional control besides pain control, which makes it an attractive therapeutic target. However, its allosteric mechanism and key residues responsible for the structural stability and signal communication are not completely clear. Here we utilize the Gaussian network model (GNM) and amino acid network (AAN) combined with perturbation methods to explore the issues. The constructed fcfGNM_MD, where the force constants are optimized with the inverse covariance estimation based on the correlated fluctuations from the available DOP molecular dynamics (MD) ensemble, shows a better performance than traditional GNM in reproducing residue fluctuations and cross-correlations and in capturing functionally low-frequency modes. Additionally, fcfGNM_MD can consider implicitly the environmental effects to some extent. The lowest mode can well divide DOP segments and identify the two sodium ion (important allosteric regulator) binding coordination shells, and from the fastest modes, the key residues important for structure stabilization are identified. Using fcfGNM_MD combined with a dynamic perturbation-response method, we explore the key residues related to the sodium ion binding. Interestingly, we identify not only the key residues in sodium ion binding shells but also the ones far away from the perturbation sites, which are involved in binding with DOP ligands, suggesting the possible long-range allosteric modulation of sodium binding for the ligand binding to DOP. Furthermore, utilizing the weighted AAN combined with attack perturbations, we identify the key residues for allosteric communication. This work helps strengthen the understanding of the allosteric communication mechanism in δ opioid receptor and can provide valuable information for drug design.

G Protein-Coupled Receptors Regulated by Membrane Potential

G protein-coupled receptors (GPCRs) are involved in a vast majority of signal transduction processes. Although they span the cell membrane, they have not been considered to be regulated by the membrane potential. Numerous studies over the last two decades have demonstrated that several GPCRs, including muscarinic, adrenergic, dopaminergic, and glutamatergic receptors, are voltage regulated. Following these observations, an effort was made to elucidate the molecular basis for this regulatory effect. In this review, we will describe the advances in understanding the voltage dependence of GPCRs, the suggested molecular mechanisms that underlie this phenomenon, and the possible physiological roles that it may play.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

Identification and functional assays of single-nucleotide variants of opsins genes in melanocytic tumors

Epidermal melanocytes sense solar light via the opsin-coupled signaling pathway which is involved in a range of biological functions, including regulating pigmentation, proliferation, apoptosis, and tumorigenesis. However, it remains unclear whether there are genetic variants within these opsins that affect opsin protein structure and function, and further melanocyte biological behaviors. Here, we examined single-nucleotide variants (SNVs) of five opsin (RGR, OPN1SW, OPN2, OPN4, and OPN5) genes in MM (malignant melanoma; n = 76) and MN (melanocytic nevi; n = 157), using next-generation sequencing. The effects of these pathogenic single-nucleotide variants (SNVs) on opsin structure and function were further investigated using molecular dynamics (MD) simulations, dynamic cross-correlation (DCC), and site-directed mutagenesis. In total, 107 SNV variants were identified. Of these variants, 14 nonsynonymous SNVs (nsSNVs) of opsin genes were detected, including three mutations in the RGR gene, three mutations in the OPN1SW gene, two mutations in the OPN2 gene, and six mutations in the OPN4 gene. The effect of these missense mutations on opsin function was then assessed using eight prediction tools to estimate the potential impact of an amino acid substitution. The impact of each nsSNV was investigated using MD simulations and DCC analysis. Furthermore, we performed in vitro fluorescence calcium imaging to assess the functional properties of nsSNV proteins using a site-directed mutagenesis method. Taken together, these results revealed that p.A103V (RGR), p.T167I (RGR), p.G141S (OPN1SW), p.R144C (OPN1SW), and p.S231F (OPN4) had more deleterious effects on protein structure and function among the 14 nsSNVs. Opsin gene alterations showed the low frequency of missense mutations in melanocytic tumors, and although rare, some mutations in these opsin genes disrupt the canonical function of opsin. Our findings provide new insight into the role of opsin variants in the loss of function.

Hydration-mediated G-protein-coupled receptor activation

Although G-protein–coupled receptors (GPCRs) control vast physiological pathways, their activation remains chemically and physically enigmatic. Our osmotic stress studies of the visual receptor rhodopsin have redefined the standard model of GPCR signaling by revealing the essential role of bulk water. We show results consistent with a large number of water molecules flooding the rhodopsin interior during activation to stabilize the effector binding conformation. These results suggest a model of GPCR activation in which the receptor becomes solvent-swollen upon formation of the active state. We thus demonstrate the mechanism whereby water acts as a powerful allosteric modulator of a pharmacologically important membrane protein family.

The Rhodopsin family of G-protein–coupled receptors (GPCRs) comprises the targets of nearly a third of all pharmaceuticals. Despite structural water present in GPCR X-ray structures, the physiological relevance of these solvent molecules to rhodopsin signaling remains unknown. Here, we show experimental results consistent with the idea that rhodopsin activation in lipid membranes is coupled to bulk water movements into the protein. To quantify hydration changes, we measured reversible shifting of the metarhodopsin equilibrium due to osmotic stress using an extensive series of polyethylene glycol (PEG) osmolytes. We discovered clear evidence that light activation entails a large influx of bulk water (∼80–100 molecules) into the protein, giving insight into GPCR activation mechanisms. Various size polymer osmolytes directly control rhodopsin activation, in which large solutes are excluded from rhodopsin and dehydrate the protein, favoring the inactive state. In contrast, small osmolytes initially forward shift the activation equilibrium until a quantifiable saturation point is reached, similar to gain-of-function protein mutations. For the limit of increasing osmolyte size, a universal response of rhodopsin to osmotic stress is observed, suggesting it adopts a dynamic, hydrated sponge-like state upon photoactivation. Our results demand a rethinking of the role of water dynamics in modulating various intermediates in the GPCR energy landscape. We propose that besides bound water, an influx of bulk water plays a necessary role in establishing the active GPCR conformation that mediates signaling.

Determination of key residues in MRGPRX2 to enhance pseudo-allergic reactions induced by fluoroquinolones

MAS-related G protein-coupled receptor X2 (MRGPRX2), expressed in human mast cells, is associated with drug-induced pseudo-allergic reactions. Dogs are highly sensitive to the anaphylactoid reactions induced by certain drugs including fluoroquinolones. Recently, dog MRGPRX2 was identified as a functional ortholog of human MRGPRX2, with dog MRGPRX2 being particularly sensitive to fluoroquinolones. The aim of this study was to determine key residues responsible for the enhanced activity of fluoroquinolone-induced histamine release associated with MRGPRX2. Firstly, a structure model of human and dog MRGPRX2 was built by homology modeling, and docking simulations with fluoroquinolones were conducted. This model indicated that E164 and D184, conserved between human and dog, are essential for the binding to fluoroquinolones. In contrast, F78 (dog: Y) and M109 (dog: W) are unconserved residues, to which the species difference in fluoroquinolone sensitivity is attributable. Intracellular calcium mobilisation assay with human MRGPRX2 mutants, in which residues at positions 78 and 109 were substituted to those of dog MRGPRX2, revealed that M109 and F78 of human MRGPRX2 are crucial residues for enhancing the fluoroquinolone-induced histamine release. In conclusion, these key residues have important clinical implications for revealing the mechanisms and predicting the risks of fluoroquinolone-mediated pseudo-allergic reactions in humans.

The Many Faces of G Protein-Coupled Receptor 143, an Atypical Intracellular Receptor

GPCRs transform extracellular stimuli into a physiological response by activating an intracellular signaling cascade initiated via binding to G proteins. Orphan G protein-coupled receptors (GPCRs) hold the potential to pave the way for development of new, innovative therapeutic strategies. In this review we will introduce G protein-coupled receptor 143 (GPR143), an enigmatic receptor in terms of classification within the GPCR superfamily and localization. GPR143 has not been assigned to any of the GPCR families due to the lack of common structural motifs. Hence we will describe the most important motifs of classes A and B and compare them to the protein sequence of GPR143. While a precise function for the receptor has yet to be determined, the protein is expressed abundantly in pigment producing cells. Many GPR143 mutations cause X-linked Ocular Albinism Type 1 (OA1, Nettleship-Falls OA), which results in hypopigmentation of the eyes and loss of visual acuity due to disrupted visual system development and function. In pigment cells of the skin, loss of functional GPR143 results in abnormally large melanosomes (organelles in which pigment is produced). Studies have shown that the receptor is localized internally, including at the melanosomal membrane, where it may function to regulate melanosome size and/or facilitate protein trafficking to the melanosome through the endolysosomal system. Numerous additional roles have been proposed for GPR143 in determining cancer predisposition, regulation of blood pressure, development of macular degeneration and signaling in the brain, which we will briefly describe as well as potential ligands that have been identified. Furthermore, GPR143 is a promiscuous receptor that has been shown to interact with multiple other melanosomal proteins and GPCRs, which strongly suggests that this orphan receptor is likely involved in many different physiological actions.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

Inline Liquid Chromatography-Fast Photochemical Oxidation of Proteins for Targeted Structural Analysis of Conformationally Heterogeneous Mixtures

Structural analysis of proteins in a conformationally heterogeneous mixture has long been a difficult problem in structural biology. In structural analysis by covalent labeling mass spectrometry, conformational heterogeneity results in data reflecting a weighted average of all conformers, complicating data analysis and potentially causing misinterpretation of results. Here, we describe a method coupling size-exclusion chromatography (SEC) with hydroxyl radical protein footprinting using inline fast photochemical oxidation of proteins (FPOP). Using a controlled synthetic mixture of holomyoglobin and apomyoglobin, we validate that we can achieve accurate footprints of each conformer using LC-FPOP when compared to offline FPOP of each pure conformer. We then applied LC-FPOP to analyze the adalimumab heat-shock aggregation process. We found that the LC-FPOP footprint of unaggregated adalimumab was consistent with a previously published footprint of the native IgG. The LC-FPOP footprint of the aggregation product indicated that heat-shock aggregation primarily protected the hinge region, suggesting that this region is involved with the heat-shock aggregation process of this molecule. LC-FPOP offers a new method to probe dynamic conformationally heterogeneous mixtures that can be separated by SEC such as biopharmaceutical aggregates and to obtain accurate information on the topography of each conformer.

Allosteric Molecular Switches in Metabotropic Glutamate Receptors

Metabotropic glutamate receptors (mGlu) are class C G protein-coupled receptors of eight subtypes that are omnipresently expressed in the central nervous system. mGlus have relevance in several psychiatric and neurological disorders, therefore they raise considerable interest as drug targets. Allosteric modulators of mGlus offer advantages over orthosteric ligands owing to their increased potential to achieve subtype selectivity, and this has prompted discovery programs that have produced a large number of reported allosteric mGlu ligands. However, the optimization of allosteric ligands into drug candidates has proved to be challenging owing to induced-fit effects, flat or steep structure-activity relationships and unexpected changes in theirpharmacology. Subtle structural changes identified as molecular switches might modulate the functional activity of allosteric ligands. Here we review these switches discovered in the metabotropic glutamate receptor family..

Membrane dynamics of γ-secretase with the anterior pharynx-defective 1B subunit

The four-subunit protease complex γ-secretase cleaves many single-pass transmembrane (TM) substrates, including Notch and β-amyloid precursor protein to generate amyloid-β (Aβ), central to Alzheimer's disease. Two of the subunits anterior pharynx-defective 1 (APH-1) and presenilin (PS) exist in two homologous forms APH1-A and APH1-B, and PS1 and PS2. The consequences of these variations are poorly understood and could affect Aβ production and γ-secretase medicine. Here, we developed the first complete structural model of the APH-1B subunit using the published cryo-electron microscopy (cryo-EM) structures of APH1-A (Protein Data Bank: 5FN2, 5A63, and 6IYC). We then performed all-atom molecular dynamics simulations at 303 K in a realistic bilayer system to understand both APH-1B alone and in γ-secretase without and with substrate C83-bound. We show that APH-1B adopts a 7TM topology with a water channel topology similar to APH-1A. We demonstrate direct transport of water through this channel, mainly via Glu84, Arg87, His170, and His196. The apo and holo states closely resemble the experimental cryo-EM structures with APH-1A, however with subtle differences: The substrate-bound APH-1B γ-secretase was quite stable, but some TM helices of PS1 and APH-1B rearranged in the membrane consistent with the disorder seen in the cryo-EM data. This produces different accessibility of water molecules for the catalytic aspartates of PS1, critical for Aβ production. In particular, we find that the typical distance between the catalytic aspartates of PS1 and the C83 cleavage sites are shorter in APH-1B, that is, it represents a more closed state, due to interactions with the C-terminal fragment of PS1. Our structural-dynamic model of APH-1B alone and in γ-secretase suggests generally similar topology but some notable differences in water accessibility which may be relevant to the protein's existence in two forms and their specific function and location.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

HomolWat: a web server tool to incorporate 'homologous' water molecules into GPCR structures

Internal water molecules play an essential role in the structure and function of membrane proteins including G protein-coupled receptors (GPCRs). However, technical limitations severely influence the number and certainty of observed water molecules in 3D structures. This may compromise the accuracy of further structural studies such as docking calculations or molecular dynamics simulations. Here we present HomolWat, a web application for incorporating water molecules into GPCR structures by using template-based modelling of homologous water molecules obtained from high-resolution structures. While there are various tools available to predict the positions of internal waters using energy-based methods, the approach of borrowing lacking water molecules from homologous GPCR structures makes HomolWat unique. The tool can incorporate water molecules into a protein structure in about a minute with around 85% of water recovery. The web server is freely available at http://lmc.uab.es/homolwat.

Mass Spectrometry-Based Protein Footprinting for Higher-Order Structure Analysis: Fundamentals and Applications

Proteins adopt different higher-order structures (HOS) to enable their unique biological functions. Understanding the complexities of protein higher-order structures and dynamics requires integrated approaches, where mass spectrometry (MS) is now positioned to play a key role. One of those approaches is protein footprinting. Although the initial demonstration of footprinting was for the HOS determination of protein/nucleic acid binding, the concept was later adapted to MS-based protein HOS analysis, through which different covalent labeling approaches "mark" the solvent accessible surface area (SASA) of proteins to reflect protein HOS. Hydrogen-deuterium exchange (HDX), where deuterium in D2O replaces hydrogen of the backbone amides, is the most common example of footprinting. Its advantage is that the footprint reflects SASA and hydrogen bonding, whereas one drawback is the labeling is reversible. Another example of footprinting is slow irreversible labeling of functional groups on amino acid side chains by targeted reagents with high specificity, probing structural changes at selected sites. A third footprinting approach is by reactions with fast, irreversible labeling species that are highly reactive and footprint broadly several amino acid residue side chains on the time scale of submilliseconds. All of these covalent labeling approaches combine to constitute a problem-solving toolbox that enables mass spectrometry as a valuable tool for HOS elucidation. As there has been a growing need for MS-based protein footprinting in both academia and industry owing to its high throughput capability, prompt availability, and high spatial resolution, we present a summary of the history, descriptions, principles, mechanisms, and applications of these covalent labeling approaches. Moreover, their applications are highlighted according to the biological questions they can answer. This review is intended as a tutorial for MS-based protein HOS elucidation and as a reference for investigators seeking a MS-based tool to address structural questions in protein science.

Light-Sensitive Membrane Proteins as Tools to Generate Precision Treatments

INTRODUCTION BY ANA-NICOLETA BONDAR, BIOPHYSICS SECTION HEAD EDITOR: This issue of the Journal of Membrane Biology inaugurates Up-and-Coming Scientist, in which investigators at early career stages are invited to present recent research in the broad context of their discipline. We inaugurate Up-and-Coming Scientist with the essay by Dr. Elena Lesca of the ETH Zürich and the Paul Scherrer Institut, Switzerland. Dr. Lesca has completed her doctoral degree at the Technical University München, Germany, in 2014, and pursued postdoctoral research at the ETH Zürich and Paul Scherrer Institut, where she is Senior Assistant since 2019. Two recent papers by Dr. Lesca et al. (references 33 and 39) have used X-ray crystallography and experimental biophysics approaches to shed light on the mechanism of action of a membrane receptor from the G Protein-Coupled Receptor (GPCR) family, Jumping Spider Rhodopsin-1 (JSR-1). JSR-1 is a visual rhodopsin activated upon absorption of light by its covalently bound retinal chromophore. Unlike the better-understood bovine rhodopsin GPCR, which is monostable, JSR-1 is bistable (i.e., in JSR-1 the Schiff base that binds retinal to the protein stays protonated throughout the reaction cycle), and absorption of a second photon resets the retinal ligand to the resting state configuration. In her essay, Dr. Lesca discusses the implications of her work on JSR-1 and, more broadly, GPCR research, for state-of-the-art applications in optogenetics and drug design.

Preferential Coupling of Dopamine D2S and D2L Receptor Isoforms with Gi1 and Gi2 Proteins-In Silico Study

The dopamine D₂ receptor belongs to rhodopsin-like G protein-coupled receptors (GPCRs) and it is an important molecular target for the treatment of many disorders, including schizophrenia and Parkinson's disease. Here, computational methods were used to construct the full models of the dopamine D₂ receptor short (D_2S) and long (D_2L) isoforms (differing with 29 amino acids insertion in the third intracellular loop, ICL3) and to study their coupling with G_i1 and G_i2 proteins. It was found that the D_2L isoform preferentially couples with the G_i2 protein and D_2S isoform with the G_i1 protein, which is in accordance with experimental data. Our findings give mechanistic insight into the interplay between isoforms of dopamine D₂ receptors and G_i proteins subtypes, which is important to understand signaling by these receptors and their mediation by pharmaceuticals, in particular psychotic and antipsychotic agents.

Development of Container Free Sample Exposure for Synchrotron X-ray Footprinting

The method of X-ray footprinting and mass spectrometry (XFMS) on large protein assemblies and membrane protein samples requires high flux density to overcome the hydroxyl radical scavenging reactions produced by the buffer constituents and the total protein content. Previously, we successfully developed microsecond XFMS using microfluidic capillary flow and a microfocused broadband X-ray source at the Advanced Light Source synchrotron beamlines, but the excessive radiation damage incurred when using capillaries prevented the full usage of a high-flux density beam. Here we present another significant advance for the XFMS method: the instrumentation of a liquid injection jet to deliver container free samples to the X-ray beam. Our preliminary experiments with a liquid jet at a bending magnet X-ray beamline demonstrate the feasibility of the approach and show a significant improvement in the effective dose for both the Alexa fluorescence assay and protein samples compared to conventional capillary flow methods. The combination of precisely controlled high dose delivery, shorter exposure times, and elimination of radiation damage due to capillary effects significantly increases the signal quality of the hydroxyl radical modification products and the dose-response data. This new approach is the first application of container free sample handling for XFMS and opens up the method for even further advances, such as high-quality microsecond time-resolved XFMS studies.

QTY code designed thermostable and water-soluble chimeric chemokine receptors with tunable ligand affinity

Chemokine receptors are of great interest as they play a critical role in many immunological and pathological processes. The ability to study chemokine receptors in aqueous solution without detergent would be significant because natural receptors require detergents to become soluble. We previously reported using the QTY code to design detergent-free chemokine receptors. We here report the design of 2 detergent-free chimeric chemokine receptors that were experimentally unattainable in detergent solution. We designed chimeric receptors by switching the N terminus and 3 extracellular (EC) loops between different receptors. Specifically, we replaced the N terminus and 3 EC loops of CCR5^QTY with the N terminus and 3 EC loops of CXCR4. The ligand for CXCR4; namely CXCL12, binds to the chimeric receptor CCR5^QTY (7TM)-CXCR4 (N terminus+3 EC loops), but with lower affinity compared to CXCR4; the CCL5 ligand of CCR5 binds the chimeric receptor with ∼20× lower affinity. The chimeric design helps to elucidate the mechanism of native receptor-ligand interaction. We also show that all detergent-free QTY-designed chemokine receptors, expressed in Escherichia coli, bind to their respective chemokines with affinities in the nanomolar (nM) range, similar to the affinities of native receptors and SF9-produced QTY variants. These QTY-designed receptors exhibit remarkable thermostability in the presence of arginine and retain ligand-binding activity after heat treatment at 60 °C for 4 h and 24 h, and at 100 °C for 10 min. Our design approach enables affordable scale-up production of detergent-free QTY variant chemokine receptors with tunable functionality for various uses.

Crystal structure of jumping spider rhodopsin-1 as a light sensitive GPCR

Light-sensitive G protein-coupled receptors (GPCRs)-rhodopsins-absorb photons to isomerize their covalently bound retinal, triggering conformational changes that result in downstream signaling cascades. Monostable rhodopsins release retinal upon isomerization as opposed to the retinal in bistable rhodopsins that "reisomerize" upon absorption of a second photon. Understanding the mechanistic differences between these light-sensitive GPCRs has been hindered by the scarcity of recombinant models of the latter. Here, we reveal the high-resolution crystal structure of a recombinant bistable rhodopsin, jumping spider rhodopsin-1, bound to the inverse agonist 9-cis retinal. We observe a water-mediated network around the ligand hinting toward the basis of their bistable nature. In contrast to bovine rhodopsin (monostable), the transmembrane bundle of jumping spider rhodopsin-1 as well that of the bistable squid rhodopsin adopts a more "activation-ready" conformation often observed in other nonphotosensitive class A GPCRs. These similarities suggest the role of jumping spider rhodopsin-1 as a potential model system in the study of the structure-function relationship of both photosensitive and nonphotosensitive class A GPCRs.

Chemokine Receptor Crystal Structures: What Can Be Learned from Them?

Chemokine receptors belong to the class A of G protein-coupled receptors (GPCRs) and are implicated in a wide variety of physiologic functions, mostly related to the homeostasis of the immune system. Chemokine receptors are also involved in multiple pathologic processes, including immune and autoimmune diseases, as well as cancer. Hence, several members of this GPCR subfamily are considered to be very relevant therapeutic targets. Since drug discovery efforts can be significantly reinforced by the availability of crystal structures, substantial efforts in the area of chemokine receptor structural biology could dramatically increase the outcome of drug discovery campaigns. This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications. SIGNIFICANCE STATEMENT: This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications.

APH-1A Component of γ-Secretase Forms an Internal Water and Ion-Containing Cavity

Anterior pharynx-defective 1A (APH-1A) is a seven transmembrane component of γ-secretase (GS), an aspartyl protease enzyme involved in the production of toxic amyloid-β peptides in Alzheimer's disease patients. Cryo-electron microscopy structures of the enzyme complex revealed a central cavity in its APH-1A component, similar to water-containing cavities in G-protein coupled receptors (GPCRs). In this work, we performed molecular dynamics and umbrella sampling simulations to understand the role of the APH-1A cavity in the GS complex. Our results suggest that APH-1A is able to store water molecules in its inner cavity and transport some of them between cell spaces. Additionally, APH-1A allows the influx of extracellular cations into a central hydrophilic cavity but cannot transport them into the intracellular space. Overall, this study seeks to describe an alternative APH-1A function in GS besides its complex stabilization role and provide novel approaches to understand the functioning of the GS enzyme.

The role of water molecules in phototransduction of retinal proteins and G protein-coupled receptors

G protein coupled receptors (GPCRs) are a key family of membrane proteins in all eukaryotes and also very important drug targets for medical intervention. The extensively studied visual pigment rhodopsin is a prime example of a family A GPCR. Its chromophore ligand retinal is covalently linked to a lysine in helix seven forming a protonated Schiff base. Interestingly, this is the same situation in other-non-GPCR-retinal proteins, like the prototype light-driven microbial proton pump bacteriorhodopsin, albeit there is no (or only a very remote) phylogenetical link. Close to the retinal ligand, several water molecules help to organise a functionally important hydrogen bond network that undergoes significant changes during photo-activation. Such water-mediated networks are also critical for ligand binding of other GPCRs and they are becoming increasingly important in drug discovery. GPCRs also contain a partially conserved water mediated hydrogen bond network that stabilises the ground state of the receptor, and rearrangement of this network leads to the stabilization of the active state. Some water molecules have a specific role in this process to appropriately orient specific residues relative to the Schiff base, and to modulate the fine structure of the transmembrane bundle, particularly near the intracellular G protein binding site. While the atomic details of these mechanisms are still missing, the recent developments in free electron lasers (FELs) are enabling us to begin to observe the changes in waters and relevant side chains shortly after photo activation at an unprecedented level of spatial and temporal resolution.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

Structural basis for ligand modulation of the CCR2 conformational landscape

CC chemokine receptor 2 (CCR2) is a part of the chemokine receptor family, an important class of therapeutic targets. These class A G-protein coupled receptors (GPCRs) are involved in mammalian signaling pathways and control cell migration toward endogenous CC chemokine ligands, named for the adjacent cysteine motif on their N terminus. Chemokine receptors and their associated ligands are involved in a wide range of diseases and thus have become important drug targets. CCR2, in particular, promotes the metastasis of cancer cells and is also implicated in autoimmunity-driven type-1 diabetes, diabetic nephropathy, multiple sclerosis, asthma, atherosclerosis, neuropathic pain, and rheumatoid arthritis. Although promising, CCR2 antagonists have been largely unsuccessful to date. Here, we investigate the effect of an orthosteric and an allosteric antagonist on CCR2 dynamics by coupling long-timescale molecular dynamics simulations with Markov-state model theory. We find that the antagonists shift CCR2 into several stable inactive conformations that are distinct from the crystal structure conformation and disrupt a continuous internal water and sodium ion pathway, preventing transitions to an active-like state. Several metastable conformations present a cryptic drug-binding pocket near the allosteric site that may be amenable to targeting with small molecules. Without antagonists, the apo dynamics reveal intermediate conformations along the activation pathway that provide insight into the basal dynamics of CCR2 and may also be useful for future drug design.

Open-Boundary Molecular Mechanics/Coarse-Grained Framework for Simulations of Low-Resolution G-Protein-Coupled Receptor-Ligand Complexes

G-protein-coupled receptors (GPCRs) constitute as much as 30% of the overall proteins targeted by FDA-approved drugs. However, paucity of structural experimental information and low sequence identity between members of the family impair the reliability of traditional docking approaches and atomistic molecular dynamics simulations for in silico pharmacological applications. We present here a dual-resolution approach tailored for such low-resolution models. It couples a hybrid molecular mechanics/coarse-grained (MM/CG) scheme, previously developed by us for GPCR-ligand complexes, with a Hamiltonian-based adaptive resolution scheme (H-AdResS) for the solvent. This dual-resolution approach removes potentially inaccurate atomistic details from the model while building a rigorous statistical ensemble-the grand canonical one-in the high-resolution region. We validate the method on a well-studied GPCR-ligand complex, for which the 3D structure is known, against atomistic simulations. This implementation paves the way for future accurate in silico studies of low-resolution ligand/GPCRs models.

Using X-ray Footprinting and Mass Spectrometry to Study the Structure and Function of Membrane Proteins

BACKGROUND

Membrane proteins are crucial for cellular sensory cascades and metabolite transport, and hence are key pharmacological targets. Structural studies by traditional highresolution techniques are limited by the requirements for high purity and stability when handled in high concentration and nonnative buffers. Hence, there is a growing requirement for the use of alternate methods in a complementary but orthogonal approach to study the dynamic and functional aspects of membrane proteins in physiologically relevant conditions. In recent years, significant progress has been made in the field of X-ray radiolytic labeling in combination with mass spectroscopy, commonly known as X-ray Footprinting and Mass Spectrometry (XFMS), which provide residue-specific information on the solvent accessibility of proteins. In combination with both lowresolution biophysical methods and high-resolution structural data, XFMS is capable of providing valuable insights into structure and dynamics of membrane proteins, which have been difficult to obtain by standalone high-resolution structural techniques. The XFMS method has also demonstrated a unique capability for identification of structural waters and their dynamics in protein cavities at both a high degree of spatial and temporal resolution, and thus capable of identifying conformational hot-spots in transmembrane proteins.

CONCLUSION

We provide a perspective on the place of XFMS amongst other structural biology methods and showcase some of the latest developments in its usage for studying conformational changes in membrane proteins.

Specificity of the chromophore-binding site in human cone opsins

The variable composition of the chromophore-binding pocket in visual receptors is essential for vision. The visual phototransduction starts with the cis-trans isomerization of the retinal chromophore upon absorption of photons. Despite sharing the common 11-cis-retinal chromophore, rod and cone photoreceptors possess distinct photochemical properties. Thus, a detailed molecular characterization of the chromophore-binding pocket of these receptors is critical to understanding the differences in the photochemistry of vision between rods and cones. Unlike for rhodopsin (Rh), the crystal structures of cone opsins remain to be determined. To obtain insights into the specific chromophore-protein interactions that govern spectral tuning in human visual pigments, here we harnessed the unique binding properties of 11-cis-6-membered-ring-retinal (11-cis-6mr-retinal) with human blue, green, and red cone opsins. To unravel the specificity of the chromophore-binding pocket of cone opsins, we applied 11-cis-6mr-retinal analog-binding analyses to human blue, green, and red cone opsins. Our results revealed that among the three cone opsins, only blue cone opsin can accommodate the 11-cis-6mr-retinal in its chromophore-binding pocket, resulting in the formation of a synthetic blue pigment (B6mr) that absorbs visible light. A combination of primary sequence alignment, molecular modeling, and mutagenesis experiments revealed the specific amino acid residue 6.48 (Tyr-262 in blue cone opsins and Trp-281 in green and red cone opsins) as a selectivity filter in human cone opsins. Altogether, the results of our study uncover the molecular basis underlying the binding selectivity of 11-cis-6mr-retinal to the cone opsins.

Diverse GPCRs exhibit conserved water networks for stabilization and activation

G protein-coupled receptors (GPCRs) represent both the largest class of drug targets and the largest family of human membrane proteins. Recent three-dimensional structures have revealed the presence of many water molecules buried inside GPCRs, but the functional role of these water molecules has been unclear. Using extensive atomic-level computer simulations, we find that, although most of these waters are highly mobile, a few are stable. These stable water molecules mediate state-dependent and state-independent polar networks that are conserved across diverse GPCRs. These findings promise to help guide the rational design of GPCR-targeted drugs.

G protein-coupled receptors (GPCRs) have evolved to recognize incredibly diverse extracellular ligands while sharing a common architecture and structurally conserved intracellular signaling partners. It remains unclear how binding of diverse ligands brings about GPCR activation, the common structural change that enables intracellular signaling. Here, we identify highly conserved networks of water-mediated interactions that play a central role in activation. Using atomic-level simulations of diverse GPCRs, we show that most of the water molecules in GPCR crystal structures are highly mobile. Several water molecules near the G protein-coupling interface, however, are stable. These water molecules form two kinds of polar networks that are conserved across diverse GPCRs: (i) a network that is maintained across the inactive and the active states and (ii) a network that rearranges upon activation. Comparative analysis of GPCR crystal structures independently confirms the striking conservation of water-mediated interaction networks. These conserved water-mediated interactions near the G protein-coupling region, along with diverse water-mediated interactions with extracellular ligands, have direct implications for structure-based drug design and GPCR engineering.

Glu2.53(90) of the GnRH receptor is part of the conserved G protein-coupled receptor structure and does not form a salt-bridge with Lys3.32(121)

GnRH receptor mutations, Glu^2.53(90)Lys and Glu^2.53(90)Asp, cause congenital hypogonadotropic hypogonadism. The Glu^2.53(90) side-chain has been proposed to form an intramolecular salt-bridge with Lys^3.32(121), but conserved intramolecular interaction networks in G protein-coupled receptor crystal structures predict that it interacts with Ser^3.35(124) and Trp^6.48(280). We investigated interhelical interactions of Glu^2.53(90) that stabilise GnRH receptor folding using functional analyses and computational modelling of mutant receptors. The Glu^2.53(90)Asp mutant was non-functional, but mutants with hydrophobic amino acids or Arg substituted for Glu^2.53(90) were functional, excluding a salt-bridge interaction. The Glu^2.53(90)Arg and Trp^6.48(280)Arg mutants had decreased affinity for GnRH. Models showed that congenital Glu^2.53(90)Lys and Glu^2.53(90)Asp mutations disrupt interactions with Ser^3.35(124) and Trp^6.48(280) respectively, whereas the Glu^2.53(90)Arg and Trp^6.48(280)Arg mutations preserve intramolecular contacts, but increase distance between the transmembrane helices. Our results show that disruption of interhelical contacts that are conserved in G protein-coupled receptors accounts for the effects of some disease-associated GnRH receptor mutations.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D^6.33 and K303E^7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D^6.33 and R307E^7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D^6.33/R307E^7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Functional characterisation of naturally occurring mutations in human melanopsin

Melanopsin is a blue light-sensitive opsin photopigment involved in a range of non-image forming behaviours, including circadian photoentrainment and the pupil light response. Many naturally occurring genetic variants exist within the human melanopsin gene (OPN4), yet it remains unclear how these variants affect melanopsin protein function and downstream physiological responses to light. Here, we have used bioinformatic analysis and in vitro expression systems to determine the functional phenotypes of missense human OPN4 variants. From 1242 human OPN4 variants collated in the NCBI Short Genetic Variation database (dbSNP), we identified 96 that lead to non-synonymous amino acid substitutions. These 96 missense mutations were screened using sequence alignment and comparative approaches to select 16 potentially deleterious variants for functional characterisation using calcium imaging of melanopsin-driven light responses in HEK293T cells. We identify several previously uncharacterised OPN4 mutations with altered functional properties, including attenuated or abolished light responses, as well as variants demonstrating abnormal response kinetics. These data provide valuable insight into the structure-function relationships of human melanopsin, including several key functional residues of the melanopsin protein. The identification of melanopsin variants with significantly altered function may serve to detect individuals with disrupted melanopsin-based light perception, and potentially highlight those at increased risk of sleep disturbance, circadian dysfunction, and visual abnormalities.

A Novel Polar Core and Weakly Fixed C-Tail in Squid Arrestin Provide New Insight into Interaction with Rhodopsin

Photoreceptors of the squid Loligo pealei contain a G-protein-coupled receptor (GPCR) signaling system that activates phospholipase C in response to light. Analogous to the mammalian visual system, signaling of the photoactivated GPCR rhodopsin is terminated by binding of squid arrestin (sArr). sArr forms a light-dependent, high-affinity complex with squid rhodopsin, which does not require prior receptor phosphorylation for interaction. This is at odds with classical mammalian GPCR desensitization where an agonist-bound phosphorylated receptor is needed to break stabilizing constraints within arrestins, the so-called "three-element interaction" and "polar core" network, before a stable receptor-arrestin complex can be established. Biophysical and mass spectrometric analysis of the squid rhodopsin-arrestin complex indicates that in contrast to mammalian arrestins, the sArr C-tail is not involved in a stable three-element interaction. We determined the crystal structure of C-terminally truncated sArr that adopts a basal conformation common to arrestins and is stabilized by a series of weak but novel polar core interactions. Unlike mammalian arrestin-1, deletion of the sArr C-tail does not influence kinetic properties of complex formation of sArr with the receptor. Hydrogen-deuterium exchange studies revealed the footprint of the light-activated rhodopsin on sArr. Furthermore, double electron-electron resonance spectroscopy experiments provide evidence that receptor-bound sArr adopts a conformation different from the one known for arrestin-1 and molecular dynamics simulations reveal the residues that account for the weak three-element interaction. Insights gleaned from studying this system add to our general understanding of GPCR-arrestin interaction.

Orthosteric and allosteric action of the C5a receptor antagonists

The C5a receptor (C5aR) is a G-protein-coupled receptor (GPCR) that can induce strong inflammatory response to the anaphylatoxin C5a. Targeting C5aR has emerged as a novel anti-inflammatory therapeutic method. However, developing potent C5aR antagonists as drugs has proven difficult. Here, we report two crystal structures of human C5aR in ternary complexes with the peptide antagonist PMX53 and a non-peptide antagonist, either avacopan or NDT9513727. The structures, together with other biophysical, computational docking and cell-based signaling data, reveal the orthosteric action of PMX53 and its effect of stabilizing the C5aR structure, as well as the allosteric action of chemically diverse non-peptide C5aR antagonists with different binding poses. Structural comparison analysis suggests the presence of similar allosteric sites in other GPCRs. We also discuss critical structural features of C5aR in activation, including a novel conformation of helix 8. On the basis of our results, we suggest novel strategies for developing C5aR-targeting drugs.

High-Resolution Hydroxyl Radical Protein Footprinting: Biophysics Tool for Drug Discovery

Hydroxyl radical footprinting (HRF) of proteins with mass spectrometry (MS) is a widespread approach for assessing protein structure. Hydroxyl radicals react with a wide variety of protein side chains, and the ease with which radicals can be generated (by radiolysis or photolysis) has made the approach popular with many laboratories. As some side chains are less reactive and thus cannot be probed, additional specific and nonspecific labeling reagents have been introduced to extend the approach. At the same time, advances in liquid chromatography and MS approaches permit an examination of the labeling of individual residues, transforming the approach to high resolution. Lastly, advances in understanding of the chemistry of the approach have led to the determination of absolute protein topologies from HRF data. Overall, the technology can provide precise and accurate measures of side-chain solvent accessibility in a wide range of interesting and useful contexts for the study of protein structure and dynamics in both academia and industry.

Serial Femtosecond Crystallography of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) represent a large superfamily of membrane proteins that mediate cell signaling and regulate a variety of physiological processes in the human body. Structure-function studies of this superfamily were enabled a decade ago by multiple breakthroughs in technology that included receptor stabilization, crystallization in a membrane environment, and microcrystallography. The recent emergence of X-ray free-electron lasers (XFELs) has further accelerated structural studies of GPCRs and other challenging proteins by overcoming radiation damage and providing access to high-resolution structures and dynamics using micrometer-sized crystals. Here, we summarize key technology advancements and major milestones of GPCR research using XFELs and provide a brief outlook on future developments in the field.

The 3D Structure of Human DP Prostaglandin G-Protein-Coupled Receptor Bound to Cyclopentanoindole Antagonist, Predicted Using the DuplexBiHelix Modification of the GEnSeMBLE Method

Prostaglandins play a critical physiological role in both cardiovascular and immune systems, acting through their interactions with 9 prostanoid G protein-coupled receptors (GPCRs). These receptors are important therapeutic targets for a variety of diseases including arthritis, allergies, type 2 diabetes, and cancer. The DP prostaglandin receptor is of interest because it has unique structural and physiological properties. Most notably, DP does not have the 3-6 ionic lock common to Class A GPCRs. However, the lack of X-ray structures for any of the 9 prostaglandin GPCRs hampers the application of structure-based drug design methods to develop more selective and active medications to specific receptors. We predict here 3D structures for the DP prostaglandin GPCR, based on the GEnSeMBLE complete sampling with hierarchical scoring (CS-HS) methodology. This involves evaluating the energy of 13 trillion packings to finally select the best 20 that are stable enough to be relevant for binding to antagonists, agonists, and modulators. To validate the predicted structures, we predict the binding site for the Merck cyclopentanoindole (CPI) selective antagonist docked to DP. We find that the CPI binds vertically in the 1-2-7 binding pocket, interacting favorably with residues R310^7.40 and K76^2.54 with additional interactions with S313^7.43, S316^7.46, S19^1.35, etc. This binding site differs significantly from that of antagonists to known Class A GPCRs where the ligand binds in the 3-4-5-6 region. We find that the predicted binding site leads to reasonable agreement with experimental Structure-Activity Relationship (SAR). We suggest additional mutation experiments including K76^2.54, E129^3.49, L123^3.43, M270^6.40, F274^6.44 to further validate the structure, function, and activation mechanism of receptors in the prostaglandin family. Our structures and binding sites are largely consistent and improve upon the predictions by Li et al. ( J. Am. Chem. Soc. 2007 , 129 ( 35 ), 10720 ) that used our earlier MembStruk prediction methodology.