How the Counterion Affects Ground- and Excited-State Properties of the Rhodopsin Chromophore

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.

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.

RubyACRs, nonalgal anion channelrhodopsins with highly red-shifted absorption

Channelrhodopsins are light-gated ion channels widely used to control neuronal firing with light (optogenetics). We report two previously unknown families of anion channelrhodopsins (ACRs), one from the heterotrophic protists labyrinthulea and the other from haptophyte algae. Four closely related labyrinthulea ACRs, named RubyACRs here, exhibit a unique retinal-binding pocket that creates spectral sensitivities with maxima at 590 to 610 nm, the most red-shifted channelrhodopsins known, long-sought for optogenetics, and more broadly the most red-shifted microbial rhodopsins thus far reported. We identified three spectral tuning residues critical for the red-shifted absorption. Photocurrents recorded from the RubyACR from Aurantiochytrium limacinum (designated AlACR1) under single-turnover excitation exhibited biphasic decay, the rate of which was only weakly voltage dependent, in contrast to that in previously characterized cryptophyte ACRs, indicating differences in channel gating mechanisms between the two ACR families. Moreover, in A. limacinum we identified three ACRs with absorption maxima at 485, 545, and 590 nm, indicating color-sensitive photosensing with blue, green, and red spectral variation of ACRs within individual species of the labyrinthulea family. We also report functional energy transfer from a cytoplasmic fluorescent protein domain to the retinal chromophore bound within RubyACRs.

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.

Theoretical Studies on the Color-Tuning Mechanism in Retinal Proteins

The excited states of the three retinal proteins, bovine rhodopsin (Rh), bacteriorhodopsin (bR), and sensory rhodopsin II (sRII) were studied using the symmetry-adapted cluster-configuration interaction (SAC-CI) and combined quantum mechanical and molecular mechanical (QM/MM) methods. The computed absorption energies are in good agreement with the experimental ones for all three proteins. The spectral tuning mechanism was analyzed in terms of three contributions:  molecular structures of the chromophore in the binding pockets, electrostatic (ES) interaction of the chromophore with the surrounding protein environment, and quantum-mechanical effect between the chromophore and the counterion group. This analysis provided an insight into the mechanism of the large blue-shifts in the absorption peak position of Rh and sRII from that of bR. Protein ES effect is primarily important both in Rh and in sRII, and the structure effect is secondary important in Rh. The quantum-mechanical interaction between the chromophore and the counterion is very important for quantitative reproduction of the excitation energy. These results indicate that the present approach is useful for studying the absorption spectra and the mechanism of the color tuning in the retinal proteins.

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.

Protein Field Effect on the Dark State of 11-cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics

The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.

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

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

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.

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.

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.

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.

Effect of polarization on the opsin shift in rhodopsins. 2. Empirical polarization models for proteins

The explicit treatment of polarization as a many-body interaction in condensed-phase systems represents a current problem in empirical force-field development. Although a variety of efficient models for molecular polarization have been suggested, polarizable force fields are still far from common use nowadays. In this work, we consider interactive polarization models employing Thole's short-range damping scheme and assess them for application on polypeptides. Despite the simplicity of the model, we find mean polarizabilities and anisotropies of amino acid side chains in excellent agreement with MP2/cc-pVQZ benchmark calculations. Combined with restrained electrostatic potential (RESP) derived atomic charges, the models are applied in a quantum-mechanical/molecular-mechanical (QM/MM) approach. An iterative scheme is used to establish a self-consistent mutual polarization between the QM and MM moieties. This ansatz is employed to study the influence of the protein polarizability on calculated optical properties of the protonated Schiff base of retinal in rhodopsin (Rh), bacterio-rhodopsin (bR), and pharaonis sensory rhodopsin II (psRII). The shifts of the excitation energy due to the instantaneous polarization response of the protein to the charge transfer on the retinal chromophore are quantified using the high level ab initio multireference spectroscopy-oriented configuration interaction (SORCI) method. The results are compared with those of previously published QM1/QM2/MM models for bR and psRII.

Effect of polarization on the opsin shift in rhodopsins. 1. A combined QM/QM/MM model for bacteriorhodopsin and pharaonis sensory rhodopsin II

The optical and IR-spectroscopic properties of the protonated Schiff base of retinal are highly sensitive to the electrostatic environment. This feature makes retinal a useful probe to study structural differences and changes in rhodopsins. It also raises an interest to theoretically predict the spectroscopic response to mutation and structural evolution. Computational models appropriate for this purpose usually combine sophisticated quantum mechanical (QM) methods with molecular mechanics (MM) force fields. In an effort to test and improve the accuracy of these QM/MM models, we consider in this article the effects of polarization and inter-residual charge transfer within the binding pocket of bacteriorhodopsin (bR) and pharaonis sensory rhodopsin II (psRII, also called pharaonis phoborhodopsin, ppR) on the excitation energy using an ab initio QM/QM/MM approach. The results will serve as reference for assessing empirical polarization models in a consecutive article.

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.

An opsin shift in rhodopsin: retinal S0-S1 excitation in protein, in solution, and in the gas phase

We considered a series of model systems for treating the photoabsorption of the 11-cis retinal chromophore in the protonated Schiff-base form in vacuum, solutions, and the protein environment. A high computational level, including the quantum mechanical-molecular mechanical (QM/MM) approach for solution and protein was utilized in simulations. The S0-S1 excitation energies in quantum subsystems were evaluated by means of an augmented version of the multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) with the ground-state geometry parameters optimized in the density functional theory PBE0/cc-pVDZ approximation. The computed positions of absorption bands lambdamax, 599(g), 448(s), and 515(p) nm for the gas phase, solution, and protein, respectively, are in excellent agreement with the corresponding experimental data, 610(g), 445(s), and 500(p) nm. Such consistency provides a support for the formulated qualitative conclusions on the role of the chromophore geometry, environmental electrostatic field, and the counterion in different media. An essentially nonplanar geometry conformation of the chromophore group in the region of the C14-C15 bond was obtained for the protein, in particular, owing to the presence of the neighboring charged amino acid residue Glu181. Nonplanarity of the C14-C15 bond region along with the influence of the negatively charged counterions Glu181 and Glu113 are found to be important to reproduce the spectroscopic features of retinal chromophore inside the Rh cavity. Furthermore, the protein field is responsible for the largest bond-order decrease at the C11-C12 double bond upon excitation, which may be the reason for the 11-cis photoisomerization specificity.

Color tuning in rhodopsins: the mechanism for the spectral shift between bacteriorhodopsin and sensory rhodopsin II

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.

Photoisomerization efficiency in UV-absorbing visual pigments: protein-directed isomerization of an unprotonated retinal Schiff base

A visual pigment consists of an opsin protein and a chromophore, 11-cis-retinal, which binds to a specific lysine residue of opsin via a Schiff base linkage. The Schiff base chromophore is protonated in pigments that absorb visible light, whereas it is unprotonated in ultraviolet-absorbing visual pigments (UV pigments). To investigate whether an unprotonated Schiff base can undergo photoisomerization as efficiently as a protonated Schiff base in the opsin environment, we measured the quantum yields of the bovine rhodopsin E113Q mutant, in which the Schiff base is unprotonated at alkaline pH, and the mouse UV pigment (mouse UV). Photosensitivities of UV pigments were measured by irradiation of the pigments followed by chromophore extraction and HPLC analysis. Extinction coefficients were estimated by comparing the maximum absorbances of the original pigments and their acid-denatured states. The quantum yield of the bovine rhodopsin E113Q mutant at pH 8.2, where the Schiff base is unprotonated, was significantly lower than that of wild-type rhodopsin, whereas the mutant gave a quantum yield almost identical to that of the wild type at pH 5.5, where the Schiff base is protonated. These results suggest that Schiff base protonation plays a role in increasing quantum yield. The quantum yield of mouse UV, which has an unprotonated Schiff base chromophore, was significantly higher than that of the unprotonated form of the rhodopsin E113Q mutant, although it was still lower than the visible-absorbing pigments. These results suggest that the mouse UV pigment has a specific mechanism for the efficient photoisomerization of its unprotonated Schiff base chromophore.

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Quantum Chemical Methods for the Investigation of Photoinitiated Processes in Biological Systems: Theory and Applications

With the advent of modern computers and advances in the development of efficient quantum chemical computer codes, the meaningful computation of large molecular systems at a quantum mechanical level became feasible. Recent experimental effort to understand photoinitiated processes in biological systems, for instance photosynthesis or vision, at a molecular level also triggered theoretical investigations in this field. In this Minireview, standard quantum chemical methods are presented that are applicable and recently used for the calculation of excited states of photoinitiated processes in biological molecular systems. These methods comprise configuration interaction singles, the complete active space self-consistent field method, and time-dependent density functional theory and its variants. Semiempirical approaches are also covered. Their basic theoretical concepts and mathematical equations are briefly outlined, and their properties and limitations are discussed. Recent successful applications of the methods to photoinitiated processes in biological systems are described and theoretical tools for the analysis of excited states are presented.

Excited-State Properties and Environmental Effects for Protonated Schiff Bases: A Theoretical Study

Complete active space self-consistent field (CASSCF), multireference configuration interaction (MRCI), density functional theory (DFT), time dependent DFT (TDDFT) and the singles and doubles coupled-cluster (CC2) methodologies have been used to study the ground state and excited states of protonated and neutral Schiff bases (PSB and SB) as models for the retinal chromophore. Systems with two to four conjugated double bonds are investigated. Geometry relaxation effects are studied in the excited pipi* state using the aforementioned methods. Taking the MRCI results as reference we find that CASSCF results are quite reliable even though overshooting of geometry changes is observed. TDDFT does not reproduce bond alternation well in the pipi* state. CC2 takes an intermediate position. Environmental effects due to solvent or protein surroundings have been studied in the excited states of the PSBs and SBs using a water molecule and solvated formate as model cases. Particular emphasis is given to the proton transfer process from the PSB to its solvent partner in the excited state. It is found that its feasibility is significantly enhanced in the excited state as compared to the ground state, which means that a proton transfer could be initiated already at an early step in the photodynamics of PSBs.

Solvent Effects on the Low-Lying Excited States of a Model of Retinal

The low-lying excited states of a solution in alcohol of a five-double-bond model of the rhodopsin protein chromophore, the protonated 11-cis-retinal Schiff base (PSB11), are studied theoretically. We combine a multireference perturbational treatment in the description of the solute molecule with molecular dynamics calculations in the description of the solvent. The geometry, charge distribution, and electronic spectra are strongly influenced by the solvent. The solvent shift values show a marked dependence on the use of relaxed geometries in solution and on the nature of the states involved in the excitation process. The dynamic correlation has a strong effect on the order of the excited states. In solution, the first two excited states almost become degenerate.

Ground and excited states of retinal schiff base chromophores by multiconfigurational perturbation theory

We have studied the wavelength dependence of retinal Schiff base absorbencies on the protonation state of the chromophore at the multiconfigurational level of theory using second order perturbation theory (CASPT2) within an atomic natural orbital basis set on MP2 optimized geometries. Quantitative agreement between calculated and experimental absorption maxima was obtained for protonated and deprotonated Schiff bases of all-trans- and 11-cis-retinal and intermediate states covering a wavelength range from 610 to 353 nm. These data will be useful as reference points for the calibration of more approximate schemes.

Origin and Consequences of Steric Strain in the Rhodopsin Binding Pocket

To study the origin and the effects of steric strain on the chromophore conformation in rhodopsin, we have performed quantum-mechanical calculations on the wild-type retinal chromophore and four retinal derivatives, 13-demethyl-, 10-methyl-13-demethyl-, 10-methyl-, and 9-demethylretinal. For the dynamics of the whole protein, a combined quantum mechanics/molecular mechanics method (DFTB/CHARMM) was used and for the calculation of excited-state properties the nonempirical CASSCF/CASPT2 method. After relaxation inside the protein, all chromophores show significant nonplanar distortions from C10 to C13, most strongly for 10-methylretinal and least pronounced for 9-demethylretinal. In all five cases, the dihedral angle of the C10-C11=C12-C13 bond is negative which attests to the strong chiral discrimination exerted by the protein pocket. The calculations show that the nonplanar distortion of the chromophore, including the sense of rotation, is caused by a combination of two effects: the fitting of both ends to the protein matrix which imposes a distance constraint and the bonding arrangement at the Schiff base terminus. With both the counterion Glu113 and Lys296 displaced off the plane of the chromophore, their binding to N16 exerts a torque on the chromophore. As a result, the polyene chain, from N16 to C13, is twisted in a clockwise manner against the remaining part of the chromophore, leading to a C11=C12 bond with the observed negative dihedral angle. Shifts of the absorption maxima are reproduced correctly, in particular, the red shift of the 10-methyl and the strong blue shift of the 9-demethyl analogue relative to the wild type. Calculated positive rotatory strengths of the alpha-CD bands are in agreement with the calculated absolute conformation of the mutant chromophores.