Quantum mechanical/molecular mechanical studies on spectral tuning mechanisms of visual pigments and other photoactive proteins

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.

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

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

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Bile Acid Binding Protein Functionalization Leads to a Fully Synthetic Rhodopsin Mimic

Rhodopsins are photoreceptive proteins using light to drive a plethora of biological functions such as vision, proton and ion pumping, cation and anion channeling, and gene and enzyme regulation. Here we combine organic synthesis, NMR structural studies, and photochemical characterization to show that it is possible to prepare a fully synthetic mimic of rhodopsin photoreceptors. More specifically, we conjugate a bile acid binding protein with a synthetic mimic of the rhodopsin protonated Schiff base chromophore to achieve a covalent complex featuring an unnatural protein host, photoswitch, and photoswitch-protein linkage with a reverse orientation. We show that, in spite of its molecular-level diversity, light irradiation of the prepared mimic fuels a photochromic cycle driven by sequential photochemical and thermal Z/E isomerizations reminiscent of the photocycles of microbial rhodopsins.

Excitation spectra of retinal by multiconfiguration pair-density functional theory

Retinal is the chromophore in proteins responsible for vision. The absorption maximum of retinal is sensitive to mutations of the protein. However, it is not easy to predict the absorption spectrum of retinal accurately, and questions remain even after intensive investigation. Retinal poses a challenge for Kohn-Sham density functional theory (KS-DFT) because of the charge transfer character in its excitations, and it poses a challenge for wave function theory because the large size of the molecule makes multiconfigurational perturbation theory methods expensive. In this study, we demonstrate that multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to predict the vertical excitation energies of 11-Z retinal, and it reproduces the experimentally determined absorption band widths and peak positions better than complete active space second-order perturbation theory (CASPT2). The consistency between complete active space self-consistent field (CASSCF) and KS-DFT dipole moments is demonstrated to be a useful criterion in selecting the active space. We also found that the nature of the terminal groups and the conformations of retinal play a significant role in the absorption spectrum. By considering a thermal distribution of conformations, we predict an absorption spectrum of retinal that is consistent with the experimental gas-phase spectrum. The location of the absorption peak and the spectral broadening based on MC-PDFT calculations agree better with experiments than those of CASPT2.

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.

Microbial rhodopsins: wide distribution, rich diversity and great potential

One of the major topics in biophysics and physicobiology is to understand and utilize biological functions using various advanced techniques. Taking advantage of the photoreactivity of the seven-transmembrane rhodopsin protein family has been actively investigated by a variety of methods. Rhodopsins serve as models for membrane-embedded proteins, for photoactive proteins and as a fundamental tool for optogenetics, a new technology to control biological activity with light. In this review, we summarize progress of microbial rhodopsin research from the viewpoint of distribution, diversity and potential.

Pre-gating conformational changes in the ChETA variant of channelrhodopsin-2 monitored by nanosecond IR spectroscopy

Light-gated ion permeation by channelrhodopsin-2 (ChR2) relies on the photoisomerization of the retinal chromophore and the subsequent photocycle, leading to the formation (on-gating) and decay (off-gating) of the conductive state. Here, we have analyzed the photocycle of a fast-cycling ChR2 variant (E123T mutation, also known as ChETA), by time-resolved UV/vis, step-scan FT-IR, and tunable quantum cascade laser IR spectroscopies with nanosecond resolution. Pre-gating conformational changes rise with a half-life of 200 ns, silent to UV/vis but detected by IR spectroscopy. They involve changes in the peptide backbone and in the H-bond of the side chain of the critical residue D156. Thus, the P1(500) intermediate must be separated into early and late states. Light-adapted ChR2 contains a mixture of all-trans and 13-cis retinal in a 70:30 ratio which are both photoactive. Analysis of ethylenic and fingerprint vibrations of retinal provides evidence that the 13-cis photocycle recovers in 1 ms. This recovery is faster than channel off-gating and most of the proton transfer reactions, implying that the 13-cis photocycle is of minor functional relevance for ChR2.

Gloeobacter rhodopsin, limitation of proton pumping at high electrochemical load

We studied the photocurrents of a cyanobacterial rhodopsin Gloeobacter violaceus (GR) in Xenopus laevis oocytes and HEK-293 cells. This protein is a light-driven proton pump with striking similarities to marine proteorhodopsins, including the D121-H87 cluster of the retinal Schiff base counterion and a glutamate at position 132 that acts as a proton donor for chromophore reprotonation during the photocycle. Interestingly, at low extracellular pH(o) and negative voltage, the proton flux inverted and directed inward. Using electrophysiological measurements of wild-type and mutant GR, we demonstrate that the electrochemical gradient limits outward-directed proton pumping and converts it into a purely passive proton influx. This conclusion contradicts the contemporary paradigm that at low pH, proteorhodopsins actively transport H(+) into cells. We identified E132 and S77 as key residues that allow inward directed diffusion. Substitution of E132 with aspartate or S77 with either alanine or cysteine abolished the inward-directed current almost completely. The proton influx is likely caused by the pK(a) of E132 in GR, which is lower than that of other microbial ion pumping rhodopsins. The advantage of such a low pK(a) is an acceleration of the photocycle and high pump turnover at high light intensities.

Spectral "Fine" Tuning in Fluorescent Proteins: The Case of the GFP-Like Chromophore in the Anionic Protonation State

Fluorescent proteins (FPs), featuring the same chromophore but different chromophore-protein interactions, display remarkable spectral variations even when the same chromophore protonation state, i.e. the anionic state, is involved. We examine the mechanisms behind this tuning by means of structural analysis, molecular dynamics simulations, and vertical excitation energy calculations using QM/MM Time-Dependent Density Functional Theory (TD-DFT), CASPT2/CASSCF, and SAC-CI. The proteins under investigation include the structurally similar, though spectrally distinct, Dronpa and mTFP0.7, with absorption peaks at 453 and 503 nm, respectively. We extend our analysis to two Green Fluorescent Protein variants, GFP-S65T (absorption peak at 484 nm), for comparison with previous computational studies, and GFP-S65G/V68L/S72A/T203Y, a yellow fluorescent protein (514 nm), in order to include one of the most red-shifted FPs containing a GFP-like chromophore. We compare different choices of the QM system, and we discuss how molecular dynamics simulations affect the calculation of excitation energies, with respect to X-ray structures. We are able to partially reproduce the spectral tuning of the FPs and correlate it to the chromophore bond-length variations, as determined by specific interactions with the chromophore environment.

Full-Quantum chemical calculation of the absorption maximum of bacteriorhodopsin: a comprehensive analysis of the amino acid residues contributing to the opsin shift

Herein, the absorption maximum of bacteriorhodopsin (bR) is calculated using our recently developed method in which the whole protein can be treated quantum mechanically at the level of INDO/S-CIS//ONIOM (B3LYP/6-31G(d,p): AMBER). The full quantum mechanical calculation is shown to reproduce the so-called opsin shift of bR with an error of less than 0.04 eV. We also apply the same calculation for 226 different bR mutants, each of which was constructed by replacing any one of the amino acid residues of the wild-type bR with Gly. This substitution makes it possible to elucidate the extent to which each amino acid contributes to the opsin shift and to estimate the inter-residue synergistic effect. It was found that one of the most important contributions to the opsin shift is the electron transfer from Tyr185 to the chromophore upon excitation. We also indicate that some aromatic (Trp86, Trp182) and polar (Ser141, Thr142) residues, located in the vicinity of the retinal polyene chain and the β-ionone ring, respectively, play an important role in compensating for the large blue-shift induced by both the counterion residues (Asp85, Asp212) and an internal water molecule (W402) located near the Schiff base linkage. In particular, the effect of Trp86 is comparable to that of Tyr185. In addition, Ser141 and Thr142 were found to contribute to an increase in the dipole moment of bR in the excited state. Finally, we provide a complete energy diagram for the opsin shift together with the contribution of the chromophore-protein steric interaction.

A single enhancer regulating the differential expression of duplicated red-sensitive opsin genes in zebrafish

A fundamental step in the evolution of the visual system is the gene duplication of visual opsins and differentiation between the duplicates in absorption spectra and expression pattern in the retina. However, our understanding of the mechanism of expression differentiation is far behind that of spectral tuning of opsins. Zebrafish (Danio rerio) have two red-sensitive cone opsin genes, LWS-1 and LWS-2. These genes are arrayed in a tail-to-head manner, in this order, and are both expressed in the long member of double cones (LDCs) in the retina. Expression of the longer-wave sensitive LWS-1 occurs later in development and is thus confined to the peripheral, especially ventral-nasal region of the adult retina, whereas expression of LWS-2 occurs earlier and is confined to the central region of the adult retina, shifted slightly to the dorsal-temporal region. In this study, we employed a transgenic reporter assay using fluorescent proteins and P1-artificial chromosome (PAC) clones encompassing the two genes and identified a 0.6-kb “LWS-activating region” (LAR) upstream of LWS-1, which regulates expression of both genes. Under the 2.6-kb flanking upstream region containing the LAR, the expression pattern of LWS-1 was recapitulated by the fluorescent reporter. On the other hand, when LAR was directly conjugated to the LWS-2 upstream region, the reporter was expressed in the LDCs but also across the entire outer nuclear layer. Deletion of LAR from the PAC clones drastically lowered the reporter expression of the two genes. These results suggest that LAR regulates both LWS-1 and LWS-2 by enhancing their expression and that interaction of LAR with the promoters is competitive between the two genes in a developmentally restricted manner. Sharing a regulatory region between duplicated genes could be a general way to facilitate the expression differentiation in duplicated visual opsins.

Among vertebrates, fish may have the most advanced color vision. They have greatly varied repertoires of color sensors called visual opsins, possibly reflecting evolutionary adaptation to their diverse photic environments in water, and are an excellent model to study the evolution of vertebrate color vision. This is achieved by multiplying opsin genes and differentiating their absorption light spectra and expression patterns. However, little is understood regarding how the opsin genes are regulated to achieve the differential expression pattern. In this study, we focused on the duplicated red-sensitive opsin genes of zebrafish to tackle this problem. We discovered an “enhancer” region near the two red opsin genes that plays a crucial role in their differential expression pattern. Our results suggest that the two red opsin genes interact with the enhancer competitively in a developmentally restricted manner. Sharing a regulatory region could be a general way to facilitate the expression differentiation in duplicated visual opsin genes.

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

We have investigated the protonation state and photoabsorption spectrum of Schiff-base (SB) nitrogen bound 11-cis-retinal in human blue and mouse UV cone visual pigments as well as in bovine rhodopsin by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have employed both multireference (MRCISD+Q, MR-SORCI+Q, and MR-DDCI2+Q) and single reference (TD-B3LYP and RI-CC2) QM methods. The calculated ground-state and vertical excitation energies show that UV-sensitive pigments have deprotonated SB nitrogen, while violet-sensitive pigments have protonated SB nitrogen, in agreement with some indirect experimental evidence. A significant blue shift of the absorption maxima of violet-sensitive pigments relative to rhodopsins arises from the increase in bond length alternation of the polyene chain of 11-cis-retinal induced by polarizing fields of these pigments. The main counterion is Glu113 in both violet-sensitive vertebrate pigments and bovine rhodopsin. Neither Glu113 nor the remaining pigment has a significant influence on the first excitation energy of 11-cis-retinal in the UV-sensitive pigments that have deprotonated SB nitrogen. There is no charge transfer between the SB and beta-ionone terminals of 11-cis-retinal in the ground and first excited states.