Channelrhodopsins (ChRs) belong to the unique class of light-gated ion channels. The structure of channelrhodopsin-2 from Chlamydomonas reinhardtii (CrChR2) has been resolved, but the mechanistic link between light-induced isomerization of the chromophore retinal and channel gating remains elusive. Replacements of residues C128 and D156 (DC gate) resulted in drastic effects in channel closure. T127 is localized close to the retinal Schiff base and links the DC gate to the Schiff base. The homologous residue in bacteriorhodopsin (T89) has been shown to be crucial for the visible absorption maximum and dark–light adaptation, suggesting an interaction with the retinylidene chromophore, but the replacement had little effect on photocycle kinetics and proton pumping activity. Here, we show that the T127A and T127S variants of CrChR2 leave the visible absorption maximum unaffected. We inferred from hybrid quantum mechanics/molecular mechanics (QM/MM) calculations and resonance Raman spectroscopy that the hydroxylic side chain of T127 is hydrogen-bonded to E123 and the latter is hydrogen-bonded to the retinal Schiff base. The C=N–H vibration of the Schiff base in the T127A variant was 1674 cm−1, the highest among all rhodopsins reported to date. We also found heterogeneity in the Schiff base ground state vibrational properties due to different rotamer conformations of E123. The photoreaction of T127A is characterized by a long-lived P2380 state during which the Schiff base is deprotonated. The conservative replacement of T127S hardly affected the photocycle kinetics. Thus, we inferred that the hydroxyl group at position 127 is part of the proton transfer pathway from D156 to the Schiff base during rise of the P3530 intermediate. This finding provides molecular reasons for the evolutionary conservation of the chemically homologous residues threonine, serine, and cysteine at this position in all channelrhodopsins known so far.

Abstract The origin of the spectral shift from a red- to a green-absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open-chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.

The solid-state photo-chemically induced dynamic nuclear polarization (photo-CIDNP) effect generates non-equilibrium nuclear spin polarization in frozen electron-transfer proteins upon illumination and radical-pair formation. The effect can be observed in various natural photosynthetic reaction center proteins using magic-angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy, and in a flavin-binding light-oxygen-voltage (LOV) domain of the blue-light receptor phototropin. In the latter system, a functionally instrumental cysteine has been mutated to interrupt the natural cysteine-involving photochemistry allowing for an electron transfer from a more distant tryptophan to the excited flavin mononucleotide chromophore. We explored the solid-state photo-CIDNP effect and its mechanisms in phototropin-LOV1-C57S from the green alga Chlamydomonas reinhardtii by using field-cycling solution NMR. We observed the 13C and, to our knowledge, for the first time, 15N photo-CIDNP signals from phototropin-LOV1-C57S. Additionally, the 1H photo-CIDNP signals of residual water in the deuterated buffer of the protein were detected. The relative strengths of the photo-CIDNP effect from the three types of nuclei, 1H, 13C and 15N were measured in dependence of the magnetic field, showing their maximum polarizations at different magnetic fields. Theoretical level crossing analysis demonstrates that anisotropic mechanisms play the dominant role at high magnetic fields.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that are promising candidates for biotechnological applications. Computational studies can contribute to an understanding at a molecular level of their wide spectral tuning and diversity. In this contribution, we benchmark methods to model a 110 nm shift in the UV/Vis absorption spectrum from a red- to a green-absorbing form of the cyanobacteriochrome Slr1393g3. Based on an assessment of semiempirical methods to describe the chromophore geometries of both forms in vacuo, we find that DFTB2+D leads to structures that are the closest to the reference method. The benchmark of the excited state calculations is based on snapshots from quantum mechanics/molecular mechanics molecular dynamics simulations. In our case, the methods RI-ADC(2) and sTD-DFT based on CAM-B3LYP ground state calculations perform the best, whereas no functional can be recommended to simulate the absorption spectra of both forms with time-dependent density functional theory. Furthermore, the difference in absorption for the lowest energy absorption maxima of both forms can already be modelled with optimized structures, but sampling is required to improve the shape of the absorption bands of both forms, in particular for the second band. This benchmark study can guide further computational studies, as it assesses essential components of a protocol to model the spectral tuning of both cyanobacteriochromes and the related phytochromes.

Abstract In this work, the electronic structure and spectroscopic properties of lumiflavin are calculated using various quantum chemical methods. The excitation energies for ten singlet and triplet states as well as the analysis of the electron density difference are assessed using various wave function-based methods and density functionals. The relative order of singlet and triplet excited states is established on the basis of the coupled cluster method CC2. We find that at least seven singlet excited states are required to assign all peaks in the UV/Vis spectrum. In addition, we have studied the solvatochromic effect on the excitation energies and found differential effects except for the first bright excited state. Vibrational frequencies as well as IR, Raman and resonance Raman intensities are simulated and compared to their experimental counterparts. We have assigned peaks, assessed the effect of anharmonicity, and confirmed the previous assignments in case of the most intense transitions. Finally, we have studied the NMR shieldings and established the effect of the solvent polarity. The present study provides data for lumiflavin in the gas phase and in implicit solvent model that can be used as a reference for the protein-embedded flavin simulations and assignment of experimental spectra.

A molecular motor potentially performing a continuous unidirectional rotation is studied by a multidisciplinary approach including organic synthesis, transient spectroscopy and excited state trajectory calculations. A stereogenic center was introduced in the N-alkylated indanylidene–pyrroline Schiff base framework of a previously investigated light-driven molecular switch in order to achieve the unidirectional CC rotary motion typical of Feringa’s motor. Here we report that the specific substitution pattern of the designed chiral molecule must critically determine the unidirectional efficiency of the light-induced rotary motion. More specifically, we find that a stereogenic center containing a methyl group and a hydrogen atom as substituents does not create a differential steric effect large enough to fully direct the motion in either the clockwise or counterclockwise direction especially along the E → Z coordinate. However, due to the documented ultrafast character and electronic circular dichroism activity of the investigated system, we find that it provides the basis for development of a novel generation of rotary motors with a biomimetic framework and operating on a picosecond time scale.