Color tuning in short wavelength-sensitive human and mouse visual pigments: ab initio quantum mechanics/molecular mechanics studies

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.

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.

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.

Adaptive Landscapes in the Age of Synthetic Biology

For nearly a century adaptive landscapes have provided overviews of the evolutionary process and yet they remain metaphors. We redefine adaptive landscapes in terms of biological processes rather than descriptive phenomenology. We focus on the underlying mechanisms that generate emergent properties such as epistasis, dominance, trade-offs and adaptive peaks. We illustrate the utility of landscapes in predicting the course of adaptation and the distribution of fitness effects. We abandon aged arguments concerning landscape ruggedness in favor of empirically determining landscape architecture. In so doing, we transform the landscape metaphor into a scientific framework within which causal hypotheses can be tested.

Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations

The initial S₁ excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S₁ energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.

Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal?

To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S₀ →S₁ and S₀ →S₂ energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.

The Dynamic Origin of Color Tuning in Proteins Revealed by a Carotenoid Pigment

Understanding the microscopic origin of the color tuning in pigment-protein complexes is a challenging yet fundamental issue in photoactive biological systems. Here, we propose a possible interpretation by using a state-of-the-art multiscale strategy based on the integration of quantum chemistry and polarizable atomistic embeddings into a dynamic description. By means of such a strategy we are able to resolve the long-standing dispute over the coloration mechanism in the crustacyanin protein. It is shown that the combination of the dynamical flexibility of the carotenoid pigments (astaxanthin) with the responsive protein environment is essential to obtain quantitative predictions of the spectral tuning. The strong linear correlation between the excitation energies and the bond length alternation in the long-chain carotenoids modulated by the dynamical protein environment is a novel finding explaining the high color tunability in crustacyanin.

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.

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.

A simple method for studying the molecular mechanisms of ultraviolet and violet reception in vertebrates

BACKGROUND

Many vertebrate species use ultraviolet (UV) reception for such basic behaviors as foraging and mating, but many others switched to violet reception and improved their visual resolution. The respective phenotypes are regulated by the short wavelength-sensitive (SWS1) pigments that absorb light maximally (λmax) at ~360 and 395-440 nm. Because of strong epistatic interactions, the biological significance of the extensive mutagenesis results on the molecular basis of spectral tuning in SWS1 pigments and the mechanisms of their phenotypic adaptations remains uncertain.

RESULTS

The magnitudes of the λmax-shifts caused by mutations in a present-day SWS1 pigment and by the corresponding forward mutations in its ancestral pigment are often dramatically different. To resolve these mutagenesis results, the A/B ratio, in which A and B are the areas formed by amino acids at sites 90, 113 and 118 and by those at sites 86, 90 and 118 and 295, respectively, becomes indispensable. Then, all critical mutations that generated the λmax of a SWS1 pigment can be identified by establishing that 1) the difference between the λmax of the ancestral pigment with these mutations and that of the present-day pigment is small (3 ~ 5 nm, depending on the entire λmax-shift) and 2) the difference between the corresponding A/B ratios is < 0.002.

CONCLUSION

Molecular adaptation has been studied mostly by using comparative sequence analyses. These statistical results provide biological hypotheses and need to be tested using experimental means. This is an opportune time to explore the currently available and new genetic systems and test these statistical hypotheses. Evaluating the λmaxs and A/B ratios of mutagenized present-day and their ancestral pigments, we now have a method to identify all critical mutations that are responsible for phenotypic adaptation of SWS1 pigments. The result also explains spectral tuning of the same pigments, a central unanswered question in phototransduction.

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.

Epistatic adaptive evolution of human color vision

Establishing genotype-phenotype relationship is the key to understand the molecular mechanism of phenotypic adaptation. This initial step may be untangled by analyzing appropriate ancestral molecules, but it is a daunting task to recapitulate the evolution of non-additive (epistatic) interactions of amino acids and function of a protein separately. To adapt to the ultraviolet (UV)-free retinal environment, the short wavelength-sensitive (SWS1) visual pigment in human (human S1) switched from detecting UV to absorbing blue light during the last 90 million years. Mutagenesis experiments of the UV-sensitive pigment in the Boreoeutherian ancestor show that the blue-sensitivity was achieved by seven mutations. The experimental and quantum chemical analyses show that 4,008 of all 5,040 possible evolutionary trajectories are terminated prematurely by containing a dehydrated nonfunctional pigment. Phylogenetic analysis further suggests that human ancestors achieved the blue-sensitivity gradually and almost exclusively by epistasis. When the final stage of spectral tuning of human S1 was underway 45-30 million years ago, the middle and long wavelength-sensitive (MWS/LWS) pigments appeared and so-called trichromatic color vision was established by interprotein epistasis. The adaptive evolution of human S1 differs dramatically from orthologous pigments with a major mutational effect used in achieving blue-sensitivity in a fish and several mammalian species and in regaining UV vision in birds. These observations imply that the mechanisms of epistatic interactions must be understood by studying various orthologues in different species that have adapted to various ecological and physiological environments.

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.

The effect of protein environment on photoexcitation properties of retinal

Retinal is the photon absorbing chromophore of rhodopsin and other visual pigments, enabling the vertebrate vision process. The effects of the protein environment on the primary photoexcitation process of retinal were studied by time-dependent density functional theory (TDDFT) and the algebraic diagrammatic construction through second order (ADC(2)) combined with our recently introduced reduction of virtual space (RVS) approximation method. The calculations were performed on large full quantum chemical cluster models of the bluecone (BC) and rhodopsin (Rh) pigments with 165-171 atoms. Absorption wavelengths of 441 and 491 nm were obtained at the B3LYP level of theory for the respective models, which agree well with the experimental values of 414 and 498 nm. Electrostatic rather than structural strain effects were shown to dominate the spectral tuning properties of the surrounding protein. The Schiff base retinal and a neighboring Glu-113 residue were found to have comparable proton affinities in the ground state of the BC model, whereas in the excited state, the proton affinity of the Schiff base is 5.9 kcal/mol (0.26 eV) higher. For the ground and excited states of the Rh model, the proton affinity of the Schiff base is 3.2 kcal/mol (0.14 eV) and 7.9 kcal/mol (0.34 eV) higher than for Glu-113, respectively. The protein environment was found to enhance the bond length alternation (BLA) of the retinyl chain and blueshift the first absorption maxima of the protonated Schiff base in the BC and Rh models relative to the chromophore in the gas phase. The protein environment was also found to decrease the intensity of the second excited state, thus improving the quantum yield of the photoexcitation process. Relaxation of the BC model on the excited state potential energy surface led to a vanishing BLA around the isomerization center of the conjugated retinyl chain, rendering the retinal accessible for cis-trans isomerization. The energy of the relaxed excited state was found to be 30 kcal/mol (1.3 eV) above the minimum ground state energy, and might be related to the transition state of the thermal activation process.

Reduction of the virtual space for coupled-cluster excitation energies of large molecules and embedded systems

We investigate how the reduction of the virtual space affects coupled-cluster excitation energies at the approximate singles and doubles coupled-cluster level (CC2). In this reduced-virtual-space (RVS) approach, all virtual orbitals above a certain energy threshold are omitted in the correlation calculation. The effects of the RVS approach are assessed by calculations on the two lowest excitation energies of 11 biochromophores using different sizes of the virtual space. Our set of biochromophores consists of common model systems for the chromophores of the photoactive yellow protein, the green fluorescent protein, and rhodopsin. The RVS calculations show that most of the high-lying virtual orbitals can be neglected without significantly affecting the accuracy of the obtained excitation energies. Omitting all virtual orbitals above 50 eV in the correlation calculation introduces errors in the excitation energies that are smaller than 0.1 eV. By using a RVS energy threshold of 50 eV, the CC2 calculations using triple-ζ basis sets (TZVP) on protonated Schiff base retinal are accelerated by a factor of 6. We demonstrate the applicability of the RVS approach by performing CC2/TZVP calculations on the lowest singlet excitation energy of a rhodopsin model consisting of 165 atoms using RVS thresholds between 20 eV and 120 eV. The calculations on the rhodopsin model show that the RVS errors determined in the gas-phase are a very good approximation to the RVS errors in the protein environment. The RVS approach thus renders purely quantum mechanical treatments of chromophores in protein environments feasible and offers an ab initio alternative to quantum mechanics/molecular mechanics separation schemes.

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.

H-bond network around retinal regulates the evolution of ultraviolet and violet vision

Ancestors of vertebrates used ultraviolet vision. Some descendants preserved ultraviolet vision, whereas some others replaced it with violet vision, and then, some of avian lineages reinvented ultraviolet vision. Ultraviolet (absorption at ∼360 nm) and violet (410-440 nm) sensitivities of visual pigments are known to be affected by around 20 amino acid substitutions. The present quantum mechanical/molecular mechanical calculations show that these substitutions modify a H-bond network formed by two waters and sites 86, 90, 113, 114, 118, and 295, which determines the protonation state of Schiff base linked 11-cis-retinal. A pigment is ultraviolet-sensitive when it is more stable with an unprotonated retinal (SBR) form than with its protonated analogue (PSBR) and is violet-sensitive when the PSBR form is more stable. These results establish for the first time the chemical basis of ultraviolet and violet vision in vertebrates.

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

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

Axial chirality of donor-donor, donor-acceptor, and tethered 1,1'-binaphthyls: a theoretical revisit with dynamics trajectories

The circular dichroism (CD) spectra of (R)-2,2'-dimethoxy-1,1'-binaphthyl (DD) and its untethered and tethered donor-acceptor analogues (DA and DA7-DA9) were investigated experimentally and theoretically. The experimental CD spectra of DD and DA resembled each other in several aspects, displaying a positive-positive-negative Cotton effect pattern in the (1)L(b)-(1)L(a) region and a strong negative couplet at the (1)B(b) band, but significantly differed in transition energy and rotatory strength. The couplet amplitude (A) of the main band was 1.6 times larger in DA than in DD, despite the comparable extinction coefficients and seemingly analogous conformations. An additional positive Cotton effect was observed at the CT (CT) band for donor-acceptor binaphthyl DA. Our theoretical prediction of the CD spectra of binaphthyls involves three sequential first principle quantum mechanics (QM) calculations. Thus, the geometry optimizations of a series of conformers with varying dihedral angles were performed by the dispersion-corrected DFT-D method using the B97-D functional and the TZV2P basis set. The potential curve as a function of the dihedral angle (θ) was obtained by using the SCS-MP2/TZVPP single-point energy calculations with and without application of the solvent correction. The CD spectrum of each conformer was independently calculated by the second-order approximate coupled cluster calculation (CC2 method) using the TZVPP basis sets and the resolution of the identity (RI-J) approximation. The (net) theoretical CD spectrum was obtained by averaging over all possible conformers, where the dynamics trajectories based on the relative SCS-MP2 energies were taken into account. By using 17 possible conformers at θ varying from 50 to 130° by 5° intervals, the experimental CD spectra were successfully reproduced in a quantitative manner, enabling us to characterize properly almost all of the important spectral features and chiroptical properties. The two-state model, reported previously, turned out to have led to the right answer with wrong reasons. The couplet sign and amplitude A are critical functions of θ and can be used not only for (qualitatively) determining the absolute configuration but also for quantitatively analyzing the binaphthyl conformations. The angle dependence of A was already argued in the classical coupled oscillator and exciton chirality theories to provide reasonable structure elucidations but only in a qualitative or semiquantitative manner. Our method is able to predict the A value quantitatively as a function of θ. For tethered binaphthyls DA7-DA9, particular care should be exercised in the conformational assessment based on the classical treatment because the amplitude A was shown to be significantly affected by the existence of the tether itself. In the present method, the couplet amplitude A was nicely related to the dihedral angle θ of DA and DD by the state-of-the-art ab initio calculations, enabling us to gain the quantitative information about the conformation of axially chiral binaphthyls. The Cotton effect at the CT band also serves as a complementary clue for elucidating the conformation of donor-acceptor binaphthyls.

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

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

Studies of the ground and excited-state surfaces of the retinal chromophore using CAM-B3LYP

The isomerization of the 11-cis isomer (PSB11) of the retinal chromophore to its all-trans isomer (PSBT) is examined. Optimized structures on both the ground state and the excited state are calculated, and the dependence on torsional angles in the carbon chain is investigated. Time-dependent density functional theory is used to produce excitation energies and the excited-state surface. To avoid problems with the description of excited states that can arise with standard DFT methods, the CAM-B3LYP functional was used. Comparing CAM-B3LYP with B3LYP results indicates that the former is significantly more accurate, as a consequence of which detailed cross sections of the retinal excited-state surface are obtained.