Origin of Spectral Tuning in Rhodopsin—It Is Not the Binding Pocket

Low-temperature FTIR spectroscopy of the L/Q switch of proteorhodopsin

Rhodopsins are photoreceptive membrane proteins containing a retinal chromophore, and the color tuning mechanism in rhodopsins is one of the important topics. Color switch is a color-determining residue at the same position, where replacement of red- and blue-shifting amino acids in two wild-type rhodopsins causes spectral blue- and red-shifts, respectively. The first and most famous color switch in microbial rhodopsins is the L/Q switch in proteorhodopsins (PRs). Green- or blue-absorbing PR (GPR or BPR) contains Leu and Gln at position 105 of the C-helix (TM3), respectively, and their replacement converted absorbing colors. The L/Q switch enables bacteria to absorb green or blue light in shallow or deep ocean waters, respectively. Although Gln and Leu are hydrophilic and hydrophobic residues, respectively, a comprehensive mutation study of position 105 in GPR revealed that the λmax correlated with the volume of residues, not the hydropathy index. To gain structural insights into the mechanism, we applied low-temperature FTIR spectroscopy of L105Q GPR, and the obtained spectra were compared with those of GPR and BPR. The difference FTIR spectra of L105Q GPR were similar to those of BPR, not GPR, implying that the L/Q switch converts the GPR structure into a BPR structure in terms of the local environments of the retinal chromophore. It includes retinal skeletal vibration, hydrogen-bonding strength of the protonated Schiff base, amide-A vibration (peptide backbone), and protein-bound water molecules. Consequently color is switched accompanying such structural alterations, and known as the L/Q switch.

Molecular Mechanism of Spectral Tuning by Chloride Binding in Monkey Green Sensitive Visual Pigment

The visual pigments of the cones perceive red, green, and blue colors. The monkey green (MG) pigment possesses a unique Cl^- binding site; however, its relationship to the spectral tuning in green pigments remains elusive. Recently, FTIR spectroscopy revealed the characteristic structural modifications of the retinal binding site by Cl^- binding. Herein, we report the computational structural modeling of MG pigments and quantum-chemical simulation to investigate its spectral redshift and physicochemical relevance when Cl^- is present. Our protein structures reflect the previously suggested structural changes. AlphaFold2 failed to predict these structural changes. Excited-state calculations successfully reproduced the experimental red-shifted absorption energies, corroborating our protein structures. Electrostatic energy decomposition revealed that the redshift results from the His197 protonation state and conformations of Glu129, Ser202, and Ala308; however, Cl^- itself contributes to the blueshift. Site-directed mutagenesis supported our analysis. These modeled structures may provide a valuable foundation for studying cone pigments.

Protein Response Effects on Cofactor Excitation Energies from First Principles: Augmenting Subsystem Time-Dependent Density-Functional Theory with Many-Body Expansion Techniques

We investigate the possibility of describing protein response effects on a chromophore excitation by means of subsystem time-dependent density-functional theory (sTDDFT) in combination with a many-body expansion (MBE) approach. While sTDDFT is in principle intrinsically able to include such contributions, addressing cofactor excitations in protein models or entire proteins with full environment-response treatments is currently out of reach. Taking different model structures of the green fluorescent protein (GFP) and bovine rhodopsin as examples, we demonstrate that an embedded-MBE approach based on sTDDFT in its simplest version leads to a good agreement of the predicted protein response effect already at second order. To reproduce reference response effects from nonsubsystem TDDFT calculations quantitatively (error ≤ 5%), however, a third- or even fourth-order MBE may be required. For the latter case, we explore a selective inclusion of fourth-order terms that drastically reduces the computational burden. In addition, we demonstrate how this sTDDFT-MBE treatment can be utilized as an analysis tool to identify residues with dominant response contributions. This, in turn, can be employed to arrive at smaller structural models for light-absorbing proteins, which still feature all of the main characteristics in terms of photoresponse properties.

Electronic Couplings and Electrostatic Interactions Behind the Light Absorption of Retinal Proteins

The photo-functional chromophore retinal exhibits a wide variety of optical absorption properties depending on its intermolecular interactions with surrounding proteins and other chromophores. By utilizing these properties, microbial and animal rhodopsins express biological functions such as ion-transport and signal transduction. In this review, we present the molecular mechanisms underlying light absorption in rhodopsins, as revealed by quantum chemical calculations. Here, symmetry-adapted cluster-configuration interaction (SAC-CI), combined quantum mechanical and molecular mechanical (QM/MM), and transition-density-fragment interaction (TDFI) methods are used to describe the electronic structure of the retinal, the surrounding protein environment, and the electronic coupling between chromophores, respectively. These computational approaches provide successful reproductions of experimentally observed absorption and circular dichroism (CD) spectra, as well as insights into the mechanisms of unique optical properties in terms of chromophore-protein electrostatic interactions and chromophore-chromophore electronic couplings. On the basis of the molecular mechanisms revealed in these studies, we also discuss strategies for artificial design of the optical absorption properties of rhodopsins.

NeoR, a near-infrared absorbing rhodopsin

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λmax value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λmax values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λmax of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λmax using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λmax values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λmax of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λmax values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

The Effective Conjugation Length Is Responsible for the Red/Green Spectral Tuning in the Cyanobacteriochrome Slr1393g3

The origin of the spectral shift from a red- to a green-absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open-chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.

Side Methyl Groups Control the Conformation and Contribute to Symmetry Breaking of Isoprenoid Chromophores

Ab initio DFT computations reveal that the essential structural and photophysical features of the conjugated π-electron system of retinal and carotenoids are dictated by "innocent" methyl substituents. These methyl groups shape the conformation and symmetry of the isoprenoid chromophores by causing a sigmoidal distortion of the polyene skeleton and increasing its flexibility, which facilitates fitting to their binding pockets in proteins. Comparison of in vacuo conformations of the chromophores with their native (protein-bound) conformations showed, surprisingly, that the peripheral groups and interactions with the protein environment are much less significant than the methyl side groups in tuning their structural features. The methyl side groups also contribute to a loss of symmetry elements specific to linear polyenes. In effect, the symmetry-imposed restrictions on the chromophore electronic properties are disabled, which is of tremendous relevance to their photophysics. This is evidenced by their non-negligible permanent dipole moments and by the simulated and experimentally measured circular dichroism spectra, which necessarily reflect the chirality of the conjugated π-electron system.

Predicting peak spectral sensitivities of vertebrate cone visual pigments using atomistic molecular simulations

Vision is the dominant sensory modality in many organisms for foraging, predator avoidance, and social behaviors including mate selection. Vertebrate visual perception is initiated when light strikes rod and cone photoreceptors within the neural retina of the eye. Sensitivity to individual colors, i.e., peak spectral sensitivities (λmax) of visual pigments, are a function of the type of chromophore and the amino acid sequence of the associated opsin protein in the photoreceptors. Large differences in peak spectral sensitivities can result from minor differences in amino acid sequence of cone opsins. To determine how minor sequence differences could result in large spectral shifts we selected a spectrally-diverse group of 14 teleost Rh2 cone opsins for which sequences and λmax are experimentally known. Classical molecular dynamics simulations were carried out after embedding chromophore-associated homology structures within explicit bilayers and water. These simulations revealed structural features of visual pigments, particularly within the chromophore, that contributed to diverged spectral sensitivities. Statistical tests performed on all the observed structural parameters associated with the chromophore revealed that a two-term, first-order regression model was sufficient to accurately predict λmax over a range of 452-528 nm. The approach was accurate, efficient and simple in that site-by-site molecular modifications or complex quantum mechanics models were not required to predict λmax. These studies identify structural features associated with the chromophore that may explain diverged spectral sensitivities, and provide a platform for future, functionally predictive opsin modeling.

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.

Beyond spectral tuning: human cone visual pigments adopt different transient conformations for chromophore regeneration

Human red and green visual pigments are seven transmembrane receptors of cone photoreceptor cells of the retina that mediate color vision. These pigments share a very high degree of homology and have been assumed to feature analogous structural and functional properties. We report on a different regeneration mechanism among red and green cone opsins with retinal analogs using UV-Vis/fluorescence spectroscopic analyses, molecular modeling and site-directed mutagenesis. We find that photoactivated green cone opsin adopts a transient conformation which regenerates via an unprotonated Schiff base linkage with its natural chromophore, whereas red cone opsin forms a typical protonated Schiff base. The chromophore regeneration kinetics is consistent with a secondary retinal uptake by the cone pigments. Overall, our findings reveal, for the first time, structural differences in the photoactivated conformation between red and green cone pigments that may be linked to their molecular evolution, and support the proposal of secondary retinal binding to visual pigments, in addition to binding to the canonical primary site, which may serve as a regulatory mechanism of dark adaptation in the phototransduction process.

Rhodopsins carrying modified chromophores--the 'making of', structural modelling and their light-induced reactivity

A series of vitamin-A aldehydes (retinals) with modified alkyl group substituents (9-demethyl-, 9-ethyl-, 9-isopropyl-, 10-methyl, 10-methyl-13-demethyl-, and 13-demethyl retinal) was synthesized and their 11-cis isomers were used as chromophores to reconstitute the visual pigment rhodopsin. Structural changes were selectively introduced around the photoisomerizing C11=C12 bond. The effect of these structural changes on rhodopsin formation and bleaching was determined. Global fit of assembly kinetics yielded lifetimes and spectral features of the assembly intermediates. Rhodopsin formation proceeds stepwise with prolonged lifetimes especially for 9-demethyl retinal (longest lifetime τ3 = 7500 s, cf., 3500 s for retinal), and for 10-methyl retinal (τ3 = 7850 s). These slowed-down processes are interpreted as either a loss of fixation (9dm) or an increased steric hindrance (10me) during the conformational adjustment within the protein. Combined quantum mechanics and molecular mechanics (QM/MM) simulations provided structural insight into the retinal analogues-assembled, full-length rhodopsins. Extinction coefficients, quantum yields and kinetics of the bleaching process (μs-to-ms time range) were determined. Global fit analysis yielded lifetimes and spectral features of bleaching intermediates, revealing remarkably altered kinetics: whereas the slowest process of wild-type rhodopsin and of bleached and 11-cis retinal assembled rhodopsin takes place with lifetimes of 7 and 3.8 s, respectively, this process for 10-methyl-13-demethyl retinal was nearly 10 h (34670 s), coming to completion only after ca. 50 h. The structural changes in retinal derivatives clearly identify the precise interactions between chromophore and protein during the light-induced changes that yield the outstanding efficiency of rhodopsin.

Role of a helix B lysine residue in the photoactive site in channelrhodopsins

In most studied microbial rhodopsins two conserved carboxylic acid residues (the homologs of Asp-85 and Asp-212 in bacteriorhodopsin) and an arginine residue (the homolog of Arg-82) form a complex counterion to the protonated retinylidene Schiff base, and neutralization of the negatively charged carboxylates causes red shifts of the absorption maximum. In contrast, the corresponding neutralizing mutations in some relatively low-efficiency channelrhodopsins (ChRs) result in blue shifts. These ChRs do not contain a lysine residue in the second helix, conserved in higher efficiency ChRs (Lys-132 in the crystallized ChR chimera). By action spectroscopy of photoinduced channel currents in HEK293 cells and absorption spectroscopy of detergent-purified pigments, we found that in tested ChRs the Lys-132 homolog controls the direction of spectral shifts in the mutants of the photoactive site carboxylic acid residues. Analysis of double mutants shows that red spectral shifts occur when this Lys is present, whether naturally or by mutagenesis, and blue shifts occur when it is replaced with a neutral residue. A neutralizing mutation of the Lys-132 homolog alone caused a red spectral shift in high-efficiency ChRs, whereas its introduction into low-efficiency ChR1 from Chlamydomonas augustae (CaChR1) caused a blue shift. Taking into account that the effective charge of the carboxylic acid residues is a key factor in microbial rhodopsin spectral tuning, these findings suggest that the Lys-132 homolog modulates their pKa values. On the other hand, mutation of the Arg-82 homolog that fulfills this role in bacteriorhodopsin caused minimal spectral changes in the tested ChRs. Titration revealed that the pKa of the Asp-85 homolog in CaChR1 lies in the alkaline region unlike in most studied microbial rhodopsins, but is substantially decreased by introduction of a Lys-132 homolog or neutralizing mutation of the Asp-212 homolog. In the three ChRs tested the Lys-132 homolog also alters channel current kinetics.

The EF loop in green proteorhodopsin affects conformation and photocycle dynamics

The proteorhodopsin family consists of retinal proteins of marine bacterial origin with optical properties adjusted to their local environments. For green proteorhodopsin, a highly specific mutation in the EF loop, A178R, has been found to cause a surprisingly large redshift of 20 nm despite its distance from the chromophore. Here, we analyze structural and functional consequences of this EF loop mutation by time-resolved optical spectroscopy and solid-state NMR. We found that the primary photoreaction and the formation of the K-like photo intermediate is almost pH-independent and slower compared to the wild-type, whereas the decay of the K-intermediate is accelerated, suggesting structural changes within the counterion complex upon mutation. The photocycle is significantly elongated mainly due to an enlarged lifetime of late photo intermediates. Multidimensional MAS-NMR reveals mutation-induced chemical shift changes propagating from the EF loop to the chromophore binding pocket, whereas dynamic nuclear polarization-enhanced (13)C-double quantum MAS-NMR has been used to probe directly the retinylidene conformation. Our data show a modified interaction network between chromophore, Schiff base, and counterion complex explaining the altered optical and kinetic properties. In particular, the mutation-induced distorted structure in the EF loop weakens interactions, which help reorienting helix F during the reprotonation step explaining the slower photocycle. These data lead to the conclusion that the EF loop plays an important role in proton uptake from the cytoplasm but our data also reveal a clear interaction pathway between the EF loop and retinal binding pocket, which might be an evolutionary conserved communication pathway in retinal proteins.

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.

The photochemical determinants of color vision: revealing how opsins tune their chromophore's absorption wavelength

The evolution of a variety of important chromophore-dependent biological processes, including microbial light sensing and mammalian color vision, relies on protein modifications that alter the spectral characteristics of a bound chromophore. Three different color opsins share the same chromophore, but have three distinct absorptions that together cover the entire visible spectrum, giving rise to trichromatic vision. The influence of opsins on the absorbance of the chromophore has been studied through methods such as model compounds, opsin mutagenesis, and computational modeling. The recent development of rhodopsin mimic that uses small soluble proteins to recapitulate the binding and wavelength tuning of the native opsins provides a new platform for studying protein-regulated spectral tuning. The ability to achieve far-red shifted absorption in the rhodopsin mimic system was attributed to a combination of the lack of a counteranion proximal to the iminium, and a uniformly neutral electrostatic environment surrounding the chromophore.

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.

Nonadiabatic photodynamics of a retinal model in polar and nonpolar environment

The nonadiabatic photodynamics of the all-trans-2,4-pentadiene-iminium cation (protonated Schiff base 3, PSB3) and the all-trans-3-methyl-2,4-pentadiene-iminium cation (MePSB3) were investigated in the gas phase and in polar (aqueous) and nonpolar (n-hexane) solutions by means of surface hopping using a multireference configuration-interaction (MRCI) quantum mechanical/molecular mechanics (QM/MM) level. Spectra, lifetimes for radiationless deactivation to the ground state, and structural and electronic parameters are compared. A strong influence of the polar solvent on the location of the crossing seam, in particular in the bond length alternation (BLA) coordinate, is found. Additionally, inclusion of the polar solvent changes the orientation of the intersection cone from sloped in the gas phase to peaked, thus enhancing considerably its efficiency for deactivation of the molecular system to the ground state. These factors cause, especially for MePSB3, a substantial decrease in the lifetime of the excited state despite the steric inhibition by the solvent.

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.

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.

Why 11-cis-retinal? Why not 7-cis-, 9-cis-, or 13-cis-retinal in the eye?

One of the basic and unresolved puzzles in the chemistry of vision concerns the natural selection of 11-cis-retinal as the light-sensing chromophore in visual pigments. A detailed computational examination of the structure, stability, energetics, and spectroscopy of 7-cis-, 9-cis-, 11-cis-, and 13-cis-retinal isomers in vertebrate (bovine, monkey) and invertebrate (squid) visual pigments was carried out using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The results show that the electrostatic interaction between retinal and opsin dominates the natural selection of 11-cis-retinal over other cis isomers in the dark state. In all of the pigments, 9-cis-retinal was found to be only slightly higher in energy than 11-cis-retinal, which provides strong evidence for the presence of 9-cis-rhodopsin in nature. 7-cis-Retinal is suggested to be an "upside-down" version of the all-trans isomer because the structural rearrangements observed for 7-cis-rhodopsin from squid were found to be very similar to those for squid bathorhodopsin. The progressive red shift in the calculated absorption wavelength (λ(max)) (431, 456, 490, and 508 nm for the 7-cis-, 9-cis-, 11-cis-, and 13-cis-retinal isomers) is due to the decrease in bond length alternation of the retinal.

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.

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.

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.

Model systems for the investigation of the opsin shift in bacteriorhodopsin

Donor-acceptor substituted styrenes and phenylbutadienes with substituents varying in donor and acceptor strength and as reconstituted chromophore-protein complexes were investigated as model compounds for the protonated Schiff base chromophore in bacteriorhodopsin (bR) both experimentally and theoretically. Charge distribution, donor-acceptor strength, and the shift of the absorption energy are correlated. The effect of the external electrostatic field was tested with a compound carrying an additional nonconjugated charge. The concept of overpolarization by the external charge, that is, the reversal of the relative importance of the two main resonance structures in S(0) and S(1), has been emphasized and related to a simple qualitative 2 x 2 interaction model. The variable donor approach is a new way for a better understanding of the Opsin shift in Bacteriorhodopsin.