Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching

Comparison of the isomerization mechanisms of human melanopsin and invertebrate and vertebrate rhodopsins

Comparative modeling and ab initio multiconfigurational quantum chemistry are combined to investigate the reactivity of the human nonvisual photoreceptor melanopsin. It is found that both the thermal and photochemical isomerization of the melanopsin 11-cis retinal chromophore occur via a space-saving mechanism involving the unidirectional, counterclockwise twisting of the =C11H-C12H= moiety with respect to its Lys340-linked frame as proposed by Warshel for visual pigments [Warshel A (1976) Nature 260(5553):679-683]. A comparison with the mechanisms documented for vertebrate (bovine) and invertebrate (squid) visual photoreceptors shows that such a mechanism is not affected by the diversity of the three chromophore cavities. Despite such invariance, trajectory computations indicate that although all receptors display less than 100 fs excited state dynamics, human melanopsin decays from the excited state ∼40 fs earlier than bovine rhodopsin. Some diversity is also found in the energy barriers controlling thermal isomerization. Human melanopsin features the highest computed barrier which appears to be ∼2.5 kcal mol(-1) higher than that of bovine rhodopsin. When assuming the validity of both the reaction speed/quantum yield correlation discussed by Warshel, Mathies and coworkers [Weiss RM, Warshel A (1979) J Am Chem Soc 101:6131-6133; Schoenlein RW, Peteanu LA, Mathies RA, Shank CV (1991) Science 254(5030):412-415] and of a relationship between thermal isomerization rate and thermal activation of the photocycle, melanopsin turns out to be a highly sensitive pigment consistent with the low number of melanopsin-containing cells found in the retina and with the extraretina location of melanopsin in nonmammalian vertebrates.

Modelling retinal chromophores photoisomerization: from minimal models in vacuo to ultimate bidimensional spectroscopy in rhodopsins

Modelling of retinal photoisomerization in different environments is reviewed and ultimate ultrafast electronic spectroscopy is proposed for obtaining new insights.

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.

Building a model of the blue cone pigment based on the wild type rhodopsin structure with QM/MM methods

Understanding the mechanism of color tuning of the retinal chromophore by its host protein became one of the key issues in the research on rhodopsins. While early mutation studies addressed its genetic origin, recent studies advanced to investigate its structural origin, based on X-ray crystallographic structures. For the human cone pigments, no crystal structures have been produced, and homology models were employed to elucidate the origin of its blue-shifted absorption. In this theoretical study, we take a different route to establish a structural model for human blue. Starting from the well-resolved structure of bovine rhodopsin, we derive multiple mutant models by stepwise mutation and equilibration using molecular dynamics simulations in a hybrid quantum mechanics/molecular mechanics framework. Our 30fold mutant reproduces the experimental UV-vis absorption shift of 0.45 eV and provides new insights about both structural and genetic factors that affect the excitation energy. Electrostatic effects of individual amino acids and collaborative structural effects are analyzed using semiempirical (OM2/MRCI) and ab initio (SORCI) multireference approaches.

Protein-bound water molecules in primate red- and green-sensitive visual pigments

Protein-bound water molecules play crucial roles in the structure and function of proteins. The functional role of water molecules has been discussed for rhodopsin, the light sensor for twilight vision, on the basis of X-ray crystallography, Fourier transform infrared (FTIR) spectroscopy, and a radiolytic labeling method, but nothing is known about the protein-bound waters in our color visual pigments. Here we apply low-temperature FTIR spectroscopy to monkey red (MR)- and green (MG)-sensitive color pigments at 77 K and successfully identify water vibrations using D(2)O and D(2)(18)O in the whole midinfrared region. The observed water vibrations are 6-8 for MR and MG, indicating that several water molecules are present near the retinal chromophore and change their hydrogen bonds upon retinal photoisomerization. In this sense, color visual pigments possess protein-bound water molecules essentially similar to those of rhodopsin. The absence of strongly hydrogen-bonded water molecules (O-D stretch at <2400 cm(-1)) is common between rhodopsin and color pigments, which greatly contrasts with the case of proton-pumping microbial rhodopsins. On the other hand, two important differences are observed in water signal between rhodopsin and color pigments. First, the water vibrations are identical between the 11-cis and 9-cis forms of rhodopsin, but different vibrational bands are observed at >2550 cm(-1) for both MR and MG. Second, strongly hydrogen-bonded water molecules (2303 cm(-1) for MR and 2308 cm(-1) for MG) are observed for the all-trans form after retinal photoisomerization, which is not the case for rhodopsin. These specific features of MR and MG can be explained by the presence of water molecules in the Cl(-)-biding site, which are located near positions C11 and C9 of the retinal chromophore. The averaged frequencies of the observed water O-D stretching vibrations for MR and MG are lower as the λ(max) is red-shifted, suggesting that water molecules are involved in the color tuning of our vision.

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.

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.

Color tuning in binding pocket models of the chlamydomonas-type channelrhodopsins

We examined the shift of absorption maxima between the chlamydomonas-type channelrhodopsins (ChRs) and bacteriorhodopsin (BR). Starting from the BR X-ray structure, we modeled the color tuning in the binding pockets of the ChRs by mutating up to 28 amino acids in the vicinity of the chromophore. By applying the efficient self-consistent charge density functional tight binding (SCC-DFTB) method in a quantum mechanical/molecular mechanical (QM/MM) framework, including explicit polarization and calculating excitation energies with the semiempirical OM2/MRCI method and the ab initio SORCI method, we have shown that multiple mutations in the binding pocket of BR causes large hypsochromic shifts that are of the same order as the experimentally observed shifts of the absorption maxima between BR and the ChRs. This study further demonstrates that mutations in the proximity of the Schiff base and complex counterion lead to a stronger but more flexible interaction with the retinal, which could serve as a possible explanation for the spectral patterns found in the ChRs.

Fundamental questions and concepts about photoreception and the case of Euglena gracilis

The ability to sense light can be considered the most fundamental and presumably the most ancient property of visual systems. This ability is the basis of phototaxis, one of the most striking behavioral responses of motile photosynthetic microorganisms (i.e. microalgae) to light stimuli, which allows them to move toward or away directional light. In order to fully exploit the information content of light (intensity, direction, distribution) microorganisms need proper perceiving devices, termed photoreceptors, which must act as sensors, to perceive wavelength and direction of light, as transducers, to convert the light signal into chemical and/or electrical information, but also as amplifiers and eventually as transmitters. This review describes the universal structural, behavioral and physiological features necessary for the proper functioning of these devices in algae, and how these features have been investigated by means of different analytical techniques such as for example microspectroscopy, digital fluorescence microscopy, two photons FLIM. The insight of the photoreceptive response mechanism is explained using the unicellular alga Euglena gracilis, in which the different structural, behavioral and physiological features combine to achieve a concerted, efficient response to light stimuli.

Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines

Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(‡) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.

Re-evaluation of rhodopsin's relaxation kinetics determined from femtosecond stimulated Raman lineshapes

This work presents a theoretical treatment of the vibrational line shape generated in a femtosecond stimulated Raman spectroscopy (FSRS) experiment under conditions in which the probed vibration undergoes a significant frequency shift during its free induction decay. This theory is applied to simulate the FSRS lineshapes previously observed in rhodopsin (Kukura et al. Science 2005, 310, 1006). The previously determined relaxation times for formation of the trans-photoproduct of rhodopsin were calculated using an incorrect equation for the time dependence of the observed frequency shifts. Here the data are reanalyzed by calculation of the corrected frequency sweep occurring during the vibrational free induction decay. It is shown that the calculated frequency shifts and general conclusions of the original work are sound but that the coherent vibrational frequency shifts of the C(10), C(11), and C(12) hydrogen-out-of-plane vibrations occur with a 140 fs time constant rather than the previously reported 325 fs time constant. This time constant provides an important constraint for models of the dynamics of the cis to trans isomerization process.

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.

Glutamic acid 181 is negatively charged in the bathorhodopsin photointermediate of visual rhodopsin

Assignment of the protonation state of the residue Glu-181 is important to our understanding of the primary event, activation processes and wavelength selection in rhodopsin. Despite extensive study, there is no general agreement on the protonation state of this residue in the literature. Electronic assignment is complicated by the location of Glu-181 near the nodal point in the electrostatic charge shift that accompanies excitation of the chromophore into the low-lying, strongly allowed ππ* state. Thus, the charge on this residue is effectively hidden from electronic spectroscopy. This situation is resolved in bathorhodopsin, because photoisomerization of the chromophore places Glu-181 well within the region of negative charge shift following excitation. We demonstrate that Glu-181 is negatively charged in bathorhodopsin on the basis of the shift in the batho absorption maxima at 10 K [λ(max) band (native) = 544 ± 2 nm, λ(max) band (E181Q) = 556 ± 3 nm] and the decrease in the λ(max) band oscillator strength (0.069 ± 0.004) of E181Q relative to that of the native protein. Because the primary event in rhodopsin does not include a proton translocation or disruption of the hydrogen-bonding network within the binding pocket, we may conclude that the Glu-181 residue in rhodopsin is also charged.

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.

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.

Modeling reaction routes from rhodopsin to bathorhodopsin

The quantum mechanical-molecular mechanical (QM/MM) theory was applied to calculate accurate structural parameters, vibrational and optical spectra of bathorhodopsin (BATHO), one of the primary photoproducts of the functional cycle of the visual pigment rhodopsin (RHO), and to characterize reaction routes from RHO to BATHO. The recently resolved crystal structure of BATHO (PDBID: 2G87) served as an initial source of coordinates of heavy atoms. Protein structures in the ground electronic state and vibrational frequencies were determined by using the density functional theory in the PBE0/cc-pVDZ approximation for the QM part and the AMBER force field parameters in the MM part. Calculated and assigned vibrational spectra of both model protein systems, BATHO and RHO, cover three main regions referring to the hydrogen-out-of-plan (HOOP) motion, the C==C ethylenic stretches, and the C--C single-bond stretches. The S(0)-S(1) electronic excitation energies of the QM part, including the chromophore group in the field of the protein matrix, were estimated by using the advanced quantum chemistry methods. The computed structural parameters as well as the spectral bands match perfectly the experimental findings. A structure of the transition state on the S(0) potential energy surface for the ground electronic state rearrangement from RHO to BATHO was located proving a possible route of the thermal protein activation to the primary photoproduct.

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.