C5a receptor activation. Genetic identification of critical residues in four transmembrane helices.

Hormones and sensory stimuli activate serpentine receptors, transmembrane switches that relay signals to heterotrimeric guanine nucleotide-binding proteins (G proteins). To understand the switch mechanism, we subjected 93 amino acids in transmembrane helices III, V, VI, and VII of the human chemoattractant C5a receptor to random saturation mutagenesis. A yeast selection identified 121 functioning mutant receptors, containing a total of 523 amino acid substitutions. Conserved hydrophobic residues are located on helix surfaces that face other helices in a modeled seven-helix bundle (Baldwin, J. M., Schertler, G. F., and Unger, V. M. (1997) J. Mol. Biol. 272, 144-164), whereas surfaces predicted to contact the surrounding lipid tolerate many substitutions. Our analysis identified 25 amino acid positions resistant to nonconservative substitutions. These appear to comprise two distinct components of the receptor switch, a surface at or near the extracellular membrane interface and a core cluster in the cytoplasmic half of the bundle. Twenty-one of the 121 mutant receptors exhibit constitutive activity. Amino acids substitutions in these activated receptors predominate in helices III and VI; other activating mutations truncate the receptor near the extracellular end of helix VI. These results identify key elements of a general mechanism for the serpentine receptor switch.     

Conformational change of the transmembrane helices II and IV of metabotropic glutamate receptor involved in G protein activation

G protein-coupled receptors are classified into several families on the basis of their amino acid sequences and the members of the same family exhibit sequence similarity but those of different families do not. In family 1 GPCRs such as rhodopsin and adrenergic receptor, extensive studies have revealed the stimulus-dependent conformational change of the receptor: the rearrangement of transmembrane helices III and VI is essential for G protein activation. In contrast, in family 3 GPCRs such as metabotropic glutamate receptor (mGluR), the inter-protomer relocation upon ligand binding has been observed but there is much less information about the structural changes of the transmsmbrane helices and the cytoplasmic domains. Here we identified constitutively active mutation sites at the cytoplasmic borders of helices II and IV of mGluR8 and successfully inhibited the G protein activation ability by engineering disulfide cross-linking between these cytoplasmic regions. The analysis of all possible single substitution mutants of these residues revealed that some steric interactions around these sites would be important to keep the receptor protein inactive. These results provided the model that the conformational changes at the cytoplasmic ends of helices II and IV of mGluR are involved in the efficient G protein coupling.     

Crystal structure of squid rhodopsin with intracellularly extended cytoplasmic region

G-protein-coupled receptors play a key step in cellular signal transduction cascades by transducing various extracellular signals via G-proteins. Rhodopsin is a prototypical G-protein-coupled receptor involved in the retinal visual signaling cascade. We determined the structure of squid rhodopsin at 3.7A resolution, which transduces signals through the G(q) protein to the phosphoinositol cascade. The structure showed seven transmembrane helices and an amphipathic helix H8 has similar geometry to structures from bovine rhodopsin, coupling to G(t), and human beta(2)-adrenergic receptor, coupling to G(s). Notably, squid rhodopsin contains a well structured cytoplasmic region involved in the interaction with G-proteins, and this region is flexible or disordered in bovine rhodopsin and human beta(2)-adrenergic receptor. The transmembrane helices 5 and 6 are longer and extrude into the cytoplasm. The distal C-terminal tail contains a short hydrophilic alpha-helix CH after the palmitoylated cysteine residues. The residues in the distal C-terminal tail interact with the neighboring residues in the second cytoplasmic loop, the extruded transmembrane helices 5 and 6, and the short helix H8. Additionally, the Tyr-111, Asn-87, and Asn-185 residues are located within hydrogen-bonding distances from the nitrogen atom of the Schiff base.     

Opioid receptor random mutagenesis reveals a mechanism for G protein-coupled receptor activation

The high resolution structure of rhodopsin has greatly enhanced current understanding of G protein-coupled receptor (GPCR) structure in the off-state, but the activation process remains to be clarified. We investigated molecular mechanisms of delta-opioid receptor activation without a preconceived structural hypothesis. Using random mutagenesis of the entire receptor, we identified 30 activating point mutations. Three-dimensional modeling revealed an activation path originating from the third extracellular loop and propagating through tightly packed helices III, VI and VII down to a VI-VII cytoplasmic switch. N- and C-terminal determinants also influence receptor activity. Findings for this therapeutically important receptor may apply to other GPCRs that respond to diffusible ligands.     

Organization of the G protein-coupled receptors rhodopsin and opsin in native membranes

G protein-coupled receptors (GPCRs), which constitute the largest and structurally best conserved family of signaling molecules, are involved in virtually all physiological processes. Crystal structures are available only for the detergent-solubilized light receptor rhodopsin. In addition, this receptor is the only GPCR for which the presumed higher order oligomeric state in native membranes has been demonstrated (Fotiadis, D., Liang, Y., Filipek, S., Saperstein, D. A., Engel, A., and Palczewski, K. (2003) Nature 421, 127-128). Here, we have determined by atomic force microscopy the organization of rhodopsin in native membranes obtained from wild-type mouse photoreceptors and opsin isolated from photoreceptors of Rpe65-/- mutant mice, which do not produce the chromophore 11-cis-retinal. The higher order organization of rhodopsin was present irrespective of the support on which the membranes were adsorbed for imaging. Rhodopsin and opsin form structural dimers that are organized in paracrystalline arrays. The intradimeric contact is likely to involve helices IV and V, whereas contacts mainly between helices I and II and the cytoplasmic loop connecting helices V and VI facilitate the formation of rhodopsin dimer rows. Contacts between rows are on the extracellular side and involve helix I. This is the first semi-empirical model of a higher order structure of a GPCR in native membranes, and it has profound implications for the understanding of how this receptor interacts with partner proteins.     

G protein-coupled receptor activation: analysis of a highly constrained, "straitjacketed" rhodopsin

G protein-coupled receptor (GPCR) activation is generally assumed to result in a significant structural rearrangement of the receptor, presumably involving the rigid body movement of transmembrane helices. We have investigated the activation of the GPCR rhodopsin by the construction and analysis of a mutant which contains a total of four disulfide bonds connecting the cytoplasmic ends of helices 1 and 7, and 3 and 5, and the extracellular ends of helices 3 and 4, and 5 and 6. Despite the constraints imposed by four disulfides, this "straitjacketed" receptor retains the ability to activate the G protein transducin and, therefore, provides insight into the molecular mechanism of the initial step in signal transduction of this important class of receptors.     

Genetic mapping of the human C5a receptor. Identification of transmembrane amino acids critical for receptor function

Many hormones and sensory stimuli signal through a superfamily of seven transmembrane-spanning receptors to activate heterotrimeric G proteins. How the seven transmembrane segments of the receptors (a molecular architecture of bundled alpha-helices conserved from yeast to man) work as "on/off" switches remains unknown. Previously, we used random saturation mutagenesis coupled with a genetic selection in yeast to determine the relative importance of amino acids in four of the seven transmembrane segments of the human C5a receptor (Baranski, T. J., Herzmark, P., Lichtarge, O., Gerber, B. O., Trueheart, J., Meng, E. C., Iiri, T., Sheikh, S. P., and Bourne, H. R. (1999) J. Biol. Chem. 274, 15757-15765). In this study, we evaluate helices I, II, and IV, thereby furnishing a complete mutational map of the seven transmembrane helices of the human C5a receptor. Our analysis identified 19 amino acid positions resistant to non-conservative substitutions. When combined with the 25 essential residues previously identified in helices III and V-VII, they delineate two distinct components of the receptor switch: a ligand-binding surface at or near the extracellular surface of the helix bundle and a core cluster in the cytoplasmic half of the bundle. In addition, we found critical amino acids in the first and second helices that are predicted to face the lipid membrane. These residues form an extended surface that might mediate interactions with lipids and other membrane proteins or function as an oligomerization domain with other receptors.     

Rhodopsin and phototransduction: a model system for G protein-linked receptors.

Rhodopsin is the photoreceptor protein in rod cells of the vertebrate retina and the first member of the class of G protein-coupled receptors for which the amino acid sequence was determined. Rhodopsin is available in greater quantities than any other receptor of its class and therefore has been studied biochemically and biophysically by methods difficult or impossible to apply to its fellow receptors. Such studies support a model in which rhodopsin consists of seven transmembrane helices that form a binding pocket for its ligand, 11-cis retinal. Insights into the structure and function of rhodopsin serve as a model for understanding the structure and function of other members of the receptor class. Rhodopsin undergoes a change in conformation upon photoexcitation and activates a G protein, transducin, and is phosphorylated by a receptor-specific kinase, rhodopsin kinase. The phosphorylated photoactivated rhodopsin is bound by arrestin, thereby terminating activity of the receptor in the signal transduction process. These auxiliary proteins that function with rhodopsin on rod cells serve as models for understanding how other members of the receptor family may function in conjunction with other G proteins, kinases, and arrestin-like proteins.     

Structure and function in rhodopsin: effects of disulfide cross-links in the cytoplasmic face of rhodopsin on transducin activation and phosphorylation by rhodopsin kinase.

Six rhodopsin mutants containing disulfide cross-links between different cytoplasmic regions were prepared: disulfide bond 1, between Cys65 (interhelical loop I-II) and Cys316 (end of helix VII); disulfide bond 2, between Cys246 (end of helix VI) and Cys312 (end of helix VII); disulfide bond 3, between Cys139 (end of helix III) and Cys248 (end of helix VI); disulfide bond 4, between Cys139 (end of helix III) and Cys250 (end of helix VI); disulfide bond 5, between Cys135 (end of helix III) and Cys250 (end of helix VI); and disulfide bond 6, between Cys245 (end of helix VI) and Cys338 (C-terminus). The effects of local restrictions caused by the cross-links on transducin (G(T)) activation and phosphorylation by rhodopsin kinase (RK) following illumination were studied. Disulfide bond 1 showed little effect on either G(T) activation or phosphorylation by RK, suggesting that the relative motion between interhelical loop I-II and helix VII is not crucial for recognition by G(T) or by RK. In contrast, disulfide bonds 2-5 abolished both G(T) activation and phosphorylation by RK. Disulfide bond 6 resulted in enhanced G(T) activation but abolished phosphorylation by RK, suggesting the structure recognized by G(T) was stabilized in this mutant by cross-linking of the C-terminus to the cytoplasmic end of helix VI. Thus, the consequences of the disulfide cross-links depended on the location of the restriction. In particular, relative motions of helix VI, with respect to both helices III and VII upon light activation, are required for recognition of rhodopsin by both G(T) and RK. Further, the conformational changes in the cytoplasmic face that are necessary for protein-protein interactions need not be cooperative, and may be segmental.     

The probable arrangement of the helices in G protein-coupled receptors.

G protein-coupled receptors form a large family of integral membrane proteins whose amino acid sequences have seven hydrophobic segments containing distinctive sequence patterns. Rhodopsin, a member of the family, is known to have transmembrane alpha-helices. The probable arrangement of the seven helices, in all receptors, was deduced from structural information extracted from a detailed analysis of the sequences. Constraints established include: (1) each helix must be positioned next to its neighbours in the sequence; (2) helices I, IV and V must be most exposed to the lipid surrounding the receptor and helix III least exposed. (1) is established from the lengths of the shortest loops. (2) is determined by considering: (i) sites of the most conserved residues; (ii) other sites where variability is restricted; (iii) sites that accommodate polar residues; (iv) sites of differences in sequence between pairs or within groups of closely related receptors. Most sites in the last category should be in unimportant positions and are most useful in determining the position and extent of lipid-facing surface in each helix. The structural constraints for the receptors are used to allocate particular helices to the peaks in the recently published projection map of rhodopsin and to propose a tentative three-dimensional arrangement of the helices in G protein-coupled receptors.     

Repulsive separation of the cytoplasmic ends of transmembrane helices 3 and 6 is linked to receptor activation in a novel thyrotropin receptor mutant (M626I)

Ligand-dependent activation of G protein-coupled receptors (GPCRs) involves repositioning of the juxtacytoplasmic ends of transmembrane helices TM3 and TM6. This concept, inferred from site-directed spin labeling studies, is supported by chemical cross-linking of the cytoplasmic ends of TM3 and TM6 blocking GPCR activation. Here we report a novel constitutive active mutation (M626I) in TM6 of the TSH receptor (TSHR), identified in affected members of a family with nonautoimmune hyperthyroidism. The specific constitutive activity of M626I, measured by its basal cAMP generation corrected for cell surface expression, was 13-fold higher than that of wild-type TSHR. Homology modeling of the TSHR serpentine domain based on the rhodopsin crystal structure suggests that M626 faces the side chain of I515 of TM3 near the membrane-cytoplasmic junction. Steric hindrance of the introduced isoleucine by I515 is consistent with the fact that shorter or more flexible side chains at position 626 did not increase constitutivity. Furthermore, a reciprocal mutation at position 515 (I515M), when introduced into the M626I background, acts as revertant mutation by allowing accommodation of the isoleucine sidechain at position 626 and fully restoring the constitutive activity to the level of wild-type TSHR. Thus, repulsive separation of the juxtacytoplasmic TM6 and TM3 in the M626I model conclusively demonstrates a direct link between the opening of this cytoplasmic face of the receptor structure and G protein coupling.     

An activation switch in the ligand binding pocket of the C5a receptor

Although agonists are thought to occupy binding pockets within the seven-helix core of serpentine receptors, the topography of these binding pockets and the conformational changes responsible for receptor activation are poorly understood. To identify the ligand binding pocket in the receptor for complement factor 5a (C5aR), we assessed binding affinities of hexapeptide ligands, each mutated at a single position, for seven mutant C5aRs, each mutated at a single position in the putative ligand binding site. In ChaW (an antagonist) and W5Cha (an agonist), the side chains at position 5 are tryptophan and cyclohexylalanine, respectively. Comparisons of binding affinities indicated that the hexapeptide residue at this position interacts with two C5aR residues, Ile-116 (helix III) and Val-286 (helix VII); in a C5aR model these two side chains point toward one another. Both the I116A and the V286A mutations markedly increased binding affinity of W5Cha but not that of ChaW. Moreover, ChaW, the antagonist hexapeptide, acted as a full agonist on the I116A mutant. These results argue that C5aR residues Ile-116 and Val-286 interact with the side chain at position 5 of the hexapeptide ligand to form an activation switch. Based on this and previous work, we present a docking model for the hexapeptide within the C5aR binding pocket. We propose that agonists induce a small change in the relative orientations of helices III and VII and that these helices work together to allow movement of helix VI away from the receptor core, thereby triggering G protein activation.     

Use of an in situ disulfide cross-linking strategy to study the dynamic properties of the cytoplasmic end of transmembrane domain VI of the M3 muscarinic acetylcholine receptor

The ligand-induced activation of G protein-coupled receptors (GPCRs) is predicted to involve pronounced conformational changes on the intracellular surface or the receptor proteins. A reorientation of the cytoplasmic end of transmembrane domain VI (TM VI) is thought to play a key role in GPCR activation and productive receptor/G protein coupling. Disulfide cross-linking studies with solubilized, Cys-substituted mutant versions of bovine rhodopsin and the M3 muscarinic acetylcholine receptor suggested that the cytoplasmic end of TM VI is conformationally highly flexible, even in the absence of activating ligands (Farrens, D. L., et al. (1996) Science 274, 768-770; Zeng, F. Y., et al. (1999) J. Biol. Chem. 274, 16629-16640). To test the hypothesis that the promiscuous disulfide cross-linking pattern observed in these studies was caused by the use of solubilized receptor proteins endowed with increased conformational flexibility, we employed a recently developed in situ disulfide cross-linking strategy that allows the detection of disulfide bonds in Cys-substituted mutant M3 muscarinic receptors present in their native membrane environment. Specifically, we used membranes prepared from transfected COS-7 cells to analyze a series of double Cys mutant M3 receptors containing one Cys residue within the sequence K484(6.29) to S493(6.38) at the cytoplasmic end of TM VI and a second Cys residue at the cytoplasmic end of TM III (I169C(3.54)). This analysis revealed a disulfide cross-linking pattern that was strikingly more restricted than that observed previously with solubilized receptor proteins, both in the absence and in the presence of the muscarinic agonist, carbachol. Carbachol stimulated the formation of disulfide bonds in only two of the 10 analyzed mutant muscarinic receptors, I169C(3.54)/K484C(6.29) and I169C(3.54)/A488C(6.33), consistent with an agonist-induced rotation of the cytoplasmic end of TM VI. These findings underline the usefulness of analyzing the structural and dynamic properties of GPCRs in their native lipid environment.     

Structural features and light-dependent changes in the sequence 59-75 connecting helices I and II in rhodopsin: a site-directed spin-labeling study.

Twenty-one single-cysteine substitution mutants were prepared in the sequence 56-75 between transmembrane helices I and II at the cytoplasmic surface of bovine rhodopsin. Each mutant was reacted with a sulfhydryl-specific reagent to produce a nitroxide side chain. The electron paramagnetic resonance of the labeled proteins in dodecyl maltoside solution was analyzed to provide the relative mobility and accessibility of the nitroxide side chain to both polar and nonpolar paramagnetic reagents. The results indicate that the hydrophobic-water interface of the micelle intersects helices I and II near residues 64 and 71, respectively. Thus, the sequence 64-71 is in the aqueous phase, while 56-63 and 72-75 lie in the transmembrane helices I and II, respectively. The lipid-facing surfaces on transmembrane helices I and II near the cytoplasmic surface correspond to approximately 180 degrees and 90 degrees of arc on the helical surfaces, respectively. Photoactivation of rhodopsin produced changes in structure in the region investigated, primarily around helix II. However, these changes are much smaller than those noted by spin labels in helix VI (Altenbach, C., Yang, K., Farrens, D., Farahbakhsh, Z., Khorana, H. G., and Hubbell, W. L. (1996) Biochemistry 35, 12470).     

Selective ligand-induced stabilization of active and desensitized parathyroid hormone type 1 receptor conformations

For many G protein-coupled receptors, agonist-induced activation is followed by desensitization, internalization, and resensitization. In most cases, these processes are dependent upon interaction of agonist-occupied receptor with cytoplasmic beta-arrestins. The ligand-induced intramolecular rearrangements of the receptor responsible for the desensitized versus active conformational states, which dictate both the pharmacological properties of ligands and the biological activity of G protein-coupled receptors, have not been fully elucidated. Here, we identify specific interactions between parathyroid hormone (PTH)-related protein and the human PTH type 1 receptor (PTH1Rc) and the related receptor conformational changes that lead to beta-arrestin-2-mediated desensitization. PTH-related protein analogs modified at position 1 induced selective stabilization of the active G protein-coupled state of the receptor, resulting in lack of beta-arrestin-2 recruitment to the cell membrane, sustained cAMP signaling, and absence of ligand-receptor complex internalization. Mechanistically, the ligands modified at position 1, interacting with the extracellular end of helix VI of PTH1Rc, produced a translocation of transmembrane helices V and VI that differed from that induced by the cognate agonist, resulting in significantly different conformations of the third intracellular loop. These results show how specific interactions between PTH1Rc and its ligands may stabilize distinct conformational states, representing either the active G protein-coupled or a desensitized beta-arrestin-coupled receptor state. In addition, they establish that sustained biological activity of PTH1Rc may be induced by appropriately designed agonist ligands.     

Effects of guanyl nucleotides and rhodopsin on ADP-ribosylation of the inhibitory GTP-binding component of adenylate cyclase by pertussis toxin

Hormonal inhibition of adenylate cyclase is mediated by a guanyl nucleotide binding protein, Gi, which is composed of alpha, beta, and gamma subunits (Gi alpha, G beta gamma). Pertussis toxin blocks hormonal inhibition by catalyzing the ADP-ribosylation of Gi alpha. With purified Gi subunits, but without nucleotides, it was observed that toxin-catalyzed ADP-ribosylation of Gi alpha was negligible in the absence of G beta gamma; ATP, previously shown to increase ADP-ribosylation in membranes, enhanced the ADP-ribosylation of Gi alpha in the absence, more than in the presence, of G beta gamma. Prior studies (Kanaho, Y., Tsai, S.-C., Adamik, R., Hewlett, E.L., Moss, J., and Vaughan, M. (1984) J. Biol. Chem. 259, 7378-7381) had demonstrated that rhodopsin, the retinal photon receptor protein, can replace inhibitory hormone receptors, and stimulate the hydrolysis of GTP by Gi alpha in the presence of G beta gamma. Photolyzed rhodopsin, but not the inactive, dark protein, inhibited ADP-ribosylation of Gi alpha in the presence of G beta gamma. ADP-ribosylation of Gi alpha, in the presence of G beta gamma and photolyzed (but not dark) rhodopsin was increased by guanosine 5'-O-(2-thiodiphosphate) or GDP, but not by (beta, gamma-methylene)guanosine triphosphate or guanosine 5'-O-(3-thiotriphosphate). Presumably, photolyzed rhodopsin and nucleoside triphosphate analogues activate Gi, whereas with dark rhodopsin and nucleoside diphosphates Gi is in the inactive state. The latter appears to be the preferred substrate for pertussis toxin. These observations are consistent with other evidence that rhodopsin and inhibitory hormone receptors are functionally similar.     

Active state of sensory rhodopsin II: structural determinants for signal transfer and proton pumping

The molecular mechanism of transmembrane signal transduction is still a pertinent question in cellular biology. Generally, a receptor can transfer an external signal via its cytoplasmic surface, as found for G-protein-coupled receptors such as rhodopsin, or via the membrane domain, such as that in sensory rhodopsin II (SRII) in complex with its transducer, HtrII. In the absence of HtrII, SRII functions as a proton pump. Here, we report on the crystal structure of the active state of uncomplexed SRII from Natronomonas pharaonis, NpSRII. The problem with a dramatic loss of diffraction quality upon loading of the active state was overcome by growing better crystals and by reducing the occupancy of the state. The conformational changes in the region comprising helices F and G are similar to those observed for the NpSRII-transducer complex but are much more pronounced. The meaning of these differences for the understanding of proton pumping and signal transduction by NpSRII is discussed.     

Control of conformational equilibria in the human B2 bradykinin receptor. Modeling of nonpeptidic ligand action and comparison to the rhodopsin structure

A prototypic study of the molecular mechanisms of activation or inactivation of peptide hormone G protein-coupled receptors was carried out on the human B2 bradykinin receptor. A detailed pharmacological analysis of receptor mutants possessing either increased constitutive activity or impaired activation or ligand recognition allowed us to propose key residues participating in intramolecular interaction networks stabilizing receptor inactive or active conformations: Asn(113) and Tyr(115) (TM III), Trp(256) and Phe(259) (TM VI), Tyr(295) (TM VII) which are homologous of the rhodopsin residues Gly(120), Glu(122), Trp(265), Tyr(268), and Lys(296), respectively. An essential experimental finding was the spatial proximity between Asn(113), which is the cornerstone of inactive conformations, and Trp(256) which plays a subtle role in controlling the balance between active and inactive conformations. Molecular modeling and mutagenesis data showed that Trp(256) and Tyr(295) constitute, together with Gln(288), receptor contact points with original nonpeptidic ligands. It provided an explanation for the ligand inverse agonist behavior on the WT receptor, with underlying restricted motions of TMs III, VI, and VII, and its agonist behavior on the Ala(113) and Phe(256) constitutively activated mutants. These data on the B2 receptor emphasize that conformational equilibria are controlled in a coordinated fashion by key residues which are located at strategic positions for several G protein-coupled receptors. They are discussed in comparison with the recently determined rhodopsin crystallographic structure.     

Agonist-induced conformational changes in thyrotropin-releasing hormone receptor type I: disulfide cross-linking and molecular modeling approaches

The conformational changes at the cytoplasmic ends of transmembrane helices 5 and 6 (TMH5 and TMH6) of thyrotropin-releasing hormone (TRH) receptor type I (TRH-R1) during activation were analyzed by cysteine-scanning mutagenesis followed by disulfide cross-linking and molecular modeling. Sixteen double cysteine mutants were constructed by substitution of one residue at the cytoplasmic end of TMH5 and the other at that of TMH6. The cross-linking experiments indicate that four mutants, Q263C/G212C, Q263C/Y211C, T265C/G212C, and T265C/Y211C, exhibited disulfide bond formation that was sensitive to TRH occupancy. We refined our previous TRH-R1 models by embedding them into a hydrated explicit lipid bilayer. Molecular dynamics simulations of the models, as well as in silico double cysteine models, generated trajectories that were in agreement with experimental results. Our findings suggest that TRH binding induces a separation of the cytoplasmic ends of TMH5 and TMH6 and a rotation of TMH6. These changes likely increase the surface accessible area at the juxtamembrane region of intracellular loop 3 that could promote interactions between G proteins and key residues within the receptor.     

Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes

As one of the best studied members of the pharmaceutically relevant family of G-protein-coupled receptors, rhodopsin serves as a prototype for understanding the mechanism of G-protein-coupled receptor activation. Here, we aim at exploring functionally relevant conformational changes and signal transmission mechanisms involved in its photoactivation brought about through a cis-trans photoisomerization of retinal. For this exploration, we propose a molecular dynamics simulation protocol that utilizes normal modes derived from the anisotropic network model for proteins. Deformations along multiple low-frequency modes of motion are used to efficiently sample collective conformational changes in the presence of explicit membrane and water environment, consistent with interresidue interactions. We identify two highly stable regions in rhodopsin, one clustered near the chromophore, the other near the cytoplasmic ends of transmembrane helices H1, H2, and H7. Due to redistribution of interactions in the neighborhood of retinal upon stabilization of the trans form, local structural rearrangements in the adjoining H3-H6 residues are efficiently propagated to the cytoplasmic end of these particular helices. In the structures obtained by our simulations, all-trans retinal interacts with Cys¹⁶⁷ on H4 and Phe²⁰³ on H5, which were not accessible in the dark state, and exhibits stronger interactions with H5, while some of the contacts made (in the cis form) with H6 are lost.