Molecular dynamics simulations of the chromophore binding site of Deinococcus radiodurans bacteriophytochrome using new force field parameters for the phytochromobilin chromophore

Restructuring the interfacial active sites to generalize the volcano curves for platinum-cobalt synergistic catalysis

Computationally derived volcano curve has become the gold standard in catalysis, whose practical application usually relies on empirical interpretations of composition or size effects by the identical active site assumption. Here, we present a proof-of-concept study on disclosing both the support- and adsorbate-induced restructuring of Pt-Co bimetallic catalysts, and the related interplays among different interfacial sites to propose the synergy-dependent volcano curves. Multiple characterizations, isotopic kinetic investigations, and multiscale simulations unravel that the progressive incorporation of Co into Pt catalysts, driven by strong Pt-C bonding (metal-support interfaces) and Co-O bonding (metal-adsorbate interfaces), initiates the formation of Pt-rich alloys accompanied by isolated Co species, then Co segregation to epitaxial CoOx overlayers and adjacent Co3O4 clusters, and ultimately structural collapse into amorphous alloys. Accordingly, three distinct synergies, involving lattice oxygen redox from Pt-Co alloy/Co3O4 clusters, dual-active sites engineering via Pt-rich alloy/CoOx overlayer, and electron coupling within exposed alloy, are identified and quantified for CO oxidation (gas-phase), ammonia borane hydrolysis (liquid-phase), and hydrogen evolution reaction (electrocatalysis), respectively. The resultant synergy-dependent volcano curves represent an advancement over traditional composition-/size-dependent ones, serving as a bridge between theoretical models and experimental observations in bimetallic catalysis.

Photoisomerization of phytochrome chromophore models: an XMS-CASPT2 study

Phytochromes are a superfamily of photoreceptors that harbor linear tetrapyrroles as chromophores. Upon light illumination, the linear tetrapyrrole chromophore undergoes a double bond isomerization which starts a photocycle. In this work, we studied the photoisomerization of chromophore models designed based on the C- and D-rings of the phycocyanobilin (PCB) chromophore. In total, five different models with varying substitutions were investigated. Firstly, the vertical excitation energies were benchmarked using different computational methods to establish the relative order of the excited states. Based on these calculations, we computed the photoisomerization profiles using the extended multi-state (XMS) version of the CASPT2 method. The profiles were obtained for both the clockwise and counterclockwise rotations of the C₁₅C₁₆ bond in the Z and E isomers using a linear interpolation of internal coordinates between the Franck-Condon and MECI geometries. In the minimal chromophore model that lacks the substitutions at the pyrrole rings, the isomerization involves both C₁₄-C₁₅ and C₁₅C₁₆ bonds of the methine bridge between the C- and D-rings, resembling the hula-twist motion. The MECIs are characterized by a partial charge transfer between the two pyrrole rings pointing towards a twisted intramolecular charge transfer. Systematic introduction of substituents leads to an increase in the steric repulsion between the two pyrrole rings causing a pretwist of the dihedral around the C₁₅C₁₆ bond, which creates a preference for the counterclockwise isomerization. An introduction of the carbonyl group at the D-ring increases the extent of charge transfer which changes the isomerization mechanism from hula-twist to one-bond flip.

Transient electronic and vibrational signatures during reversible photoswitching of a cyanobacteriochrome photoreceptor

Cyanobacteriochromes (CBCRs) are an emerging class of photoreceptors that are distant relatives of the phytochromes family. Unlike phytochromes, CBCRs have gained popularity in optogenetics due to their highly diverse spectral properties spanning the UV to near-IR region and only needing a single compact binding domain. AnPixJg2 is a CBCR that can reversibly photoswitch between its red-absorbing (^15ZP_r) and green-absorbing (^15EP_g) forms of the phycocyanobilin (PCB) cofactor. To reveal primary events of photoconversion, we implemented femtosecond transient absorption spectroscopy with a homemade LED box and a miniature peristaltic pump flow cell to track transient electronic responses of the photoexcited AnPixJg2 on molecular time scales. The 525 nm laser-induced P_g-to-P_r reverse conversion exhibits a ~3 ps excited-state lifetime before reaching the conical intersection (CI) and undergoing further relaxation on the 30 ps time scale to generate a long-lived Lumi-G ground state intermediate en route to P_r. The 650 nm laser-induced P_r-to-P_g forward conversion is less efficient than reverse conversion, showing a longer-lived excited state which requires two steps with ~13 and 217 ps time constants to enter the CI region. Furthermore, using a tunable ps Raman pump with broadband Raman probe on both the Stokes and anti-Stokes sides, we collected the pre-resonance ground-state femtosecond stimulated Raman spectroscopy (GS-FSRS) data with mode assignments aided by quantum calculations. Key vibrational marker bands at ~850, 1050, 1615, and 1649 cm^-1 of the P_r conformer exhibit a notable blueshift to those of the P_g conformer inside AnPixJg2, reflecting the PCB chromophore terminal D (major) and A (minor) ring twist along the primary photoswitching reaction coordinate. This integrated ultrafast spectroscopy and computational platform has the potential to elucidate photochemistry and photophysics of more CBCRs and photoactive proteins in general, providing the highly desirable mechanistic insights to facilitate the rational design of functional molecular sensors and devices.

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

Dynamic Properties of the Photosensory Domain of Deinococcus radiodurans Bacteriophytochrome

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.

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.

Structural heterogeneity in a parent ground-state structure of AnPixJg2 revealed by theory and spectroscopy

We investigated the red absorbing, dark stable state (Pr state) of the second GAF domain of the cyanobacteriochrome AnPixJ (AnPixJg2) by a molecular dynamics simulation of 1 μs duration. Our results reveal two distinct conformational isoforms of the chromophore, from which only one was known from crystallographic experiments. The interconversion between both isoforms is accompanied by alterations in the hydrogen bond pattern between the chromophore and the protein and the solvation structure of the chromophore binding pocket. The existence of sub-states in the Pr form of AnPixJg2 is supported by the results from experimental ¹³C MAS NMR spectroscopy. Our finding is consistent with the observation of structural heterogeneity in other cyanobacteriochromes and phytochromes.

Light-induced structural changes in a monomeric bacteriophytochrome

Phytochromes sense red light in plants and various microorganism. Light absorption causes structural changes within the protein, which alter its biochemical activity. Bacterial phytochromes are dimeric proteins, but the functional relevance of this arrangement remains unclear. Here, we use time-resolved X-ray scattering to reveal the solution structural change of a monomeric variant of the photosensory core module of the phytochrome from Deinococcus radiodurans. The data reveal two motions, a bend and a twist of the PHY domain with respect to the chromophore-binding domains. Infrared spectroscopy shows the refolding of the PHY tongue. We conclude that a monomer of the phytochrome photosensory core is sufficient to perform the light-induced structural changes. This implies that allosteric cooperation with the other monomer is not needed for structural activation. The dimeric arrangement may instead be intrinsic to the biochemical output domains of bacterial phytochromes.

Deciphering solution scattering data with experimentally guided molecular dynamics simulations

Time-resolved X-ray solution scattering is an increasingly popular method to measure conformational changes in proteins. Extracting structural information from the resulting difference X-ray scattering data is a daunting task. We present a method in which the limited but precious information encoded in such scattering curves is combined with the chemical knowledge of molecular force fields. The molecule of interest is then refined toward experimental data using molecular dynamics simulation. Therefore, the energy landscape is biased toward conformations that agree with experimental data. We describe and verify the method, and we provide an implementation in GROMACS.

Light-induced Changes in the Dimerization Interface of Bacteriophytochromes

Phytochromes are dimeric photoreceptor proteins that sense red light levels in plants, fungi, and bacteria. The proteins are structurally divided into a light-sensing photosensory module consisting of PAS, GAF, and PHY domains and a signaling output module, which in bacteriophytochromes typically is a histidine kinase (HK) domain. Existing structural data suggest that two dimerization interfaces exist between the GAF and HK domains, but their functional roles remain unclear. Using mutational, biochemical, and computational analyses of the Deinococcus radiodurans phytochrome, we demonstrate that two dimerization interfaces between sister GAF and HK domains stabilize the dimer with approximately equal contributions. The existence of both dimerization interfaces is critical for thermal reversion back to the resting state. We also find that a mutant in which the interactions between the GAF domains were removed monomerizes under red light. This implies that the interactions between the HK domains are significantly altered by photoconversion. The results suggest functional importance of the dimerization interfaces in bacteriophytochromes.

Structural parameters controlling the fluorescence properties of phytochromes

Phytochromes constitute a class of photoreceptors that can be photoconverted between two stable states. The tetrapyrrole chromophore absorbs in the red spectral region and displays fluorescence maxima above 700 nm, albeit with low quantum yields. Because this wavelength region is particularly advantageous for fluorescence-based deep tissue imaging, there is a strong interest to engineer phytochrome variants with increased fluorescence yields. Such targeted design efforts would substantially benefit from a deeper understanding of those structural parameters that control the photophysical properties of the protein-bound chromophore. Here we have employed resonance Raman (RR) spectroscopy and molecular dynamics simulations for elucidating the chromophore structural changes in a fluorescence-optimized mutant (iRFP) derived from the PAS-GAF domain of the bacteriophytochrome RpBphP2 from Rhodopseudomas palustris . Both methods consistently reveal the structural consequences of the amino acid substitutions in the vicinity of the biliverdin chromophore that may account for lowering the propability of nonradiative excited state decays. First, compared to the wild-type protein, the tilt angle of the terminal ring D with respect to ring C is increased in iRFP, accompanied by the loss of hydrogen bond interactions of the ring D carbonyl function and the reduction of the number of water molecules in that part of the chromophore pocket. Second, the overall flexibility of the chromophore is significantly reduced, particularly in the region of rings D and A, thereby reducing the conformational heterogeneity of the methine bridge between rings A and B and the ring A carbonyl group, as concluded from the RR spectra of the wild-type proteins.

Photoconversion mechanism of the second GAF domain of cyanobacteriochrome AnPixJ and the cofactor structure of its green-absorbing state

Cyanobacteriochromes are members of the phytochrome superfamily. In contrast to classical phytochromes, these small photosensors display a considerable variability of electronic absorption maxima. We have studied the light-induced conversions of the second GAF domain of AnPixJ, AnPixJg2, a phycocyanobilin-binding protein from the cyanobacterium Anabaena PCC 7120, using low-temperature resonance Raman spectroscopy combined with molecular dynamics simulations. AnPixJg2 is formed biosynthetically as a red-absorbing form (Pr) and can be photoconverted into a green-absorbing form (Pg). Forward and backward phototransformations involve the same reaction sequences and intermediates of similar cofactor structures as the corresponding processes in canonical phytochromes, including a transient cofactor deprotonation. Whereas the cofactor of the Pr state shows far-reaching similarities to the Pr states of classical phytochromes, the Pg form displays significant upshifts of the methine bridge stretching frequencies concomitant to the hypsochromically shifted absorption maximum. However, the cofactor in Pg is protonated and adopts a conformation very similar to the Pfr state of classical phytochromes. The spectral differences are probably related to an increased solvent accessibility of the chromophore which may reduce the π-electron delocalization in the phycocyanobilin and thus raise the energies of the first electronic transition and the methine bridge stretching modes. Molecular dynamics simulations suggest that the Z → E photoisomerization of the chromophore at the C-D methine bridge alters the interactions with the nearby Trp90 which in turn may act as a gate, allowing the influx of water molecules into the chromophore pocket. Such a mechanism of color tuning AnPixJg2 is unique among the cyanobacteriochromes studied so far.

Chemical inhomogeneity in the ultrafast dynamics of the DXCF cyanobacteriochrome Tlr0924

Cyanobacteriochromes (CBCRs) are diverse biliprotein photosensors distantly related to the red/far-red photoreceptors of the phytochrome family. There are several subfamilies of CBCRs, displaying varied spectral responses spanning the entire visible region. Tlr0924 belongs to the DXCF subfamily that utilizes the Cys residue in a conserved Asp-Xaa-Cys-Phe (DXCF) motif to form a second covalent linkage to the chromophore, resulting in a blue-absorbing dark state. Photoconversion leads to elimination of this linkage, resulting in a green-absorbing photoproduct. Tlr0924 initially incorporates phycocyanobilin (PCB) as a chromophore, exhibiting a blue/orange photocycle, but slowly isomerizes PCB to phycoviolobilin (PVB) to yield a blue/green photocycle. Ultrafast transient absorption spectroscopy was used to study both forward and reverse reaction photodynamics of the recombinant GAF domain of Tlr0924. Primary photoproducts were identified, as were subsequent intermediates at 1 ms. PCB and PVB population photodynamics were decomposed using global target analysis. PCB and PVB populations exhibit similar and parallel photocycles in Tlr0924, but the PVB population exhibits faster excited-state decay in both reaction directions. On the basis of longer time analysis, we show that the photochemical coordinate (15,16-isomerization) and second-linkage coordinate (elimination or bond formation at C10) are separate processes in both directions.

Femtosecond photodynamics of the red/green cyanobacteriochrome NpR6012g4 from Nostoc punctiforme. 2. reverse dynamics

Phytochromes are red/far-red photosensory proteins that utilize photoisomerization of a linear tetrapyrrole (bilin) chromophore to photoconvert reversibly between red- and far-red-absorbing forms (P(r) and P(fr), respectively). Cyanobacteriochromes (CBCRs) are related photosensory proteins with more diverse spectral sensitivity. The mechanisms that underlie this spectral diversity have not yet been fully elucidated. One of the main CBCR subfamilies photoconverts between a red-absorbing 15Z ground state, like the familiar P(r) state of phytochromes, and a green-absorbing photoproduct ((15E)P(g)). We have previously used the red/green CBCR NpR6012g4 from the cyanobacterium Nostoc punctiforme to examine ultrafast photodynamics of the forward photoreaction. Here, we examine the reverse reaction. Using excitation-interleaved transient absorption spectroscopy with broadband detection and multicomponent global analysis, we observed multiphasic excited-state dynamics. Interleaved excitation allowed us to identify wavelength-dependent shifts in the ground-state bleach that equilibrated on a 200 ps time scale, indicating ground-state heterogeneity. Compared to the previously studied forward reaction, the reverse reaction has much faster excited-state decay time constants and significantly higher photoproduct yield. This work thus demonstrates striking differences between the forward and reverse reactions of NpR6012g4 and provides clear evidence of ground-state heterogeneity in the phytochrome superfamily.

Bacterial phytochromes: more than meets the light

Phytochromes are environmental sensors, historically thought of as red/far-red photoreceptors in plants. Their photoperception occurs through a covalently linked tetrapyrrole chromophore, which undergoes a light-dependent conformational change propagated through the protein to a variable output domain. The phytochrome composition is modular, typically consisting of a PAS-GAF-PHY architecture for the N-terminal photosensory core. A collection of three-dimensional structures has uncovered key features, including an unusual figure-of-eight knot, an extension reaching from the PHY domain to the chromophore-binding GAF domain, and a centrally located, long α-helix hypothesized to be crucial for intramolecular signaling. Continuing identification of phytochromes in microbial systems has expanded the assigned sensory abilities of this family out of the red and into the yellow, green, blue, and violet portions of the spectrum. Furthermore, phytochromes acting not as photoreceptors but as redox sensors have been recognized. In addition, architectures other than PAS-GAF-PHY are known, thus revealing phytochromes to be a varied group of sensory receptors evolved to utilize their modular design to perceive a signal and respond accordingly. This review focuses on the structures of bacterial phytochromes and implications for signal transmission. We also discuss the small but growing set of bacterial phytochromes for which a physiological function has been ascertained.

On the origin of fluorescence in bacteriophytochrome infrared fluorescent proteins

Tsien et al. (Science, 2009, 324, 804-807) recently reported the creation of the first infrared fluorescent protein (IFP). It was engineered from bacterial phytochrome by removing the PHY and histidine kinase-related domains, by optimizing the protein to prevent dimerization, and by limiting the biliverdins conformational freedom, especially around its D ring. We have used database analyses and molecular dynamics simulations with freely rotating chromophoric dihedrals in order to model the dihedral freedom available to the biliverdin D ring in the excited state and to show that the tetrapyrrole ligands in phytochromes are flexible and can adopt many conformations; however, their conformational space is limited/defined by the chemospatial characteristics of the protein cavity. Our simulations confirm that the reduced accessibility to conformations geared to an excited state proton transfer may be responsible for the fluorescence in IFP, just as has been suggested by Kennis et al. (Proc. Natl. Acad. Sci. U.S.A., 2010, 107, 9170-9175) for fluorescent bacteriophytochrome from Rhodopseudomonas palustris.