Low-barrier hydrogen bond in photoactive yellow protein

Direct evidence for a deprotonated lysine serving as a H-bond "acceptor" in a photoreceptor protein

Deprotonation or suppression of the pKa of the amino group of a lysine sidechain is a widely recognized phenomenon whereby the sidechain amino group transiently can act as a nucleophile at the active site of enzymatic reactions. However, a deprotonated lysine and its molecular interactions have not been directly experimentally detected. Here, we demonstrate a deprotonated lysine stably serving as an "acceptor" in a H-bond between the photosensor protein RcaE and its chromophore. Signal splitting and trans-H-bond J coupling observed by NMR spectroscopy provide direct evidence that Lys261 is deprotonated and serves as a H-bond acceptor for the chromophore NH group. Quantum mechanical/molecular mechanical calculations also indicate that this H-bond exists stably. Interestingly, the sidechain amino group of the lysine can act as both donor and acceptor. The remarkable shift in the H-bond characteristics arises from a decrease in solvation, triggered by photoisomerization. Our results provide insights into the dual role of this lysine. This mechanism has broad implications for other biological reactions in which lysine plays a role.

Very Strong Hydrogen Bond in Nitrophthalic Cocrystals

This work presents the studies of a very strong hydrogen bond (VSHB) in biologically active phthalic acids. Research on VSHB comes topical due to its participation in many biological processes. The studies cover the modelling of intermolecular interactions and phthalic acids with 2,4,6-collidine and N,N-dimethyl-4-pyridinamine complexes with aim to obtain a VSHB. The four synthesized complexes were studied by experimental X-ray, IR, and Raman methods, as well as theoretical Car-Parrinello Molecular Dynamics (CP-MD) and Density Functional Theory (DFT) simulations. By variation of the steric repulsion and basicity of the complex' components, a very short intramolecular hydrogen bond was achieved. The potential energy curves calculated by the DFT method were characterized by a low barrier (0.7 and 0.9 kcal/mol) on proton transfer in the OHN intermolecular hydrogen bond for 3-nitrophthalic acid with either 2,4,6-collidine or N,N-dimethyl-4-pyridinamine cocrystals. Moreover, the CP-MD simulations exposed very strong bridging proton dynamics in the intermolecular hydrogen bonds. The accomplished crystallographic and spectroscopic studies indicate that the OHO intramolecular hydrogen bond in 4-nitrophthalic cocrystals is VSHB. The influence of a strong steric effect on the geometry of the studied cocrystals and the stretching vibration bands of the carboxyl and carboxylate groups was elaborated.

gmXtal: Cooking Crystals with GROMACS

Molecular dynamics (MD) simulations are routinely performed of biomolecules in solution, because this is their native environment. However, the structures used in such simulations are often obtained with X-ray crystallography, which provides the atomic coordinates of the biomolecule in a crystal environment. With the advent of free electron lasers and time-resolved techniques, X-ray crystallography can now also access metastable states that are intermediates in a biochemical process. Such experiments provide additional data, which can be used, for example, to optimize MD force fields. Doing so requires that the simulation of the biomolecule is also performed in the crystal environment. However, in contrast to simulations of biomolecules in solution, setting up a crystal is challenging. In particular, because not all solvent molecules are resolved in X-ray crystallography, adding a suitable number of solvent molecules, such that the properties of the crystallographic unit cell are preserved in the simulation, can be difficult and typically is a trial-and-error based procedure requiring manual interventions. Such interventions preclude high throughput applications. To overcome this bottleneck, we introduce gmXtal, a tool for setting up crystal simulations for MD simulations with GROMACS. With the information from the protein data bank (rcsb.org) gmXtal automatically (i) builds the crystallographic unit cell; (ii) sets the protonation of titratable residues; (iii) builds missing residues that were not resolved experimentally; and (iv) adds an appropriate number of solvent molecules to the system. gmXtal is available as a standalone tool https://gitlab.com/pbuslaev/gmxtal .

Watching a signaling protein function: What has been learned over four decades of time-resolved studies of photoactive yellow protein

Photoactive yellow protein (PYP) is a signaling protein whose internal p-coumaric acid chromophore undergoes reversible, light-induced trans-to-cis isomerization, which triggers a sequence of structural changes that ultimately lead to a signaling state. Since its discovery nearly 40 years ago, PYP has attracted much interest and has become one of the most extensively studied proteins found in nature. The method of time-resolved crystallography, pioneered by Keith Moffat, has successfully characterized intermediates in the PYP photocycle at near atomic resolution over 12 decades of time down to the sub-picosecond time scale, allowing one to stitch together a movie and literally watch a protein as it functions. But how close to reality is this movie? To address this question, results from numerous complementary time-resolved techniques including x-ray crystallography, x-ray scattering, and spectroscopy are discussed. Emerging from spectroscopic studies is a general consensus that three time constants are required to model the excited state relaxation, with a highly strained ground-state cis intermediate formed in less than 2.4 ps. Persistent strain drives the sequence of structural transitions that ultimately produce the signaling state. Crystal packing forces produce a restoring force that slows somewhat the rates of interconversion between the intermediates. Moreover, the solvent composition surrounding PYP can influence the number and structures of intermediates as well as the rates at which they interconvert. When chloride is present, the PYP photocycle in a crystal closely tracks that in solution, which suggests the epic movie of the PYP photocycle is indeed based in reality.

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.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Measurement of charges and chemical bonding in a cryo-EM structure

Hydrogen bonding, bond polarity, and charges in protein molecules play critical roles in the stabilization of protein structures, as well as affecting their functions such as enzymatic catalysis, electron transfer, and ligand binding. These effects can potentially be measured in Coulomb potentials using cryogenic electron microscopy (cryo-EM). We here present charges and bond properties of hydrogen in a sub-1.2 Å resolution structure of a protein complex, apoferritin, by single-particle cryo-EM. A weighted difference map reveals positive densities for most hydrogen atoms in the core region of the complex, while negative densities around acidic amino-acid side chains are likely related to negative charges. The former positive densities identify the amino- and oxo-termini of asparagine and glutamine side chains. The latter observations were verified by spatial-resolution selection and a dose-dependent frame series. The average position of the hydrogen densities depends on the parent bonded-atom type, and this is validated by the estimated level of the standard uncertainties in the bond lengths.

Chemical Features and Machine Learning Assisted Predictions of Protein-Ligand Short Hydrogen Bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 A closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Brønsted Acids Promote Olefin Oxidations by Bioinspired Nonheme CoIII(PhIO)(OH) Complexes: A Role for Low-Barrier Hydrogen Bonds

Introduction of Brønsted acids into biomimetic nonheme reactions promotes the oxidative ability of metal-oxygen complexes significantly. However, the molecular machinery of the promoted effects is missing. Herein, a comprehensive investigation of styrene oxidation by a cobalt(III)-iodosylbenzene complex, [(TQA)Co^III(OIPh)(OH)]²⁺ (1, TQA = tris(2-quinolylmethyl)amine), in the presence and absence of triflic acid (HOTf) was performed using density functional theory calculations. Results revealed for the first time that there is a low-barrier hydrogen bond (LBHB) between HOTf and the hydroxyl ligand of 1, which forms two valence-resonance structures [(TQA)Co^III(OIPh)(HO^---HOTf)]²⁺ (1_LBHB) and [(TQA)Co^III(OIPh)(H₂O--OTf^-)]²⁺ (1'_LBHB). Due to the oxo-wall, these complexes (1_LBHB and 1'_LBHB) cannot convert to high-valent cobalt-oxyl species. Instead, styrene oxidation by these oxidants (1_LBHB and 1'_LBHB) shows novel spin-state selectivity, i.e., on the ground closed-shell singlet state, styrene is oxidized to an epoxide, whereas on the excited triplet and quintet states, an aldehyde product, phenylacetaldehyde, is formed. The preferred pathway is styrene oxidation by 1'_LBHB, which is initiated by a rate-limiting bond-formation-coupled electron transfer process with an energy barrier of 12.2 kcal mol^-1. The nascent PhIO-styrene-radical-cation intermediate undergoes an intramolecular rearrangement to produce an aldehyde. The halogen bond between the OH^-/H₂O ligand and the iodine of PhIO modulates the activity of the cobalt-iodosylarene complexes 1_LBHB and 1'_LBHB. These new mechanistic findings enrich our knowledge of nonheme chemistry and hypervalent iodine chemistry and will play a positive role in the rational design of new catalysts.

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Absorption wavelength along chromophore low-barrier hydrogen bonds

In low-barrier hydrogen bonds (H-bonds), the pK_a values for the H-bond donor and acceptor moieties are nearly equal, whereas the redox potential values depend on the H⁺ position. Spectroscopic details of low-barrier H-bonds remain unclear. Here, we report the absorption wavelength along low-barrier H-bonds in protein environments, using a quantum mechanical/molecular mechanical approach. Low-barrier H-bonds form between Glu46 and p-coumaric acid (pCA) in the intermediate pR_CW state of photoactive yellow protein and between Asp116 and the retinal Schiff base in the intermediate M-state of the sodium-pumping rhodopsin KR2. The H⁺ displacement of only ∼0.4 Å, which does not easily occur without low-barrier H-bonds, is responsible for the ∼50 nm-shift in the absorption wavelength. This may be a basis of how photoreceptor proteins have evolved to proceed photocycles using abundant protons.

•

The low-barrier H-bond formation is a prerequisite for proton transfer

•

How the absorption wavelength changes as H⁺ moves is an open question

•

The H⁺ displacement of ∼0.4 Å leads to the absorption wavelength shift of ∼50 nm

•

The localization of the molecular orbitals plays a key role in the wavelength shift

Re-evaluation of protein neutron crystallography with and without X-ray/neutron joint refinement

Protein neutron crystallography is a powerful technique to determine the positions of H atoms, providing crucial biochemical information such as the protonation states of catalytic groups and the geometry of hydrogen bonds. Recently, the crystal structure of a bacterial copper amine oxidase was determined by joint refinement using X-ray and neutron diffraction data sets at resolutions of 1.14 and 1.72 Å, respectively [Murakawa et al. (2020 ▸). Proc. Natl Acad. Sci. USA, 117, 10818-10824]. While joint refinement is effective for the determination of the accurate positions of heavy atoms on the basis of the electron density, the structural information on light atoms (hydrogen and deuterium) derived from the neutron diffraction data might be affected by the X-ray data. To unravel the information included in the neutron diffraction data, the structure determination was conducted again using only the neutron diffraction data at 1.72 Å resolution and the results were compared with those obtained in the previous study. Most H and D atoms were identified at essentially the same positions in both the neutron-only and the X-ray/neutron joint refinements. Nevertheless, neutron-only refinement was found to be less effective than joint refinement in providing very accurate heavy-atom coordinates that lead to significant improvement of the neutron scattering length density map, especially for the active-site cofactor. Consequently, it was confirmed that X-ray/neutron joint refinement is crucial for determination of the real chemical structure of the catalytic site of the enzyme.

Effective prediction of short hydrogen bonds in proteins via machine learning method

Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.

Newly identified C–H⋯O hydrogen bond in histidine

Histidine C–H bonds observed in protein structures include (clockwise from top left): myoglobin, β-lactamase, and photoactive yellow protein; calculations indicate that tautomeric/protonation state influences H-bonding ability (bottom left).

Protein Response Effects on Cofactor Excitation Energies from First Principles: Augmenting Subsystem Time-Dependent Density-Functional Theory with Many-Body Expansion Techniques

We investigate the possibility of describing protein response effects on a chromophore excitation by means of subsystem time-dependent density-functional theory (sTDDFT) in combination with a many-body expansion (MBE) approach. While sTDDFT is in principle intrinsically able to include such contributions, addressing cofactor excitations in protein models or entire proteins with full environment-response treatments is currently out of reach. Taking different model structures of the green fluorescent protein (GFP) and bovine rhodopsin as examples, we demonstrate that an embedded-MBE approach based on sTDDFT in its simplest version leads to a good agreement of the predicted protein response effect already at second order. To reproduce reference response effects from nonsubsystem TDDFT calculations quantitatively (error ≤ 5%), however, a third- or even fourth-order MBE may be required. For the latter case, we explore a selective inclusion of fourth-order terms that drastically reduces the computational burden. In addition, we demonstrate how this sTDDFT-MBE treatment can be utilized as an analysis tool to identify residues with dominant response contributions. This, in turn, can be employed to arrive at smaller structural models for light-absorbing proteins, which still feature all of the main characteristics in terms of photoresponse properties.

Aspartate or arginine? Validated redox state X-ray structures elucidate mechanistic subtleties of FeIV = O formation in bacterial dye-decolorizing peroxidases

Structure determination of proteins and enzymes by X-ray crystallography remains the most widely used approach to complement functional and mechanistic studies. Capturing the structures of intact redox states in metalloenzymes is critical for assigning the chemistry carried out by the metal in the catalytic cycle. Unfortunately, X-rays interact with protein crystals to generate solvated photoelectrons that can reduce redox active metals and hence change the coordination geometry and the coupled protein structure. Approaches to mitigate such site-specific radiation damage continue to be developed, but nevertheless application of such approaches to metalloenzymes in combination with mechanistic studies are often overlooked. In this review, we summarize our recent structural and kinetic studies on a set of three heme peroxidases found in the bacterium Streptomyces lividans that each belong to the dye decolourizing peroxidase (DyP) superfamily. Kinetically, each of these DyPs has a distinct reactivity with hydrogen peroxide. Through a combination of low dose synchrotron X-ray crystallography and zero dose serial femtosecond X-ray crystallography using an X-ray free electron laser (XFEL), high-resolution structures with unambiguous redox state assignment of the ferric and ferryl (Fe^IV = O) heme species have been obtained. Experiments using stopped-flow kinetics, solvent-isotope exchange and site-directed mutagenesis with this set of redox state validated DyP structures have provided the first comprehensive kinetic and structural framework for how DyPs can modulate their distal heme pocket Asp/Arg dyad to use either the Asp or the Arg to facilitate proton transfer and rate enhancement of peroxide heterolysis.

On-off transition and ultrafast decay of amino acid luminescence driven by modulation of supramolecular packing

Luminescence of biomolecules in the visible range of the spectrum has been experimentally observed upon aggregation, contrary to their monomeric state. However, the physical basis for this phenomenon is still elusive. Here, we systematically examine all coded amino acids to provide non-biased empirical insights. Several amino acids, including non-aromatic, show intense visible luminescence. Lysine crystals display the highest signal, whereas the very chemically similar non-coded ornithine does not, implying a role for molecular packing rather than the chemical characteristics. Furthermore, cysteine shows luminescence that is indeed crystal packing dependent as repeated rearrangements between two crystal structures result in a reversible on-off optical transition. In addition, ultrafast lifetime decay is experimentally validated, corroborating a recently raised hypothesis regarding the governing role of nπ∗ states in the emission formation. Collectively, our study supports that electronic interactions between non-fluorescent, non-absorbing molecules at the monomeric state may result in reversible optically active states by the formation of supramolecular fluorophores.

Neutron crystallography for the elucidation of enzyme catalysis

Hydrogen atoms and hydration water molecules in proteins are indispensable for many biochemical processes, especially enzymatic catalysis. The locations of hydrogen atoms in proteins are usually predicted based on X-ray structures, but it is still very difficult to know the ionization states of the catalytic residues, the hydration structure of the protein, and the characteristics of hydrogen-bonding interactions. Neutron crystallography allows the direct observation of hydrogen atoms that play crucial roles in molecular recognition and the catalytic reactions of enzymes. In this review, we present the current status of neutron crystallography in structural biology and recent neutron structural analyses of three enzymes: ascorbate peroxidase, the main protease of severe acute respiratory syndrome coronavirus 2, and copper-containing nitrite reductase.

Combined NMR and UV-Vis Spectroscopic Studies of Models for the Hydrogen Bond System in the Active Site of Photoactive Yellow Protein: H-Bond Cooperativity and Medium Effects

Intramolecular hydrogen bonds in aprotic media were studied by combined (simultaneous) NMR and UV-vis spectroscopy. The species under investigation were anionic and featured single or coupled H-bonds between, for example, carboxylic groups and phenolic oxygen atoms (COO···H···OC)^-, among phenolic oxygen atoms (CO···H···OC)^-, and hydrogen bond chains between a carboxylic group and two phenolic oxygen atoms (COO···H···(OC)···H···OC)^-. The last anion may be regarded as a small molecule model for the hydrogen bond system in the active site of wild-type photoactive yellow protein (PYP) while the others mimic the corresponding H-bonds in site-selective mutants. Proton positions in isolated hydrogen bonds and hydrogen bond chains were assessed by calculations for vacuum conditions and spectroscopically for the two media, CD₂Cl₂ and the liquefied gas mixture CDClF₂/CDF₃ at low temperatures. NMR parameters allow for the estimation of time-averaged H-bond geometries, and optical spectra give additional information about geometry distributions. Comparison of the results from the various systems revealed the effects of the formation of hydrogen bond chains and changes of medium conditions on the geometry of individual H-bonds. In particular, the proton in a hydrogen bond to a carboxylic group shifts from the phenolic oxygen atom in the system COO^-···H-OC to the carboxylic group in COO-H···(OC)^-···H-OC as a result of hydrogen bond formation to the additional phenolic donor. Increase in medium polarity may, however, induce the conversion of a structure of a type COO-H···(OC)^-···H-OC to the type COO^-···H-(OC)···H-OC. Application of these results obtained from the model systems to PYP suggests that both cooperative effects within the hydrogen bond chain and a low-polarity protein environment are prerequisites for the stabilization of negative charge on the cofactor and hence for the spectral tuning of the photoreceptor.

Short hydrogen bonds enhance nonaromatic protein-related fluorescence

Intrinsic fluorescence of nonaromatic amino acids is a puzzling phenomenon with an enormous potential in biophotonic applications. The physical origins of this effect, however, remain elusive. Herein, we demonstrate how specific hydrogen bond networks can modulate fluorescence. We highlight the key role played by short hydrogen bonds, present in the protein structure, on the ensuing fluorescence. We provide detailed experimental and molecular evidence to explain these unusual nonaromatic optical properties. Our findings should benefit the design of novel optically active biomaterials for applications in biosensing and imaging.

Fluorescence in biological systems is usually associated with the presence of aromatic groups. Here, by employing a combined experimental and computational approach, we show that specific hydrogen bond networks can significantly affect fluorescence. In particular, we reveal that the single amino acid L-glutamine, by undergoing a chemical transformation leading to the formation of a short hydrogen bond, displays optical properties that are significantly enhanced compared with L-glutamine itself. Ab initio molecular dynamics simulations highlight that these short hydrogen bonds prevent the appearance of a conical intersection between the excited and the ground states and thereby significantly decrease nonradiative transition probabilities. Our findings open the door to the design of new photoactive materials with biophotonic applications.

Controlling Light-Induced Proton Transfer from the GFP Chromophore

Green Fluorescent Protein (GFP) is known to undergo excited-state proton transfer (ESPT). Formation of a short H-bond favors ultrafast ESPT in GFP-like proteins, such as the GFP S65T/H148D mutant, but the detailed mechanism and its quantum nature remain to be resolved. Here we study in vacuo, light-induced proton transfer from the GFP chromophore in hydrogen-bonded complexes with two anionic proton acceptors, I^- and deprotonated trichloroacetic acid (TCA^- ). We address the role of the strong H-bond and the quantum mechanical proton-density distribution in the excited state, which determines the proton-transfer probability. Our study shows that chemical modifications to the molecular network drastically change the proton-transfer probability and it can become strongly wavelength dependent. The proton-transfer branching ratio is found to be 60 % for the TCA complex and 10 % for the iodide complex, being highly dependent on the photon energy in the latter case. Using high-level ab initio calculations, we show that light-induced proton transfer takes place in S₁ , revealing intrinsic photoacid properties of the isolated GFP chromophore in strongly bound H-bonded complexes. ESPT is found to be very sensitive to the topography of the highly anharmonic potential in S₁ , depending on the quantum-density distribution upon vibrational excitation. We also show that the S₁ potential-energy surface, and hence excited-state proton transfer, can be controlled by altering the chromophore microenvironment.

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.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Serial Femtosecond Zero Dose Crystallography Captures a Water-Free Distal Heme Site in a Dye-Decolorising Peroxidase to Reveal a Catalytic Role for an Arginine in FeIV =O Formation

Obtaining structures of intact redox states of metal centers derived from zero dose X-ray crystallography can advance our mechanistic understanding of metalloenzymes. In dye-decolorising heme peroxidases (DyPs), controversy exists regarding the mechanistic role of the distal heme residues aspartate and arginine in the heterolysis of peroxide to form the catalytic intermediate compound I (FeIV =O and a porphyrin cation radical). Using serial femtosecond X-ray crystallography (SFX), we have determined the pristine structures of the FeIII and FeIV =O redox states of a B-type DyP. These structures reveal a water-free distal heme site that, together with the presence of an asparagine, imply the use of the distal arginine as a catalytic base. A combination of mutagenesis and kinetic studies corroborate such a role. Our SFX approach thus provides unique insight into how the distal heme site of DyPs can be tuned to select aspartate or arginine for the rate enhancement of peroxide heterolysis.

Electron transfer and spin dynamics of the radical-pair in the cryptochrome from Chlamydomonas reinhardtii by computational analysis

In an effort to elucidate the origin of avian magnetoreception, it was postulated that a radical-pair formed in a cryptochrome upon light activation provided the basis for the mechanism that enables an inclination compass sensitive to the geomagnetic field. Photoreduction in this case involves formation of a flavin adenine dinucleotide (FAD)-tryptophan (TRP) radical-pair, following electron transfer within a conserved TRP triad in the cryptochrome. Recently, an animal-like cryptochrome from Chlamydomonas reinhardtii (CraCRY) was analyzed, demonstrating the role of a fourth aromatic residue, which serves as a terminal electron donor in the photoreduction pathway, resulting in the creation of a more distal radical-pair and exhibiting fast electron transfer. In this work, we investigated the electron transfer in CraCRY with a combination of free energy molecular dynamics (MD) simulations, frozen density functional theory, and QM/MM MD simulations, supporting the suggestion of a proton coupled electron transfer mechanism. Spin dynamics simulations discerned details on the dependence of the singlet yield on the direction of the external magnetic field for the [FAD^•- TYRH^•+] and [FAD^•- TYR^•] radical-pairs in CraCRY, in comparison with the previously modeled [FAD^•- TRPH^•+] radical-pair.

Opportunities and challenges in neutron crystallography

Neutron and X-ray crystallography are complementary to each other. While X-ray scattering is directly proportional to the number of electrons of an atom, neutrons interact with the atomic nuclei themselves. Neutron crystallography therefore provides an excellent alternative in determining the positions of hydrogens in a biological molecule. In particular, since highly polarized hydrogen atoms (H+) do not have electrons, they cannot be observed by X-rays. Neutron crystallography has its own limitations, mainly due to inherent low flux of neutrons sources, and as a consequence, the need for much larger crystals and for different data collection and analysis strategies. These technical challenges can however be overcome to yield crucial structural insights about protonation states in enzyme catalysis, ligand recognition, as well as the presence of unusual hydrogen bonds in proteins.

Neutron crystallography of copper amine oxidase reveals keto/enolate interconversion of the quinone cofactor and unusual proton sharing

Recent advances in neutron crystallographic studies have provided structural bases for quantum behaviors of protons observed in enzymatic reactions. Thus, we resolved the neutron crystal structure of a bacterial copper (Cu) amine oxidase (CAO), which contains a prosthetic Cu ion and a protein-derived redox cofactor, topa quinone (TPQ). We solved hitherto unknown structures of the active site, including a keto/enolate equilibrium of the cofactor with a nonplanar quinone ring, unusual proton sharing between the cofactor and the catalytic base, and metal-induced deprotonation of a histidine residue that coordinates to the Cu. Our findings show a refined active-site structure that gives detailed information on the protonation state of dissociable groups, such as the quinone cofactor, which are critical for catalytic reactions.

Visualizing the protons in a metalloenzyme electron proton transfer pathway

In redox metalloenzymes, the process of electron transfer often involves the concerted movement of a proton. These processes are referred to as proton-coupled electron transfer, and they underpin a wide variety of biological processes, including respiration, energy conversion, photosynthesis, and metalloenzyme catalysis. The mechanisms of proton delivery are incompletely understood, in part due to an absence of information on exact proton locations and hydrogen bonding structures in a bona fide metalloenzyme proton pathway. Here, we present a 2.1-Å neutron crystal structure of the complex formed between a redox metalloenzyme (ascorbate peroxidase) and its reducing substrate (ascorbate). In the neutron structure of the complex, the protonation states of the electron/proton donor (ascorbate) and all of the residues involved in the electron/proton transfer pathway are directly observed. This information sheds light on possible proton movements during heme-catalyzed oxygen activation, as well as on ascorbate oxidation.

Direct Observation of the Protonation States in the Mutant Green Fluorescent Protein

Neutron crystallography has been used to elucidate the protonation states for the enhanced green fluorescent protein, which has revolutionized imaging technologies. The structure has a deprotonated hydroxyl group in the fluorescent chromophore. Also, the protonation states of His148 and Thr203, as well as the orientation of a critical water molecule in direct contact with the chromophore, could be determined. The results demonstrate that the deprotonated hydroxyl group in the chromophore and the nitrogen atom ND1 in His148 are charged negatively and positively, respectively, forming an ion pair. The position of the two deuterium atoms in the critical water molecule appears to be displaced slightly toward the acceptor oxygen atoms according to their omit maps. This displacement implies the formation of an intriguing electrostatic potential realized inside of the protein. Our findings provide new insights into future protein design strategies along with developments in quantum chemical calculations.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in Phenix

Diffraction (X-ray, neutron and electron) and electron cryo-microscopy are powerful methods to determine three-dimensional macromolecular structures, which are required to understand biological processes and to develop new therapeutics against diseases. The overall structure-solution workflow is similar for these techniques, but nuances exist because the properties of the reduced experimental data are different. Software tools for structure determination should therefore be tailored for each method. Phenix is a comprehensive software package for macromolecular structure determination that handles data from any of these techniques. Tasks performed with Phenix include data-quality assessment, map improvement, model building, the validation/rebuilding/refinement cycle and deposition. Each tool caters to the type of experimental data. The design of Phenix emphasizes the automation of procedures, where possible, to minimize repetitive and time-consuming manual tasks, while default parameters are chosen to encourage best practice. A graphical user interface provides access to many command-line features of Phenix and streamlines the transition between programs, project tracking and re-running of previous tasks.

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.

Fifth-order time-domain Raman spectroscopy of photoactive yellow protein for visualizing vibrational coupling in its excited state

We report fifth-order time-domain Raman spectroscopy of photoactive yellow protein (PYP), with the aim to visualize vibrational coupling in its excited state. After the ultrashort actinic pump pulse prepared the vibrational coherence and population in the excited state, the evolving vibrational structure was tracked by time-resolved impulsive stimulated Raman spectroscopy using sub-7-fs pulses. The obtained fifth-order time-domain Raman data were translated to a two-dimensional (2D) frequency-frequency correlation map, which visualizes the correlation between low- and high-frequency vibrational modes of the excited state. The 2D map of PYP reveals a cross peak, indicating the coupling between the phenolic C─O stretch mode of the chromophore and the low-frequency modes (~160 cm^-1), assignable to the intermolecular motions involving the surrounding hydrogen-bonded amino acids. The unveiled coupling suggests the importance of the low-frequency vibrational motion in the primary photoreaction of PYP, highlighting the unique capability of this spectroscopic approach for studying ultrafast reaction dynamics.

Evidence of proton-coupled mixed-valency by electrochemical behavior on transition metal complex dimers bridged by two Ag+ ions

H-Bonded metal complex dimers with reversible redox behaviour, which are connected by a low-barrier hydrogen bond (LBHB) with a very low energy barrier for proton transfer, can provide a unique mixed-valency state stabilized by the proton-coupled electron transfer (PCET) phenomenon. Using cyclic voltammetry measurements, newly prepared [ReIIICl2(PnPr3)2(Hbim)]2 (2) and [OsIIICl2(PnPr3)2(Hbim)]2 (3) existing as H-bonded dimers in a CH2Cl2 solution showed a four-step and four-electron transfer containing two mixed-valency states of ReIIReIII and ReIIIReIV, and OsIIOsIII and OsIIIOsVI, respectively. Furthermore, [ReIIICl2(PnPr3)2(Agbim)]2 (4) and [OsIIICl2(PnPr3)2(Agbim)]2 (5), bridged by two Ag+ ions instead of two H-bonding protons, were prepared, and their electrochemical behaviours changed to a two-step and four-electron transfer. It is clear that the H-bonded complex dimers 2 and 3, connected by an LBHB, can be electrochemically stabilized into unique pairs of mixed-valency states by PCET, and the H-bonding proton transfer also controls the electrochemical redox behaviour.

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.