Protonation states of membrane-embedded carboxylic acid groups in rhodopsin and metarhodopsin II: a Fourier-transform infrared spectroscopy study of site-directed mutants

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.

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.

Quantum mechanical/molecular mechanical studies on spectral tuning mechanisms of visual pigments and other photoactive proteins

The protein environments surrounding the retinal tune electronic absorption maximum from 350 to 630 nm. Hybrid quantum mechanical/molecular mechanical (QM/MM) methods can be used in calculating excitation energies of retinal in its native protein environments and in studying the molecular basis of spectral tuning. We hereby review recent QM/MM results on the phototransduction of bovine rhodopsin, bacteriorhodopsin, sensory rhodopsin II, nonretinal photoactive yellow protein and their mutants.

An ion-responsive motif in the second transmembrane segment of rhodopsin-like receptors

A L(M)xxxD(N, E) motif (x=a non-ionic amino acid residue, most frequently A, S, L or F; small capitals indicating a minor representation) is found in the second transmembrane (tm2) segment of most G-protein coupling metazoan receptors of the rhodopsin family (Rh-GPCRs). Changes in signal transduction, agonist binding and receptor cycling are known for numerous receptors bearing evolved or experimentally introduced mutations in this tm2 motif, especially of its aspartate residue. The [Na(+)] sensitivity of the receptor-agonist interaction relates to this aspartate in a number of Rh-GPCRs. Native non-conservative mutations in the tm2 motif only rarely coincide with significant changes in two other ubiquitous features of the rhodopsin family, the seventh transmembrane N(D)PxxY(F) motif and the D(E)RY(W,F) or analogous sequence at the border of the third transmembrane helix and the second intracellular loop. Native tm2 mutations with Rh-GPCRs frequently result in constitutive signaling, and with visual opsins also in shifts to short-wavelength sensitivity. Substitution of a strongly basic residue for the tm2 aspartate in Taste-2 receptors could be connected to a lack of sodium sensing by these receptors. These properties could be consistent with ionic interactions, and even of ion transfer, that involve the tm2 motif. A decrease in cation sensing by this motif is usually connected to an enhanced constitutive interaction of the mutated receptors with cognate G- proteins, and also relates to both the constitutive and the overall activity of the short-wavelength opsins.

Generation of a homology model of the human histamine H3 receptor for ligand docking and pharmacophore-based screening

The human histamine H(3) receptor (hH(3)R) is a G-protein coupled receptor (GPCR), which modulates the release of various neurotransmitters in the central and peripheral nervous system and therefore is a potential target in the therapy of numerous diseases. Although ligands addressing this receptor are already known, the discovery of alternative lead structures represents an important goal in drug design. The goal of this work was to study the hH(3)R and its antagonists by means of molecular modelling tools. For this purpose, a strategy was pursued in which a homology model of the hH(3)R based on the crystal structure of bovine rhodopsin was generated and refined by molecular dynamics simulations in a dipalmitoylphosphatidylcholine (DPPC)/water membrane mimic before the resulting binding pocket was used for high-throughput docking using the program GOLD. Alternatively, a pharmacophore-based procedure was carried out where the alleged bioactive conformations of three different potent hH(3)R antagonists were used as templates for the generation of pharmacophore models. A pharmacophore-based screening was then carried out using the program Catalyst. Based upon a database of 418 validated hH(3)R antagonists both strategies could be validated in respect of their performance. Seven hits obtained during this screening procedure were commercially purchased, and experimentally tested in a [(3)H]N(alpha)-methylhistamine binding assay. The compounds tested showed affinities at hH(3)R with K ( i ) values ranging from 0.079 to 6.3 muM.

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.

Molecular Dynamics Simulations of Rhodopsin in Different One-Component Lipid Bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.

Molecular dynamics of rhodopsin and free opsin: computer simulation

Computer simulation was used to perform a comparative study of the molecular dynamics of rhodopsin containing the chromophore group (11-cis-retinal) and free opsin. The molecular dynamics were followed over a time interval of 3000 psec; a total of 3 x 10(6) discrete conformational states of rhodopsin and opsin. The presence of the chromophore group in the chromophore center of opsin was shown to have significant effects on the immediate protein environment of the chromophore and the conformational state of the cytoplasmic domain, but to have virtually no effect on the conformational state of the intradisk domain. The simulation results are used to discuss the possible intramolecular mechanism by which rhodopsin is maintained as a G-protein-coupled receptor in the inactive state, i.e., the function of the chromophore as an effective antagonist ligand.

Coupling of Protonation Switches During Rhodopsin Activation†

Recent studies of the activation mechanism of rhodopsin involving Fourier-transform infrared spectroscopy and a combination of chromophore modifications and site-directed mutagenesis reveal an allosteric coupling between two protonation switches. In particular, the ring and the 9-methyl group of the all-trans retinal chromophore serve to couple two proton-dependent activation steps: proton uptake by a cytoplasmic network between transmembrane (TM) helices 3 and 6 around the conserved ERY (Glu-Arg-Tyr) motif and disruption of a salt bridge between the retinal protonated Schiff base (PSB) and a protein counterion in the TM core of the receptor. Retinal analogs lacking the ring or 9-methyl group are only partial agonists--the conformational equilibrium between inactive Meta I and active Meta II photoproduct states is shifted to Meta I. An artificial pigment was engineered, in which the ring of retinal was removed and the PSB salt bridge was weakened by fluorination of C14 of the retinal polyene. These modifications abolished allosteric coupling of the proton switches and resulted in a stabilized Meta I state with a deprotonated Schiff base (Meta I(SB)). This state had a partial Meta II-like conformation due to disruption of the PSB salt bridge, but still lacked the cytoplasmic proton uptake reaction characteristic of the final transition to Meta II. As activation of native rhodopsin is known to involve deprotonation of the retinal Schiff base prior to formation of Meta II, this Meta I(SB) state may serve as a model for the structural characterization of a key transient species in the activation pathway of a prototypical G protein-coupled receptor.

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.

Critical role of electrostatic interactions of amino acids at the cytoplasmic region of helices 3 and 6 in rhodopsin conformational properties and activation

The cytoplasmic sides of transmembrane helices 3 and 6 of G-protein-coupled receptors are connected by a network of ionic interactions that play an important role in maintaining its inactive conformation. To investigate the role of such a network in rhodopsin structure and function, we have constructed single mutants at position 134 in helix 3 and at positions 247 and 251 in helix 6, as well as combinations of these to obtain double mutants involving the two helices. These mutants have been expressed in COS-1 cells, immunopurified using the rho-1D4 antibody, and studied by UV-visible spectrophotometry. Most of the single mutations did not affect chromophore formation, but double mutants, especially those involving the T251K mutant, resulted in low yield of protein and impaired 11-cis-retinal binding. Single mutants E134Q, E247Q, and E247A showed the ability to activate transducin in the dark, and E134Q and E247A enhanced activation upon illumination, with regard to wild-type rhodopsin. Mutations E247A and T251A (in E134Q/E247A and E134Q/T251A double mutants) resulted in enhanced activation compared with the single E134Q mutant in the dark. A role for Thr(251) in this network is proposed for the first time in rhodopsin. As a result of these mutations, alterations in the hydrogen bond interactions between the amino acid side chains at the cytoplasmic region of transmembrane helices 3 and 6 have been observed using molecular dynamics simulations. Our combined experimental and modeling results provide new insights into the details of the structural determinants of the conformational change ensuing photoactivation of rhodopsin.

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.

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.

Methods to probe protein transitions with ATR infrared spectroscopy

We describe techniques that can be used in conjunction with modern attenuated total reflection (ATR) infrared micro-prisms to allow proteins to be manipulated cyclically between different states whilst simultaneously monitoring both mid-IR and UV/visible/near IR changes. These methods provide increased flexibility of the types of changes that can be induced in proteins in comparison to transmission methods. Quantitative measurements can be made of vibrational changes associated with conversion between stable catalytic reaction intermediates, ligand binding and oxidation-reduction. Both hydrophobic and soluble proteins can be analysed and the ability to induce transitions repetitively allows IR difference spectra to be acquired at a signal/noise sufficient to resolve changes due to specific cofactors or amino acids. Such spectra can often be interpreted at the atomic level by standard IR methods of comparisons with model compounds, by isotope and mutation effects and, increasingly, by ab initio simulations. Combination of such analyses with atomic 3D structural models derived from X-ray and NMR studies can lead to a deeper understanding of molecular mechanisms of enzymatic reactions.

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

A methyl group at C7 of 11-cis-retinal allows chromophore formation but affects rhodopsin activation

The newly synthesized 11-cis-7-methylretinal can form an artificial visual pigment with kinetic and spectroscopic properties similar to the native pigment in the dark-state. However, its photobleaching behavior is altered, showing a Meta I-like photoproduct. This behavior reflects a steric constraint imposed by the 7-methyl group that affects the conformational change in the binding pocket as a result of retinal photoisomerization. Transducin activation is reduced, when compared to the native pigment with 11-cis-retinal. Molecular dynamics simulations suggest coupling of the C7 methyl group and the beta-ionone ring with Met207 in transmembrane helix 5 in agreement with recent experimental results.

The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution

We demonstrate that "brute force" quantum-mechanics/molecular-mechanics computations based on ab initio (i.e., first principles) multiconfigurational perturbation theory can reproduce the absorption maxima of a set of modified bovine rhodopsins with an accuracy allowing for the analysis of the factors determining their colors. In particular, we show that the theory accounts for the changes in excitation energy even when the proteins display the same charge distribution. Three color-tuning mechanisms, leading to changes of close magnitude, are demonstrated to operate in these conditions. The first is based on the change of the conformation of the conjugated backbone of the retinal chromophore. The second operates through the control of the distance between the positive charge residing on the chromophore and the carboxylate counterion. Finally, the third mechanism operates through the changes in orientation of the chromophore relative to the protein. These results offer perspectives for the unbiased computational design of mutants or chemically modified proteins with wanted optical properties.

Time-resolved rapid-scan Fourier transform infrared difference spectroscopy on a noncyclic photosystem: rhodopsin photointermediates from Lumi to Meta II

The visual pigment rhodopsin has been extensively studied for the kinetics of its photointermediates by various spectroscopic methods. Unlike such archaeal retinal proteins as bacteriorhodopsin, visual rhodopsin does not thermally recover its dark state after photoexcitation, which precludes repeated excitation of a single sample and thereby complicates time-resolved experiments. Kinetic data on the late rhodopsin photointermediates have so far been available mainly from time-resolved ultraviolet (UV)-visible spectroscopy, but not from Fourier transform infrared (FTIR) spectroscopy. The latter has the advantage of being informative of structural changes of both chromophore and protein, but does not allow the highly reproducible, automated sample exchange procedures available to UV-visible spectroscopy. Using rapid-scan FTIR difference spectroscopy, we obtained time-resolved data sets that were analyzed by a maximum entropy inverse Laplace-transform. Covering the time range from 8 ms to 15 s at temperatures of 0 and -7 degrees C, the transitions from the Lumi to the Meta I and from the Meta I to the Meta II photoproduct states could be resolved. In the transition from Meta I to Meta II, our data reveal a partial deprotonation of the retinal Schiff base preceding the conformational change of the receptor protein to Meta II. The technique and the results are discussed in regard to its advantages as well as its limitations.

Proton movement and photointermediate kinetics in rhodopsin mutants

The role of ionizable amino acid side chains in the bovine rhodopsin activation mechanism was studied in mutants E134Q, E134R/R135E, H211F, and E122Q. All mutants exhibited bathorhodopsin stability on the 30 ns to 1 micros time scale similar to that of the wild type. Lumirhodopsin decay was also similar to that of the wild type except for the H211F mutant where early decay (20 micros) to a second form of lumirhodopsin was seen, followed by formation of an extremely long-lived Meta I(480) product (34 ms), an intermediate which forms to a much reduced extent, if at all, in dodecyl maltoside suspensions of wild-type rhodopsin. A smaller amount of a similar long-lived Meta I(480) product was seen after photolysis of E122Q, but E134Q and E134R/R135Q displayed kinetics much more similar to those of the wild type under these conditions (i.e., no Meta I(480) product). These results support the idea that specific interaction of His211 and Glu122 plays a significant role in deprotonation of the retinylidene Schiff base and receptor activation. Proton uptake measurements using bromcresol purple showed that E122Q was qualitatively similar to wild-type rhodopsin, with at least one proton being released during lumirhodopsin decay per Meta I(380) intermediate formed, followed by uptake of at least two protons per rhodopsin bleached on a time scale of tens of milliseconds. Different results were obtained for H211F, E134Q, and E134R/R135E, which all released approximately two protons per rhodopsin bleached. These results show that several ionizable groups besides the Schiff base imine are affected by the structural changes involved in rhodopsin activation. At least two proton uptake groups and probably at least one proton release group in addition to the Schiff base are present in rhodopsin.

Computational studies of the primary phototransduction event in visual rhodopsin

This Account addresses recent advances in the elucidation of the detailed molecular rearrangements due to the primary photochemical event in rhodopsin, a prototypical G-protein-coupled receptor (GPCR) responsible for the signal transmission cascade in the vertebrate vision process. The reviewed studies provide fundamental insight on long-standing problems regarding the assembly and function of the individual residues and bound water molecules that form the rhodopsin active site, a center that catalyzes the 11-cis/all-trans isomerization of the retinyl chromophore in the primary step of the phototransduction mechanism. Emphasis is placed on the authors' recent computational studies, based on state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, addressing the structural refinement of the retinyl chromophore binding site in high-resolution X-ray structures of bovine visual rhodopsin, the energy storage mechanism, and the molecular origin of spectroscopic changes due to the primary photochemical event.

Location of Trp265 in metarhodopsin II: implications for the activation mechanism of the visual receptor rhodopsin

Isomerization of the 11-cis retinal chromophore in the visual pigment rhodopsin is coupled to motion of transmembrane helix H6 and receptor activation. We present solid-state magic angle spinning NMR measurements of rhodopsin and the metarhodopsin II intermediate that support the proposal that interaction of Trp265(6.48) with the retinal chromophore is responsible for stabilizing an inactive conformation in the dark, and that motion of the beta-ionone ring allows Trp265(6.48) and transmembrane helix H6 to adopt active conformations in the light. Two-dimensional dipolar-assisted rotational resonance NMR measurements are made between the C19 and C20-methyl groups of the retinal and uniformly 13C-labeled Trp265(6.48). The retinal C20-Trp265(6.48) contact present in the dark-state of rhodopsin is lost in metarhodopsin II, and a new contact is formed with the C19 methyl group. We have previously shown that the retinal translates 4-5 A toward H5 in metarhodopsin II. This motion, in conjunction with the Trp-C19 contact, implies that the Trp265(6.48) side-chain moves significantly upon rhodopsin activation. NMR measurements also show that a packing interaction in rhodopsin between Trp265(6.48) and Gly121(3.36) is lost in metarhodopsin II, consistent with H6 motion away from H3. However, a close contact between Gly120(3.35) on H3 and Met86(2.53) on H2 is observed in both rhodopsin and metarhodopsin II, suggesting that H3 does not change orientation significantly upon receptor activation.

The role of Glu181 in the photoactivation of rhodopsin

The visual pigment rhodopsin is a prototypical seven transmembrane helical G protein-coupled receptor. Photoisomerization of its protonated Schiff base (PSB) retinylidene chromophore initiates a progression of metastable intermediates. We studied the structural dynamics of receptor activation by FTIR spectroscopy of recombinant pigments. Formation of the active state, Meta II, is characterized by neutralization of the PSB and its counterion Glu113. We focused on testing the hypothesis of a PSB counterion switch from Glu113 to Glu181 during the transition of rhodopsin to the still inactive Meta I photointermediate. Our results, especially from studies of the E181Q mutant, support the view that both Glu113 and Glu181 are deprotonated, forming a complex counterion to the PSB in rhodopsin, and that the function of the primary counterion shifts from Glu113 to Glu181 during the transition to Meta I. The Meta I conformation in the E181Q mutant is less constrained compared with that of wild-type Meta I. In particular, the hydrogen bonded network linking transmembrane helices 1, 2, and 7, adopts a conformation that is already Meta II-like, while other parts of the receptor appear to be in a Meta I-like conformation similar to wild-type. We conclude that Glu181 is responsible, in part, for controlling the extraordinary high pK(a) of the chromophore PSB in the dark state, which very likely decreases upon transition to Meta I in a stepwise weakening of the interaction between PSB and its complex counterion during the course of receptor activation. A model for the specific role in coupling chromophore isomerization to protein conformational changes concomitant with receptor activation is presented.

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.

Solvent and Protein Effects on the Structure and Dynamics of the Rhodopsin Chromophore

The structure and dynamics of the retinal chromophore of rhodopsin are investigated systematically in different environments (vacuum, methanol solution, and protein binding pocket) and with different computational approaches (classical, quantum, and hybrid quantum mechanics/molecular mechanics (QM/MM) descriptions). Finite temperature effects are taken into account by molecular dynamics simulations. The different components that determine the structure and dynamics of the chromophore in the protein are dissected, both in the dark state and in the early photointermediates. In vacuum and in solution the chromophore displays a very high flexibility, which is significantly reduced by the protein environment. In the 11-cis chromophore, the bond-length alternation, which is correlated with the dipole moment, is found to be similar in solution and in the protein, while it differs greatly with respect to minimum-energy vacuum structures. In the model of the earliest protein photointermediate, the highly twisted chromophore shows a very reduced bond-length alternation.

Structure, Spectroscopy, and Spectral Tuning of the Gas-Phase Retinal Chromophore: The β-Ionone “Handle” and Alkyl Group Effect

The low-lying singlet states (i.e. S0, S1, and S2) of the chromophore of rhodopsin, the protonated Schiff base of 11-cis-retinal (PSB11), and of its all-trans photoproduct have been studied in isolated conditions by using ab initio multiconfigurational second-order perturbation theory. The computed spectroscopic features include the vertical excitation, the band origin, and the fluorescence maximum of both isomers. On the basis of the S0-->S1 vertical excitation, the gas-phase absorption maximum of PSB11 is predicted to be 545 nm (2.28 eV). Thus, the predicted absorption maximum appears to be closer to that of the rhodopsin pigment (2.48 eV) and considerably red-shifted with respect to that measured in solution (2.82 eV in methanol). In addition, the absorption maxima associated with the blue, green, and red cone visual pigments are tentatively rationalized in terms of the spectral changes computed for PSB11 structures featuring differently twisted beta-ionone rings. More specifically, a blue-shifted absorption maximum is explained in terms of a large twisting of the beta-ionone ring (with respect to the main conjugated chain) in the visual S-cone (blue) pigment chromophore. In contrast, the chromophore of the visual L-cone (red) pigment is expected to have a nearly coplanar beta-ionone ring yielding a six double bond fully conjugated framework. Finally, the M-cone (green) chromophore is expected to feature a twisting angle between 10 and 60 degrees. The spectroscopic effects of the alkyl substituents on the PSB11 spectroscopic properties have also been investigated. It is found that they have a not negligible stabilizing effect on the S1-S0 energy gap (and, thus, cause a red shift of the absorption maximum) only when the double bond of the beta-ionone ring conjugates significantly with the rest of the conjugated chain.

Changes in interhelical hydrogen bonding upon rhodopsin activation

Hydrogen bonding interactions between transmembrane helices stabilize the visual pigment rhodopsin in an inactive conformation in the dark. The crystal structure of rhodopsin has previously revealed that Glu122 and Trp126 on transmembrane helix H3 form a complex hydrogen bonding network with Tyr206 and His211 on H5, while the indole nitrogen of Trp265 on H6 forms a water-mediated hydrogen bond with Asn302 on H7. Here, we use solid-state magic angle spinning NMR spectroscopy to probe the changes in hydrogen bonding upon rhodopsin activation. The NMR chemical shifts of 15N-labeled tryptophan are consistent with the indole nitrogens of Trp126 and Trp265 becoming more weakly hydrogen bonded between rhodopsin and metarhodopsin II. The NMR chemical shifts of 15N-labeled histidine show that His211 is neutral; the unprotonated imidazole nitrogen is not coordinated to zinc in rhodopsin and becomes more strongly hydrogen bonded in metarhodopsin II. Moreover, measurements of rhodopsin containing 13C-labeled histidine show that a strong hydrogen bond between the side-chain of Glu122 and the backbone carbonyl of His211 is disrupted in metarhodopsin II. The implications of these observations for the activation mechanism of rhodopsin are discussed.

A vibrational spectral maker for probing the hydrogen-bonding status of protonated Asp and Glu residues

Hydrogen bonding is a fundamental element in protein structure and function. Breaking a single hydrogen bond may impair the stability of a protein. We report an infrared vibrational spectral marker for probing the hydrogen-bond number for buried, protonated Asp or Glu residues in proteins. Ab initio computational studies were performed on hydrogen-bonding interactions of a COOH group with a variety of side-chain model compounds of polar and charged amino acids in vacuum using density function theory. For hydrogen-bonding interactions with polar side-chain groups, our results show a strong correlation between the C=O stretching frequency and the hydrogen bond number of a COOH group: approximately 1759-1776 cm(-1) for zero, approximately 1733-1749 cm(-1) for one, and 1703-1710 cm(-1) for two hydrogen bonds. Experimental evidence for this correlation will be discussed. In addition, we show an approximate linear correlation between the C=O stretching frequency and the hydrogen-bond strength. We propose that a two-dimensional infrared spectroscopy, C=O stretching versus O-H stretching, may be employed to identify the specific type of hydrogen-bonding interaction. This vibrational spectral marker for hydrogen-bonding interaction is expected to enhance the power of time-resolved Fourier transform infrared spectroscopy for structural characterization of functionally important intermediates of proteins.

Structure, initial excited-state relaxation, and energy storage of rhodopsin resolved at the multiconfigurational perturbation theory level

We demonstrate that a "brute force" quantum chemical calculation based on an ab initio multiconfigurational second order perturbation theory approach implemented in a quantum mechanics/molecular mechanics strategy can be applied to the investigation of the excited state of the visual pigment rhodopsin (Rh) with a computational error <5 kcal.mol(-1). As a consequence, the simulation of the absorption and fluorescence of Rh and its retinal chromophore in solution allows for a nearly quantitative analysis of the factors determining the properties of the protein environment. More specifically, we demonstrate that the Rh environment is more similar to the "gas phase" than to the solution environment and that the so-called "opsin shift" originates from the inability of the solvent to effectively "shield" the chromophore from its counterion. The same strategy is used to investigate three transient structures involved in the photoisomerization of Rh under the assumption that the protein cavity does not change shape during the reaction. Accordingly, the analysis of the initially relaxed excited-state structure, the conical intersection driving the excited-state decay, and the primary isolable bathorhodopsin intermediate supports a mechanism where the photoisomerization coordinate involves a "motion" reminiscent of the so-called bicycle-pedal reaction coordinate. Most importantly, it is shown that the mechanism of the approximately 30 kcal.mol(-1) photon energy storage observed for Rh is not consistent with a model based exclusively on the change of the electrostatic interaction of the chromophore with the protein/counterion environment.

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.

Counterion Controlled Photoisomerization of Retinal Chromophore Models: a Computational Investigation

CASPT2//CASSCF photoisomerization path computations have been used to unveil the effects of an acetate counterion on the photochemistry of two retinal protonated Schiff base (PSB) models: the 2-cis-penta-2,4-dieniminium and the all-trans-epta-2,4,6-trieniminium cations. Different positions/orientations of the counterion have been investigated and related to (i) the spectral tuning and relative stability of the S0, S1, and S2 singlet states; (ii) the selection of the photochemically relevant excited state; (iii) the control of the radiationless decay and photoisomerization rates; and, finally, (iv) the control of the photoisomerization stereospecificity. A rationale for the results is given on the basis of a simple (electrostatic) qualitative model. We show that the model readily explains the computational results providing a qualitative explanation for different aspects of the experimentally observed "environment" dependent PSB photochemistry. Electrostatic effects likely involved in controlling retinal photoisomerization stereoselectivity in the protein are also discussed under the light of these results, and clues for a stereocontrolled electrostatically driven photochemical process are presented. These computations provide a rational basis for the formulation of a mechanistic model for photoisomerization electrostatic catalysis.

Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Rhodopsin is the only member of the pharmacologically important superfamily of G-protein-coupled receptors with a known structure at atomic resolution. A molecular dynamics model of rhodopsin in a POPC phospholipid bilayer was simulated for 15 ns, revealing a conformation significantly different from the recent crystal structures. The structure of the bilayer compared with a protein-free POPC control indicated hydrophobic matching with the nonpolar interface of the receptor, in agreement with deuterium NMR experiments. A new generalized molecular surface method, based on a three-dimensional Voronoi cell construction for atoms with different radii, was developed to quantify cross-sectional area profiles for the protein, lipid acyl chains and headgroups, and water. Thus, it was possible to investigate the bilayer deformation due to curvature of the individual lipid monolayers. Moreover, the generalized molecular surface derived hydrophobic interface allowed benchmarking of the hydropathy sequence analysis, an important structural genomics tool. Five water molecules diffused into internal hydration sites during the simulation, yielding a total of 12 internal waters. The cytoplasmic loops and the C-terminal tail, containing the G-protein recognition and protein sorting sequences, exhibited a high mobility, in marked contrast to the extracellular and transmembrane domains. The proposed functional coupling of the highly conserved ERY motif to the lipid-water interface via the cytoplasmic loops provides insight into lipid effects on G-protein-coupled receptor activation in terms of a flexible surface model, involving the spontaneous monolayer curvature.

The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

A new high-resolution structure is reported for bovine rhodopsin, the visual pigment in rod photoreceptor cells. Substantial improvement of the resolution limit to 2.2 A has been achieved by new crystallization conditions, which also reduce significantly the probability of merohedral twinning in the crystals. The new structure completely resolves the polypeptide chain and provides further details of the chromophore binding site including the configuration about the C6-C7 single bond of the 11-cis-retinal Schiff base. Based on both an earlier structure and the new improved model of the protein, a theoretical study of the chromophore geometry has been carried out using combined quantum mechanics/force field molecular dynamics. The consistency between the experimental and calculated chromophore structures is found to be significantly improved for the 2.2 A model, including the angle of the negatively twisted 6-s-cis-bond. Importantly, the new crystal structure refinement reveals significant negative pre-twist of the C11-C12 double bond and this is also supported by the theoretical calculation although the latter converges to a smaller value. Bond alternation along the unsaturated chain is significant, but weaker in the calculated structure than the one obtained from the X-ray data. Other differences between the experimental and theoretical structures in the chromophore binding site are discussed with respect to the unique spectral properties and excited state reactivity of the chromophore.

QM/MM study of energy storage and molecular rearrangements due to the primary event in vision

The energy storage and the molecular rearrangements due to the primary photochemical event in rhodopsin are investigated by using quantum mechanics/molecular mechanics hybrid methods in conjunction with high-resolution structural data of bovine visual rhodopsin. The analysis of the reactant and product molecular structures reveals the energy storage mechanism as determined by the detailed molecular rearrangements of the retinyl chromophore, including rotation of the (C11-C12) dihedral angle from -11 degrees in the 11-cis isomer to -161 degrees in the all-trans product, where the preferential sense of rotation is determined by the steric interactions between Ala-117 and the polyene chain at the C13 position, torsion of the polyene chain due to steric constraints in the binding pocket, and stretching of the salt bridge between the protonated Schiff base and the Glu-113 counterion by reorientation of the polarized bonds that localize the net positive charge at the Schiff-base linkage. The energy storage, computed at the ONIOM electronic-embedding approach (B3LYP/6-31G*:AMBER) level of theory and the S0-->S1 electronic-excitation energies for the dark and product states, obtained at the ONIOM electronic-embedding approach (TD-B3LYP/6-31G*//B3LYP/6-31G*:AMBER) level of theory, are in very good agreement with experimental data. These results are particularly relevant to the development of a first-principles understanding of the structure-function relations in prototypical G-protein-coupled receptors.

Rhodopsin photoproducts in 2D crystals

The published electron microscope and X-ray structures of rhodopsin have made available a detailed picture of the inactive dark state of rhodopsin. Yet, the photointermediates of rhodopsin that ultimately lead to the activated receptor species still await a similar analysis. Such an analysis first requires the generation and characterization of the photoproducts that can be obtained in crystals of rhodopsin. We therefore studied with Fourier-transform infrared (FTIR) difference spectroscopy the photoproducts in 2D crystals of bovine rhodopsin in a p22(1)2(1) crystal form. The spectra obtained by cryotrapping revealed that in this crystal form the still inactive early intermediates batho, lumi, and meta I are similar to those obtained from rhodopsin in native disk membranes, although the transition from lumi to meta I is shifted to a higher temperature. However, at room temperature, the formation of the active state, meta II, is blocked in the crystalline environment. Instead, an intermediate state is formed that bears some features of meta II but lacks the specific conformational changes required for activity. Despite being unable to activate its cognate G protein, transducin, to a significant extent, this intermediate state is capable of interacting with functional transducin-derived peptides to a limited extent. Therefore, while unable to support formation of rhodopsin's active state meta II, 2D p22(1)2(1) crystals proved to be very suitable for determining 3D structures of its still inactive precursors, batho, lumi, and meta I. In future studies, FTIR spectroscopy may serve as a sensitive assay to screen crystals grown under altered conditions for potential formation of the active state, meta II.

Structure of the intersection space associated with Z/E photoisomerization of retinal in rhodopsin proteins

In this paper we employ a CASSCF/AMBER quantum-mechanics/molecular-mechanics tool to map the intersection space (IS) of a protein. In particular, we provide evidence that the S1 excited-state potential-energy surface of the visual photoreceptor rhodopsin is spanned by an IS segment located right at the bottom of the surface. Analysis of the molecular structures of the protein chromophore (a protonated Schiff base of retinal) along IS reveals a type of geometrical deformation not observed in vacuo. Such a structure suggests that conical intersections mediating different photochemical reactions reside along the same intersection space. This conjecture is investigated by mapping the intersection space of the rhodopsin chromophore model 2-Z-hepta-2,4,6-trieniminium cation and of the conjugated hydrocarbon 3-Z-deca-1,3,5,6,7-pentaene.

Fourier transform IR spectroscopy study for new insights into molecular properties and activation mechanisms of visual pigment rhodopsin

Fourier transform IR (FTIR) spectroscopy has been successfully applied in recent years to examine the functional and structural properties of the membrane protein rhodopsin, a prototype G protein coupled receptor. Unlike UV-visible spectroscopy, FTIR spectroscopy is structurally sensitive. It may give us both global information about the conformation of the protein and very detailed information about the retinal chromophore and all other functional groups, even when these are not directly related to the chromophore. Furthermore, it can be successfully applied to the photointermediates of rhodopsin, including the active receptor species, metarhodopsin II, and its decay products, which is not expected presently or even in the near future from crystallographic approaches. In this review we show how FTIR spectroscopy has significantly contributed to the understanding of very different aspects of rhodopsin, comprising both structural properties and the mechanisms leading to receptor activation and deactivation.

Probing the rhodopsin cavity with reduced retinal models at the CASPT2//CASSCF/AMBER level of theory

We show that the ab initio CASPT2//CASSCF strategy previously used to investigate the ground and excited states of the chromophore of the vision receptor rhodopsin (Rh) in vacuo can be successfully implemented in a QM/MM scheme allowing for CASPT2//CASSCF/AMBER geometry optimization and excited state property evaluation in proteins. Two receptor models (Rh-1 and Rh-2) incorporating different reduced chromophores are investigated. It is shown that Rh-2 features a chromophore equilibrium structure with the correct helicity and a lambdamax that is only 52 nm blue-shifted from the observed value. This result should open the way to a qualitatively correct ab initio QM/MM modeling of the early excited state transient species involved in the vision process.

Molecular dynamics investigation of primary photoinduced events in the activation of rhodopsin

Retinal cis-trans isomerization and early relaxation steps have been studied in a 10-ns molecular dynamics simulation of a fully hydrated model of membrane-embedded rhodopsin. The isomerization, induced by transiently switching the potential energy function governing the C(11)==C(12) dihedral angle of retinal, completes within 150 fs and yields a strongly distorted retinal. The most significant conformational changes in the binding pocket are straightening of retinal's polyene chain and separation of its beta-ionone ring from Trp-265. In the following 500 ps, transition of 6s-cis to 6s-trans retinal and dramatic changes in the hydrogen bonding network of the binding pocket involving the counterion for the protonated Schiff base, Glu-113, occur. Furthermore, the energy initially stored internally in the distorted retinal is transformed into nonbonding interactions of retinal with its environment. During the following 10 ns, increased mobilities of some parts of the protein, such as the kinked regions of the helices, mainly helix VI, and the intracellular loop I2, were observed, as well as transient structural changes involving the conserved salt bridge between Glu-134 and Arg-135. These features prepare the protein for major structural transformations achieved later in the photocycle. Retinal's motion, in particular, can be compared to an opening turnstile freeing the way for the proposed rotation of helix VI. This was demonstrated by a steered molecular dynamics simulation in which an applied torque enforced the rotation of helix VI.

Rhodopsin with 11-cis-locked chromophore is capable of forming an active state photoproduct

The visual pigment rhodopsin is characterized by an 11-cis retinal chromophore bound to Lys-296 via a protonated Schiff base. Following light absorption the C(11)=C(12) double bond isomerizes to trans configuration and triggers protein conformational alterations. These alterations lead to the formation of an active intermediate (Meta II), which binds and activates the visual G protein, transducin. We have examined by UV-visible and Fourier transform IR spectroscopy the photochemistry of a rhodopsin analogue with an 11-cis-locked chromophore, where cis to trans isomerization around the C(11)=C(12) double bond is prevented by a 6-member ring structure (Rh(6.10)). Despite this lock, the pigment was found capable of forming an active photoproduct with a characteristic protein conformation similar to that of native Meta II. This intermediate is further characterized by a protonated Schiff base and protonated Glu-113, as well as by its ability to bind a transducin-derived peptide previously shown to interact efficiently with native Meta II. The yield of this active photointermediate is pH-dependent and decreases with increasing pH. This study shows that with the C(11)=C(12) double bond being locked, isomerization around the C(9)=C(10) or the C(13)=C(14) double bonds may well lead to an activation of the receptor. Additionally, prolonged illumination at pH 7.5 produces a new photoproduct absorbing at 385 nm, which, however, does not exhibit the characteristic active protein conformation.

Rhodopsin: insights from recent structural studies

The recent report of the crystal structure of rhodopsin provides insights concerning structure-activity relationships in visual pigments and related G protein-coupled receptors (GPCRs). The seven transmembrane helices of rhodopsin are interrupted or kinked at multiple sites. An extensive network of interhelical interactions stabilizes the ground state of the receptor. The ligand-binding pocket of rhodopsin is remarkably compact, and several chromophore-protein interactions were not predicted from mutagenesis or spectroscopic studies. The helix movement model of receptor activation, which likely applies to all 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 includes a helical domain extending from the seventh transmembrane segment parallel to the bilayer surface. The cytoplasmic surface appears to be approximately large enough to bind to the transducin heterotrimer in a one-to-one complex. The structural basis for several unique biophysical properties of rhodopsin, including its extremely low dark noise level and high quantum efficiency, can now be addressed using a combination of structural biology and various spectroscopic methods. Future high-resolution structural studies of rhodopsin and other GPCRs will form the basis to elucidate the detailed molecular mechanism of GPCR-mediated signal transduction.

Receptor-dependent G-protein activation in lipidic cubic phase

The primary step in cellular signaling by G-protein-coupled receptors (GPCRs) is the interaction of the agonist-activated transmembrane receptor with an intracellular G-protein. Understanding the underlying molecular mechanisms requires the structural determination of receptor G-protein complexes that are not yet achieved. The crystal structure of the bovine photoreceptor rhodopsin, a prototypical GPCR, was solved recently and the structures of different states of engineered G-proteins were reported. Posttranslational hydrophobic modifications of G-proteins are in most cases removed for crystallization but play functional roles for interactions among G-protein subunits with receptors, as well as membranes. Bovine rhodopsin is reconstituted into lipidic cubic phases to assess their potential for crystallization of receptor G-protein complexes under conditions that may preserve the structural and functional roles of hydrophobic protein modifications. Three-dimensional bilayers of a bicontinuous lipidic cubic phase are successfully employed for crystallization of membrane and soluble proteins. UV-visible absorption and attenuated total reflection Fourier transform IR difference spectroscopy reveal that light activation of cubic phase reconstituted rhodopsin results in the generation of a metarhodopsin II-like state. Via diffusion along aqueous channels, transducin couples efficiently to this photoproduct as evidenced by the nucleotide-dependent increase of transducin fluorescence. Thus, rhodopsin transducin interactions do not crucially depend on the presence of sn1 and sn2 acyl chains, phospholipid head groups, or membrane planarity. Because lipidic cubic phases preserve the essential functional and structural properties of native rhodopsin and transducin, they appear suitable for the detergent-free crystallization of receptor G-protein complexes carrying a normal pattern of hydrophobic modifications.