Water structural changes involved in the activation process of photoactive yellow protein

Potential Second-Harmonic Ghost Bands in Fourier Transform Infrared (FT-IR) Difference Spectroscopy of Proteins

Fourier transform infrared (FT-IR) difference absorption spectroscopy is a common method for studying the structural and dynamical aspects behind protein function. In particular, the 2800-1800 cm^-1 spectral range has been used to obtain information about internal (deuterated) water molecules, as well as site-specific details about cysteine residues and chemically modified and artificial amino acids. Here, we report on the presence of ghost bands in cryogenic light-induced FT-IR difference spectra of the protein bacteriorhodopsin. The presence of these ghost bands can be particularly problematic in the 2800-1900 cm^-1 region, showing intensities similar to O-D vibrations from water molecules. We demonstrate that they arise from second harmonics from genuine chromophore bands located in the 1400-850 cm^-1 region, generated by double-modulation artifacts caused from reflections of the IR beam at the sample and at the cryostat windows back to the interferometer (inter-reflections). The second-harmonic ghost bands can be physically removed by placing an optical filter of suitable cutoff in the beam path, but at the cost of losing part of the multiplexing advantage of FT-IR spectroscopy. We explored alternatives to the use of optical filters. Tilting the cryostat windows was effective in reducing the intensity of the second harmonic artifacts but tilting the sample windows was not, presumably by their close proximity to the focal point of the IR beam. We also introduce a simple numerical post-processing approach that can partially, but not fully, correct for second-harmonic ghost bands in FT-IR difference spectra.

Changes in the hydrogen-bonding strength of internal water molecules and cysteine residues in the conductive state of channelrhodopsin-1

Water plays an essential role in the structure and function of proteins, particularly in the less understood class of membrane proteins. As the first of its kind, channelrhodopsin is a light-gated cation channel and paved the way for the new and vibrant field of optogenetics, where nerve cells are activated by light. Still, the molecular mechanism of channelrhodopsin is not understood. Here, we applied time-resolved FT-IR difference spectroscopy to channelrhodopsin-1 from Chlamydomonas augustae. It is shown that the (conductive) P2(380) intermediate decays with τ ≈ 40 ms and 200 ms after pulsed excitation. The vibrational changes between the closed and the conductive states were analyzed in the X-H stretching region (X = O, S, N), comprising vibrational changes of water molecules, sulfhydryl groups of cysteine side chains and changes of the amide A of the protein backbone. The O-H stretching vibrations of "dangling" water molecules were detected in two different states of the protein using H2 (18)O exchange. Uncoupling experiments with a 1:1 mixture of H2O:D2O provided the natural uncoupled frequencies of the four O-H (and O-D) stretches of these water molecules, each with a very weakly hydrogen-bonded O-H group (3639 and 3628 cm(-1)) and with the other O-H group medium (3440 cm(-1)) to moderately strongly (3300 cm(-1)) hydrogen-bonded. Changes in amide A and thiol vibrations report on global and local changes, respectively, associated with the formation of the conductive state. Future studies will aim at assigning the respective cysteine group(s) and at localizing the "dangling" water molecules within the protein, providing a better understanding of their functional relevance in CaChR1.

Multiscale approach to the determination of the photoactive yellow protein signaling state ensemble

The nature of the optical cycle of photoactive yellow protein (PYP) makes its elucidation challenging for both experiment and theory. The long transition times render conventional simulation methods ineffective, and yet the short signaling-state lifetime makes experimental data difficult to obtain and interpret. Here, through an innovative combination of computational methods, a prediction and analysis of the biological signaling state of PYP is presented. Coarse-grained modeling and locally scaled diffusion map are first used to obtain a rough bird's-eye view of the free energy landscape of photo-activated PYP. Then all-atom reconstruction, followed by an enhanced sampling scheme; diffusion map-directed-molecular dynamics are used to focus in on the signaling-state region of configuration space and obtain an ensemble of signaling state structures. To the best of our knowledge, this is the first time an all-atom reconstruction from a coarse grained model has been performed in a relatively unexplored region of molecular configuration space. We compare our signaling state prediction with previous computational and more recent experimental results, and the comparison is favorable, which validates the method presented. This approach provides additional insight to understand the PYP photo cycle, and can be applied to other systems for which more direct methods are impractical.

Many protein systems of biological interest undergo dynamical changes on a time scale too long to be modeled using standard computational methods. One example is photoactive yellow protein (PYP), found in several bacterial species. Blue light, potentially harmful for DNA, triggers several structural changes in PYP, eventually resulting in a conformation that changes the swimming behavior of bacteria. This conformation is difficult to investigate, as it is too short lived. In addition, understanding this “signaling state” is computationally difficult because of the long timescale of the transition. We overcome this by constructing a coarse-grained model to rapidly induce transitions to the signaling state. We then reconstruct and further sample the all-atom configurations from these coarse-grained representations. Our results are consistent with all available experimental and computational evidence.

Proton transfer reactions and hydrogen-bond networks in protein environments

In protein environments, proton transfer reactions occur along polar or charged residues and isolated water molecules. These species consist of H-bond networks that serve as proton transfer pathways; therefore, thorough understanding of H-bond energetics is essential when investigating proton transfer reactions in protein environments. When the pKa values (or proton affinity) of the H-bond donor and acceptor moieties are equal, significantly short, symmetric H-bonds can be formed between the two, and proton transfer reactions can occur in an efficient manner. However, such short, symmetric H-bonds are not necessarily stable when they are situated near the protein bulk surface, because the condition of matching pKa values is opposite to that required for the formation of strong salt bridges, which play a key role in protein-protein interactions. To satisfy the pKa matching condition and allow for proton transfer reactions, proteins often adjust the pKa via electron transfer reactions or H-bond pattern changes. In particular, when a symmetric H-bond is formed near the protein bulk surface as a result of one of these phenomena, its instability often results in breakage, leading to large changes in protein conformation.

Effects of dehydration on light-induced conformational changes in bacterial photosynthetic reaction centers probed by optical and differential FTIR spectroscopy

Following light-induced electron transfer between the primary donor (P) and quinone acceptor (Q(A)) the bacterial photosynthetic reaction center (RC) undergoes conformational relaxations which stabilize the primary charge separated state P(+)Q(A)(-). Dehydration of RCs from Rhodobacter sphaeroides hinders these conformational dynamics, leading to acceleration of P(+)Q(A)(-) recombination kinetics [Malferrari et al., J. Phys. Chem. B 115 (2011) 14732-14750]. To clarify the structural basis of the conformational relaxations and the involvement of bound water molecules, we analyzed light-induced P(+)Q(A)(-)/PQ(A) difference FTIR spectra of RC films at two hydration levels (relative humidity r=76% and r=11%). Dehydration reduced the amplitude of bands in the 3700-3550cm(-1) region, attributed to water molecules hydrogen bonded to the RC, previously proposed to stabilize the charge separation by dielectric screening [Iwata et al., Biochemistry 48 (2009) 1220-1229]. Other features of the FTIR difference spectrum were affected by partial depletion of the hydration shell (r=11%), including contributions from modes of P (9-keto groups), and from NH or OH stretching modes of amino acidic residues, absorbing in the 3550-3150cm(-1) range, a region so far not examined in detail for bacterial RCs. To probe in parallel the effects of dehydration on the RC conformational relaxations, we analyzed by optical absorption spectroscopy the kinetics of P(+)Q(A)(-) recombination following the same photoexcitation used in FTIR measurements (20s continuous illumination). The results suggest a correlation between the observed FTIR spectral changes and the conformational rearrangements which, in the hydrated system, strongly stabilize the P(+)Q(A)(-) charge separated state over the second time scale.

H atom positions and nuclear magnetic resonance chemical shifts of short H bonds in photoactive yellow protein

Recent neutron diffraction studies on photoactive yellow protein (PYP) proposed that the H bond between protonated Glu46 and the chromophore-ionized p-coumaric acid (pCA) is a low-barrier H bond (LBHB) mainly because the H atom position was assigned at the midpoint of the O(Glu46)-O(pCA) bond. However, the (1)H nuclear magnetic resonance (NMR) chemical shift (δ(H)) was 15.2 ppm, which is lower than the values of 17-19 ppm for typical LBHBs. We evaluated the dependence of δ(H) on an H atom position in the O(Glu46)-O(pCA) bond in the PYP ground state by using a quantum mechanical/molecular mechanical (QM/MM) approach. The calculated chemical shift unambiguously suggested that a δ(H) of 15.2 ppm for the O(Glu46)-O(pCA) bond in NMR studies should correspond to the QM/MM geometry (δ(H) = 14.5 ppm), where the H atom belongs to the Glu moiety, rather than the neutron diffraction geometry (δ(H) = 19.7 ppm), where the H atom is near the midpoint of the donor and acceptor atoms.

Conformational heterogeneity and propagation of structural changes in the LOV2/Jalpha domain from Avena sativa phototropin 1 as recorded by temperature-dependent FTIR spectroscopy

Phototropins control phototropism, chloroplast movement, stomatal opening, and leaf expansion in plants. Phototropin 1 (phot1) is composed of a kinase domain linked to two blue light-sensing domains, LOV2 and LOV1, which bind flavin mononucleotide. Disruption of the interaction between the LOV2 domain and a helical segment named Jalpha, joining LOV to the kinase domain, induces the subsequent kinase activity of phototropin 1 and further-downstream signal transduction. Here we study the effects of temperature and hydration on the light-triggered signal propagation in the phot1 LOV2 domain of Avena sativa (AsLOV2/Jalpha), using Fourier transform infrared spectroscopy to unravel part of the molecular mechanism of phototropin 1. We report that AsLOV2/Jalpha shows an intense signal in the amide I and II regions, arising mainly from beta-sheet changes and the unbinding of the Jalpha helix from the Per-ARNT-Sim core and its subsequent partial unfolding. Importantly, these structural changes only occur under conditions of full hydration and at temperatures above 280 K. We characterized a newly isolated low-hydration intermediate that shows a downshift of high-frequency amide I signals and that possibly corresponds to loop tightening, without large beta-sheet or Jalpha structural changes. In addition, we report a heterogeneity in AsLOV2/Jalpha involving two different C(4)=O conformer populations, coexisting in the dark state and characterized by C(4)=O carbonyl frequencies at 1712 cm(-1) and 1694 cm(-1) that are attributable to a single H-bond and two H-bonds at this site, respectively. Such conformers display slightly shifted absorption spectra and cause a splitting of the 475-nm band in the ultraviolet/visible spectra of LOV domains at low temperature.

Hydration-Induced Far-Infrared Absorption Increase in Myoglobin

The interaction of proteins with an aqueous environment leads to a thin region of "biological water", the molecules of which have properties that differ from those of bulk water, in particular, reduced absorption of far-infrared radiation caused by protein-induced hindrance of the water rotational and vibrational degrees of freedom. New results at terahertz (THz) frequencies, however, show that absorption per protein molecule is increased by the presence of biological water. Absorption measurements were made of the heme protein myoglobin mixed with water from 3.6 to 98 wt % in the frequency range of 0.1-1.2 THz, using THz time-domain spectroscopy. Analysis shows greater THz absorption when compared to a non-interacting protein-water model. Including the suppressed absorption of biological water leads to a substantial hydration-dependent increase in absorption per protein molecule over a wide range of concentration and frequencies, meaning that water increases the protein's polarizability.

Structural changes during the photocycle of photoactive yellow protein monitored by ultraviolet resonance raman spectra of tyrosine and tryptophan

Photoactive yellow protein (PYP) is a bacterial blue light photoreceptor, and photoexcitation of dark-state PYP (PYP(dark)) triggers a photocycle that involves several intermediate states. We report the ultraviolet resonance Raman spectra of PYP with 225-250 nm excitations and investigate protein structural changes accompanying the formation of the putative signaling state denoted PYP(M). The PYP(M)-PYP(dark) difference spectra show several features of tyrosine and tryptophan, indicating environmental changes for these amino acid residues. The tyrosine difference signals show small upshifts with intensity changes in Y8a and Y9a bands. Although there are five tyrosine residues in PYP, Tyr42 and Tyr118 are suggested to be responsible for the difference signals on the basis of a global fitting analysis of the difference spectra at different excitation wavelengths and the crystal structure of PYP(dark). A further experiment on the Thr50-->Val mutant supports environmental changes in Tyr42. The observed upshift of the Y8a band suggests a weaker or broken hydrogen bond between Tyr42 and the chromophore in PYP(M). In addition, a reorientation of the OH group in Tyr42 is suggested from the upshift of the Y9a band. For tryptophan, the Raman bands of W3, W16, and W18 modes diminish in intensity upon formation of PYP(M). The loss of intensities is attributable to an exposure of tryptophan in PYP(M). PYP contains only one tryptophan (Trp119) that is located more than 10 A from the active site. Thus the observed changes are indicative of global conformational changes in protein during the transition from PYP(dark) to PYP(M). These results are in line with the currently proposed photocycle mechanism of PYP.

pH-dependent Equilibrium between Long Lived Near-UV Intermediates of Photoactive Yellow Protein

The long lived intermediate (signaling state) of photoactive yellow protein (PYP(M)), which is formed in the photocycle, was characterized at various pHs. PYP(M) at neutral pH was in equilibrium between two spectroscopically distinct states. Absorption maxima of the acidic form (PYP(M)(acid)) and alkaline form (PYP(M)(alkali)) were located at 367 and 356 nm, respectively. Equilibrium was represented by the Henderson-Hasselbalch equation, in which apparent pK(a) was 6.4. Content of alpha- and/or beta-structure of PYP(M)(acid) was significantly greater than PYP(M)(alkali) as demonstrated by the molar ellipticity at 222 nm. In addition, changes in amide I and II modes of beta-structure in the difference Fourier transform infrared spectra for formation of PYP(M)(acid) was smaller than that of PYP(M)(alkali). The vibrational mode at 1747 cm(-1) of protonated Glu-46 was found as a small band for PYP(M)(acid) but not for PYP(M)(alkali), suggesting that Glu-46 remains partially protonated in PYP(M)(acid), whereas it is fully deprotonated in PYP(M)(alkali). Small angle x-ray scattering measurements demonstrated that the radius of gyration of PYP(M)(acid) was 15.7 Angstroms, whereas for PYP(M)(alkali) it was 16.2 Angstroms. These results indicate that PYP(M)(acid) assumes a more ordered and compact structure than PYP(M)(alkali). Binding of citrate shifts this equilibrium toward PYP(M)(alkali). UV-visible absorption spectra and difference infrared spectra of the long lived intermediate formed from E46Q mutant was consistent with those of PYP(M)(acid), indicating that the mutation shifts this equilibrium toward PYP(M)(acid). Alterations in the nature of PYP(M) by pH, citrate, and mutation of Glu-46 are consistently explained by the shift of the equilibrium between PYP(M)(acid) and PYP(M)(alkali).

Predicting the signaling state of photoactive yellow protein

As a bacterial blue light sensor the photoactive yellow protein (PYP) undergoes conformational changes upon signal transduction. The absorption of a photon triggers a series of events that are initially localized around the protein chromophore, extends to encompass the whole protein within microseconds, and leads to the formation of the transient pB signaling state. We study the formation of this signaling state pB by molecular simulation and predict its solution structure. Conventional straightforward molecular dynamics is not able to address this formation process due to the long (microsecond) timescales involved, which are (partially) caused by the presence of free energy barriers between the metastable states. To overcome these barriers, we employed the parallel tempering (or replica exchange) method, thus enabling us to predict qualitatively the formation of the PYP signaling state pB. In contrast to the receptor state pG of PYP, the characteristics of this predicted pB structure include a wide open chromophore-binding pocket, with the chromophore and Glu(46) fully solvent-exposed. In addition, loss of alpha-helical structure occurs, caused by the opening motion of the chromophore-binding pocket and the disruptive interaction of the negatively charged Glu(46) with the backbone atoms in the hydrophobic core of the N-terminal cap. Recent NMR experiments agree very well with these predictions.

Deuterium isotope effects in the photocycle transitions of the photoactive yellow protein

The Photoactive Yellow Protein (PYP) from Halorhodospira halophila (formerly Ectothiorhodospira halophila) is increasingly used as a model system. As such, a thorough understanding of the photocycle of PYP is essential. In this study we have combined information from pOH- (or pH-) dependence and (kinetic) deuterium isotope effects to elaborate on existing photocycle models. For several characteristics of PYP we were able to make a distinction between pH- and pOH-dependence, a nontrivial distinction when comparing data from samples dissolved in H(2)O and D(2)O. It turns out that most characteristics of PYP are pOH-dependent. We confirmed the existence of a pB' intermediate in the pR to pB transition of the photocycle. In addition, we were able to show that the pR to pB' transition is reversible, which explains the previously observed biexponential character of the pR-to-pB photocycle step. Also, the absorption spectrum of pB' is slightly red-shifted with respect to pB. The recovery of the pG state is accompanied by an inverse kinetic deuterium isotope effect. Our interpretation of this is that before the chromophore can be isomerized, it is deprotonated by a hydroxide ion from solution. From this we propose a new photocycle intermediate, pB(deprot), from which pG is recovered and which is in equilibrium with pB. This is supported in our data through the combination of the observed pOH and pH dependence, together with the kinetic deuterium isotope effect.

Photoreaction of the cysteine S-H group in the LOV2 domain of Adiantum phytochrome3

Phy3-LOV2 is the chromophore domain of a blue-light photoreceptor for tropic responses of plants. FMN is noncovalently bound to phy3-LOV2, and the protein structure is characteristic of the LOV (light-oxygen-voltage) domain. Primary photoreaction is considered to be adduct formation between FMN and a cysteine, while deprotonation of the cysteine S-H group was suggested. On the basis of the infrared spectral analysis, however, we have shown that the cysteine of phy3-LOV2 is in the protonated S-H form, and not in the thiolate form in the ground state. Upon formation of S390, the S-H group disappears, presumably forming the adduct between FMN and Cys966. S390 can be trapped at 150 K, and the protein structure of S390 may not be changed drastically at 295 K.

Structural change of site-directed mutants of PYP: new dynamics during pR state

The energetics, protein dynamics, and diffusion coefficients of three mutants of photoactive yellow protein, R52Q, P68A, and W119G, were studied by the transient grating and pulsed laser-induced photoacoustic method. We observed a new dynamics with a lifetime of approximately 1 micro s in the transient grating signal, which is silent by the light absorption technique. This fact indicates that, after the structure change around the chromophore is completed (pR(1)), the protein part located far from the chromophore is still moving to finally create another pR (pR(2)) species, which can transform to the next intermediate, pB. Although the kinetics of pR(2)-->pB-->pG are very different depending on the mutants, the enthalpies of the first long-lived (in micro seconds, 100-micro s range) intermediate species (pR(2)) are similar and very high for all mutants. The diffusion coefficients of the parent (pG) and pB species of the mutants are also similar to that of the wild-type photoactive yellow protein. From the temperature dependence of the volume change, the difference in the thermal expansion coefficients taken as indicator of the flexibility of the structure between pG and pR(2) is measured. They are also similar to that of the wild-type photoactive yellow protein. These results suggest that the protein structures of pR(2) and pB in these mutants are globally different from that of pG, and this structural change is not altered so much by the single amino acid residue mutation. This is consistent with the partially unfolded nature of these intermediate species. On the other hand, the volume changes during pR(1)-->pR(2) are sensitive to the mutations, which may suggest that the volume change reflects a rather local character of the structure, such as the chromophore-protein interaction.

The role of the N-terminal domain of photoactive yellow protein in the transient partial unfolding during signalling state formation

It is shown that the N-terminal domain of photoactive yellow protein (PYP), which appears relatively independently folded in the ground state of the protein, plays a key role in the transient unfolding during signalling state formation: genetic truncation of the N-terminal domain of PYP significantly decreases the extent of cooperativity of the titration curve that describes chromophore protonation in the ground state of PYP, which is in agreement with the notion that the N-terminal domain is linked through a hydrogen-bonding network with the chromophore-containing domain of the protein. Furthermore, deletion of the N-terminal domain completely abolishes the non-linearity of the Arrhenius plot of the rate of ground state recovery.

Tight Asp-85--Thr-89 association during the pump switch of bacteriorhodopsin

Unidirectional proton transport in bacteriorhodopsin is enforced by the switching machinery of the active site. Threonine 89 is located in this region, with its O--H group forming a hydrogen bond with Asp-85, the acceptor for proton transfer from the Schiff base of the retinal chromophore. Previous IR spectroscopy of [3-(18)O]threonine-labeled bacteriorhodopsin showed that the hydrogen bond of the O--D group of Thr-89 in D(2)O is strengthened in the K photocycle intermediate. Here, we show that the strength and orientation of this hydrogen bond remains unchanged in the L intermediate and through the M intermediate. Furthermore, a strong interaction between Asp-85 and the O--H (O--D) group of Thr-89 in M is indicated by a shift in the C==O stretching vibration of the former because of (18)O substitution in the latter. Thus, the strong hydrogen bond between Asp-85 and Thr-89 in K persists through M, contrary to structural models based on x-ray crystallography of the photocycle intermediates. We propose that, upon photoisomerization of the chromophore, Thr-89 forms a tight, persistent complex with one of the side-chain oxygens of Asp-85 and is thereby precluded from participating in the switching process. On the other hand, the loss of hydrogen bonding at the other oxygen of Asp-85 in M may be related to the switching event.