Unveiling the mechanism of photoinduced isomerization of the photoactive yellow protein (PYP) chromophore

Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response

In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.

Investigation of Ultrafast Configurational Photoisomerization of Bilirubin Using Femtosecond Stimulated Raman Spectroscopy

Phototherapy is an efficient and safe way to reduce high levels of free 4Z,15Z-bilirubin (ZZ-BR) in the serum of newborns. The success of BR phototherapy lies in photoinduced configurational and structural isomerization processes that form excretable isomers. However, the physical picture of photoinduced photoisomerization of ZZ-BR is still unclear. Here, we strategically implement tunable femtosecond stimulated Raman spectroscopy and several time-resolved electronic spectroscopies, assisted by quantum chemical calculations, to dissect the detailed primary configurational isomerization dynamics of free ZZ-BR in organic solvents. The results of this study demonstrate that upon photoexcitation, ultrafast configurational isomerization proceeds by a volume-conserving "hula twist", followed by intramolecular hydrogen-bond distortion and large-scale rotation of the two dipyrrinone halves of the ZZ-BR isomer in a few picoseconds. After that, most of the population recovers back to ZZ-BR, and a very small amount is converted into stable BR isomers via structural isomerization.

Excited-State Barrier Controls E → Z Photoisomerization in p-Hydroxycinnamate Biochromophores

Molecules based on the deprotonated p-hydroxycinnamate moiety are widespread in nature, including serving as UV filters in the leaves of plants and as the biochromophore in photoactive yellow protein. The photophysical behavior of these chromophores is centered around a rapid E → Z photoisomerization by passage through a conical intersection seam. Here, we use photoisomerization and photodissociation action spectroscopies with deprotonated 4-hydroxybenzal acetone (pCK^-) to characterize a wavelength-dependent bifurcation between electron autodetachment (spontaneous ejection of an electron from the S₁ state because it is situated in the detachment continuum) and E → Z photoisomerization. While autodetachment occurs across the entire S₁(ππ*) band (370-480 nm), E → Z photoisomerization occurs only over a blue portion of the band (370-430 nm). No E → Z photoisomerization is observed when the ketone functional group in pCK^- is replaced with an ester or carboxylic acid. The wavelength-dependent bifurcation is consistent with potential energy surface calculations showing that a barrier separates the Franck-Condon region from the E → Z isomerizing conical intersection. The barrier height, which is substantially higher in the gas phase than in solution, depends on the functional group and governs whether E → Z photoisomerization occurs more rapidly than autodetachment.

Four resonance structures elucidate double-bond isomerisation of a biological chromophore

Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.

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.

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.

Controlling electron emission from the photoactive yellow protein chromophore by substitution at the coumaric acid group

Understanding how the interactions between a chromophore and its surrounding protein control the function of a photoactive protein remains a challenge. Here, we present the results of photoelectron spectroscopy measurements and quantum chemistry calculations aimed at investigating how substitution at the coumaryl tail of the photoactive yellow protein chromophore controls competing relaxation pathways following photoexcitation of isolated chromophores in the gas phase with ultraviolet light in the range 350-315 nm. The photoelectron spectra are dominated by electrons resulting from direct detachment and fast detachment from the 2(1)ππ* state but also have a low electron kinetic energy component arising from autodetachment from lower lying electronically excited states or thermionic emission from the electronic ground state. We find that substituting the hydrogen atom of the carboxylic acid group with a methyl group lowers the threshold for electron detachment but has very little effect on the competition between the different relaxation pathways, whereas substituting with a thioester group raises the threshold for electron detachment and appears to 'turn off' the competing electron emission processes from lower lying electronically excited states. This has potential implications in terms of tuning the light-induced electron donor properties of photoactive yellow protein.

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.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

Influence of a chromophore analogue in the protein cage of a photoactive yellow protein

Time-resolved spectra of a photoactive yellow protein (PYP) containing cyano-p-coumaric acid (CHCA) were recorded. To understand the mechanism of photo-isomerization, an electron-withdrawing CN group was introduced into the PYP to alter the C[double bond, length as m-dash]C double bond character. Free CHCA chromophores in aqueous solution underwent photo-isomerization whereas PYP with a bound CHCA (PYP-CN) exhibited no photocycle at acidic or alkaline pH or in urea and other solutions. Furthermore, no photocycle was observed with PYP mutants after illumination. This phenomenon cannot be fully explained by the electron-withdrawing properties of the CN group. We conclude that the CHCA chromophore in PYP was locked in the protein cage and that the CN group interacted with the protein residues.