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.photoisomerization | quantum mechanics/molecular mechanics (QM/MM) | retinal proteins T he light-activated proton pump bacteriorhodopsin (bR) is an Archaea receptor contained in the purple membrane of Halobacterium salinarium. As shown in Scheme 1, upon photoexcitation, the all-trans retinal chromophore (PSBT) bounded to Lys216 via a protonated Schiff base linkage is converted to the 13-cis isomer (PSB13). This event triggers a series of conformational changes that ultimately result in a proton translocation from the cytoplasmic to the extracellular domain (1).The bR protein environment affects both the spectroscopic and photochemical properties of PSBT (2-4). First, the absorption maximum (568 nm) is red-shifted with respect to the one observed in solution (440 nm in methanol). Second, timeresolved spectroscopy reveals a dominant excited state lifetime component of 450 fs (5, 6) almost matching the 500-fs time scale for formation of the vibrationally hot primary photoproduct J (7,8). In contrast, PSBT in solution features a biexponential decay dynamics with a dominant (ca. 10-fold longer) 2-ps component (9-11). Finally, although bR photoisomerizes stereoselectively (leading exclusively to PSB13) with high (ca. 67%) quantum yield (9, 12), irradiation of PSBT in solution leads to a mixture of different stereoisomers with a smaller (ca. 25%) total quantum yield (9, 13).Low-temperature spectroscopic studies provided evidence for the existence of a tiny (≤1 kcal mol −1 ) energy barrier on the excited state (S 1 ) potential energy surface of bR (14,15). Such a barrier would explain the 450-fs short-lived quasistationary state observed by Ruhman et al. (5) that is assigned to the fluorescent state I (16-18) and precedes decay to the ground state (S 0 ). Recent experiments by Léonard et al. (19) have also revealed that the photoinduced changes in the permanent dipole moment of PSBT are mirrored by the absorption of Trp86 up to the first stable (i.e., cryogenically trapped) photoproduct K (7). Relative to bR, K bears 11.7-...