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

Light activation of the isomerization and deprotonation of the protonated Schiff base retinal

We perform an ab initio analysis of the photoisomerization of the protonated Schiff base of retinal (PSB-retinal) from 11-cis to 11-trans rotating the C10-C11=C12-C13 dihedral angle from 0° (cis) to -180° (trans). We find that the retinal molecule shows the lowest rotational barrier (0.22 eV) when its charge state is zero as compared to the barrier for the protonated molecule which is ∼0.89 eV. We conclude that rotation most likely takes place in the excited state of the deprotonated retinal. The addition of a proton creates a much larger barrier implying a switching behavior of retinal that might be useful for several applications in molecular electronics. All conformations of the retinal compound absorb in the green region with small shifts following the dihedral angle rotation; however, the Schiff base of retinal (SB-retinal) at trans-conformation absorbs in the violet region. The rotation of the dihedral angle around the C11=C12 π-bond affects the absorption energy of the retinal and the binding energy of the SB-retinal with the proton at the N-Schiff; the binding energy is slightly lower at the trans-SB-retinal than at other conformations of the retinal.

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.

Anabaena sensory rhodopsin is a light-driven unidirectional rotor

The implementation of multiconfigurational quantum chemistry methods into a quantum-mechanics/molecular-mechanics protocol has allowed the construction of a realistic computer model for the sensory rhodopsin of the cyanobacterium Anabaena PCC 7120. The model, which reproduces the absorption spectra of both the all-trans and 13-cis forms of the protein and their associated K and L intermediates, is employed to investigate the light-driven steps of the photochromic cycle exhibited by the protein. It is found that the photoisomerizations of the all-trans and 13-cis retinal chromophores occur through unidirectional, counterclockwise 180° rotations of the =C14-C15= moiety with respect to the Lys210-linked end of the chromophore axis. Thus, the sequential interconversions of the all-trans and 13-cis forms during a single photochromic cycle yield a complete (360°) unidirectional rotation of the =C14-C15= moiety. This finding implies that Anabaena sensory rhodopsin is a biological realization of a light-driven molecular rotor.

Aborted double bicycle-pedal isomerization with hydrogen bond breaking is the primary event of bacteriorhodopsin proton pumping

Quantum mechanics/molecular mechanics calculations based on ab initio multiconfigurational second order perturbation theory are employed to construct a computer model of Bacteriorhodopsin that reproduces the observed static and transient electronic spectra, the dipole moment changes, and the energy stored in the photocycle intermediate K. The computed reaction coordinate indicates that the isomerization of the retinal chromophore occurs via a complex motion accounting for three distinct regimes: (i) production of the excited state intermediate I, (ii) evolution of I toward a conical intersection between the excited state and the ground state, and (iii) formation of K. We show that, during stage ii, a space-saving mechanism dominated by an asynchronous double bicycle-pedal deformation of the C10═C11─C12═C13─C14═N moiety of the chromophore dominates the isomerization. On this same stage a N─H/water hydrogen bond is weakened and initiates a breaking process that is completed during stage iii.

Modeling, preparation, and characterization of a dipole moment switch driven by Z/E photoisomerization

We report the results of a multidisciplinary research effort where the methods of computational photochemistry and retrosynthetic analysis/synthesis have contributed to the preparation of a novel N-alkylated indanylidene-pyrroline Schiff base featuring an exocyclic double bond and a permanent zwitterionic head. We show that, due to its large dipole moment and efficient photoisomerization, such a system may constitute the prototype of a novel generation of electrostatic switches achieving a reversible light-induced dipole moment change on the order of 30 D. The modeling of a peptide fragment incorporating the zwitterionic head into a conformationally rigid side chain shows that the switch can effectively modulate the fluorescence of a tryptophan probe.

Multiple functions of Schiff base counterion in rhodopsins

In rhodopsins, visible-light absorption is achieved by the protonation of the chromophore Schiff base. The Schiff base proton is stabilized by the negative charge of an amino acid residue called the Schiff base counterion. Since E113 was identified as the counterion in bovine rhodopsin, there has been growing evidence that the counterion has multiple functions besides proton stabilization. Here, we first introduce generally accepted findings as well as some controversial theories about the identity of the Schiff base counterion in the dark and in intermediate states and then review multiple functions of the counterion in vertebrate visual pigments. Special focus is placed on the recently demonstrated role in photoisomerization efficiency. Finally, differences in the position of the counterion between vertebrate visual pigments and other opsins and its relevance to the molecular evolution of opsins are discussed.

Studies of the ground and excited-state surfaces of the retinal chromophore using CAM-B3LYP

The isomerization of the 11-cis isomer (PSB11) of the retinal chromophore to its all-trans isomer (PSBT) is examined. Optimized structures on both the ground state and the excited state are calculated, and the dependence on torsional angles in the carbon chain is investigated. Time-dependent density functional theory is used to produce excitation energies and the excited-state surface. To avoid problems with the description of excited states that can arise with standard DFT methods, the CAM-B3LYP functional was used. Comparing CAM-B3LYP with B3LYP results indicates that the former is significantly more accurate, as a consequence of which detailed cross sections of the retinal excited-state surface are obtained.

Modeling reaction routes from rhodopsin to bathorhodopsin

The quantum mechanical-molecular mechanical (QM/MM) theory was applied to calculate accurate structural parameters, vibrational and optical spectra of bathorhodopsin (BATHO), one of the primary photoproducts of the functional cycle of the visual pigment rhodopsin (RHO), and to characterize reaction routes from RHO to BATHO. The recently resolved crystal structure of BATHO (PDBID: 2G87) served as an initial source of coordinates of heavy atoms. Protein structures in the ground electronic state and vibrational frequencies were determined by using the density functional theory in the PBE0/cc-pVDZ approximation for the QM part and the AMBER force field parameters in the MM part. Calculated and assigned vibrational spectra of both model protein systems, BATHO and RHO, cover three main regions referring to the hydrogen-out-of-plan (HOOP) motion, the C==C ethylenic stretches, and the C--C single-bond stretches. The S(0)-S(1) electronic excitation energies of the QM part, including the chromophore group in the field of the protein matrix, were estimated by using the advanced quantum chemistry methods. The computed structural parameters as well as the spectral bands match perfectly the experimental findings. A structure of the transition state on the S(0) potential energy surface for the ground electronic state rearrangement from RHO to BATHO was located proving a possible route of the thermal protein activation to the primary photoproduct.

Color tuning in short wavelength-sensitive human and mouse visual pigments: ab initio quantum mechanics/molecular mechanics studies

We have investigated the protonation state and photoabsorption spectrum of Schiff-base (SB) nitrogen bound 11-cis-retinal in human blue and mouse UV cone visual pigments as well as in bovine rhodopsin by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have employed both multireference (MRCISD+Q, MR-SORCI+Q, and MR-DDCI2+Q) and single reference (TD-B3LYP and RI-CC2) QM methods. The calculated ground-state and vertical excitation energies show that UV-sensitive pigments have deprotonated SB nitrogen, while violet-sensitive pigments have protonated SB nitrogen, in agreement with some indirect experimental evidence. A significant blue shift of the absorption maxima of violet-sensitive pigments relative to rhodopsins arises from the increase in bond length alternation of the polyene chain of 11-cis-retinal induced by polarizing fields of these pigments. The main counterion is Glu113 in both violet-sensitive vertebrate pigments and bovine rhodopsin. Neither Glu113 nor the remaining pigment has a significant influence on the first excitation energy of 11-cis-retinal in the UV-sensitive pigments that have deprotonated SB nitrogen. There is no charge transfer between the SB and beta-ionone terminals of 11-cis-retinal in the ground and first excited states.

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

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

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.

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

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

Acceleration of the Z to E photoisomerization of penta-2,4-dieniminium by hydrogen out-of-plane motion: theoretical study on a model system of retinal protonated Schiff base

We report the result of comparison between two reaction coordinates [on the potential energy surface of the first excited state (S(1))] produced by CASSCF and these energies recalculated by MRMP2 in the Z to E photoisomerization of penta-2,4-dieniminium (PDI) as the minimal model of the retinal protonated Schiff base (RPSB). One coordinate is the S(1) state minimum-energy-path (MEP) in mass-weighted coordinates from the S(1) vertically excited point, where a strong hydrogen-out-of plane (HOOP) motion is not exhibited. The energy profile of the S(1) MEP at the MRMP2//CASSCF level shows a barrier for the rotation around the reactive C-C and hits the S(1)/S(0) degeneracy space where the central C-C-C-C dihedral angle is distorted by 65 degrees . The other coordinate is an S(1) coordinate obtained by the relaxed scan strategy. The relaxed coordinate along the central C-C-C-C dihedral angle, which we call the HOOP coordinate, shows strong HOOP motion. According to the MRMP2//CASSCF calculation, there is no barrier on the HOOP coordinate. Furthermore, the S(1) to S(0) transition may be possible without the large skeletal deformation by HOOP motion because the HOOP coordinate encounters the S(1)/S(0) degeneracy space where the central C-C-C-C dihedral angle is distorted by only 40 degrees . Consequently, if PDI is a suitable model molecule for the RPSB as often assumed, the 11-cis to all-trans photoisomerization is predicted to be accelerated by the HOOP motion.

Photochemical reaction dynamics of the primary event of vision studied by means of a hybrid molecular simulation

The photoisomerization reaction dynamics of a retinal chromophore in the visual receptor rhodopsin was investigated by means of hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations. The photoisomerization reaction of retinal constitutes the primary step of vision and is known as one of the fastest reactions in nature. To elucidate the molecular mechanism of the high efficiency of the reaction, we carried out hybrid ab initio QM/MM MD simulations of the complete reaction process from the vertically excited state to the photoproduct via electronic transition in the entire chromophore-protein complex. An ensemble of reaction trajectories reveal that the excited-state dynamics is dynamically homogeneous and synchronous even in the presence of thermal fluctuation of the protein, giving rise to the very fast formation of the photoproduct. The synchronous nature of the reaction dynamics in rhodopsin is found to originate from weak perturbation of the protein surroundings and from dynamic regulation of volume-conserving motions of the chromophore. The simulations also provide a detailed view of time-dependent modulations of hydrogen-out-of-plane vibrations during the reaction process, and identify molecular motions underlying the experimentally observed dynamic spectral modulations.

An artificial molecular switch that mimics the visual pigment and completes its photocycle in picoseconds

Single molecules that act as light-energy transducers (e.g., converting the energy of a photon into atomic-level mechanical motion) are examples of minimal molecular devices. Here, we focus on a molecular switch designed by merging a conformationally locked diarylidene skeleton with a retinal-like Schiff base and capable of mimicking, in solution, different aspects of the transduction of the visual pigment Rhodopsin. Complementary ab initio multiconfigurational quantum chemistry-based computations and time-resolved spectroscopy are used to follow the light-induced isomerization of the switch in methanol. The results show that, similar to rhodopsin, the isomerization occurs on a 0.3-ps time scale and is followed by <10-ps cooling and solvation. The entire (2-photon-powered) switch cycle was traced by following the evolution of its infrared spectrum. These measurements indicate that a full cycle can be completed within 20 ps.

Dielectric relaxation in proteins: the computational perspective

In photoexcitation and electron transfer, a new dipole or charge is introduced, and the structure is adjusted. This adjustment represents dielectric relaxation, which is the focus of this review. We concentrate on a few selected topics. We discuss linear response theory, as a unifying framework and a tool to describe non-equilibrium states. We review recent, molecular dynamics simulation studies that illustrate the calculation of dynamic and thermodynamic properties, such as Stokes shifts or reorganization free energies. We then turn to the macroscopic, continuum electrostatic view. We recall the physical definition of a dielectric constant and revisit the decomposition of the free energy into a reorganization and a static term. We review some illustrative continuum studies and discuss some difficulties that can arise with the continuum approach. In conclusion, we consider recent developments that will increase the accuracy and broaden the scope of all these methods.

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.

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.

Spectral tuning of photoactive yellow protein

We report a theoretical study on the optical properties of a small, water-soluble photosensory receptor, photoactive yellow protein (PYP). A hierarchical ab initio molecular orbital calculation accurately evaluated the optical absorption maximum of the wild-type, as well as the lambda(max) values of 12 mutants. Electronic excitation of the chromophore directly affects the electronic state of nearby atoms in the protein environment. This effect is explicitly considered in the present study. Furthermore, the spectral tuning mechanism of PYP was investigated at the atomic level. The static disorder of a protein molecule is intimately related to the complex nature of its energy landscape. By using molecular dynamics simulation and quantum mechanical structure optimization, we obtained multiple minimum energy conformations of PYP. The statistical distribution of electronic excitation energies of these minima was compared with the hole-burning experiment (Masciangioli, T. [2000] Photochem. Photobiol. 72, 639), a direct observation of the distribution of excitation energies.

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

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

Stairway to the conical intersection: a computational study of the retinal isomerization

The potential-energy surface of the first excited state of the 11-cis-retinal protonated Schiff base (PSB11) chromophore has been studied at the density functional theory (DFT) level using the time-dependent perturbation theory approach (TDDFT) in combination with Becke's three-parameter hybrid functional (B3LYP). The potential-energy curves for torsion motions around single and double bonds of the first excited state have also been studied at the coupled-cluster approximate singles and doubles (CC2) level. The corresponding potential-energy curves for the ground state have been calculated at the B3LYP DFT and second-order Møller-Plesset (MP2) levels. The TDDFT study suggests that the electronic excitation initiates a turn of the beta-ionone ring around the C6-C7 bond. The torsion is propagating along the retinyl chain toward the cis to trans isomerization center at the C11=C12 double bond. The torsion twist of the C10-C11 single bond leads to a significant reduction in the deexcitation energy indicating that a conical intersection is being reached by an almost barrierless rotation around the C10-C11 single bond. The energy released when passing the conical intersection can assist the subsequent cis to trans isomerization of the C11=C12 double bond. The CC2 calculations also show that the torsion barrier for the twist of the retinyl C10-C11 single bond adjacent to the isomerization center almost vanishes for the excited state. Because of the reduced torsion barriers of the single bonds, the retinyl chain can easily deform in the excited state. Thus, the CC2 and TDDFT calculations suggest similar reaction pathways on the potential-energy surface of the excited state leading toward the conical intersection and resulting in a cis to trans isomerization of the retinal chromophore. According to the CC2 calculations the cis to trans isomerization mechanism does not involve any significant torsion motion of the beta-ionone ring.

Tracking the excited-state time evolution of the visual pigment with multiconfigurational quantum chemistry

The primary event that initiates vision is the photoinduced isomerization of retinal in the visual pigment rhodopsin (Rh). Here, we use a scaled quantum mechanics/molecular mechanics potential that reproduces the isomerization path determined with multiconfigurational perturbation theory to follow the excited-state evolution of bovine Rh. The analysis of a 140-fs trajectory provides a description of the electronic and geometrical changes that prepare the system for decay to the ground state. The data uncover a complex change of the retinal backbone that, at approximately 60-fs delay, initiates a space saving "asynchronous bicycle-pedal or crankshaft" motion, leading to a conical intersection on a 110-fs time scale. It is shown that the twisted structure achieved at decay features a momentum that provides a natural route toward the photoRh structure recently resolved by using femtosecond-stimulated Raman spectroscopy.

Femtosecond Stimulated Raman Spectroscopy

Femtosecond stimulated Raman spectroscopy (FSRS) is a new ultrafast spectroscopic technique that provides vibrational structural information with high temporal (50-fs) and spectral (10-cm(1)) resolution. As a result of these unique capabilities, FSRS studies of chemical and biochemical reaction dynamics are expected to grow rapidly, giving previously unattainable insight into the structural dynamics of reactively evolving systems with atomic spatial and femtosecond temporal resolution. This review discusses the experimental and theoretical concepts behind FSRS, with an emphasis on the origins of its unique temporal and spectral capabilities. We illustrate these capabilities with vibrational studies of ultrafast electronic dynamics, as well as the direct structural observation of nonstationary vibrational wave-packet motion in small molecules and in complex biochemical reaction dynamics.

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.

Photochemistry of Visual Pigment Chromophore Models by Ab Initio Molecular Dynamics

Ab initio excited-state molecular dynamics calculations have been performed to study the effect of methyl substitution and chromophore distortion on the photoreaction of different four-double-bond retinal model chromophores. Randomly distributed starting geometries were generated by zero-point energy sampling; after Franck-Condon excitation the reaction was followed on the S1 surface. For determining the photoproduct and its configuration, a simplified approach--torsion angle following--is discussed and applied. We find that chromophore distortion significantly affects the outcome of the photoreaction: with dihedral angles taken from the rhodopsin-embedded 11-cis-retinal chromophore, the reaction rate of the model chromophore is increased by a factor of 3 compared to that of the relaxed chromophore. Also, the reaction proceeds in a completely stereoselective manner involving only the cis double bond and with a minimum quantum yield of 72%. Bond torsion is more effective than methyl substitution for fast and selective photochemistry, which is in agreement with photophysical measurements on rhodopsin analogues. We conclude that apart from the geometric distortions caused by the protein pocket it is not necessary to postulate other specific interactions between the protein and the chromophore to effect the selective and ultrafast photoreaction in rhodopsin.