Ground and excited states of retinal schiff base chromophores by multiconfigurational perturbation theory

Ultrafast terahertz Stark spectroscopy reveals the excited-state dipole moments of retinal in bacteriorhodopsin

The photoinduced all-trans to 13-cis isomerization of the retinal Schiff base represents the ultrafast first step in the reaction cycle of bacteriorhodopsin (BR). Extensive experimental and theoretical work has addressed excited-state dynamics and isomerization via a conical intersection with the ground state. In conflicting molecular pictures, the excited state potential energy surface has been modeled as a pure S[Formula: see text] state that intersects with the ground state, or in a 3-state picture involving the S[Formula: see text] and S[Formula: see text] states. Here, the photoexcited system passes two crossing regions to return to the ground state. The electric dipole moment of the Schiff base in the S[Formula: see text] and S[Formula: see text] state differs strongly and, thus, its measurement allows for assessing the character of the excited-state potential. We apply the method of ultrafast terahertz (THz) Stark spectroscopy to measure electric dipole changes of wild-type BR and a BR D85T mutant upon electronic excitation. A fully reversible transient broadening and spectral shift of electronic absorption is induced by a picosecond THz field of several megavolts/cm and mapped by a 120-fs optical probe pulse. For both BR variants, we derive a moderate electric dipole change of 5 [Formula: see text] 1 Debye, which is markedly smaller than predicted for a neat S[Formula: see text]-character of the excited state. In contrast, S[Formula: see text]-admixture and temporal averaging of excited-state dynamics over the probe pulse duration gives a dipole change in line with experiment. Our results support a picture of electronic and nuclear dynamics governed by the interaction of S[Formula: see text] and S[Formula: see text] states in a 3-state model.

Multistate Multiconfiguration Quantum Chemical Computation of the Two-Photon Absorption Spectra of Bovine Rhodopsin

Recently, progress in IR sources has led to the discovery that humans can detect infrared (IR) light. This is hypothesized to be due to the two-photon absorption (TPA) events promoting the retina dim-light rod photoreceptor rhodopsin to the same excited state populated via one-photon absorption (OPA). Here, we combine quantum mechanics/molecular mechanics and extended multiconfiguration quasi-degenerate perturbation theory calculations to simulate the TPA spectrum of bovine rhodopsin (Rh) as a model for the human photoreceptor. The results show that the TPA spectrum of Rh has an intense S₀ → S₁ band but shows also S₀ → S₂ and S₀ → S₃ transitions whose intensities, relative to the S₀ → S₁ band, are significantly increased when compared to the corresponding bands of the OPA spectrum. In conclusion, we show that IR light in the 950 nm region can be perceived by rod photoreceptors, thus supporting the two-photon origin of the IR perception. We also found that the same photoreceptor can perceive red (i.e., close to 680 nm) light provided that TPA induces population of S₂.

pH-dependent absorption spectra of rhodopsin mutant E113Q: On the role of counterions and protein

The absorption spectra of bovine rhodopsin mutant E113Q in solutions were investigated at the molecular level by using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The calculations suggest the mechanism of the absorption variations of E113Q at different pH values. The results indicate that the polarizations of the counterions in the vicinity of Schiff base under protonation and unprotonation states of the mutant E113Q would be a crucial factor to change the energy gap of the retinal to tune the absorption spectra. Glu-181 residue, which is close to the chromophore, cannot serve as the counterion of the protonated Schiff base of E113Q in dark state. Moreover, the results of the absorption maximum in mutant E113Q with the various anions (Cl^-, Br^-, I^- and NO₃^-) manifested that the mutant E113Q could have the potential for use as a template of anion biosensors at visible wavelength.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Protein Field Effect on the Dark State of 11-cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics

The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.

Gas-Phase Retinal Spectroscopy: Temperature Effects Are But a Mirage

We employ state-of-the-art first-principle approaches to investigate whether temperature effects are responsible for the unusually broad and flat spectrum of protonated Schiff base retinal observed in photodissociation spectroscopy, as has recently been proposed. We first carefully calibrate how to construct a realistic geometrical model of retinal and show that the exchange-correlation M06-2X functional yields an accurate description while the commonly used complete active space self-consistent field method (CASSCF) is not adequate. Using modern multiconfigurational perturbative methods (NEVPT2) to compute the excitations, we then demonstrate that conformations with different orientations of the β-ionone ring are characterized by similar excitations. Moreover, other degrees of freedom identified as active in room-temperature molecular dynamics simulations do not yield the shift required to explain the anomalous spectral shape. Our findings indicate that photodissociation experiments are not representative of the optical spectrum of retinal in the gas phase and call for further experimental characterization of the dissociation spectra.

The effect of protein environment on photoexcitation properties of retinal

Retinal is the photon absorbing chromophore of rhodopsin and other visual pigments, enabling the vertebrate vision process. The effects of the protein environment on the primary photoexcitation process of retinal were studied by time-dependent density functional theory (TDDFT) and the algebraic diagrammatic construction through second order (ADC(2)) combined with our recently introduced reduction of virtual space (RVS) approximation method. The calculations were performed on large full quantum chemical cluster models of the bluecone (BC) and rhodopsin (Rh) pigments with 165-171 atoms. Absorption wavelengths of 441 and 491 nm were obtained at the B3LYP level of theory for the respective models, which agree well with the experimental values of 414 and 498 nm. Electrostatic rather than structural strain effects were shown to dominate the spectral tuning properties of the surrounding protein. The Schiff base retinal and a neighboring Glu-113 residue were found to have comparable proton affinities in the ground state of the BC model, whereas in the excited state, the proton affinity of the Schiff base is 5.9 kcal/mol (0.26 eV) higher. For the ground and excited states of the Rh model, the proton affinity of the Schiff base is 3.2 kcal/mol (0.14 eV) and 7.9 kcal/mol (0.34 eV) higher than for Glu-113, respectively. The protein environment was found to enhance the bond length alternation (BLA) of the retinyl chain and blueshift the first absorption maxima of the protonated Schiff base in the BC and Rh models relative to the chromophore in the gas phase. The protein environment was also found to decrease the intensity of the second excited state, thus improving the quantum yield of the photoexcitation process. Relaxation of the BC model on the excited state potential energy surface led to a vanishing BLA around the isomerization center of the conjugated retinyl chain, rendering the retinal accessible for cis-trans isomerization. The energy of the relaxed excited state was found to be 30 kcal/mol (1.3 eV) above the minimum ground state energy, and might be related to the transition state of the thermal activation process.

Quantum mechanical/molecular mechanical structure, enantioselectivity, and spectroscopy of hydroxyretinals and insights into the evolution of color vision in small white butterflies

Since Vogt's discovery of A(3)-retinal or 3-hydroxyretinal in insects in 1983 and Matsui's discovery of A(4)-retinal or 4-hydroxyretinal in firefly squid in 1988, hydroxyretinal-protein interactions mediating vision have remained largely unexplored. In the present study, A(3)- and A(4)-retinals are theoretically incorporated into squid and bovine visual pigments by use of the hybrid quantum mechanics/molecular mechanics [SORCI+Q//B3LYP/6-31G(d):Amber96] method, and insights into structure, enantioselectivity, and spectroscopy are gathered and presented for the first time. Contrary to general perception, our findings rule out the formation of a hydrogen bond between the hydroxyl-bearing β-ionone ring portion of retinal and opsin. Compared to A(1)-pigments, A(3)- and A(4)-pigments exhibit slightly blue-shifted absorption maxima due to increase in bond-length alternation of the hydroxyretinal. We suggest that (i) the binding site of firefly squid (Watasenia scintillans) opsin is very similar to that of the Japanese common squid (Todarodes pacificus) opsin; (ii) the molecular mechanism of spectral tuning in small white butterflies involve sites S116 and T185 and breaking of a hydrogen bond between sites E180 and T185; and finally (iii) A(3)-retinal may have occurred during the conversion of A(1)- to A(2)-retinal and insects may have acquired them, in order to absorb light in the blue-green wavelength region and to speed up the G-protein signaling cascade.

The opsin shift and mechanism of spectral tuning in rhodopsin

Molecular dynamics simulations and combined quantum mechanical and molecular mechanical calculations have been performed to investigate the mechanism of the opsin shift and spectral tuning in rhodopsin. A red shift of -980 cm(-1) was estimated in the transfer of the chromophore from methanol solution environment to the protonated Schiff base (PSB)-binding site of the opsin. The conformational change from a 6-s-cis-all-trans configuration in solution to the 6-s-cis-11-cis conformer contributes additional -200 cm(-1), and the remaining effects were attributed to dispersion interactions with the aromatic residues in the binding site. An opsin shift of 2100 cm(-1) was obtained, in reasonable accord with experiment (2730 cm(-1)). Dynamics simulations revealed that the 6-s-cis bond can occupy two main conformations for the β-ionone ring, resulting in a weighted average dihedral angle of about -50°, which may be compared with the experimental estimate of -28° from solid-state NMR and Raman data. We investigated a series of four single mutations, including E113D, A292S, T118A, and A269T, which are located near the PSB, along the polyene chain of retinal and close to the ionone ring. The computational results on absorption energy shift provided insights into the mechanism of spectral tuning, which involves all means of electronic structural effects, including the stabilization or destabilization of either the ground or the electronically excited state of the retinal PSB.

QM/MM study of the structure, energy storage, and origin of the bathochromic shift in vertebrate and invertebrate bathorhodopsins

By comparing the results from a hybrid quantum mechanics/molecular mechanics method (SORCI+Q//B3LYP/6-31G*:Amber) between vertebrate (bovine) and invertebrate (squid) visual pigments, the mechanism of molecular rearrangements, energy storage, and origin of the bathochromic shift accompanying the transformation of rhodopsin to bathorhodopsin have been evaluated. The analysis reveals that, in the presence of an unrelaxed binding site, bathorhodopsin was found to carry almost 27 kcal/mol energy in both visual pigments and absorb (λ(max)) at 528 nm in bovine and 554 nm in squid. However, when the residues within 4.0 Å radius of the retinal are relaxed during the isomerization event, almost ∼16 kcal/mol energy is lost in squid compared to only ∼8 kcal/mol in bovine. Loss of a larger amount of energy in squid is attributed to the presence of a flexible binding site compared to a rigid binding site in bovine. Structure of the squid bathorhodopsin is characterized by formation of a direct H-bond between the Schiff base and Asn87.

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.

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.

Modeling the Excited States of Biological Chromophores within Many-Body Green's Function Theory

First-principle many-body Green's function theory (MBGFT) has been successfully used to describe electronic excitations in many materials, from bulk crystals to nanoparticles. Here we assess its performance for the calculations of the excited states of biological chromophores. MBGFT is based on a set of Green's function equations, whose key ingredients are the electron's self-energy Σ, which is obtained by Hedin's GW approach, and the electron-hole interaction, which is described by the Bethe-Salpeter equation (BSE). The GW approach and the BSE predict orbital energies and excitation energies with high accuracy, respectively. We have calculated the low-lying excited states of a series of model biological chromophores, related to the photoactive yellow protein (PYP), rhodopsin, and the green fluorescent protein (GFP), obtaining a very good agreement with the available experimental and accurate theoretical data; the order of the excited states is also correctly predicted. MBGFT bridges the gap between time-dependent density functional theory and high-level quantum chemistry methods, combining the efficiency of the former with the accuracy of the latter: this makes MBGFT a promising method for studying excitations in complex biological systems.

Water-mediated spectral shifts in rhodopsin and bathorhodopsin

The role of water molecules in spectral tuning of proteins has been left largely unexplored. This topic is important because changing hydrogen bond patterns during the activation process may lead to spectral shifts which can be of diagnostic value for the underlying structures. Arguments put forward in this article are based on spectral shift calculations of the rhodopsin and bathorhodopsin chromophore due to wat2a and 2b in the presence and absence of the counterion and of the amino acids lining the rhodopsin binding pocket. They show, among others, that a single water molecule can shift the absorbance by up to 0.1 eV or 34 nm depending on the environment of the chromophore.

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.

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.

Probing the Sub-microsecond Photodissociation Dynamics in Gas-Phase Retinal Chromophores

The photoinduced fragmentation of a retinal model chromophore (all-trans-n-butyl protonated Schiff-base retinal) was studied in vacuo using a new experimental technique. The apparatus is able to record the photodissociation yield of gas-phase biomolecular ions in the first microseconds after absorption. Together with the existing ion storage ring ELISA, which operates on the millisecond to second time scale, the complete decay dynamics of such molecules can now be followed. In the case of retinal, the time-dependent fragmentation yield observed after irradiation with a 410 nm laser pulse exhibits contributions from one- and two-photon absorption, which decay non-exponentially with lifetimes on the order of 1 ms and 1 micros, respectively. The decay can be simulated using a statistical model, yielding good agreement with the experimental findings on both the millisecond and the microsecond time scales. No indication for nonstatistical processes is found for this molecule, the upper limit for a possible direct rate being a factor of 10(4) below the observed statistical dissociation rate.

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.

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.