Red Shift in the Absorption Spectrum of Phototropin LOV1 upon the Formation of a Semiquinone Radical: Reconstructing the Orbital Architecture.

Flavin mononucleotide (FMN) is a ubiquitous blue-light pigment due to its ability to drive one- and two-electron transfer reactions. In both light-oxygen-voltage (LOV) domains of phototropin from the green algae Chlamydomonas reinhardtii, FMN is noncovalently bound. In the LOV1 cysteine-to-serine mutant (C57S), light-induced electron transfer from a nearby tryptophan occurs, and a transient spin-correlated radical pair (SCRP) is formed. Within this photocycle, nuclear hyperpolarization is created by the solid-state photochemically induced dynamic nuclear polarization (photo-CIDNP) effect. In a side reaction, a stable protonated semiquinone radical (FMNH·) forms undergoing a significant bathochromic shift of the first electronic transition from 445 to 591 nm. The incorporation of phototropin LOV1-C57S into an amorphous trehalose matrix, stabilizing the radical, allows for application of various magnetic resonance experiments at ambient temperatures, which are combined with quantum-chemical calculations. As a result, the bathochromic shift of the first absorption band is explained by lifting the degeneracy of the molecular orbital energy levels for electrons with alpha and beta spins in FMNH· due to the additional electron.

Electronic Structures of Radical-Pair-Forming Cofactors in a Heliobacterial Reaction Center.

Photosynthetic reaction centers (RCs) are membrane proteins converting photonic excitations into electric gradients. The heliobacterial RCs (HbRCs) are assumed to be the precursors of all known RCs, making them a compelling subject for investigating structural and functional relationships. A comprehensive picture of the electronic structure of the HbRCs is still missing. In this work, the combination of selective isotope labelling of 13C and 15N nuclei and the utilization of photo-CIDNP MAS NMR (photochemically induced dynamic nuclear polarization magic-angle spinning nuclear magnetic resonance) allows for highly enhanced signals from the radical-pair-forming cofactors. The remarkable magnetic-field dependence of the solid-state photo-CIDNP effect allows for observation of positive signals of the electron donor cofactor at 4.7 T, which is interpreted in terms of a dominant contribution of the differential relaxation (DR) mechanism. Conversely, at 9.4 T, the emissive signals mainly originate from the electron acceptor, due to the strong activation of the three-spin mixing (TSM) mechanism. Consequently, we have utilized two-dimensional homonuclear photo-CIDNP MAS NMR at both 4.7 T and 9.4 T. These findings from experimental investigations are corroborated by calculations based on density functional theory (DFT). This allows us to present a comprehensive investigation of the electronic structure of the cofactors involved in electron transfer (ET).

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.

Stabilization of a flavoprotein for solid-state photo-CIDNP MAS NMR at room temperature by embedding in a glassy sugar matrix.

Hyperpolarization via the solid-state photochemically induced dynamic nuclear polarization (photo-CIDNP) effect can be detected in frozen solutions of electron transfer proteins generating a radical-pair upon illumination. The effect has been observed in various natural photosynthetic reaction centers and in light-oxygen-voltage (LOV) sensing domains incorporating a flavin mononucleotide (FMN) as chromophore. In LOV domains, where a highly conserved cysteine is mutated to a flavin to interrupt its natural photochemistry, a radical-pair is generated by electron transfer from a nearby tryptophan to the photoexcited triplet state of FMN. During the photocycle, both the LOV domain and the chromophore are photochemically degraded, e.g., by the formation of singlet oxygen. This limits the time for collection of hyperpolarized nuclear magnetic resonance (NMR) data. We show that embedding of the protein into a trehalose sugar glass matrix stabilizes the protein for 13C solid-state photo-CIDNP NMR experiments which can be conducted at room temperature in a powder sample. Additionally, this preparation allows for incorporation of high amounts of protein further boosting the intensity of the detected signals from FMN and tryptophan at natural abundance. Signal assignment is aided by quantum chemical calculations of absolute shieldings. The underlying mechanism for the surprising absorption-only signal pattern is not yet understood. Comparison to calculated isotropic hyperfine couplings imply that the enhancement is not due to the classical radical-pair mechanism (RPM). Analysis of the anisotropic hyperfine couplings associated with solid-state photo-CIDNP mechanisms also show no simple correlation, suggesting a more complex underlying mechanism.

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.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State.

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

Gas Phase Reactivity of [Mo6X14]2- Dianions (X = Cl - I).

We investigate collision-induced dissociation (CID) of [Mo6X14]2- (X = Cl, Br, I) and the reactivity of fragment ions of these precursors with background gases. Ion mobility measurements and theoretical calculations provide structural information for some of the observed ions. Sequential losses of MoX2 units dominate the dissociation pathways of [Mo6Cl14]2-. Meanwhile, loss of X radicals is the main channel for X = Br and I. Ion mobility measurements and computational investigations indicate minor structural changes in the octahedral Mo6 unit for [Mo6Im]- (m = 6-13) fragments. We observe that mass spectra obtained using CID substantially vary among mass spectrometers: Specifically, ions with molecular formula [Mo6Xm(O2)n]- (X = Br and I) are observed as dominant species produced through reactions with O2 in several mass spectrometers, but also adduct free fragment ions were observed in other instruments, depending on the background conditions. Ion-trap fragmentation combined with theoretical investigations indicates that spontaneous losses of X radicals occur upon binding of O2 to [Mo6Im]- fragments (m ≤ 12). Theoretical investigations indicate that both oxygen atoms are bound to the vacant sites of the Mo6 units. This study opens up a new vista to generate and study a large variety of hexanuclear Mo6Xm(O2)n anions.

Spectroscopic insight into post-synthetic surface modification of porous glass beads as a silica model system.

Applications in catalysis, adsorption and separation require high surface areas as provided by mesoporous materials. Particularly attractive is the class of silica-based mesoporous glasses, which are mechanically and chemically very stable and post-synthetically modifiable allowing specific surface properties to be introduced. One of the catalytically relevant moieties is the sulfonic acid group. To optimize the performance of mesoporous glass systems, analytical methods are required to determine the state of surface modification and its effect on the porosity. To this end, we here propose a specific combination of spectroscopic methods: The porosity during the introduction of thiol functionalities and subsequent oxidation into sulfonic acid groups on the surface of porous micro glass beads is investigated using hyperpolarized 129Xe NMR, revealing that during the two modification steps the textural properties are preserved. The grafting mode as well as the surface coverage are determined using 29Si MAS NMR. The oxidation step is demonstrated to be complete as probed by Raman spectroscopy and 13C MAS NMR. Our combined analysis demonstrates the successful and complete surface modification as well as the maintenance of the free accessibility of the mesopore system.

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.

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 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.

Correction: Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2.

Correction for 'Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2' by Aditya G. Rao et al., Phys. Chem. Chem. Phys., 2021, DOI: 10.1039/d0cp05314g.

Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2.

Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that bind a linear tetrapyrrole as a chromophore. They show photochromicity by having two stable states that can be interconverted by the photoisomerization of the chromophore. These photochemical properties make them an attractive target for biotechnological applications. However, their application is impeded by structural heterogeneity that reduces the yield of the photoconversion. The heterogeneity can originate either from the chromophore structure or the protein environment. Here, we study the origin of the heterogeneity in AnPixJg2, a representative member of the red/green cyanobacteriochrome family, that has a red absorbing parental state and a green absorbing photoproduct state. Using molecular dynamics simulations and umbrella sampling we have identified the protonation state of a conserved histidine residue as a trigger for structural heterogeneity. When the histidine is in a neutral form, the chromophore structure is homogenous, while in a positively charged form, the chromophore is heterogeneous with two different conformations. We have identified a correlation between the protonation of the histidine and the structural heterogeneity of the chromophore by detailed characterization of the interactions in the protein binding site. Our findings reconcile seemingly contradicting spectroscopic studies that attribute the heterogeneity to different sources. Furthermore, we predict that circular dichroism can be used as a diagnostic tool to distinguish different substates.

Hydrogen Bond between a Tyrosine Residue and the C-Ring Propionate Has a Direct Influence on Conformation and Absorption of the Bilin Cofactor in Red/Green Cyanobacteriochromes.

Cyanobacteriochromes (CBCRs) are photoreceptors of the phytochrome superfamily showing remarkable variability in the wavelengths of the first electronic transition-sometimes denoted as Q band-compared to canonical phytochromes. Both classes carry the same cofactor, a bilin, but the molecular basis for the wide variation of their absorption properties is still a matter of debate. The interaction between the cofactor and the surrounding protein moiety has been proposed as a possible tuning factor. Here, we address the impact of hydrogen-bonding interaction between the covalently bound tetrapyrrole cofactor (phycocyanobilin, PCB) and a conserved tyrosine residue (Y302) in the second GAF (cGMP-specific phosphodiesterase, adenylyl cyclases, and FhlA) domain of the red-/green-switching CBCR AnPixJ (AnPixJg2). In the wild type, AnPixJg2 shows absorption maxima of 648 and 543 nm for the dark-adapted (Pr) and photoproduct (Pg) states, respectively. The Y302F mutation leads to the occurrence of an additional absorption band at 687 nm, which is assigned to a new spectroscopically identified sub-state called PIII. Similar spectral changes result upon mutating the Y302F-homologue in another representative red-/green-switching CBCR, Slr1393g3. Molecular dynamics simulations on the dark-adapted state suggest that the removal of the hydrogen bond leads to an additional PCB sub-state differing in its A- and D-ring geometries. The origin of the Q band satellite in the dark-adapted state is discussed.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

The interplay between chromophore and protein determines the extended excited state dynamics in a single-domain phytochrome.

Phytochromes are a diverse family of bilin-binding photoreceptors that regulate a wide range of physiological processes. Their photochemical properties make them attractive for applications in optogenetics and superresolution microscopy. Phytochromes undergo reversible photoconversion triggered by the Z ⇄ E photoisomerization about the double bond in the bilin chromophore. However, it is not fully understood at the molecular level how the protein framework facilitates the complex photoisomerization dynamics. We have studied a single-domain bilin-binding photoreceptor All2699g1 (Nostoc sp. PCC 7120) that exhibits photoconversion between the red light-absorbing (Pr) and far red-absorbing (Pfr) states just like canonical phytochromes. We present the crystal structure and examine the photoisomerization mechanism of the Pr form as well as the formation of the primary photoproduct Lumi-R using time-resolved spectroscopy and hybrid quantum mechanics/molecular mechanics simulations. We show that the unusually long excited state lifetime (broad lifetime distribution centered at ∼300 picoseconds) is due to the interactions between the isomerizing pyrrole ring D and an adjacent conserved Tyr142. The decay kinetics shows a strongly distributed character which is imposed by the nonexponential protein dynamics. Our findings offer a mechanistic insight into how the quantum efficiency of the bilin photoisomerization is tuned by the protein environment, thereby providing a structural framework for engineering bilin-based optical agents for imaging and optogenetics applications.

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.

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.

The Effective Conjugation Length Is Responsible for the Red/Green Spectral Tuning in the Cyanobacteriochrome Slr1393g3.

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.

A QM/MM study of the initial excited state dynamics of green-absorbing proteorhodopsin.

The primary photochemical reaction of the green-absorbing proteorhodopsin is studied by means of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The simulations are based on a homology model derived from the blue-absorbing proteorhodopsin crystal structure. The geometry of retinal and the surrounding sidechains in the protein binding pocket were optimized using the QM/MM method. Starting from this geometry the isomerization was studied with a relaxed scan along the C13[double bond, length as m-dash]C14 dihedral. It revealed an "aborted bicycle pedal" mechanism of isomerization that was originally proposed by Warshel for bovine rhodopsin and bacteriorhodopsin. However, the isomerization involved the concerted rotation about C13[double bond, length as m-dash]C14 and C15[double bond, length as m-dash]N, with the latter being highly twisted but not isomerized. Further, the simulation showed an increased steric interaction between the hydrogen at the C14 of the isomerizing bond and the hydroxyl group at the neighbouring tyrosine 200. In addition, we have simulated a nonadiabatic trajectory which showed the timing of the isomerization. In the first 20 fs upon excitation the order of the conjugated double and single bonds is inverted, consecutively the C13[double bond, length as m-dash]C14 rotation is activated for 200 fs until the S1-S0 transition is detected. However, the isomerization is reverted due to the specific interaction with the tyrosine as observed along the relaxed scan calculation. Our simulations indicate that the retinal - tyrosine 200 interaction plays an important role in the outcome of the photoisomerization.