The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution

Ultrafast photochemistry of anabaena sensory rhodopsin: experiment and theory

Light induced isomerization of the retinal chromophore activates biological function in all retinal protein (RP) driving processes such as ion-pumping, vertebrate vision and phototaxis in organisms as primitive as archea, or as complex as mammals. This process and its consecutive reactions have been the focus of experimental and theoretical research for decades. The aim of this review is to demonstrate how the experimental and theoretical research efforts can now be combined to reach a more comprehensive understanding of the excited state process on the molecular level. Using the Anabaena Sensory Rhodopsin as an example we will show how contemporary time-resolved spectroscopy and recently implemented excited state QM/MM methods consistently describe photochemistry in retinal proteins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Electrostatic spectral tuning mechanism of the green fluorescent protein

Understanding the mechanism of spectral tuning of biological chromophores is a major challenge in photobiology. We show here using large-scale full quantum chemical calculations of the green fluorescent protein that state-of-the-art coupled-cluster calculations provide accurate excitation energies and detailed insight about specific environmental effects. We obtain vertical excitation energies of 3.13 eV (396 nm) and 2.68 eV (463 nm), which are in quantitative agreement with the experimental absorption energies of 3.12 eV (397 nm) and 2.61 eV (475 nm) for the A- and B-forms of the protein. We find that the protein environment redshifts the absorption spectra by ∼0.56 eV and ∼0.22 eV for the two states, which can be attributed to ∼80% electrostatic effects and ∼20% steric effects.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

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.

Tuning the electronic absorption of protein-embedded all-trans-retinal

Protein-chromophore interactions are a central component of a wide variety of critical biological processes such as color vision and photosynthesis. To understand the fundamental elements that contribute to spectral tuning of a chromophore inside the protein cavity, we redesigned human cellular retinol binding protein II (hCRBPII) to fully encapsulate all-trans-retinal and form a covalent bond as a protonated Schiff base. This system, using rational mutagenesis designed to alter the electrostatic environment within the binding pocket of the host protein, enabled regulation of the absorption maximum of the pigment in the range of 425 to 644 nanometers. With only nine point mutations, the hCRBPII mutants induced a systematic shift in the absorption profile of all-trans-retinal of more than 200 nanometers across the visible spectrum.

Color vision: "OH-site" rule for seeing red and green

Eyes gather information, and color forms an extremely important component of the information, more so in the case of animals to forage and navigate within their immediate environment. By using the ONIOM (QM/MM) (ONIOM = our own N-layer integrated molecular orbital plus molecular mechanics) method, we report a comprehensive theoretical analysis of the structure and molecular mechanism of spectral tuning of monkey red- and green-sensitive visual pigments. We show that interaction of retinal with three hydroxyl-bearing amino acids near the β-ionone ring part of the retinal in opsin, A164S, F261Y, and A269T, increases the electron delocalization, decreases the bond length alternation, and leads to variation in the wavelength of maximal absorbance of the retinal in the red- and green-sensitive visual pigments. On the basis of the analysis, we propose the "OH-site" rule for seeing red and green. This rule is also shown to account for the spectral shifts obtained from hydroxyl-bearing amino acids near the Schiff base in different visual pigments: at site 292 (A292S, A292Y, and A292T) in bovine and at site 111 (Y111) in squid opsins. Therefore, the OH-site rule is shown to be site-specific and not pigment-specific and thus can be used for tracking spectral shifts in any visual pigment.

Anomalous pH Effect of Blue Proteorhodopsin

Proteorhodopsin (PR) is a light-driven proton pump found in marine bacteria, and thousands of PRs are classified into blue-absorbing PR (B-PR; λmax ≈ 490 nm) and green-absorbing PR (G-PR; λmax ≈ 525 nm). In this report, we present conversion of B-PR into G-PR using anomalous pH effect. B-PR in LC1-200, marine γ-proteobacteria, absorbs 497 and 513 nm maximally at pH 7 and 4, respectively, whose pH titration was reversible (pKa = 4.8). When pH was lowered from 4, the λmax was further red-shifted (528 nm at pH 2). This is unusual because blue shift occurs by chloride binding in the case of bacteriorhodopsin. Surprisingly, when pH was increased from 2 to 7, the λmax of this B-PR was further red-shifted to 540 nm, indicating that green-absorbing PR (PR540) is created only by changing pH. The present study reports the conformational flexibility of microbial rhodopsins, leading to the switch of absorbing color by a simple pH change.

Color Tuning in rhodopsins: the origin of the spectral shift between the chloride-bound and anion-free forms of halorhodopsin

Detailed knowledge of the molecular mechanisms that control the spectral properties in the rhodopsin protein family is important for understanding the functions of these photoreceptors and for the rational design of artificial photosensitive proteins. Here we used a high-level ab initio QM/MM method to investigate the mechanism of spectral tuning in the chloride-bound and anion-free forms of halorhodopsin from Natronobacterium pharaonis (phR) and the interprotein spectral shift between them. We demonstrate that the chloride ion tunes the spectral properties of phR via two distinct mechanisms: (i) electrostatic interaction with the chromophore, which results in a 95 nm difference between the absorption maxima of the two forms, and (ii) induction of a structural reorganization in the protein, which changes the positions of charged and polar residues and reduces this difference to 29 nm. The present study expands our knowledge concerning the role of the reorganization of the internal H-bond network for color tuning in general and provides a detailed investigation of the tuning mechanism in phR in particular.

Color tuning in photofunctional proteins

Depending on protein environment, a single photofunctional chromophore shows a wide variation of photoabsorption/emission energies. This photobiological phenomenon, known as color tuning, is observed in human visual cone pigments, firefly luciferase, and red fluorescent protein. We investigate the origin of color tuning by quantum chemical calculations on the excited states: symmetry-adapted cluster-configuration interaction (SAC-CI) method for excited states and a combined quantum mechanical (QM)/molecular mechanical (MM) method for protein environments. This Minireview summarizes our theoretical studies on the above three systems and explains a common feature of their color-tuning mechanisms. It also discuss the possibility of artificial color tuning toward a rational design of photoabsorption/emission properties.

Intramolecular charge transfer and dual fluorescence of 4-(dimethylamino)benzonitrile: ultrafast branching followed by a two-fold decay mechanism

In this contribution we present new experimental and theoretical results for the intramolecular charge transfer (ICT) reaction underlying the dual fluorescence of 4-(dimethylamino)benzonitrile (DMABN), which indicate that the fully twisted ICT (TICT) state is responsible for the time-resolved transient absorption spectrum while a distinct partially twisted ICT (pTICT) structure is suggested for the fluorescent ICT state.

Conical intersections in solution: formulation, algorithm, and implementation with combined quantum mechanics/molecular mechanics method

The significance of conical intersections in photophysics, photochemistry, and photodissociation of polyatomic molecules in gas phase has been demonstrated by numerous experimental and theoretical studies. Optimization of conical intersections of small- and medium-size molecules in gas phase has currently become a routine optimization process, as it has been implemented in many electronic structure packages. However, optimization of conical intersections of small- and medium-size molecules in solution or macromolecules remains inefficient, even poorly defined, due to large number of degrees of freedom and costly evaluations of gradient difference and nonadiabatic coupling vectors. In this work, based on the sequential quantum mechanics and molecular mechanics (QM/MM) and QM/MM-minimum free energy path methods, we have designed two conical intersection optimization methods for small- and medium-size molecules in solution or macromolecules. The first one is sequential QM conical intersection optimization and MM minimization for potential energy surfaces; the second one is sequential QM conical intersection optimization and MM sampling for potential of mean force surfaces, i.e., free energy surfaces. In such methods, the region where electronic structures change remarkably is placed into the QM subsystem, while the rest of the system is placed into the MM subsystem; thus, dimensionalities of gradient difference and nonadiabatic coupling vectors are decreased due to the relatively small QM subsystem. Furthermore, in comparison with the concurrent optimization scheme, sequential QM conical intersection optimization and MM minimization or sampling reduce the number of evaluations of gradient difference and nonadiabatic coupling vectors because these vectors need to be calculated only when the QM subsystem moves, independent of the MM minimization or sampling. Taken together, costly evaluations of gradient difference and nonadiabatic coupling vectors in solution or macromolecules can be reduced significantly. Test optimizations of conical intersections of cyclopropanone and acetaldehyde in aqueous solution have been carried out successfully.

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.

The protonation state of Glu181 in rhodopsin revisited: interpretation of experimental data on the basis of QM/MM calculations

The structure and spectroscopy of rhodopsin have been intensely studied in the past decade both experimentally and theoretically; however, important issues still remain unresolved. Of central interest is the protonation state of Glu181, where controversial and contradictory experimental evidence has appeared. While FTIR measurements indicate this residue to be unprotonated, preresonance Raman and UV-vis spectra have been interpreted in favor of a protonated Glu181. Previous computational approaches were not able to resolve this issue, providing contradicting data as well. Here, we perform hybrid QM/MM calculations using DFT methods for the electronic ground state, MRCI methods for the electronically excited states, and a polarization model for the MM part in order to investigate this issue systematically. We constructed various active-site models for protonated as well as unprotonated Glu181, which were evaluated by computing NMR, IR, Raman, and UV-vis spectroscopic data. The resulting differences in the UV-vis and Raman spectra between protonated and unprotonated models are very subtle, which has two major consequences. First, the common interpretation of prior Raman and UV-vis experiments in favor of a neutral Glu181 appears questionable, as it is based on the assumption that a charge at the Glu181 location would have a sizable impact. Second, also theoretical results should be interpreted with care. Spectroscopic differences between the structural models must be related to modeling uncertainties and intrinsic methodological errors. Despite a detailed comparison of various rhodopsins and mutants and consistently favorite results with charged Glu181 models, we find merely weak evidence from UV-vis and Raman calculations. On the contrary, difference FTIR and NMR chemical shift measurements on Rh mutants are indicative of the protonation state of Glu181. Supported by our results, they provide strong and independent evidence for a charged Glu181.

Anabaena sensory rhodopsin is a light-driven unidirectional rotor

The implementation of multiconfigurational quantum chemistry methods into a quantum-mechanics/molecular-mechanics protocol has allowed the construction of a realistic computer model for the sensory rhodopsin of the cyanobacterium Anabaena PCC 7120. The model, which reproduces the absorption spectra of both the all-trans and 13-cis forms of the protein and their associated K and L intermediates, is employed to investigate the light-driven steps of the photochromic cycle exhibited by the protein. It is found that the photoisomerizations of the all-trans and 13-cis retinal chromophores occur through unidirectional, counterclockwise 180° rotations of the =C14-C15= moiety with respect to the Lys210-linked end of the chromophore axis. Thus, the sequential interconversions of the all-trans and 13-cis forms during a single photochromic cycle yield a complete (360°) unidirectional rotation of the =C14-C15= moiety. This finding implies that Anabaena sensory rhodopsin is a biological realization of a light-driven molecular rotor.

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.

Color-tuning mechanism of firefly investigated by multi-configurational perturbation method

This is the first report on a multiconfigurational reference second-order perturbation theory-molecular mechanics study of the color modulation of the observed bioluminescence of the oxyluciferin-luciferase complex of the Japanese genji-botaru firefly using structures according to recent X-ray data. Our theoretical results do not support the experimentally deduced conclusion that the color modulation of the emitted light primarily depends on the size of the compact luciferase protein cavity embedding the excited oxyluciferin molecule. Rather, we find, in agreement with recent experimental observations, that the wavelength of the emitted light depends on the polarity of the microenvironment at the phenol/phenolate terminal of the benzothiazole fragment in oxyluciferin.

Excited states of fluorescent proteins, mKO and DsRed: chromophore-protein electrostatic interaction behind the color variations

The emitting states of green fluorescent protein (GFP), monomeric Kusabira orange (mKO), and Discosoma red (DsRed) were studied using QM/MM and SAC-CI methods. By comparing the electronic structures among the green-, orange-, and red-emitting states as well as their electrostatic and quantum mechanical interactions within the protein cavity, the basic mechanisms for determining emission colors have been clarified. We found that the orange and red emissions of mKO and DsRed, respectively, result from cancellation between two effects, the pi skeleton extension (red shift) and protein electrostatic potential (blue shift). The extension of the pi skeleton enhances the intramolecular charge-transfer character of the transition, which makes the fluorescence energy more sensitive to the protein's electrostatic potential. On the basis of this mechanism, we predicted amino acid mutations that could red shift the emission energy of DsRed. A novel single amino acid mutation, which was examined computationally, reduced the DsRed emission energy from 2.14 (579 nm) to 1.95 eV (636 nm), which is approaching near-infrared fluorescence.

MOLCAS 7: the next generation

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.

Color tuning in short wavelength-sensitive human and mouse visual pigments: ab initio quantum mechanics/molecular mechanics studies

We have investigated the protonation state and photoabsorption spectrum of Schiff-base (SB) nitrogen bound 11-cis-retinal in human blue and mouse UV cone visual pigments as well as in bovine rhodopsin by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have employed both multireference (MRCISD+Q, MR-SORCI+Q, and MR-DDCI2+Q) and single reference (TD-B3LYP and RI-CC2) QM methods. The calculated ground-state and vertical excitation energies show that UV-sensitive pigments have deprotonated SB nitrogen, while violet-sensitive pigments have protonated SB nitrogen, in agreement with some indirect experimental evidence. A significant blue shift of the absorption maxima of violet-sensitive pigments relative to rhodopsins arises from the increase in bond length alternation of the polyene chain of 11-cis-retinal induced by polarizing fields of these pigments. The main counterion is Glu113 in both violet-sensitive vertebrate pigments and bovine rhodopsin. Neither Glu113 nor the remaining pigment has a significant influence on the first excitation energy of 11-cis-retinal in the UV-sensitive pigments that have deprotonated SB nitrogen. There is no charge transfer between the SB and beta-ionone terminals of 11-cis-retinal in the ground and first excited states.

Molecular mechanism of long-range synergetic color tuning between multiple amino acid residues in conger rhodopsin

The synergetic effects of multiple rhodopsin mutations on color tuning need to be completely elucidated. Systematic genetic studies and spectroscopy have demonstrated an interesting example of synergetic color tuning between two amino acid residues in conger rhodopsin's ancestral pigment (p501): -a double mutation at one nearby and one distant residue led to a significant λ(max) blue shift of 13 nm, whereas neither of the single mutations at these two sites led to meaningful shifts.To analyze the molecular mechanisms of this synergetic color tuning, we performed homology modeling, molecular simulations, and electronic state calculations. For the double mutant, N195A/A292S, in silico mutation analysis demonstrated conspicuous structural changes in the retinal chromophore, whereas that of the single mutant, A292S, was almost unchanged. Using statistical ensembles of QM/MM optimized structures, the excitation energy of retinal chromophore was evaluated for the three visual pigments. As a result, the λ(max) shift of double mutant (DM) from p501 was -8 nm, while that of single mutant (SM) from p501 was +1 nm. Molecular dynamics simulation for DM demonstrated frequent isomerization between 6-s-cis and 6-s-trans conformers. Unexpectedly, however, the two conformers exhibited almost identical excitation energy, whereas principal component analysis (PCA) identified the retinal-counterion cooperative change of BLA (bond length alternation) and retinal-counterion interaction lead to the shift.

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.

Color-changing mutation in the E-F loop of proteorhodopsin

It is usually assumed that only amino acids located near the retinal chromophore are responsible for color tuning of rhodopsins. However, we recently found that replacement of Ala178 with Arg in the E-F loop of proteorhodopsin (PR), an archaeal-type rhodopsin in marine bacteria, shifts the lambda(max) from 525 to 545 nm at neutral pH [Yoshitsugu, M., Shibata, M., Ikeda, D., Furutani, Y., and Kandori, H. (2008) Angew. Chem., Int. Ed. 47, 3923-3926]. Since the location of Ala178 is distant from the retinal chromophore (approximately 25 A), the molecular mechanism of the unusual mutation effect on color tuning is intriguing. Here we studied this mechanism by using additional mutations and some analytical methods. Introduction of Arg into the corresponding amino acid in bacteriorhodopsin (BR, M163R mutant) does not change the absorption spectra, indicating that the effect is specific to PR. Introduction of Arg into the A-B or C-D loop yields little (3 nm) or no color change, respectively. T177R and P180R mutants exhibited absorption spectra identical to that of the wild type, while N176R and S179R mutants exhibit lambda(max) values of 528 and 535 nm, respectively. Therefore, the observed color change is position-specific, being fully effective at position 178 and half-effective at position 179. Salt affects the absorption spectra of wild-type and A178R PR similarly. FTIR spectroscopy at 77 K indicated similar chromophore structures for wild-type and A178R PR, and A178R PR pumps protons normally. We infer that the E-F loop has a unique structure in PR and the mutation of Ala178 disrupts the structure that includes the transmembrane region, leading to the observed changes in color and pK(a).

DNA nucleobase properties and photoreactivity: Modeling environmental effects

The accurate ab initio quantum chemical (QM) method multiconfigurational second-order perturbation (CASSPT2)/complete active space self-consistent field (CASSCF) has been used in conjunction with molecular mechanics (MM) procedures to compute molecular properties and photoinduced reactivity of DNA/RNA nucleobases (NABs) in isolation and within a realistic environment, in which the double helix strand, the aqueous media, and the external counterions are included. It is illustrated that the use of an MM model is helpful both to account for short- and long-range effects of the system surrounding the QM molecular core and to provide the proper structural constraints that allow more accurate QM geometry determinations.