Mechanism of Spectral Tuning Going from Retinal in Vacuo to Bovine Rhodopsin and its Mutants: Multireference ab Initio Quantum Mechanics/Molecular Mechanics Studies

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.

Intramolecular hydrogen bonding plays a crucial role in the photophysics and photochemistry of the GFP chromophore

In commonly studied GFP chromophore analogues such as 4-(4-hydroxybenzylidene)-1,2-dimethyl-1H-imidazol-5(4H)-one (PHBDI), the dominant photoinduced processes are cis-trans isomerization and subsequent S(1) → S(0) decay via a conical intersection characterized by a highly twisted double bond. The recently synthesized 2-hydroxy-substituted isomer (OHBDI) shows an entirely different photochemical behavior experimentally, since it mainly undergoes ultrafast intramolecular excited-state proton transfer, followed by S(1) → S(0) decay and ground-state reverse hydrogen transfer. We have chosen 4-(2-hydroxybenzylidene)-1H-imidazol-5(4H)-one (OHBI) to model the gas-phase photodynamics of such 2-hydroxy-substituted chromophores. We first use various electronic structure methods (DFT, TDDFT, CC2, DFT/MRCI, OM2/MRCI) to explore the S(0) and S(1) potential energy surfaces of OHBI and to locate the relevant minima, transition state, and minimum-energy conical intersection. These static calculations suggest the following decay mechanism: upon photoexcitation to the S(1) state, an ultrafast adiabatic charge-transfer induced excited-state intramolecular proton transfer (ESIPT) occurs that leads to the S(1) minimum-energy structure. Nearby, there is a S(1)/S(0) minimum-energy conical intersection that allows for an efficient nonadiabatic S(1) → S(0) internal conversion, which is followed by a fast ground-state reverse hydrogen transfer (GSHT). This mechanism is verified by semiempirical OM2/MRCI surface-hopping dynamics simulations, in which the successive ESIPT-GSTH processes are observed, but without cis-trans isomerization (which is a minor path experimentally with less than 5% yield). These gas-phase simulations of OHBI give an estimated first-order decay time of 476 fs for the S(1) state, which is larger but of the same order as the experimental values measured for OHBDI in solution: 270 fs in CH(3)CN and 230 fs in CH(2)Cl(2). The differences between the photoinduced processes of the 2- and 4-hydroxy-substituted chromophores are attributed to the presence or absence of intramolecular hydrogen bonding between the two rings.

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.

[Recent progress in the theoretical studies of structure, function, and reaction of biological molecules]

Essential biomolecular functions often involve electron-related events such as chemical reactions and photoluminescence phenomena. Theoretical description of such electronic processes requires the use of quantum mechanics (QM), but the number of atoms that can be handled with QM is usually smaller than the number of atoms present in a single protein. A reasonable strategy is therefore to give priority to a few tens or hundreds of atoms in the system and deal with them quantum mechanically. Lower-priority atoms influence the event occurring in the higher-priority area; therefore, their effect should also be taken into account. Under these circumstances, a reasonable approach is to apply two or more different theoretical methods to differently prioritized subsystems. QM can be combined, for example, with less accurate yet much less demanding molecular mechanics (MM). Our own N-layered integrated molecular orbital and molecular mechanics (ONIOM) method allows for such hybrid calculations, and our group has been applying it to a wide range of biology-related problems. In this paper, we briefly explain the theoretical background and the procedure for the theoretical investigation of biological systems. Subsequently, we provide an overview of some of our recent studies of metalloenzymes and photobiology-related problems.

Light activation of the isomerization and deprotonation of the protonated Schiff base retinal

We perform an ab initio analysis of the photoisomerization of the protonated Schiff base of retinal (PSB-retinal) from 11-cis to 11-trans rotating the C10-C11=C12-C13 dihedral angle from 0° (cis) to -180° (trans). We find that the retinal molecule shows the lowest rotational barrier (0.22 eV) when its charge state is zero as compared to the barrier for the protonated molecule which is ∼0.89 eV. We conclude that rotation most likely takes place in the excited state of the deprotonated retinal. The addition of a proton creates a much larger barrier implying a switching behavior of retinal that might be useful for several applications in molecular electronics. All conformations of the retinal compound absorb in the green region with small shifts following the dihedral angle rotation; however, the Schiff base of retinal (SB-retinal) at trans-conformation absorbs in the violet region. The rotation of the dihedral angle around the C11=C12 π-bond affects the absorption energy of the retinal and the binding energy of the SB-retinal with the proton at the N-Schiff; the binding energy is slightly lower at the trans-SB-retinal than at other conformations of the retinal.

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.

QM/MM study of dehydro and dihydro β-ionone retinal analogues in squid and bovine rhodopsins: implications for vision in salamander rhodopsin

Visual pigment rhodopsin provides a decisive crossing point for interaction between organisms and environment. Naturally occurring visual pigments contain only PSB11 and 3,4-dehydro-PSB11 as chromophores. Therefore, the ability of visual opsin to discriminate between the retinal geometries is investigated by means of QM/MM incorporation of PSB11, 6-s-cis and 6-s-trans forms of 3,4-dehydro-PSB11, and 3,4-dehydro-5,6-dihydro-PSB11 and 5,6-dihydro-PSB11 analogues into squid and bovine rhodopsin environments. The analogue-protein interaction reveals the binding site of squid rhodopsin to be malleable and ductile, while that of bovine rhodopsin is rigid and stiff. On the basis of these studies, a tentative model of the salamander rhodopsin binding site is also proposed.

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.

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.

Drawing the Retinal Out of Its Comfort Zone: An ONIOM(QM/MM) Study of Mutant Squid Rhodopsin

Engineering squid rhodopsin with modified retinal analogues is essential for understanding the conserved steric and electrostatic interaction networks that govern the architecture of the Schiff base binding site. Depriving the retinal of its steric and electrostatic contacts affects the positioning of the Schiff-base relative to the key residues Asn87, Tyr111, and Glu180. Displacement of the W1 and W2 positions and the impact on the structural rearrangements near the Schiff base binding region reiterates the need for the presence of internal water molecules and the accessibility of binding sites to them. Also, the dominant role of the Glu180 counterion in inducing the S(1)/S(2) state reversal for SBR is shown for the first time in squid rhodopsin.