The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

pH-dependent absorption spectra of rhodopsin mutant E113Q: On the role of counterions and protein

The absorption spectra of bovine rhodopsin mutant E113Q in solutions were investigated at the molecular level by using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The calculations suggest the mechanism of the absorption variations of E113Q at different pH values. The results indicate that the polarizations of the counterions in the vicinity of Schiff base under protonation and unprotonation states of the mutant E113Q would be a crucial factor to change the energy gap of the retinal to tune the absorption spectra. Glu-181 residue, which is close to the chromophore, cannot serve as the counterion of the protonated Schiff base of E113Q in dark state. Moreover, the results of the absorption maximum in mutant E113Q with the various anions (Cl^-, Br^-, I^- and NO₃^-) manifested that the mutant E113Q could have the potential for use as a template of anion biosensors at visible wavelength.

From Homology Models to a Set of Predictive Binding Pockets-a 5-HT1A Receptor Case Study

Despite its remarkable importance in the arena of drug design, serotonin 1A receptor (5-HT1A) has been elusive to the X-ray crystallography community. This lack of direct structural information not only hampers our knowledge regarding the binding modes of many popular ligands (including the endogenous neurotransmitter-serotonin), but also limits the search for more potent compounds. In this paper we shed new light on the 3D pharmacological properties of the 5-HT1A receptor by using a ligand-guided approach (ALiBERO) grounded in the Internal Coordinate Mechanics (ICM) docking platform. Starting from a homology template and set of known actives, the method introduces receptor flexibility via Normal Mode Analysis and Monte Carlo sampling, to generate a subset of pockets that display enriched discrimination of actives from inactives in retrospective docking. Here, we thoroughly investigated the repercussions of using different protein templates and the effect of compound selection on screening performance. Finally, the best resulting protein models were applied prospectively in a large virtual screening campaign, in which two new active compounds were identified that were chemically distinct from those described in the literature.

7×7 RMSD matrix: A new method for quantitative comparison of the transmembrane domain structures in the G-protein coupled receptors

The G-protein coupled receptors (GPCRs) share a conserved heptahelical fold in the transmembrane (TM) region, but the exact arrangements of the seven TM helices vary with receptors and their activation states. The differences or the changes have been observed in the experimentally solved structures, but have not been systematically and quantitatively investigated due to lack of suitable methods. In this work, we describe a novel method, called 7×7 RMSD matrix that is proposed specifically for comparing the characteristic 7TM bundle structures of GPCRs. Compared to the commonly used overall TM bundle RMSD as a single parameter, a 7×7 RMSD matrix contains 49 parameters, which reveal changes of the relative orientations of the seven TMs. We demonstrate the novelty and advantages of this method by tackling two problems that are challenging for the existing methods. With this method, we are able to identify and quantify the helix movements in the activated receptor structures and reveal structural conservation and divergence as well as the structural relationships of different GPCRs in terms of the relative orientations of the seven TMs.

Point charge embedding for ONIOM excited states calculations

Hybrid quantum mechanical methods can assist in the interpretation and prediction of the electronic spectra of large molecular structures. In this work, we study the performance of the ONIOM (Our own N-layered Integrated molecular Orbital molecular Mechanics) hybrid method for the calculation of transition energies and oscillator strengths by embedding the core region in a field of fixed point charges. These charges introduce polarization effects from the substituent groups to the core region. We test various charge definitions, with particular attention to the issue of overpolarization near the boundary between layers. To minimize this issue, we fit the charges on the electrostatic potential of the entire structure in the presence of the link atoms used to cap dangling bonds. We propose two constrained fitting strategies: one that produces an average set of charges common to both model system calculations, EE(L1), and one that produces two separate sets of embedding charges, EE(L2). The results from our tests show that indeed electronic embedding with constrained-fitted charges tends to improve the performance of ONIOM compared to non-embedded calculations. However, the EE(L2) charges work best for transition energies, and the EE(L1) charges work best for oscillator strengths. This may be an indication that fixed point charges do not have enough flexibility to adapt to each system, and other effects (e.g., polarization of the embedding field) may be necessary.

Model-Based Discovery of Synthetic Agonists for the Zn2+-Sensing G-Protein-Coupled Receptor 39 (GPR39) Reveals Novel Biological Functions

The G-protein-coupled receptor 39 (GPR39) is a G-protein-coupled receptor activated by Zn²⁺. We used a homology model-based approach to identify small-molecule pharmacological tool compounds for the receptor. The method focused on a putative binding site in GPR39 for synthetic ligands and knowledge of ligand binding to other receptors with similar binding pockets to select iterative series of minilibraries. These libraries were cherry-picked from all commercially available synthetic compounds. A total of only 520 compounds were tested in vitro, making this method broadly applicable for tool compound development. The compounds of the initial library were inactive when tested alone, but lead compounds were identified using Zn²⁺ as an allosteric enhancer. Highly selective, highly potent Zn²⁺-independent GPR39 agonists were found in subsequent minilibraries. These agonists identified GPR39 as a novel regulator of gastric somatostatin secretion.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Molecular-Docking-Based Drug Design and Discovery

Currently 30-50% of drug targets are G Protein-Coupled Receptors (GPCRs). However, the clinical useful drugs for targeting GPCR have been limited by the lack of subtype selectivity or efficacy, leading to undesirable side effects. To develop subtype-selective GPCR ligands with desired molecular properties, better understanding is needed of the pharmacophore elements and of the binding mechanism required for subtype selectivity. To illustrate these issues, we describe here three successful applications to understand the binding mechanism associated with subtype selectivity: 5-HT2B (5-Hydroxytryptamine, 5-HT) serotonin receptor (HT2BR), H3 histamine receptor (H3HR) and A3 adenosine receptor (A3AR). The understanding of structure-function relationships among individual types and subtypes of GPCRs gained from such computational predictions combined with experimental validation and testing is expected the development of new highly selective and effective ligands to address such diseases while minimizing side-effects.

Modulation of thermal noise and spectral sensitivity in Lake Baikal cottoid fish rhodopsins

Lake Baikal is the deepest and one of the most ancient lakes in the world. Its unique ecology has resulted in the colonization of a diversity of depth habitats by a unique fauna that includes a group of teleost fish of the sub-order Cottoidei. This relatively recent radiation of cottoid fishes shows a gradual blue-shift in the wavelength of the absorption maximum of their visual pigments with increasing habitat depth. Here we combine homology modeling and quantum chemical calculations with experimental in vitro measurements of rhodopsins to investigate dim-light adaptation. The calculations, which were able to reproduce the trend of observed absorption maxima in both A1 and A2 rhodopsins, reveal a Barlow-type relationship between the absorption maxima and the thermal isomerization rate suggesting a link between the observed blue-shift and a thermal noise decrease. A Nakanishi point-charge analysis of the electrostatic effects of non-conserved and conserved amino acid residues surrounding the rhodopsin chromophore identified both close and distant sites affecting simultaneously spectral tuning and visual sensitivity. We propose that natural variation at these sites modulate both the thermal noise and spectral shifting in Baikal cottoid visual pigments resulting in adaptations that enable vision in deep water light environments.

Structure and activation of the TSH receptor transmembrane domain

PURPOSE

The thyroid-stimulating hormone receptor (TSHR) is the target autoantigen for TSHR-stimulating autoantibodies in Graves' disease. The TSHR is composed of: a leucine-rich repeat domain (LRD), a hinge region or cleavage domain (CD) and a transmembrane domain (TMD). The binding arrangements between the TSHR LRD and the thyroid-stimulating autoantibody M22 or TSH have become available from the crystal structure of the TSHR LRD-M22 complex and a comparative model of the TSHR LRD in complex with TSH, respectively. However, the mechanism by which the TMD of the TSHR and the other glycoprotein hormone receptors (GPHRs) becomes activated is unknown.

METHODS

We have generated comparative models of the structures of the inactive (TMD_In) and active (TMD_Ac) conformations of the TSHR, follicle-stimulating hormone receptor (FSHR) and luteinizing hormone receptor (LHR) TMDs. The structures of TMD_Ac and TMD_In were obtained using class A GPCR crystal structures for which fully active and inactive conformations were available.

RESULTS

Most conserved motifs observed in GPCR TMDs are also observed in the amino acid sequences of GPHR TMDs. Furthermore, most GPCR TMD conserved helix distortions are observed in our models of the structures of GPHR TMDs. Analysis of these structures has allowed us to propose a mechanism for activation of GPHR TMDs.

CONCLUSIONS

Insight into the mechanism of activation of the TSHR by both TSH and TSHR autoantibodies is likely to be useful in the development of new treatments for Graves' disease.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Structural and Functional Effect of an Oscillating Electric Field on the Dopamine-D3 Receptor: A Molecular Dynamics Simulation Study

Dopamine as a neurotransmitter plays a critical role in the functioning of the central nervous system. The structure of D3 receptor as a member of class A G-protein coupled receptors (GPCRs) has been reported. We used MD simulation to investigate the effect of an oscillating electric field, with frequencies in the range 0.6–800 GHz applied along the z-direction, on the dopamine-D3R complex. The simulations showed that at some frequencies, the application of an external oscillating electric field along the z-direction has a considerable effect on the dopamine-D3R. However, there is no enough evidence for prediction of changes in specific frequency, implying that there is no order in changes. Computing the correlation coefficient parameter showed that increasing the field frequency can weaken the interaction between dopamine and D3R and may decrease the Arg128{3.50}-Glu324{6.30} distance. Because of high stability of α helices along the z-direction, applying an oscillating electric field in this direction with an amplitude 10-time higher did not have a considerable effect. However, applying the oscillating field at the frequency of 0.6 GHz along other directions, such as X-Y and Y-Z planes, could change the energy between the dopamine and the D3R, and the number of internal hydrogen bonds of the protein. This can be due to the effect of the direction of the electric field vis-à-vis the ligands orientation and the interaction of the oscillating electric field with the dipole moment of the protein.

Coupled HOOP signature correlates with quantum yield of isorhodopsin and analog pigments

With a quantum yield of 0.66±0.03 the photoisomerization efficiency of the visual pigment rhodopsin (11-cis⇒all-trans chromophore) is exceptionally high. This is currently explained by coherent coupling of the excited state electronic wavepacket with local vibrational nuclear modes, facilitating efficient cross-over at a conical intersection onto the photoproduct energy surface. The 9-cis counterpart of rhodopsin, dubbed isorhodopsin, has a much lower quantum yield (0.26±0.03), which, however, can be markedly enhanced by modification of the retinal chromophore (7,8-dihydro and 9-cyclopropyl derivatives). The coherent coupling in the excited state is promoted by torsional skeletal and coupled HOOP vibrational modes, in combination with a twisted conformation around the isomerization region. Since such torsion will strongly enhance the infrared intensity of coupled HOOP modes, we investigated FTIR difference spectra of rhodopsin, isorhodopsin and several analog pigments in the spectral range of isolated and coupled HCCH wags. As a result we propose that the coupled HOOP signature in these retinal pigments correlates with the distribution of torsion over counteracting segments in the retinylidene polyene chain. As such the HOOP signature can act as an indicator for the photoisomerization efficiency, and can explain the higher quantum yield of the 7,8-dihydro and 9-cyclopropyl-isorhodopsin analogs.

Powdered G-Protein-Coupled Receptors

Preparation and storage of functional membrane proteins such as G-protein-coupled receptors (GPCRs) are crucial to the processes of drug delivery and discovery. Here, we describe a method of preparing powdered GPCRs using rhodopsin as the prototype. We purified rhodopsin in CHAPS detergent with low detergent to protein ratio so the bulk of the sample represented protein (ca. 72% w/w). Our new method for generating powders of membrane proteins followed by rehydration paves the way for conducting functional and biophysical experiments. As an illustrative application, powdered rhodopsin was prepared with and without the cofactor 11-cis-retinal to enable partial rehydration of the protein with D₂O in a controlled manner. Quasi-elastic neutron scattering studies using both spatial motion and energy landscape models form the basis for crucial insights into structural fluctuations and thermodynamics of GPCR activation.

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy

With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the β₂ AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/template sequence identity. This suggests that the likely accuracy of the homology models that can be built for a given receptor can be generally forecasted on the basis of the available templates. We also probed whether sequence alignments that allow for the presence of gaps within the transmembrane domains to account for structural irregularities afford better models than the classical alignment procedures that do not allow for the presence of gaps within such domains. As our results indicated, although the overall differences are very subtle, the inclusion of internal gaps within the transmembrane domains has a noticeable a beneficial effect on the local structural accuracy of the domain in question.

Retinal orientation and interactions in rhodopsin reveal a two-stage trigger mechanism for activation

The 11-cis retinal chromophore is tightly packed within the interior of the visual receptor rhodopsin and isomerizes to the all-trans configuration following absorption of light. The mechanism by which this isomerization event drives the outward rotation of transmembrane helix H6, a hallmark of activated G protein-coupled receptors, is not well established. To address this question, we use solid-state NMR and FTIR spectroscopy to define the orientation and interactions of the retinal chromophore in the active metarhodopsin II intermediate. Here we show that isomerization of the 11-cis retinal chromophore generates strong steric interactions between its β-ionone ring and transmembrane helices H5 and H6, while deprotonation of its protonated Schiff's base triggers the rearrangement of the hydrogen-bonding network involving residues on H6 and within the second extracellular loop. We integrate these observations with previous structural and functional studies to propose a two-stage mechanism for rhodopsin activation.

Update on melatonin receptors: IUPHAR Review 20

Melatonin receptors are seven transmembrane-spanning proteins belonging to the GPCR superfamily. In mammals, two melatonin receptor subtypes exist - MT1 and MT2 - encoded by the MTNR1A and MTNR1B genes respectively. The current review provides an update on melatonin receptors by the corresponding subcommittee of the International Union of Basic and Clinical Pharmacology. We will highlight recent developments of melatonin receptor ligands, including radioligands, and give an update on the latest phenotyping results of melatonin receptor knockout mice. The current status and perspectives of the structure of melatonin receptor will be summarized. The physiological importance of melatonin receptor dimers and biologically important and type 2 diabetes-associated genetic variants of melatonin receptors will be discussed. The role of melatonin receptors in physiology and disease will be further exemplified by their functions in the immune system and the CNS. Finally, antioxidant and free radical scavenger properties of melatonin and its relation to melatonin receptors will be critically addressed.

A cure for the blues: opsin duplication and subfunctionalization for short-wavelength sensitivity in jewel beetles (Coleoptera: Buprestidae)

Background

Arthropods have received much attention as a model for studying opsin evolution in invertebrates. Yet, relatively few studies have investigated the diversity of opsin proteins that underlie spectral sensitivity of the visual pigments within the diverse beetles (Insecta: Coleoptera). Previous work has demonstrated that beetles appear to lack the short-wavelength-sensitive (SWS) opsin class that typically confers sensitivity to the “blue” region of the light spectrum. However, this is contrary to established physiological data in a number of Coleoptera. To explore potential adaptations at the molecular level that may compensate for the loss of the SWS opsin, we carried out an exploration of the opsin proteins within a group of beetles (Buprestidae) where short-wave sensitivity has been demonstrated. RNA-seq data were generated to identify opsin proteins from nine taxa comprising six buprestid species (including three male/female pairs) across four subfamilies. Structural analyses of recovered opsins were conducted and compared to opsin sequences in other insects across the main opsin classes—ultraviolet, short-wavelength, and long-wavelength.

Results

All nine buprestids were found to express two opsin copies in each of the ultraviolet and long-wavelength classes, contrary to the single copies recovered in all other molecular studies of adult beetle opsin expression. No SWS opsin class was recovered. Furthermore, the male Agrilus planipennis (emerald ash borer—EAB) expressed a third LWS opsin at low levels that is presumed to be a larval copy. Subsequent homology and structural analyses identified multiple amino acid substitutions in the UVS and LWS copies that could confer short-wavelength sensitivity.

Conclusions

This work is the first to compare expressed opsin genes against known electrophysiological data that demonstrate multiple peak sensitivities in Coleoptera. We report the first instance of opsin duplication in adult beetles, which occurs in both the UVS and LWS opsin classes. Through structural comparisons of known insect opsins, we suggest that opsin duplication and amino acid variation within the chromophore binding pocket explains sensitivity in the short-wavelength portion of the visible light spectrum in these species. These findings are the first to reveal molecular complexity of the color vision system within beetles.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0674-4) contains supplementary material, which is available to authorized users.

Retinal Conformation Changes Rhodopsin's Dynamic Ensemble

G protein-coupled receptors are vital membrane proteins that allosterically transduce biomolecular signals across the cell membrane. However, the process by which ligand binding induces protein conformation changes is not well understood biophysically. Rhodopsin, the mammalian dim-light receptor, is a unique test case for understanding these processes because of its switch-like activity; the ligand, retinal, is bound throughout the activation cycle, switching from inverse agonist to agonist after absorbing a photon. By contrast, the ligand-free opsin is outside the activation cycle and may behave differently. We find that retinal influences rhodopsin dynamics using an ensemble of all-atom molecular dynamics simulations that in aggregate contain 100 μs of sampling. Active retinal destabilizes the inactive state of the receptor, whereas the active ensemble was more structurally homogenous. By contrast, simulations of an active-like receptor without retinal present were much more heterogeneous than those containing retinal. These results suggest allosteric processes are more complicated than a ligand inducing protein conformational changes or simply capturing a shifted ensemble as outlined in classic models of allostery.

A comparative study of rhodopsin function in the great bowerbird (Ptilonorhynchus nuchalis): Spectral tuning and light-activated kinetics

Rhodopsin is the visual pigment responsible for initiating the phototransduction cascade in vertebrate rod photoreceptors. Although well-characterized in a few model systems, comparative studies of rhodopsin function, particularly for nonmammalian vertebrates are comparatively lacking. Bowerbirds are rare among passerines in possessing a key substitution, D83N, at a site that is otherwise highly conserved among G protein-coupled receptors. While this substitution is present in some dim-light adapted vertebrates, often accompanying another unusual substitution, A292S, its functional relevance in birds is uncertain. To investigate functional effects associated with these two substitutions, we use the rhodopsin gene from the great bowerbird (Ptilonorhynchus nuchalis) as a background for site-directed mutagenesis, in vitro expression and functional characterization. We also mutated these sites in two additional rhodopsins that do not naturally possess N83, chicken and bovine, for comparison. Both sites were found to contribute to spectral blue-shifts, but had opposing effects on kinetic rates. Substitutions at site 83 were found to primarily affect the kinetics of light-activated rhodopsin, while substitutions at site 292 had a larger impact on spectral tuning. The contribution of substitutions at site 83 to spectral tuning in particular depended on genetic background, but overall, the effects of substitutions were otherwise surprisingly additive, and the magnitudes of functional shifts were roughly similar across all three genetic backgrounds. By employing a comparative approach with multiple species, our study provides new insight into the joint impact of sites 83 and 292 on rhodopsin structure-function as well as their evolutionary significance for dim-light vision across vertebrates.

In Silico Modeling of Novel Drug Ligands for Treatment of Concussion Associated Tauopathy

The objective of this study was to develop an in silico screening model for characterization of potential novel ligands from commercial drug libraries able to functionally activate certain olfactory receptors (ORs), which are members of the class A rhodopsin-like family of G protein couple receptors (GPCRs), in the brain of murine models of concussion. We previously found that concussions may significantly influence expression of certain ORs, for example, OR4M1 in subjects with a history of concussion/traumatic brain injury (TBI). In this study, we built a 3-D OR4M1 model and used it in in silico screening of potential novel ligands from commercial drug libraries. We report that in vitro activation of OR4M1 with the commercially available ZINC library compound 10915775 led to a significant attenuation of abnormal tau phosphorylation in embryonic cortico-hippocampal neuronal cultures derived from NSE-OR4M1 transgenic mice, possibly through modulation of the JNK signaling pathway. The attenuation of abnormal tau phosphorylation was rather selective since ZINC10915775 significantly decreased tau phosphorylation on tau Ser202/T205 (AT8 epitope) and tau Thr212/Ser214 (AT100 epitope), but not on tau Ser396/404 (PHF-1 epitope). Moreover, no response of ZINC10915775 was found in control hippocampal neuronal cultures derived from wild type littermates. Our in silico model provides novel means to pharmacologically modulate select ubiquitously expressed ORs in the brain through high affinity ligand activation to prevent and eventually to treat concussion induced down regulation of ORs and subsequent cascade of tau pathology. J. Cell. Biochem. 117: 2241-2248, 2016. © 2016 Wiley Periodicals, Inc.

Computational analysis of the CB1 carboxyl-terminus in the receptor-G protein complex

Despite the important role of the carboxyl-terminus (Ct) of the activated brain cannabinoid receptor one (CB1) in the regulation of G protein signaling, a structural understanding of interactions with G proteins is lacking. This is largely due to the highly flexible nature of the CB1 Ct that dynamically adapts its conformation to the presence of G proteins. In the present study, we explored how the CB1 Ct can interact with the G protein by building on our prior modeling of the CB1-Gi complex (Shim, Ahn, and Kendall, The Journal of Biological Chemistry 2013;288:32449-32465) to incorporate a complete CB1 Ct (Glu416(Ct)-Leu472(Ct)). Based on the structural constraints from NMR studies, we employed ROSETTA to predict tertiary folds, ZDOCK to predict docking orientation, and molecular dynamics (MD) simulations to obtain two distinct plausible models of CB1 Ct in the CB1-Gi complex. The resulting models were consistent with the NMR-determined helical structure (H9) in the middle region of the CB1 Ct. The CB1 Ct directly interacted with both Gα and Gβ and stabilized the receptor at the Gi interface. The results of site-directed mutagenesis studies of Glu416(Ct), Asp423(Ct), Asp428(Ct), and Arg444(Ct) of CB1 Ct suggested that the CB1 Ct can influence receptor-G protein coupling by stabilizing the receptor at the Gi interface. This research provided, for the first time, models of the CB1 Ct in contact with the G protein.

CONDENSED-MATTER SPECTROSCOPY SPECTRAL METHODS FOR STUDY OF THE G-PROTEIN-COUPLED RECEPTOR RHODOPSIN. II. MAGNETIC RESONANCE METHODS

This article continues our review of spectroscopic studies of G-protein-coupled receptors. Magnetic resonance methods including electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) provide specific structural and dynamical data for the protein in conjunction with optical methods (vibrational, electronic spectroscopy) as discussed in the accompanying article. An additional advantage is the opportunity to explore the receptor proteins in the natural membrane lipid environment. Solid-state ²H and ¹³C NMR methods yield information about the both local structure and dynamics of the cofactor bound to the protein and its light induced changes. Complementary site-directed spin labeling studies monitor the structural alterations over larger distances and correspondingly longer time scales. A multi-scale reaction mechanism describes how local changes of the retinal cofactor unlock the receptor to initiate large-scale conformational changes of rhodopsin. Activation of the G-protein-coupled receptor involves an ensemble of conformational substates within the rhodopsin manifold that characterize the dynamically active receptor.

Rhodopsins carrying modified chromophores--the 'making of', structural modelling and their light-induced reactivity

A series of vitamin-A aldehydes (retinals) with modified alkyl group substituents (9-demethyl-, 9-ethyl-, 9-isopropyl-, 10-methyl, 10-methyl-13-demethyl-, and 13-demethyl retinal) was synthesized and their 11-cis isomers were used as chromophores to reconstitute the visual pigment rhodopsin. Structural changes were selectively introduced around the photoisomerizing C11=C12 bond. The effect of these structural changes on rhodopsin formation and bleaching was determined. Global fit of assembly kinetics yielded lifetimes and spectral features of the assembly intermediates. Rhodopsin formation proceeds stepwise with prolonged lifetimes especially for 9-demethyl retinal (longest lifetime τ3 = 7500 s, cf., 3500 s for retinal), and for 10-methyl retinal (τ3 = 7850 s). These slowed-down processes are interpreted as either a loss of fixation (9dm) or an increased steric hindrance (10me) during the conformational adjustment within the protein. Combined quantum mechanics and molecular mechanics (QM/MM) simulations provided structural insight into the retinal analogues-assembled, full-length rhodopsins. Extinction coefficients, quantum yields and kinetics of the bleaching process (μs-to-ms time range) were determined. Global fit analysis yielded lifetimes and spectral features of bleaching intermediates, revealing remarkably altered kinetics: whereas the slowest process of wild-type rhodopsin and of bleached and 11-cis retinal assembled rhodopsin takes place with lifetimes of 7 and 3.8 s, respectively, this process for 10-methyl-13-demethyl retinal was nearly 10 h (34670 s), coming to completion only after ca. 50 h. The structural changes in retinal derivatives clearly identify the precise interactions between chromophore and protein during the light-induced changes that yield the outstanding efficiency of rhodopsin.

Computational Prediction and Biochemical Analyses of New Inverse Agonists for the CB1 Receptor

Human cannabinoid type 1 (CB1) G-protein coupled receptor is a potential therapeutic target for obesity. The previously predicted and experimentally validated ensemble of ligand-free conformations of CB1 [Scott, C. E. et al. Protein Sci. 2013 , 22 , 101 - 113 ; Ahn, K. H. et al. Proteins 2013 , 81 , 1304 - 1317] are used here to predict the binding sites for known CB1-selective inverse agonists including rimonabant and its seven known derivatives. This binding pocket, which differs significantly from previously published models, is used to identify 16 novel compounds expected to be CB1 inverse agonists by exploiting potential new interactions. We show experimentally that two of these compounds exhibit inverse agonist properties including inhibition of basal and agonist-induced G-protein coupling activity, as well as an enhanced level of CB1 cell surface localization. This demonstrates the utility of using the predicted binding sites for an ensemble of CB1 receptor structures for designing new CB1 inverse agonists.

Comparative sequence analyses of rhodopsin and RPE65 reveal patterns of selective constraint across hereditary retinal disease mutations

Retinitis pigmentosa (RP) comprises several heritable diseases that involve photoreceptor, and ultimately retinal, degeneration. Currently, mutations in over 50 genes have known links to RP. Despite advances in clinical characterization, molecular characterization of RP remains challenging due to the heterogeneous nature of causal genes, mutations, and clinical phenotypes. In this study, we compiled large datasets of two important visual genes associated with RP: rhodopsin, which initiates the phototransduction cascade, and the retinoid isomerase RPE65, which regenerates the visual cycle. We used a comparative evolutionary approach to investigate the relationship between interspecific sequence variation and pathogenic mutations that lead to degenerative retinal disease. Using codon-based likelihood methods, we estimated evolutionary rates (dN/dS) across both genes in a phylogenetic context to investigate differences between pathogenic and nonpathogenic amino acid sites. In both genes, disease-associated sites showed significantly lower evolutionary rates compared to nondisease sites, and were more likely to occur in functionally critical areas of the proteins. The nature of the dataset (e.g., vertebrate or mammalian sequences), as well as selection of pathogenic sites, affected the differences observed between pathogenic and nonpathogenic sites. Our results illustrate that these methods can serve as an intermediate step in understanding protein structure and function in a clinical context, particularly in predicting the relative pathogenicity (i.e., functional impact) of point mutations and their downstream phenotypic effects. Extensions of this approach may also contribute to current methods for predicting the deleterious effects of candidate mutations and to the identification of protein regions under strong constraint where we expect pathogenic mutations to occur.

Theoretical Studies on the Color-Tuning Mechanism in Retinal Proteins

The excited states of the three retinal proteins, bovine rhodopsin (Rh), bacteriorhodopsin (bR), and sensory rhodopsin II (sRII) were studied using the symmetry-adapted cluster-configuration interaction (SAC-CI) and combined quantum mechanical and molecular mechanical (QM/MM) methods. The computed absorption energies are in good agreement with the experimental ones for all three proteins. The spectral tuning mechanism was analyzed in terms of three contributions:  molecular structures of the chromophore in the binding pockets, electrostatic (ES) interaction of the chromophore with the surrounding protein environment, and quantum-mechanical effect between the chromophore and the counterion group. This analysis provided an insight into the mechanism of the large blue-shifts in the absorption peak position of Rh and sRII from that of bR. Protein ES effect is primarily important both in Rh and in sRII, and the structure effect is secondary important in Rh. The quantum-mechanical interaction between the chromophore and the counterion is very important for quantitative reproduction of the excitation energy. These results indicate that the present approach is useful for studying the absorption spectra and the mechanism of the color tuning in the retinal proteins.

Protein Influence on Electronic Spectra Modeled by Multipoles and Polarizabilities

We have developed automatic methods to calculate multipoles and anisotropic polarizabilities for all atoms and bond centers in a protein and to include such a model in the calculation of electronic properties at any level of quantum mechanical theory. This approach is applied for the calculation of the electronic spectra of retinal in rhodopsin at the CASPT2//CASSCF level (second-order multiconfigurational perturbation theory) for the wild-type protein, as well as two mutants and isorhodopsin in QM/MM structures based on two crystal structures. We also perform a detailed investigation of the importance and distance dependence of the multipoles and the polarizabilities for both the absolute and the relative absorption energies. It is shown that the model of the surrounding protein strongly influences the spectrum and that different models give widely different results. For example, the Amber 1994 and 2003 force fields give excitation energies that differ by up to 16 kJ/mol. For accurate excitation energies, multipoles up to quadrupoles and anisotropic polarizabilities are needed. However, interactions with residues more than 10 Å from the chromophore can be treated with a standard polarizable force field without any dipoles or quadrupoles.

Pharmacophore-Map-Pick: A Method to Generate Pharmacophore Models for All Human GPCRs

GPCR-based drug discovery is hindered by a lack of effective screening methods for most GPCRs that have neither ligands nor high-quality structures. With the aim to identify lead molecules for these GPCRs, we developed a new method called Pharmacophore-Map-Pick to generate pharmacophore models for all human GPCRs. The model of ADRB2 generated using this method not only predicts the binding mode of ADRB2-ligands correctly but also performs well in virtual screening. Findings also demonstrate that this method is powerful for generating high-quality pharmacophore models. The average enrichment for the pharmacophore models of the 15 targets in different GPCR families reached 15-fold at 0.5 % false-positive rate. Therefore, the pharmacophore models can be applied in virtual screening directly with no requirement for any ligand information or shape constraints. A total of 2386 pharmacophore models for 819 different GPCRs (99 % coverage (819/825)) were generated and are available at http://bsb.kiz.ac.cn/GPCRPMD.

Responsiveness of G protein-coupled odorant receptors is partially attributed to the activation mechanism

Mammals detect and discriminate numerous odors via a large family of G protein-coupled odorant receptors (ORs). However, little is known about the molecular and structural basis underlying OR response properties. Using site-directed mutagenesis and computational modeling, we studied ORs sharing high sequence homology but with different response properties. When tested in heterologous cells by diverse odorants, MOR256-3 responded broadly to many odorants, whereas MOR256-8 responded weakly to a few odorants. Out of 36 mutant MOR256-3 ORs, the majority altered the responses to different odorants in a similar manner and the overall response of an OR was positively correlated with its basal activity, an indication of ligand-independent receptor activation. Strikingly, a single mutation in MOR256-8 was sufficient to confer both high basal activity and broad responsiveness to this receptor. These results suggest that broad responsiveness of an OR is at least partially attributed to its activation likelihood.

Transmembrane signal transduction by peptide hormones via family B G protein-coupled receptors

Although family B G protein-coupled receptors (GPCRs) contain only 15 members, they play key roles in transmembrane signal transduction of hormones. Family B GPCRs are drug targets for developing therapeutics for diseases ranging from metabolic to neurological disorders. Despite their importance, the molecular mechanism of activation of family B GPCRs remains largely unexplored due to the challenges in expression and purification of functional receptors to the quantity for biophysical characterization. Currently, there is no crystal structure available of a full-length family B GPCR. However, structures of key domains, including the extracellular ligand binding regions and seven-helical transmembrane regions, have been solved by X-ray crystallography and NMR, providing insights into the mechanisms of ligand recognition and selectivity, and helical arrangements within the cell membrane. Moreover, biophysical and biochemical methods have been used to explore functions, key residues for signaling, and the kinetics and dynamics of signaling processes. This review summarizes the current knowledge of the signal transduction mechanism of family B GPCRs at the molecular level and comments on the challenges and outlook for mechanistic studies of family B GPCRs.

Eukaryotic G protein-coupled receptors as descendants of prokaryotic sodium-translocating rhodopsins

Abstract

Microbial rhodopsins and G-protein coupled receptors (GPCRs, which include animal rhodopsins) are two distinct (super) families of heptahelical (7TM) membrane proteins that share obvious structural similarities but no significant sequence similarity. Comparison of the recently solved high-resolution structures of the sodium-translocating bacterial rhodopsin and various Na⁺-binding GPCRs revealed striking similarity of their sodium-binding sites. This similarity allowed us to construct a structure-guided sequence alignment for the two (super)families, which highlighted their evolutionary relatedness. Our analysis supports a common underlying molecular mechanism for both families that involves a highly conserved aromatic residue playing a pivotal role in rotation of the 6th transmembrane helix.

Reviewers

This article was reviewed by Oded Beja, G. P. S. Raghava and L. Aravind.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0091-4) contains supplementary material, which is available to authorized users.

Explicit spatiotemporal simulation of receptor-G protein coupling in rod cell disk membranes

Dim-light vision is mediated by retinal rod cells. Rhodopsin (R), a G-protein-coupled receptor, switches to its active form (R(∗)) in response to absorbing a single photon and activates multiple copies of the G-protein transducin (G) that trigger further downstream reactions of the phototransduction cascade. The classical assumption is that R and G are uniformly distributed and freely diffusing on disk membranes. Recent experimental findings have challenged this view by showing specific R architectures, including RG precomplexes, nonuniform R density, specific R arrangements, and immobile fractions of R. Here, we derive a physical model that describes the first steps of the photoactivation cascade in spatiotemporal detail and single-molecule resolution. The model was implemented in the ReaDDy software for particle-based reaction-diffusion simulations. Detailed kinetic in vitro experiments are used to parametrize the reaction rates and diffusion constants of R and G. Particle diffusion and G activation are then studied under different conditions of R-R interaction. It is found that the classical free-diffusion model is consistent with the available kinetic data. The existence of precomplexes between inactive R and G is only consistent with the data if these precomplexes are weak, with much larger dissociation rates than suggested elsewhere. Microarchitectures of R, such as dimer racks, would effectively immobilize R but have little impact on the diffusivity of G and on the overall amplification of the cascade at the level of the G protein.

The role of the non-covalent β-ionone-ring binding site in rhodopsin: historical and physiological perspective

Bleached rhodopsin regenerates by way of the Schiff base formation between the 11-cis retinal and opsin. Recovery of human vision from light adapted states follows biphasic kinetics and each adaptive phase is assigned to two distinct classes of visual pigments in cones and rods, respectively, suggesting that the speed of Schiff base formation differs between iodopsin and rhodopsin. Matsumoto and Yoshizawa predicted the existence of a β-ionone ring-binding site in rhodopsin, which has been proven by structural studies. They postulated that rhodopsin regeneration starts with a non-covalent binding of the β-ionone ring moiety of 11-cis-retinal, followed by the Schiff base formation. Recent physiological investigation revealed that non-covalent occupation of the β-ionone ring binding site transiently activates the visual transduction cascade in the dark. In order to understand the role of non-covalent binding of 11-cis-retinal to opsin during regeneration, we studied the kinetics of rhodopsin regeneration from opsin and 11-cis-retinal and found that the Schiff base formation is accelerated ∼10(7) times compared to that between retinal and free amine. According to Cordes and Jencks, Schiff base formation in solution exhibits a bell-shaped pH dependence. However, we discovered that the rhodopsin formation is independent of pH over a wide pH range, suggesting that aqueous solvents do not have access to the Schiff base milieu during its formation. According to Hecht et al. the regeneration of iodopsin must be significantly faster than that of rhodopsin. Does this suggest that the Schiff base formation in iodopsin is favored due to its structural architecture? The iodopsin structure once solved would answer such a question as how molecular fine-tuning of retinal proteins realizes their dark adaptive functions. In contrast, bacteriorhodopsin does not require occupancy of a distinct β-ionone ring-binding site, enabling an aldehyde without the cyclohexene ring to form a pigment. Studies of regeneration reaction of other retinal proteins, which are scarcely available, would clarify the molecular structure-phenotype relationships and their physiological roles.

Analysis of Conserved Glutamate and Aspartate Residues in Drosophila Rhodopsin 1 and Their Influence on Spectral Tuning

The molecular mechanisms that regulate invertebrate visual pigment absorption are poorly understood. Studies of amphioxus Go-opsin have demonstrated that Glu-181 functions as the counterion in this pigment. This finding has led to the proposal that Glu-181 may function as the counterion in other invertebrate visual pigments as well. Here we describe a series of mutagenesis experiments to test this hypothesis and to also test whether other conserved acidic amino acids in Drosophila Rhodopsin 1 (Rh1) may serve as the counterion of this visual pigment. Of the 5 Glu and Asp residues replaced by Gln or Asn in our experiments, none of the mutant pigments shift the absorption of Rh1 by more than 6 nm. In combination with prior studies, these results suggest that the counterion in Drosophila Rh1 may not be located at Glu-181 as in amphioxus, or at Glu-113 as in bovine rhodopsin. Conversely, the extremely low steady state levels of the E194Q mutant pigment (bovine opsin site Glu-181), and the rhabdomere degeneration observed in flies expressing this mutant demonstrate that a negatively charged residue at this position is essential for normal rhodopsin function in vivo. This work also raises the possibility that another residue or physiologic anion may compensate for the missing counterion in the E194Q mutant.

Comparative Mutagenesis Studies of Retinal Release in Light-Activated Zebrafish Rhodopsin Using Fluorescence Spectroscopy

Rhodopsin is the visual pigment responsible for initiating scotopic (dim-light) vision in vetebrates. Once activated by light, release of all-trans-retinal from rhodopsin involves hydrolysis of the Schiff base linkage, followed by dissociation of retinal from the protein moiety. This kinetic process has been well studied in model systems such as bovine rhodopsin, but not in rhodopsins from cold-blooded animals, where physiological temperatures can vary considerably. Here, we characterize the rate of retinal release from light-activated rhodopsin in an ectotherm, zebrafish (Danio rerio), demonstrating in a fluorescence assay that this process occurs more than twice as fast as bovine rhodopsin at similar temperatures in 0.1% dodecyl maltoside. Using site-directed mutagenesis, we found that differences in retinal release rates can be attributed to a series of variable residues lining the retinal channel in three key structural motifs: an opening in metarhodopsin II between transmembrane helix 5 (TM5) and TM6, in TM3 near E122, and in the "retinal plug" formed by extracellular loop 2 (EL2). The majority of these sites are more proximal to the β-ionone ring of retinal than the Schiff base, indicating their influence on retinal release is more likely due to steric effects during retinal dissociation, rather than alterations to Schiff base stability. An Arrhenius plot of zebrafish rhodopsin was consistent with this model, inferring that the activation energy for Schiff base hydrolysis is similar to that of bovine rhodopsin. Functional variation at key sites identified in this study is consistent with the idea that retinal release might be an adaptive property of rhodopsin in vertebrates. Our study is one of the few investigating a nonmammalian rhodopsin, which will help establish a better understanding of the molecular mechanisms contributing to vision in cold-blooded vertebrates.

Conserved Residues Control Activation of Mammalian G Protein-Coupled Odorant Receptors

Odorant receptor (OR) genes and proteins represent more than 2% of our genome and 4% of our proteome and constitute the largest subgroup of G protein-coupled receptors (GPCRs). The mechanism underlying OR activation remains poorly understood, as they do not share some of the highly conserved motifs critical for activation of non-olfactory GPCRs. By combining site-directed mutagenesis, heterologous expression, and molecular dynamics simulations that capture the conformational change of constitutively active mutants, we tentatively identified crucial residues for the function of these receptors using the mouse MOR256-3 (Olfr124) as a model. The toggle switch for sensing agonists involves a highly conserved tyrosine residue in helix VI. The ionic lock is located between the "DRY" motif in helix III and a positively charged "R/K" residue in helix VI. This study provides an unprecedented model that captures the main mechanisms of odorant receptor activation.