Ultrafast hydrogen-bonding dynamics in the electronic excited state of photoactive yellow protein revealed by femtosecond stimulated Raman spectroscopy

Photogrammetry of Ultrafast Excited-State Intramolecular Proton Transfer Pathways in the Fungal Pigment Draconin Red

Proton transfer processes of organic molecules are key to charge transport and photoprotection in biological systems. Among them, excited-state intramolecular proton transfer (ESIPT) reactions are characterized by quick and efficient charge transfer within a molecule, resulting in ultrafast proton motions. The ESIPT-facilitated interconversion between two tautomers (PS and PA) comprising the tree fungal pigment Draconin Red in solution was investigated using a combination of targeted femtosecond transient absorption (fs-TA) and excited-state femtosecond stimulated Raman spectroscopy (ES-FSRS) measurements. Transient intensity (population and polarizability) and frequency (structural and cooling) dynamics of -COH rocking and -C=C, -C=O stretching modes following directed stimulation of each tautomer elucidate the excitation-dependent relaxation pathways, particularly the bidirectional ESIPT progression out of the Franck-Condon region to the lower-lying excited state, of the intrinsically heterogeneous chromophore in dichloromethane solvent. A characteristic overall excited-state PS-to-PA transition on the picosecond timescale leads to a unique "W"-shaped excited-state Raman intensity pattern due to dynamic resonance enhancement with the Raman pump-probe pulse pair. The ability to utilize quantum mechanics calculations in conjunction with steady-state electronic absorption and emission spectra to induce disparate excited-state populations in an inhomogeneous mixture of similar tautomers has broad implications for the modeling of potential energy surfaces and delineation of reaction mechanisms in naturally occurring chromophores. Such fundamental insights afforded by in-depth analysis of ultrafast spectroscopic datasets are also beneficial for future development of sustainable materials and optoelectronics.

Mapping the Complete Photocycle that Powers a Large Stokes Shift Red Fluorescent Protein

Large Stokes shift (LSS) red fluorescent proteins (RFPs) are highly desirable for bioimaging advances. The RFP mKeima, with coexisting cis- and trans-isomers, holds significance as an archetypal system for LSS emission due to excited-state proton transfer (ESPT), yet the mechanisms remain elusive. We implemented femtosecond stimulated Raman spectroscopy (FSRS) and various time-resolved electronic spectroscopies, aided by quantum calculations, to dissect the cis- and trans-mKeima photocycle from ESPT, isomerization, to ground-state proton transfer in solution. This work manifests the power of FSRS with global analysis to resolve Raman fingerprints of intermediate states. Importantly, the deprotonated trans-isomer governs LSS emission at 620 nm, while the deprotonated cis-isomer's 520 nm emission is weak due to an ultrafast cis-to-trans isomerization. Complementary spectroscopic techniques as a table-top toolset are thus essential to study photochemistry in physiological environments.

Skeletal Structure of the Chromophore of Photoactive Yellow Protein in the Excited State Investigated by Ultraviolet Femtosecond Stimulated Raman Spectroscopy

We studied ultrafast structural dynamics of photoactive yellow protein (PYP) using ultraviolet femtosecond stimulated Raman spectroscopy. By employing the Raman pump and probe pulses in the ultraviolet region, resonantly enhanced, rich vibrational features of the excited-state chromophore were observed in the fingerprint region. In contrast to the marked spectral change reported for the excited-state chromophore in solution, in the protein, all of the observed Raman bands in the fingerprint region did not show any noticeable spectral shifts nor band shape changes during the excited-state lifetime of PYP. This indicates that the significant skeletal change does not occur on the chromophore in the excited state of PYP and that the trans conformation is retained in its lifetime. Based on the femtosecond Raman data of PYP obtained so far, we discuss a comprehensive picture of the excited-state structural dynamics of PYP.

Twisted Intramolecular Charge Transfer State of a "Push-Pull" Emitter

The excited state Raman spectra of 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran (DCM) in the locally-excited (LE) and the intramolecular charge transfer (ICT) states have been separately measured by time-resolved stimulated Raman spectroscopy. In a polar dimethylsulfoxide solution, the ultrafast ICT of DCM with a time constant of 1.0 ps was observed in addition to the vibrational relaxation in the ICT state of 4–7 ps. On the other hand, the energy of the ICT state of DCM becomes higher than that of the LE state in a less polar chloroform solution, where the initially-photoexcited ICT state with the LE state shows the ultrafast internal conversion to the LE state with a time constant of 300 fs. The excited-state Raman spectra of the LE and ICT state of DCM showed several major vibrational modes of DCM in the LE and ICT conformer states coexisting in the excited state. Comparing to the time-dependent density functional theory simulations and the experimental results of similar push-pull type molecules, a twisted geometry of the dimethylamino group is suggested for the structure of DCM in the S₁/ICT state.

Mapping Structural Dynamics of Proteins with Femtosecond Stimulated Raman Spectroscopy

The structure-function relationships of biomolecules have captured the interest and imagination of the scientific community and general public since the field of structural biology emerged to enable the molecular understanding of life processes. Proteins that play numerous functional roles in cellular processes have remained in the forefront of research, inspiring new characterization techniques. In this review, we present key theoretical concepts and recent experimental strategies using femtosecond stimulated Raman spectroscopy (FSRS) to map the structural dynamics of proteins, highlighting the flexible chromophores on ultrafast timescales. In particular, wavelength-tunable FSRS exploits dynamic resonance conditions to track transient-species-dependent vibrational motions, enabling rational design to alter functions. Various ways of capturing excited-state chromophore structural snapshots in the time and/or frequency domains are discussed. Continuous development of experimental methodologies, synergistic correlation with theoretical modeling, and the expansion to other nonequilibrium, photoswitchable, and controllable protein systems will greatly advance the chemical, physical, and biological sciences.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Unveiling coupled electronic and vibrational motions of chromophores in condensed phases

The quest for capturing molecular movies of functional systems has motivated scientists and engineers for decades. A fundamental understanding of electronic and nuclear motions, two principal components of the molecular Schrödinger equation, has the potential to enable the de novo rational design for targeted functionalities of molecular machines. We discuss the development and application of a relatively new structural dynamics technique, femtosecond stimulated Raman spectroscopy with broadly tunable laser pulses from the UV to near-IR region, in tracking the coupled electronic and vibrational motions of organic chromophores in solution and protein environments. Such light-sensitive moieties hold broad interest and significance in gaining fundamental knowledge about the intramolecular and intermolecular Hamiltonian and developing effective strategies to control macroscopic properties. Inspired by recent experimental and theoretical advances, we focus on the in situ characterization and spectroscopy-guided tuning of photoacidity, excited state proton transfer pathways, emission color, and internal conversion via a conical intersection.

Does the wavelength dependent photoisomerization process of the p‑coumaric acid come out from the electronic state dependent pathways?

Similar to the anion photoactive yellow protein (PYP) chromophore, the neutral form of the PYP chromophore was also found to exhibit a the wavelength-dependent photoisomerization quantum yield. The isomerization quantum yield increases with the increasing excitation energy on the S₁ state, while decreases when being excited to the S₂ state. Does this wavelength dependent product yield come out from the specific reaction pathways of the S₁ and S₂ states? This would mean that, the relaxation pathway of the S₂ state is distinct from that of the S₁ state and does not involve twisting motion. Does it break Kasha's rule by exhibiting a direct transition from the S₂ state to the ground state? The underlying mechanism needs further in. In this article, we employed the on-the-fly dynamics simulations and static electronic structure calculations to reveal the deactivation mechanism of the neutral form of the PYP chromophore. Our results indicated that the CC twisting motion dominates the S₁ state decay process. In contrast, for the decay process of the S₂ state, an ultrafast transition from the S₂ to the S₁ state through a planar conical intersection is observed, and the excess energy activates a new reaction channel to the ground state characterized by a puckering distortion of the ring. This pathway competes with the photoisomerization channel. No direct transition from S₂ to S₀ is observed, hence Kasha's rule is valid for this process. Our calcualtions can provide a reasonable explanation of the wavelength-dependent isomerization quantum yield of neutral PYP chromophore, and we hope it can provide theoretical foundations for comparing the effect of protonation state on the dynamcal behaviors of PYP chromophore.

Excited State Structural Evolution of a GFP Single-Site Mutant Tracked by Tunable Femtosecond-Stimulated Raman Spectroscopy

Tracking vibrational motions during a photochemical or photophysical process has gained momentum, due to its sensitivity to the progression of reaction and change of environment. In this work, we implemented an advanced ultrafast vibrational technique, femtosecond-stimulated Raman spectroscopy (FSRS), to monitor the excited state structural evolution of an engineered green fluorescent protein (GFP) single-site mutant S205V. This mutation alters the original excited state proton transfer (ESPT) chain. By strategically tuning the Raman pump to different wavelengths (i.e., 801, 539, and 504 nm) to achieve pre-resonance with transient excited state electronic bands, the characteristic Raman modes of the excited protonated (A*) chromophore species and intermediate deprotonated (I*) species can be selectively monitored. The inhomogeneous distribution/population of A* species go through ESPT with a similar ~300 ps time constant, confirming that bridging a water molecule to protein residue T203 in the ESPT chain is the rate-limiting step. Some A* species undergo vibrational cooling through high-frequency motions on the ~190 ps time scale. At early times, a portion of the largely protonated A* species could also undergo vibrational cooling or return to the ground state with a ~80 ps time constant. On the photoproduct side, a ~1330 cm⁻¹ delocalized motion is observed, with dispersive line shapes in both the Stokes and anti-Stokes FSRS with a pre-resonance Raman pump, which indicates strong vibronic coupling, as the mode could facilitate the I* species to reach a relatively stable state (e.g., the main fluorescent state) after conversion from A*. Our findings disentangle the contributions of various vibrational motions active during the ESPT reaction, and offer new structural dynamics insights into the fluorescence mechanisms of engineered GFPs and other analogous autofluorescent proteins.

Excitation-Wavelength-Dependent Photocycle Initiation Dynamics Resolve Heterogeneity in the Photoactive Yellow Protein from Halorhodospira halophila

Photoactive yellow proteins (PYPs) make up a diverse class of blue-light-absorbing bacterial photoreceptors. Electronic excitation of the p-coumaric acid chromophore covalently bound within PYP results in triphasic quenching kinetics; however, the molecular basis of this behavior remains unresolved. Here we explore this question by examining the excitation-wavelength dependence of the photodynamics of the PYP from Halorhodospira halophila via a combined experimental and computational approach. The fluorescence quantum yield, steady-state fluorescence emission maximum, and cryotrapping spectra are demonstrated to depend on excitation wavelength. We also compare the femtosecond photodynamics in PYP at two excitation wavelengths (435 and 475 nm) with a dual-excitation-wavelength-interleaved pump-probe technique. Multicompartment global analysis of these data demonstrates that the excited-state photochemistry of PYP depends subtly, but convincingly, on excitation wavelength with similar kinetics with distinctly different spectral features, including a shifted ground-state beach and altered stimulated emission oscillator strengths and peak positions. Three models involving multiple excited states, vibrationally enhanced barrier crossing, and inhomogeneity are proposed to interpret the observed excitation-wavelength dependence of the data. Conformational heterogeneity was identified as the most probable model, which was supported with molecular mechanics simulations that identified two levels of inhomogeneity involving the orientation of the R52 residue and different hydrogen bonding networks with the p-coumaric acid chromophore. Quantum calculations were used to confirm that these inhomogeneities track to altered spectral properties consistent with the experimental results.

Combined TDDFT and AIM Insights into Photoinduced Excited State Intramolecular Proton Transfer (ESIPT) Mechanism in Hydroxyl- and Amino-Anthraquinone Solution

Time-dependent density functional theory (TDDFT) and atoms in molecules (AIM) theory are combined to study the photoinduced excited state intramolecular proton transfer (ESIPT) dynamics for eight anthraquinones (AQs) derivatives in solution. The calculated absorption and emission spectra are consistent with the available experimental data, verifying the suitability of the theory selected. The systems with the excited-state exothermic proton transfer, such as 1-HAQ, 1,5-DHAQ and TFAQ, emit completely from transfer structure (T), while the reactions for those without ESIPT including 1,4-DHAQ and AAAQ appear to be endothermic. Three reaction properties of three systems (1,8-DHAQ, DCAQ and CAAQ) are between the exothermic and endothermic, sensitive to the solvent. Energy scanning shows that 1,4-DHAQ and AAAQ exhibit the higher ESIPT energy barriers compared to 1-HAQ, 1,5-DHAQ and TFAQ with the "barrierless" ESIPT process. The ESIPT process is facilitated by the strengthening of hydrogen bonds in excited state. With AIM theory, it is observed that the change in electrons density ρ(r) and potential energy density V(r) at BCP position between ground state and excited state are crucial factors to quantitatively elucidate the ESIPT.

Tracking Ultrafast Vibrational Cooling during Excited-State Proton Transfer Reaction with Anti-Stokes and Stokes Femtosecond Stimulated Raman Spectroscopy

Energy dissipation following photoexcitation is foundational to photophysics and chemistry. Consequently, understanding such processes on molecular time scales holds paramount importance. Femtosecond stimulated Raman spectroscopy (FSRS) has been used to study the molecular structure-function relationships but usually on the Stokes side. Here, we perform both Stokes and anti-Stokes FSRS to track energy dissipation and excited-state proton transfer (ESPT) for the photoacid pyranine in aqueous solution. We reveal biphasic vibrational cooling on fs-ps time scales during ESPT. Characteristic low-frequency motions (<800 cm^-1) exhibit initial energy dissipation (∼2 ps) that correlates with functional events of forming contact ion pairs via H-bonds between photoacid and water, which lengthens to ∼9 ps in methanol where ESPT is inhibited. The interplay between photoinduced dissipative and reactive channels is implied. Thermal cooling to bulk solvent occurs on the ∼50 ps time scale. These results demonstrate the combined Stokes and anti-Stokes FSRS as a powerful toolset to elucidate structural dynamics.

QM/MM Study of Substituent and Solvent Effects on the Excited State Dynamics of the Photoactive Yellow Protein Chromophore

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

Probing the early stages of photoreception in photoactive yellow protein with ultrafast time-domain Raman spectroscopy

Unveiling the nuclear motions of photoreceptor proteins in action is a crucial goal in protein science in order to understand their elaborate mechanisms and how they achieve optimal selectivity and efficiency. Previous studies have provided detailed information on the structures of intermediates that appear during the later stages (>ns) of such photoreception cycles, yet the initial events immediately after photoabsorption remain unclear because of experimental challenges in monitoring nuclear rearrangements on ultrafast timescales, including protein-specific low-frequency motions. Using time-domain Raman probing with sub-7-fs pulses, we obtain snapshot vibrational spectra of photoactive yellow protein and a mutant with high sensitivity, providing insights into the key responses that drive photoreception. Our data show a drastic intensity drop of the excited-state marker band at 135 cm^-1 within a few hundred femtoseconds, suggesting a rapid weakening of the hydrogen bond that anchors the chromophore. We also track formation of the first ground-state intermediate over the first few picoseconds and fully characterize its vibrational structure, revealing a substantially-twisted cis conformation.

Solvent effects on de-excitation channels in the p-coumaric acid methyl ester anion, an analogue of the photoactive yellow protein (PYP) chromophore

Environmental effects on the deactivation channels of the PYP chromophore in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level.

How fissors works: observing vibrationally adiabatic conformational change through femtosecond stimulated Raman spectroscopy

With the help of a two-dimensional model system comprising a slow conformational degree of freedom and a higher-frequency vibration, we investigate the molecular-level origin and dynamical information content of femtosecond stimulated Raman spectroscopy (fissors) signals. Our treatment avails itself of the time scale separation between conformational and vibrational modes by incorporating a vibrationally adiabatic approximation to the conformational dynamics. We derive an expression for the fissors signal without resort to the macroscopic concepts of light- and phonon-wave propagation employed in prior coupled-wave analyses. Numerical calculations of fissors spectra illustrate the case of relatively small conformational mass (still large enough that conformational motion does not induce any change in the vibrational quantum number) in which conformational sidebands accompany a central peak in the Raman gain at a conformationally averaged vibrational transition frequency, and the case of a larger conformational mass in which the sidebands merge with the central peak and the frequency of the latter tracks the time-evolving conformational coordinate.