Molecular Modeling Clarifies the Mechanism of Chromophore Maturation in the Green Fluorescent Protein

Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins

The photochemically active sites of the proteins sfGFP66azF and Venus66azF, members of the green fluorescent protein (GFP) family, contain a non-canonical amino acid residue p-azidophenylalanine (azF) instead of Tyr66. The light-induced decomposition of azF at these sites leads to the formation of reactive arylnitrene (nF) intermediates followed by the formation of phenylamine-containing chromophores. We report the first study of the reaction mechanism of the reduction of the arylnitrene intermediates in sfGFP66nF and Venus66nF using molecular modeling methods. The Gibbs energy profiles for the elementary steps of the chemical reaction in sfGFP66nF are computed using molecular dynamics simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. Structures and energies along the reaction pathway in Venus66nF are evaluated using a QM/MM approach. According to the results of the simulations, arylnitrene reduction is coupled with oxidation of the histidine side chain on the His148 residue located near the chromophore.

Computational Modeling of the Interaction of Molecular Oxygen with the miniSOG Protein—A Light Induced Source of Singlet Oxygen

Interaction of molecular oxygen 3O2 with the flavin-dependent protein miniSOG after light illumination results in creation of singlet oxygen 1O2 and superoxide O2●−. Despite the recently resolved crystal structures of miniSOG variants, oxygen-binding sites near the flavin chromophore are poorly characterized. We report the results of computational studies of the protein−oxygen systems using molecular dynamics (MD) simulations with force-field interaction potentials and quantum mechanics/molecular mechanics (QM/MM) potentials for the original miniSOG and the mutated protein. We found several oxygen-binding pockets and pointed out possible tunnels bridging the bulk solvent and the isoalloxazine ring of the chromophore. These findings provide an essential step toward understanding photophysical properties of miniSOG—an important singlet oxygen photosensitizer.

Applications of fluorescence spectroscopy in protein conformational changes and intermolecular contacts

Emission fluorescence is one of the most versatile and powerful biophysical techniques used in several scientific subjects. It is extensively applied in the studies of proteins, their conformations, and intermolecular contacts, such as in protein-ligand and protein-protein interactions, allowing qualitative, quantitative, and structural data elucidation. This review, aimed to outline some of the most widely used fluorescence techniques in this area, illustrate their applications and display a few examples. At first, the data on the intrinsic fluorescence of proteins is disclosed, mainly on the tryptophan side chain. Predominantly, research to study protein conformational changes, protein interactions, and changes in intensities and shifts of the fluorescence emission maximums were discussed. Fluorescence anisotropy or fluorescence polarization is a measurement of the changing orientation of a molecule in space, concerning the time between the absorption and emission events. Absorption and emission indicate the spatial alignment of the molecule's dipoles relative to the electric vector of the electromagnetic wave of excitation and emitted light, respectively. In other words, if the fluorophore population is excited with vertically polarized light, the emitted light will retain some polarization based on how fast it rotates in solution. Therefore, fluorescence anisotropy can be successfully used in protein-protein interaction investigations. Then, green fluorescent proteins (GFPs), photo-transformable fluorescent proteins (FPs) such as photoswitchable and photoconvertible FPs, and those with Large Stokes Shift (LSS) are disclosed in more detail. FPs are potent tools for the study of biological systems. Their versatility and wide range of colours and properties allow many applications. Finally, the application of fluorescence in life sciences is exposed, especially the application of FPs in fluorescence microscopy techniques with super-resolution that enables precise in vivo photolabeling to monitor the movement and interactions of target proteins.

AlphaFold2 and RoseTTAFold predict posttranslational modifications. Chromophore formation in GFP-like proteins

AlphaFold2 and RoseTTAfold are able to predict, based solely on their sequence whether GFP-like proteins will post-translationally form a chromophore (the part of the protein responsible for fluorescence) or not. Their training has not only taught them protein structure and folding, but also chemistry. The structures of 21 sequences of GFP-like fluorescent proteins that will post-translationally form a chromophore and of 23 GFP-like non-fluorescent proteins that do not have the residues required to form a chromophore were determined by AlphaFold2 and RoseTTAfold. The resultant structures were mined for a series of geometric measurements that are crucial to chromophore formation. Statistical analysis of these measurements showed that both programs conclusively distinguished between chromophore forming and non-chromophore forming proteins. A clear distinction between sequences capable of forming a chromophore and those that do not have the residues required for chromophore formation can be obtained by examining a single measurement—the RMSD of the overlap of the central alpha helices of the crystal structure of S65T GFP and the AlphaFold2 determined structure. Only 10 of the 578 GFP-like proteins in the pdb have no chromophore, yet when AlphaFold2 and RoseTTAFold are presented with the sequences of 44 GFP-like proteins that are not in the pdb they fold the proteins in such a way that one can unequivocally distinguish between those that can and cannot form a chromophore.

Temperature-Dependent Fluorescence of mPlum Fluorescent Protein from 295 to 20 K

The development of bright fluorescent proteins (FPs) emitting beyond 600 nm continues to be of interest both from a fundamental perspective in understanding protein-chromophore interactions and from a practical perspective as these FPs would be valuable for cellular imaging. We previously reported ultrafast spectral observations of the excited-state dynamics in mPlum resulting from interconversion between direct hydrogen bonding and water-mediated hydrogen bonding between the chromophore acylimine carbonyl and the Glu16 side chain. Here, we report temperature-dependent steady-state and time-resolved fluorescence measurements of mPlum and its E16H variant, which does not contain a side-chain permitting hydrogen bonding with the acylimine carbonyl. Lowering the temperature of the system freezes interconversion between the hydrogen-bonding states, thus revealing the spectral signatures of the two states. Analysis of the temperature-dependent spectra assuming Boltzmann populations of the two states yields a 205 cm^-1 energy difference. This value agrees with the predictions from a quantum mechanics/molecular mechanics study of mPlum (198 cm^-1). This study demonstrates the first use of cryogenic spectroscopy to quantify the energetics and timescales of FP chromophore structural states that were only previously obtained from computational methods and further confirms the importance of acylimine hydrogen-bonding dynamics to the fluorescence spectral shifts of red FPs.

A novel violet fluorescent protein contains a unique oxidized tyrosine as the simplest chromophore ever reported in fluorescent proteins

We describe an engineered violet fluorescent protein from the lancelet Branchiostoma floridae (bfVFP). This is the first example of a GFP-like fluorescent protein with a stable fluorescent chromophore lacking an imidazolinone ring; instead, it consists of oxidized tyrosine 68 flanked by glycine 67 and alanine 69. bfVFP contains the simplest chromophore reported in fluorescent proteins and was generated from the yellow protein lanFP10A2 by two synergetic mutations, S148H and C166I. The chromophore structure was confirmed crystallographically and by high-resolution mass spectrometry. The photophysical characteristics of bfVFP (323/430 nm, quantum yield 0.33, and Ec 14,300 M-1  cm-1 ) make it potentially useful for multicolor experiments to expand the excitation range of available FP biomarkers and Förster resonance energy transfer with blue and cyan fluorescent protein acceptors.

A novel fluorescent reporter sensitive to serine mis-incorporation

High-fidelity translation was considered a requirement for living cells. The frozen accident theory suggested that any deviation from the standard genetic code should result in the production of so much mis-made and non-functional proteins that cells cannot remain viable. Studies in bacterial, yeast, and mammalian cells show that significant levels of mistranslation (1–10% per codon) can be tolerated or even beneficial under conditions of oxidative stress. Single tRNA mutants, which occur naturally in the human population, can lead to amino acid mis-incorporation at a codon or set of codons. The rate or level of mistranslation can be difficult or impossible to measure in live cells. We developed a novel red fluorescent protein reporter that is sensitive to serine (Ser) mis-incorporation at proline (Pro) codons. The mCherry Ser151Pro mutant is efficiently produced in Escherichia coli but non-fluorescent. We demonstrated in cells and with purified mCherry protein that the fluorescence of mCherry Ser151Pro is rescued by two different tRNA^Ser gene variants that were mutated to contain the Pro (UGG) anticodon. Ser mis-incorporation was confirmed by mass spectrometry. Remarkably, E. coli tolerated mistranslation rates of ~10% per codon with negligible reduction in growth rate. Conformational sampling simulations revealed that the Ser151Pro mutant leads to significant changes in the conformational freedom of the chromophore precursor, which is indicative of a defect in chromophore maturation. Together our data suggest that the mCherry Ser151 mutants may be used to report Ser mis-incorporation at multiple other codons, further expanding the ability to measure mistranslation in living cells.

Why Are Gly31, Gly33, and Gly35 Highly Conserved in All Fluorescent Proteins?

Green fluorescent protein (GFP)-like fluorescent proteins have been found in more than 120 species. Although the proteins have little sequence identity, Gly31, 33, and 35 are 87, 100, and 95% conserved across all species, respectively. All GFP-like proteins have a β-barrel structure composed of 11 β-sheets, and the 3 conserved glycines are located in the second β-sheet. Molecular dynamics (MD) simulations have shown that mutating one or more of the glycines to alanines most likely does not reduce chromophore formation in correctly folded immature fluorescent proteins. MD and protein characterization of alanine mutants indicate that mutation of the conserved glycines leads to misfolding. Gly31, 33, and 35 are essential to maintain the integrity of the β_1-3 triad that is the last structural element to slot in place in the formation of the canonical fluorescent protein β-barrel. Glycines located in β-sheets may have a similar role in the formation of other non-GFP β-barrels.

Beyond protein tagging: Rewiring the genetic code of fluorescent proteins - A review

Fluorescent proteins (FP) are an integral part of modern biology due to its diverse biochemical and photophysical properties. The boundaries of FP have been extended through conventional mutagenesis and directed evolution approaches. Engineering of FP based on the standard genetic code consisting of 20 amino acids with limited functional groups restrict its diversification. Degeneracy of genetic code has helped in covering this substantial gap through genetic code engineering, wherein introduction of unnatural amino acid (UAA) analogues resulted in a collection of FP with varying properties. This review features the work carried till date in the area of FP incorporated with UAAs and explores strategies employed for incorporation, impact of UAAs in chromophore and surrounding residues and changes in inherent properties of FP. The long-standing association of FP as a tool for high throughput screening of orthogonal aaRS/tRNA pairs used in site specific incorporation of UAAs is expounded. Insertion of UAAs in FP has enabled their use in contemporary fields such as biophotovoltaics, bioremediation, biosensors, biomaterials and imaging of acidic vesicles. Thus, expansion of genetic code of FP is envisaged to rejig the existing spectra of colors and future research initiative in this direction is expected to glow brighter and brighter.

A genetically encoded sensor with improved fluorescence intensity for opioid detection at cellular resolution

The mu-opioid receptor (MOR) regulates the neuronal pathways involved in pain, reward, and respiration. To increase our understanding of MOR's roles in these pathways, there is a need to detect opioids at cellular resolution. Here, we engineered an improved opioid-sensor, called M-SPOTIT2, which is 11x brighter than our previously engineered M-SPOTIT1.1. We engineered M-SPOTIT2 by adding the amino acids YNSH, located near the fluorophore of the enhanced green fluorescent protein, to the circular permuted green fluorescent protein in M-SPOTIT2. M-SPOTIT2 is 11x brighter than our previously engineered M-SPOTIT1.1 in HEK293T cell culture and 2.7x brighter in neuronal culture. M-SPOTIT2 will potentially be useful for the detection of opioids in cell culture for drug screening and the detection of opioids at cellular resolution in animal tissues. By using M-SPOTIT2, researchers can gain more understanding about the mechanisms of addiction, respiratory suppression, and pain-modulation involved in opioid signaling.

Amino acid residue at the 165th position tunes EYFP chromophore maturation. A structure-based design

For the whole GFP family, a few cases, when a single mutation in the chromophore environment strongly inhibits maturation, were described. Here we study EYFP-F165G - a variant of the enhanced yellow fluorescent protein - obtained by a single F165G replacement, and demonstrated multiple fluorescent states represented by the minor emission peaks in blue and yellow ranges (~470 and ~530 nm), and the major peak at ~330 nm. The latter has been assigned to tryptophan fluorescence, quenched due to excitation energy transfer to the mature chromophore in the parental EYFP protein. EYFP-F165G crystal structure revealed two general independent routes of post-translational chemistry, resulting in two main states of the polypeptide chain with the intact chromophore forming triad (~85%) and mature chromophore (~15%). Our experiments thus highlighted important stereochemical role of the 165th position strongly affecting spectral characteristics of the protein. On the basis of the determined EYFP-F165G three-dimensional structure, new variants with ~ 2-fold improved brightness were engineered.

Stalling chromophore synthesis of the fluorescent protein Venus reveals the molecular basis of the final oxidation step

Fluorescent proteins (FPs) have revolutionised the life sciences, but the mechanism of chromophore maturation is still not fully understood. Here we show that incorporation of a photo-responsive non-canonical amino acid within the chromophore stalls maturation of Venus, a yellow FP, at an intermediate stage; a crystal structure indicates the presence of O2 located above a dehydrated enolate form of the imidazolone ring, close to the strictly conserved Gly67 that occupies a twisted conformation. His148 adopts an "open" conformation so forming a channel that allows O2 access to the immature chromophore. Absorbance spectroscopy supported by QM/MM simulations suggests that the first oxidation step involves formation of a hydroperoxyl intermediate in conjunction with dehydrogenation of the methylene bridge. A fully conjugated mature chromophore is formed through release of H2O2, both in vitro and in vivo. The possibility of interrupting and photochemically restarting chromophore maturation and the mechanistic insights open up new approaches for engineering optically controlled fluorescent proteins.

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.

Interplay between Locally Excited and Charge Transfer States Governs the Photoswitching Mechanism in the Fluorescent Protein Dreiklang

We present the results of high-level electronic structure and dynamics simulations of the photoactive protein Dreiklang. With the goal of understanding the details of the Dreiklang photocycle, we carefully characterize the excited states of the ON- and OFF-forms of Dreiklang. The key finding of our study is the existence of a low-lying excited state of a charge-transfer character in the neutral ON form and that population of this state, which is nearly isoenergetic with the locally excited bright state, initiates a series of steps that ultimately lead to the formation of the hydrated dark chromophore (OFF state). These results allow us to refine the mechanistic picture of Dreiklang's photocycle and photoactivation.

Evidence for rate-dependent filtering of global extrinsic noise by biochemical reactions in mammalian cells

Recent studies have revealed that global extrinsic noise arising from stochasticity in the intracellular biochemical environment plays a critical role in heterogeneous cell physiologies. However, it remains largely unclear how such extrinsic noise dynamically influences downstream reactions and whether it could be neutralized by cellular reactions. Here, using fluorescent protein (FP) maturation as a model biochemical reaction, we explored how cellular reactions might combat global extrinsic noise in mammalian cells. We developed a novel single-cell assay to systematically quantify the maturation rate and the associated noise for over a dozen FPs. By exploiting the variation in the maturation rate for different FPs, we inferred that global extrinsic noise could be temporally filtered by maturation reactions, and as a result, the noise levels for slow-maturing FPs are lower compared to fast-maturing FPs. This mechanism is validated by directly perturbing the maturation rates of specific FPs and measuring the resulting noise levels. Together, our results revealed a potentially general principle governing extrinsic noise propagation, where timescale separation allows cellular reactions to cope with dynamic global extrinsic noise.

Influence of the First Chromophore-Forming Residue on Photobleaching and Oxidative Photoconversion of EGFP and EYFP

Enhanced green fluorescent protein (EGFP)—one of the most widely applied genetically encoded fluorescent probes—carries the threonine-tyrosine-glycine (TYG) chromophore. EGFP efficiently undergoes green-to-red oxidative photoconversion (“redding”) with electron acceptors. Enhanced yellow fluorescent protein (EYFP), a close EGFP homologue (five amino acid substitutions), has a glycine-tyrosine-glycine (GYG) chromophore and is much less susceptible to redding, requiring halide ions in addition to the oxidants. In this contribution we aim to clarify the role of the first chromophore-forming amino acid in photoinduced behavior of these fluorescent proteins. To that end, we compared photobleaching and redding kinetics of EGFP, EYFP, and their mutants with reciprocally substituted chromophore residues, EGFP-T65G and EYFP-G65T. Measurements showed that T65G mutation significantly increases EGFP photostability and inhibits its excited-state oxidation efficiency. Remarkably, while EYFP-G65T demonstrated highly increased spectral sensitivity to chloride, it is also able to undergo redding chloride-independently. Atomistic calculations reveal that the GYG chromophore has an increased flexibility, which facilitates radiationless relaxation leading to the reduced fluorescence quantum yield in the T65G mutant. The GYG chromophore also has larger oscillator strength as compared to TYG, which leads to a shorter radiative lifetime (i.e., a faster rate of fluorescence). The faster fluorescence rate partially compensates for the loss of quantum efficiency due to radiationless relaxation. The shorter excited-state lifetime of the GYG chromophore is responsible for its increased photostability and resistance to redding. In EYFP and EYFP-G65T, the chromophore is stabilized by π-stacking with Tyr203, which suppresses its twisting motions relative to EGFP.

Design of a chromogenic substrate for elastase based on split GFP system-Proof of concept for colour switch sensors

Recent studies have demonstrated that human neutrophil elastase (HNE) can be used as marker for inflammation/infection of chronic wounds since it was found to be present in high concentration in exudate collected from chronic wounds. Biosensors used in wound care benefit from a chromogenic signalling due to the readiness of signal interpretation, but the most common use faint yellow chromogenic molecules such as p-nitroaniline (pNa). In addition, if to be converted into smart dressings, the colour of the detection system should not be masked by the exudate's colour. In this work, we designed a chromogenic substrate for HNE aiming to be incorporated in a smart dressing as a colour switch sensor. The substrate was developed using the GFP-like chromoprotein ultramarine (UM), following the split GFP technology. The cleavage sequence for HNE (Ala-Ala-Pro-Val) was embedded into the sensing moiety of the substrate corresponding to the 11th β-sheet. In the presence of HNE, the 11th β-sheet is able to interact to the signalling moiety composed of the β1-β10 incomplete barrel, allowing the re-establishment of the chromophore environment and, hence, the colour production. Structural homology and molecular dynamics simulations were conducted to aid on the disclosure of the structural changes that are the base of the mechanism of action of this HNE switch substrate. Our findings explore the possible application of GFP-like chromogenic sensors in point-of-care devices for the evaluation of the wounds status, representing a major step in the medical field.

Inspecting fluctuation and coordination around chromophore inside green fluorescent protein from water to nonpolar solvent

Green fluorescent protein (GFP) is a widely used biomarker that demands systematical rational approaches to its structure function redesign. In this work, we mainly utilized atomistic molecular dynamics simulations to inspect and visualize internal fluctuation and coordination around chromophore inside GFP, from water to nonpolar octane solvent. We found that GFP not only maintains its β-barrel structure well into the octane, but also sustains internal residue and water coordination to position the chromophore stably while suppress dihedral fluctuations of the chromophore, so that functional robustness of GFP is achieved. Our accompanied fluorescence microscope measurements accordingly confirmed the GFP functioning into the octane. Furthermore, we identified that crucial water sites inside GFP along with permeable pores on the β-barrel of the protein are largely preserved from the water to the octane solvent, which allows sufficiently fast exchanges of internal water with the bulk or with the water layer kept on the surface of the protein. By additionally pulling GFP from bulk water to octane, we suggest that the GFP function can be well maintained into the nonpolar solvent as long as, first, the protein does not denature in the nonpolar solvent nor across the polar-nonpolar solvent interface; second, a minimal set of water molecules are in accompany with the protein; third, the nonpolar solvent molecules may need to be large enough to be nonpermeable via the water pores on the β-barrel.

Assessment of Functionals for TDDFT Calculations of One- and Two-Photon Absorption Properties of Neutral and Anionic Fluorescent Proteins Chromophores

Performance of DFT functionals with different percentages of exact Hartree-Fock exchange energy (EX) is assessed for recovery of the CC2 reference one- (OPA) and two-photon absorption (TPA) spectra of fluorescent proteins chromophores in vacuo. The investigated DFT functionals, together with their EX contributions are BLYP (0%), B3LYP (20%), B1LYP (25%), BHandHLYP (50%), and CAM-B3LYP (19% at short range and 65% at long range). Our test set consists of anionic and neutral chromophores as naturally occurring in the fluorescent proteins. For the first time, we compare TDDFT and CC2 methods for higher excited states than the S₁ state, exhibiting relatively large TPA intensity. Our TDDFT results for neutral chromophores reveal an increase in excitation energies as well as TPA and OPA intensities errors, compared to CC2-derived results, as the DFT functional contains less exact exchange. The long-range-corrected CAM-B3LYP functional performs the best, closely followed by BHandHLYP, while BLYP usually significantly underestimates all investigated spectral properties, hence being the worst in reproducing the reference CC2 results. The hybrid B3LYP and B1LYP functionals can be roughly placed in between. We propose that TDDFT may underestimate the TPA intensities for neutral chromophores of fluorescent proteins due to underestimated oscillator strengths between some excited states. In the case of anionic chromophores, we find that B3LYP and B1LYP functionals overcome others in terms of reproducing CC2 excitation energies. On the other hand, however, TPA intensity is usually significantly underestimated, and in this respect, CAM-B3LYP functional seems to be again superior. In contrast to the case of neutral chromophores, it seems that a large magnitude of excited-state dipole moments or changes in dipole moments upon excitation may be the driving force behind high TPA transition moments.

Engineering 'Golden' Fluorescence by Selective Pressure Incorporation of Non-canonical Amino Acids and Protein Analysis by Mass Spectrometry and Fluorescence

Fluorescent proteins are fundamental tools for the life sciences, in particular for fluorescence microscopy of living cells. While wild-type and engineered variants of the green fluorescent protein from Aequorea victoria (avGFP) as well as homologs from other species already cover large parts of the optical spectrum, a spectral gap remains in the near-infrared region, for which avGFP-based fluorophores are not available. Red-shifted fluorescent protein (FP) variants would substantially expand the toolkit for spectral unmixing of multiple molecular species, but the naturally occurring red-shifted FPs derived from corals or sea anemones have lower fluorescence quantum yield and inferior photo-stability compared to the avGFP variants. Further manipulation and possible expansion of the chromophore's conjugated system towards the far-red spectral region is also limited by the repertoire of 20 canonical amino acids prescribed by the genetic code. To overcome these limitations, synthetic biology can achieve further spectral red-shifting via insertion of non-canonical amino acids into the chromophore triad. We describe the application of SPI to engineer avGFP variants with novel spectral properties. Protein expression is performed in a tryptophan-auxotrophic E. coli strain and by supplementing growth media with suitable indole precursors. Inside the cells, these precursors are converted to the corresponding tryptophan analogs and incorporated into proteins by the ribosomal machinery in response to UGG codons. The replacement of Trp-66 in the enhanced "cyan" variant of avGFP (ECFP) by an electron-donating 4-aminotryptophan results in GdFP featuring a 108 nm Stokes shift and a strongly red-shifted emission maximum (574 nm), while being thermodynamically more stable than its predecessor ECFP. Residue-specific incorporation of the non-canonical amino acid is analyzed by mass spectrometry. The spectroscopic properties of GdFP are characterized by time-resolved fluorescence spectroscopy as one of the valuable applications of genetically encoded FPs in life sciences.

Non-fluorescent mutant of green fluorescent protein sheds light on the mechanism of chromophore formation

The mechanism of green fluorescent protein (GFP) chromophore formation is still not clearly defined. Two mechanisms have been proposed: cyclisation-dehydration-oxidation (Mechanism A) and cyclisation-oxidation-dehydration (Mechanism B). To distinguish between these mechanisms, we generated a non-fluorescent mutant of GFP, S65T/G67A-GFP. This mutant folds to a stable, native-like structure but lacks fluorescence due to interruption of the chromophore maturation process. Mass spectrometric analysis of peptides derived from this mutant reveal that chromophore formation follows only mechanism A, but that the final oxidation reaction is suppressed. This result is unexpected within the pool of examined GFP mutants, since for the wild-type GFP, there is strong support for mechanism B.

The mechanism of a green fluorescent protein proton shuttle unveiled in the time-resolved frequency domain by excited state ab initio dynamics

We simulated an excited state proton transfer in green fluorescent protein by excited state ab initio dynamics, and examined the reaction mechanism in both the time and the frequency domain through a multi resolution wavelet analysis. This original approach allowed us, for the first time, to directly compare the trends of photoactivated vibrations to femtosecond stimulated Raman spectroscopy results, and to give an unequivocal interpretation of the role played by low frequency modes in promoting the reaction. We could attribute the main driving force of the reaction to an important photoinduced softening of the ring-ring orientational motion of the chromophore, thus permitting the tightening of the hydrogen bond network and the opening of the reaction pathway. We also found that both the chromophore (in terms of its inter-ring dihedral angle and phenolic C-O and imidazolinone C-N bond distances) and its pocket (in terms of the inter-molecular oxygen's dihedral angle of the chromophore pocket) relaxations are modulated by low frequency (about 120 cm-1) modes involving the oxygen atoms of the network. This is in agreement with the femtosecond Raman spectroscopy findings in the time-frequency domain. Moreover, the rate in proximity to the Franck Condon region involves a picosecond time scale, with a significant influence from fluctuations of nearby hydrogen bonded residues such as His148. This approach opens a new scenario with ab initio simulations as routinely used tools to understand photoreactivity and the results of advanced time resolved spectroscopy techniques.

Competition between two cysteines in covalent binding of biliverdin to phytochrome domains

In this work, we disclose a mechanism of competing chemical reactions of protein assembly for a bacterial phytochrome using modern methods of molecular modeling.