The retinal chromophore/chloride ion pair: structure of the photoisomerization path and interplay of charge transfer and covalent states

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.

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.

Computational Photobiology and Beyond

In this paper we review the results of a group of computational studies of the spectroscopy and photochemistry of light-responsive proteins. We focus on the use of quantum mechanics/molecular mechanics protocols based on a multiconfigurational quantum chemical treatment. More specifically, we discuss the use, limitations, and application of the ab initio CASPT2//CASSCF protocol that, presently, constitutes the method of choice for the investigation of excited state organic molecules, most notably, biological chromophores and fluorophores. At the end of this Review we will also see how the computational investigation of the visual photoreceptor rhodopsin is providing the basis for the design of light-driven artificial molecular devices.

Evidence of disruption of conjugation involving delta bonds in intramolecular electronic coupling

A dimer of dimers containing two quadruply bonded [Mo(2)(DAniF)(3)](+) units (DAniF = N,N'-di(p-anisyl)formamidinate) linked by the S-donor linker, dimethyldithiooxamidate was synthesized, structurally characterized, and electronic communication was probed. The core of [Mo(2)(DAniF)(3)](2)(C(2)S(2)N(2)Me(2)), 1, formed by the Mo(2)NSC(2)SNMo(2) atoms shows two fused but non planar six-membered rings, which differs from that of the beta form of dimethyloxamidate analogue that has a heteronaphthalene-type structure (Cotton, F. A.; Liu, C. Y.; Murillo, C. A.; Villagran, D.; Wang, X. J. Am. Chem. Soc. 2004, 126, 14822). For these two analogous compounds electronic coupling between the two [Mo(2)] units, as determined by electrochemical measurements, diminishes considerably upon replacement of O-donor by S-donor atoms (DeltaE(1/2) = 531 mV and 440 mV, respectively). This suggests that the non planar conformation of the linker in 1 hampers a pathway leading to pi conjugation. Density functional theory (DFT) calculations show that the highest occupied molecular orbitals HOMO-HOMO-1 energy gap of 0.12 eV for 1 is much smaller than that of 0.61 eV for the O-donor analogue, which is consistent with the electrochemical data.

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.

Pump-shaped dump optimal control reveals the nuclear reaction pathway of isomerization of a photoexcited cyanine dye

Using optimal control as a spectroscopic tool we decipher the details of the molecular dynamics of the essential multidimensional excited-state photoisomerization - a fundamental chemical reaction of key importance in biology. Two distinct nuclear motions are identified in addition to the overall bond-twisting motion: Initially, the reaction is dominated by motion perpendicular to the torsion coordinate. At later times, a second optically active vibration drives the system along the reaction path to the bottom of the excited-state potential. The time scales of the wavepacket motion on a different part of the excited-state potential are detailed by pump-shaped dump optimal control. This technique offers new means to control a chemical reaction far from the Franck-Condon point of absorption and to map details of excited-state reaction pathways revealing unique insights into the underlying reaction mechanism.

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.

Isomerization Through Conical Intersections

The standard model for photoinduced cis-trans isomerization about carbon double bonds is framed in terms of two electronic states and a one-dimensional reaction coordinate. We review recent work that suggests that a minimal picture of the reaction mechanism requires the consideration of at least two molecular coordinates and three electronic states. In this chapter, we emphasize the role of conical intersections and charge transfer in the photoisomerization mechanism.

S1 and S2 Excited States of Gas-Phase Schiff-Base Retinal Chromophores: A Time-Dependent Density Functional Theoretical Investigation

In concert with the recent photoabsorption experiments of gas-phase Schiff-base retinal chromophores (Nielsen et al. Phys. Rev. Lett. 2006, 96, 018304), quantum chemical calculations using time-dependent density functional theory coupled with different functionals and under the Tamm-Dancoff approximation were made on the first two excited states (S1 and S2) of two retinal chromophores: 11-cis and all-trans protonated Schiff bases. The calculated vertical excitation energies (Tv) and oscillator strengths (f) are consistent with the experimental absorption bands. The experimentally observed phenomenon that the transition dipole moment (mu) of S2 is much smaller that of S1 was interpreted by 3D representation of transition densities. The different optical behaviors (linear and nonlinear optical responds) of the excited states were investigated by considering different strengths of external electric fields.

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Photochemistry of a Retinal Protonated Schiff-Base Analogue Mimicking the Opsin Shift of Bacteriorhodopsin

A retinal Schiff base analogue which artificially mimics the protein-induced red shifting of absorption in bacteriorhodopsin (BR) has been investigated with femtosecond multichannel pump probe spectroscopy. The objective is to determine if the catalysis of retinal internal conversion in the native protein BR, which absorbs at 570 nm, is directly correlated with the protein-induced Stokes shifting of this absorption band otherwise known as the "opsin shift". Results demonstrate that the red shift afforded in the model system does not hasten internal conversion relative to that taking place in a free retinal-protonated Schiff base (RPSB) in methanol solution, and stimulated emission takes place with biexponential kinetics and characteristic timescales of approximately 2 and 10.5 ps. This shows that interactions between the prosthetic group and the protein that lead to the opsin shift in BR are not directly involved in reducing the excited-state lifetime by nearly an order of magnitude. A sub-picosecond phase of spectral evolution, analogues of which are detected in photoexcited retinal proteins and RPSBs in solution, is observed after excitation anywhere within the intense visible absorption band. It consists of a large and discontinuous spectral shift in excited-state absorption and is assigned to electronic relaxation between excited states, a scenario which might also be relevant to those systems as well. Finally, a transient excess bleach component that tunes with the excitation wavelength is detected in the data and tentatively assigned to inhomogeneous broadening in the ground state absorption band. Possible sources of such inhomogeneity and its relevance to native RPSB photochemistry are discussed.

Coupled-cluster studies of the lowest excited states of the 11-cis-retinal chromophore

The first few excited states of the 11-cis-retinal (PSB11) chromophore have been studied at the coupled-cluster approximative singles and doubles (CC2) level using triple-zeta quality basis sets augmented with double sets of polarisation functions. The two lowest vertical excitation energies of 2.14 and 3.21 eV are in good agreement with recently reported experimental values of 2.03 and 3.18 eV obtained in molecular beam measurements. Calculations at the time-dependent density functional theory (TDDFT) level using the B3LYP hybrid functional yield vertical excitation energies of 2.34 and 3.10 eV for the two lowest states. Zero-point vibrational energy (ZPVE) corrections of -0.09 and -0.17 eV were deduced from the harmonic vibrational frequencies for the ground and excited states calculated at the density functional theory (DFT) and TDDFT level, respectively, using the B3LYP hybrid functional.

The Role of the β-Ionone Ring in the Photochemical Reaction of Rhodopsin

Time-dependent density functional theory (TDDFT) calculations on the photoabsorption process of the 11-cis retinal protonated Schiff base (PSB) chromophore show that the Franck-Condon relaxation of the first excited state of the chromophore involves a torsional twist motion of the beta-ionone ring relative to the conjugated retinyl chain. For the ground state, the beta-ionone ring and the retinyl chain of the free retinal PSB chromophore form a -40 degrees dihedral angle as compared to -94 degrees for the first excited state. The double bonds of the retinal are shorter for the fully optimized structure of the excited state than for the ground state suggesting a higher cis-trans isomerization barrier for the excited state than for the ground state. According to the present TDDFT calculations, the excitation of the retinal PSB chromophore does not primarily lead to a reaction along the cis-trans torsional coordinate at the C11-C12 bond. The activation of the isomerization center seems to occur at a later stage of the photo reaction. The results obtained at the TDDFT level are supported by second-order Møller-Plesset (MP2) and approximate singles and doubles-coupled cluster (CC2) calculations on retinal chromophore models; the MP2 and CC2 calculations yield for them qualitatively the same ground state and excited-state structures as obtained in the density functional theory and TDDFT calculations.

Quantum Chemical Methods for the Investigation of Photoinitiated Processes in Biological Systems: Theory and Applications

With the advent of modern computers and advances in the development of efficient quantum chemical computer codes, the meaningful computation of large molecular systems at a quantum mechanical level became feasible. Recent experimental effort to understand photoinitiated processes in biological systems, for instance photosynthesis or vision, at a molecular level also triggered theoretical investigations in this field. In this Minireview, standard quantum chemical methods are presented that are applicable and recently used for the calculation of excited states of photoinitiated processes in biological molecular systems. These methods comprise configuration interaction singles, the complete active space self-consistent field method, and time-dependent density functional theory and its variants. Semiempirical approaches are also covered. Their basic theoretical concepts and mathematical equations are briefly outlined, and their properties and limitations are discussed. Recent successful applications of the methods to photoinitiated processes in biological systems are described and theoretical tools for the analysis of excited states are presented.

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.

Ion pairs of crystal violet in sodium bis(2-ethylhexyl)sulfosuccinate reverse micelles

The interfacial localization and the ion pair formation of the positively charged dye crystal violet (CV) in sodium bis(2-ethylhexyl)sulfosuccinate reverse micelles (AOT RMs) were studied by several structural and spectroscopic techniques and by quantum chemical calculations. The size and shape of the AOT RMs in the presence of CV were investigated by small-angle X-ray scattering, showing that CV does not significantly change the RM structure. CV localization as a function of the water to surfactant molar ratio (w(0)) was characterized by H(1) and (13)C NMR, indicating the close proximity of CV to the sulfosuccinate group of AOT at small and large w(0) values. These results were confirmed by calculation of magnetic shielding constants using the gauge-independent atomic orbital method with the HF/6-31G(d) basis set. Two different types of ion pairs between AOT and CV, i.e., contact ion pair (CIPs) and solvent-separated ion pair (SSIPs), were characterized by UV-vis spectroscopy and quantum chemical calculations using the semiempirical ZINDO-CI method. In nonpolar isotropic solvents CIPs are formed with an association constant (K(ASSOC)) of 2 x 10(4) mol(-1) L in isooctane and 750 mol(-1) L in chloroform. In AOT RMs at low w(0), CV-AOT CIPs are also formed. By increasing w(0), there is a sharp decrease in the CIP association free energy, and SSIPs are formed. (CV(+))(H(2)O)(AOT(-)) SSIPs are stable in the AOT RM up to the largest w(0) tested (w(0) = 33).

Methyl Substituents at the 11 or 12 Position of Retinal Profoundly and Differentially Affect Photochemistry and Signalling Activity of Rhodopsin

The C-11=C-12 double bond of the retinylidene chromophore of rhodopsin holds a central position in its light-induced photoisomerization and hence the photosensory function of this visual pigment. To probe the local environment of the HC-11=C-12H element we have prepared the 11-methyl and 12-methyl derivatives of 11-Z retinal and incorporated these into opsin to generate the rhodopsin analogs 11-methyl and 12-methyl rhodopsin. These analog pigments form with much slower kinetics and lower efficiency than the native pigment. The initial photochemistry and the signaling activity of the analog pigments were investigated by UV-vis and FTIR spectroscopy, and by a G protein activation assay. Our data indicate that the ultrafast formation of the first photointermediate is strongly perturbed by the presence of an 11-methyl substituent, but much less by a 12-methyl substituent. These results support the current concept of the mechanism of the primary photoisomerization event in rhodopsin. An important stronghold of this concept is an out-of-plane movement of the C-12H element, which is facilitated by torsion as well as extended positive charge delocalization into the C-10-C-13 segment of the chromophore. We argue that this mechanism is maintained principally with a methyl substituent at C-12. In addition, we show that both an 11-methyl and a 12-methyl substitutent perturb the photointermediate cascade and finally yield a low-activity state of the receptor. The 11-methyl pigment retains about 30% of the G protein activation rate of native rhodopsin, while the 12-methyl chromophore behaves like an inverse agonist up to at least 20 degrees C, trapping the protein in a perturbed Meta-I-like conformation. We conclude that the isomerization region of the chromophore and the spatial structure of the binding site are finely tuned, in order to achieve a high photosensory potential with an efficient pathway to a high-activity state.

Excited-State Charge Transfer at a Conical Intersection: Effects of an Environment

The influence of a polar and polarizable environment on charge transfer processes at a conical intersection (CI) can be described by a diabatic free energy model yielding coupled surfaces as a function of both molecular coordinates and a solvent coordinate. We extend and apply this model for the S1-S0 CI in protonated Schiff bases, representing a model for retinal isomerization (Faraday Discuss. 2004, 127, 395, 2004). A dielectric continuum description of the solvent is combined with a minimal, two-electron-two-orbital electronic structure model according to Bonacić-Koutecký, Koutecký, and Michl (Angew. Chem. 1987, 26, 170), which characterizes the charge translocation effects at the CI. The model predicts that the nonequilibrium solvent state resulting from the S0-->S1 Franck-Condon transition can entail the disappearance of the CI, such that solvent motion is necessary to reach the CI seam. The concerted evolution of the intramolecular coordinates and the solvent coordinate is illustrated by an excited-state minimum energy path.

First steps of retinal photoisomerization in proteorhodopsin

The early steps (<1 ns) in the photocycle of the detergent solubilized proton pump proteorhodopsin are analyzed by ultrafast spectroscopic techniques. A comparison to the first primary events in reconstituted proteorhodopsin as well as to the well known archaeal proton pump bacteriorhodopsin is given. A dynamic Stokes shift observed in fs-time-resolved fluorescence experiments allows a direct observation of early motions on the excited state potential energy surface. The initial dynamics is dominated by sequentially emerging stretching (<150 fs) and torsional (approximately 300 fs) modes of the retinal. The different protonation states of the primary proton acceptor Asp-97 drastically affect the reaction rate and the overall quantum efficiencies of the isomerization reactions, mainly evidenced for time scales above 1 ps. However, no major influence on the fast time scales (approximately 150 fs) could be seen, indicating that the movement out of the Franck-Condon region is fairly robust to electrostatic changes in the retinal binding pocket. Based on fs-time-resolved absorption and fluorescence spectra, ground and exited state contributions can be disentangled and allow to construct a reaction model that consistently explains pH-dependent effects in solubilized and reconstituted proteorhodopsin.

S1 and S2 excited States of gas-phase Schiff-base retinal chromophores

Photoabsorption studies of 11-cis and all-trans Schiff-base retinal chromophore cations in the gas phase have been performed at the electrostatic ion storage ring in Aarhus. A broad absorption band due to the optically allowed excitation to the first electronically excited singlet state (S1) is observed at around 600 nm. A second "dark" excited state (S2) just below 400 nm is reported for the first time. It is located approximately 1.2 eV above S1 for both chromophores. The S2 state was not visible in a solution measurement where only one highly blueshifted absorption band corresponding to the first excited state was visible. Knowledge of the position of the excited states in retinal is essential for the understanding of the fast photoisomerization in, for example, visual pigments.