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...