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 Ϸ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.photoisomerization ͉ rhodopsin ͉ vision T he visual pigment rhodopsin (Rh) (1, 2) is a G protein-coupled receptor containing a 11-cis retinal chromophore (PSB11) bounded to a lysine residue (Lys-296) via a protonated Schiff base linkage (see Fig. 1). While the biological activity of Rh is triggered by the light-induced 11-cis all-trans isomerization of PSB11, this reaction owes its efficiency (e.g., short time scale and high quantum yields) to the protein cavity (1). Recently, the mechanism of the isomerization of retinal in Rh has been investigated by using femtosecond-stimulated Raman spectroscopy (FSRS) (3). Kukura et al. (3) have reported on experimentally derived structures of photoRh and bathoRh, namely the first and second ground-state intermediates of the Rh photocycle.While such progress has provided information on the structural changes achieved 200 fs after light absorption, the faster structural changes prompting the excited-state decay of PSB11 (i.e., the central event of the isomerization mechanism) remain to be established. Indeed, it has been suggested that such decay may occur on a 60-fs time scale through fast hydrogen out-ofplane (HOOP) motion (3), whereas the traditional view points to a slower Ϸ150-fs decay driven by cis-trans isomerization motion (4). In principle, molecular dynamics simulations featuring a quantum chemical description of the chromophore can be used to address such issues. This fact was shown by Warshel (5) using semiempirical quantum chemistry to describe PSB11 and geometrical constraints to account for the protein environment. Later, Birge and Hubbard (6) reported a different semiempirical study of an explicit chromophore-counterion pair evolving along a single coordinate. While the first simulation of the retinal photoisomerization using a full atomic-level protein model (7) was reported for the related receptor bacterio-Rh (bR), attempts to simulate the PSB11 excited-state motion in a complete Rh model are more recent (8-10). On the other hand, a quantitative evaluation of the isomerization coordinate and time scale requires, as a prerequisite, an accurate excited-st...

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 excitedstate 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 Ϸ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.photoisomerization ͉ quantum mechanics ͉ molecular mechanics ͉ retinal ͉ vision T he visual pigment rhodopsin (Rh) (1, 2) is a G proteincoupled receptor containing an 11-cis retinal chromophore (PSB11) bounded to a lysine residue (Lys-296) via a protonated Schiff base linkage (see Scheme 1). While the biological activity of Rh is triggered by the light-induced 11-cis 3 all-trans isomerization of PSB11, this reaction owes its efficiency (e.g., short time scale and quantum yields) to the protein cavity (1). Accordingly, investigation of the environment-dependent properties of PSB11 is a prerequisite for understanding the Rh ''catalytic'' effect. The equilibrium geometry, absorption maxima ( max a ), and fluorescence maxima ( max f ) are indicators of the environment effect. In fact, whereas the geometry of PSB11 is nearly planar in a crystal (3), in bovine Rh it has a helical conformation (4). Similarly, the 445-nm max a observed for PSB11 in methanol (5) is red-shifted to 498 nm in Rh (1, 2): an effect known as the opsin shift.The Rh fluorescence band ranges from 530 to 780 nm (6). The max f has been reported (6) to be excitation wavelength-dependent, shifting from 595 to 704 nm when the excitation wavelength is shifted from 472 to 568 nm. This observation is consistent with the idea that the emission arises from a nonstationary unrelaxed excited-state population. In methanol solution the PSB11...

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

A "brute-force" ab initio CASPT2//CASSCF/CHARMM computational approach is employed to investigate the properties of the emitting state of the wild-type green fluorescence protein. The results indicate that the emitting moiety corresponds to a slightly perturbed H2O- - -chromophore complex. Thus, the protein matrix seems to be designed in such a way to mimic an environment that is more similar to gas-phase than water solution.