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

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.

Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment

Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.

Calcium-Sensitive Microbial Rhodopsin VirChR1: A Femtosecond to Second Photocycle Study

Viral rhodopsins are light-gated cation channels representing a novel class of microbial rhodopsins. For viral rhodopsin 1 subfamily members VirChR1 and OLPVR1, channel activity is abolished above a certain calcium concentration. Here we present a calcium-dependent spectroscopic analysis of VirChR1 on the femtosecond to second time scale. Unlike channelrhodopsin-2, VirChR1 possesses two intermediate states P1 and P2 on the ultrafast time scale, similar to J and K in ion-pumping rhodopsins. Subsequently, we observe multifaceted photocycle kinetics with up to seven intermediate states. Calcium predominantly affects the last photocycle steps, including the appearance of additional intermediates P6Ca and P7 representing the blocked channel. Furthermore, the photocycle of the counterion variant D80N is drastically altered, yielding intermediates with different spectra and kinetics compared to those of the wt. These findings demonstrate the central role of the counterion within the defined reaction sequence of microbial rhodopsins that ultimately defines the protein function.

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.

Excited-state dynamics of all-trans protonated retinal Schiff base in CRABPII-based rhodopsin mimics

The rhodopsin mimic is a chemically synthetized complex with retinyl Schiff base (RSB) formed between protein and the retinal chromophore that can mimic the natural rhodopsin-like protein. The artificial rhodopsin mimic is more stable and designable than the natural protein and hence has wider uses in photon detection devices. The mimic structure RSB, like the case in the actual rhodopsin-like protein, undergoes isomerization and protonation throughout the photoreaction process. As a result, understanding the dynamics of the RSB in the photoreaction process is critical. In this study, the ultrafast transient absorption spectra of three mutants of the cellular retinoic acid-binding protein II-based rhodopsin mimic at acidic environment were recorded, from which the related excited-state dynamics of the all-trans protonated RSB (AT-PRSB) were investigated. The transient fluorescence spectra measurements are used to validate some of the dynamic features. We find that the excited-state dynamics of AT-PRSB in three mutants share a similar pattern that differs significantly from the dynamics of 15-cis PRSB of the rhodopsin mimic in neutral solution. By comparing the dynamics across the three mutants, we discovered that the aromatic residues near the β-ionone ring structure of the retinal may help stabilize the AT-PRSB and hence slow down its isomerization rate. The experimental results provide implications on designing a rhodopsin-like protein with significant infrared fluorescence, which can be particularly useful in the applications in biosensing or bioimaging in deeper tissues.

Mechanism of the Irreversible Transition from Pentamer to Monomer at pH 2 in a Blue Proteorhodopsin

Proteorhodopsin (PR) is a light-driven proton pump found in marine bacteria, and thousands of PRs are classified as blue-absorbing PRs (BPR; λ_max ∼ 490 nm) and green-absorbing PRs (GPR; λ_max ∼ 525 nm). We previously converted BPR into GPR using an anomalous pH effect, which was achieved by an irreversible process at around pH 2. Recent size-exclusion chromatography (SEC) and atomic force microscopy (AFM) analyses of BPR from Vibrio califitulae (VcBPR) revealed the anomalous pH effect owing to the irreversible transition from pentamer to monomer. Different pK_a values of the Schiff base counterion between pentamer and monomer lead to different colors at the same pH. Here, we incorporate systematic mutation into VcBPR and examine the anomalous pH effect. The anomalous pH effect was observed for the mutants of key residues near the retinal chromophore such as D76N, D206N, and Q84L, indicating that the Schiff base counterions and the L/Q switch do not affect the irreversible transition from pentamer to monomer at pH ∼ 2. We then focus on the two specific interactions at the intermonomer interface in a pentamer, E29/R30/D31 and W13/H54. Single mutants such as E29Q, R30A, W13A, and H54A and the wild type (WT) exhibited an anomalous pH effect. In contrast, the anomalous pH effect was lost for E29Q/H54A, R30A/H54A, and W13A/E29Q. Size-exclusion chromatography (SEC) and atomic force microscopy (AFM) measurements showed monomer forms in the original states of the double mutants, being a clear contrast to the pentamer forms of all single mutants in the original states. It was concluded that the pentamer structure of VcBPR was stabilized by an electrostatic interaction in the E29/R30/D31 region and a hydrogen-bonding interaction in the W13/H54 region, which was disrupted at pH 2 and converted into monomers.

Additive and epistatic effects influence spectral tuning in molluscan retinochrome opsin

The relationship between genotype and phenotype is nontrivial due to often complex molecular pathways that make it difficult to unambiguously relate phenotypes to specific genotypes. Photopigments, an opsin apoprotein bound to a light-absorbing chromophore, present an opportunity to directly relate the amino acid sequence to an absorbance peak phenotype (λmax). We examined this relationship by conducting a series of site-directed mutagenesis experiments of retinochrome, a non-visual opsin, from two closely related species: the common bay scallop, Argopecten irradians, and the king scallop, Pecten maximus. Using protein folding models, we identified three amino acid sites of likely functional importance and expressed mutated retinochrome proteins in vitro. Our results show that the mutation of amino acids lining the opsin binding pocket are responsible for fine spectral tuning, or small changes in the λmax of these light sensitive proteins Mutations resulted in a blue or red shift as predicted, but with dissimilar magnitudes. Shifts ranged from a 16 nm blue shift to a 12 nm red shift from the wild-type λmax. These mutations do not show an additive effect, but rather suggests the presence of epistatic interactions. This work highlights the importance of binding pocket shape in the evolution of spectral tuning and builds on our ability to relate genotypic changes to phenotypes in an emerging model for opsin functional analysis.

Molecular Origin of the Anomalous pH Effect in Blue Proteorhodopsin

Proteorhodopsin (PR) is a light-driven proton pump found in marine bacteria, and thousands of PRs are classified into blue-absorbing PR (BPR; λ_max ∼ 490 nm) and green-absorbing PR (GPR; λ_max ∼ 525 nm). We previously presented conversion of BPR into GPR using the anomalous pH effect. When we lowered the pH of a BPR to pH 2 and returned to pH 7, the protein absorbs green light. This suggests the existence of the critical point of the irreversible process at around pH 2, but the mechanism of anomalous pH effect was fully unknown. The present size exclusion chromatography (SEC) and atomic force microscope (AFM) analysis of BPR from Vibrio califitulae (VcBPR) revealed the anomalous pH effect because of the conversion from pentamer to monomer. The different pK_a of the Schiff base counterion between pentamer and monomer leads to different colors at the same pH.

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.

Simple Models to Study Spectral Properties of Microbial and Animal Rhodopsins: Evaluation of the Electrostatic Effect of Charged and Polar Residues on the First Absorption Band Maxima

A typical feature of proteins from the rhodopsin family is the sensitivity of their absorption band maximum to protein amino acid composition. For this reason, studies of these proteins often require methodologies that determine spectral shift caused by amino acid substitutions. Generally, quantum mechanics/molecular mechanics models allow for the calculation of a substitution-induced spectral shift with high accuracy, but their application is not always easy and requires special knowledge. In the present study, we propose simple models that allow us to estimate the direct effect of a charged or polar residue substitution without extensive calculations using only rhodopsin three-dimensional structure and plots or tables that are provided in this article. The models are based on absorption maximum values calculated at the SORCI+Q level of theory for cis- and trans-forms of retinal protonated Schiff base in an external electrostatic field of charges and dipoles. Each value corresponds to a certain position of a charged or polar residue relative to the retinal chromophore. The proposed approach was evaluated against an example set consisting of twelve bovine rhodopsin and sodium pumping rhodopsin mutants. The limits of the applicability of the models are also discussed. The results of our study can be useful for the interpretation of experimental data and for the rational design of rhodopsins with required spectral properties.

Color-tuning of natural variants of heliorhodopsin

Microbial rhodopsins are distributed through many microorganisms. Heliorhodopsins are newly discovered but have an unclear function. They have seven transmembrane helices similar to type-I and type-II rhodopsins, but they are different in that the N-terminal region of heliorhodopsin is cytoplasmic. We chose 13 representative heliorhodopsins from various microorganisms, expressed and purified with an N-terminal His tag, and measured the absorption spectra. The 13 natural variants had an absorption maximum (λmax) in the range 530–556 nm similar to proteorhodopsin (λmax = 490–525 nm). We selected several candidate residues that influence rhodopsin color-tuning based on sequence alignment and constructed mutants via site-directed mutagenesis to confirm the spectral changes. We found two important residues located near retinal chromophore that influence λmax. We also predict the 3D structure via homology-modeling of Thermoplasmatales heliorhodopsin. The results indicate that the color-tuning mechanism of type-I rhodopsin can be applied to understand the color-tuning of heliorhodopsin.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Benchmarking the Performance of Time-Dependent Density Functional Theory Methods on Biochromophores

Quantum chemical calculations are important for elucidating light-capturing mechanisms in photobiological systems. The time-dependent density functional theory (TDDFT) has become a popular methodology because of its balance between accuracy and computational scaling, despite its problems in describing, for example, charge transfer states. As a step toward systematically understanding the performance of TDDFT calculations on biomolecular systems, we study here 17 commonly used density functionals, including seven long-range separated functionals, and compare the obtained results with excitation energies calculated at the approximate second order coupled-cluster theory level (CC2). The benchmarking set includes the first five singlet excited states of 11 chemical analogues of biochromophores from the green fluorescent protein, rhodopsin/bacteriorhodopsin (Rh/bR), and the photoactive yellow protein. We find that commonly used pure density functionals such as BP86, PBE, M11-L, and hybrid functionals with 20-25% of Hartree-Fock (HF) exchange (B3LYP, PBE0) have a tendency to consistently underestimate vertical excitation energies (VEEs) relative to the CC2 values, whereas hybrid density functionals with around 50% HF exchange such as BHLYP, PBE50, and M06-2X and long-range corrected functionals such as CAM-B3LYP, ωPBE, ωPBEh, ωB97X, ωB97XD, BNL, and M11 overestimate the VEEs. We observe that calculations using the CAM-B3LYP and ωPBEh functionals with 65% and 100% long-range HF exchange, respectively, lead to an overestimation of the VEEs by 0.2-0.3 eV for the benchmarking set. To reduce the systematic error, we introduce here two new empirical functionals, CAMh-B3LYP and ωhPBE0, for which we adjusted the long-range HF exchange to 50%. The introduced parameterization reduces the mean signed average (MSA) deviation to 0.07 eV and the root mean square (rms) deviation to 0.17 eV as compared to the CC2 values. In the present study, TDDFT calculations using the aug-def2-TZVP basis sets, the best performing functionals relative to CC2 are ωhPBE0 (rms = 0.17, MSA = 0.06 eV); CAMh-B3LYP (rms = 0.16, MSA = 0.07 eV); and PBE0 (rms = 0.23, MSA = -0.14 eV). For the popular range-separated CAM-B3LYP functional, we obtain an rms value of 0.31 eV and an MSA value of 0.25 eV, which can be compared with the rms and MSA values of 0.37 and -0.31 eV, respectively, as obtained at the B3LYP level.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Red-shifting mutation of light-driven sodium-pump rhodopsin

Microbial rhodopsins are photoreceptive membrane proteins that transport various ions using light energy. While they are widely used in optogenetics to optically control neuronal activity, rhodopsins that function with longer-wavelength light are highly demanded because of their low phototoxicity and high tissue penetration. Here, we achieve a 40-nm red-shift in the absorption wavelength of a sodium-pump rhodopsin (KR2) by altering dipole moment of residues around the retinal chromophore (KR2 P219T/S254A) without impairing its ion-transport activity. Structural differences in the chromophore of the red-shifted protein from that of the wildtype are observed by Fourier transform infrared spectroscopy. QM/MM models generated with an automated protocol show that the changes in the electrostatic interaction between protein and chromophore induced by the amino-acid replacements, lowered the energy gap between the ground and the first electronically excited state. Based on these insights, a natural sodium pump with red-shifted absorption is identified from Jannaschia seosinensis.

Microbial rhodopsins are photoreceptive and widely used in optogenetics for which they should preferable function with longer-wavelength light. Here, authors achieve a 40-nm red-shift in the absorption wavelength of a sodium-pump rhodopsin (KR2) by altering the distribution of the retinal chromophore.

Expression of Anabaena sensory rhodopsin is influenced by different codons of seven residues at the N-terminal region

Microbial rhodopsins are well-known seven-transmembrane proteins that have been extensively studied for their structure and function. These retinal-binding proteins can be divided into two types. Type I is microbial rhodopsin, and type II (visual pigment) is expressed mostly in mammalian eyes. The two primary functions of type I rhodopsin are ion pumping activity and sensory transduction. Anabaena sensory rhodopsin (ASR) is a microbial rhodopsin with a specific function of photosensory transduction. ASR is expressed at moderate levels in Escherichia coli, but its expression level is lower compared to the general green light absorbing proteorhodopsin (GPR). In this study, full-length ASR was used to test the influence of codon usage on expression E. coli. Seven amino acids at the N-terminal region of ASR after the Met start codon were changed randomly using designed primers, which allowed for 8192 different nucleotide combinations. The codon changes were screened for the preferable codons that resulted in higher expression yield. Among the 57 selected mutations, 24 color-enhanced E. coli colonies contained ASR proteins, and they expressed ASR at a higher level than the bacteria with wild-type ASR codon usage. This result strongly suggests that the specific codon usage of only the N-terminal portion of a protein can increase the expression level of the entire protein.

Impact of Electronic State Mixing on the Photoisomerization Time Scale of the Retinal Chromophore

Spectral data show that the photoisomerization of retinal protonated Schiff base (rPSB) chromophores occurs on a 100 fs time scale or less in vertebrate rhodopsins, it is several times slower in microbial rhodopsins and it is between one and 2 orders of magnitude slower in solution. These time scale variations have been attributed to specific modifications of the topography of the first excited state potential energy surface of the chromophore. However, it is presently not clear which specific environment effects (e.g., electrostatic, electronic, or steric) are responsible for changing the surface topography. Here, we use QM/MM models and excited state trajectory computations to provide evidence for an increase in electronic mixing between the first and the second excited state of the chromophore when going from vertebrate rhodopsin to the solution environments. Ultimately, we argue that a correlation between the lifetime of the first excited state and electronic mixing between such state and its higher neighbor, may have been exploited to evolve rhodopsins toward faster isomerization and, possibly, light-sensitivity.

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

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

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

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

pH-Dependent absorption spectrum of a protein: a minimal electrostatic model of Anabaena sensory rhodopsin

A minimal electrostatic model is introduced which aims at reproducing and analyzing the visible-light absorption energy shift of a protein with pH.

Electrostatic versus Resonance Interactions in Photoreceptor Proteins: The Case of Rhodopsin

Light sensing in photoreceptor proteins is subtly modulated by the multiple interactions between the chromophoric unit and its binding pocket. Many theoretical and experimental studies have tried to uncover the fundamental origin of these interactions but reached contradictory conclusions as to whether electrostatics, polarization, or intrinsically quantum effects prevail. Here, we select rhodopsin as a prototypical photoreceptor system to reveal the molecular mechanism underlying these interactions and regulating the spectral tuning. Combining a multireference perturbation method and density functional theory with a classical but atomistic and polarizable embedding scheme, we show that accounting for electrostatics only leads to a qualitatively wrong picture, while a responsive environment can successfully capture both the classical and quantum dominant effects. Several residues are found to tune the excitation by both differentially stabilizing ground and excited states and through nonclassical "inductive resonance" interactions. The results obtained with such a quantum-in-classical model are validated against both experimental data and fully quantum calculations.

Polarizable Embedding Density Matrix Renormalization Group

The polarizable embedding (PE) approach is a flexible embedding model where a preselected region out of a larger system is described quantum mechanically, while the interaction with the surrounding environment is modeled through an effective operator. This effective operator represents the environment by atom-centered multipoles and polarizabilities derived from quantum mechanical calculations on (fragments of) the environment. Thereby, the polarization of the environment is explicitly accounted for. Here, we present the coupling of the PE approach with the density matrix renormalization group (DMRG). This PE-DMRG method is particularly suitable for embedded subsystems that feature a dense manifold of frontier orbitals which requires large active spaces. Recovering such static electron-correlation effects in multiconfigurational electronic structure problems, while accounting for both electrostatics and polarization of a surrounding environment, allows us to describe strongly correlated electronic structures in complex molecular environments. We investigate various embedding potentials for the well-studied first excited state of water with active spaces that correspond to a full configuration-interaction treatment. Moreover, we study the environment effect on the first excited state of a retinylidene Schiff base within a channelrhodopsin protein. For this system, we also investigate the effect of dynamical correlation included through short-range density functional theory.

Space and Time Evolution of the Electrostatic Potential During the Activation of a Visual Pigment

Animal and microbial retinal proteins employ the Schiff base of retinal as their chromophore. Here, the possible consequences of the charge translocation associated with the light-induced dynamics of the chromophore of a visual opsin are investigated along a representative semiclassical trajectory. We show that the evolution of the electrostatic potential projected by the chromophore onto the surrounding protein displays intense but topographically localized sudden variations in proximity of the decay region. pKa calculations carried out on selected snapshots used as probes, indicate that the only residue which may be sensitive to the electrostatic potential shift is Glu181. Accordingly, our results suggest that the frail Tyr191/268-Glu181-Wat2-Ser186 hydrogen bond network may be perturbed by the transient variations of the electrostatic potential.

Theoretical Investigation of the Geometries and UV-vis Spectra of Poly(l-glutamic acid) Featuring a Photochromic Azobenzene Side Chain

The geometries and UV-vis spectra of azobenzene dyes grafted as a side chain on poly(l-glutamic acid) have been investigated using a combination of quantum mechanics/molecular mechanics (QM/MM) and time-dependent density functional theory (TD-DFT) methods at the TD-PBE0/6-311+G(d,p)//B3LYP/6-311G(d,p):Amber ff99 level of theory. The influence of the secondary structure of the polypeptide on the electronic properties of both the trans and cis conformations of azobenzene dyes has been studied. It turns out that the grafted dyes exhibit a red-shift of the π → π* absorption energies mainly due to the auxochromic shift induced by the peptidic group used to link the chromophoric unit to the polypeptide and that specific interactions between the glutamic side chain and the azobenzene moiety lead to a large blue-shift of the n → π* transition.

Protein Influence on Electronic Spectra Modeled by Multipoles and Polarizabilities

We have developed automatic methods to calculate multipoles and anisotropic polarizabilities for all atoms and bond centers in a protein and to include such a model in the calculation of electronic properties at any level of quantum mechanical theory. This approach is applied for the calculation of the electronic spectra of retinal in rhodopsin at the CASPT2//CASSCF level (second-order multiconfigurational perturbation theory) for the wild-type protein, as well as two mutants and isorhodopsin in QM/MM structures based on two crystal structures. We also perform a detailed investigation of the importance and distance dependence of the multipoles and the polarizabilities for both the absolute and the relative absorption energies. It is shown that the model of the surrounding protein strongly influences the spectrum and that different models give widely different results. For example, the Amber 1994 and 2003 force fields give excitation energies that differ by up to 16 kJ/mol. For accurate excitation energies, multipoles up to quadrupoles and anisotropic polarizabilities are needed. However, interactions with residues more than 10 Å from the chromophore can be treated with a standard polarizable force field without any dipoles or quadrupoles.

Molecular bases for the selection of the chromophore of animal rhodopsins

The functions of microbial and animal rhodopsins are triggered by the isomerization of their all-trans and 11-cis retinal chromophores, respectively. To lay the molecular basis driving the evolutionary transition from the all-trans to the 11-cis chromophore, multiconfigurational quantum chemistry is used to compare the isomerization mechanisms of the sensory rhodopsin from the cyanobacterium Anabaena PCC 7120 (ASR) and of the bovine rhodopsin (Rh). It is found that, despite their evolutionary distance, these eubacterial and vertebrate rhodopsins start to isomerize via distinct implementations of the same bicycle-pedal mechanism originally proposed by Warshel [Warshel A (1976) Nature 260:678-683]. However, by following the electronic structure changes of ASR (featuring the all-trans chromophore) during the isomerization, we find that ASR enters a region of degeneracy between the first and second excited states not found in Rh (featuring the 11-cis chromophore). We show that such degeneracy is modulated by the preorganized structure of the chromophore and by the position of the reactive double bond. It is argued that the optimization of the electronic properties of the chromophore, which affects the photoisomerization efficiency and the thermal isomerization barrier, provided a key factor for the emergence of the striking amino acid sequence divergence observed between the microbial and animal rhodopsins.

Advances in understanding the molecular basis of the first steps in color vision

Serving as one of our primary environmental inputs, vision is the most sophisticated sensory system in humans. Here, we present recent findings derived from energetics, genetics and physiology that provide a more advanced understanding of color perception in mammals. Energetics of cis-trans isomerization of 11-cis-retinal accounts for color perception in the narrow region of the electromagnetic spectrum and how human eyes can absorb light in the near infrared (IR) range. Structural homology models of visual pigments reveal complex interactions of the protein moieties with the light sensitive chromophore 11-cis-retinal and that certain color blinding mutations impair secondary structural elements of these G protein-coupled receptors (GPCRs). Finally, we identify unsolved critical aspects of color tuning that require future investigation.

Atomistic design of microbial opsin-based blue-shifted optogenetics tools

Microbial opsins with a bound chromophore function as photosensitive ion transporters and have been employed in optogenetics for the optical control of neuronal activity. Molecular engineering has been utilized to create colour variants for the functional augmentation of optogenetics tools, but was limited by the complexity of the protein-chromophore interactions. Here we report the development of blue-shifted colour variants by rational design at atomic resolution, achieved through accurate hybrid molecular simulations, electrophysiology and X-ray crystallography. The molecular simulation models and the crystal structure reveal the precisely designed conformational changes of the chromophore induced by combinatory mutations that shrink its π-conjugated system which, together with electrostatic tuning, produce large blue shifts of the absorption spectra by maximally 100 nm, while maintaining photosensitive ion transport activities. The design principle we elaborate is applicable to other microbial opsins, and clarifies the underlying molecular mechanism of the blue-shifted action spectra of microbial opsins recently isolated from natural sources.

Identical Hydrogen-Bonding Strength of the Retinal Schiff Base between Primate Green- and Red-Sensitive Pigments: New Insight into Color Tuning Mechanism

Three aspects are generally considered in the color-tuning mechanism of vision: (I) chromophore distortion, (II) electrostatic interaction between the protonated Schiff base and counterion, and (III) polarity around the β-ionone ring and polyene chain. Primate green- and red-sensitive proteins are highly homologous but display maximum absorption at 530 and 560 nm, respectively. In the present study, the N-D stretching frequency of monkey green-sensitive protein was identified by using C15-D retinal. The hydrogen-bonding strength between monkey green and red was identical. Together with a previous resonance Raman study, we conclude that the 30 nm difference originates exclusively from the polarity around the β-ionone ring and polyene chain. Three amino acids (Ala, Phe, and Ala in monkey green and Ser, Tyr, and Thr in monkey red, respectively) may be responsible for color tuning together with protein-bound water molecules around the β-ionone ring and polyene chain but not at the Schiff base region.

A color-determining amino acid residue of proteorhodopsin

Proteorhodopsin (PR) is a light-driven proton pump found in marine bacteria. More than 1000 PRs are classified as blue-absorbing (λmax ∼ 490 nm) and green-absorbing (λmax ∼ 525 nm) PRs. The color determinant is known to be at position 105, where blue-absorbing and green-absorbing PRs possess Gln and Leu, respectively. This suggests hydrophobicity at position 105 plays a key role in color tuning. Here we successfully introduced 19 amino acid residues into position 105 of green-absorbing PR in the membrane environment and investigated the absorption properties. High-performance liquid chromatography analysis shows that the isomeric composition of the all-trans form is >70% for all mutants, indicating little influence of different isomers on color tuning. Absorption spectra of the wild-type and 19 mutant proteins were well-characterized by the pH-dependent equilibria of the protonated and deprotonated counterion (Asp97) of the Schiff base, whereas the λmax values of these two states and the pKa value differed significantly among mutants. Although Gln and Leu are hydrophilic and hydrophobic residues, respectively, the λmax values of the two states and the pKa value did not correlate with the hydropathy index of residues. In contrast, the λmax and pKa were correlated with the volume of residues, though Gln and Leu possess similar volumes. This observation concludes that the λmax and pKa of Asp97 are determined by local and specific interactions in the Schiff base moiety, in which the volume of the residue at position 105 is more influential than its hydrophobicity. We suggest that the hydrogen-bonding network in the Schiff base moiety plays a key role in the λmax and pKa of Asp97, and the hydrogen-bonding network is significantly perturbed by large amino acid residues but may be preserved by additional water molecule(s) for small amino acid residues at position 105.

Directed evolution of Gloeobacter violaceus rhodopsin spectral properties

Proton-pumping rhodopsins (PPRs) are photoactive retinal-binding proteins that transport ions across biological membranes in response to light. These proteins are interesting for light-harvesting applications in bioenergy production, in optogenetics applications in neuroscience, and as fluorescent sensors of membrane potential. Little is known, however, about how the protein sequence determines the considerable variation in spectral properties of PPRs from different biological niches or how to engineer these properties in a given PPR. Here we report a comprehensive study of amino acid substitutions in the retinal-binding pocket of Gloeobacter violaceus rhodopsin (GR) that tune its spectral properties. Directed evolution generated 70 GR variants with absorption maxima shifted by up to ±80nm, extending the protein's light absorption significantly beyond the range of known natural PPRs. While proton-pumping activity was disrupted in many of the spectrally shifted variants, we identified single tuning mutations that incurred blue and red shifts of 42nm and 22nm, respectively, that did not disrupt proton pumping. Blue-shifting mutations were distributed evenly along the retinal molecule while red-shifting mutations were clustered near the residue K257, which forms a covalent bond with retinal through a Schiff base linkage. Thirty eight of the identified tuning mutations are not found in known microbial rhodopsins. We discovered a subset of red-shifted GRs that exhibit high levels of fluorescence relative to the WT (wild-type) protein.

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.

Learning from photobiology how to design molecular devices using a computer

Learning how to model photo-responsive proteins may open the way to the design of lightpowered biomimetic molecular devices.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.