QM/MM benchmarking of cyanobacteriochrome Slr1393g3 absorption spectra

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

A Massively Parallel Algorithm for the Approximate Calculation of Inverse p-th Roots of Large Sparse Matrices

We present the submatrix method, a highly parallelizable method for the approximate calculation of inverse p-th roots of large sparse symmetric matrices which are required in different scientific applications. We follow the idea of Approximate Computing, allowing imprecision in the final result in order to be able to utilize the sparsity of the input matrix and to allow massively parallel execution. For an n x n matrix, the proposed algorithm allows to distribute the calculations over n nodes with only little communication overhead. The approximate result matrix exhibits the same sparsity pattern as the input matrix, allowing for efficient reuse of allocated data structures. We evaluate the algorithm with respect to the error that it introduces into calculated results, as well as its performance and scalability. We demonstrate that the error is relatively limited for well-conditioned matrices and that results are still valuable for error-resilient applications like preconditioning even for ill-conditioned matrices. We discuss the execution time and scaling of the algorithm on a theoretical level and present a distributed implementation of the algorithm using MPI and OpenMP. We demonstrate the scalability of this implementation by running it on a high-performance compute cluster comprised of 1024 CPU cores, showing a speedup of 665x compared to single-threaded execution

Optical spectrum, perceived color, refractive index, and non-adiabatic dynamics of the photochromic diarylethene CMTE

Photochromism allows for reversible light-induced conversion of a molecular species into a different form with significantly altered optical properties. One promising compound that excels with high fatigue resistance and shows its photochromic functionality both in solution and in molecular solid films is the diarylethene derivative CMTE. Here we present a comprehensive study of its photophysical properties with density-functional theory based methods and benchmark the results against higher-level quantum-chemical approaches and experiments. In addition to static properties such as optical absorption, perceived color, and refractive index, we also investigate reaction dynamics based on non-adiabatic ab initio molecular dynamics. This gives detailed insight into the molecules' ultrafast reaction dynamics and enables us to extract reaction time scales and quantum yields for the observed electrocyclic reaction following photoexcitation

Photocycle of a cyanobacteriochrome: a charge defect on ring C impairs conjugation in chromophore

A large number of novel phytochromes named cyanobacteriochromes (CBCRs) have been recently identified. CBCRs appear to be attractive for further in-depth studies as paradigms for phytochromes because of their related photochemistry, but simpler domain architecture. Elucidating the mechanisms of spectral tuning for the bilin chromophore down to the molecular/atomic level is a prerequisite to design fine-tuned photoswitches for optogenetics. Several explanations for the blue shift during photoproduct formation associated with the red/green CBCRs represented by Slr1393g3 have been developed. There are, however, only sparse mechanistic data concerning the factors controlling stepwise absorbance changes along the reaction pathways from the dark state to the photoproduct and vice versa in this subfamily. Conventional cryotrapping of photocycle intermediates of phytochromes has proven experimentally intractable for solid-state NMR spectroscopy within the probe. Here, we have developed a simple method to circumvent this hindrance by incorporating proteins into trehalose glasses which allows four photocycle intermediates of Slr1393g3 to be isolated for NMR use. In addition to identifying the chemical shifts and chemical shift anisotropy principal values of selective chromophore carbons in various photocycle states, we generated QM/MM models of the dark state and photoproduct as well as of the primary intermediate of the backward-reaction. We find the motion of all three methine bridges in both reaction directions but in different orders. These molecular events channel light excitation to drive distinguishable transformation processes. Our work also suggests that polaronic self-trapping of a conjugation defect by displacement of the counterion during the photocycle would play a role in tuning the spectral properties of both the dark state and photoproduct

Unraveling the electrochemical and spectroscopic properties of neutral and negatively charged perylene tetraethylesters

A detailed investigation of the energy levels of perylene-3,4,9,10-tetracarboxylic tetraethylester as a representative compound for the whole family of perylene esters was performed. It was revealed via electrochemical measurements that one oxidation and two reductions take place. The bandgaps determined via the electrochemical approach are in good agreement with the optical bandgap obtained from the absorption spectra via a Tauc plot. In addition, absorption spectra in dependence of the electrochemical potential were the basis for extensive quantum-chemical calculations of the neutral, monoanionic, and dianionic molecules. For this purpose, calculations based on density functional theory were compared with post-Hartree–Fock methods and the CAM-B3LYP functional proved to be the most reliable choice for the calculation of absorption spectra. Furthermore, spectral features found experimentally could be reproduced with vibronic calculations and allowed to understand their origins. In particular, the two lowest energy absorption bands of the anion are not caused by absorption of two distinct electronic states, which might have been expected from vertical excitation calculations, but both states exhibit a strong vibronic progression resulting in contributions to both bands

Structural elements regulating the photochromicity in a cyanobacteriochrome

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

Light- and pH-dependent structural changes in cyanobacteriochrome AnPixJg2

Cyanobacteriochromes (CBCRs) are phytochrome-related photosensory proteins that play an essential role in regulating phototaxis, chromatic acclimation, and cell aggregation in cyanobacteria. Here, we apply solid-state NMR spectroscopy to the red/green GAF2 domain of the CBCR AnPixJ assembled in vitro with a uniformly 13C- and 15N-labeled bilin chromophore, tracking changes in electronic structure, geometry, and structural heterogeneity of the chromophore as well as intimate contacts between the chromophore and protein residues in the photocycle. Our data confirm that the bilin ring D is strongly twisted with respect to the B–C plane in both dark and photoproduct states. We also identify a greater structural heterogeneity of the bilin chromophore in the photoproduct than in the dark state. In addition, the binding pocket is more hydrated in the photoproduct. Observation of interfacial 1H contacts of the photoproduct chromophore, together with quantum mechanics/molecular mechanics (QM/MM)-based structural models for this photoproduct, clearly suggests the presence of a biprotonated (cationic) imidazolium side-chain for a conserved histidine residue (322) at a distance of ~2.7 Å, generalizing the recent theoretical findings that explicitly link the structural heterogeneity of the dark-state chromophore to the protonation of this specific residue. Moreover, we examine pH effects on this in vitro assembled holoprotein, showing a substantially altered electronic structure and protonation of the photoproduct chromophore even with a small pH drop from 7.8 to 7.2. Our studies provide further information regarding the light- and pH-induced changes of the chromophore and the rearrangements of the hydrogen-bonding and electrostatic interaction network around it. Possible correlations between structural heterogeneity of the chromophore, protonation of the histidine residue nearby, and hydration of the pocket in both photostates are discussed

Impact of Chlorine on the Internal Transition Rates and Excited States of the Thermally Delayed Activated Fluorescence Molecule 3CzClIPN

We analyze internal transition rates and the singlet-triplet energy gap of the thermally activated delayed fluorescence (TADF) molecule 3CzClIPN, which recently was introduced as an efficient photocatalyst. Distribution and origin of the non-monoexponential decays, which are commonly observed in TADF films, are revealed by analysis of transient fluorescence with an inverse Laplace transform. A numerically robust global rate fit routine, which extracts all relevant TADF parameters by modeling the complete set of data, is introduced. To compare and verify the results, all methods are also applied to the well-known 4CzIPN. The influence of the molecular matrix is discussed by embedding low concentrations of TADF molecules in polystyrene films. Finally, quantum chemical calculations are compared to the experimental results to demonstrate that the chlorine atom increases the charge transfer character of the relevant states, resulting in a reduction of the singlet-triplet energy gap

Ultrafast Electronic Energy Transfer in an orthogonal molecular dyad

The St Andrews group acknowledges support from the European Research Council (grant number 321305) and the Engineering and Physical Sciences Research Council (grant EP/L017008/1). I.D.W.S. also acknowledges support from a Royal Society Wolfson Research Merit Award.Understanding electronic energy transfer (EET) is an important ingredient in the development of artificial photosynthetic systems and photovoltaic technologies. Although EET is at the heart of these applications and crucially influences their light-harvesting efficiency, the nature of EET over short distances for covalently bound donor and acceptor units is often not well understood. Here we investigate EET in an orthogonal molecular dyad (BODT4) in which simple models fail to explain the very origin of EET. Based on nonadiabatic ab initio molecular dynamics calculations and fluorescence depolarization experiments we gain detailed microscopic insights into the ultrafast electro-vibrational dynamics following photoexcitation. Our analysis offers molecular-level insights into these processes and reveals that it takes place on timescales ≲ 100 fs and occurs through an intermediate charge-transfer state.PostprintPeer reviewe

Structure and function of the metagenomic plastic-degrading polyester hydrolase PHL7 bound to its product

The recently discovered metagenomic-derived polyester hydrolase PHL7 is able to efficiently degrade amorphous polyethylene terephthalate (PET) in post-consumer plastic waste. We present the cocrystal structure of this hydrolase with its hydrolysis product terephthalic acid and elucidate the influence of 17 single mutations on the PET-hydrolytic activity and thermal stability of PHL7. The substrate-binding mode of terephthalic acid is similar to that of the thermophilic polyester hydrolase LCC and deviates from the mesophilic IsPETase. The subsite I modifications L93F and Q95Y, derived from LCC, increased the thermal stability, while exchange of H185S, derived from IsPETase, reduced the stability of PHL7. The subsite II residue H130 is suggested to represent an adaptation for high thermal stability, whereas L210 emerged as the main contributor to the observed high PET-hydrolytic activity. Variant L210T showed significantly higher activity, achieving a degradation rate of 20 µm h−1 with amorphous PET films