Spectral tuning in photoactive yellow protein by modulation of the shape of the excited state energy surface

From femtoseconds to minutes: Spectroscopic study of optically induced thermal diffusion in aqueous solution of rhodamine B

Given that non-equilibrium molecular motion in thermal gradients is influenced by both solute and solvent, the application of spectroscopic methods that probe each component in a binary mixture can provide insights into the molecular mechanisms of thermal diffusion for a large class of systems. In the present work, we use an all-optical setup whereby near-infrared excitation of the solvent leads to a steady-state thermal gradient in solution, followed by characterization of the non-equilibrium system with electronic spectroscopy, imaging, and intensity. Using rhodamine B in water as a case study, we perform measurements as a function of solute concentration, temperature, wavelength, time, near-infrared laser power, visible excitation wavelength, and isotope effect. We find that spectroscopic excitation of solute and solvent as a function of temperature can be used to understand the thermal diffusion signal. Non-equilibrium molecular dynamics simulations and analysis of infrared spectra and heat diffusion, and stochastic dynamics simulations of the coupled Brownian/spectroscopic non-equilibrium dynamics. The stochastic/spectroscopic model shows how near-infrared excitation of the solvent influences lifetime dynamics on the picosecond timescale as well as Brownian dynamics in real time. Overall, the results presented here exemplify how spectroscopic probing of solute and solvent can be useful for understanding molecular mechanisms of optically induced thermal diffusion in aqueous solution of rhodamine B.

Not All Photoactive Yellow Proteins Are Built Alike: Surprises and Insights into Chromophore Photoisomerization, Protonation, and Thermal Reisomerization of the Photoactive Yellow Protein Isolated from Salinibacter ruber

We demonstrate that the Halorhodospira halophila (Hhal) photoactive yellow protein (PYP) is not representative of the greater PYP family. The photodynamics of the PYP isolated from Salinibacter ruber (Srub) is characterized with a comprehensive range of spectroscopic techniques including ultrafast transient absorption, photostationary light titrations, Fourier transform infrared, and cryokinetics spectroscopies. We demonstrate that the dark-adapted pG state consists of two subpopulations differing in the protonation state of the chromophore and that both are photoactive, with the protonated species undergoing excited-state proton transfer. However, the primary I₀ photoproduct observed in the Hhal PYP photocycle is absent in the Srub PYP photodynamics, which indicates that this intermediate, while important in Hhal photodynamics, is not a critical intermediate in initiating all PYP photocycles. The excited-state lifetime of Srub PYP is the longest of any PYP resolved to date (∼30 ps), which we ascribe to the more constrained chromophore binding pocket of Srub PYP and the absence of the critical Arg52 residue found in Hhal PYP. The final stage of the Srub PYP photocycle involves the slowest known thermal dark reversion of a PYP (∼40 min vs 350 ms in Hhal PYP). This property allowed the characterization of a pH-dependent equilibrium between the light-adapted pB state with a protonated cis chromophore and a newly resolved pG' intermediate with a deprotonated cis chromophore and pG-like protein conformation. This result demonstates that protein conformational changes and chromophore deprotonation precede chromophore reisomerization during the thermal recovery of the PYP photocycle.

Orthogonal fluorescent chemogenetic reporters for multicolor imaging

Spectrally separated fluorophores allow the observation of multiple targets simultaneously inside living cells, leading to a deeper understanding of the molecular interplay that regulates cell function and fate. Chemogenetic systems combining a tag and a synthetic fluorophore provide certain advantages over fluorescent proteins since there is no requirement for chromophore maturation. Here, we present the engineering of a set of spectrally orthogonal fluorogen activating tags based on the Fluorescence Activating and absorption Shifting Tag (FAST), that are compatible with two-color, live cell imaging. The resulting tags, greenFAST and redFAST, demonstrate orthogonality not only in their fluorogen recognition capabilities, but also in their one- and two-photon absorption profiles. This pair of orthogonal tags allowed the creation of a two-color cell cycle sensor capable of detecting very short, early cell cycles in zebrafish development, and the development of split complementation systems capable of detecting multiple protein-protein interactions by live cell fluorescence microscopy.

The Origin of Ultrafast Multiphasic Dynamics in Photoisomerization of Bacteriophytochrome

Red-light bacteriophytochromes regulate many physiological functions through photoisomerization of a linear tetrapyrrole chromophore. In this work, we mapped out femtosecond-resolved fluorescence spectra of the excited Pr state and observed unique active-site relaxations on the picosecond time scale with unusual spectral tuning of rises on the blue side and decays on the red side of the emission. We also observed initial wavepacket dynamics in femtoseconds with two low-frequency modes of 38 and 181 cm^-1 as well as the intermediate product formation after isomerization in hundreds of picoseconds. With critical mutations at the active site, we observed similar dynamic patterns with different times for both relaxation and isomerization, consistent with the structural and chemical changes induced by the mutations. The observed multiphasic dynamics clearly represents the active-site relaxation, not different intermediate reactions or excitation of heterogeneous ground states. The active-site relaxation must be considered in understanding overall isomerization reactions in phytochromes, and such a molecular mechanism should be general.

Electrostatic control of photoisomerization pathways in proteins

Rotation around a specific bond after photoexcitation is central to vision and emerging opportunities in optogenetics, super-resolution microscopy, and photoactive molecular devices. Competing roles for steric and electrostatic effects that govern bond-specific photoisomerization have been widely discussed, the latter originating from chromophore charge transfer upon excitation. We systematically altered the electrostatic properties of the green fluorescent protein chromophore in a photoswitchable variant, Dronpa2, using amber suppression to introduce electron-donating and electron-withdrawing groups to the phenolate ring. Through analysis of the absorption (color), fluorescence quantum yield, and energy barriers to ground- and excited-state isomerization, we evaluate the contributions of sterics and electrostatics quantitatively and demonstrate how electrostatic effects bias the pathway of chromophore photoisomerization, leading to a generalized framework to guide protein design.

Hydrogen-Bonding Interaction Regulates Photoisomerization of a Single-Bond-Rotation Locked Photoactive Yellow Protein Chromophore in Protein

We have employed the QM(CASPT2//CASSCF)/MM method to explore the excited-state isomerization and decay mechanism of a single-bond-rotation locked photoactive yellow protein (PYP) chromophore in wild-type and mutant proteins. The S₁ state is a spectroscopically bright state in the Franck-Condon region. In this state, there exist two excited-state isomerization pathways separately related to the clockwise and anticlockwise rotations of the C=C bond. The clockwise path is favorable because of a small barrier of 2 kcal/mol and uses a novel bicycle-pedal unidirectional photoisomerization mechanism in which the involved two dihedral angles rotate asynchronously because of the reinforced hydrogen-bonding interaction between the chromophore and Cys69. Near the twisted S₁ minimum, the chromophore hops to the S₀ state via the S₁/S₀ conical intersection. Finally, the R52A mutation has small effects on the excited-state properties and photoisomerization of the locked PYP chromophore. The present work provides new insights for understanding the photochemistry of PYP chromophores in protein surroundings.

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.

Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins

Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.

The Dynamic Origin of Color Tuning in Proteins Revealed by a Carotenoid Pigment

Understanding the microscopic origin of the color tuning in pigment-protein complexes is a challenging yet fundamental issue in photoactive biological systems. Here, we propose a possible interpretation by using a state-of-the-art multiscale strategy based on the integration of quantum chemistry and polarizable atomistic embeddings into a dynamic description. By means of such a strategy we are able to resolve the long-standing dispute over the coloration mechanism in the crustacyanin protein. It is shown that the combination of the dynamical flexibility of the carotenoid pigments (astaxanthin) with the responsive protein environment is essential to obtain quantitative predictions of the spectral tuning. The strong linear correlation between the excitation energies and the bond length alternation in the long-chain carotenoids modulated by the dynamical protein environment is a novel finding explaining the high color tunability in crustacyanin.

Mechanism and bottlenecks in strand photodissociation of split green fluorescent proteins (GFPs)

Split GFPs have been widely applied for monitoring protein-protein interactions by expressing GFPs as two or more constituent parts linked to separate proteins that only fluoresce on complementing with one another. Although this complementation is typically irreversible, it has been shown previously that light accelerates dissociation of a noncovalently attached β-strand from a circularly permuted split GFP, allowing the interaction to be reversible. Reversible complementation is desirable, but photodissociation has too low of an efficiency (quantum yield <1%) to be useful as an optogenetic tool. Understanding the physical origins of this low efficiency can provide strategies to improve it. We elucidated the mechanism of strand photodissociation by measuring the dependence of its rate on light intensity and point mutations. The results show that strand photodissociation is a two-step process involving light-activated cis-trans isomerization of the chromophore followed by light-independent strand dissociation. The dependence of the rate on temperature was then used to establish a potential energy surface (PES) diagram along the photodissociation reaction coordinate. The resulting energetics-function model reveals the rate-limiting process to be the transition from the electronic excited-state to the ground-state PES accompanying cis-trans isomerization. Comparisons between split GFPs and other photosensory proteins, like photoactive yellow protein and rhodopsin, provide potential strategies for improving the photodissociation quantum yield.

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.

A PYP chromophore acts as a 'photoacid' in an isolated hydrogen bonded complex

The light-induced response of a neutral Photoactive Yellow Protein chromophore in a hydrogen-bonded complex with a proton acceptor has been studied by dual-detection action absorption spectroscopy and density functional theory. We show that the chromophore is a 'photoacid' and that ultrafast excited-state proton transfer might be operative in an isolated complex.

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.

Dynamics and mechanisms of DNA repair by photolyase

Photolyases, a class of flavoproteins, use blue light to repair two types of ultraviolet-induced DNA damage, a cyclobutane pyrimidine dimer (CPD) and a pyrimidine-pyrimidone (6-4) photoproduct (6-4PP). In this perspective, we review the recent progress in the repair dynamics and mechanisms of both types of DNA restoration by photolyases. We first report the spectroscopic characterization of flavin in various redox states and the active-site solvation dynamics in photolyases. We then systematically summarize the detailed repair dynamics of damaged DNA by photolyases and a biomimetic system through resolving all elementary steps on ultrafast timescales, including multiple intermolecular electron- and proton-transfer reactions and bond-breaking and -making processes. We determined the unique electron tunneling pathways, identified the key functional residues and revealed the molecular origin of high repair efficiency, and thus elucidate the molecular mechanisms and repair photocycles at the most fundamental level. We finally conclude that the active sites of photolyases, unlike the aqueous solution for the biomimetic system, provide a unique electrostatic environment and local flexibility and thus a dedicated synergy for all elementary dynamics to maximize the repair efficiency. This repair photomachine is the first enzyme that the entire functional evolution is completely mapped out in real time.

Strong ionic hydrogen bonding causes a spectral isotope effect in photoactive yellow protein

Standard hydrogen bonds are of great importance for protein structure and function. Ionic hydrogen bonds often are significantly stronger than standard hydrogen bonds and exhibit unique properties, but their role in proteins is not well understood. We report that hydrogen/deuterium exchange causes a redshift in the visible absorbance spectrum of photoactive yellow protein (PYP). We expand the range of interpretable isotope effects by assigning this spectral isotope effect (SIE) to a functionally important hydrogen bond at the active site of PYP. The inverted sign and extent of this SIE is explained by the ionic nature and strength of this hydrogen bond. These results show the relevance of ionic hydrogen bonding for protein active sites, and reveal that the inverted SIE is a novel, to our knowledge, tool to probe ionic hydrogen bonds. Our results support a classification of hydrogen bonds that distinguishes the properties of ionic hydrogen bonds from those of both standard and low barrier hydrogen bonds, and show how this classification helps resolve a recent debate regarding active site hydrogen bonding in PYP.

Towards in situ fluorescence spectroscopy and microscopy investigations of asphaltene precipitation kinetics

We perform a spectroscopic analysis of asphaltene in solution and in crude oil with the goal of designing an optical probe of asphaltene precipitation inside high-pressure cells. Quantitative analysis of steady-state spectroscopic data is employed to identify fluorescence and Raman contributions to the observed signals. Time-resolved fluorescence spectroscopy indicates that fluorescence lifetime can be used as a spectroscopic probe of asphaltene in crude oil. Quantitative confocal laser-scanning microscopy studies of asphaltene in n-heptane are used to calculate particle-size distributions as a function of time, both at the sample surface and asphaltene interior. The resulting precipitation kinetics is well described by stochastic numerical simulations of diffusion-limited aggregation. Based on these results, we present the design and construction of an apparatus to optically probe the in situ precipitation of asphaltene suitable for studies inside high pressure cells. Design considerations include the use of a spatial light modulator for aberration correction in microscopy measurements, together with the design of epi-fluorescence spectrometer, both fiber-based and for remote sensing fluorescence spectroscopy.

Nonlinear spectroscopic theory of displaced harmonic oscillators with differing curvatures: a correlation function approach

We present a theory for a bath model in which we approximate the adiabatic nuclear potential surfaces on the ground and excited electronic states by displaced harmonic oscillators that differ in curvature. Calculations of the linear and third-order optical response functions employ an effective short-time approximation coupled with the cumulant expansion. In general, all orders of correlation contribute to the optical response, indicating that the solvation process cannot be described as Gaussian within the model. Calculations of the linear absorption and fluorescence spectra resulting from the theory reveal a stronger temperature dependence of the Stokes shift along with a general asymmetry between absorption and fluorescence line shapes, resulting purely from the difference in the phonon side band. We discuss strategies for controlling spectral tuning and energy-transfer dynamics through the manipulation of the excited-state and ground-state curvature. Calculations of the nonlinear response also provide insights into the dynamics of the system-bath interactions and reveal that multidimensional spectroscopies are sensitive to a difference in curvature between the ground- and excited-state adiabatic surfaces. This extension allows for the elucidation of short-time dynamics of dephasing that are accessible in nonlinear spectroscopic methods.

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

The ultrafast structural dynamics in the electronic excited state of photoactive yellow protein (PYP) is studied by femtosecond stimulated Raman spectroscopy. Stimulated Raman spectra in the electronic excited state, S(1), can be obtained by using a Raman pump pulse in resonance with the S(1)-S(0) transition. This is confirmed by comparing the experimental results with numerical calculations based on the density matrix treatment. We also investigate the hydrogen-bonding network surrounding the wild-type (WT)-PYP chromophore in the ground and excited states by comparing its stimulated Raman spectra with those of the E46Q-PYP mutant. We focus on the relative intensity of the Raman band at 1555 cm(-1), which includes both vinyl bond C═C stretching and ring vibrations and is sensitive to the hydrogen-bonding network around the phenolic oxygen of the chromophore. The relative intensity for the WT-PYP decreases after actinic excitation within the 150 fs time resolution and reaches a similar intensity to that for E46Q-PYP. These observations indicate that the WT-PYP hydrogen-bonding network is immediately rearranged in the electronic excited state to form a structure similar to that of E46Q-PYP.

Electronic Absorption Spectra from MM and ab initio QM/MM Molecular Dynamics: Environmental Effects on the Absorption Spectrum of Photoactive Yellow Protein

We describe a new interface of the GPU parallelized TeraChem electronic structure package and the Amber molecular dynamics package for quantum mechanical (QM) and mixed QM and molecular mechanical (MM) molecular dynamics simulations. This QM/MM interface is used for computation of the absorption spectra of the photoactive yellow protein (PYP) chromophore in vacuum, aqueous solution, and protein environments. The computed excitation energies of PYP require a very large QM region (hundreds of atoms) covalently bonded to the chromophore in order to achieve agreement with calculations that treat the entire protein quantum mechanically. We also show that 40 or more surrounding water molecules must be included in the QM region in order to obtain converged excitation energies of the solvated PYP chromophore. These results indicate that large QM regions (with hundreds of atoms) are a necessity in QM/MM calculations.

On the involvement of single-bond rotation in the primary photochemistry of photoactive yellow protein

Prior experimental observations, as well as theoretical considerations, have led to the proposal that C(4)-C(7) single-bond rotation may play an important role in the primary photochemistry of photoactive yellow protein (PYP). We therefore synthesized an analog of this protein's 4-hydroxy-cinnamic acid chromophore, (5-hydroxy indan-(1E)-ylidene)acetic acid, in which rotation across the C(4)-C(7) single bond has been locked with an ethane bridge, and we reconstituted the apo form of the wild-type protein and its R52A derivative with this chromophore analog. In PYP reconstituted with the rotation-locked chromophore, 1), absorption spectra of ground and intermediate states are slightly blue-shifted; 2), the quantum yield of photochemistry is ∼60% reduced; 3), the excited-state dynamics of the chromophore are accelerated; and 4), dynamics of the thermal recovery reaction of the protein are accelerated. A significant finding was that the yield of the transient ground-state intermediate in the early phase of the photocycle was considerably higher in the rotation-locked samples than in the corresponding samples reconstituted with p-coumaric acid. In contrast to theoretical predictions, the initial photocycle dynamics of PYP were observed to be not affected by the charge of the amino acid residue at position 52, which was varied by 1), varying the pH of the sample between 5 and 10; and 2), site-directed mutagenesis to construct R52A. These results imply that C(4)-C(7) single-bond rotation in PYP is not an alternative to C(7)=C(8) double-bond rotation, in case the nearby positive charge of R52 is absent, but rather facilitates, presumably with a compensatory movement, the physiological Z/E isomerization of the blue-light-absorbing chromophore.

Modulation of the photocycle of a LOV domain photoreceptor by the hydrogen-bonding network

An extended hydrogen-bonding (HB) network stabilizes the isoalloxazine ring of the flavin mononucleotide (FMN) chromophore within the photosensing LOV domain of blue-light protein receptors, via interactions between the C(2)═O, N(3)H, C(4)═O, and N(5) groups and conserved glutamine and asparagine residues. In this work we studied the influence of the HB network on the efficiency, kinetics, and energetics of a LOV protein photocycle, involving the reversible formation of a FMN-cysteine covalent adduct. The following results were found for mutations of the conserved amino acids N94, N104, and Q123 in the Bacillus subtilis LOV protein YtvA: (i) Increased (N104D, N94D) or strongly reduced (N94A) rate of adduct formation; this latter mutation extends the lifetime of the flavin triplet state, i.e., adduct formation, more than 60-fold, from 2 μs for the wild-type (WT) protein to 129 μs. (ii) Acceleration of the overall photocycle for N94S, N94A, and Q123N, with recovery lifetimes 20, 45, and 85 times faster than for YtvA-WT, respectively. (iii) Slight modifications of FMN spectral features, correlated with the polarization of low-energy transitions. (iv) Strongly reduced (N94S) or suppressed (Q123N) structural volume changes accompanying adduct formation, as determined by optoacoustic spectroscopy. (v) Minor effects on the quantum yield, with the exception of a considerable reduction for Q123N, i.e., 0.22 vs 0.49 for YtvA-WT. The data stress the importance of the HB network in modulating the photocycle of LOV domains, while at the same time establishing a link with functional responses.

A conserved helical capping hydrogen bond in PAS domains controls signaling kinetics in the superfamily prototype photoactive yellow protein

PAS domains form a divergent protein superfamily with more than 20 000 members that perform a wide array of sensing and regulatory functions in all three domains of life. Only nine residues are well-conserved in PAS domains, with an Asn residue at the start of α-helix 3 showing the strongest conservation. The molecular functions of these nine conserved residues are unknown. We use static and time-resolved visible and FTIR spectroscopy to investigate receptor activation in the photosensor photoactive yellow protein (PYP), a PAS domain prototype. The N43A and N43S mutants allow an investigation of the role of side-chain hydrogen bonding at this conserved position. The mutants exhibit a blue-shifted visible absorbance maximum and up-shifted chromophore pK(a). Disruption of the hydrogen bonds in N43A PYP causes both a reduction in protein stability and a 3400-fold increase in the lifetime of the signaling state of this photoreceptor. A significant part of this increase in lifetime can be attributed to the helical capping interaction of Asn43. This extends the known importance of helical capping for protein structure to regulating functional protein kinetics. A model for PYP activation has been proposed in which side-chain hydrogen bonding of Asn43 is critical for relaying light-induced conformational changes. However, FTIR spectroscopy shows that both Asn43 mutants retain full allosteric transmission of structural changes. Analysis of 30 available high-resolution structures of PAS domains reveals that the side-chain hydrogen bonding of residue 43 but not residue identity is highly conserved and suggests that its helical cap affects signaling kinetics in other PAS domains.

Mapping solvation dynamics at the function site of flavodoxin in three redox states

Flavoproteins are unique redox coenzymes, and the dynamic solvation at their function sites is critical to the understanding of their electron-transfer properties. Here, we report our complete characterization of the function-site solvation of holoflavodoxin in three redox states and of the binding-site solvation of apoflavodoxin. Using intrinsic flavin cofactor and tryptophan residue as the local optical probes with two site-specific mutations, we observed distinct ultrafast solvation dynamics at the function site in the three states and at the related recognition site of the cofactor, ranging from a few to hundreds of picoseconds. The initial ultrafast motion in 1-2.6 ps reflects the local water-network relaxation around the shallow, solvent-exposed function site. The second relaxation in 20-40 ps results from the coupled local water-protein fluctuation. The third dynamics in hundreds of picoseconds is from the intrinsic fluctuation of the loose loops flanking the cofactor at the function site. These solvation dynamics with different amplitudes well correlate with the redox states from the oxidized form, to the more rigid semiquinone and to the much looser hydroquinone. This observation of the redox control of local protein conformation plasticity and water network flexibility is significant, and such an intimate relationship is essential to the biological function of interprotein electron transfer.