The role of Glu181 in the photoactivation of rhodopsin

Methodology of pulsed photoacoustics and its application to probe photosystems and receptors

We review recent advances in the methodology of pulsed time-resolved photoacoustics and its application to studies of photosynthetic reaction centers and membrane receptors such as the G protein-coupled receptor rhodopsin. The experimental parameters accessible to photoacoustics include molecular volume change and photoreaction enthalpy change. Light-driven volume change secondary to protein conformational changes or electrostriction is directly related to the photoreaction and thus can be a useful measurement of activity and function. The enthalpy changes of the photochemical reactions observed can be measured directly by photoacoustics. With the measurement of enthalpy change, the reaction entropy can also be calculated when free energy is known. Dissecting the free energy of a photoreaction into enthalpic and entropic components may provide critical information about photoactivation mechanisms of photosystems and photoreceptors. The potential limitations and future applications of time-resolved photoacoustics are also discussed.

6-s-cis Conformation and polar binding pocket of the retinal chromophore in the photoactivated state of rhodopsin

The visual pigment rhodopsin is unique among the G protein-coupled receptors in having an 11-cis retinal chromophore covalently bound to the protein through a protonated Schiff base linkage. The chromophore locks the visual receptor in an inactive conformation through specific steric and electrostatic interactions. This efficient inverse agonist is rapidly converted to an agonist, the unprotonated Schiff base of all-trans retinal, upon light activation. Here, we use magic angle spinning NMR spectroscopy to obtain the (13)C chemical shifts (C5-C20) of the all-trans retinylidene chromophore and the (15)N chemical shift of the Schiff base nitrogen in the active metarhodopsin II intermediate. The retinal chemical shifts are sensitive to the conformation of the chromophore and its molecular interactions within the protein-binding site. Comparison of the retinal chemical shifts in metarhodopsin II with those of retinal model compounds reveals that the Schiff base environment is polar. In particular, the (13)C15 and (15)Nepsilon chemical shifts indicate that the C horizontal lineN bond is highly polarized in a manner that would facilitate Schiff base hydrolysis. We show that a strong perturbation of the retinal (13)C12 chemical shift observed in rhodopsin is reduced in wild-type metarhodopsin II and in the E181Q mutant of rhodopsin. On the basis of the T(1) relaxation time of the retinal (13)C18 methyl group and the conjugated retinal (13)C5 and (13)C8 chemical shifts, we have determined that the conformation of the retinal C6-C7 single bond connecting the beta-ionone ring and the retinylidene chain is 6-s-cis in both the inactive and the active states of rhodopsin. These results are discussed within the general framework of ligand-activated G protein-coupled receptors.

Functional interaction structures of the photochromic retinal protein rhodopsin

We studied functional interaction structures of the vertebrate membrane photoreceptor rhodopsin containing retinal as a chromophore. Using time-resolved fluorescence depolarization we analyzed real-time dynamics and conformational changes of the cytoplasmic helix 8 (H8) preceding the long C-terminal tail of rhodopsin. H8 runs parallel to the membrane surface and extends from transmembrane helix 7 whose highly conserved NPxxY(x)F motif connects that region of rhodopsin with the retinal binding pocket. Our measurements indicate that photo-induced retinal isomerization from 11-cis to all-trans provokes conformational changes of H8, including slower motion and reduced flexibility, that are specific for the active metarhodopsin-II photo-intermediate. These conformational changes are absent in the retinal-devoid state opsin and in the phosphorylated metarhodopsin-II state upon receptor deactivation. Furthermore we show that membrane rim effects can influence interfacial reactions at the cytoplasmic rhodopsin surface such as proton transfer reactions between surface and aqueous bulk phase or binding of the signaling protein transducin visualized with single-molecule widefield microscopy. These findings are important for an understanding of the effects of membrane structure on the photo-transduction mechanism.

Molecular mechanism of long-range synergetic color tuning between multiple amino acid residues in conger rhodopsin

The synergetic effects of multiple rhodopsin mutations on color tuning need to be completely elucidated. Systematic genetic studies and spectroscopy have demonstrated an interesting example of synergetic color tuning between two amino acid residues in conger rhodopsin's ancestral pigment (p501): -a double mutation at one nearby and one distant residue led to a significant λ(max) blue shift of 13 nm, whereas neither of the single mutations at these two sites led to meaningful shifts.To analyze the molecular mechanisms of this synergetic color tuning, we performed homology modeling, molecular simulations, and electronic state calculations. For the double mutant, N195A/A292S, in silico mutation analysis demonstrated conspicuous structural changes in the retinal chromophore, whereas that of the single mutant, A292S, was almost unchanged. Using statistical ensembles of QM/MM optimized structures, the excitation energy of retinal chromophore was evaluated for the three visual pigments. As a result, the λ(max) shift of double mutant (DM) from p501 was -8 nm, while that of single mutant (SM) from p501 was +1 nm. Molecular dynamics simulation for DM demonstrated frequent isomerization between 6-s-cis and 6-s-trans conformers. Unexpectedly, however, the two conformers exhibited almost identical excitation energy, whereas principal component analysis (PCA) identified the retinal-counterion cooperative change of BLA (bond length alternation) and retinal-counterion interaction lead to the shift.

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

Structural restraints provided by solid-state NMR measurements of the metarhodopsin II intermediate are combined with molecular dynamics simulations to help visualize structural changes in the light activation of rhodopsin. Since the timescale for the formation of the metarhodopsin II intermediate (>1 ms) is beyond that readily accessible by molecular dynamics, we use NMR distance restraints derived from 13C dipolar recoupling measurements to guide the simulations. The simulations yield a working model for how photoisomerization of the 11-cis retinylidene chromophore bound within the interior of rhodopsin is coupled to transmembrane helix motion and receptor activation. The mechanism of activation that emerges is that multiple switches on the extracellular (or intradiscal) side of rhodopsin trigger structural changes that converge to disrupt the ionic lock between helices H3 and H6 on the intracellular side of the receptor.

A G protein-coupled receptor at work: the rhodopsin model

G protein-coupled receptors (GPCRs) are ubiquitous signal transducers in cell membranes, as well as important drug targets. Interaction with extracellular agonists turns the seven transmembrane helix (7TM) scaffold of a GPCR into a catalyst for GDP and GTP exchange in heterotrimeric Galphabetagamma proteins. Activation of the model GPCR, rhodopsin, is triggered by photoisomerization of its retinal ligand. From the augmentation of biochemical and biophysical studies by recent high-resolution 3D structures, its activation intermediates can now be interpreted as the stepwise engagement of protein domains. Rearrangement of TM5-TM6 opens a crevice at the cytoplasmic side of the receptor into which the C terminus of the Galpha subunit can bind. The Galpha C-terminal helix is used as a transmission rod to the nucleotide binding site. The mechanism relies on dynamic interactions between conserved residues and could therefore be common to other GPCRs.

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state (2)H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific (2)H labels have been introduced into the methyl groups of retinal and solid-state (2)H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent (2)H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the beta-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the beta4 strand of the E2 loop and the side chains of Glu(122) and Trp(265) within the binding pocket. The solid-state (2)H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.

Rhodopsin activation switches in a native membrane environment

The elucidation of structure-function relationships of membrane proteins still poses a considerable challenge due to the sometimes profound influence of the lipid bilayer on the functional properties of the protein. The visual pigment rhodopsin is a prototype of the family of G protein-coupled transmembrane receptors and a considerable part of our knowledge on its activation mechanisms has been derived from studies on detergent-solubilized proteins. This includes in particular the events associated with the conformational transitions of the receptor from the still inactive Meta I to the Meta II photoproduct states, which are involved in signaling. These events involve disruption of an internal salt bridge of the retinal protonated Schiff base, movement of helices and proton uptake from the solvent by the conserved cytoplasmic E(D)RY network around Glu134. As the equilibria associated with these events are considerably altered by the detergent environment, we set out to investigate these equilibria in the native membrane environment and to develop a coherent thermodynamic model of these activating steps using UV-visible and Fourier-transform infrared spectroscopy as complementary techniques. Particular emphasis is put on the role of protonation of Glu134 from the solvent, which is a thermodynamic prerequisite for full receptor activation in membranes, but not in detergent. In view of the conservation of this carboxylate group in family A G protein-coupled receptors, it may also play a similar role in the activation of other family members.

Retinal conformation and dynamics in activation of rhodopsin illuminated by solid-state H NMR spectroscopy

Solid-state NMR spectroscopy gives a powerful avenue for investigating G protein-coupled receptors and other integral membrane proteins in a native-like environment. This article reviews the use of solid-state (2)H NMR to study the retinal cofactor of rhodopsin in the dark state as well as the meta I and meta II photointermediates. Site-specific (2)H NMR labels have been introduced into three regions (methyl groups) of retinal that are crucially important for the photochemical function of rhodopsin. Despite its phenomenal stability (2)H NMR spectroscopy indicates retinal undergoes rapid fluctuations within the protein binding cavity. The spectral lineshapes reveal the methyl groups spin rapidly about their three-fold (C(3)) axes with an order parameter for the off-axial motion of SC(3) approximately 0.9. For the dark state, the (2)H NMR structure of 11-cis-retinal manifests torsional twisting of both the polyene chain and the beta-ionone ring due to steric interactions of the ligand and the protein. Retinal is accommodated within the rhodopsin binding pocket with a negative pretwist about the C11=C12 double bond. Conformational distortion explains its rapid photochemistry and reveals the trajectory of the 11-cis to trans isomerization. In addition, (2)H NMR has been applied to study the retinylidene dynamics in the dark and light-activated states. Upon isomerization there are drastic changes in the mobility of all three methyl groups. The relaxation data support an activation mechanism whereby the beta-ionone ring of retinal stays in nearly the same environment, without a large displacement of the ligand. Interactions of the beta-ionone ring and the retinylidene Schiff base with the protein transmit the force of the retinal isomerization. Solid-state (2)H NMR thus provides information about the flow of energy that triggers changes in hydrogen-bonding networks and helix movements in the activation mechanism of the photoreceptor.

Dynamics of the internal water molecules in squid rhodopsin

Understanding the mechanism of G-protein coupled receptors action is of major interest for drug design. The visual rhodopsin is the prototype structure for the family A of G-protein coupled receptors. Upon photoisomerization of the covalently bound retinal chromophore, visual rhodopsins undergo a large-scale conformational change that prepares the receptor for a productive interaction with the G-protein. The mechanism by which the local perturbation of the retinal cis-trans isomerization is transmitted throughout the protein is not well understood. The crystal structure of the visual rhodopsin from squid solved recently suggests that a chain of water molecules extending from the retinal toward the cytoplasmic side of the protein may play a role in the signal transduction from the all-trans retinal geometry to the activated receptor. As a first step toward understanding the role of water in rhodopsin function, we performed a molecular dynamics simulation of squid rhodopsin embedded in a hydrated bilayer of polyunsaturated lipid molecules. The simulation indicates that the water molecules present in the crystal structure participate in favorable interactions with side chains in the interhelical region and form a persistent hydrogen-bond network in connecting Y315 to W274 via D80.

Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophore's optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Towards an interpretation of 13C chemical shifts in bathorhodopsin, a functional intermediate of a G-protein coupled receptor

Photoisomerization of the membrane-bound light receptor protein rhodopsin leads to an energy-rich photostate called bathorhodopsin, which may be trapped at temperatures of 120 K or lower. We recently studied bathorhodopsin by low-temperature solid-state NMR, using in situ illumination of the sample in a purpose-built NMR probe. In this way we acquired (13)C chemical shifts along the retinylidene chain of the chromophore. Here we compare these results with the chemical shifts of the dark state chromophore in rhodopsin, as well as with the chemical shifts of retinylidene model compounds in solution. An earlier solid-state NMR study of bathorhodopsin found only small changes in the (13)C chemical shifts upon isomerization, suggesting only minor perturbations of the electronic structure in the isomerized retinylidene chain. This is at variance with our recent measurements which show much larger perturbations of the (13)C chemical shifts. Here we present a tentative interpretation of our NMR results involving an increased charge delocalization inside the polyene chain of the bathorhodopsin chromophore. Our results suggest that the bathochromic shift of bathorhodopsin is due to modified electrostatic interactions between the chromophore and the binding pocket, whereas both electrostatic interactions and torsional strain are involved in the energy storage mechanism of bathorhodopsin.

Two protonation switches control rhodopsin activation in membranes

Activation of the G protein-coupled receptor (GPCR) rhodopsin is initiated by light-induced isomerization of the retinal ligand, which triggers 2 protonation switches in the conformational transition to the active receptor state Meta II. The first switch involves disruption of an interhelical salt bridge by internal proton transfer from the retinal protonated Schiff base (PSB) to its counterion, Glu-113, in the transmembrane domain. The second switch consists of uptake of a proton from the solvent by Glu-134 of the conserved E(D)RY motif at the cytoplasmic terminus of helix 3, leading to pH-dependent receptor activation. By using a combination of UV-visible and FTIR spectroscopy, we study the activation mechanism of rhodopsin in different membrane environments and show that these 2 protonation switches become partially uncoupled at physiological temperature. This partial uncoupling leads to approximately 50% population of an entropy-stabilized Meta II state in which the interhelical PSB salt bridge is broken and activating helix movements have taken place but in which Glu-134 remains unprotonated. This partial activation is converted to full activation only by coupling to the pH-dependent protonation of Glu-134 from the solvent, which stabilizes the active receptor conformation by lowering its enthalpy. In a membrane environment, protonation of Glu-134 is therefore a thermodynamic rather than a structural prerequisite for activating helix movements. In light of the conservation of the E(D)RY motif in rhodopsin-like GPCRs, protonation of this carboxylate also may serve a similar function in signal transduction of other members of this receptor family.

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Crystal structure of the ligand-free G-protein-coupled receptor opsin

In the G-protein-coupled receptor (GPCR) rhodopsin, the inactivating ligand 11-cis-retinal is bound in the seven-transmembrane helix (TM) bundle and is cis/trans isomerized by light to form active metarhodopsin II. With metarhodopsin II decay, all-trans-retinal is released, and opsin is reloaded with new 11-cis-retinal. Here we present the crystal structure of ligand-free native opsin from bovine retinal rod cells at 2.9 ångström (A) resolution. Compared to rhodopsin, opsin shows prominent structural changes in the conserved E(D)RY and NPxxY(x)(5,6)F regions and in TM5-TM7. At the cytoplasmic side, TM6 is tilted outwards by 6-7 A, whereas the helix structure of TM5 is more elongated and close to TM6. These structural changes, some of which were attributed to an active GPCR state, reorganize the empty retinal-binding pocket to disclose two openings that may serve the entry and exit of retinal. The opsin structure sheds new light on ligand binding to GPCRs and on GPCR activation.

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 micros, which predict that substantial changes in internal hydration play a functional role in rhodopsin activation. We confirm with (1)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.

Structural impact of the E113Q counterion mutation on the activation and deactivation pathways of the G protein-coupled receptor rhodopsin

Disruption of an interhelical salt bridge between the retinal protonated Schiff base linked to H7 and Glu113 on H3 is one of the decisive steps during activation of rhodopsin. Using previously established stabilization strategies, we engineered a stabilized E113Q counterion mutant that converted rhodopsin to a UV-absorbing photoreceptor with deprotonated Schiff base and allowed reconstitution into native-like lipid membranes. Fourier-transform infrared difference spectroscopy reveals a deprotonated Schiff base in the photoproducts of the mutant up to the active state Meta II, the absence of the classical pH-dependent Meta I/Meta II conformational equilibrium in favor of Meta II, and an anticipation of active state features under conditions that stabilize inactive photoproduct states in wildtype rhodopsin. Glu181 on extracellular loop 2, is found to be unable to maintain a counterion function to the Schiff base on the activation pathway of rhodopsin in the absence of the primary counterion, Glu113. The Schiff base becomes protonated in the transition to Meta III. This protonation is, however, not associated with a deactivation of the receptor, in contrast to wildtype rhodopsin. Glu181 is suggested to be the counterion in the Meta III state of the mutant and appears to be capable of stabilizing a protonated Schiff base in Meta III, but not of constraining the receptor in an inactive conformation.

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.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Roof and floor of the muscarinic binding pocket: variations in the binding modes of orthosteric ligands

Alanine substitution mutagenesis has been used to investigate residues that make up the roof and floor of the muscarinic binding pocket and regulate ligand access. We mutated the amino acids in the second extracellular loop of the M1 muscarinic acetylcholine receptor that are homologous to the cis-retinal contact residues in rhodopsin, the disulfide-bonded Cys178 and Cys98 that anchor the loop to transmembrane helix 3, the adjoining acidic residue Asp99, and the conserved aromatic residues Phe197 and Trp378 in the transmembrane domain. The effects on ligand binding, kinetics, and receptor function suggest that the second extracellular loop does not provide primary contacts for orthosteric ligands, including acetylcholine, but that it does contribute to microdomains that are important for the conformational changes that accompany receptor activation. Kinetic studies suggest that the disulfide bond between Cys98 and Cys178 may contribute to structures that regulate the access of positively charged ligands such as N-methyl scopolamine to the binding pocket. Asp99 may act as a gatekeeper residue to this channel. In contrast, the bulkier lipophilic ligand 3-quinuclidinyl benzilate may require breathing motions of the receptor to access the binding site. Trp378 is a key residue for receptor activation as well as binding, whereas Phe197 represents the floor of the N-methyl scopolamine binding pocket but does not interact with acetylcholine or 3-quinuclidinyl benzilate. Differences between the binding modes of N-methyl scopolamine, 3-quinuclidinyl benzilate, and acetylcholine have been modeled. Although the head groups of these ligands occupy overlapping volumes within the binding site, their side chains may follow significantly different directions.

An opsin shift in rhodopsin: retinal S0-S1 excitation in protein, in solution, and in the gas phase

We considered a series of model systems for treating the photoabsorption of the 11-cis retinal chromophore in the protonated Schiff-base form in vacuum, solutions, and the protein environment. A high computational level, including the quantum mechanical-molecular mechanical (QM/MM) approach for solution and protein was utilized in simulations. The S0-S1 excitation energies in quantum subsystems were evaluated by means of an augmented version of the multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) with the ground-state geometry parameters optimized in the density functional theory PBE0/cc-pVDZ approximation. The computed positions of absorption bands lambdamax, 599(g), 448(s), and 515(p) nm for the gas phase, solution, and protein, respectively, are in excellent agreement with the corresponding experimental data, 610(g), 445(s), and 500(p) nm. Such consistency provides a support for the formulated qualitative conclusions on the role of the chromophore geometry, environmental electrostatic field, and the counterion in different media. An essentially nonplanar geometry conformation of the chromophore group in the region of the C14-C15 bond was obtained for the protein, in particular, owing to the presence of the neighboring charged amino acid residue Glu181. Nonplanarity of the C14-C15 bond region along with the influence of the negatively charged counterions Glu181 and Glu113 are found to be important to reproduce the spectroscopic features of retinal chromophore inside the Rh cavity. Furthermore, the protein field is responsible for the largest bond-order decrease at the C11-C12 double bond upon excitation, which may be the reason for the 11-cis photoisomerization specificity.

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.

Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations

Rhodopsin is currently the only available atomic-resolution template for understanding biological functions of the G protein-coupled receptor (GPCR) family. The structural basis for the phenomenal dark state stability of 11-cis-retinal bound to rhodopsin and its ultrafast photoreaction are active topics of research. In particular, the beta-ionone ring of the retinylidene inverse agonist is crucial for the activation mechanism. We analyzed a total of 23 independent, 100 ns all-atom molecular dynamics simulations of rhodopsin embedded in a lipid bilayer in the microcanonical (N,V,E) ensemble. Analysis of intramolecular fluctuations predicts hydrogen-out-of-plane (HOOP) wagging modes of retinal consistent with those found in Raman vibrational spectroscopy. We show that sampling and ergodicity of the ensemble of simulations are crucial for determining the distribution of conformers of retinal bound to rhodopsin. The polyene chain is rigidly locked into a single, twisted conformation, consistent with the function of retinal as an inverse agonist in the dark state. Most surprisingly, the beta-ionone ring is mobile within its binding pocket; interactions are non-specific and the cavity is sufficiently large to enable structural heterogeneity. We find that retinal occupies two distinct conformations in the dark state, contrary to most previous assumptions. The beta-ionone ring can rotate relative to the polyene chain, thereby populating both positively and negatively twisted 6-s-cis enantiomers. This result, while unexpected, strongly agrees with experimental solid-state (2)H NMR spectra. Correlation analysis identifies the residues most critical to controlling mobility of retinal; we find that Trp265 moves away from the ionone ring prior to any conformational transition. Our findings reinforce how molecular dynamics simulations can challenge conventional assumptions for interpreting experimental data, especially where existing models neglect conformational fluctuations.

Regulation of photoactivation in vertebrate short wavelength visual pigments: protonation of the retinylidene Schiff base and a counterion switch

Xenopus violet cone opsin (VCOP) and its counterion variant (VCOP-D108A) are expressed in mammalian COS1 cells and regenerated with 11-cis-retinal. The phototransduction process in VCOP-D108A is investigated via cryogenic electronic spectroscopy, homology modeling, molecular dynamics, and molecular orbital theory. The VCOP-D108A variant is a UV-like pigment that displays less efficient photoactivation than the mouse short wavelength sensitive visual pigment (MUV) and photobleaching properties that are significantly different. Theoretical calculations trace the difference to the protonation state of the nearby glutamic acid residue E176, which is the homology equivalent of E181 in rhodopsin. We find that E176 is negatively charged in MUV but neutral (protonated) in VCOP-D108A. In the dark state, VCOP-D108A has an unprotonated Schiff base (SB) chromophore (lambdamax = 357 nm). Photolysis of VCOP-D108A at 70 K generates a bathochromic photostationary state (lambdamax = 380 nm). We identify two lumi intermediates, wherein the transitions from batho to the lumi intermediates are temperature- and pH-dependent. The batho intermediate decays to a more red-shifted intermediate called lumi I. The SB becomes protonated during the lumi I to lumi II transition. Decay of lumi II forms meta I, followed by the formation of meta II. We conclude that even in the absence of a primary counterion in VCOP-D108A, the SB becomes protonated during the photoactivation cascade. We examine the relevance of this observation to the counterion switch mechanism of visual pigment activation.

Structural analysis and dynamics of retinal chromophore in dark and meta I states of rhodopsin from 2H NMR of aligned membranes

Rhodopsin is a prototype for G protein-coupled receptors (GPCRs) that are implicated in many biological responses in humans. A site-directed (2)H NMR approach was used for structural analysis of retinal within its binding cavity in the dark and pre-activated meta I states. Retinal was labeled with (2)H at the C5, C9, or C13 methyl groups by total synthesis, and was used to regenerate the opsin apoprotein. Solid-state (2)H NMR spectra were acquired for aligned membranes in the low-temperature lipid gel phase versus the tilt angle to the magnetic field. Data reduction assumed a static uniaxial distribution, and gave the retinylidene methyl bond orientations plus the alignment disorder (mosaic spread). The dark-state (2)H NMR structure of 11-cis-retinal shows torsional twisting of the polyene chain and the beta-ionone ring. The ligand undergoes restricted motion, as evinced by order parameters of approximately 0.9 for the spinning C-C(2)H(3) groups, with off-axial fluctuations of approximately 15 degrees . Retinal is accommodated within the rhodopsin binding pocket with a negative pre-twist about the C11=C12 double bond that explains its rapid photochemistry and the trajectory of 11-cis to trans isomerization. In the cryo-trapped meta I state, the (2)H NMR structure shows a reduction of the polyene strain, while torsional twisting of the beta-ionone ring is maintained. Distortion of the retinal conformation is interpreted through substituent control of receptor activation. Steric hindrance between trans retinal and Trp265 can trigger formation of the subsequent activated meta II state. Our results are pertinent to quantum and molecular mechanics simulations of ligands bound to GPCRs, and illustrate how (2)H NMR can be applied to study their biological mechanisms of action.

Linkage Between the Intramembrane H-bond Network Around Aspartic Acid 83 and the Cytosolic Environment of Helix 8 in Photoactivated Rhodopsin

Understanding the coupling between conformational changes in the intramembrane domain and at the membrane-exposed surface of the bovine photoreceptor rhodopsin, a prototypical G protein-coupled receptor (GPCR), is crucial for the elucidation of molecular mechanisms in GPCR activation. Here, we have combined Fourier transform infrared (FTIR) and fluorescence spectroscopy to address the coupling between conformational changes in the intramembrane region around the retinal and the environment of helix 8, a putative cytosolic surface switch region in class-I GPCRs. Using FTIR/fluorescence cross-correlation we show specifically that surface alterations monitored by emission changes of fluorescein bound to Cys316 in helix 8 of rhodopsin are highly correlated with (i) H-bonding to Asp83 proximal of the retinal Schiff base but not to Glu122 close to the beta-ionone and (ii) with a metarhodopsin II (MII)-specific 1643 cm(-1) IR absorption change, indicative of a partial loss of secondary structure in helix 8 upon MII formation. These correlations are disrupted by limited C-terminal proteolysis but are maintained upon binding of a transducin alpha-subunit (G(talpha))-derived peptide, which stabilizes the MII state. Our results suggest that additional C-terminal cytosolic loop contacts monitored by an amide II absorption at 1557 cm(-1) play a functionally crucial role in keeping helix 8 in the position in which its environment is strongly coupled to the retinal-binding site near the Schiff base. In the intramembrane region, this coupling is mediated by the H-bonding network that connects Asp83 to the NPxxY(x)F motif preceding helix 8.

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Deactivation and proton transfer in light-induced metarhodopsin II/metarhodopsin III conversion: a time-resolved fourier transform infrared spectroscopic study

Vertebrate rhodopsin shares with other retinal proteins the 11-cis-retinal chromophore and the light-induced 11-cis/trans isomerization triggering its activation pathway. However, only in rhodopsin the retinylidene Schiff base bond to the apoprotein is eventually hydrolyzed, making a complex regeneration pathway necessary. Metabolic regeneration cannot be short-cut, and light absorption in the active metarhodopsin (Meta) II intermediate causes anti/syn isomerization around the retinylidene linkage rather than reversed trans/cis isomerization. A new deactivating pathway is thereby triggered, which ends in the Meta III "retinal storage" product. Using time-resolved Fourier transform infrared spectroscopy, we show that the identified steps of receptor activation, including Schiff base deprotonation, protein structural changes, and proton uptake by the apoprotein, are all reversed. However, Schiff base reprotonation is much faster than the activating deprotonation, whereas the protein structural changes are slower. The final proton release occurs with pK approximately 4.5, similar to the pK of a free Glu residue and to the pK at which the isolated opsin apoprotein becomes active. A forced deprotonation, equivalent to the forced protonation in the activating pathway, which occurs against the unfavorable pH of the medium, is not observed. This explains properties of the final Meta III product, which displays much higher residual activity and is less stable than rhodopsin arising from regeneration with 11-cis-retinal. We propose that the anti/syn conversion can only induce a fast reorientation and distance change of the Schiff base but fails to build up the full set of dark ground state constraints, presumably involving the Glu(134)/Arg(135) cluster.

Structural and functional properties of metarhodopsin III: recent spectroscopic studies on deactivation pathways of rhodopsin

The activation of rhodopsin has been the focus of researchers over the past decades, revealing many aspects of the activation pathways of this prototypical G protein-coupled receptor on a molecular level, starting with the light-dependent isomerization of its retinal chromophore from 11-cis to all-trans and leading eventually to the large scale helix movements in the transition to the active receptor state, Meta II. Comparatively little is known, however, on the deactivation pathways of the light receptor, which represent essential steps in maintaining a functional photoreceptor cell. Rhodopsin's active receptor species, Meta II, decays by two fundamentally different pathways, either forming the apoprotein opsin by release of the activating all-trans retinal ligand from its binding pocket, or by a thermal isomerization of this ligand to a less activating species in the transition to metarhodopsin III (Meta III). Both decay products, opsin and Meta III, are largely inactive under physiological conditions, yet they do not restore the complete inactivity of the dark state. Although some properties of Meta III have been described already in the 1960s, its molecular nature and the pathways of its formation have remained rather obscure. In this review, we focus on recent studies from our laboratories, which have provided a major progress in our understanding of the Meta III deactivation pathway and its potential physiological roles. Using Fourier-transform infrared (FTIR) difference spectroscopy in combination with a variety of other spectroscopic and biochemical techniques and quantum chemical calculations, we have developed a general picture of the interplay between the retinal ligand and the receptor protein, which is compared to similar reaction mechanisms in invertebrate photoreceptors and microbial retinal proteins.

Curvature and hydrophobic forces drive oligomerization and modulate activity of rhodopsin in membranes

G protein-coupled receptors (GPCRs) are essential components of cellular signaling pathways. They are the targets of many current pharmaceuticals and are postulated to dimerize or oligomerize in cellular membranes in conjunction with their functional mechanisms. We demonstrate using fluorescence resonance energy transfer how association of rhodopsin occurs by long-range lipid-protein interactions due to geometrical forces, yielding greater receptor crowding. Constitutive association of rhodopsin is promoted by a reduction in membrane thickness (hydrophobic mismatch), but also by an increase in protein/lipid molar ratio, showing the importance of interactions extending well beyond a single annulus of boundary lipids. The fluorescence data correlate with the pK(a) for the MI-to-MII transition of rhodopsin, where deprotonation of the retinylidene Schiff base occurs in conjunction with helical movements leading to activation of the photoreceptor. A more dispersed membrane environment optimizes formation of the MII conformation that results in visual function. A flexible surface model explains both the dispersal and activation of rhodopsin in terms of bilayer curvature deformation (strain) and hydrophobic solvation energy. The bilayer stress is related to the lateral pressure profile in terms of the spontaneous curvature and associated bending rigidity. Transduction of the strain energy (frustration) of the bilayer drives protein oligomerization and conformational changes in a coupled manner. Our findings illuminate the physical principles of membrane protein association due to chemically nonspecific interactions in fluid lipid bilayers. Moreover, they yield a conceptual framework for understanding how the tightly regulated lipid compositions of cellular membranes influence their protein-mediated functions.

Retinal Counterion Switch Mechanism in Vision Evaluated by Molecular Simulations

Photoisomerization of the retinylidene chromophore of rhodopsin is the starting point in the vision cascade. A counterion switch mechanism that stabilizes the retinal protonated Schiff base (PSB) has been proposed to be an essential step in rhodopsin activation. On the basis of vibrational and UV-visible spectroscopy, two counterion switch models have emerged. In the first model, the PSB is stabilized by Glu181 in the meta I state, while in the most recent proposal, it is stabilized by Glu113 as well as Glu181. We assess these models by conducting a pair of microsecond scale, all-atom molecular dynamics simulations of rhodopsin embedded in a 99-lipid bilayer of SDPC, SDPE, and cholesterol (2:2:1 ratio) varying the starting protonation state of Glu181. Theoretical simulations gave different orientations of retinal for the two counterion switch mechanisms, which were used to simulate experimental 2H NMR spectra for the C5, C9, and C13 methyl groups. Comparison of the simulated 2H NMR spectra with experimental data supports the complex-counterion mechanism. Hence, our results indicate that Glu113 and Glu181 stabilize the retinal PSB in the meta I state prior to activation of rhodopsin.

The all-trans-15-syn-retinal chromophore of metarhodopsin III is a partial agonist and not an inverse agonist

Meta III is formed during the decay of rhodopsin's active receptor state at neutral to alkaline pH by thermal isomerization of the retinal Schiff base C15=N bond, converting the ligand from all-trans 15-anti to all-trans 15-syn. The thereby induced change of ligand geometry switches the receptor to an inactive conformation, such that the decay pathway to Meta III contributes to the deactivation of the signaling state at higher pH values. We have examined the conformation of Meta III over a wider pH range and found that Meta III exists in a pH-dependent conformational equilibrium between this inactive conformation at neutral to alkaline pH and an active conformation similar to that of Meta II, which, however, is assumed at very acidic pH only. The apparent pKa of this transition is around 5.1 and thus several units lower than that of the Meta I/Meta II photoproduct equilibrium with its all-trans 15-anti ligand, but still about 1 unit higher than that of the opsin conformational equilibrium in the absence of ligand. The all-trans-15-syn-retinal chromophore is therefore not an inverse agonist like 11-cis- or 9-cis-retinal, which lock the receptor in an inactive conformation, but a classical partial agonist, which is capable of activating the receptor, yet with an efficiency considerably lower than the full agonist all-trans 15-anti. As the Meta III chromophore differs structurally from this full agonist only in the isomeric state of the C15=N bond, this ligand represents an excellent model system to study principal mechanisms of partial agonism which are helpful to understand the partial agonist behavior of other ligands.

Development and application of hybrid structure based method for efficient screening of ligands binding to G-protein coupled receptors

G-protein coupled receptors (GPCRs) comprise a large superfamily of proteins that are targets for nearly 50% of drugs in clinical use today. In the past, the use of structure-based drug design strategies to develop better drug candidates has been severely hampered due to the absence of the receptor's three-dimensional structure. However, with recent advances in molecular modeling techniques and better computing power, atomic level details of these receptors can be derived from computationally derived molecular models. Using information from these models coupled with experimental evidence, it has become feasible to build receptor pharmacophores. In this study, we demonstrate the use of the Hybrid Structure Based (HSB) method that can be used effectively to screen and identify prospective ligands that bind to GPCRs. Essentially; this multi-step method combines ligand-based methods for building enriched libraries of small molecules and structure-based methods for screening molecules against the GPCR target. The HSB method was validated to identify retinal and its analogues from a random dataset of approximately 300,000 molecules. The results from this study showed that the 9 top-ranking molecules are indeed analogues of retinal. The method was also tested to identify analogues of dopamine binding to the dopamine D2 receptor. Six of the ten top-ranking molecules are known analogues of dopamine including a prodrug, while the other thirty-four molecules are currently being tested for their activity against all dopamine receptors. The results from both these test cases have proved that the HSB method provides a realistic solution to bridge the gap between the ever-increasing demand for new drugs to treat psychiatric disorders and the lack of efficient screening methods for GPCRs.

Modulating Rhodopsin Receptor Activation by Altering the pKa of the Retinal Schiff Base

The visual pigment rhodopsin is a seven-transmembrane (7-TM) G protein-coupled receptor (GPCR). Activation of rhodopsin involves two pH-dependent steps: proton uptake at a conserved cytoplasmic motif between TM helices 3 and 6, and disruption of a salt bridge between a protonated Schiff base (PSB) and its carboxylate counterion in the transmembrane core of the receptor. Formation of an artificial pigment with a retinal chromophore fluorinated at C14 decreases the intrinsic pKa of the PSB and thereby destabilizes this salt bridge. Using Fourier transform infrared difference and UV-visible spectroscopy, we characterized the pH-dependent equilibrium between the active photoproduct Meta II and its inactive precursor, Meta I, in the 14-fluoro (14-F) analogue pigment. The 14-F chromophore decreases the enthalpy change of the Meta I-to-Meta II transition and shifts the Meta I/Meta II equilibrium toward Meta II. Combining C14 fluorination with deletion of the retinal beta-ionone ring to form a 14-F acyclic artificial pigment uncouples disruption of the Schiff base salt bridge from transition to Meta II and in particular from the cytoplasmic proton uptake reaction, as confirmed by combining the 14-F acyclic chromophore with the E134Q mutant. The 14-F acyclic analogue formed a stable Meta I state with a deprotonated Schiff base and an at least partially protonated protein counterion. The combination of retinal modification and site-directed mutagenesis reveals that disruption of the protonated Schiff base salt bridge is the most important step thermodynamically in the transition from Meta I to Meta II. This finding is particularly important since deprotonation of the retinal PSB is known to precede the transition to the active state in rhodopsin activation and is consistent with models of agonist-dependent activation of other GPCRs.