Vibrational analysis on the revised potential energy curve of the low-barrier hydrogen bond in photoactive yellow protein

Nuclear quantum effects on the intramolecular hydrogen bonds in biuret and biguanide

We focus on the unique aspects of biuret and biguanide, which form six-membered ring structures via intramolecular hydrogen bonds. The proton donor and acceptor atoms differ between biuret and biguanide, leading to varying energy barrier heights for proton transfer. We performed path integral molecular dynamics (PIMD) simulations for biuret and biguanide to investigate the correlation between proton transfer and the degree of the delocalization of π-electrons in the six-membered ring framework structure. The results indicate that the π-electrons in the framework structure are delocalized regardless of the ease of intramolecular proton transfer.

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.

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.

Neutron crystallography of photoactive yellow protein reveals unusual protonation state of Arg52 in the crystal

Because of its high pK_a, arginine (Arg) is believed to be protonated even in the hydrophobic environment of the protein interior. However, our neutron crystallographic structure of photoactive yellow protein, a light sensor, demonstrated that Arg52 adopts an electrically neutral form. We also showed that the hydrogen bond between the chromophore and Glu46 is a so-called low barrier hydrogen bond (LBHB). Because both the neutral Arg and LBHB are unusual in proteins, these observations remain controversial. To validate our findings, we carried out neutron crystallographic analysis of the E46Q mutant of PYP. The resultant structure revealed that the proportion of the cationic form is higher in E46Q than in WT, although the cationic and neutral forms of Arg52 coexist in E46Q. These observations were confirmed by the occupancy of the deuterium atom bound to the N _η1 atom combined with an alternative conformation of the N_(η2)D₂ group comprising sp² hybridisation. Based on these results, we propose that the formation of the LBHB decreases the proton affinity of Arg52, stabilizing the neutral form in the crystal.

The Low Barrier Hydrogen Bond in the Photoactive Yellow Protein: A Vacuum Artifact Absent in the Crystal and Solution

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large ¹H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK_a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical ¹H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.