Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs

G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.     

Modeling GPCR active state conformations: the β(2)-adrenergic receptor

The recent publication of several G protein-coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β(2)-adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site-directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site-directed mutagenesis information. Virtual screening experiments show that the models are selective for β-adrenergic agonists over other GPCR ligands, for (R)- over (S)-β-hydroxy agonists and for β(2)-selective agonists over β(1)-selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C-terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C-terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR-G protein model has been used to explain site-directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β(2)-adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs.     

Relevance of rhodopsin studies for GPCR activation

Rhodopsin, the dim-light photoreceptor present in the rod cells of the retina, is both a retinal-binding protein and a G protein-coupled receptor (GPCR). Due to this conjunction, it benefits from an arsenal of spectroscopy techniques that can be used for its characterization, while being a model system for the important family of Class A (also referred to as “rhodopsin-like”) GPCRs. For instance, rhodopsin has been a crucial player in the field of GPCR structural biology. Until 2007, it was the only GPCR for which a high-resolution crystal structure was available, so all structure–activity analyses on GPCRs, from structure-based drug discovery to studies of structural changes upon activation, were based on rhodopsin. At present, about a third of currently available GPCR structures are still from rhodopsin. In this review, I show some examples of how these structures can still be used to gain insight into general aspects of GPCR activation. First, the analysis of the third intracellular loop in rhodopsin structures allows us to gain an understanding of the structural and dynamic properties of this region, which is absent (due to protein engineering or poor electron density) in most of the currently available GPCR structures. Second, a detailed analysis of the structure of the transmembrane domains in inactive, intermediate and active rhodopsin structures allows us to detect early conformational changes in the process of ligand-induced GPCR activation. Finally, the analysis of a conserved ligand-activated transmission switch in the transmembrane bundle of GPCRs in the context of the rhodopsin activation cycle, allows us to suggest that the structures of many of the currently available agonist-bound GPCRs may correspond to intermediate active states. While the focus in GPCR structural biology is inevitably moving away from rhodopsin, in other aspects rhodopsin is still at the forefront. For instance, the first studies of the structural basis of disease mutants in GPCRs, or the most detailed analysis of cellular GPCR signal transduction networks using a systems biology approach, have been carried out in rhodopsin. Finally, due again to its unique properties among GPCRs, rhodopsin will likely play an important role in the application of X-ray free electron laser crystallography to time-resolved structural biology in membrane proteins. Rhodopsin, thus, still remains relevant as a model system to study the molecular mechanisms of GPCR activation. This article is part of a Special Issue entitled: Retinal Proteins—You can teach an old dog new tricks.     

Structural insights into agonist-induced activation of G-protein-coupled receptors

Recent years have seen tremendous breakthroughs in structure determination of G-protein-coupled receptors (GPCRs). In 2011, two agonist-bound active-state structures of rhodopsin have been published. Together with structures of several rhodopsin activation intermediates and a wealth of biochemical and spectroscopic information, they provide a unique structural framework on which to understand GPCR activation. Here we use this framework to compare the recent crystal structures of the agonist-bound active states of the β(2) adrenergic receptor (β(2)AR) and the A(2A) adenosine receptor (A(2A)AR). While activation of these three GPCRs results in rearrangements of TM5 and TM6, the extent of this conformational change varies considerably. Displacements of the cytoplasmic side of TM6 ranges between 3 and 8? depending on whether selective stabilizers of the active conformation are used (i.e. a G-protein peptide in the case of rhodopsin or a conformationally selective nanobody in the case of the β(2)AR) or not (A(2A)AR). The agonist-induced conformational changes in the ligand-binding pocket are largely receptor specific due to the different chemical nature of the agonists. However, several similarities can be observed, including a relocation of conserved residues W6.48 and F6.44 towards L5.51 and P5.50, and of I/L3.40 away from P5.50. This transmission switch links agonist binding to the movement of TM5 and TM6 through the rearrangement of the TM3-TM5-TM6 interface, and possibly constitutes a common theme of GPCR activation.     

Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors

Activation of G-protein-coupled receptors (GPCRs) is initiated by conformational changes in the transmembrane (TM) helices and the intra- and extracellular loops induced by ligand binding. Understanding the conformational changes in GPCRs leading to activation is imperative in deciphering the role of these receptors in the pathology of diseases. Since the crystal structures of activated GPCRs are not yet available, computational methods and biophysical techniques have been used to predict the structures of GPCR active states. We have recently applied the computational method LITiCon to understand the ligand-induced conformational changes in β₂-adrenergic receptor by ligands of varied efficacies. Here we report a study of the conformational changes associated with the activation of bovine rhodopsin for which the crystal structure of the inactive state is known. Starting from the inactive (dark) state, we have predicted the TM conformational changes that are induced by the isomerization of 11-cis retinal to all-trans retinal leading to the fully activated state, metarhodopsin II. The predicted active state of rhodopsin satisfies all of the 30 known experimental distance constraints. The predicted model also correlates well with the experimentally observed conformational switches in rhodopsin and other class A GPCRs, namely, the breaking of the ionic lock between R135^3.50 at the intracellular end of TM3 (part of the DRY motif) and E247^6.30 on TM6, and the rotamer toggle switch on W265^6.48 on TM6. We observe that the toggling of the W265^6.48 rotamer modulates the bend angle of TM6 around the conserved proline. The rotamer toggling is facilitated by the formation of a water wire connecting S298^7.45, W265^6.48 and H211^5.46. As a result, the intracellular ends of TMs 5 and 6 move outward from the protein core, causing large conformational changes at the cytoplasmic interface. The predicted outward movements of TM5 and TM6 are in agreement with the recently published crystal structure of opsin, which is proposed to be close to the active-state structure. In the predicted active state, several residues in the intracellular loops, such as R69, V139^3.54, T229, Q237, Q239, S240, T243 and V250^6.33, become more water exposed compared to the inactive state. These residues may be involved in mediating the conformational signal from the receptor to the G protein. From mutagenesis studies, some of these residues, such as V139^3.54, T229 and V250^6.33, are already implicated in G-protein activation. The predicted active state also leads to the formation of new stabilizing interhelical hydrogen-bond contacts, such as those between W265^6.48 and H211^5.46 and E122^3.37 and C167^4.56. These hydrogen-bond contacts serve as potential conformational switches offering new opportunities for future experimental investigations. The calculated retinal binding energy surface shows that binding of an agonist makes the receptor dynamic and flexible and accessible to many conformations, while binding of an inverse agonist traps the receptor in the inactive state and makes the other conformations inaccessible.     

Conserved activation pathways in G-protein-coupled receptors

GPCRs (G-protein-coupled receptors) are seven-transmembrane helix proteins that transduce exogenous and endogenous signals to modulate the activity of downstream effectors inside the cell. Despite the relevance of these proteins in human physiology and pharmaceutical research, we only recently started to understand the structural basis of their activation mechanism. In the period 2008-2011, nine active-like structures of GPCRs were solved. Among them, we have determined the structure of light-activated rhodopsin with all the features of the active metarhodopsin-II, which represents so far the most native-like model of an active GPCR. This structure, together with the structures of other inactive, intermediate and active states of rhodopsin constitutes a unique structural framework on which to understand the conserved aspects of the activation mechanism of GPCRs. This mechanism can be summarized as follows: retinal isomerization triggers a series of local structural changes in the binding site that are amplified into three intramolecular activation pathways through TM (transmembrane helix) 5/TM3, TM6 and TM7/TM2. Sequence analysis strongly suggests that these pathways are conserved in other GPCRs. Differential activation of these pathways by ligands could be translated into the stabilization of different active states of the receptor with specific signalling properties.     

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.     

Molecular dynamics simulations reveal insights into key structural elements of adenosine receptors

The crystal structure of the human A(2A) adenosine receptor, a member of the G protein-coupled receptor (GPCR) family, is used as a starting point for the structural characterization of the conformational equilibrium around the inactive conformation of the human A(2) (A(2A) and A(2B)) adenosine receptors (ARs). A homology model of the closely related A(2B)AR is reported, and the two receptors were simulated in their apo form through all-atom molecular dynamics (MD) simulations. Different conditions were additionally explored in the A(2A)AR, including the protonation state of crucial histidines or the presence of the cocrystallized ligand. Our simulations reveal the role of several conserved residues in the ARs in the conformational equilibrium of the receptors. The "ionic lock" absent in the crystal structure of the inactive A(2A)AR is rapidly formed in the two simulated receptors, and a complex network of interacting residues is presented that further stabilizes this structural element. Notably, the observed rotameric transition of Trp6.48 ("toggle switch"), which is thought to initiate the activation process in GPCRs, is accompanied by a concerted rotation of the conserved residue of the A(2)ARs, His6.52. This new conformation is further stabilized in the two receptors under study by a novel interaction network involving residues in transmembrane (TM) helices TM5 (Asn5.42) and TM3 (Gln3.37), which resemble the conformational changes recently observed in the agonist-bound structure of β-adrenoreceptors. Finally, the interaction between Glu1.39 and His7.43, a pair of conserved residues in the family of ARs, is found to be weaker than previously thought, and the role of this interaction in the structure and dynamics of the receptor is thoroughly examined. All these findings suggest that, despite the commonalities with other GPCRs, the conformational equilibrium of ARs is also modulated by specific residues of the family.     

Protein‐based virtual screening of chemical databases. II. Are homology models of g‐protein coupled receptors suitable targets?

The aim of the current study is to investigate whether homology models of G-Protein-Coupled Receptors (GPCRs) that are based on bovine rhodopsin are reliable enough to be used for virtual screening of chemical databases. Starting from the recently described 2.8 A-resolution X-ray structure of bovine rhodopsin, homology models of an "antagonist-bound" form of three human GPCRs (dopamine D3 receptor, muscarinic M1 receptor, vasopressin V1a receptor) were constructed. The homology models were used to screen three-dimensional databases using three different docking programs (Dock, FlexX, Gold) in combination with seven scoring functions (ChemScore, Dock, FlexX, Fresno, Gold, Pmf, Score). Rhodopsin-based homology models turned out to be suitable, indeed, for virtual screening since known antagonists seeded in the test databases could be distinguished from randomly chosen molecules. However, such models are not accurate enough for retrieving known agonists. To generate receptor models better suited for agonist screening, we developed a new knowledge- and pharmacophore-based modeling procedure that might partly simulate the conformational changes occurring in the active site during receptor activation. Receptor coordinates generated by this new procedure are now suitable for agonist screening. We thus propose two alternative strategies for the virtual screening of GPCR ligands, relying on a different set of receptor coordinates (antagonist-bound and agonist-bound states).     

Formation and Decay of the Arrestin·Rhodopsin Complex in Native Disc Membranes

In the G protein-coupled receptor rhodopsin, light-induced cis/trans isomerization of the retinal ligand triggers a series of distinct receptor states culminating in the active Metarhodopsin II (Meta II) state, which binds and activates the G protein transducin (G_t). Long before Meta II decays into the aporeceptor opsin and free all-trans-retinal, its signaling is quenched by receptor phosphorylation and binding of the protein arrestin-1, which blocks further access of G_t to Meta II. Although recent crystal structures of arrestin indicate how it might look in a precomplex with the phosphorylated receptor, the transition into the high affinity complex is not understood. Here we applied Fourier transform infrared spectroscopy to monitor the interaction of arrestin-1 and phosphorylated rhodopsin in native disc membranes. By isolating the unique infrared signature of arrestin binding, we directly observed the structural alterations in both reaction partners. In the high affinity complex, rhodopsin adopts a structure similar to G_t-bound Meta II. In arrestin, a modest loss of β-sheet structure indicates an increase in flexibility but is inconsistent with a large scale structural change. During Meta II decay, the arrestin-rhodopsin stoichiometry shifts from 1:1 to 1:2. Arrestin stabilizes half of the receptor population in a specific Meta II protein conformation, whereas the other half decays to inactive opsin. Altogether these results illustrate the distinct binding modes used by arrestin to interact with different functional forms of the receptor.     

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.     

Differences in Allosteric Communication Pipelines in the Inactive and Active States of a GPCR

G-protein-coupled receptors (GPCRs) are membrane proteins that allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to proteins in the intracellular (IC) domain. However, the complete pathway of allosteric communication from the EC to the IC domain, including the role of individual amino acids in the pathway is not known. Using the correlation in torsion angle movements calculated from microseconds-long molecular-dynamics simulations, we elucidated the allosteric pathways in three different conformational states of β₂-adrenergic receptor (β₂AR): 1), the inverse-agonist-bound inactive state; 2), the agonist-bound intermediate state; and (3), the agonist- and G-protein-bound fully active state. The inactive state is less dynamic compared with the intermediate and active states, showing dense clusters of allosteric pathways (allosteric pipelines) connecting the EC with the IC domain. The allosteric pipelines from the EC domain to the IC domain are weakened in the intermediate state, thus decoupling the EC domain from the IC domain and making the receptor more dynamic compared with the other states. Also, the orthosteric ligand-binding site becomes the initiator region for allosteric communication in the intermediate state. This finding agrees with the paradigm that the nature of the agonist governs the specific signaling state of the receptor. These results provide an understanding of the mechanism of allosteric communication in class A GPCRs. In addition, our analysis shows that mutations that affect the ligand efficacy, but not the binding affinity, are located in the allosteric pipelines. This clarifies the role of such mutations, which has hitherto been unexplained.     

Differences in Allosteric Communication Pipelines in the Inactive and Active States of a GPCR

G-protein-coupled receptors (GPCRs) are membrane proteins that allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to proteins in the intracellular (IC) domain. However, the complete pathway of allosteric communication from the EC to the IC domain, including the role of individual amino acids in the pathway is not known. Using the correlation in torsion angle movements calculated from microseconds-long molecular-dynamics simulations, we elucidated the allosteric pathways in three different conformational states of β₂-adrenergic receptor (β₂AR): 1), the inverse-agonist-bound inactive state; 2), the agonist-bound intermediate state; and (3), the agonist- and G-protein-bound fully active state. The inactive state is less dynamic compared with the intermediate and active states, showing dense clusters of allosteric pathways (allosteric pipelines) connecting the EC with the IC domain. The allosteric pipelines from the EC domain to the IC domain are weakened in the intermediate state, thus decoupling the EC domain from the IC domain and making the receptor more dynamic compared with the other states. Also, the orthosteric ligand-binding site becomes the initiator region for allosteric communication in the intermediate state. This finding agrees with the paradigm that the nature of the agonist governs the specific signaling state of the receptor. These results provide an understanding of the mechanism of allosteric communication in class A GPCRs. In addition, our analysis shows that mutations that affect the ligand efficacy, but not the binding affinity, are located in the allosteric pipelines. This clarifies the role of such mutations, which has hitherto been unexplained.     

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation

Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs)-the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state ((2))H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Structural basis for ligand binding and specificity in adrenergic receptors: implications for GPCR-targeted drug discovery

Crystal structures of engineered human beta 2-adrenergic receptors (ARs) in complex with an inverse agonist ligand, carazolol, provide three-dimensional snapshots of the disposition of seven transmembrane helices and the ligand-binding site of an important G protein-coupled receptor (GPCR). As expected, beta 2-AR shares substantial structural similarities with rhodopsin, the dim-light photoreceptor of the rod cell. However, although carazolol and the 11- cis-retinylidene moiety of rhodopsin are situated in the same general binding pocket, the second extracellular (E2) loop structures are quite distinct. E2 in rhodopsin shows beta-sheet structure and forms part of the chromophore-binding site. In the beta 2-AR, E2 is alpha-helical and seems to be distinct from the receptor's active site, allowing a potential entry pathway for diffusible ligands. The structures, together with extensive structure-activity relationship (SAR) data from earlier studies, provide insight about possible structural determinants of ligand specificity and how the binding of agonist ligands might alter receptor conformation. We review key features of the new beta 2-AR structures in the context of recent complementary work on the conformational dynamics of GPCRs. We also report 600 ns molecular dynamics simulations that quantified beta 2-AR receptor mobility in a membrane bilayer environment and show how the binding of an agonist ligand, adrenaline (epinephrine), causes conformational changes to the ligand-binding pocket and neighboring helices.     

Mechanism of activation of a G protein-coupled receptor, the human cholecystokinin-2 receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.     

Structural biology of G protein‐coupled receptor signaling complexes

G protein‐coupled receptors (GPCRs) constitute the largest family of cell surface receptors that mediate numerous cell signaling pathways, and are targets of more than one‐third of clinical drugs. Thanks to the advancement of novel structural biology technologies, high‐resolution structures of GPCRs in complex with their signaling transducers, including G‐protein and arrestin, have been determined. These 3D complex structures have significantly improved our understanding of the molecular mechanism of GPCR signaling and provided a structural basis for signaling‐biased drug discovery targeting GPCRs. Here we summarize structural studies of GPCR signaling complexes with G protein and arrestin using rhodopsin as a model system, and highlight the key features of GPCR conformational states in biased signaling including the sequence motifs of receptor TM6 that determine selective coupling of G proteins, and the phosphorylation codes of GPCRs for arrestin recruitment. We envision the future of GPCR structural biology not only to solve more high‐resolution complex structures but also to show stepwise GPCR signaling complex assembly and disassembly and dynamic process of GPCR signal transduction.     

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.     

Internal hydration increases during activation of the G-protein-coupled receptor rhodopsin

Rhodopsin, the membrane protein responsible for dim-light vision, until recently was the only G-protein-coupled receptor (GPCR) with a known crystal structure. As a result, there is enormous interest in studying its structure, dynamics, and function. Here we report the results of three all-atom molecular dynamics simulations, each at least 1.5 μs, which predict that substantial changes in internal hydration play a functional role in rhodopsin activation. We confirm with ¹H magic angle spinning NMR that the increased hydration is specific to the metarhodopsin-I intermediate. The internal water molecules interact with several conserved residues, suggesting that changes in internal hydration may be important during the activation of other GPCRs. The results serve to illustrate the synergism of long-time-scale molecular dynamics simulations and NMR in enhancing our understanding of GPCR function.     

A Constitutively Activating Mutation Alters the Dynamics and Energetics of a Key Conformational Change in a Ligand-free G Protein-coupled Receptor

G protein-coupled receptors (GPCRs) undergo dynamic transitions between active and inactive conformations. Usually, these conversions are triggered when the receptor detects an external signal, but some so-called constitutively activating mutations, or CAMs, induce a GPCR to bind and activate G proteins in the absence of external stimulation, in ways still not fully understood. Here, we investigated how a CAM alters the structure of a GPCR and the dynamics involved as the receptor transitions between different conformations. Our approach used site-directed fluorescence labeling (SDFL) spectroscopy to compare opsin, the ligand-free form of the GPCR rhodopsin, with opsin containing the CAM M257Y, focusing specifically on key movements that occur in the sixth transmembrane helix (TM6) during GPCR activation. The site-directed fluorescence labeling data indicate opsin is constrained to an inactive conformation both in detergent micelles and lipid membranes, but when it contains the M257Y CAM, opsin is more dynamic and can interact with a G protein mimetic. Further study of these receptors using tryptophan-induced quenching (TrIQ) methods indicates that in detergent, the CAM significantly increases the population of receptors in the active state, but not in lipids. Subsequent Arrhenius analysis of the TrIQ data suggests that, both in detergent and lipids, the CAM lowers the energy barrier for TM6 movement, a key transition required for conversion between the inactive and active conformations. Together, these data suggest that the lowered energy barrier is a primary effect of the CAM on the receptor dynamics and energetics.