Color Tuning in rhodopsins: the origin of the spectral shift between the chloride-bound and anion-free forms of halorhodopsin

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.

Fluorescence of the Retinal Chromophore in Microbial and Animal Rhodopsins

Fluorescence of the vast majority of natural opsin-based photoactive proteins is extremely low, in accordance with their functions that depend on efficient transduction of absorbed light energy. However, several recently proposed classes of engineered rhodopsins with enhanced fluorescence, along with the discovery of a new natural highly fluorescent rhodopsin, NeoR, opened a way to exploit these transmembrane proteins as fluorescent sensors and draw more attention to studies on this untypical rhodopsin property. Here, we review the available data on the fluorescence of the retinal chromophore in microbial and animal rhodopsins and their photocycle intermediates, as well as different isomers of the protonated retinal Schiff base in various solvents and the gas phase.

Fluorescence Imaging of Cell Membrane Potential: From Relative Changes to Absolute Values

Membrane potential is a fundamental property of biological cells. Changes in membrane potential characterize a vast number of vital biological processes, such as the activity of neurons and cardiomyocytes, tumorogenesis, cell-cycle progression, etc. A common strategy to record membrane potential changes that occur in the process of interest is to utilize organic dyes or genetically-encoded voltage indicators with voltage-dependent fluorescence. Sensors are introduced into target cells, and alterations of fluorescence intensity are recorded with optical methods. Techniques that allow recording relative changes of membrane potential and do not take into account fluorescence alterations due to factors other than membrane voltage are already widely used in modern biological and biomedical studies. Such techniques have been reviewed previously in many works. However, in order to investigate a number of processes, especially long-term processes, the measured signal must be corrected to exclude the contribution from voltage-independent factors or even absolute values of cell membrane potential have to be evaluated. Techniques that enable such measurements are the subject of this review.

Azobenzene/Tetraethyl Ammonium Photochromic Potassium Channel Blockers: Scope and Limitations for Design of Para-Substituted Derivatives with Specific Absorption Band Maxima and Thermal Isomerization Rate

Azobenzene/tetraethyl ammonium photochromic ligands (ATPLs) are photoactive compounds with a large variety of photopharmacological applications such as nociception control or vision restoration. Absorption band maximum and lifetime of the less stable isomer are important characteristics that determine the applicability of ATPLs. Substituents allow to adjust these characteristics in a range limited by the azobenzene/tetraethyl ammonium scaffold. The aim of the current study is to find the scope and limitations for the design of ATPLs with specific spectral and kinetic properties by introducing para substituents with different electronic effects. To perform this task we synthesized ATPLs with various electron acceptor and electron donor functional groups and studied their spectral and kinetic properties using flash photolysis and conventional spectroscopy techniques as well as quantum chemical modeling. As a result, we obtained diagrams that describe correlations between spectral and kinetic properties of ATPLs (absorption maxima of E and Z isomers of ATPLs, the thermal lifetime of their Z form) and both the electronic effect of substituents described by Hammett constants and structural parameters obtained from quantum chemical calculations. The provided results can be used for the design of ATPLs with properties that are optimal for photopharmacological applications.

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.

Insights into the Protein Functions and Absorption Wavelengths of Microbial Rhodopsins

Using a quantum mechanical/molecular mechanical approach, the absorption wavelength of the retinal Schiff base was calculated based on 13 microbial rhodopsin crystal structures. The results showed that the protein electrostatic environment decreases the absorption wavelength significantly in the cation-conducting rhodopsin but only slightly in the sensory rhodopsin. Among the microbial rhodopsins with different functions, the differences in the absorption wavelengths are caused by differences in the arrangement of the charged residues at the retinal Schiff base binding moiety, namely, one or two counterions at the three common positions. Among the microbial rhodopsins with similar functions, the differences in the polar residues at the retinal Schiff base binding site are responsible for the differences in the absorption wavelengths. Counterions contribute to an absorption wavelength shift of 50-120 nm, whereas polar groups contribute to a shift of up to ∼10 nm. It seems likely that protein function is directly associated with the absorption wavelength in microbial rhodopsins.

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.

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.

A Comparative Study of Modern Homology Modeling Algorithms for Rhodopsin Structure Prediction

Rhodopsins are seven α-helical membrane proteins that are of great importance in chemistry, biology, and modern biotechnology. Any in silico study on rhodopsin properties and functioning requires a high-quality three-dimensional structure. Due to particular difficulties with obtaining membrane protein structures from the experiment, in silico prediction of the three-dimensional rhodopsin structure based only on its primary sequence is an especially important task. For the last few years, significant progress was made in the field of protein structure prediction, especially for methods based on comparative modeling. However, the majority of this progress was made for soluble proteins and further investigations are needed to achieve similar progress for membrane proteins. In this paper, we evaluate the performance of modern protein structure prediction methodologies (implemented in the Medeller, I-TASSER, and Rosetta packages) for their ability to predict rhodopsin structures. Three widely used methodologies were considered: two general methodologies that are commonly applied to soluble proteins and a methodology that uses constraints that are specific for membrane proteins. The test pool consisted of 36 target-template pairs with different sequence similarities that was constructed on the basis of 24 experimental rhodopsin structures taken from the RCSB database. As a result, we showed that all three considered methodologies allow obtaining rhodopsin structures with the quality that is close to the crystallographic one (root mean square deviation (RMSD) of the predicted structure from the corresponding X-ray structure up to 1.5 Å) if the target-template sequence identity is higher than 40%. Moreover, all considered methodologies provided structures of average quality (RMSD < 4.0 Å) if the target-template sequence identity is higher than 20%. Such structures can be subsequently used for further investigation of molecular mechanisms of protein functioning and for the development of modern protein-based biotechnologies.

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.

A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling

It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.

Identification of Specific Effect of Chloride on the Spectral Properties and Structural Stability of Multiple Extracellular Glutamic Acid Mutants of Bacteriorhodopsin

In the present work we combine spectroscopic, DSC and computational approaches to examine the multiple extracellular Glu mutants E204Q/E194Q, E204Q/E194Q/E9Q and E204Q/E194Q/E9Q/E74Q of bacteriorhodopsin by varying solvent ionic strength and composition. Absorption spectroscopy data reveal that the absorption maxima of multiple EC Glu mutants can be tuned by the chloride concentration in the solution. Visible Circular dichroism spectra imply that the specific binding of Cl- can modulate weakened exciton chromophore coupling and reestablish wild type-like bilobe spectral features of the mutants. The DSC data display reappearance of the reversible thermal transition, higher Tm of denaturation and an increase in the enthalpy of unfolding of the mutants in 1 M KCl solutions. Molecular dynamics simulations indicate high affinity binding of Cl- to Arg82 and to Gln204 and Gln194 residues in the mutants. Analysis of the experimental data suggests that simultaneous elimination of the negatively charged side chain of Glu194 and Glu204 is the major cause for mutants' alterations. Specific Cl- binding efficiently coordinates distorted hydrogen bonding interactions of the EC region and reconstitutes the conformation and structure stability of mutated bR in WT-like fashion.

Excited-State Proton Transfer Controls Irreversibility of Photoisomerization in Mononuclear Ruthenium(II) Monoaquo Complexes: A DFT Study

The detailed DFT investigation clears the working mechanism of the irreversible photoisomerization of trans-[Ru(tpy)(pynp)(OH2)](2+) (TA) and cis-[Ru(tpy)(pynp)(OH2)](2+) (CA) complexes. Both TA and CA complexes present two types of low lying triplet states, one resulting from a triplet metal-ligand charge-transfer (TMLCT) and the other from a triplet metal-centered d-d transition (TMC). The vertical excitation of the singlet ground state of the complexes leads to a singlet excited state, which undergoes ultrafast decay to the corresponding TMLCT. For TA, this TMLCT transforms with a low barrier to a TMC state. The dissociative nature of the TMC state leads to easy water removal to produce a five-coordinate intermediate that can isomerize via rotation of a pynp ligand and proceed towards the CA product. For CA, however, during this excitation and intersystem crossing process, an excited-state proton transfer (ESPT) occurs and the resultant TMLCT is very much stabilized with a very strong Ru(II)-OH bond; the high barrier from this TMLCT blocks conversion to a TMC state and thus prevents isomerization from the cis to the trans isomer. This high barrier also prevents the possibility of the isomerization process from TA to CA solely on the adiabatic triplet pathway. Instead, crossing points (XMC-CB, XMC-CA) near the minimum of the triplet metal-centered state of the cis isomer provide nonadiabatic decay channels to the ground-state S0--CA, which completes the photoisomerization pathway from TA to CA.

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.

Protein-Induced Color Shift of Carotenoids in β-Crustacyanin

β-Crustacyanin (β-CR) is a pigment protein responsible for the blue color of lobsters. We show using correlated ab initio calculations how the protein environment tunes the chromophores of β-CR through electrostatic and steric effects.

Synthetic retinal analogues modify the spectral and kinetic characteristics of microbial rhodopsin optogenetic tools

Optogenetic tools have become indispensable in neuroscience to stimulate or inhibit excitable cells by light. Channelrhodopsin-2 (ChR2) variants have been established by mutating the opsin backbone or by mining related algal genomes. As an alternative strategy, we surveyed synthetic retinal analogues combined with microbial rhodopsins for functional and spectral properties, capitalizing on assays in C. elegans, HEK cells and larval Drosophila. Compared with all-trans retinal (ATR), Dimethylamino-retinal (DMAR) shifts the action spectra maxima of ChR2 variants H134R and H134R/T159C from 480 to 520 nm. Moreover, DMAR decelerates the photocycle of ChR2(H134R) and (H134R/T159C), thereby reducing the light intensity required for persistent channel activation. In hyperpolarizing archaerhodopsin-3 and Mac, naphthyl-retinal and thiophene-retinal support activity alike ATR, yet at altered peak wavelengths. Our experiments enable applications of retinal analogues in colour tuning and altering photocycle characteristics of optogenetic tools, thereby increasing the operational light sensitivity of existing cell lines or transgenic animals.