Phytochromes are biological photoreceptors found in all kingdoms of life. Numerous physicochemical and spectroscopic studies of phytochromes have been carried out for many decades, both experimentally and computationally, with the main focus on the photoconversion mechanism involving a tetrapyrrole chromophore. In this computational work, we concentrate on the long-scale dynamic motion of the photosensory domain of Deinococcus radiodurans by means of classical all-atom molecular dynamics (MD) simulations. Conventional and accelerated MD methods in combination with two different force fields, CHARMM27 and AMBER ff14SB, are tested in long atomistic simulations to confront the dynamics of monomer and dimer forms. These calculations highlight dissimilar equilibrium conformations in aqueous solutions and, in turn, different large-scale dynamic behaviors of the monomer form vs the dimer form. While the phytochrome in a monomer form tends to close the cavity entailed between the GAF and PHY domains, the opposite trend is predicted for the phytochrome dimer, which opens up as a consequence of the formation of strong salt bridges between the PHY domains of two molecules in water.Phytochromes are biological photoreceptors found in all kingdoms of life. Numerous physicochemical and spectroscopic studies of phytochromes have been carried out for many decades, both experimentally and computationally, with the main focus on the photoconversion mechanism involving a tetrapyrrole chromophore. In this computational work, we concentrate on the long-scale dynamic motion of the photosensory domain of Deinococcus radiodurans by means of classical all-atom molecular dynamics (MD) simulations. Conventional and accelerated MD methods in combination with two different force fields, CHARMM27 and AMBER ff14SB, are tested in long atomistic simulations to confront the dynamics of monomer and dimer forms. These calculations highlight dissimilar equilibrium conformations in aqueous solutions and, in turn, different large-scale dynamic behaviors of the monomer form vs the dimer form. While the phytochrome in a monomer form tends to close the cavity entailed between the GAF and PHY domains, the opposite trend is predicted for the phytochrome dimer, which opens up as a consequence of the formation of strong salt bridges between the PHY domains of two molecules in water.

ABC importers are membrane proteins responsible for the transport of nutrients into the cells of prokaryotes. Although the structures of ABC importers vary, all contain four conserved domains: two nucleotide-binding domains (NBDs), which bind and hydrolyze ATP, and two transmembrane domains (TMDs), which help translocate the substrate. ABC importers are also dependent on an additional protein component, a high-affinity substrate-binding protein (SBP) that specifically binds the target ligand for delivery to the appropriate ABC transporter. AbnE is a SBP belonging to the ABC importer for arabino-oligosaccharides in the Gram-positive thermophilic bacterium Geobacillus stearothermophilus. Using isothermal titration calorimetry (ITC), purified AbnE was shown to bind medium-sized arabino-oligosaccharides, in the range of arabino-triose (A3) to arabino-octaose (A8), all with Kd values in the nanomolar range. We describe herein the 3D structure of AbnE in its closed conformation in complex with a wide range of arabino-oligosaccharide substrates (A2-A8). These structures provide the basis for the detailed structural analysis of the AbnE-sugar complexes, and together with complementary quantum chemical calculations, site-specific mutagenesis, and isothermal titration calorimetry (ITC) experiments, provide detailed insights into the AbnE-substrate interactions involved. Small-angle X-ray scattering (SAXS) experiments and normal mode analysis (NMA) are used to study the conformational changes of AbnE, and these data, taken together, suggest clues regarding its binding mode to the full ABC importer.

Channelrhodopsins (ChR) are light-gated ion-channels heavily used in optogenetics. Upon light excitation an ultrafast all-trans to 13-cis isomerization of the retinal chromophore takes place. It is still uncertain by what means this reaction leads to further protein changes and channel conductivity. Channelrhodopsin-1 in Chlamydomonas augustae exhibits a 100 fs photoisomerization and a protonated counterion complex. By polarization resolved ultrafast spectroscopy in the mid-IR we show that the initial reaction of the retinal is accompanied by changes in the protein backbone and ultrafast protonation changes at the counterion complex comprising Asp299 and Glu169. In combination with homology modelling and quantum mechanics/molecular mechanics (QM/MM) geometry optimization we assign the protonation dynamics to ultrafast deprotonation of Glu169, and transient protonation of the Glu169 backbone, followed by a proton transfer from the backbone to the carboxylate group of Asp299 on a timescale of tens of picoseconds. The second proton transfer is not related to retinal dynamics and reflects pure protein changes in the first photoproduct. We assume these protein dynamics to be the first steps in a cascade of protein-wide changes resulting in channel conductivity.
 

Phytochromes are biological photoswitches that interconvert between two parent states (Pr and Pfr). The transformation is initiated by photoisomerization of the tetrapyrrole chromophore, followed by a sequence of chromophore and protein structural changes. In the last step, a phytochrome-specific peptide segment (tongue) undergoes a secondary structure change, which in prokaryotic phytochromes is associated with the (de)activation of the output module. The focus of this work is the Pfr-to-Pr photoconversion of the bathy bacteriophytochrome Agp2 in which Pfr is the thermodynamically stable state. Using spectroscopic techniques, we studied the structural and functional consequences of substituting Arg211, Tyr165, His278, and Phe192 close to the biliverdin (BV) chromophore. In Pfr, substitutions of these residues do not affect the BV structure. The characteristic Pfr properties of bathy phytochromes, including the protonated propionic side chain of ring C (propC) of BV, are preserved. However, replacing Arg211 or Tyr165 blocks the photoconversion in the Meta-F state, prior to the secondary structure transition of the tongue and without deprotonation of propC. The Meta-F state of these variants displays low photochemical activity, but electronic excitation causes ultrafast alterations of the hydrogen bond network surrounding the chromophore. In all variants studied here, thermal back conversion from the photoproducts to Pfr is decelerated but substitution of His278 or Phe192 is not critical for the Pfr-to-Pr photoconversion. These variants do not impair deprotonation of propC or the α-helix/β-sheet transformation of the tongue during the Meta-F-to-Pr decay. Thus, we conclude that propC deprotonation is essential for restructuring of the tongue.Phytochromes are biological photoswitches that interconvert between two parent states (Pr and Pfr). The transformation is initiated by photoisomerization of the tetrapyrrole chromophore, followed by a sequence of chromophore and protein structural changes. In the last step, a phytochrome-specific peptide segment (tongue) undergoes a secondary structure change, which in prokaryotic phytochromes is associated with the (de)activation of the output module. The focus of this work is the Pfr-to-Pr photoconversion of the bathy bacteriophytochrome Agp2 in which Pfr is the thermodynamically stable state. Using spectroscopic techniques, we studied the structural and functional consequences of substituting Arg211, Tyr165, His278, and Phe192 close to the biliverdin (BV) chromophore. In Pfr, substitutions of these residues do not affect the BV structure. The characteristic Pfr properties of bathy phytochromes, including the protonated propionic side chain of ring C (propC) of BV, are preserved. However, replacing Arg211 or Tyr165 blocks the photoconversion in the Meta-F state, prior to the secondary structure transition of the tongue and without deprotonation of propC. The Meta-F state of these variants displays low photochemical activity, but electronic excitation causes ultrafast alterations of the hydrogen bond network surrounding the chromophore. In all variants studied here, thermal back conversion from the photoproducts to Pfr is decelerated but substitution of His278 or Phe192 is not critical for the Pfr-to-Pr photoconversion. These variants do not impair deprotonation of propC or the α-helix/β-sheet transformation of the tongue during the Meta-F-to-Pr decay. Thus, we conclude that propC deprotonation is essential for restructuring of the tongue.

Phytochromes and related photoreceptors distinguish themselves for their long-wavelength absorption and large spectral shift between parental state and photoproduct. Both features are not well understood, partly due to lack of high-resolution structural data and insufficient support from quantum-chemical calculations. The red–green switching cyanobacteriochrome Slr1393g3 shows an absorption shift larger than 110 nm. Both parental state and photoproduct could be crystallized with high resolution, together with a “hybrid” form carrying the chromophore in parental state geometry, whereas the protein remained in the photoproduct conformation. The crystal structures reveal how chromophore and protein mutually regulate their conformational changes, yielding the observed spectral shift. Quantum-chemical calculations, based on these structures, provide a deeper understanding of the spectral tuning mechanisms.The three-dimensional (3D) crystal structures of the GAF3 domain of cyanobacteriochrome Slr1393 (Synechocystis PCC6803) carrying a phycocyanobilin chromophore could be solved in both 15-Z dark-adapted state, Pr, λmax = 649 nm, and 15-E photoproduct, Pg, λmax = 536 nm (resolution, 1.6 and 1.86 Å, respectively). The structural data allowed identifying the large spectral shift of the Pr-to-Pg conversion as resulting from an out-of-plane rotation of the chromophore’s peripheral rings and an outward movement of a short helix formed from a formerly unstructured loop. In addition, a third structure (2.1-Å resolution) starting from the photoproduct crystals allowed identification of elements that regulate the absorption maxima. In this peculiar form, generated during X-ray exposition, protein and chromophore conformation still resemble the photoproduct state, except for the D-ring already in 15-Z configuration and tilted out of plane akin the dark state. Due to its formation from the photoproduct, it might be considered an early conformational change initiating the parental state-recovering photocycle. The high quality and the distinct features of the three forms allowed for applying quantum-chemical calculations in the framework of multiscale modeling to rationalize the absorption maxima changes. A systematic analysis of the PCB chromophore in the presence and absence of the protein environment showed that the direct electrostatic effect is negligible on the spectral tuning. However, the protein forces the outer pyrrole rings of the chromophore to deviate from coplanarity, which is identified as the dominating factor for the color regulation.

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 ¹³C and, to our knowledge, for the first time, ¹⁵N photo-CIDNP signals from phototropin-LOV1-C57S. Additionally, the ¹H 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, ¹H, ¹³C and ¹⁵N 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.

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.

Color vision is based on the ability of different opsin proteins to tune the absorption band of their chromophore, the retinal protonated Schiff base (RPSB). Two main mechanisms proposed for this tunability are geometric and electrostatic. Here we probe the latter effect experimentally and by a quantum chemical calculation of the absorption by an isolated complex containing the retinal chromophore and molecules with a strong dipole moment. Betaine complexation causes an anomalously large blue shift. The shift provides direct evidence that the strong charge-transfer character of the electronic transition is the cause of the opsin shift, and shows that the electric field of the counterion is responsible for the color tuning, which allows absorption of light in the blue region of the visible spectrum by opsin proteins.

The first event of the channelrhodopsin-2 (ChR2) photocycle, i.e. trans-to-cis photoisomerization, is studied by means of quantum mechanics/molecular mechanics, taking into account the flexible retinal environment in the ground state. By treating the chromophore at the ab initio multiconfigurational level of theory, we can rationalize the experimental findings based on pump–probe spectroscopy, explaining the different and more complex scenario found for ChR2 in comparison to other rhodopsins. In particular, we find that depending on the hydrogen bonding pattern, different excited states are involved, hence making it possible to suggest one pattern as the most productive. Moreover, after photoisomerization the structure of the first photocycle intermediate, P⁵⁰⁰₁, is characterized by simulating the infrared spectrum and compared to available experimental data. This was obtained by extensive molecular dynamics, where the chromophore is described by a semi-empirical method based on density functional theory. The results clearly identify which counterion is responsible for accepting the proton from the retinal Schiff base: the side chain of the glutamic acid E123.

ABSTRACT The performance of the extended multi-state (XMS)-complete active space second-order perturbation theory (CASPT2) method has been assessed for the benchmark of a truncated retinal model, the penta-2,4-dieniminium cation (PSB3). This benchmark presents a challenge for multireference electronic structure methods because the wave function character is changing considerably. The assessment comprises ground and excited state pathways of the isomerisation, including transition states and conical intersection (CI) points. It also includes circular paths centred around different CIs, and 2D potential energy scans located in the branching planes. In this work, we compare the performance of the previous formulations of CASPT2, the single-state and the multi-state, with the recently developed XMS-CASPT2. Besides, we have also tested two variants of internal contraction in XMS-CASPT2, namely, the single-state single reference (SS-SR) and multi-state multireference (MS-MR) schemes. In our study, we find that XMS-CASPT2 corrects the artefacts and discontinuities present in other CASPT2 variants. The investigation of a circular loop and 2D potential energy surfaces around the surface crossing point shows that XMS-CASPT2 exhibits a smooth topology at the CI with the correct degeneracy. It also agrees better with the reference method MRCISD+Q in regions of the potential energy surfaces further away from CIs. Another observation is the close agreement between the results from the SS-SR contraction scheme and the more demanding MS-MR scheme.

Photochemically induced dynamic nuclear polarization (photo-CIDNP) has been observed in the homodimeric, type-1 photochemical reaction centers (RCs) of the acidobacterium, Chloracidobacterium (Cab.) thermophilum, by ¹⁵N magic-angle spinning (MAS) solid-state NMR under continuous white-light illumination. Three light-induced emissive (negative) signals are detected. In the RCs of Cab. thermophilum, three types of (bacterio)chlorophylls have previously been identified: bacteriochlorophyll a (BChl a), chlorophyll a (Chl a), and Zn-bacteriochlorophyll a′ (Zn-BChl a′) (Tsukatani et al. in J Biol Chem 287:5720–5732, 2012). Based upon experimental and quantum chemical ¹⁵N NMR data, we assign the observed signals to a Chl a cofactor. We exclude Zn-BChl because of its measured spectroscopic properties. We conclude that Chl a is the primary electron acceptor, which implies that the primary donor is most likely Zn-BChl a′. Chl a and 8¹-OH Chl a have been shown to be the primary electron acceptors in green sulfur bacteria and heliobacteria, respectively, and thus a Chl a molecule serves this role in all known homodimeric type-1 RCs.