Low-barrier hydrogen bond in photoactive yellow protein

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.

Evaluation of models determined by neutron diffraction and proposed improvements to their validation and deposition

The Protein Data Bank (PDB) contains a growing number of models that have been determined using neutron diffraction or a hybrid method that combines X-ray and neutron diffraction. The advantage of neutron diffraction experiments is that the positions of all atoms can be determined, including H atoms, which are hardly detectable by X-ray diffraction. This allows the determination of protonation states and the assignment of H atoms to water molecules. Because neutrons are scattered differently by hydrogen and its isotope deuterium, neutron diffraction in combination with H/D exchange can provide information on accessibility, dynamics and chemical lability. In this study, the deposited data, models and model-to-data fit for all PDB entries that used neutron diffraction as the source of experimental data have been analysed. In many cases, the reported Rwork and Rfree values were not reproducible. In such cases, the model and data files were analysed to identify the reasons for this mismatch. The issues responsible for the discrepancies are summarized and explained. The analysis unveiled limitations to the annotation, deposition and validation of models and data, and a lack of community-wide accepted standards for the description of neutron models and data, as well as deficiencies in current model refinement tools. Most of the issues identified concern the handling of H atoms. Since the primary use of neutron macromolecular crystallography is to locate and directly visualize H atoms, it is important to address these issues, so that the deposited neutron models allow the retrieval of the maximum amount of information with the smallest effort of manual intervention. A path forward to improving the annotation, validation and deposition of neutron models and hybrid X-ray and neutron models is suggested.

Excitation-Wavelength-Dependent Photocycle Initiation Dynamics Resolve Heterogeneity in the Photoactive Yellow Protein from Halorhodospira halophila

Photoactive yellow proteins (PYPs) make up a diverse class of blue-light-absorbing bacterial photoreceptors. Electronic excitation of the p-coumaric acid chromophore covalently bound within PYP results in triphasic quenching kinetics; however, the molecular basis of this behavior remains unresolved. Here we explore this question by examining the excitation-wavelength dependence of the photodynamics of the PYP from Halorhodospira halophila via a combined experimental and computational approach. The fluorescence quantum yield, steady-state fluorescence emission maximum, and cryotrapping spectra are demonstrated to depend on excitation wavelength. We also compare the femtosecond photodynamics in PYP at two excitation wavelengths (435 and 475 nm) with a dual-excitation-wavelength-interleaved pump-probe technique. Multicompartment global analysis of these data demonstrates that the excited-state photochemistry of PYP depends subtly, but convincingly, on excitation wavelength with similar kinetics with distinctly different spectral features, including a shifted ground-state beach and altered stimulated emission oscillator strengths and peak positions. Three models involving multiple excited states, vibrationally enhanced barrier crossing, and inhomogeneity are proposed to interpret the observed excitation-wavelength dependence of the data. Conformational heterogeneity was identified as the most probable model, which was supported with molecular mechanics simulations that identified two levels of inhomogeneity involving the orientation of the R52 residue and different hydrogen bonding networks with the p-coumaric acid chromophore. Quantum calculations were used to confirm that these inhomogeneities track to altered spectral properties consistent with the experimental results.

Neutron macromolecular crystallography

Neutron diffraction techniques permit direct determination of the hydrogen (H) and deuterium (D) positions in crystal structures of biological macromolecules at resolutions of ∼1.5 and 2.5 Å, respectively. In addition, neutron diffraction data can be collected from a single crystal at room temperature without radiation damage issues. By locating the positions of H/D-atoms, protonation states and water molecule orientations can be determined, leading to a more complete understanding of many biological processes and drug-binding. In the last ca. 5 years, new beamlines have come online at reactor neutron sources, such as BIODIFF at Heinz Maier-Leibnitz Zentrum and IMAGINE at Oak Ridge National Laboratory (ORNL), and at spallation neutron sources, such as MaNDi at ORNL and iBIX at the Japan Proton Accelerator Research Complex. In addition, significant improvements have been made to existing beamlines, such as LADI-III at the Institut Laue-Langevin. The new and improved instrumentations are allowing sub-mm3 crystals to be regularly used for data collection and permitting the study of larger systems (unit-cell edges >100 Å). Owing to this increase in capacity and capability, many more studies have been performed and for a wider range of macromolecules, including enzymes, signalling proteins, transport proteins, sugar-binding proteins, fluorescent proteins, hormones and oligonucleotides; of the 126 structures deposited in the Protein Data Bank, more than half have been released since 2013 (65/126, 52%). Although the overall number is still relatively small, there are a growing number of examples for which neutron macromolecular crystallography has provided the answers to questions that otherwise remained elusive.

Unique Hydrogen Bonds in Membrane Protein Monitored by Whole Mid-IR ATR Spectroscopy in Aqueous Solution

Protein function is coupled to its structural changes, for which stimulus-induced difference Fourier-transform infrared (FTIR) spectroscopy is a powerful method. By optimizing the attenuated total reflection (ATR)-FTIR analysis on sodium-pumping rhodopsin KR2 in aqueous solution, we first measured the accurate difference spectra upon sodium binding in the whole IR region (4000-1000 cm^-1). The new spectral window allows the analysis of not only the fingerprint region (1800-1000 cm^-1) but also the hydrogen-bonding donor region (4000-1800 cm^-1), revealing an unusually strong hydrogen bond of Tyr located in the sodium binding site of KR2. Progress in ATR-FTIR difference spectroscopy provides an approach to investigating stimulus-induced structural changes of membrane proteins under physiological aqueous conditions.

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.

Computerized implementation of higher-order electron-correlation methods and their linear-scaling divide-and-conquer extensions

We have implemented a linear-scaling divide-and-conquer (DC)-based higher-order coupled-cluster (CC) and Møller-Plesset perturbation theories (MPPT) as well as their combinations automatically by means of the tensor contraction engine, which is a computerized symbolic algebra system. The DC-based energy expressions of the standard CC and MPPT methods and the CC methods augmented with a perturbation correction were proposed for up to high excitation orders [e.g., CCSDTQ, MP4, and CCSD(2)_TQ ]. The numerical assessment for hydrogen halide chains, polyene chains, and first coordination sphere (C1) model of photoactive yellow protein has revealed that the DC-based correlation methods provide reliable correlation energies with significantly less computational cost than that of the conventional implementations. © 2017 Wiley Periodicals, Inc.

Active-Site pKa Determination for Photoactive Yellow Protein Rationalizes Slow Ground-State Recovery

The ability to avoid blue-light radiation is crucial for bacteria to survive. In Halorhodospira halophila, the putative receptor for this response is known as photoactive yellow protein (PYP). Its response to blue light is mediated by changes in the optical properties of the chromophore para-coumaric acid (pCA) in the protein active site. PYP displays photocycle kinetics with a strong pH dependence for ground-state recovery, which has remained enigmatic. To resolve this problem, a comprehensive pK_a determination of the active-site residues of PYP is required. Herein, we show that Glu-46 stays protonated from pH 3.4 to pH 11.4 in the ground (pG) state. This conclusion is supported by the observed hydrogen-bonded protons between Glu-46 and pCA and Tyr-42 and pCA, which are persistent over the entire pH range. Our experimental results show that none of the active-site residues of PYP undergo pH-induced changes in the pG state. Ineluctably, the pH dependence of pG recovery is linked to conformational change that is dependent upon the population of the relevant protonation state of Glu-46 and the pCA chromophore in the excited state, collaterally explaining why pG recovery is slow.

Neutron structure of the T26H mutant of T4 phage lysozyme provides insight into the catalytic activity of the mutant enzyme and how it differs from that of wild type

T4 phage lysozyme is an inverting glycoside hydrolase that degrades the murein of bacterial cell walls by cleaving the β‐1,4‐glycosidic bond. The substitution of the catalytic Thr26 residue to a histidine converts the wild type from an inverting to a retaining enzyme, which implies that the original general acid Glu11 can also act as an acid/base catalyst in the hydrolysis. Here, we have determined the neutron structure of the perdeuterated T26H mutant to clarify the protonation states of Glu11 and the substituted His26, which are key in the retaining reaction. The 2.09‐Å resolution structure shows that the imidazole group of His26 is in its singly protonated form in the active site, suggesting that the deprotonated Nɛ2 atom of His26 can attack the anomeric carbon of bound substrate as a nucleophile. The carboxyl group of Glu11 is partially protonated and interacts with the unusual neutral state of the guanidine moiety of Arg145, as well as two heavy water molecules. Considering that one of the water‐binding sites has the potential to be occupied by a hydronium ion, the bulk solvent could be the source for the protonation of Glu11. The respective protonation states of Glu11 and His26 are consistent with the bond lengths determined by an unrestrained refinement of the high‐resolution X‐ray structure of T26H at 1.04‐Å resolution. The detail structural information, including the coordinates of the deuterium atoms in the active site, provides insight into the distinctively different catalytic activities of the mutant and wild type enzymes.

PDB Code(s): 5XPE; 5XPF

Neutron Crystallography for the Study of Hydrogen Bonds in Macromolecules

The hydrogen bond (H bond) is one of the most important interactions that form the foundation of secondary and tertiary protein structure. Beyond holding protein structures together, H bonds are also intimately involved in solvent coordination, ligand binding, and enzyme catalysis. The H bond by definition involves the light atom, H, and it is very difficult to study directly, especially with X-ray crystallographic techniques, due to the poor scattering power of H atoms. Neutron protein crystallography provides a powerful, complementary tool that can give unambiguous information to structural biologists on solvent organization and coordination, the electrostatics of ligand binding, the protonation states of amino acid side chains and catalytic water species. The method is complementary to X-ray crystallography and the dynamic data obtainable with NMR spectroscopy. Also, as it gives explicit H atom positions, it can be very valuable to computational chemistry where exact knowledge of protonation and solvent orientation can make a large difference in modeling. This article gives general information about neutron crystallography and shows specific examples of how the method has contributed to structural biology, structure-based drug design; and the understanding of fundamental questions of reaction mechanisms.

Probing the early stages of photoreception in photoactive yellow protein with ultrafast time-domain Raman spectroscopy

Unveiling the nuclear motions of photoreceptor proteins in action is a crucial goal in protein science in order to understand their elaborate mechanisms and how they achieve optimal selectivity and efficiency. Previous studies have provided detailed information on the structures of intermediates that appear during the later stages (>ns) of such photoreception cycles, yet the initial events immediately after photoabsorption remain unclear because of experimental challenges in monitoring nuclear rearrangements on ultrafast timescales, including protein-specific low-frequency motions. Using time-domain Raman probing with sub-7-fs pulses, we obtain snapshot vibrational spectra of photoactive yellow protein and a mutant with high sensitivity, providing insights into the key responses that drive photoreception. Our data show a drastic intensity drop of the excited-state marker band at 135 cm^-1 within a few hundred femtoseconds, suggesting a rapid weakening of the hydrogen bond that anchors the chromophore. We also track formation of the first ground-state intermediate over the first few picoseconds and fully characterize its vibrational structure, revealing a substantially-twisted cis conformation.

Photocycle of Photoactive Yellow Protein in Cell-Mimetic Environments: Molecular Volume Changes and Kinetics

Using various spectroscopic techniques such as UV-visible spectroscopy, circular dichroism spectroscopy, NMR spectroscopy, small-angle X-ray scattering, transient grating, and transient absorption techniques, we investigated how cell-mimetic environments made by crowding influence the photocycle of photoactive yellow protein (PYP) in terms of the molecular volume change and kinetics. Upon addition of molecular crowding agents, the ratio of the diffusion coefficient of the blue-shifted intermediate (pB) to that of the ground species (pG) significantly changes from 0.92 and approaches 1.0. This result indicates that the molecular volume change accompanied by the photocycle of PYP in molecularly crowded environments is much smaller than that which occurs in vitro and that the pB intermediate under crowded environments favors a compact conformation due to the excluded volume effect. The kinetics of the photocycle of PYP in cell-mimetic environments is greatly decelerated by the dehydration, owing to the interaction between the protein and small crowding agents, but is barely affected by the excluded volume effect. The results lead to the inference that the signaling transducer of PYP may not necessarily utilize the conformational change of PYP to sense the signaling state.

An extended N-H bond, driven by a conserved second-order interaction, orients the flavin N5 orbital in cholesterol oxidase

The protein microenvironment surrounding the flavin cofactor in flavoenzymes is key to the efficiency and diversity of reactions catalysed by this class of enzymes. X-ray diffraction structures of oxidoreductase flavoenzymes have revealed recurrent features which facilitate catalysis, such as a hydrogen bond between a main chain nitrogen atom and the flavin redox center (N5). A neutron diffraction study of cholesterol oxidase has revealed an unusual elongated main chain nitrogen to hydrogen bond distance positioning the hydrogen atom towards the flavin N5 reactive center. Investigation of the structural features which could cause such an unusual occurrence revealed a positively charged lysine side chain, conserved in other flavin mediated oxidoreductases, in a second shell away from the FAD cofactor acting to polarize the peptide bond through interaction with the carbonyl oxygen atom. Double-hybrid density functional theory calculations confirm that this electrostatic arrangement affects the N-H bond length in the region of the flavin reactive center. We propose a novel second-order partial-charge interaction network which enables the correct orientation of the hydride receiving orbital of N5. The implications of these observations for flavin mediated redox chemistry are discussed.

Short Carboxylic Acid-Carboxylate Hydrogen Bonds Can Have Fully Localized Protons

Short hydrogen bonds (H-bonds) have been proposed to play key functional roles in several proteins. The location of the proton in short H-bonds is of central importance, as proton delocalization is a defining feature of low-barrier hydrogen bonds (LBHBs). Experimentally determining proton location in H-bonds is challenging. Here, bond length analysis of atomic (1.15-0.98 Å) resolution X-ray crystal structures of the human protein DJ-1 and its bacterial homologue, YajL, was used to determine the protonation states of H-bonded carboxylic acids. DJ-1 contains a buried, dimer-spanning 2.49 Å H-bond between Glu15 and Asp24 that satisfies standard donor-acceptor distance criteria for a LBHB. Bond length analysis indicates that the proton is localized on Asp24, excluding a LBHB at this location. However, similar analysis of the Escherichia coli homologue YajL shows both residues may be protonated at the H-bonded oxygen atoms, potentially consistent with a LBHB. A Protein Data Bank-wide screen identifies candidate carboxylic acid H-bonds in approximately 14% of proteins, which are typically short [⟨d_O-O⟩ = 2.542(2) Å]. Chemically similar H-bonds between hydroxylated residues (Ser/Thr/Tyr) and carboxylates show a trend of lengthening O-O distance with increasing H-bond donor pK_a. This trend suggests that conventional electronic effects provide an adequate explanation for short, charge-assisted carboxylic acid-carboxylate H-bonds in proteins, without the need to invoke LBHBs in general. This study demonstrates that bond length analysis of atomic resolution X-ray crystal structures provides a useful experimental test of certain candidate LBHBs.

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.

Mechanisms for Two-Step Proton Transfer Reactions in the Outward-Facing Form of MATE Transporter

Bacterial pathogens or cancer cells can acquire multidrug resistance, which causes serious clinical problems. In cells with multidrug resistance, various drugs or antibiotics are extruded across the cell membrane by multidrug transporters. The multidrug and toxic compound extrusion (MATE) transporter is one of the five families of multidrug transporters. MATE from Pyrococcus furiosus uses H(+) to transport a substrate from the cytoplasm to the outside of a cell. Crystal structures of MATE from P. furiosus provide essential information on the relevant H(+)-binding sites (D41 and D184). Hybrid quantum mechanical/molecular mechanical simulations and continuum electrostatic calculations on the crystal structures predict that D41 is protonated in one structure (Straight) and, both D41 and D184 protonated in another (Bent). All-atom molecular dynamics simulations suggest a dynamic equilibrium between the protonation states of the two aspartic acids and that the protonation state affects hydration in the substrate binding cavity and lipid intrusion in the cleft between the N- and C-lobes. This hypothesis is examined in more detail by quantum mechanical/molecular mechanical calculations on snapshots taken from the molecular dynamics trajectories. We find the possibility of two proton transfer (PT) reactions in Straight: the 1st PT takes place between side-chains D41 and D184 through a transient formation of low-barrier hydrogen bonds and the 2nd through another H(+) from the headgroup of a lipid that intrudes into the cleft resulting in a doubly protonated (both D41 and D184) state. The 1st PT affects the local hydrogen bond network and hydration in the N-lobe cavity, which would impinge on the substrate-binding affinity. The 2nd PT would drive the conformational change from Straight to Bent. This model may be applicable to several prokaryotic H(+)-coupled MATE multidrug transporters with the relevant aspartic acids.

Unambiguous Determination of Protein Arginine Ionization States in Solution by NMR Spectroscopy

Because arginine residues in proteins are expected to be in their protonated form almost without exception, reports demonstrating that a protein arginine residue is charge-neutral are rare and potentially controversial. Herein, we present a ¹³ C-detected NMR experiment for probing individual arginine residues in proteins notwithstanding the presence of chemical and conformational exchange effects. In the experiment, the ¹⁵ N^η and ¹⁵ N^ϵ chemical shifts of an arginine head group are correlated with that of the directly attached ¹³ C^ζ . In the resulting spectrum, the number of protons in the arginine head group can be obtained directly from the ¹⁵ N-¹ H scalar coupling splitting pattern. We applied this method to unambiguously determine the ionization state of the R52 side chain in the photoactive yellow protein from Halorhodospira halophila. Although only three H^η atoms were previously identified by neutron crystallography, we show that R52 is predominantly protonated in solution.

Transient low-barrier hydrogen bond in the photoactive state of green fluorescent protein

In this paper, we have analyzed the feasibility of spontaneous proton transfer in GFP at the Franck-Condon region directly after photoexcitation. Computation of a sizeable portion of the potential energy surface at the Franck-Condon region of A the structure shows the process of proton transfer to be unfavorable by 3 kcal mol(-1) in S1 if no further structural relaxation is permitted. The ground vibrational state is found to lie above the potential energy barrier of the proton transfer in both S0 and S1. Expectation values of the geometry reveal that the proton shared between the chromophore and W22, and the proton shared between Ser205 and Glu222 are very close to the center of the respective hydrogen bonds, giving support to the claim that the first transient intermediate detected after photoexcitation (I0*) has characteristics similar to those of a low-barrier hydrogen bond [Di Donato et al., Phys. Chem. Chem. Phys., 2012, 13, 16295]. A quantum dynamical calculation of the evolution in the excited state shows an even larger probability of finding those two protons close to the center compared to in the ground state, but no formation of the proton-transferred product is observed. A QM/MM photoactive state geometry optimization, initiated using a configuration obtained by taking the A minimum and moving the protons to the product side, yields a minimum energy structure with the protons transferred and in which the His148 residue is substantially closer to the now anionic chromophore. These results indicate that: (1) proton transfer is not possible if structural relaxation of the surroundings of the chromophore is prevented; (2) protons H1 and H3 especially are found very close to the point halfway between the donor and acceptor after photoexcitation when the zero-point energy is considered; (3) a geometrical parameter exists (the His148-Cro distance) under which the structure with the protons transferred is not a minimum, and that, if included, should lead to the fluorescing I* structure. The existence of an oscillating stationary state between the reactants and products of the triple proton transfer reaction can explain the dual emission reported for the I0* intermediate of wtGFP.

Inverse Ubbelohde effect in the short hydrogen bond of photosystem II: Relation between H/D isotope effect and symmetry in potential energy profile

The short hydrogen bond between tyrosine Yz and D1-His190 of photosystem II (PSII) was investigated using multicomponent quantum mechanics, where the quantum fluctuation of a hydrogen nucleus was incorporated into electronic structure calculation. Our computation demonstrated that the deuteration for hydrogen in the short hydrogen bond of PSII led to the reduction of the O…N distance. It indicated an inverse Ubbelohde effect typically recognized in strong and symmetric hydrogen-bonding clusters such as FHF(-) and H3O2-. We confirmed that the relation between the geometric isotope effect and the symmetry of the potential energy profile of FHF(-) was reasonably agreed with that of PSII. According to this agreement, the short hydrogen bond in PSII can be regarded as a short strong hydrogen bond. © 2016 Wiley Periodicals, Inc.

The use of neutron scattering to determine the functional structure of glycoside hydrolase

Neutron diffraction provides different information from X-ray diffraction, because neutrons are scattered by atomic nuclei, whereas X-rays are scattered by electrons. One of the key advantages of neutron crystallography is the ability to visualize hydrogen and deuterium atoms, making it possible to observe the protonation state of amino acid residues, hydrogen bonds, networks of water molecules and proton relay pathways in enzymes. But, because of technical difficulties, less than 100 enzyme structures have been evaluated by neutron crystallography to date. In this review, we discuss the advantages and disadvantages of neutron crystallography as a tool to investigate the functional structure of glycoside hydrolases, with some examples.

Femtosecond structural dynamics drives the trans/cis isomerization in photoactive yellow protein

A variety of organisms have evolved mechanisms to detect and respond to light, in which the response is mediated by protein structural changes after photon absorption. The initial step is often the photoisomerization of a conjugated chromophore. Isomerization occurs on ultrafast time scales and is substantially influenced by the chromophore environment. Here we identify structural changes associated with the earliest steps in the trans-to-cis isomerization of the chromophore in photoactive yellow protein. Femtosecond hard x-ray pulses emitted by the Linac Coherent Light Source were used to conduct time-resolved serial femtosecond crystallography on photoactive yellow protein microcrystals over a time range from 100 femtoseconds to 3 picoseconds to determine the structural dynamics of the photoisomerization reaction.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Conversion of light-energy into molecular strain in the photocycle of the photoactive yellow protein

The Photoactive Yellow Protein (PYP) converts light energy into molecular strain, stored in the early pR₀-photocycle intermediate.

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

Photoactive yellow protein (PYP) has a characteristic hydrogen bond (H bond) between p-coumaric acid chromophore and Glu46, whose OH bond length has been observed to be 1.21 Å by the neutron diffraction technique [Proc. Natl. Acad. Sci. 106, 440-4]. Although it has been expected that such a drastic elongation of the OH bond could be caused by the quantum effect of the hydrogen nucleus, previous theoretical computations including the nuclear quantum effect have so far underestimated the bond length by more than 0.07 Å. To elucidate the origin of the difference, we performed a vibrational analysis of the H bond on potential energy curve with O…O distance of 2.47 Å on the equilibrium structure, and that with O…O distance of 2.56 Å on the experimental crystal structure. While the vibrationally averaged OH bond length for equilibrium structure was underestimated, the corresponding value for crystal structure was in reasonable agreement with the corresponding experimental values. The elongation of the O…O distance by the quantum mechanical or thermal fluctuation would be indispensable for the formation of a low-barrier hydrogen bond in PYP.

Theoretical analysis of geometry and NMR isotope shift in hydrogen-bonding center of photoactive yellow protein by combination of multicomponent quantum mechanics and ONIOM scheme

Multicomponent quantum mechanical (MC_QM) calculation has been extended with ONIOM (our own N-layered integrated molecular orbital + molecular mechanics) scheme [ONIOM(MC_QM:MM)] to take account of both the nuclear quantum effect and the surrounding environment effect. The authors have demonstrated the first implementation and application of ONIOM(MC_QM:MM) method for the analysis of the geometry and the isotope shift in hydrogen-bonding center of photoactive yellow protein. ONIOM(MC_QM:MM) calculation for a model with deprotonated Arg52 reproduced the elongation of O-H bond of Glu46 observed by neutron diffraction crystallography. Among the unique isotope shifts in different conditions, the model with protonated Arg52 with solvent effect reasonably provided the best agreement with the corresponding experimental values from liquid NMR measurement. Our results implied the availability of ONIOM(MC_QM:MM) to distinguish the local environment around hydrogen bonds in a biomolecule.

Influence of a chromophore analogue in the protein cage of a photoactive yellow protein

Time-resolved spectra of a photoactive yellow protein (PYP) containing cyano-p-coumaric acid (CHCA) were recorded. To understand the mechanism of photo-isomerization, an electron-withdrawing CN group was introduced into the PYP to alter the C[double bond, length as m-dash]C double bond character. Free CHCA chromophores in aqueous solution underwent photo-isomerization whereas PYP with a bound CHCA (PYP-CN) exhibited no photocycle at acidic or alkaline pH or in urea and other solutions. Furthermore, no photocycle was observed with PYP mutants after illumination. This phenomenon cannot be fully explained by the electron-withdrawing properties of the CN group. We conclude that the CHCA chromophore in PYP was locked in the protein cage and that the CN group interacted with the protein residues.

(15)N and (1)H Solid-State NMR Investigation of a Canonical Low-Barrier Hydrogen-Bond Compound: 1,8-Bis(dimethylamino)naphthalene

Strong or low-barrier hydrogen bonds have often been proposed in proteins to explain enzyme catalysis and proton-transfer reactions. So far (1)H chemical shifts and scalar couplings have been used as the main NMR spectroscopic signatures for strong H-bonds. In this work, we report simultaneous measurements of (15)N and (1)H chemical shifts and N-H bond lengths by solid-state NMR in (15)N-labeled 1,8-bis(dimethylamino)naphthalene (DMAN), which contains a well-known strong NHN H-bond. We complexed DMAN with three different counteranions to examine the effects of the chemical environment on the H-bond lengths and chemical shifts. All three DMAN compounds exhibit significantly elongated N-H distances compared to the covalent bond length, and the (1)H(N) chemical shifts are larger than ∼17 ppm, consistent with strong NHN H-bonds in the DMAN cation. However, the (15)N and (1)H chemical shifts and the precise N-H distances differ among the three compounds, and the (15)N chemical shifts show opposite dependences on the proton localization from the general trend in organic compounds, indicating the significant effects of the counteranions on the electronic structure of the H-bond. These data provide useful NMR benchmarks for strong H-bonds and caution against the sole reliance on chemical shifts for identifying strong H-bonds in proteins since neighboring side chains can exert influences on chemical shifts similar to those of the bulky organic anions in DMAN. Instead, N-H bond lengths should be measured, in conjunction with chemical shifts, as a more fundamental parameter of H-bond strength.

Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography

Ligand binding can change the pKa of protein residues and influence enzyme catalysis. Herein, we report three ultrahigh resolution X-ray crystal structures of CTX-M β-lactamase, directly visualizing protonation state changes along the enzymatic pathway: apo protein at 0.79 Å, precovalent complex with nonelectrophilic ligand at 0.89 Å, and acylation transition state (TS) analogue at 0.84 Å. Binding of the noncovalent ligand induces a proton transfer from the catalytic Ser70 to the negatively charged Glu166, and the formation of a low-barrier hydrogen bond (LBHB) between Ser70 and Lys73, with a length of 2.53 Å and the shared hydrogen equidistant from the heteroatoms. QM/MM reaction path calculations determined the proton transfer barrier to be 1.53 kcal/mol. The LBHB is absent in the other two structures although Glu166 remains neutral in the covalent complex. Our data represents the first X-ray crystallographic example of a hydrogen engaged in an enzymatic LBHB, and demonstrates that desolvation of the active site by ligand binding can provide a protein microenvironment conducive to LBHB formation. It also suggests that LBHBs may contribute to stabilization of the TS in general acid/base catalysis together with other preorganized features of enzyme active sites. These structures reconcile previous experimental results suggesting alternatively Glu166 or Lys73 as the general base for acylation, and underline the importance of considering residue protonation state change when modeling protein-ligand interactions. Additionally, the observation of another LBHB (2.47 Å) between two conserved residues, Asp233 and Asp246, suggests that LBHBs may potentially play a special structural role in proteins.

New insights into the enzymatic mechanism of human chitotriosidase (CHIT1) catalytic domain by atomic resolution X-ray diffraction and hybrid QM/MM

Chitotriosidase (CHIT1) is a human chitinase belonging to the highly conserved glycosyl hydrolase family 18 (GH18). GH18 enzymes hydrolyze chitin, an N-acetylglucosamine polymer synthesized by lower organisms for structural purposes. Recently, CHIT1 has attracted attention owing to its upregulation in immune-system disorders and as a marker of Gaucher disease. The 39 kDa catalytic domain shows a conserved cluster of three acidic residues, Glu140, Asp138 and Asp136, involved in the hydrolysis reaction. Under an excess concentration of substrate, CHIT1 and other homologues perform an additional activity, transglycosylation. To understand the catalytic mechanism of GH18 chitinases and the dual enzymatic activity, the structure and mechanism of CHIT1 were analyzed in detail. The resolution of the crystals of the catalytic domain was improved from 1.65 Å (PDB entry 1waw) to 0.95-1.10 Å for the apo and pseudo-apo forms and the complex with chitobiose, allowing the determination of the protonation states within the active site. This information was extended by hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. The results suggest a new mechanism involving changes in the conformation and protonation state of the catalytic triad, as well as a new role for Tyr27, providing new insights into the hydrolysis and transglycosylation activities.

Short Hydrogen Bonds and Proton Delocalization in Green Fluorescent Protein (GFP)

Short hydrogen bonds and specifically low-barrier hydrogen bonds (LBHBs) have been the focus of much attention and controversy for their possible role in enzymatic catalysis. The green fluorescent protein (GFP) mutant S65T, H148D has been found to form a very short hydrogen bond between Asp148 and the chromophore resulting in significant spectral perturbations. Leveraging the unique autocatalytically formed chromophore and its sensitivity to this interaction we explore the consequences of proton affinity matching across this putative LBHB. Through the use of noncanonical amino acids introduced through nonsense suppression or global incorporation, we systematically modify the acidity of the GFP chromophore with halogen substituents. X-ray crystal structures indicated that the length of the interaction with Asp148 is unchanged at ∼2.45 Å while the absorbance spectra demonstrate an unprecedented degree of color tuning with increasing acidity. We utilized spectral isotope effects, isotope fractionation factors, and a simple 1D model of the hydrogen bond coordinate in order to gain insight into the potential energy surface and particularly the role that proton delocalization may play in this putative short hydrogen bond. The data and model suggest that even with the short donor-acceptor distance (∼2.45 Å) and near perfect affinity matching there is not a LBHB, that is, the barrier to proton transfer exceeds the H zero-point energy.

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.