Anticipatory active-site motions and chromophore distortion prime photoreceptor PYP for light activation

Spectroscopic Validation of Crystallographic Structures of a Protein Active Site by Chiroptical Spectroscopy

Out-of-plane distortions of a cofactor molecule in a protein active site are functionally important, and in photoreceptors, it has been proposed that they are crucial for spectral tuning and energy storage in photocycle intermediates. However, these subtle structural features are often beyond the grasp of structural biology. This issue is strikingly exemplified by photoactive yellow protein: its 14 independently determined crystal structures exhibit considerable differences in the dihedral angles defining the chromophore geometry, even though most of these are at excellent resolution. Here we developed a strategy to verify cofactor distortions in crystal structures by using quantum chemical calculations and chiroptical spectroscopy, particularly Raman optical activity and electronic circular dichroism spectroscopies. Based on this approach, we identify seven crystal structures with the chromophore geometries inconsistent with the experimentally observed data. The strategy implemented here promises to be widely applicable to uncovering cofactor distortions at active sites and to studies of reaction intermediates.

Photoactive Yellow Protein Represents a Distinct, Evolutionarily Novel Family of PAS Domains

Photoactive yellow protein (PYP) is a model photoreceptor. It binds a p-coumaric acid as a chromophore, thus enabling blue light sensing. The first discovered single-domain PYP from Halorhodospira halophila has been studied thoroughly in terms of its structural dynamics and photochemical properties. However, the evolutionary origins and biological role of PYP homologs are not well understood. Here, we show that PYP is an evolutionarily novel domain family of the ubiquitous PAS (Per-Arnt-Sim) superfamily. It likely originated from the phylum Myxococcota and was then horizontally transferred to representatives of a few other bacterial phyla. We show that PYP is associated with signal transduction either by domain fusion or by genome context. Key cellular functions modulated by PYP-initiated signal transduction pathways likely involve gene expression, motility, and biofilm formation. We identified three clades of the PYP family, one of which is poorly understood and potentially has novel functional properties. The Tyr42, Glu46, and Cys69 residues that are involved in p-coumaric acid binding in the model PYP from H. halophila are well conserved in the PYP family. However, we also identified cases where substitutions in these residues might have led to neofunctionalization, such as the proposed transition from light to redox sensing. Overall, this study provides definition, a newly built hidden Markov model, and the current genomic landscape of the PYP family and sets the stage for the future exploration of its signaling mechanisms and functional diversity. IMPORTANCE Photoactive yellow protein is a model bacterial photoreceptor. For many years, it was considered a prototypical model of the ubiquitous PAS domain superfamily. Here, we show that, in fact, the PYP family is evolutionarily novel, restricted to a few bacterial phyla and distinct from other PAS domains. We also reveal the diversity of PYP-containing signal transduction proteins and their potential mechanisms.

Photo-isomerization of the isolated photoactive yellow protein chromophore: what comes before the primary step?

Photoactive proteins typically rely on structural changes in a small chromophore to initiate a biological response. While these changes often involve isomerization as the "primary step", preceding this is an ultrafast relaxation of the molecular framework caused by the sudden change in electronic structure upon photoexcitation. Here, we capture this motion for an isolated model chromophore of the photoactive yellow protein using time-resolved photoelectron imaging. It occurs in <150 fs and is apparent from a spectral shift of ∼70 meV and a change in photoelectron anisotropy. Electronic structure calculations enable the quantitative assignment of the geometric and electronic structure changes to a planar intermediate from which the primary step can then proceed.

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.

NanoFAST: structure-based design of a small fluorogen-activating protein with only 98 amino acids

One of the essential characteristics of any tag used in bioscience and medical applications is its size. The larger the label, the more it may affect the studied object, and the more it may distort its behavior. In this paper, using NMR spectroscopy and X-ray crystallography, we have studied the structure of fluorogen-activating protein FAST both in the apo form and in complex with the fluorogen. We showed that significant change in the protein occurs upon interaction with the ligand. While the protein is completely ordered in the complex, its apo form is characterized by higher mobility and disordering of its N-terminus. We used structural information to design the shortened FAST (which we named nanoFAST) by truncating 26 N-terminal residues. Thus, we created the shortest genetically encoded tag among all known fluorescent and fluorogen-activating proteins, which is composed of only 98 amino acids.

Protonation Equilibrium in the Active Site of the Photoactive Yellow Protein

The role and existence of low-barrier hydrogen bonds (LBHBs) in enzymatic and protein activity has been largely debated. An interesting case is that of the photoactive yellow protein (PYP). In this protein, two short HBs adjacent to the chromophore, p-coumaric acid (pCA), have been identified by X-ray and neutron diffraction experiments. However, there is a lack of agreement on the chemical nature of these H-bond interactions. Additionally, no consensus has been reached on the presence of LBHBs in the active site of the protein, despite various experimental and theoretical studies having been carried out to investigate this issue. In this work, we perform a computational study that combines classical and density functional theory (DFT)-based quantum mechanical/molecular mechanical (QM/MM) simulations to shed light onto this controversy. Furthermore, we aim to deepen our understanding of the chemical nature and dynamics of the protons involved in the two short hydrogen bonds that, in the dark state of PYP, connect pCA with the two binding pocket residues (E46 and Y42). Our results support the existence of a strong LBHB between pCA and E46, with the H fully delocalized and shared between both the carboxylic oxygen of E46 and the phenolic oxygen of pCA. Additionally, our findings suggest that the pCA interaction with Y42 can be suitably described as a typical short ionic H-bond of moderate strength that is fully localized on the phenolic oxygen of Y42.

Structural and spectroscopic characterization of photoactive yellow protein and photoswitchable fluorescent protein constructs containing heavy atoms

Photo-induced structural rearrangements of chromophore-containing proteins are essential for various light-dependent signaling pathways and optogenetic applications. Ultrafast structural and spectroscopic methods have offered insights into these structural rearrangements across many timescales. However, questions still remain about exact mechanistic details, especially regarding photoisomerization of the chromophore within these proteins femtoseconds to picoseconds after photoexcitation. Instrumentation advancements for time-resolved crystallography and ultrafast electron diffraction provide a promising opportunity to study these reactions, but achieving enough signal-to-noise is a constant challenge. Here we present four new photoactive yellow protein constructs and one new fluorescent protein construct that contain heavy atoms either within or around the chromophore and can be expressed with high yields. Structural characterization of these constructs, most at atomic resolution, show minimal perturbation caused by the heavy atoms compared to wild-type structures. Spectroscopic studies report the effects of the heavy atom identity and location on the chromophore's photophysical properties. None of the substitutions prevent photoisomerization, although certain rates within the photocycle may be affected. Overall, these new proteins containing heavy atoms are ideal samples for state-of-theart time-resolved crystallography and electron diffraction experiments to elucidate crucial mechanistic information of photoisomerization.

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.

Local dynamics of the photo-switchable protein PYP in ground and signalling state probed by 2D-IR spectroscopy of –SCN labels

We employ 2D-IR spectroscopy of the protein label –SCN to describe the local dynamics in the photo-switchable protein PYP in its dark state (pG) and after photoactivation, concomitant with vast structural rearrangements, in its signalling state (pB).

Vibrational Lifetime of the SCN Protein Label in H2O and D2O Reports Site-Specific Solvation and Structure Changes During PYP's Photocycle

The application of vibrational labels such as thiocyanate  (-S-C≡N) for studying protein structure and dynamics is thriving. Absorption spectroscopy is usually employed to obtain wavenumber and line shape of the label. An observable of great significance might be the vibrational lifetime, which can be obtained by pump probe or 2D-IR spectroscopy. Due to the insulating effect of the heavy sulfur atom in the case of the SCN label, the lifetime of the C≡N oscillator is expected to be particularly sensitive to its surrounding as it is not dominated by through-bond relaxation. We therefore investigate the vibrational lifetime of the SCN label at various positions in the blue light sensor protein Photoactive Yellow Protein (PYP) in the ground state and signaling state of the photoreceptor. We find that the vibrational lifetime of the C≡N stretching mode is strongly affected both by its protein environment and by the degree of exposure to the solvent. Even for label positions where the line shape and wavenumber observed by FTIR are barely changing upon activation of the photoreceptor, we find that the lifetime can change considerably. To obtain an unambiguous measure for the solvent exposure of the labeled site, we show that it is imperative to compare the lifetimes in H₂O and D₂O. Importantly, the lifetimes shorten in H₂O as compared to D₂O for water exposed labels, while they stay largely the same for buried labels. We quantify this effect by defining a solvent exclusion coefficient (SEC). The response of the label's vibrational lifetime to its solvent exposure renders it a suitable universal probe for protein investigations. This applies even to systems that are otherwise hard to address, such as transient or short-lived states, which could be created during a protein's working cycle (as here in PYP) or during protein folding. It is also applicable to flexible systems (intrinsically disordered proteins), protein-protein and protein-membrane interactions.

Visualization of H atoms in the X-ray crystal structure of photoactive yellow protein: Does it contain low-barrier hydrogen bonds?

The hydrogen bond (HB) between 4-hydroxycinnamic acid (HC4) and glutamic acid E46 of photoactive yellow protein is exceptionally strong. In the 0.82-å resolution X-ray structure for this protein (PDB ID: 1NWZ), the OH…O distance is only 2.57 å. The position of the H atom between these two O atoms has not been determined in that structure, and in the absence of that information, it is impossible to determine whether or not this HB is a low-barrier HB (LBHB), as was proposed recently based on neutron structures of this protein (Yamaguchi et al., Proceedings of the National Academy of Sciences of the United States of America, 2009, 106: 440-444). Residual electron density maps computed using the 1NWZ data reveal that this H atom is 0.92 å from the O_ε2 atom of E46 and 1.67 å from the O₄ ' of HC4, and that the OH…O bond angle is 167°. These observations indicate that E46 is protonated, and HC4 is deprotonated, as was originally suggested, and that the HB in question is not an LBHB.

Unraveling the structural and chemical features of biological short hydrogen bonds

Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.

The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.

Perturbation of Short Hydrogen Bonds in Photoactive Yellow Protein via Noncanonical Amino Acid Incorporation

Photoactive yellow protein (PYP) is a small photoreceptor protein that has two unusually short hydrogen bonds between the deprotonated p-coumaric acid chromophore and two amino acids, a tyrosine and a glutamic acid. This has led to considerable debate as to whether the glutamic acid-chromophore hydrogen bond is a low barrier hydrogen bond, with conflicting results in the literature. We have modified the p Ka of the tyrosine by amber suppression and of the chromophore by chemical substitution. X-ray crystal structures of these modified proteins are nearly identical to the wild-type protein, so the heavy atom distance between proton donor and acceptor is maintained, even though these modifications change the relative proton affinity between donor and acceptor. Despite a considerable change in relative proton affinity, the NMR chemical shifts of the hydrogen-bonded protons are only moderately affected. QM/MM calculations were used to explore the protons' potential energy surface and connect the calculated proton position with empirically measured proton chemical shifts. The results are inconsistent with a low barrier hydrogen bond but in all cases are consistent with a localized proton, suggesting an ionic hydrogen bond rather than a low barrier hydrogen bond.

Following local light-induced structure changes and dynamics of the photoreceptor PYP with the thiocyanate IR label

Photoactive Yellow Protein (PYP) is a bacterial blue light receptor that enters a photocycle after excitation. The intermediate states are formed on time scales ranging from femtoseconds up to hundreds of milliseconds, after which the signaling state with a lifetime of about 1 s is reached. To investigate structural changes and dynamics, we incorporated the SCN IR label at distinct positions of the photoreceptor via cysteine mutation and cyanylation. FT-IR measurements of the SCN label at different sites of the well-established dark state structure of PYP characterized the spectral response of the label to differences in the environment. Under constant blue light irradiation, we observed the formation of the signaling state with significant changes of wavenumber and lineshape of the SCN bands. Thereby we deduced light-induced structural changes in the local environment of the labels. These results were supported by molecular dynamics simulations on PYP providing the solvent accessible surface area (SASA) at the different positions. To follow protein dynamics via the SCN label during the photocycle, we performed step-scan FT-IR measurements with a time resolution of 10 μs. Global analysis yielded similar time constants of τ1 = 70 μs, τ2 = 640 μs, and τ3 > 20 ms for the wild type and τ1 = 36 μs, τ2 = 530 μs, and τ3 > 20 ms for the SCN-labeled mutant PYP-A44C*, a mutant which provided a sufficiently large SCN difference signal to measure step-scan FT-IR spectra. In comparison to the protein (amide, E46) and chromophore bands the dynamics of the SCN label show a different behavior. This result indicates that the local kinetics sensed by the label are different from the global protein kinetics.

Improved Chemical-Genetic Fluorescent Markers for Live Cell Microscopy

Inducible chemical-genetic fluorescent markers are promising tools for live cell imaging requiring high spatiotemporal resolution and low background fluorescence. The fluorescence-activating and absorption shifting tag (FAST) was recently developed to form fluorescent molecular complexes with a family of small, synthetic fluorogenic chromophores (so-called fluorogens). Here, we use rational design to modify the binding pocket of the protein and screen for improved fluorescence performances with four different fluorogens. The introduction of a single mutation results in improvements in both quantum yield and dissociation constant with nearly all fluorogens tested. Our improved FAST (iFAST) allowed the generation of a tandem iFAST (td-iFAST) that forms green and red fluorescent reporters 1.6-fold and 2-fold brighter than EGFP and mCherry, respectively, while having a comparable size.

Spectroscopic ruler for measuring active-site distortions based on Raman optical activity of a hydrogen out-of-plane vibration

Photoactive yellow protein (PYP), from the phototrophic bacterium Halorhodospira halophila, is a small water-soluble photoreceptor protein and contains p-coumaric acid (pCA) as a chromophore. PYP has been an attractive model for studying the physical chemistry of protein active sites. Here, we explore how Raman optical activity (ROA) can be used to extract quantitative information on distortions of the pCA chromophore at the active site in PYP. We use ¹³C8-pCA to assign an intense signal at 826 cm^-1 in the ROA spectrum of PYP to a hydrogen out-of-plane vibration of the ethylenic moiety of the chromophore. Quantum-chemical calculations based on density functional theory demonstrate that the sign of this ROA band reports the direction of the distortion in the dihedral angle about the ethylenic C=C bond, while its amplitude is proportional to the dihedral angle. These results document the ability of ROA to quantify structural deformations of a cofactor molecule embedded in a protein moiety.

High-throughput, image-based screening of pooled genetic-variant libraries

We report a high-throughput screening method that allows diverse genotypes and corresponding phenotypes to be imaged in individual cells. We achieve genotyping by introducing barcoded genetic variants into cells as pooled libraries and reading the barcodes out using massively multiplexed fluorescence in situ hybridization. To demonstrate the power of image-based pooled screening, we identified brighter and more photostable variants of the fluorescent protein YFAST among 60,000 variants.

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.

Molecular modeling of the AhR structure and interactions can shed light on ligand-dependent activation and transformation mechanisms

Molecular modeling has given important contributions to elucidation of the main stages in the AhR signal transduction pathway. Despite the lack of experimentally determined structures of the AhR functional domains, information derived from homologous systems has been exploited for modeling their structure and interactions. Homology models of the AhR PASB domain have provided information on the binding cavity and contributed to elucidate species-specific differences in ligand binding. Molecular Docking simulations of the ligand binding process have given insights into differences in binding of diverse agonists, antagonists, and selective AhR modulators, and their application to virtual screening of large databases of compounds have allowed identification of novel AhR ligands. Recently available structural information on protein-protein and protein-DNA complexes of other bHLH-PAS systems has opened the way for modeling the AhR:ARNT dimer structure and investigating the mechanisms of AhR transformation and DNA binding. Future research directions should include simulation of the protein dynamics to obtain a more reliable description of intermolecular interactions involved in signal transmission.

On the reliability of peptide nonplanarity seen in ultra-high resolution crystal structures

Ultra-high resolution protein crystal structures have been considered as relatively reliable sources for defining details of protein geometry, such as the extent to which the peptide unit deviates from planarity. Chellapa and Rose (Proteins 2015; 83:1687) recently called this into question, reporting that for a dozen representative protein structures determined at ∼ 1 Å resolution, the diffraction data could be equally well fit with models restrained to have highly planar peptides, i.e. having a standard deviation of the ω torsion angles of only ∼ 1° instead of the typically observed value of ∼ 6°. Here, we document both conceptual and practical shortcomings of that study and show that the more tightly restrained models are demonstrably incorrect and do not fit the diffraction data equally well. We emphasize the importance of inspecting electron density maps when investigating the agreement between a model and its experimental data. Overall, this report reinforces that modern standard refinement protocols have been well-conceived and that ultra-high resolution protein crystal structures, when evaluated carefully and used with an awareness of their levels of coordinate uncertainty, are powerful sources of information for providing reliable information about the details of protein geometry.

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.

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.

Photo-control of DNA binding by an engrailed homeodomain-photoactive yellow protein hybrid

A photo-controlled version of the engrailed homeodomain (zENG) was created by inserting the homeodomain into a surface loop of a circularly permuted version of the photoactive yellow protein (cPYP). The two proteins fold independently as judged by NMR and fluorescence denaturation measurements. In the dark, the affinity of the zENG domain for its cognate DNA is inhibited >100-fold compared to wild-type zENG. Blue-light irradiation of the hybrid protein leads to enhanced conformational dynamics of the cPYP portion and a two-fold enhancement of the DNA binding affinity of the zENG domain. These results suggest that insertion into a surface loop of cPYP can be a general approach for conferring an initial level of photo-control on a given target protein. Focussed mutation/selection strategies may then be used to enhance the degree of photo-control.

Room temperature structures beyond 1.5 Å by serial femtosecond crystallography

About 2.5 × 10(6) snapshots on microcrystals of photoactive yellow protein (PYP) from a recent serial femtosecond crystallographic (SFX) experiment were reanalyzed to maximum resolution. The resolution is pushed to 1.46 Å, and a PYP structural model is refined at that resolution. The result is compared to other PYP models determined at atomic resolution around 1 Å and better at the synchrotron. By comparing subtleties such as individual isotropic temperature factors and hydrogen bond lengths, we were able to assess the quality of the SFX data at that resolution. We also show that the determination of anisotropic temperature factor ellipsoids starts to become feasible with the SFX data at resolutions better than 1.5 Å.

Origins of the Intermediate Spectral Form in M100 Mutants of Photoactive Yellow Protein

Numerous single-site mutants of photoactive yellow protein (PYP) from Halorhodospira halophila and as well as PYP homologs from other species exhibit a shoulder on the short wavelength side of the absorbance maximum in their dark-adapted states. The structural basis for the occurrence of this shoulder, called the "intermediate spectral form," has only been investigated in detail for the Y42F mutation. Here we explore the structural basis for occurrence of the intermediate spectral form in a M121E derivative of a circularly permuted H. halophila PYP (M121E-cPYP). The M121 site in M121E-cPYP corresponds to the M100 site in wild-type H. halophila PYP. High-resolution NMR measurements with a salt-tolerant cryoprobe enabled identification of those residues directly affected by increasing concentrations of ammonium chloride, a salt that greatly enhances the fraction of the intermediate spectra form. Residues in the surface loop containing the M121E (M100E) mutation were found to be affected by ammonium chloride as well as a discrete set of residues that link this surface loop to the buried hydroxyl group of the chromophore via a hydrogen bond network. Localized changes in the conformational dynamics of a surface loop can thereby produce structural rearrangements near the buried hydroxyl group chromophore while leaving the large majority of residues in the protein unaffected.

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.

Are there really low-barrier hydrogen bonds in proteins? The case of photoactive yellow protein

For a long time, low-barrier hydrogen bonds (LBHBs) have been proposed to exist in many enzymes and to play an important role in their catalytic function, but the proof of their existence has been elusive. The transient formation of an LBHB in a protein system has been detected for the first time using neutron diffraction techniques on a photoactive yellow protein (PYP) crystal in a study published in 2009 (Yamaguchi, S.; et al. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 440-444). However, very recent theoretical studies based on electronic structure calculations and NMR resonance experiments on PYP in solution (Saito, K.; et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 167-172) strongly indicate that there is not such an LBHB. By means of electronic structure calculations combined with the solution of the nuclear Schrödinger equation, we analyze here under which conditions an LBHB can exist in PYP, thus leading to a more reasonable and conciliating understanding of the above-mentioned studies.

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.

Simulations of Two-dimensional Infrared and Stimulated Resonance Raman Spectra of Photoactive Yellow Protein

We present simulations of one and two-dimensional infrared (2DIR) and stimulated resonance Raman (SRR) spectra of the dark state (pG) and early red-shifted intermediate (pR) of photoactive yellow protein (PYP). Shifts in the amide I and Glu46 COOH stretching bands distinguish between pG and pR in the IR absorption and 2DIR spectra. The one-dimensional SRR spectra are similar to the spontaneous RR spectra. The two-dimensional SRR spectra show large changes in cross peaks involving the C=O stretch of the two species and are more sensitive to the chromophore structure than 2DIR spectra.

A circularly permuted photoactive yellow protein as a scaffold for photoswitch design

Upon blue light irradiation, photoactive yellow protein (PYP) undergoes a conformational change that involves large movements at the N-terminus of the protein. We reasoned that this conformational change might be used to control other protein or peptide sequences if these were introduced as linkers connecting the N- and C-termini of PYP in a circular permutant. For such a design strategy to succeed, the circularly permuted PYP (cPYP) would have to fold normally and undergo a photocycle similar to that of the wild-type protein. We created a test cPYP by connecting the N- and C-termini of wild-type PYP (wtPYP) with a GGSGGSGG linker polypeptide and introducing new N- and C-termini at G115 and S114, respectively. Biophysical analysis indicated that this cPYP adopts a dark-state conformation much like wtPYP and undergoes wtPYP-like photoisomerization driven by blue light. However, thermal recovery of dark-state cPYP is ∼10-fold faster than that of wtPYP, so that very bright light is required to significantly populate the light state. Targeted mutations at M121E (M100 in wtPYP numbering) were found to enhance the light sensitivity substantially by lengthening the lifetime of the light state to ∼10 min. Nuclear magnetic resonance (NMR), circular dichroism, and UV-vis analysis indicated that the M121E-cPYP mutant also adopts a dark-state structure like that of wtPYP, although protonated and deprotonated forms of the chromophore coexist, giving rise to a shoulder near 380 nm in the UV-vis absorption spectrum. Fluorine NMR studies with fluorotryptophan-labeled M121E-cPYP show that blue light drives large changes in conformational dynamics and leads to solvent exposure of Trp7 (Trp119 in wtPYP numbering), consistent with substantial rearrangement of the N-terminal cap structure. M121E-cPYP thus provides a scaffold that may allow a wider range of photoswitchable protein designs via replacement of the linker polypeptide with a target protein or peptide sequence.

Raman Optical Activity Probing Structural Deformations of the 4-Hydroxycinnamyl Chromophore in Photoactive Yellow Protein

Many biological cofactors, such as light-absorbing chromophores in photoreceptors, contain a π-electron system and are planar molecules. These cofactors are, however, usually nonplanar within a protein environment, and such structural distortions have been shown to be functionally important. Because the nonplanar structure makes the molecule chiral, Raman optical activity (ROA) provides a wealth of stereochemical information about the structural and conformational details of cofactors. The present study applied a near-infrared excited ROA to photoactive yellow protein, a blue light receptor. We successfully obtained the ROA spectra of the 4-hydroxycinnamyl chromophore embedded in a protein environment. Furthermore, calculations of the ROA spectra utilizing density functional theory provide detailed structural information, such as data on out-of-plane distortions of the chromophore. The structural information obtained from the ROA spectra includes the positions of hydrogen atoms, which are usually not detected in the crystal structures of biological samples.

Detection of alternative conformations:Shift_plotandAC_predictionprograms

Increased mobility of poorly ordered atoms compared to well ordered ones in unrestrained refinement was adopted as the basis for the development of a method for predicting alternative conformations. Two programs for the analysis of atomic mobility are presented.Shift_plotmay be useful for visual analysis of magnitudes of atomic shifts in refinement.AC_predictionimplements automatic decision-making procedures based on atomic shift magnitudes in unrestrained refinement to highlight the specific residues that are likely to harbour alternative conformations and that hence require checkingviaelectron density maps.

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.

Formation of an unusually short hydrogen bond in photoactive yellow protein

The photoactive chromophore of photoactive yellow protein (PYP) is p-coumaric acid (pCA). In the ground state, the pCA chromophore exists as a phenolate anion, which is H-bonded by protonated Glu46 (O(Glu46)-O(pCA)=~2.6Å) and protonated Tyr42. On the other hand, the O(Glu46)-O(pCA) H-bond was unusually short (O(Glu46)-O(pCA)=2.47Å) in the intermediate pR(CW) state observed in time-resolved Laue diffraction studies. To understand how the existence of the unusually short H-bond is energetically possible, we analyzed the H-bond energetics adopting a quantum mechanical/molecular mechanical (QM/MM) approach based on the atomic coordinates of the PYP crystal structures. In QM/MM calculations, the O(Glu46)-O(pCA) bond is 2.60Å in the ground state, where Tyr42 donates an H-bond to pCA. In contrast, when the hydroxyl group of Tyr42 is flipped away from pCA, the H-bond was significantly shortened to 2.49Å in the ground state. The same H-bond pattern reproduced the unusually short H-bond in the pR(CW) structure (O(Glu46)-O(pCA)=2.49Å). Intriguingly, the potential-energy profile resembles that of a single-well H-bond, suggesting that the pK(a) values of the donor (Glu46) and acceptor (pCA) moieties are nearly equal. The present results indicate that the "equal pK(a)" requirement for formation of single-well or low-barrier H-bond (LBHB) is satisfied only when Tyr42 does not donate an H-bond to pCA, and argue against the possibility that the O(Glu46)-O(pCA) bond is an LBHB in the ground state, where Tyr42 donates an H-bond to pCA.

Detection of alternative conformations by unrestrained refinement

Unrestrained refinement is stable for the vast majority of atoms when working at atomic resolution. Nevertheless, geometrical restraints should be retained in refinement for residues that are present in several (alternative) conformations in the crystal used for the X-ray experiment; otherwise, such residues deteriorate significantly. The authors believe that a large distortion of a residue in unrestrained refinement may hint at the presence of alternative conformations of this residue. To obtain these hints in a routine way, two methods of analyzing the shifts of atomic centres resulting from several cycles of unrestrained refinement are described. A simple diagram plotting the values of the atomic shifts against the residue number may give an idea of the crystallographic order of different parts of the structure at a qualitative level. To put the analysis on a more quantitative basis, several decision-making procedures were developed and tested which compose a list of residues that are likely to be present in alternative conformations or to be disordered and so should be checked thoroughly using Fourier syntheses and included in the model with alternative conformations when necessary. The parameters and performance of the suggested procedures were estimated by the use of 203 PDB structures refined at resolutions better than 1.2 Å. Decision-making procedures based on analysis of atomic shifts were found to be more reliable than similar procedures based on atomic displacement parameters or density values calculated at atomic centres.

Protein conformational switches: from nature to design

Protein conformational switches alter their shape upon receiving an input signal, such as ligand binding, chemical modification, or change in environment. The apparent simplicity of this transformation--which can be carried out by a molecule as small as a thousand atoms or so--belies its critical importance to the life of the cell as well as its capacity for engineering by humans. In the realm of molecular switches, proteins are unique because they are capable of performing a variety of biological functions. Switchable proteins are therefore of high interest to the fields of biology, biotechnology, and medicine. These molecules are beginning to be exploited as the core machinery behind a new generation of biosensors, functionally regulated enzymes, and "smart" biomaterials that react to their surroundings. As inspirations for these designs, researchers continue to analyze existing examples of allosteric proteins. Recent years have also witnessed the development of new methodologies for introducing conformational change into proteins that previously had none. Herein we review examples of both natural and engineered protein switches in the context of four basic modes of conformational change: rigid-body domain movement, limited structural rearrangement, global fold switching, and folding-unfolding. Our purpose is to highlight examples that can potentially serve as platforms for the design of custom switches. Accordingly, we focus on inducible conformational changes that are substantial enough to produce a functional response (e.g., in a second protein to which it is fused), yet are relatively simple, structurally well-characterized, and amenable to protein engineering efforts.

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.

Energetics of short hydrogen bonds in photoactive yellow protein

Recent neutron diffraction studies of photoactive yellow protein (PYP) proposed that the H bond between protonated Glu46 and the chromophore [ionized p-coumaric acid (pCA)] was a low-barrier H bond (LBHB). Using the atomic coordinates of the high-resolution crystal structure, we analyzed the energetics of the short H bond by two independent methods: electrostatic pK(a) calculations and a quantum mechanical/molecular mechanical (QM/MM) approach. (i) In the QM/MM optimized geometry, we reproduced the two short H-bond distances of the crystal structure: Tyr42-pCA (2.50 Å) and Glu46-pCA (2.57 Å). However, the H atoms obviously belonged to the Tyr or Glu moieties, and were not near the midpoint of the donor and acceptor atoms. (ii) The potential-energy curves of the two H bonds resembled those of standard asymmetric double-well potentials, which differ from those of LBHB. (iii) The calculated pK(a) values for Glu46 and pCA were 8.6 and 5.4, respectively. The pK(a) difference was unlikely to satisfy the prerequisite for LBHB. (iv) The LBHB in PYP was originally proposed to stabilize the ionized pCA because deprotonated Arg52 cannot stabilize it. However, the calculated pK(a) of Arg52 and QM/MM optimized geometry suggested that Arg52 was protonated on the protein surface. The short H bond between Glu46 and ionized pCA in the PYP ground state could be simply explained by electrostatic stabilization without invoking LBHB.