Color tuning in binding pocket models of the chlamydomonas-type channelrhodopsins

Structures of channelrhodopsin paralogs in peptidiscs explain their contrasting K+ and Na+ selectivities

Kalium channelrhodopsin 1 from Hyphochytrium catenoides (HcKCR1) is a light-gated channel used for optogenetic silencing of mammalian neurons. It selects K+ over Na+ in the absence of the canonical tetrameric K+ selectivity filter found universally in voltage- and ligand-gated channels. The genome of H. catenoides also encodes a highly homologous cation channelrhodopsin (HcCCR), a Na+ channel with >100-fold larger Na+ to K+ permeability ratio. Here, we use cryo-electron microscopy to determine atomic structures of these two channels embedded in peptidiscs to elucidate structural foundations of their dramatically different cation selectivity. Together with structure-guided mutagenesis, we show that K+ versus Na+ selectivity is determined at two distinct sites on the putative ion conduction pathway: in a patch of critical residues in the intracellular segment (Leu69/Phe69, Ile73/Ser73 and Asp116) and within a cluster of aromatic residues in the extracellular segment (primarily, Trp102 and Tyr222). The two filters are on the opposite sides of the photoactive site involved in channel gating.

Structural Factors Determining the Absorption Spectrum of Channelrhodopsins: A Case Study of the Chimera C1C2

Channelrhodopsins are photosensitive proteins that trigger flagella motion in single-cell algae and have been successfully utilized in optogenetic applications. In optogenetics, light is used to activate neural cells in living organisms, which can be achieved by exploiting the ion channel signaling of channelrhodopsins. Tailoring channelrhodopsins for such applications includes the tuning of the absorption maximum. In order to establish rational design and to obtain a desired spectral shift, a basic understanding of the absorption spectrum is required. We have studied the chimera C1C2 as a representative of this protein family and the first member with an available crystal structure. For this purpose, we sampled the conformations of C1C2 using quantum mechanical/molecular mechanical molecular dynamics and subjected the resulting snapshots of the trajectory to excitation energy calculations using ADC(2) and simplified time-dependent density functional theory. In contrast to previous reports, we found that different hydrogen-bonding networks-involving the retinal protonated Schiff base, the putative counterions E162 and D292, and water molecules-had only a small impact on the absorption spectrum. However, in the case of deprotonated E162, increasing the distance to the Schiff base hydrogen-bonding partner led to a systematic blue shift. The β-ionone ring rotation was identified as another important contributor. Yet the most important factors were found to be the bond length alternation and bond order alternation that were linearly correlated to the absorption maximum by up to 62 and 82%, respectively. We ascribe this novel insight into the structural basis of the absorption spectrum to our enhanced protein setup that includes membrane embedding as well as long and extensive sampling.

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.

Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes

The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.

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

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

Active site structure and absorption spectrum of channelrhodopsin-2 wild-type and C128T mutant

We show by extensive ground state and absorption spectra simulations that the channelrhodopsin-2 active site samples three different hydrogen-bonding patterns.

In spite of considerable interest, the active site of channelrhodopsin still lacks a detailed atomistic description, the understanding of which could strongly enhance the development of novel optogenetics tools. We present a computational study combining different state-of-the-art techniques, including hybrid quantum mechanics/molecular mechanics schemes and high-level quantum chemical methods, to properly describe the hydrogen-bonding pattern between the retinal chromophore and its counterions in channelrhodopsin-2 Wild-Type and C128T mutant. Especially, we show by extensive ground state dynamics that the active site, containing a glutamic acid (E123) and a water molecule, is highly dynamic, sampling three different hydrogen-bonding patterns. This results in a broad absorption spectrum that is representative of the different structural motifs found. A comparison with bacteriorhodopsin, characterized by a pentagonal hydrogen-bonded active site structure, elucidates their different absorption properties.

A biomimetic molecular switch at work: coupling photoisomerization dynamics to peptide structural rearrangement

A combined experimental and computational study of a peptide-linked retinal-like molecular switch shows the effects on photoreactivity and the α-helix structure.

Early formation of the ion-conducting pore in channelrhodopsin-2

Channelrhodopsins (ChRs) are light-gated ion channels that are widely used in optogenetics. They allow precise control of neuronal activity with light, but a detailed understanding of how the channel is gated and the ions are conducted is still lacking. The recent determination of the X-ray structural model in the closed state marks an important milestone. Herein the open state structure is presented and the early formation of the ion conducting pore is elucidated in atomic detail using time-resolved FTIR spectroscopy. Photo-isomerization of the retinal-chromophore causes a downward movement of the highly conserved E90, which opens the pore. Molecular dynamic (MD) simulations show that water molecules invade through this opened pore, Helix 2 tilts and the channel fully opens within ms. Since E90 is a highly conserved residue, the proposed E90-Helix2-tilt (EHT) model might describe a general activation mechanism and provides a new avenue for further mechanistic studies and engineering.

Role of a helix B lysine residue in the photoactive site in channelrhodopsins

In most studied microbial rhodopsins two conserved carboxylic acid residues (the homologs of Asp-85 and Asp-212 in bacteriorhodopsin) and an arginine residue (the homolog of Arg-82) form a complex counterion to the protonated retinylidene Schiff base, and neutralization of the negatively charged carboxylates causes red shifts of the absorption maximum. In contrast, the corresponding neutralizing mutations in some relatively low-efficiency channelrhodopsins (ChRs) result in blue shifts. These ChRs do not contain a lysine residue in the second helix, conserved in higher efficiency ChRs (Lys-132 in the crystallized ChR chimera). By action spectroscopy of photoinduced channel currents in HEK293 cells and absorption spectroscopy of detergent-purified pigments, we found that in tested ChRs the Lys-132 homolog controls the direction of spectral shifts in the mutants of the photoactive site carboxylic acid residues. Analysis of double mutants shows that red spectral shifts occur when this Lys is present, whether naturally or by mutagenesis, and blue shifts occur when it is replaced with a neutral residue. A neutralizing mutation of the Lys-132 homolog alone caused a red spectral shift in high-efficiency ChRs, whereas its introduction into low-efficiency ChR1 from Chlamydomonas augustae (CaChR1) caused a blue shift. Taking into account that the effective charge of the carboxylic acid residues is a key factor in microbial rhodopsin spectral tuning, these findings suggest that the Lys-132 homolog modulates their pKa values. On the other hand, mutation of the Arg-82 homolog that fulfills this role in bacteriorhodopsin caused minimal spectral changes in the tested ChRs. Titration revealed that the pKa of the Asp-85 homolog in CaChR1 lies in the alkaline region unlike in most studied microbial rhodopsins, but is substantially decreased by introduction of a Lys-132 homolog or neutralizing mutation of the Asp-212 homolog. In the three ChRs tested the Lys-132 homolog also alters channel current kinetics.

Effect of protein environment within cytochrome P450cam evaluated using a polarizable-embedding QM/MM method

Metalloenzymes accommodate cofactors and substrates in their active sites, thereby exerting powerful catalytic effects. Understanding the key elements of the mechanism via which such binding is accomplished using a number of atoms in a protein is a fundamental challenge. To address this issue computationally, here we used mechanical-embedding (ME), electronic-embedding (EE), and polarizable-embedding (PE) hybrid quantum mechanics and molecular mechanics (QM/MM) methods and performed an energy decomposition analysis (EDA) of the nonbonding protein environmental effect in the "compound I" intermediate state of cytochrome P450cam. The B3LYP and AMBER99/QP302 methods were used to deal with the QM and MM subsystems, respectively, and the nonbonding interaction energy between these subsystems was decomposed into electrostatic, van der Waals, and polarization contributions. The PE-QM/MM calculation was performed using polarizable force fields that were capable of describing induced dipoles within the MM subsystem, which arose in response to the electric field generated by QM electron density, QM nuclei, and MM point charges. The present QM/MM EDA revealed that the electrostatic term constituted the largest stabilizing interaction between the QM and MM subsystems. When proper adjustment was made for the point charges of the MM atoms located at the QM-MM boundary, EE-QM/MM and PE-QM/MM calculations yielded similar QM electron density distributions, indicating that the MM polarization effect does not have a significant influence on the extent of QM polarization in this particular enzyme system.

Gas-Phase Retinal Spectroscopy: Temperature Effects Are But a Mirage

We employ state-of-the-art first-principle approaches to investigate whether temperature effects are responsible for the unusually broad and flat spectrum of protonated Schiff base retinal observed in photodissociation spectroscopy, as has recently been proposed. We first carefully calibrate how to construct a realistic geometrical model of retinal and show that the exchange-correlation M06-2X functional yields an accurate description while the commonly used complete active space self-consistent field method (CASSCF) is not adequate. Using modern multiconfigurational perturbative methods (NEVPT2) to compute the excitations, we then demonstrate that conformations with different orientations of the β-ionone ring are characterized by similar excitations. Moreover, other degrees of freedom identified as active in room-temperature molecular dynamics simulations do not yield the shift required to explain the anomalous spectral shape. Our findings indicate that photodissociation experiments are not representative of the optical spectrum of retinal in the gas phase and call for further experimental characterization of the dissociation spectra.