Light activation of rhodopsin: insights from molecular dynamics simulations guided by solid-state NMR distance restraints

Interaction of Synthetic Cannabinoid Receptor Agonists with Cannabinoid Receptor I: Insights into Activation Molecular Mechanism

Synthetic cannabinoid receptor agonists (SCRAs) have become a wide group of new psychoactive substances since the 2010s. For the last few years, the X-ray structures of the complexes of cannabinoid receptor I (CB1) with SCRAs as well as the complexes of CB1 with its antagonist have been published. Based on those data, SCRA-CB1 interactions are analyzed in detail, using molecular modeling and molecular dynamics simulations. The molecular mechanism of the conformational transformation of the transmembrane domain of CB1 caused by its interaction with SCRA is studied. These conformational changes allosterically modulate the CB1-Gi complex, providing activation of the Gi protein. Based on the X-ray-determined structures of the CB1-ligand complexes, a stable apo conformation of inactive CB1 with a relatively low potential barrier of receptor activation was modeled. For that model, molecular dynamic simulations of SCRA binding to CB1 led to the active state of CB1, which allowed us to explore the key features of this activation and the molecular mechanism of the receptor's structural transformation. The simulated CB1 activation is in accordance with the previously published experimental data for the activation at protein mutations or structural changes of ligands. The key feature of the suggested activation mechanism is the determination of the stiff core of the CB1 transmembrane domain and the statement that the entire conformational transformation of the receptor to the active state is caused by a shift of alpha helix TM7 relative to this core. The shift itself is caused by protein-ligand interactions. It was verified via steered molecular dynamics simulations of the X-ray-determined structures of the inactive receptor, which resulted in the active conformation of CB1 irrespective of the placement of agonist ligand in the receptor's active site.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Residue-Residue Mutual Work Analysis of Retinal-Opsin Interaction in Rhodopsin: Implications for Protein-Ligand Binding

Energetic contributions at the single-residue level for retinal-opsin interactions in rhodopsin were studied by combining molecular dynamics simulations, transition path sampling, and a newly developed energy decomposition approach. The virtual work at an infinitesimal time interval was decomposed into the work components on one residue due to its interaction with another residue, which were then averaged over the transition path ensemble along a proposed reaction coordinate. Such residue-residue mutual work analysis on 62 residues within the active center of rhodopsin resulted in a very sparse interaction matrix, which is generally not symmetric but antisymmetric to some extent. Fourteen residues were identified to be major players in retinal relaxation along a plausible pathway from bathorhodopsin to the blue-shifted intermediate, which is in good agreement with an existing NMR study. Based on the matrix of mutual work, a comprehensive network was constructed to provide detailed insights into the chromophore-protein interaction from a viewpoint of energy flow.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Insight into the chromophore of rhodopsin and its Meta-II photointermediate by 19F solid-state NMR and chemical shift tensor calculations

19F nuclei are useful labels in solid-state NMR studies, since their chemical shift and tensor elements are very sensitive to the electrostatic and space-filling properties of their local environment. In this study we have exploited a fluorine substituent, strategically placed at the C-12-position of 11-cis retinal, the chromophore of visual rhodopsins. This label was used to explore the local environment of the chromophore in the ground state of bovine rhodopsin and its active photo-intermediate Meta II. In addition, the chemical shift and tensor elements of the chromophore in the free state in a membrane environment and the bound state in the protein were determined. Upon binding of the chromophore into rhodopsin and Meta II, the isotropic chemical shift changes in the opposite direction by +9.7 and -8.4 ppm, respectively. An unusually large isotropic shift difference of 35.9 ppm was observed between rhodopsin and Meta II. This partly originates in the light-triggered 11-cis to all-trans isomerization of the chromophore. The other part reflects the local conformational rearrangements in the chromophore and the binding pocket. These NMR data were correlated with the available X-ray structures of rhodopsin and Meta II using bond polarization theory. For this purpose hydrogen atoms have to be inserted and hereto a family of structures were derived that best correlated with the well-established 13C chemical shifts. Based upon these structures, a 12-F derivative was obtained that best corresponded with the experimentally determined 19F chemical shifts and tensor elements. The combined data indicate strong changes in the local environment of the C-12 position and a substantially different interaction pattern with the protein in Meta II as compared to rhodopsin.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

New views on phototransduction from atomic force microscopy and single molecule force spectroscopy on native rods

By combining atomic force microscopy (AFM) imaging and single-molecule force spectroscopy (SMFS), we analyzed membrane proteins of the rod outer segments (OS). With this combined approach we were able to study the membrane proteins in their natural environment. In the plasma membrane we identified native cyclic nucleotide-gated (CNG) channels which are organized in single file strings. We also identified rhodopsin located both in the discs and in the plasma membrane. SMFS reveals strikingly different mechanical properties of rhodopsin unfolding in the two environments. Molecular dynamic simulations suggest that this difference is likely to be related to the higher hydrophobicity of the plasma membrane, due to the higher cholesterol concentration. This increases rhodopsin mechanical stability lowering the rate of transition towards its active form, hindering, in this manner, phototransduction.

Photoactivation Intermediates of a G-Protein Coupled Receptor Rhodopsin Investigated by a Hybrid Molecular Simulation

Rhodopsin is a G-protein coupled receptor functioning as a photoreceptor for vision through photoactivation of a covalently bound ligand of a retinal protonated Schiff base chromophore. Despite the availability of structural information on the inactivated and activated forms of the receptor, the transition processes initiated by the photoabsorption have not been well understood. Here we theoretically examined the photoactivation processes by means of molecular dynamics (MD) simulations and ab initio quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimizations which enabled accurate geometry determination of the ligand molecule in ample statistical conformational samples of the protein. Structures of the intermediate states of the activation process, blue-shifted intermediate and Lumi, as well as the dark state first generated by MD simulations and then refined by the QM/MM free energy geometry optimizations were characterized by large displacement of the β-ionone ring of retinal along with change in the hydrogen bond of the protonated Schiff base. The ab initio calculations of vibrational and electronic spectroscopic properties of those states well reproduced the experimental observations and successfully identified the molecular origins underlying the spectroscopic features. The structural evolution in the formation of the intermediates provides a molecular insight into the efficient activation processes of the receptor.

Propagation of the Allosteric Modulation Induced by Sodium in the δ-Opioid Receptor

Allosteric sodium in the helix bundle of a G protein-coupled receptor (GPCR) can modulate the receptor activation on the intracellular side. This phenomenon has confounded the GPCR community for decades. In this work, we present a theoretical model that reveals the mechanism of the allosteric modulation induced by sodium in the δ-opioid receptor. We found that the allosteric sodium ion exploits a distinct conformation of the key residue Trp274^6.48 to propagate the modulation to helices 5 and 6, which further transmits along the helices and regulates their positions on the intracellular side. This mechanism is supported by subsequent functional assays. Remarkably, our results highlight the contrast between the allosteric effects towards two GPCR partners, the G protein and β-arrestin, as indicated by the fact that the allosteric modulation initiated by the sodium ion significantly affects the β-arrestin recruitment, while it alters the G protein signaling only moderately. We believe that the mechanism revealed in this work can be used to explain allosteric effects initiated by sodium in other GPCRs since the allosteric sodium is highly conserved across GPCRs.

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.

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.

Photosynthetic Proteins in Supported Lipid Bilayers: Towards a Biokleptic Approach for Energy Capture

In nature, plants and some bacteria have evolved an ability to convert solar energy into chemical energy usable by the organism. This process involves several proteins and the creation of a chemical gradient across the cell membrane. To transfer this process to a laboratory environment, several conditions have to be met: i) proteins need to be reconstituted into a lipid membrane, ii) the proteins need to be correctly oriented and functional and, finally, iii) the lipid membrane should be capable of maintaining chemical and electrical gradients. Investigating the processes of photosynthesis and energy generation in vivo is a difficult task due to the complexity of the membrane and its associated proteins. Solid, supported lipid bilayers provide a good model system for the systematic investigation of the different components involved in the photosynthetic pathway. In this review, the progress made to date in the development of supported lipid bilayer systems suitable for the investigation of membrane proteins is described; in particular, there is a focus on those used for the reconstitution of proteins involved in light capture.

Construction of stable mammalian cell lines for inducible expression of G protein-coupled receptors

The large-scale expression of many membrane proteins, including the members of the G protein-coupled receptor superfamily, in a correctly folded and fully functional form remains a formidable challenge. In this chapter, we focus on the construction of stable mammalian cell lines to overcome this hurdle. First, we will outline the steps for establishing a tightly regulated gene expression system in human HEK293S cells. This system utilizes separate plasmids containing components of well-defined genetic control elements from the Escherichia coli tetracycline operon to control the powerful cytomegalovirus immediate early enhancer/promoter. Next, we describe the assembly of this expression system into HEK293S cells and a derivative cell line devoid of complex N-glycosylation. Finally, we describe methods for the growth of these cells lines in scalable suspension culture for the preparation of milligram amounts of recombinant protein.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

Structure and dynamics of G-protein coupled receptors

G-protein coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell-to-cell communication. They also form the largest superfamily of drug targets. Hence detailed studies of the three dimensional structure and dynamics are critical to understanding the functional role of GPCRs in signal transduction pathways, and for drug design. In this chapter we compare the features of the crystal structures of various biogenic amine receptors, such as β1 and β2 adrenergic receptors, dopamine D3 receptor, M2 and M3 muscarinic acetylcholine receptors. This analysis revealed that conserved residues are located facing the inside of the transmembrane domain in these GPCRs improving the efficiency of packing of these structures. The NMR structure of the chemokine receptor CXCR1 without any ligand bound, shows significant dynamics of the transmembrane domain, especially the helical kink angle on the transmembrane helix6. The activation mechanism of the β2-adrenergic receptor has been studied using multiscale computational methods. The results of these studies showed that the receptor without any ligand bound, samples conformations that resemble some of the structural characteristics of the active state of the receptor. Ligand binding stabilizes some of the conformations already sampled by the apo receptor. This was later observed in the NMR study of the dynamics of human β2-adrenergic receptor. The dynamic nature of GPCRs leads to a challenge in obtaining purified receptors for biophysical studies. Deriving thermostable mutants of GPCRs has been a successful strategy to reduce the conformational heterogeneity and stabilize the receptors. This has lead to several crystal structures of GPCRs. However, the cause of how these mutations lead to thermostability is not clear. Computational studies are beginning to shed some insight into the possible structural basis for the thermostability. Molecular Dynamics simulations studying the conformational ensemble of thermostable mutants have shown that the stability could arise from both enthalpic and entropic factors. There are regions of high stress in the wild type GPCR that gets relieved upon mutation conferring thermostability.

Structural approaches to understanding retinal proteins needed for vision

The past decade has witnessed an impressive expansion of our knowledge of retinal photoreceptor signal transduction and the regulation of the visual cycle required for normal eyesight. Progress in human genetics and next generation sequencing technologies have revealed the complexity behind many inherited retinal diseases. Structural studies have markedly increased our understanding of the visual process. Moreover, technical innovations and improved methodologies in proteomics, macromolecular crystallization and high resolution imaging at different levels set the scene for even greater advances. Pharmacology combined with structural biology of membrane proteins holds great promise for developing innovative accessible therapies for millions robbed of their sight or progressing toward blindness.

Isotropic magnetic shielding constants of retinal derivatives in aprotic and protic solvents

We investigate the nuclear isotropic shielding constants σ((13)C) and σ((17)O) of isomers of retinoic acid and retinal in gas-phase and in chloroform, acetonitrile, methanol, and water solutions via Monte Carlo simulation and quantum mechanics calculations using the GIAO-B3LYP∕6-311++G(2d,2p) approach. Electronic solute polarization effects due to protic and aprotic solvents are included iteratively and play an important role in the quantitative determination of oxygen shielding constants. Our MP2∕6-31G+(d) results show substantial increases of the dipole moment of both retinal derivatives in solution as compared with the gas-phase results (between 22% and 26% in chloroform and between 55% and 99% in water). For the oxygen atoms the influence of the solute polarization is mild for σ((17)O) of hydroxyl group, even in protic solvents, but it is particularly important for σ((17)O) of carbonyl group. For the latter, there is a sizable increase in the magnitude with increasing solvent polarity. For the carbon atoms, the solvent effects on the σ((13)C) values are in general small, being more appreciable in carbon atoms of the polyene chain than in the carbon atoms of the β-ionone ring and methyl groups. The results also show that isomeric changes on the backbones of the polyene chains have marked influence on the (13)C chemical shifts of carbon atoms near to the structural distortions, in good agreement with the experimental results measured in solution.

Deconstructing activation events in rhodopsin

Activation of class-A G-protein-coupled receptors (GPCRs) involves large-scale reorganization of the H3/H6 interhelical network. In rhodopsin (Rh), this process is coupled to a change in the protonation state of a key residue, E134, whose exact role in activation is not well understood. Capturing this millisecond pH-dependent process is a well-appreciated challenge. We have developed a scheme combining the harmonic Fourier beads (HFB) method and constant-pH molecular dynamics with pH-based replica exchange (pH-REX) to gain insight into the structural changes that occur along the activation pathway as a function of the protonation state of E134. Our results indicate that E134 is protonated as a consequence of tilting of H6 by ca. 4.0° with respect to its initial position and simultaneous rotation by ca. 23° along its principal axis. The movement of H6 is associated with breakage of the E247-R135 and R135-E134 salt bridges and concomitant release of the E134 side chain, which results in an increase in its pKa value above physiological pH. An increase in the hydrophobicity of the environment surrounding E134 leads to further tilting and rotation of H6 and upshift of the E134 pKa. Such atomic-level information, which is not accessible through experiments, refines the earlier proposed sequential model of Rh activation (see: Zaitseva, E.; et al. Sequential Rearrangement of Interhelical Networks Upon Rhodopsin Activation in Membranes: The Meta IIa Conformational Substate . J. Am. Chem. Soc. 2010, 132, 4815) and argues that the E134 protonation switch is both a cause and a consequence of the H6 motion.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Thermostabilization of the β1-adrenergic receptor correlates with increased entropy of the inactive state

The dynamic nature of GPCRs is a major hurdle in their purification and crystallization. Thermostabilization can facilitate GPCR structure determination, as has been shown by the structure of the thermostabilized β1-adrenergic receptor (β1AR) mutant, m23-β1AR, which has been thermostabilized in the inactive state. However, it is unclear from the structure how the six thermostabilizing mutations in m23-β1AR affect receptor dynamics. We have used molecular dynamics simulations in explicit solvent to compare the conformational ensembles for both wild type β1AR (wt-β1AR) and m23-β1AR. Thermostabilization results in an increase in the number of accessible microscopic conformational states within the inactive state ensemble, effectively increasing the side chain entropy of the inactive state at room temperature, while suppressing large-scale main chain conformational changes that lead to activation. We identified several diverse mechanisms of thermostabilization upon mutation. These include decrease of long-range correlated movement between residues in the G-protein coupling site to the extracellular region (Y227A(5.58), F338M(7.48)), formation of new hydrogen bonds (R68S), and reduction of local stress (Y227(5.58), F327(7.37), and F338(7.48)). This study provides insights into microscopic mechanisms underlying thermostability that leads to an understanding of the effect of these mutations on the structure of the receptor.

G-protein-coupled receptor structure, ligand binding and activation as studied by solid-state NMR spectroscopy

GPCRs (G-protein-coupled receptors) are versatile signalling molecules at the cell surface and make up the largest and most diverse family of membrane receptors in the human genome. They convert a large variety of extracellular stimuli into intracellular responses through the activation of heterotrimeric G-proteins, which make them key regulatory elements in a broad range of normal and pathological processes, and are therefore one of the most important targets for pharmaceutical drug discovery. Knowledge of a GPCR structure enables us to gain a mechanistic insight into its function and dynamics, and further aid rational drug design. Despite intensive research carried out over the last three decades, resolving the structural basis of GPCR function is still a major activity. The crystal structures obtained in the last 5 years provide the first opportunity to understand how protein structure dictates the unique functional properties of these complex signalling molecules. However, owing to the intrinsic hydrophobicity, flexibility and instability of membrane proteins, it is still a challenge to crystallize GPCRs, and, when this is possible, it is no longer in its native membrane environment and no longer without modification. Furthermore, the conformational change of the transmembrane α-helices associated with the structure activation increases the difficulty of capturing the activation state of a GPCR to a higher resolution by X-ray crystallography. On the other hand, solid-state NMR may offer a unique opportunity to study membrane protein structure, ligand binding and activation at atomic resolution in the native membrane environment, as well as described functionally significant dynamics. In the present review, we discuss some recent achievements of solid-state NMR for understanding GPCRs, the largest mammalian proteome at ~1% of the total expressed proteins. Structural information, details of determination, details of ligand conformations and the consequences of ligand binding to initiate activation can all be explored with solid-state NMR.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Modulation of constitutive activity and signaling bias of the ghrelin receptor by conformational constraint in the second extracellular loop

Based on a rare, natural Glu for Ala-204(C+6) variant located six residues after the conserved Cys residue in extracellular loop 2b (ECL2b) associated with selective elimination of the high constitutive signaling of the ghrelin receptor, this loop was subjected to a detailed structure functional analysis. Introduction of Glu in different positions demonstrated that although the constitutive signaling was partly reduced when introduced in position 205(C+7) it was only totally eliminated in position 204(C+6). No charge-charge interaction partner could be identified for the Glu(C+6) variant despite mutational analysis of a number of potential partners in the extracellular loops and outer parts of the transmembrane segments. Systematic probing of position 204(C+6) with amino acid residues of different physicochemical properties indicated that a positively charged Lys surprisingly provided phenotypes similar to those of the negatively charged Glu residue. Computational chemistry analysis indicated that the propensity for the C-terminal segment of extracellular loop 2b to form an extended α-helix was increased from 15% in the wild type to 89 and 82% by introduction in position 204(C+6) of a Glu or a Lys residue, respectively. Moreover, the constitutive activity of the receptor was inhibited by Zn(2+) binding in an engineered metal ion site, stabilizing an α-helical conformation of this loop segment. It is concluded that the high constitutive activity of the ghrelin receptor is dependent upon flexibility in the C-terminal segment of extracellular loop 2 and that mutations or ligand binding that constrains this segment and thereby conceivably the movements of transmembrane domain V relative to transmembrane domain III inhibits the high constitutive signaling.

Insights into the activation mechanism of the visual receptor rhodopsin

Recent advances in the structural biology of GPCRs (G-protein-coupled receptors) have provided insights into their structure and function. Comparisons of the visual and ligand-activated receptors highlight the unique elements of rhodopsin that allow it to function as a highly sensitive dim-light photoreceptor in vertebrates, as well as the common elements that it shares with the large class A GPCR family. However, despite progress, a number of questions remain unanswered about how these receptors are activated.

The significance of G protein-coupled receptor crystallography for drug discovery

Crucial as molecular sensors for many vital physiological processes, seven-transmembrane domain G protein-coupled receptors (GPCRs) comprise the largest family of proteins targeted by drug discovery. Together with structures of the prototypical GPCR rhodopsin, solved structures of other liganded GPCRs promise to provide insights into the structural basis of the superfamily's biochemical functions and assist in the development of new therapeutic modalities and drugs. One of the greatest technical and theoretical challenges to elucidating and exploiting structure-function relationships in these systems is the emerging concept of GPCR conformational flexibility and its cause-effect relationship for receptor-receptor and receptor-effector interactions. Such conformational changes can be subtle and triggered by relatively small binding energy effects, leading to full or partial efficacy in the activation or inactivation of the receptor system at large. Pharmacological dogma generally dictates that these changes manifest themselves through kinetic modulation of the receptor's G protein partners. Atomic resolution information derived from increasingly available receptor structures provides an entrée to the understanding of these events and practically applying it to drug design. Supported by structure-activity relationship information arising from empirical screening, a unified structural model of GPCR activation/inactivation promises to both accelerate drug discovery in this field and improve our fundamental understanding of structure-based drug design in general. This review discusses fundamental problems that persist in drug design and GPCR structural determination.

Structural transitions of transmembrane helix 6 in the formation of metarhodopsin I

Absorption of light by the visual pigment rhodopsin triggers a rapid cis-trans photoisomerization of its retinal chromophore and a series of conformational changes in both the retinal and protein. The largest structural change is an outward tilt of transmembrane helix H6 that increases the separation of the intracellular ends of H6 and H3 and opens up the G-protein binding site. In the dark state of rhodopsin, Glu247 at the intracellular end of H6 forms a salt bridge with Arg135 on H3 to tether H6 in an inactive conformation. The Arg135-Glu247 interaction is broken in the active state of the receptor, and Arg135 is then stabilized by interactions with Tyr223, Met257, and Tyr306 on helices H5, H6, and H7, respectively. To address the mechanism of H6 motion, solid-state NMR measurements are undertaken of Metarhodopsin I (Meta I), the intermediate preceding the active Metarhodopsin II (Meta II) state of the receptor. (13)C NMR dipolar recoupling measurements reveal an interhelical contact of (13)Cζ-Arg135 with (13)Cε-Met257 in Meta I but not with (13)Cζ-Tyr223 or (13)Cζ-Tyr306. These observations suggest that helix H6 has rotated in the formation of Meta I but that structural changes involving helices H5 and H7 have not yet occurred. Together, our results provide insights into the sequence of events leading up to the outward motion of H6, a hallmark of G protein-coupled receptor activation.

Structural intermediates during α-synuclein fibrillogenesis on phospholipid vesicles

α-Synuclein (AS) fibrils are the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease and other related disorders. AS forms helices that bind phospholipid membranes with high affinity, but no atomic level data for AS aggregation in the presence of lipids is yet available. Here, we present direct evidence of a conversion from α-helical conformation to β-sheet fibrils in the presence of anionic phospholipid vesicles and direct conversion to β-sheet fibrils in their absence. We have trapped intermediate states throughout the fibril formation pathways to examine the structural changes using solid-state NMR spectroscopy and electron microscopy. The comparison between mature AS fibrils formed in aqueous buffer and those derived in the presence of anionic phospholipids demonstrates no major changes in the overall fibril fold. However, a site-specific comparison of these fibrillar states demonstrates major perturbations in the N-terminal domain with a partial disruption of the long β-strand located in the 40s and small perturbations in residues located in the "non-β amyloid component" (NAC) domain. Combining all these results, we propose a model for AS fibrillogenesis in the presence of phospholipid vesicles.

Identification of essential cannabinoid-binding domains: structural insights into early dynamic events in receptor activation

The classical cannabinoid agonist HU210, a structural analog of (-)-Δ(9)-tetrahydrocannabinol, binds to brain cannabinoid (CB1) receptors and activates signal transduction pathways. To date, an exact molecular description of the CB1 receptor is not yet available. Utilizing the minor binding pocket of the CB1 receptor as the primary ligand interaction site, we explored HU210 binding using lipid bilayer molecular dynamics (MD) simulations. Among the potential ligand contact residues, we identified residues Phe-174(2.61), Phe-177(2.64), Leu-193(3.29), and Met-363(6.55) as being critical for HU210 binding by mutational analysis. Using these residues to guide the simulations, we determined essential cannabinoid-binding domains in the CB1 receptor, including the highly sought after hydrophobic pocket important for the binding of the C3 alkyl chain of classical and nonclassical cannabinoids. Analyzing the simulations of the HU210-CB1 receptor complex, the CP55940-CB1 receptor complex, and the (-)-Δ(9)-tetrahydrocannabinol-CB1 receptor complex, we found that the positioning of the C3 alkyl chain and the aromatic stacking between Trp-356(6.48) and Trp-279(5.43) is crucial for the Trp-356(6.48) rotamer change toward receptor activation through the rigid-body movement of H6. The functional data for the mutant receptors demonstrated reductions in potency for G protein activation similar to the reductions seen in ligand binding affinity for HU210.

The roles of transmembrane domain helix-III during rhodopsin photoactivation

BACKGROUND

Rhodopsin, the prototypic member of G protein-coupled receptors (GPCRs), undergoes isomerization of 11-cis-retinal to all-trans-retinal upon photoactivation. Although the basic mechanism by which rhodopsin is activated is well understood, the roles of whole transmembrane (TM) helix-III during rhodopsin photoactivation in detail are not completely clear.

PRINCIPAL FINDINGS

We herein use single-cysteine mutagenesis technique to investigate conformational changes in TM helices of rhodopsin upon photoactivation. Specifically, we study changes in accessibility and reactivity of cysteine residues introduced into the TM helix-III of rhodopsin. Twenty-eight single-cysteine mutants of rhodopsin (P107C-R135C) were prepared after substitution of all natural cysteine residues (C140/C167/C185/C222/C264/C316) by alanine. The cysteine mutants were expressed in COS-1 cells and rhodopsin was purified after regeneration with 11-cis-retinal. Cysteine accessibility in these mutants was monitored by reaction with 4, 4'-dithiodipyridine (4-PDS) in the dark and after illumination. Most of the mutants except for T108C, G109C, E113C, I133C, and R135C showed no reaction in the dark. Wide variation in reactivity was observed among cysteines at different positions in the sequence 108-135 after photoactivation. In particular, cysteines at position 115, 119, 121, 129, 131, 132, and 135, facing 11-cis-retinal, reacted with 4-PDS faster than neighboring amino acids. The different reaction rates of mutants with 4-PDS after photoactivation suggest that the amino acids in different positions in helix-III are exposed to aqueous environment to varying degrees.

SIGNIFICANCE

Accessibility data indicate that an aqueous/hydrophobic boundary in helix-III is near G109 and I133. The lack of reactivity in the dark and the accessibility of cysteine after photoactivation indicate an increase of water/4-PDS accessibility for certain cysteine-mutants at Helix-III during formation of Meta II. We conclude that photoactivation resulted in water-accessible at the chromophore-facing residues of Helix-III.

Light activation of the isomerization and deprotonation of the protonated Schiff base retinal

We perform an ab initio analysis of the photoisomerization of the protonated Schiff base of retinal (PSB-retinal) from 11-cis to 11-trans rotating the C10-C11=C12-C13 dihedral angle from 0° (cis) to -180° (trans). We find that the retinal molecule shows the lowest rotational barrier (0.22 eV) when its charge state is zero as compared to the barrier for the protonated molecule which is ∼0.89 eV. We conclude that rotation most likely takes place in the excited state of the deprotonated retinal. The addition of a proton creates a much larger barrier implying a switching behavior of retinal that might be useful for several applications in molecular electronics. All conformations of the retinal compound absorb in the green region with small shifts following the dihedral angle rotation; however, the Schiff base of retinal (SB-retinal) at trans-conformation absorbs in the violet region. The rotation of the dihedral angle around the C11=C12 π-bond affects the absorption energy of the retinal and the binding energy of the SB-retinal with the proton at the N-Schiff; the binding energy is slightly lower at the trans-SB-retinal than at other conformations of the retinal.

Identification of residue-to-residue contact between a peptide ligand and its G protein-coupled receptor using periodate-mediated dihydroxyphenylalanine cross-linking and mass spectrometry

Fundamental knowledge about how G protein-coupled receptors and their ligands interact is important for understanding receptor-ligand binding and the development of new drug discovery strategies. We have used cross-linking and tandem mass spectrometry analyses to investigate the interaction of the N terminus of the Saccharomyces cerevisiae tridecapeptide pheromone, α-factor (WHWLQLKPGQPMY), and Ste2p, its cognate G protein-coupled receptor. The Trp(1) residue of α-factor was replaced by 3,4-dihydroxyphenylalanine (DOPA) for periodate-mediated chemical cross-linking, and biotin was conjugated to Lys(7) for detection purposes to create the peptide [DOPA(1),Lys(7)(BioACA),Nle(12)]α-factor, called Bio-DOPA(1)-α-factor. This ligand analog was a potent agonist and bound to Ste2p with ∼65 nanomolar affinity. Immunoblot analysis of purified Ste2p samples that were treated with Bio-DOPA(1)-α-factor showed that the peptide analog cross-linked efficiently to Ste2p. The cross-linking was inhibited by the presence of either native α-factor or an α-factor antagonist. MALDI-TOF and immunoblot analyses revealed that Bio-DOPA(1)-α-factor cross-linked to a fragment of Ste2p encompassing residues Ser(251)-Met(294). Fragmentation of the cross-linked fragment and Ste2p using tandem mass spectrometry pinpointed the cross-link point of the DOPA(1) of the α-factor analog to the Ste2p Lys(269) side chain near the extracellular surface of the TM6-TM7 bundle. This conclusion was confirmed by a greatly diminished cross-linking of Bio-DOPA(1)-α-factor into a Ste2p(K269A) mutant. Based on these and previously obtained binding contact data, a mechanism of α-factor binding to Ste2p is proposed. The model for bound α-factor shows how ligand binding leads to conformational changes resulting in receptor activation of the signal transduction pathway.

Putative active states of a prototypic g-protein-coupled receptor from biased molecular dynamics

A major current focus of structural work on G-protein-coupled receptors (GPCRs) pertains to the investigation of their active states. However, for virtually all GPCRs, active agonist-bound intermediate states have been difficult to characterize experimentally owing to their higher conformational flexibility, and thus intrinsic instability, as compared to inactive inverse agonist-bound states. In this work, we explored possible activation pathways of the prototypic GPCR bovine rhodopsin by means of biased molecular dynamics simulations. Specifically, we used an explicit atomistic representation of the receptor and its environment, and sampled the conformational transition from the crystal structure of a photoactivated deprotonated state of rhodopsin to the low pH crystal structure of opsin in the presence of 11-trans-retinal, using adiabatic biased molecular dynamics simulations. We then reconstructed the system free-energy landscape along the predetermined transition trajectories using a path collective variable approach based on metadynamics. Our results suggest that the two experimental endpoints of rhodopsin/opsin are connected by at least two different pathways, and that the conformational transition is populated by at least four metastable states of the receptor, characterized by a different amplitude of the outward movement of transmembrane helix 6.

What can solid state NMR contribute to our understanding of protein folding?

Complete understanding of the folding process that connects a structurally disordered state of a protein to an ordered, biochemically functional state requires detailed characterization of intermediate structural states with high resolution and site specificity. While the intrinsically inhomogeneous and dynamic nature of unfolded and partially folded states limits the efficacy of traditional X-ray diffraction and solution NMR in structural studies, solid state NMR methods applied to frozen solutions can circumvent the complications due to molecular motions and conformational exchange encountered in unfolded and partially folded states. Moreover, solid state NMR methods can provide both qualitative and quantitative structural information at the site-specific level, even in the presence of structural inhomogeneity. This article reviews relevant solid state NMR methods and their initial applications to protein folding studies. Using either chemical denaturation to prepare unfolded states at equilibrium or a rapid freezing apparatus to trap non-equilibrium, transient structural states on a sub-millisecond time scale, recent results demonstrate that solid state NMR can contribute essential information about folding processes that is not available from more familiar biophysical methods.

Transmembrane signaling by GPCRs: insight from rhodopsin and opsin structures

G-protein-coupled receptors (GPCRs), also known as seven-transmembrane (7TM) receptors, are the largest family of membrane proteins in the human genome. As versatile signaling molecules, they mediate cellular responses to extracellular signals. Diffusible ligands like hormones and neurotransmitters bind to GPCRs to modulate GPCR activity. An extraordinary and highly specialized GPCR is the photoreceptor rhodopsin which contains the chromophore retinal as its covalently bound ligand. For receptor activation the configuration of retinal is altered by photon absorption. To date, rhodopsin is the only GPCR for which crystal structures of inactive, active and ligand-free conformations are known. Although the photochemical activation is unique to rhodopsin, many mechanistic insights from this receptor can be generalized for GPCRs.

The properties of thermostabilised G protein-coupled receptors (StaRs) and their use in drug discovery

G protein-coupled receptors (GPCRs) are one of the most important target classes in the central nervous system (CNS) drug discovery, however the fact they are integral membrane proteins and are unstable when purified out of the cell precludes them from a wide range of structural and biophysical techniques that are used for soluble proteins. In this study we demonstrate how protein engineering methods can be used to identify mutations which can both increase the thermostability of receptors, when purified in detergent, as well as biasing the receptor towards a specific physiologically relevant conformational state. We demonstrate this method for the adenosine A(2A) receptor and muscarinic M(1) receptor. The resultant stabilised receptors (known as StaRs) have a pharmacological profile consistent with the inverse agonist conformation. The stabilised receptors can be purified in large quantities, whilst retaining correct folding, thus generating reagents suitable for a broad range of structural and biophysical studies. In the case of the A(2A)-StaR we demonstrate that surface plasmon resonance can be used to profile the association and dissociation rates of a range of antagonists, a technique that can be used to improve the in vivo efficacy of receptor antagonists.

Sequential rearrangement of interhelical networks upon rhodopsin activation in membranes: the Meta II(a) conformational substate

Photon absorption by rhodopsin is proposed to lead to an activation pathway that is described by the extended reaction scheme Meta I <==>Meta II(a) <==> Meta II(b) <==> Meta II(b)H(+), where Meta II(b)H(+) is thought to be the conformational substate that activates the G protein transducin. Here we test this extended scheme for rhodopsin in a membrane bilayer environment by investigating lipid perturbation of the activation mechanism. We found that symmetric membrane lipids having two unsaturated acyl chains, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), selectively stabilize the Meta II(a) substate in the above mechanism. By combining FTIR and UV-visible difference spectroscopy, we characterized the structural and functional changes involved in the transition to the Meta II(a) intermediate, which links the inactive Meta I intermediate with the Meta II(b) states formed by helix rearrangement. Besides the opening of the Schiff base ionic lock, the Meta II(a) substate is characterized by an activation switch in a conserved water-mediated hydrogen-bonded network involving transmembrane helices H1/H2/H7, which is sensed by its key residue Asp83. On the other hand, movement of retinal toward H5 and its interaction with another interhelical H3/H5 network mediated by His211 and Glu122 is absent in Meta II(a). The latter rearrangement takes place only in the subsequent transition to Meta II(b), which has been previously associated with movement of H6. Our results imply that activating structural changes in the H1/H2/H7 network are triggered by disruption of the Schiff base salt bridge and occur prior to other chromophore-induced changes in the H3/H5 network and the outward tilt of H6 in the activation process.