Structure and function in rhodopsin: topology of the C-terminal polypeptide chain in relation to the cytoplasmic loops

Conformational dynamics of activation for the pentameric complex of dimeric G protein-coupled receptor and heterotrimeric G protein

Photoactivation of rhodopsin (Rho), a G protein-coupled receptor, causes conformational changes that provide a specific binding site for the rod G protein, G(t). In this work we employed structural mass spectrometry techniques to elucidate the structural changes accompanying transition of ground state Rho to photoactivated Rho (Rho(∗)) and in the pentameric complex between dimeric Rho(∗) and heterotrimeric G(t). Observed differences in hydroxyl radical labeling and deuterium uptake between Rho(∗) and the (Rho(∗))(2)-G(t) complex suggest that photoactivation causes structural relaxation of Rho following its initial tightening upon G(t) coupling. In contrast, nucleotide-free G(t) in the complex is significantly more accessible to deuterium uptake allowing it to accept GTP and mediating complex dissociation. Thus, we provide direct evidence that in the critical step of signal amplification, Rho(∗) and G(t) exhibit dissimilar conformational changes when they are coupled in the (Rho(∗))(2)-G(t) complex.

Identification of motions in membrane proteins by elastic network models and their experimental validation

Identifying the functional motions of membrane proteins is difficult because they range from large-scale collective dynamics to local small atomic fluctuations at different timescales that are difficult to measure experimentally due to the hydrophobic nature of these proteins. Elastic Network Models, and in particular their most widely used implementation, the Anisotropic Network Model (ANM), have proven to be useful computational methods in many recent applications to predict membrane protein dynamics. These models are based on the premise that biomolecules possess intrinsic mechanical characteristics uniquely defined by their particular architectures. In the ANM, interactions between residues in close proximity are represented by harmonic potentials with a uniform spring constant. The slow mode shapes generated by the ANM provide valuable information on the global dynamics of biomolecules that are relevant to their function. In its recent extension in the form of ANM-guided molecular dynamics (MD), this coarse-grained approach is augmented with atomic detail. The results from ANM and its extensions can be used to guide experiments and thus speedup the process of quantifying motions in membrane proteins. Testing the predictions can be accomplished through (a) direct observation of motions through studies of structure and biophysical probes, (b) perturbation of the motions by, e.g., cross-linking or site-directed mutagenesis, and (c) by studying the effects of such perturbations on protein function, typically through ligand binding and activity assays. To illustrate the applicability of the combined computational ANM-experimental testing framework to membrane proteins, we describe-alongside the general protocols-here the application of ANM to rhodopsin, a prototypical member of the pharmacologically relevant G-protein coupled receptor family.

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.

Generating properly weighted ensemble of conformations of proteins from sparse or indirect distance constraints

Inferring three-dimensional structural information of biomacromolecules such as proteins from limited experimental data is an important and challenging task. Nuclear Overhauser effect measurements based on nucleic magnetic resonance, disulfide linking, and electron paramagnetic resonance labeling studies can all provide useful partial distance constraint characteristic of the conformations of proteins. In this study, we describe a general approach for reconstructing conformations of biomolecules that are consistent with given distance constraints. Such constraints can be in the form of upper bounds and lower bounds of distances between residue pairs, contact maps based on specific contact distance cutoff values, or indirect distance constraints such as experimental phi-value measurement. Our approach is based on the framework of sequential Monte Carlo method, a chain growth-based method. We have developed a novel growth potential function to guide the generation of conformations that satisfy given distance constraints. This potential function incorporates not only the distance information of current residue during growth but also the distance information of future residues by introducing global distance upper bounds between residue pairs and the placement of reference points. To obtain protein conformations from indirect distance constraints in the form of experimental phi-values, we first generate properly weighted contact maps satisfying phi-value constraints, we then generate conformations from these contact maps. We show that our approach can faithfully generate conformations that satisfy the given constraints, which approach the native structures when distance constraints for all residue pairs are given.

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(167) on H4 and Phe(203) 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.

Conformational changes in G-protein-coupled receptors-the quest for functionally selective conformations is open

The G-protein-coupled receptors (GPCRs) represent one the largest families of drug targets. Upon agonist binding a receptor undergoes conformational rearrangements that lead to a novel protein conformation which in turn can interact with effector proteins. During the last decade significant progress has been made to prove that different conformational changes occur. Today it is mostly accepted that individual ligands can induce different receptor conformations. However, the nature or molecular identity of the different conformations is still ill-known. Knowledge of the potential functionally selective conformations will help to develop drugs that select specific conformations of a given GPCR which couple to specific signalling pathways and may, ultimately, lead to reduced side effects. In this review we will summarize recent progress in biophysical approaches that have led to the current understanding of conformational changes that occur during GPCR activation.

Modulation of the absorption maximum of rhodopsin by amino acids in the C-terminus

Vision begins when light is absorbed by visual pigments. It is commonly believed that the absorption spectra of visual pigments are modulated by interactions between the retinal and amino acids within or near 4.5 angstroms of the retinal in the transmembrane (TM) segments. However, this dogma has not been rigorously tested. In this study, we show that the retinal-opsin interactions extend well beyond the retinal binding pocket. We found that, although it is positioned outside of TM segments, the C-terminus of the rhodopsin in the rockfish longspine thornyhead (Sebastolobus altivelis) modulates its lambda(max) by interacting mainly with the last TM segment. Our results illustrate how amino acids in the C-terminus are likely to interact with the retinal. We anticipate our analyses to be a starting point for viewing the spectral tuning of visual pigments as interactions between the retinal and key amino acids that are distributed throughout the entire pigment.

Predisposition of the dark state of rhodopsin to functional changes in structure

As the only member of the family of G-protein-coupled receptors for which atomic coordinates are available, rhodopsin is widely studied for insight into the molecular mechanism of G-protein-coupled receptor activation. The currently available structures refer to the inactive, dark state, of rhodopsin, rather than the light-activated metarhodopsin II (Meta II) state. A model for the Meta II state is proposed here by analyzing elastic network normal modes in conjunction with experimental data. Key mechanical features and interactions broken/formed in the proposed model are found to be consistent with the experimental data. The model is further tested by using a set of Meta II fluorescence decay rates measured to empirically characterize the deactivation of rhodopsin mutants. The model is found to correctly predict 93% of the experimentally observed effects in 119 rhodopsin mutants for which the decay rates and misfolding data have been measured, including a systematic analysis of Cys-->Ser replacements reported here. Based on the detailed comparison between model and experiments, a cooperative activation mechanism is deduced that couples retinal isomerization to concerted changes in conformation, facilitated by the intrinsic dynamics of rhodopsin. A global hinge site is identified near the retinal-binding pocket that ensures the efficient propagation of signals from the central transmembrane region to both cytoplasmic and extracellular ends. The predicted activation mechanism opens the transmembrane helices at the critical G-protein binding cytoplasmic domain. This model provides a detailed, mechanistic description of the activation process, extending experimental observations and yielding new insights for further tests.

Dynamics of the nitroxide side chain in spin-labeled proteins

The dynamics of the tether linking methanethiosulfonate (MTSSL) spin probes to alpha-helices has been investigated with the purpose of rationalizing its effects on ESR line shapes. Torsional profiles for the chain bonds have been calculated ab initio, and steric interactions with the alpha-helix and the neighboring residues have been introduced at the excluded-volume level. As a consequence of the restrictions deriving from chain geometry and local constraints, a limited number of allowed conformers has been identified that undergo torsional oscillations and conformational jumps. Torsional fluctuations are described as damped oscillations, while transition rates between conformers are calculated according to the Langer multidimensional extension of the Kramers theory. The time scale and amplitude of the different motions are compared; the major role played by rotations of the outermost bonds of the side chain emerges, along with the effects of substituents in the pyrroline ring on the conformer distribution and dynamics. The extent and symmetry of magnetic tensor averaging produced by the side chain motions are estimated, the implications for the ESR spectra of spin-labeled proteins are discussed, and suggestions for the introduction of realistic features of the spin probe dynamics into the line shape simulation are presented.

G-protein coupled receptor structure

Because of their central role in regulation of cellular function, structure/function relationships for G-protein coupled receptors (GPCR) are of vital importance, yet only recently have sufficient data been obtained to begin mapping those relationships. GPCRs regulate a wide range of cellular processes, including the senses of taste, smell, and vision, and control a myriad of intracellular signaling systems in response to external stimuli. Many diseases are linked to GPCRs. A critical need exists for structural information to inform studies on mechanism of receptor action and regulation. X-ray crystal structures of only one GPCR, in an inactive state, have been obtained to date. However considerable structural information for a variety of GPCRs has been obtained using non-crystallographic approaches. This review begins with a review of the very earliest GPCR structural information, mostly derived from rhodopsin. Because of the difficulty in crystallizing GPCRs for X-ray crystallography, the extensive published work utilizing alternative approaches to GPCR structure is reviewed, including determination of three-dimensional structure from sparse constraints. The available X-ray crystallographic analyses on bovine rhodopsin are reviewed as the only available high-resolution structures for any GPCR. Structural information available on ligand binding to several receptors is included. The limited information on excited states of receptors is also reviewed. It is concluded that while considerable basic structural information has been obtained, more data are needed to describe the molecular mechanism of activation of a GPCR.

Rhodopsin-transducin coupling: role of the Galpha C-terminus in nucleotide exchange catalysis

In the early steps of visual signal transduction, light-activated rhodopsin (R*) catalyzes GDP/GTP exchange in the heterotrimeric G protein (Galphabetagamma) transducin. We recently reported that the catalytic interaction involves two sequential steps. An initial docking between R* and Gbetagamma leads to conformational changes which make the C-terminus of Galpha (CTalpha) available for binding to R*. Binding of CTalpha by R* then triggers GDP/GTP exchange in the Galpha subunit. To further study this two-step mechanism, we investigated different single amino acid substitutions within CTalpha and discuss the effects of high affinity mutations on nucleotide exchange catalysis.

Structure and dynamics of dark-state bovine rhodopsin revealed by chemical cross-linking and high-resolution mass spectrometry

Recent work using chemical cross-linking to define interresidue distance constraints in proteins has shown that these constraints are useful for testing tertiary structural models. We applied this approach to the G-protein-coupled receptor bovine rhodopsin in its native membrane using lysine- and cysteine-targeted bifunctional cross-linking reagents. Cross-linked proteolytic peptides of rhodopsin were identified by combined liquid chromatography and FT-ICR mass spectrometry with automated data-reduction and assignment software. Tandem mass spectrometry was used to verify cross-link assignments and locate the exact sites of cross-link attachment. Cross-links were observed to form between 10 pairs of residues in dark-state rhodopsin. For each pair, cross-linkers with a range of linker lengths were tested to determine an experimental distance-of-closest-approach (DCA) between reactive side-chain atoms. In all, 28 cross-links were identified using seven different cross-linking reagents. Molecular mechanics procedures were applied to published crystal structure data to calculate energetically achievable theoretical DCAs between reactive atoms without altering the position of the protein backbone. Experimentally measured DCAs are generally in good agreement with the theoretical DCAs. However, a cross-link between C316 and K325 in the C-terminal region cannot be rationalized by DCA simulations and suggests that backbone reorientation relative to the crystal coordinates occurs on the timescale of cross-linking reactions. Biochemical and spectroscopic data from other studies have found that the C-terminal region is highly mobile in solution and not fully represented by X-ray crystallography data. Our results show that chemical cross-linking can provide reliable three-dimensional structural information and insight into local conformational dynamics in a membrane protein.

Molecular motion of spin labeled side chains in the C-terminal domain of RGL2 protein: A SDSL-EPR and MD study

Five singly spin labeled side chains at surface sites in the C-terminal domain of RGL2 protein have been analyzed to investigate the general relationship between nitroxide side chain mobility and protein structure. At these sites, the structural perturbation produced by replacement of a native residue with a nitroxide side chain appears to be very slight at the level of the backbone fold. The primary determinants of the nitroxide side chain mobility are backbone dynamics and tertiary interactions. On the exposed surfaces of alpha-helices, the side chain mobility is not restricted by tertiary interactions but appears to be determined by backbone dynamics, while in loop sites, the side chain mobility is even higher. For a better understanding of the changes in the EPR spectral line shape, molecular dynamics simulations were performed and found in agreement with EPR spectral data.

Molecular mechanism of 7TM receptor activation--a global toggle switch model

The multitude of chemically highly different agonists for 7TM receptors apparently do not share a common binding mode or active site but nevertheless act through induction of a common molecular activation mechanism. A global toggle switch model is proposed for this activation mechanism to reconcile the accumulated biophysical data supporting an outward rigid-body movement of the intracellular segments, as well as the recent data derived from activating metal ion sites and tethered ligands, which suggests an opposite, inward movement of the extracellular segments of the transmembrane helices. According to this model, a vertical see-saw movement of TM-VI-and to some degree TM-VII-around a pivot corresponding to the highly conserved prolines will occur during receptor activation, which may involve the outer segment of TM-V in an as yet unclear fashion. Small-molecule agonists can stabilize such a proposed active conformation, where the extracellular segments of TM-VI and -VII are bent inward toward TM-III, by acting as molecular glue deep in the main ligand-binding pocket between the helices, whereas larger agonists, peptides, and proteins can stabilize a similar active conformation by acting as Velcro at the extracellular ends of the helices and the connecting loops.

Interacting Residues in an Activated State of a G Protein-coupled Receptor

Ste2p, the G protein-coupled receptor (GPCR) for the tridecapeptide pheromone alpha-factor of Saccharomyces cerevisiae, was used as a model GPCR to investigate the role of specific residues in the resting and activated states of the receptor. Using a series of biological and biochemical analyses of wild-type and site-directed mutant receptors, we identified Asn(205) as a potential interacting partner with the Tyr(266) residue. An N205H/Y266H double mutant showed pH-dependent functional activity, whereas the N205H receptor was non-functional and the Y266H receptor was partially active indicating that the histidine 205 and 266 residues interact in an activated state of the receptor. The introduction of N205K or Y266D mutations into the P258L/S259L constitutively active receptor suppressed the constitutive activity; in contrast, the N205K/Y266D/P258L/S259L quadruple mutant was fully constitutively active, again indicating an interaction between residues at the 205 and 206 positions in the receptor-active state. To further test this interaction, we introduced the N205C/Y266C, F204C/Y266C, and N205C/A265C double mutations into wild-type and P258L/S259L constitutively active receptors. After trypsin digestion, we found that a disulfide-cross-linked product, with the molecular weight expected for a receptor fragment with a cross-link between N205C and Y266C, formed only in the N205C/Y266C constitutively activated receptor. This study represents the first experimental demonstration of an interaction between specific residues in an active state, but not the resting state, of Ste2p. The information gained from this study should contribute to an understanding of the conformational differences between resting and active states in GPCRs.

A sequence alignment-independent method for protein classification

Annotation of the rapidly accumulating body of sequence data relies heavily on the detection of remote homologues and functional motifs in protein families. The most popular methods rely on sequence alignment. These include programs that use a scoring matrix to compare the probability of a potential alignment with random chance and programs that use curated multiple alignments to train profile hidden Markov models (HMMs). Related approaches depend on bootstrapping multiple alignments from a single sequence. However, alignment-based programs have limitations. They make the assumption that contiguity is conserved between homologous segments, which may not be true in genetic recombination or horizontal transfer. Alignments also become ambiguous when sequence similarity drops below 40%. This has kindled interest in classification methods that do not rely on alignment. An approach to classification without alignment based on the distribution of contiguous sequences of four amino acids (4-grams) was developed. Interest in 4-grams stemmed from the observation that almost all theoretically possible 4-grams (20(4)) occur in natural sequences and the majority of 4-grams are uniformly distributed. This implies that the probability of finding identical 4-grams by random chance in unrelated sequences is low. A Bayesian probabilistic model was developed to test this hypothesis. For each protein family in Pfam-A and PIR-PSD, a feature vector called a probe was constructed from the set of 4-grams that best characterised the family. In rigorous jackknife tests, unknown sequences from Pfam-A and PIR-PSD were compared with the probes for each family. A classification result was deemed a true positive if the probe match with the highest probability was in first place in a rank-ordered list. This was achieved in 70% of cases. Analysis of false positives suggested that the precision might approach 85% if selected families were clustered into subsets. Case studies indicated that the 4-grams in common between an unknown and the best matching probe correlated with functional motifs from PRINTS. The results showed that remote homologues and functional motifs could be identified from an analysis of 4-gram patterns.

Assessing protein disorder and induced folding

Intrinsically disordered proteins (IDPs) defy the structure-function paradigm as they fulfill essential biological functions while lacking well-defined secondary and tertiary structures. Conformational and spectroscopic analyses showed that IDPs do not constitute a uniform family, and can be divided into subfamilies as a function of their residual structure content. Residual intramolecular interactions are thought to facilitate binding to a partner and then induced folding. Comprehensive information about experimental approaches to investigate structural disorder and induced folding is still scarce. We herein provide hints to readily recognize features typical of intrinsic disorder and review the principal techniques to assess structural disorder and induced folding. We describe their theoretical principles and discuss their respective advantages and limitations. Finally, we point out the necessity of using different approaches and show how information can be broadened by the use of multiples techniques.

Inter- and intra-molecular distances determined by EPR spectroscopy and site-directed spin labeling reveal protein-protein and protein-oligonucleotide interaction

Recent developments including pulse and multi-frequency techniques make the combination of site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy an attractive approach for the study of protein-protein or protein-oligonucleotide interaction. Analysis of the spin label side chain mobility, its solvent accessibility, the polarity of the spin label micro-environment and distances between spin label side chains allow the modeling of protein domains or protein-protein interaction sites and their conformational changes with a spatial resolution at the level of the backbone fold. Structural changes can be detected with millisecond time resolution. Inter- and intra-molecular distances are accessible in the range from approximately 0.5 to 8 nm by the combination of continuous wave and pulse EPR methods. Recent applications include the study of transmembrane substrate transport, membrane channel gating, gene regulation and signal transfer.

Structure of bovine rhodopsin in a trigonal crystal form

We have determined the structure of bovine rhodopsin at 2.65 A resolution using untwinned native crystals in the space group P3(1), by molecular replacement from the 2.8 A model (1F88) solved in space group P4(1). The new structure reveals mechanistically important details unresolved previously, which are considered in the membrane context by docking the structure into a cryo-electron microscopy map of 2D crystals. Kinks in the transmembrane helices facilitate inter-helical polar interactions. Ordered water molecules extend the hydrogen bonding networks, linking Trp265 in the retinal binding pocket to the NPxxY motif near the cytoplasmic boundary, and the Glu113 counterion for the protonated Schiff base to the extracellular surface. Glu113 forms a complex with a water molecule hydrogen bonded between its main chain and side-chain oxygen atoms. This can be expected to stabilise the salt-bridge with the protonated Schiff base linking the 11-cis-retinal to Lys296. The cytoplasmic ends of helices H5 and H6 have been extended by one turn. The G-protein interaction sites mapped to the cytoplasmic ends of H5 and H6 and a spiral extension of H5 are elevated above the bilayer. There is a surface cavity next to the conserved Glu134-Arg135 ion pair. The cytoplasmic loops have the highest temperature factors in the structure, indicative of their flexibility when not interacting with G protein or regulatory proteins. An ordered detergent molecule is seen wrapped around the kink in H6, stabilising the structure around the potential hinge in H6. These findings provide further explanation for the stability of the dark state structure. They support a mechanism for the activation, initiated by photo-isomerisation of the chromophore to its all-trans form, that involves pivoting movements of kinked helices, which, while maintaining hydrophobic contacts in the membrane interior, can be coupled to amplified translation of the helix ends near the membrane surfaces.

Optimal bundling of transmembrane helices using sparse distance constraints

We present a two-step approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only sparse distance constraints, such as those derived from chemical cross-linking, dipolar EPR and FRET experiments. In Step 1, using an algorithm, we developed, the conformational space of membrane protein folds matching a set of distance constraints is explored to provide initial structures for local conformational searches. In Step 2, these structures refined against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. We begin by describing the statistical analysis of the solved membrane protein structures from which the theoretical portion of the penalty function was derived. We then describe the penalty function, and, using a set of six test cases, demonstrate that it is capable of distinguishing helical bundles that are close to the native bundle from those that are far from the native bundle. Finally, using a set of only 27 distance constraints extracted from the literature, we show that our method successfully recovers the structure of dark-adapted rhodopsin to within 3.2 A of the crystal structure.

Electron crystallography reveals the structure of metarhodopsin I

Rhodopsin is the prototypical G protein-coupled receptor, responsible for detection of dim light in vision. Upon absorption of a photon, rhodopsin undergoes structural changes, characterised by distinct photointermediates. Currently, only the ground-state structure has been described. We have determined a density map of a photostationary state highly enriched in metarhodopsin I, to a resolution of 5.5 A in the membrane plane, by electron crystallography. The map shows density for helix 8, the cytoplasmic loops, the extracellular plug, all tryptophan residues, an ordered cholesterol molecule and the beta-ionone ring. Comparison of this map with X-ray structures of the ground state reveals that metarhodopsin I formation does not involve large rigid-body movements of helices, but there is a rearrangement close to the bend of helix 6, at the level of the retinal chromophore. There is no gradual build-up of the large conformational change known to accompany metarhodopsin II formation. The protein remains in a conformation similar to that of the ground state until late in the photobleaching process.

Molecular Dynamics Simulation of Dark-adapted Rhodopsin in an Explicit Membrane Bilayer: Coupling between Local Retinal and Larger Scale Conformational Change

The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature of this change, involving a coupling between a local process and larger changes throughout the protein, may be important for many membrane proteins. In addition, functional work has shown that this coupling occurs with different efficiency in different lipid settings. To begin to understand the nature of the efficiency of this coupling in different lipid settings, we present a molecular dynamics study of rhodopsin in an explicit dioleoyl-phosphatidylcholine bilayer. Our system was simulated for 40 ns and provides insights into the very early events of the visual cascade, before the full transition and activation have occurred. In particular, we see an event near 10 ns that begins with a change in hydrogen bonding near the retinal and that leads through a series of coupled changes to a shift in helical tilt. This type of event, though rare on the molecular dynamics time-scale, could be an important clue to the types of coupling that occur between local and large-scale conformational change in many membrane proteins.

Exploring the conformational space of membrane protein folds matching distance constraints

Herein we present a computational technique for generating helix-membrane protein folds matching a predefined set of distance constraints, such as those obtained from NMR NOE, chemical cross-linking, dipolar EPR, and FRET experiments. The purpose of the technique is to provide initial structures for local conformational searches based on either energetic considerations or ad-hoc scoring criteria. In order to properly screen the conformational space, the technique generates an exhaustive list of conformations within a specified root-mean-square deviation (RMSD) where the helices are positioned in order to match the provided distances. Our results indicate that the number of structures decreases exponentially as the number of distances increases, and increases exponentially as the errors associated with the distances increases. We also found the number of solutions to be smaller when all the distances share one helix in common, compared to the case where the distances connect helices in a daisy-chain manner. We found that for 7 helices, at least 15 distances with errors up to 8 A are needed to produce a number of solutions that is not too large to be processed by local search refinement procedures. Finally, without energetic considerations, our enumeration technique retrieved the transmembrane domains of Bacteriorhodopsin (PDB entry1c3w), Halorhodopsin (1e12), Rhodopsin (1f88), Aquaporin-1 (1fqy), Glycerol uptake facilitator protein (1fx8), Sensory Rhodopsin (1jgj), and a subunit of Fumarate reductase flavoprotein (1qlaC) with Calpha level RMSDs of 3.0 A, 2.3 A, 3.2 A, 4.6 A, 6.0 A, 3.7 A, and 4.4 A, respectively.

Role of the conserved NPxxY(x)5,6F motif in the rhodopsin ground state and during activation

In the G protein-coupled receptor rhodopsin, the conserved NPxxY(x)(5,6)F motif connects the transmembrane helix VII and the cytoplasmic helix 8. The less geometrically constrained retinal analogue 9-demethyl-retinal prevents efficient transformation of rhodopsin to signaling metarhodopsin (Meta) II after retinal photoisomerization. Here, we demonstrate that Ala replacement mutations within the NPxxY(x)(5,6)F domain, which eliminate an interaction between aromatic residues Y306 and F313, allow formation of Meta II despite the presence of 9-demethyl-retinal. Also a disulfide bond linking residues 306 and 313 in the 9-demethyl-retinal-reconstituted mutant Y306C/F313C/C316S prevented Meta II formation, whereas the reduced form of the mutant readily transformed to Meta II after illumination. These observations suggest that the interaction between residues 306 and 313 is disrupted during the Meta I/Meta II transition. However, this enhancement in Meta II formation is not reflected in the G protein activation, which is dramatically reduced for these mutants, suggesting that changes in the Y306-F313 interaction also lead to a proper realigning of helix 8 after photoisomerization. The E134Q mutation, located in the second conserved motif, D(E)RY, rescues activity in 9-demethyl-retinal-reconstituted mutants to different degrees, depending on the position of the Ala replacement in the NPxxY(x)(5,6)F motif, thus revealing distinct roles for the NP and Y(x)(5,6)F portions. Our studies underscore the importance of the NPxxY(x)(5,6)F and D(E)RY motifs in providing structural constraints in rhodopsin that rearrange in response to photoisomerization during formation of the G protein-activating Meta II. The dual control of the structural rearrangements secures reliable transformation of quiescent rhodopsin to activating Meta II.

Structural studies on rhodopsin

Bovine rhodopsin is the prototypical G protein coupled receptor (GPCR). It was the first GPCR to be obtained in quantity and studied in detail. It is also the first GPCR for which detailed three dimensional structural information has been obtained. Reviewed here are the experiments leading up to the high resolution structure determination of rhodopsin and the most recent structural information on the activation and stability of this integral membrane protein.

Conformational changes that occur during M3 muscarinic acetylcholine receptor activation probed by the use of an in situ disulfide cross-linking strategy

The structural changes involved in ligand-dependent activation of G protein-coupled receptors are not well understood at present. To address this issue, we developed an in situ disulfide cross-linking strategy using the rat M(3) muscarinic receptor, a prototypical G(q)-coupled receptor, as a model system. It is known that a tyrosine residue (Tyr(254)) located at the C terminus of transmembrane domain (TM) V and several primarily hydrophobic amino acids present within the cytoplasmic portion of TM VI play key roles in determining the G protein coupling selectivity of the M(3) receptor subtype. To examine whether M3 receptor activation involves changes in the relative orientations of these functionally critical residues, pairs of cysteine residues were substituted into a modified version of the M(3) receptor that contained a factor Xa cleavage site within the third intracellular loop and lacked most endogenous cysteine residues. All analyzed mutant receptors contained a Y254C point mutation and a second cysteine substitution within the segment Lys(484)-Ser(493) at the intracellular end of TM VI. Following their transient expression in COS-7 cells, mutant receptors present in their native membrane environment (in situ) were subjected to mild oxidizing conditions, either in the absence or in the presence of the muscarinic agonist, carbachol. The successful formation of disulfide cross-links was monitored by studying changes in the electrophoretic mobility of oxidized, factor Xa-treated receptors on SDS gels. The observed cross-linking patterns indicated that M(3) receptor activation leads to structural changes that allow the cytoplasmic ends of TM V and TM VI to move closer to each other and that also appear to involve a major change in secondary structure at the cytoplasmic end of TM VI. This is the first study employing an in situ disulfide cross-linking strategy to examine agonist-dependent dynamic structural changes in a G protein-coupled receptor.

Receptor activation: what does the rhodopsin structure tell us?

G-protein-coupled receptors (GPCRs) are a large family of seven-transmembrane-helix proteins that mediate responses to hormones, neurotransmitters and, in the case of rhodopsin, photons. The recent determination of the structure of rhodopsin at atomic resolution opens avenues to a deeper understanding of GPCR activation and transmembrane signaling. Data from previous crosslinking, spin labeling and scanning accessibility experiments on rhodopsin have been mapped onto the high-resolution structure. These data correlate well and are consistent with the structure, and suggest that activation by light opens a cleft at the cytoplasmic end of the seven-helix bundle of rhodopsin. Furthermore, lessons learned from rhodopsin might also apply to other members of this essential family of receptors. (For an animation of the crystal structure of rhodopsin see http://archive.bmn.com/supp/tips/tips2211a.html)

Rhodopsin: structural basis of molecular physiology

The crystal structure of rod cell visual pigment rhodopsin was recently solved at 2.8-A resolution. A critical evaluation of a decade of structure-function studies is now possible. It is also possible to begin to explain the structural basis for several unique physiological properties of the vertebrate visual system, including extremely low dark noise levels as well as high gain and color detection. The ligand-binding pocket of rhodopsin is remarkably compact, and several apparent chromophore-protein interactions were not predicted from extensive mutagenesis or spectroscopic studies. The transmembrane helices are interrupted or kinked at multiple sites. An extensive network of interhelical interactions stabilizes the ground state of the receptor. The helix movement model of receptor activation, which might apply to all G protein-coupled receptors (GPCRs) of the rhodopsin family, is supported by several structural elements that suggest how light-induced conformational changes in the ligand-binding pocket are transmitted to the cytoplasmic surface. The cytoplasmic domain of the receptor is remarkable for a carboxy-terminal helical domain extending from the seventh transmembrane segment parallel to the bilayer surface. Thus the cytoplasmic surface appears to be approximately the right size to bind to the transducin heterotrimer in a one-to-one complex. Future high-resolution structural studies of rhodopsin and other GPCRs will form a basis to elucidate the detailed molecular mechanism of GPCR-mediated signal transduction.

G protein-coupled receptor drug discovery: implications from the crystal structure of rhodopsin

G protein-coupled receptors (GPCRs) are a functionally diverse group of membrane proteins that play a critical role in signal transduction. Because of the lack of a high-resolution structure, the heptahelical transmembrane bundle within the N-terminal extracellular and C-terminal intracellular region of these receptors has initially been modeled based on the high-resolution structure of bacterial retinal-binding protein, bacteriorhodopsin. However, the low-resolution structure of rhodopsin, a prototypical GPCR, revealed that there is a minor relationship between GPCRs and bacteriorhodopsins. The high-resolution crystal structure of the rhodopsin ground state and further refinements of the model provide the first structural information about the entire organization of the polypeptide chain and post-translational moieties. These studies provide a structural template for Family 1 GPCRs that has the potential to significantly improve structure-based approaches to GPCR drug discovery.

Mapping of contact sites in complex formation between transducin and light-activated rhodopsin by covalent crosslinking: use of a photoactivatable reagent

Interaction of light-activated rhodopsin with transducin (T) is the first event in visual signal transduction. We use covalent crosslinking approaches to map the contact sites in interaction between the two proteins. Here we use a photoactivatable reagent, N-[(2-pyridyldithio)-ethyl], 4-azido salicylamide. The reagent is attached to the SH group of cytoplasmic monocysteine rhodopsin mutants by a disulfide-exchange reaction with the pyridylthio group, and the derivatized rhodopsin then is complexed with T by illumination at lambda >495 nm. Subsequent irradiation of the complex at lambda310 nm generates covalent crosslinks between the two proteins. Crosslinking was demonstrated between T and a number of single cysteine rhodopsin mutants. However, sites of crosslinks were investigated in detail only between T and the rhodopsin mutant S240C (cytoplasmic loop V-VI). Crosslinking occurred predominantly with T(alpha). For identification of the sites of crosslinks in T(alpha), the strategy used involved: (i) derivatization of all of the free cysteines in the crosslinked proteins with N-ethylmaleimide; (ii) reduction of the disulfide bond linking the two proteins and isolation of all of the T(alpha) species carrying the crosslinked moiety with a free SH group; (iii) adduct formation of the latter with the N-maleimide moiety of the reagent, maleimido-butyryl-biocytin, containing a biotinyl group; (iv) trypsin degradation of the resulting T(alpha) derivatives and isolation of T(alpha) peptides carrying maleimido-butyryl-biocytin by avidin-agarose chromatography; and (v) identification of the isolated peptides by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. We found that crosslinking occurred mainly to two C-terminal peptides in T(alpha) containing the amino acid sequences 310-313 and 342-345.

NMR spectroscopy in studies of light-induced structural changes in mammalian rhodopsin: applicability of solution (19)F NMR

We report high resolution solution (19)F NMR spectra of fluorine-labeled rhodopsin mutants in detergent micelles. Single cysteine substitution mutants in the cytoplasmic face of rhodopsin were labeled by attachment of the trifluoroethylthio (TET), CF(3)-CH(2)-S, group through a disulfide linkage. TET-labeled cysteine mutants at amino acid positions 67, 140, 245, 248, 311, and 316 in rhodopsin were thus prepared. Purified mutant rhodopsins (6-10 mg), in dodecylmaltoside, were analyzed at 20 degrees C by solution (19)F NMR spectroscopy. The spectra recorded in the dark showed the following chemical shifts relative to trifluoroacetate: Cys-67, 9.8 ppm; Cys-140, 10.6 ppm; Cys-245, 9.9 ppm; Cys-248, 9.5 ppm; Cys-311, 9.9 ppm; and Cys-316, 10.0 ppm. Thus, all mutants showed chemical shifts downfield that of free TET (6.5 ppm). On illumination to form metarhodopsin II, upfield changes in chemical shift were observed for (19)F labels at positions 67 (-0.2 ppm) and 140 (-0.4 ppm) and downfield changes for positions 248 (+0.1 ppm) and 316 (+0.1 ppm) whereas little or no change was observed at positions 311 and 245. On decay of metarhodopsin II, the chemical shifts reverted largely to those originally observed in the dark. The results demonstrate the applicability of solution (19)F NMR spectroscopy to studies of the tertiary structures in the cytoplasmic face of intact rhodopsin in the dark and on light activation.

Nitroxide spin-spin interactions: applications to protein structure and dynamics

Measurement of the distance between two spin label probes in proteins permits the spatial orientation of elements of defined secondary structure. By using site-directed spin labeling, it is possible to determine multiple distance constraints and thereby build tertiary and quaternary structural models as well as measure the kinetics of structural changes. New analytical methods for determining interprobe distances and relative orientations for uniquely oriented spin labels have been developed using global analysis of multifrequency electron paramagnetic resonance data. New methods have also been developed for determining interprobe distances for randomly oriented spin labels. These methods are being applied to a wide range of structural problems, including peptides, soluble proteins, and membrane proteins, that are not readily characterized by other structural techniques.

Use of a disulfide cross-linking strategy to study muscarinic receptor structure and mechanisms of activation

To gain insight into the molecular architecture of the cytoplasmic surface of G protein-coupled receptors, we have developed a disulfide cross-linking strategy using the m3 muscarinic receptor as a model system. To facilitate the interpretation of disulfide cross-linking data, we initially generated a mutant m3 muscarinic receptor (referred to as m3'(3C)-Xa) in which most native Cys residues had been deleted or substituted with Ala or Ser (remaining Cys residues Cys-140, Cys-220, and Cys-532) and in which the central portion of the third intracellular loop had been replaced with a factor Xa cleavage site. Radioligand binding and second messenger assays showed that the m3'(3C)-Xa mutant receptor was fully functional. In the next step, pairs of Cys residues were reintroduced into the m3'(3C)-Xa construct, thus generating 10 double Cys mutant receptors. All 10 mutant receptors contained a Cys residue at position 169 at the beginning of the second intracellular loop and a second Cys within the C-terminal portion of the third intracellular loop, at positions 484-493. Radioligand binding studies and phosphatidylinositol assays indicated that all double Cys mutant receptors were properly folded. Membrane lysates prepared from COS-7 cells transfected with the different mutant receptor constructs were incubated with factor Xa protease and the oxidizing agent Cu(II)-(1,10-phenanthroline)3, and the formation of intramolecular disulfide bonds between juxtaposed Cys residues was monitored by using a combined immunoprecipitation/immunoblotting strategy. To our surprise, efficient disulfide cross-linking was observed with 8 of the 10 double Cys mutant receptors studied (Cys-169/Cys-484 to Cys-491), suggesting that the intracellular m3 receptor surface is characterized by pronounced backbone fluctuations. Moreover, [35S]guanosine 5'-3-O-(thio)triphosphate binding assays indicated that the formation of intramolecular disulfide cross-links prevented or strongly inhibited receptor-mediated G protein activation, suggesting that the highly dynamic character of the cytoplasmic receptor surface is a prerequisite for efficient receptor-G protein interactions. This is the first study using a disulfide mapping strategy to examine the three-dimensional structure of a hormone-activated G protein-coupled receptor.

Structure and function in rhodopsin: further elucidation of the role of the intradiscal cysteines, Cys-110, -185, and -187, in rhodopsin folding and function

The disulfide bond between Cys-110 and Cys-187 in the intradiscal domain is required for correct folding in vivo and function of mammalian rhodopsin. Misfolding in rhodopsin, characterized by the loss of ability to bind 11-cis-retinal, has been shown to be caused by an intradiscal disulfide bond different from the above native disulfide bond. Further, naturally occurring single mutations of the intradiscal cysteines (C110F, C110Y, and C187Y) are associated with retinitis pigmentosa (RP). To elucidate further the role of every one of the three intradiscal cysteines, mutants containing single-cysteine replacements by alanine residues and the above three RP mutants have been studied. We find that C110A, C110F, and C110Y all form a disulfide bond between C185 and C187 and cause loss of retinal binding. C185A allows the formation of a C110-C187 disulfide bond, with wild-type-like rhodopsin phenotype. C187A forms a disulfide bond between C110 and C185 and binds retinal, and the pigment formed has markedly altered bleaching behavior. However, the opsin from the RP mutant C187Y forms no rhodopsin chromophore.

Magic angle spinning NMR of the protonated retinylidene Schiff base nitrogen in rhodopsin: expression of 15N-lysine- and 13C-glycine-labeled opsin in a stable cell line

The apoprotein corresponding to the mammalian photoreceptor rhodopsin has been expressed by using suspension cultures of HEK293S cells in defined media that contained 6-15N-lysine and 2-13C-glycine. Typical yields were 1.5-1.8 mg/liter. Incorporation of 6-15N-lysine was quantitative, whereas that of 2-13C-glycine was about 60%. The rhodopsin pigment formed by binding of 11-cis retinal was spectrally indistinguishable from native bovine rhodopsin. Magic angle spinning (MAS) NMR spectra of labeled rhodopsin were obtained after its incorporation into liposomes. The 15N resonance corresponding to the protonated retinylidene Schiff base nitrogen was observed at 156.8 ppm in the MAS spectrum of 6-15N-lysine-labeled rhodopsin. This chemical shift corresponds to an effective Schiff base-counterion distance of greater than 4 A, consistent with structural water in the binding site hydrogen bonded with the Schiff base nitrogen and the Glu-113 counterion. The present study demonstrates that structural studies of rhodopsin and other G protein-coupled receptors by using MAS NMR are feasible.

Conformational changes in rhodopsin. Movement of helix f detected by site-specific chemical labeling and fluorescence spectroscopy

A recent proposal for the formation of functionally active rhodopsin has placed critical importance on a movement of one of its transmembrane helices (Farrens, D. L., Altenbach, C., Yang, K., Hubbell, W. L., and Khorana, H. G. (1996) Science 274, 768-770). We investigated this hypothesis using a series of eight rhodopsin mutants containing single reactive cysteine residues in the region (helix F) where movement was previously detected. The cysteine mutants were studied in two ways, by measuring their reactivity to a cysteine-specific reagent (PyMPO-maleimide), and by labeling the cysteines with a fluorescent label (monobromobimane) followed by fluorescence spectroscopic analysis. The chemical reactivity data showed sequence-specific variations in reactivity for the mutants in the dark state, resulting in a pattern suggestive of an alpha helix. Interestingly, only upon photoactivation to the MII form did residues found on the inner "face" of this helix react with the PyMPO-maleimide. The ability of the dark state mutants to react with the fluorescent label monobromobimane followed a similar pattern. Furthermore, fluorescence measurements indicate that a bimane label on the inner face of the helix (at V250C) detects changes in the polarity of its environment and accessibility to a fluorescence quenching agent upon MII formation. Viewed together, the data provide further direct evidence that rhodopsin activation involves a conformational change at helix F.

Light-induced exposure of the cytoplasmic end of transmembrane helix seven in rhodopsin

A key step in signal transduction in the visual cell is the light-induced conformational change of rhodopsin that triggers the binding and activation of the guanine nucleotide-binding protein. Site-directed mAbs against bovine rhodopsin were produced and used to detect and characterize these conformational changes upon light activation. Among several antibodies that bound exclusively to the light-activated state, an antibody (IgG subclass) with the highest affinity (Ka approximately 6 x 10(-9) M) was further purified and characterized. The epitope of this antibody was mapped to the amino acid sequence 304-311. This epitope extends from the central region to the cytoplasmic end of the seventh transmembrane helix and incorporates a part of a highly conserved NPXXY motif, a critical region for signaling and agonist-induced internalization of several biogenic amine and peptide receptors. In the dark state, no binding of the antibody to rhodopsin was detected. Accessibility of the epitope to the antibody correlated with formation of the metarhodopsin II photointermediate and was reduced significantly at the metarhodopsin III intermediate. Further, incubation of the antigen-antibody complex with 11-cis-retinal failed to regenerate the native rhodopsin chromophore. These results suggest significant and reversible conformational changes in close proximity to the cytoplasmic end of the seventh transmembrane helix of rhodopsin that might be important for folding and signaling.

Recent advances in site-directed spin labeling of proteins

Site-directed spin labeling of proteins is experiencing a phase of rapid technical evolution, application and evaluation. New strategies have been introduced for determining membrane protein topography, electrostatic potentials, the orientation of proteins at membrane surfaces and inter-residue distances. New applications include studies of beta strands, structure mapping using spin-spin interactions, domain motions in soluble proteins and extensive structural analysis of a number of membrane and soluble proteins.

Structure and function in rhodopsin: rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state

The Glu-134-Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as "molecular sensors" of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.