Solid-state NMR studies of the mechanism of the opsin shift in the visual pigment rhodopsin

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

pH-dependent absorption spectra of rhodopsin mutant E113Q: On the role of counterions and protein

The absorption spectra of bovine rhodopsin mutant E113Q in solutions were investigated at the molecular level by using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The calculations suggest the mechanism of the absorption variations of E113Q at different pH values. The results indicate that the polarizations of the counterions in the vicinity of Schiff base under protonation and unprotonation states of the mutant E113Q would be a crucial factor to change the energy gap of the retinal to tune the absorption spectra. Glu-181 residue, which is close to the chromophore, cannot serve as the counterion of the protonated Schiff base of E113Q in dark state. Moreover, the results of the absorption maximum in mutant E113Q with the various anions (Cl^-, Br^-, I^- and NO₃^-) manifested that the mutant E113Q could have the potential for use as a template of anion biosensors at visible wavelength.

Opsin Effect on the Electronic Structure of the Retinylidene Chromophore in Rhodopsin

Direct examination of experimental NMR parameters combined with electronic structure analysis was used to provide a first-principle interpretation of NMR experiments and give a precise evaluation of how the electronic perturbation of the protein environment affects the electronic properties of the retinylidene chromophere in rhodopsin. To this end, we pursued a theoretical analysis using a combination of tools including quantum mechanics/molecular mechanics (QM/MM) at the Density Functional Theory (DFT) level, in conjunction with gauge independent atomic orbital (GIAO) calculations of (13)C NMR chemical shieldings and (1)J(CC) spin-spin coupling constants obtained with the Coupled Perturbed DFT (CPDFT) method. The opsin effect on the retinylidene chromophere is interpreted as an inductive effect of Glu-113 which readjusts the weighting factors of resonance substructures of the conjugated chain of the chromophere. These changes give a rationalization to the alternating effect of the (13)C chemical shifts magnitudes when comparing the retinylidene chromophere in the presence and absence of the protein environment. Conversely, perturbation of π orbitals has little to no effect over (1)J (13)C-(13)C spin-spin coupling constants, as they are mainly dominated by the Fermi contact term, and hence the counteraion effect is restricted to the vicinity of the perturbation. Thus, the apparent contradiction between experimental findings based on chemical shifts (deep penetration) and one-bond J-couplings (localized effects of the protonated Schiff base at the chain terminus) is in fact a consequence of different properties responding differently to the same external perturbation.

Low-Temperature Trapping of Photointermediates of the Rhodopsin E181Q Mutant

Three active-site components in rhodopsin play a key role in the stability and function of the protein: 1) the counter-ion residues which stabilize the protonated Schiff base, 2) water molecules, and 3) the hydrogen-bonding network. The ionizable residue Glu-181, which is involved in an extended hydrogen-bonding network with Ser-186, Tyr-268, Tyr-192, and key water molecules within the active site of rhodopsin, has been shown to be involved in a complex counter-ion switch mechanism with Glu-113 during the photobleaching sequence of the protein. Herein, we examine the photobleaching sequence of the E181Q rhodopsin mutant by using cryogenic UV-visible spectroscopy to further elucidate the role of Glu-181 during photoactivation of the protein. We find that lower temperatures are required to trap the early photostationary states of the E181Q mutant compared to native rhodopsin. Additionally, a Blue Shifted Intermediate (BSI, λ_max = 498 nm, 100 K) is observed after the formation of E181Q Bathorhodopsin (Batho, λ_max = 556 nm, 10 K) but prior to formation of E181Q Lumirhodopsin (Lumi, λ_max = 506 nm, 220 K). A potential energy diagram of the observed photointermediates suggests the E181Q Batho intermediate has an enthalpy value 7.99 KJ/mol higher than E181Q BSI, whereas in rhodopsin, the BSI is 10.02 KJ/mol higher in enthalpy than Batho. Thus, the Batho to BSI transition is enthalpically driven in E181Q and entropically driven in native rhodopsin. We conclude that the substitution of Glu-181 with Gln-181 results in a significant perturbation of the hydrogen-bonding network within the active site of rhodopsin. In addition, the removal of a key electrostatic interaction between the chromophore and the protein destabilizes the protein in both the dark state and Batho intermediate conformations while having a stabilizing effect on the BSI conformation. The observed destabilization upon this substitution further supports that Glu-181 is negatively charged in the early intermediates of the photobleaching sequence of rhodopsin.

Isotropic magnetic shielding constants of retinal derivatives in aprotic and protic solvents

We investigate the nuclear isotropic shielding constants σ((13)C) and σ((17)O) of isomers of retinoic acid and retinal in gas-phase and in chloroform, acetonitrile, methanol, and water solutions via Monte Carlo simulation and quantum mechanics calculations using the GIAO-B3LYP∕6-311++G(2d,2p) approach. Electronic solute polarization effects due to protic and aprotic solvents are included iteratively and play an important role in the quantitative determination of oxygen shielding constants. Our MP2∕6-31G+(d) results show substantial increases of the dipole moment of both retinal derivatives in solution as compared with the gas-phase results (between 22% and 26% in chloroform and between 55% and 99% in water). For the oxygen atoms the influence of the solute polarization is mild for σ((17)O) of hydroxyl group, even in protic solvents, but it is particularly important for σ((17)O) of carbonyl group. For the latter, there is a sizable increase in the magnitude with increasing solvent polarity. For the carbon atoms, the solvent effects on the σ((13)C) values are in general small, being more appreciable in carbon atoms of the polyene chain than in the carbon atoms of the β-ionone ring and methyl groups. The results also show that isomeric changes on the backbones of the polyene chains have marked influence on the (13)C chemical shifts of carbon atoms near to the structural distortions, in good agreement with the experimental results measured in solution.

Glutamic acid 181 is negatively charged in the bathorhodopsin photointermediate of visual rhodopsin

Assignment of the protonation state of the residue Glu-181 is important to our understanding of the primary event, activation processes and wavelength selection in rhodopsin. Despite extensive study, there is no general agreement on the protonation state of this residue in the literature. Electronic assignment is complicated by the location of Glu-181 near the nodal point in the electrostatic charge shift that accompanies excitation of the chromophore into the low-lying, strongly allowed ππ* state. Thus, the charge on this residue is effectively hidden from electronic spectroscopy. This situation is resolved in bathorhodopsin, because photoisomerization of the chromophore places Glu-181 well within the region of negative charge shift following excitation. We demonstrate that Glu-181 is negatively charged in bathorhodopsin on the basis of the shift in the batho absorption maxima at 10 K [λ(max) band (native) = 544 ± 2 nm, λ(max) band (E181Q) = 556 ± 3 nm] and the decrease in the λ(max) band oscillator strength (0.069 ± 0.004) of E181Q relative to that of the native protein. Because the primary event in rhodopsin does not include a proton translocation or disruption of the hydrogen-bonding network within the binding pocket, we may conclude that the Glu-181 residue in rhodopsin is also charged.

The protonation state of Glu181 in rhodopsin revisited: interpretation of experimental data on the basis of QM/MM calculations

The structure and spectroscopy of rhodopsin have been intensely studied in the past decade both experimentally and theoretically; however, important issues still remain unresolved. Of central interest is the protonation state of Glu181, where controversial and contradictory experimental evidence has appeared. While FTIR measurements indicate this residue to be unprotonated, preresonance Raman and UV-vis spectra have been interpreted in favor of a protonated Glu181. Previous computational approaches were not able to resolve this issue, providing contradicting data as well. Here, we perform hybrid QM/MM calculations using DFT methods for the electronic ground state, MRCI methods for the electronically excited states, and a polarization model for the MM part in order to investigate this issue systematically. We constructed various active-site models for protonated as well as unprotonated Glu181, which were evaluated by computing NMR, IR, Raman, and UV-vis spectroscopic data. The resulting differences in the UV-vis and Raman spectra between protonated and unprotonated models are very subtle, which has two major consequences. First, the common interpretation of prior Raman and UV-vis experiments in favor of a neutral Glu181 appears questionable, as it is based on the assumption that a charge at the Glu181 location would have a sizable impact. Second, also theoretical results should be interpreted with care. Spectroscopic differences between the structural models must be related to modeling uncertainties and intrinsic methodological errors. Despite a detailed comparison of various rhodopsins and mutants and consistently favorite results with charged Glu181 models, we find merely weak evidence from UV-vis and Raman calculations. On the contrary, difference FTIR and NMR chemical shift measurements on Rh mutants are indicative of the protonation state of Glu181. Supported by our results, they provide strong and independent evidence for a charged Glu181.

Correlation between absorption maxima and thermal isomerization rates in bacteriorhodopsin

The reported rates of thermal 13-cis to all-trans isomerization of the protonated Schiff base of retinal (PSBR) in solution and in bacteriorhodopsin (BR) are shown to be correlated with the red shift in the absorption maximum of the chromophore, though the linear fit is different for BR and for a model PSBR in solution. Because the red shift in the absorption has been previously shown to be correlated with pi-electron delocalization in the chromophore, this suggests that the thermal isomerization rate is largely regulated by the amount of double bond character in the chromophore. Because the linear fit of isomerization rates with absorption maxima is different for BR and the model PSBR, specific interactions of the protein with the chromophore must also be a factor in determining thermal isomerization rates in BR. A model of the later steps in the photocycle of BR is presented in which the 13-cis to all-trans thermal isomerization occurs during the O intermediate.

Synthesis and use of stable isotope enriched retinals in the field of vitamin A

The role of vitamin A and its metabolites in the life processes starting with the historical background and its up to date information is discussed in the introduction. Also the role of 11Z-retinal in vision and retinoic acid in the biological processes is elucidated. The essential role of isotopically enriched systems in the progress of vision research, nutrition research etc. is discussed. In part B industrial commercial syntheses of vitamin A by the two leading companies Hoffmann-La Roche (now DSM) and BASF are discussed. The knowledge obtained via these pioneering syntheses has been essential for the further synthetic efforts in vitamin A field by other scientific groups. The rest of the paper is devoted to the synthetic efforts of the Leiden group that gives an access to the preparation of site directed high level isotope enrichment in retinals. First the synthesis of the retinals with deuterium incorporation in the conjugated side chain is reviewed. Then, 13C-labeled retinals are discussed. This is followed by the discussion of a convergent synthetic scheme that allows a rational access to prepare any isotopomer of retinals. The schemes that provide access to prepare any possible isotope enriched chemically modified systems are discussed. Finally, nor-retinals and bridged retinals that give access to a whole (as yet incomplete) library of possible isotopomers are reviewed.

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.

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.

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.

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.

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.

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.

G-protein coupled receptor structure

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

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.

A 3-D structural model of solid self-assembled chlorophyll a/H(2)O from multispin labeling and MAS NMR 2-D dipolar correlation spectroscopy in high magnetic field

Magic angle spinning (MAS) NMR with Lee-Goldburg cross-polarization (LG-CP) is used to promote long-range heteronuclear transfer of magnetization and to constrain a structural model for uniformly labeled chlorophyll a/H(2)O. An effective maximum transfer range d(max) can be determined experimentally from the detection of a gradually decreasing series of intramolecular correlations with the (13)C along the molecular skeleton. To probe intermolecular contacts, d(max) can be set to approximately 4.2 A by choosing an LG-CP contact time of 2 ms. Long-range (1)H-(13)C correlations are used in conjunction with carbon and proton aggregation shifts to establish the stacking of the chlorophyll a (Chl a) molecules. First, high-field (14.1 T) 2-D MAS NMR homonuclear ((13)C-(13)C) dipolar correlation spectra provide a complete assignment of the carbon chemical shifts. Second, proton chemical shifts are obtained from (1)H-(13)C heteronuclear dipolar correlation spectroscopy in high magnetic field. The shift constraints and long-range (1)H-(13)C intermolecular correlations reveal a 2-D stacking homologous to the molecular arrangement in crystalline solid ethyl-chlorophyllide a. A doubling of a small subset of the carbon resonances, in the 7-methyl region of the molecule, provides evidence for two marginally different well-defined molecular environments. Evidence is found for the presence of neutral structural water molecules forming a hydrogen-bonded network to stabilize Chl a sheets. In line with the microcrystalline order observed for the rings, the long T(1)'s, and absence of conformational shifts for the (13)C in the phytyl tails, it is proposed that the Chl a form a rigid 3-D space-filling structure. Probably the only way this can be realized with the sheets is by forming bilayers with interpenetration of elongated tails. Such a 3-D space-filling organization of the aggregated Chl a from MAS NMR would match existing models inferred from electron microscopy and low-resolution X-ray powder diffraction, while a micellar model based on neutron diffraction and antiparallel stacking observed in solution can be discarded.

Retinylidene ligand structure in bovine rhodopsin, metarhodopsin-I, and 10-methylrhodopsin from internuclear distance measurements using 13C-labeling and 1-D rotational resonance MAS NMR

Rhodopsin is the G-protein coupled photoreceptor that initiates the rod phototransduction cascade in the vertebrate retina. Using specific isotope enrichment and magic angle spinning (MAS) NMR, we examine the spatial structure of the C10-C11=C12-C13-C20 motif in the native retinylidene chromophore, its 10-methyl analogue, and the predischarge photoproduct metarhodopsin-I. For the rhodopsin study 11-Z-[10,20-(13)C(2)]- and 11-Z-[11,20-(13)C(2)]-retinal were synthesized and incorporated into bovine opsin while maintaining a natural lipid environment. The ligand is covalently bound to Lys(296) in the photoreceptor. The C10-C20 and C11-C20 distances were measured using a novel 1-D CP/MAS NMR rotational resonance experimental procedure that was specifically developed for the purpose of these measurements [Verdegem, P. J. E., Helmle, M., Lugtenburg, J., and de Groot, H. J. M. (1997) J. Am. Chem. Soc. 119, 169]. We obtain r(10,20) = 0.304 +/- 0.015 nm and r(11,20) = 0.293 +/- 0.015 nm, which confirms that the retinylidene is 11-Z and shows that the C10-C13 unit is conformationally twisted. The corresponding torsional angle is about 44 degrees as indicated by Car-Parrinello modeling studies. To increase the nonplanarity in the chromophore, 11-Z-[10,20-(13)C(2)]-10-methylretinal and 11-Z-[(10-CH(3)), 13-(13)C(2)]-10-methylretinal were prepared and incorporated in opsin. For the resulting analogue pigment r(10,20) = 0.347 +/- 0.015 nm and r((10)(-)(CH)()3())(,)(13) = 0.314 +/- 0.015 nm were obtained, consistent with a more distorted chromophore. The analogue data are in agreement with the induced fit principle for the interaction of opsin with modified retinal chromophores. Finally, we determined the intraligand distances r(10,20) and r(11,20) also for the photoproduct metarhodopsin-I, which has a relaxed all-E structure. The results (r(10,20) >/= 0.435 nm and r(11,20) = 0.283 +/- 0.015 nm) fully agree with such a relaxed all-E structure, which further validates the 1-D rotational resonance technique for measuring intraligand distances and probing ligand structure. As far as we are aware, these results represent the first highly precise distance determinations in a ligand at the active site of a membrane protein. Overall, the MAS NMR data indicate a tight binding pocket, well defined to bind specifically only one enantiomer out of four possibilities and providing a steric complement to the chromophore in an ultrafast ( approximately 200 fs) isomerization process.

BUNDLE: a program for building the transmembrane domains of G-protein-coupled receptors

The only information available at present about the structural features of G-protein-coupled receptors (GPCRs) comes from low resolution electron density maps of rhodopsin obtained from electron microscopy studies on 2D crystals. Despite their low resolution, maps can be used to extract information about transmembrane helix relative positions and their tilt. This information, together with a reliable algorithm to assess the residues involved in each of the membrane spanning regions, can be used to construct a 3D model of the transmembrane domains of rhodopsin at atomic resolution. In the present work, we describe an automated procedure applicable to generate such a model and, in general, to construct a 3D model of any given GPCR with the only assumption that it adopts the same helix arrangement as in rhodopsin. The present approach avoids uncertainties associated with other procedures available for constructing models of GPCRs based on a template, since sequence identity among GPCRs of different families in most of the cases is not significant. The steps involved in the construction of the model are: (i) locate the centers of the helices according to the low-resolution electron density map; (ii) compute the tilt of each helix based on the elliptical shape observed by each helix in the map; (iii) define a local coordinate system for each of the helices; (iv) bring them together in an antiparallel orientation; (v) rotate each helix through the helical axis in such a way that its hydrophobic moment points in the same direction of the bisector formed between three consecutive helices in the bundle; (vi) rotate each helix through an axis perpendicular to the helical one to assign a proper tilt; and (vii) translate each helix to its center deduced from the projection map.

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.

Automated method for modeling seven-helix transmembrane receptors from experimental data

A rule-based automated method is presented for modeling the structures of the seven transmembrane helices of G-protein-coupled receptors. The structures are generated by using a simulated annealing Monte Carlo procedure that positions and orients rigid helices to satisfy structural restraints. The restraints are derived from analysis of experimental information from biophysical studies on native and mutant proteins, from analysis of the sequences of related proteins, and from theoretical considerations of protein structure. Calculations are presented for two systems. The method was validated through calculations using appropriate experimental information for bacteriorhodopsin, which produced a model structure with a root mean square (rms) deviation of 1.87 A from the structure determined by electron microscopy. Calculations are also presented using experimental and theoretical information available for bovine rhodopsin to assign the helices to a projection density map and to produce a model of bovine rhodopsin that can be used as a template for modeling other G-protein-coupled receptors.

High-resolution structural studies of the retinal--Glu113 interaction in rhodopsin

The key to understanding the reaction mechanism of rhodopsin lies in determining the structure of the retinal binding site and in defining the charge interactions between Glu113 and the retinal protonated Schiff base chromophore. We have been using 13C-NMR chemical shift data to determine the location of the Glu113 carboxyl side chain in relation to the retinal. The NMR data constrain one of the carboxylate oxygens of Glu113 to be ca. 3 A from the C12 position of the retinal with the second oxygen oriented away from the conjugated chain. A water molecule forming a hydrogen bond with the Schiff base is incorporated into the model to account for the high C = N stretching frequency [Han et al., Biophys. J., 65 (1993) 899]. In this study, we have refined the counterion position and have shown that it can reproduce the observed chemical shift data as well as the red-shifted absorption maximum of rhodopsin. Furthermore, the retinal binding site geometry derived from the NMR constraints can be readily incorporated into a recent structural model of the apoprotein.

Membrane protein structure: the contribution and potential of novel solid state NMR approaches

Alternative methods for describing molecular detail for large integral membrane proteins are required in the absence of routine crystallographic approaches. Novel solid state NMR methods, devised for the study of large molecular assemblies, are now finding applications in biological systems, including integral membrane proteins. Wild-type and genetically engineered proteins can be investigated and detailed information about side chains, prosthetic groups, ligands (e.g. drugs) and binding sites can be deduced. The molecular structure and dynamics of selected parts of the proteins are accessible by a range of different solid state NMR approaches. Inter- and intra-atomic distances can be determined rather accurately (within ångströms) and the orientation of molecular bonds (within 2 degrees) can be measured in ideal cases. Here, a brief description of the methods is given and then some specific examples described with an indication of the future potential for the approaches in studying membrane proteins. It is anticipated that this emerging NMR methodology will be more widely used in the future, not only for resolving local structure, but also for more expansive descriptions of membrane protein structure at atomic resolution.

Protein-chromophore interactions in alpha-crustacyanin, the major blue carotenoprotein from the carapace of the lobster, Homarus gammarus. A study by 13C magic angle spinning NMR

MAS (magic angle spinning) 13C NMR has been used to study protein-chromophore interactions in alpha-crustacyanin, the blue astaxanthin-binding carotenoprotein of the lobster, Homarus gammarus, reconstituted with astaxanthins labelled with 13C at the 14,14' or 15,15' positions. Two signals are seen for alpha-crustacyanin containing [14,14'-13C2]astaxanthin, shifted 6.9 and 4.0 ppm downfield from the 134.1 ppm signal of uncomplexed astaxanthin in the solid state. With alpha-crustacyanin containing [15,15'-13C2]astaxanthin, one essentially unshifted broad signal is seen. Hence binding to the protein causes a decrease in electronic charge density, providing the first experimental evidence that a charge redistribution mechanism contributes to the bathochromic shift of the astaxanthin in alpha-crustacyanin, in agreement with inferences based on resonance Raman data [Salares, et al. (1979) Biochim. Biophys. Acta 576, 176-191]. The splitting of the 14 and 14' signals provides evidence for asymmetric binding of each astaxanthin molecule by the protein.

Analyzing the red-shift characteristics of azulenic, naphthyl, other ring-fused and retinyl pigment analogs of bacteriorhodopsin

Prompted by the near infrared-absorbing properties of some of the azulenic bacteriorhodopsin (bR) analogs, we have analyzed their absorption characteristics along with 11 new related ring-fused analogs and the corresponding Schiff bases (SB) and protonated Schiff bases (PSB). The following three factors are believed to contribute to the total red shift of each of the pigment analogs (sigma RS): perturbation of the basic chromophore (SB shift, delta SB), protonation of the SB (PSB shift, PSBS) and protein perturbation (the opsin shift, OS). For each factor, effects of structural modifications were examined. For the red-shifted pigments, percent OS has been suggested as an alternate way of measuring protein perturbation. Computer-simulated chromophores provided evidence against any explanation involving altered shapes of the binding pocket as a major cause for absorption differences. Implications of the current bR results on preparation of further red-shifted bR and possible application to visual pigment analogs are discussed.

Localization of the retinal protonated Schiff base counterion in rhodopsin

Semiempirical molecular orbital calculations are combined with 13C NMR chemical shifts to localize the counterion in the retinal binding site of vertebrate rhodopsin. Charge densities along the polyene chain are calculated for an 11-cis-retinylidene protonated Schiff base (11-cis-RPSB) chromophore with 1) a chloride counterion at various distances from the Schiff base nitrogen, 2) one or two chloride counterions at different positions along the retinal chain from C10 to C15 and at the Schiff base nitrogen, and 3) a carboxylate counterion out of the retinal plane near C12. Increasing the distance of the negative counterion from the Schiff base results in an enhancement of alternating negative and positive partial charge on the even- and odd-numbered carbons, respectively, when compared to the 11-cis-RPSB chloride model compound. In contrast, the observed 13C NMR data of rhodopsin exhibit downfield chemical shifts from C8 to C13 relative to the 11-cis-RPSB.Cl corresponding to a net increase of partial positive or decrease of partial negative charge at these positions (Smith, S. O., I. Palings, M. E. Miley, J. Courtin, H. de Groot, J. Lugtenburg, R. A. Mathies, and R. G. Griffin. 1990. Biochemistry. 29:8158-8164). The anomalous changes in charge density reflected in the rhodopsin NMR chemical shifts can be qualitatively modeled by placing a single negative charge above C12. The calculated fit improves when a carboxylate counterion is used to model the retinal binding site. Inclusion of water in the model does not alter the fit to the NMR data, although it is consistent with observations based on other methods. These data constrain the location and the orientation of the Glu113 side chain, which is known to be the counterion in rhodopsin, and argue for a strong interaction centered at C12 of the retinylidene chain.

Magic angle spinning NMR studies on the metarhodopsin II intermediate of bovine rhodopsin: evidence for an unprotonated Schiff base

Magic angle spinning (MAS)13C-NMR spectra of the metarhodopsin II intermediate have been obtained using bovine rhodopsin regenerated with retinal 13C-labeled at the C-13 and C-15 positions to investigate the protonation state of the retinal Schiff base linkage. The 13C-labeled rhodopsin was reconstituted into 1,2-dipalmitoleoylphosphatidylcholine bilayers to increase the amount of meta II trapped at low temperature. Both the 13C-15 (159.2 ppm) and 13C-13 (144.0 ppm) isotropic chemical shifts are characteristic of an unprotonated Schiff base, while the 13C-15 shift is significantly different from that of retinal (191 ppm) or a tetrahedral carbinolamine group (70-90 ppm) previously proposed as an intermediate in the hydrolysis of the Schiff base at the meta II stage. This rules out the possibility that meta II non-covalently binds retinal or is a carbinolamine intermediate and provides convincing evidence that Schiff base deprotonation occurs in the meta I-meta II transition, an event that is likely to be important in triggering the activation of transducin.

9,11-Dicis-12-fluororhodopsin. Chromophore properties from 19F-NMR studies

From a 19F-NMR study of 9,11-dicis-12-fluororhodopsin and its photobleached product, we concluded that the initially formed chromophore retained its configuration and the photoproduct corresponded to the two-bond isomerized all-trans. Upon standing, it slowly isomerized to the 9-cis isomer. The method represents a direct, non-destructive procedure for determining configuration purity of the pigment formed. Its unique fluorine opsin shift value is consistent with the expected different orientation of the fluoro-substituent in a dicis pigment.

NMR studies of retinal proteins

A review is given of the use of nuclear magnetic resonance (NMR) spectroscopy to study bacteriorhodopsin and bovine rhodopsin. Solution and solid-state approaches are included. The studies of the bacterial proton pump examine the chromophore, the peptide backbone, and the protein side chains. The studies of the bovine visual pigment are limited to the chromophore. Various forms of each pigment are considered. Both structural and dynamic features are addressed.

Retinal analog restoration of photophobic responses in a blind Chlamydomonas reinhardtii mutant. Evidence for an archaebacterial like chromophore in a eukaryotic rhodopsin

The strain CC-2359 of the unicellular eukaryotic alga Chlamydomonas reinhardtii originally described as a low pigmentation mutant is found to be devoid of photophobic stop responses to photostimuli over a wide range of light intensities. Photophobic responses of the mutant are restored by exogenous addition of all-trans retinal. We have combined computer-based cell-tracking and motion analysis with retinal isomer and retinal analog reconstitution of CC-2359 to investigate properties of the photophobic response receptor. Most rapid and most complete reconstitution is obtained with all-trans retinal compared to 13-cis, 11-cis, and 9-cis retinal. An analog locked by a carbon bridge in a 6-s-trans conformation reconstitutes whereas the corresponding 6-s-cis locked analog does not. Retinal analogs prevented from isomerization around the 13-14 double bond by a five-membered ring in the polyene chain (locked in either the 13-trans or 13-cis configuration) do not restore the response, but enter the chromophore binding pocket as evidenced by their inhibition of all-trans retinal regeneration of the response. Results of competition experiments between all-trans and each of the 13-locked analogs fit a model in which each chromophore exhibits reversible binding to the photoreceptor apoprotein. A competitive inhibition scheme closely fits the data and permits calculation of apparent dissociation constants for the in vivo reconstitution process of 2.5 x 10(-11) M, 5.2 x 10(-10) M, and 5.4 x 10(-9) M, for all-trans, 13-trans-locked and 13-cis-locked analogs, respectively. The chromophore requirement for the trans configuration and 6-s-trans conformation, and the lack of signaling function from analogs locked at the 13 position, are characteristic of archaebacterial rhodopsins, rather than the previously studied eukaryotic rhodopsins (i.e., visual pigments).

NMR studies of fluorinated visual pigment analogs

The 19F-nmr chemical shift data of isomeric pigments (11-cis and 9-cis) of four vinyl fluororhodopsins and two trifluororhodopsins have been recorded. When compared with model protonated Schiff bases, a set of F-nmr opsin shift parameter (FOS) was obtained. The data revealed regiospecific protein perturbations on the F-resonances. They can be interpreted in terms of specific protein interactions such as the postulated second point charge and other polar interactions as well as the common hydrophobic protein perturbation.

Molecular determinants of visual pigment function

Mutagenesis studies and comparisons of natural variants of rhodopsin and related visual pigments have led to new insights concerning photoreceptor function. The studies identify domains important for receptor folding, the residues that set the wavelength of absorption for the ligand 11-cis retinal, and residues, that when mutated, trigger the cell death of photoreceptors.

13C magic-angle spinning NMR studies of bathorhodopsin, the primary photoproduct of rhodopsin

Magic-angle spinning NMR spectra have been obtained of the bathorhodopsin photointermediate trapped at low temperature (less than 130 K) by using isorhodopsin samples regenerated with retinal specifically 13C-labeled at positions 8, 10, 11, 12, 13, 14, and 15. Comparison of the chemical shifts of the bathorhodopsin resonances with those of an all-trans-retinal protonated Schiff base (PSB) chloride salt show the largest difference (6.2 ppm) at position 13 of the protein-bound retinal. Small differences in chemical shift between bathorhodopsin and the all-trans PSB model compound are also observed at positions 10, 11, and 12. The effects are almost equal in magnitude to those previously observed in rhodopsin and isorhodopsin. Consequently, the energy stored in the primary photoproduct bathorhodopsin does not give rise to any substantial change in the average electron density at the labeled positions. The data indicate that the electronic and structural properties of the protein environment are similar to those in rhodopsin and isorhodopsin. In particular, a previously proposed perturbation near position 13 of the retinal appears not to change its position significantly with respect to the chromophore upon isomerization. The data effectively exclude charge separation between the chromophore and a protein residue as the main mechanism for energy storage in the primary photoproduct and argue that the light energy is stored in the form of distortions of the bathorhodopsin chromophore.

Structure of the retinal chromophore in 7,9-dicis-rhodopsin

Bovine rhodopsin was bleached and regenerated with 7,9-dicis-retinal to form 7,9-dicis-rhodopsin, which was purified on a concanavalin A affinity column. The absorption maximum of the 7,9-dicis pigment is 453 nm, giving an opsin shift of 1600 cm-1 compared to 2500 cm-1 for 11-cis-rhodopsin and 2400 cm-1 for 9-cis-rhodopsin. Rapid-flow resonance Raman spectra have been obtained of 7,9-dicis-rhodopsin in H2O and D2O at room temperature. The shift of the 1654-cm-1 C = N stretch to 1627 cm-1 in D2O demonstrates that the Schiff base nitrogen is protonated. The absence of any shift in the 1201-cm-1 mode, which is assigned as the C14-C15 stretch, or of any other C-C stretching modes in D2O indicates that the Schiff base C = N configuration is trans (anti). Assuming that the cyclohexenyl ring binds with the same orientation in 7,9-dicis-, 9-cis-, and 11-cis-rhodopsins, the presence of two cis bonds requires that the N-H bond of the 7,9-dicis chromophore points in the opposite direction from that in the 9-cis or 11-cis pigment. However, the Schiff base C = NH+ stretching frequency and its D2O shift in 7,9-dicis-rhodopsin are very similar to those in 11-cis- and 9-cis-rhodopsin, indicating that the Schiff base electrostatic/hydrogen-bonding environments are effectively the same. The C = N trans (anti) Schiff base geometry of 7,9-dicis-rhodopsin and the insensitivity of its Schiff base vibrational properties to orientation are rationalized by examining the binding site specificity with molecular modeling.