Direct determination of protonation states and visualization of hydrogen bonding in a glycoside hydrolase with neutron crystallography

Development of Tunable Mechanism-Based Carbasugar Ligands that Stabilize Glycoside Hydrolases through the Formation of Transient Covalent Intermediates

Mutations in many members of the set of human lysosomal glycoside hydrolases cause a wide range of lysosomal storage diseases. As a result, much effort has been directed toward identifying pharmacological chaperones of these lysosomal enzymes. The majority of the candidate chaperones are active site-directed competitive iminosugar inhibitors but these have met with limited success. As a first step toward an alternative class of pharmacological chaperones we explored the potential of small molecule mechanism-based reversible covalent inhibitors to form transient enzyme-inhibitor adducts. By serial synthesis and kinetic analysis of candidate molecules, we show that rational tuning of the chemical reactivity of glucose-configured carbasugars delivers cyclohexenyl-based allylic carbasugar that react with the lysosomal enzyme β-glucocerebrosidase (GCase) to form covalent enzyme-adducts with different half-lives. X-ray structural analysis of these compounds bound noncovalently to GCase, along with the structures of the covalent adducts of compounds that reacted with the catalytic nucleophile of GCase, reveal unexpected reactivities of these compounds. Using differential scanning fluorimetry, we show that formation of a transient covalent intermediate stabilizes the folded enzyme against thermal denaturation. In addition, these covalent adducts break down to liberate the active enzyme and a product that is no longer inhibitory. We further show that the one compound, which reacts through an unprecedented SN1'-like mechanism, exhibits exceptional reactivity-illustrated by this compound also covalently labeling an α-glucosidase. We anticipate that such carbasugar-based single turnover covalent ligands may serve as pharmacological chaperones for lysosomal glycoside hydrolases and other disease-associated retaining glycosidases. The unusual reactivity of these molecules should also open the door to creation of new chemical biology probes to explore the biology of this important superfamily of glycoside hydrolases.

Temperature-dependent structural changes in xylanase II from Trichoderma longibrachiatum

Endo-β-1,4-xylanases degrade heteroxylans that constitute the lignocellulosic plant cell wall. This enzyme is widely used in the food, paper, textile, and biorefinery industries. Temperature affects the optimum activity of xylanase and is an important factor in its application. Various structural analyses of xylanase have been performed, but its structural influence by temperature is not fully elucidated. To better understand the structural influence of xylanase due to temperature, the crystal structure of xylanase II from Trichoderma longibrachiatum (TloXynII) at room and cryogenic temperatures was determined at 2.1 and 1.9 Å resolution, respectively. The room-temperature structure of TloXynII (TloXynIIRT) showed a B-factor value 2.09 times higher than that of the cryogenic-temperature structure of TloXynII (TloXynIICryo). Subtle movement of the catalytic and substrate binding residues was observed between TloXynIIRT and TloXynIICryo. In TloXynIIRT, the thumb domain exhibited high flexibility, whereas in TloXynIICryo, the finger domain exhibited high flexibility. The substrate binding cleft of TloXynIIRT was narrower than that of TloXynIICryo, indicating a distinct finger domain conformation. Numerous water molecule networks were observed in the substrate binding cleft of TloXynIICryo, whereas only a few water molecules were observed in TloXynIIRT. These structural analyses expand our understanding of the temperature-dependent conformational changes in xylanase.

Experimental Evidence for the Role of Dynamics in pH-Dependent Enzymatic Activity

Enzymatic activity is heavily influenced by pH, but the rationale for the dynamical mechanism of pH-dependent enzymatic activity has not been fully understood. In this work, combined neutron scattering techniques, including quasielastic neutron scattering (QENS) and small angle neutron scattering (SANS), are used to study the structural and dynamic changes of a model enzyme, xylanase, under different pH and temperature environments. The QENS results reveal that xylanase at optimal pH exhibits faster relaxational dynamics and a lower energy barrier between conformational substates. The SANS results demonstrate that pH affects both xylanase's stability and monodispersity. Our findings indicate that enzymes have optimized stability and function under their optimal pH conditions, with both structure and dynamics being affected. The current study offers valuable insights into enzymatic functionality mechanisms, allowing for broad industrial applications.

Neutron diffraction from a microgravity-grown crystal reveals the active site hydrogens of the internal aldimine form of tryptophan synthase

Pyridoxal 5'-phosphate (PLP), the biologically active form of vitamin B6, is an essential cofactor in many biosynthetic pathways. The emergence of PLP-dependent enzymes as drug targets and biocatalysts, such as tryptophan synthase (TS), has underlined the demand to understand PLP-dependent catalysis and reaction specificity. The ability of neutron diffraction to resolve the positions of hydrogen atoms makes it an ideal technique to understand how the electrostatic environment and selective protonation of PLP regulates PLP-dependent activities. Facilitated by microgravity crystallization of TS with the Toledo Crystallization Box, we report the 2.1 Å joint X-ray/neutron (XN) structure of TS with PLP in the internal aldimine form. Positions of hydrogens were directly determined in both the α- and β-active sites, including PLP cofactor. The joint XN structure thus provides insight into the selective protonation of the internal aldimine and the electrostatic environment of TS necessary to understand the overall catalytic mechanism.

Protomer of Imipramine Captured in Cucurbit[7]uril

Small molecules possessing multiple proton-accessible sites are important not only to many biological systems but also to host-guest chemistry; their protonation states are causal to boosting or hindering specific host-guest interactions. However, determining the protonation site is not trivial. Herein, we conduct electrospray ionization ion mobility spectrometry-mass spectrometry to imipramine, a known molecule with two protonation sites, based on the introduction of cucurbit[7]uril as a host molecule. For protonated imipramine, the proposed strategy allows clear distinction of the two protomers as host-guest complex ions and can be leveraged to capture the energetically less preferable protomer of the protonated imipramine.

Neutron crystallographic refinement with REFMAC5 from the CCP4 suite

Hydrogen (H) atoms are abundant in macromolecules and often play critical roles in enzyme catalysis, ligand-recognition processes and protein-protein interactions. However, their direct visualization by diffraction techniques is challenging. Macromolecular X-ray crystallography affords the localization of only the most ordered H atoms at (sub-)atomic resolution (around 1.2 Å or higher). However, many H atoms of biochemical significance remain undetectable by this method. In contrast, neutron diffraction methods enable the visualization of most H atoms, typically in the form of deuterium (2H) atoms, at much more common resolution values (better than 2.5 Å). Thus, neutron crystallography, although technically demanding, is often the method of choice when direct information on protonation states is sought. REFMAC5 from the Collaborative Computational Project No. 4 (CCP4) is a program for the refinement of macromolecular models against X-ray crystallographic and cryo-EM data. This contribution describes its extension to include the refinement of structural models obtained from neutron crystallographic data. Stereochemical restraints with accurate bond distances between H atoms and their parent atom nuclei are now part of the CCP4 Monomer Library, the source of prior chemical information used in the refinement. One new feature for neutron data analysis in REFMAC5 is refinement of the protium/deuterium (1H/2H) fraction. This parameter describes the relative 1H/2H contribution to neutron scattering for hydrogen isotopes. The newly developed REFMAC5 algorithms were tested by performing the (re-)refinement of several entries available in the PDB and of one novel structure (FutA) using either (i) neutron data only or (ii) neutron data supplemented by external restraints to a reference X-ray crystallographic structure. Re-refinement with REFMAC5 afforded models characterized by R-factor values that are consistent with, and in some cases better than, the originally deposited values. The use of external reference structure restraints during refinement has been observed to be a valuable strategy, especially for structures at medium-low resolution.

Investigating the effect of substrate binding on the catalytic activity of xylanase

XynAF1 from Aspergillus fumigatus Z5 is an efficient thermophilic xylanase belonging to glycoside hydrolase family 10 (GH10). The non-catalytic amino acids N179 and R246 in its catalytic center formed one and three intermolecular H-bonds with the substrate in the aglycone region, respectively. Here we purified XynAF1-N179S and XynAF1-R246K, and obtained the protein-product complex structures by X-ray diffraction. The snapshots indicated that mutations at N179 and R246 had decreased the substrate-binding ability in the aglycone region. XynAF1-N179S, XynAF1-R246K, and XynAF1-N179S-R246K lost one, three, and four H-bonds with the substrate in comparison with the wild-type XynAF1, respectively, but this had little influence on the protein structure. As expected, N179S, R246K, and N179S-R246K led to a gradual decrease of substrate affinity of XynAF1. Interestingly, the enzyme assay showed that N179S increased catalytic efficiency, while both R246K and N179S-R246K had decreased catalytic efficiency. KEY POINTS: • The non-catalytic amino acids of XynAF1 could form H-bonds with the substrate. • The protein-product complex structures were obtained by X-ray diffraction. • The enzyme-substrate-binding capacity could affect enzyme catalytic efficiency.

Description of peptide bond planarity from high-resolution neutron crystallography

Neutron crystallography is a highly effective method for visualizing hydrogen atoms in proteins. In our recent study, we successfully determined the high-resolution (1.2 Å) neutron structure of high-potential iron-sulfur protein, refining the coordinates of some amide protons without any geometric restraints. Interestingly, we observed that amide protons are deviated from the peptide plane due to electrostatic interactions. Moreover, the difference in the position of the amide proton of Cys75 between reduced and oxidized states is possibly attributed to the electron storage capacity of the iron-sulfur cluster. Additionally, we have discussed about the rigidity of the iron-sulfur cluster based on the results of the hydrogen-deuterium exchange. Our research underscores the significance of neutron crystallography in protein structure elucidation, enriching our understanding of protein functions at an atomic resolution.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Structural and biophysical characterization of the multidomain xylanase Xyl

The depletion of fossil fuels, associated pollution, and resulting health hazards are of concern worldwide. Woody biomass constitutes an alternative source of cleaner and renewable energy. The efficient use of woody biomass depends on xylan depolymerisation as the endo-β-1,4-xylopyranosyl homopolymer is the main component of hemicellulose, the second most abundant component of wood. Xylan depolymerisation is achieved by hemicellulolytic xylanases of glycoside hydrolase (GH) families 5, 8, 10, 11, 30 and 43 of the CAZY database. We analysed a multidomain xylanase (Xyl) from the hindgut metagenome of the snouted harvester termite Trinervitermes trinervoides that releases xylobiose and xylotriose from beech and birch xylan and wheat arabinoxylan. The four domains of Xyl include an N-terminal GH11 xylanase domain, two family 36-like carbohydrate-binding domains CBM36-1 and 2, and a C-terminal CE4 esterase domain. Previous analyses indicated that CBM36-1 deletion slightly increased GH11 catalysis at low pH whereas removal of both CBMs decreased xylanase activity at 60°C from 90 to 56%. Possible cooperativity between the domains suggested by these observations was explored. A crystal structure of the two-domain construct, GH11-CBM36-1, confirmed the structure of the GH11 domain whereas the CBM36-1 domain lacked electron density, possibly indicating a random orientation of the CBM36-1 domain around the GH11 domain. Isothermal titration calorimetry (ITC) experiments similarly did not indicate specific interactions between the individual domains of Xyl supporting a "beads-on-a-string" model for Xyl domains.

Engineering the optimum pH of β-galactosidase from Aspergillus oryzae for efficient hydrolysis of lactose

β-Galactosidase (lacA) from Aspergillus oryzae is widely used in the dairy industry. Its acidic pH optimum and severe product inhibition limit its application for lactose hydrolysis in milk. In the present study, structure-based sequence alignment was conducted to determine the candidate mutations to shift the pH optimum of lacA to the neutral range. The Y138F and Y364F mutants shifted the pH optimum of lacA from 4.5 to 5.5 and 6.0, respectively. The acid dissociation constant (pKa) values of catalytic acid/base residues with upwards shift were consistent with the increased pH optimum. All variants in the present study also alleviated galactose inhibition to various extents. Molecular dynamics demonstrated that the less rigid tertiary structures and lower galactose-binding free energy of Y138F and Y364F might facilitate the release of the end product. Both Y138F and Y364F mutants exhibited better hydrolytic ability than lacA in milk lactose hydrolysis. The higher pH optimum and lower galactose inhibition of Y138F and Y364F may explain their superiority over lacA. The Y138F and Y364F mutants in the present study showed potential in producing low-lactose milk, and our studies provide a novel strategy for engineering the pH optimum of glycoside hydrolase.

Protonation State of an Important Histidine from High Resolution Structures of Lytic Polysaccharide Monooxygenases

Lytic Polysaccharide Monooxygenases (LPMOs) oxidatively cleave recalcitrant polysaccharides. The mechanism involves (i) reduction of the Cu, (ii) polysaccharide binding, (iii) binding of different oxygen species, and (iv) glycosidic bond cleavage. However, the complete mechanism is poorly understood and may vary across different families and even within the same family. Here, we have investigated the protonation state of a secondary co-ordination sphere histidine, conserved across AA9 family LPMOs that has previously been proposed to be a potential proton donor. Partial unrestrained refinement of newly obtained higher resolution data for two AA9 LPMOs and re-refinement of four additional data sets deposited in the PDB were carried out, where the His was refined without restraints, followed by measurements of the His ring geometrical parameters. This allowed reliable assignment of the protonation state, as also validated by following the same procedure for the His brace, for which the protonation state is predictable. The study shows that this histidine is generally singly protonated at the Nε2 atom, which is close to the oxygen species binding site. Our results indicate robustness of the method. In view of this and other emerging evidence, a role as proton donor during catalysis is unlikely for this His.

The human gut symbiont Ruminococcus gnavus shows specificity to blood group A antigen during mucin glycan foraging: Implication for niche colonisation in the gastrointestinal tract

The human gut symbiont Ruminococcus gnavus displays strain-specific repertoires of glycoside hydrolases (GHs) contributing to its spatial location in the gut. Sequence similarity network analysis identified strain-specific differences in blood-group endo-β-1,4-galactosidase belonging to the GH98 family. We determined the substrate and linkage specificities of GH98 from R. gnavus ATCC 29149, RgGH98, against a range of defined oligosaccharides and glycoconjugates including mucin. We showed by HPAEC-PAD and LC-FD-MS/MS that RgGH98 is specific for blood group A tetrasaccharide type II (BgA II). Isothermal titration calorimetry (ITC) and saturation transfer difference (STD) NMR confirmed RgGH98 affinity for blood group A over blood group B and H antigens. The molecular basis of RgGH98 strict specificity was further investigated using a combination of glycan microarrays, site-directed mutagenesis, and X-ray crystallography. The crystal structures of RgGH98 in complex with BgA trisaccharide (BgAtri) and of RgGH98 E411A with BgA II revealed a dedicated hydrogen network of residues, which were shown by site-directed mutagenesis to be critical to the recognition of the BgA epitope. We demonstrated experimentally that RgGH98 is part of an operon of 10 genes that is overexpresssed in vitro when R. gnavus ATCC 29149 is grown on mucin as sole carbon source as shown by RNAseq analysis and RT-qPCR confirmed RgGH98 expression on BgA II growth. Using MALDI-ToF MS, we showed that RgGH98 releases BgAtri from mucin and that pretreatment of mucin with RgGH98 confered R. gnavus E1 the ability to grow, by enabling the E1 strain to metabolise BgAtri and access the underlying mucin glycan chain. These data further support that the GH repertoire of R. gnavus strains enable them to colonise different nutritional niches in the human gut and has potential applications in diagnostic and therapeutics against infection.

Differences in Gluco and Galacto Substrate-Binding Interactions in a Dual 6Pβ-Glucosidase/6Pβ-Galactosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis

Bacterial glycoside hydrolase 1 (GH1) enzymes with 6-phospho-β-galactosidase and 6-phospho-β-glucosidase activities have the important task of releasing phosphorylated and nonphosphorylated monosaccharides into the cytoplasm. Curiously, dual 6-phospho-β-galactosidase/6-phospho-β-glucosidase (dual-phospho) enzymes have broad specificity and are able to hydrolyze galacto- and gluco-derived substrates. This study investigates the structure and substrate specificity of a GH family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglC. The enzyme structure has been solved, and sequence analysis, molecular dynamics simulations, and binding free energy calculations offered evidence of dual-phospho activity. Both test ligands p-nitrophenyl-β-d-galactoside-6-phosphate (PNP6Pgal) and p-nitrophenyl-β-d-glucoside-6-phosphate (PNP6Pglc) demonstrated strong binding to BlBglC although the pose and interactions of the PNP6Pglc triplicates were slightly more consistent. Interestingly, known specificity-inducing residues, Gln23 and Trp433, bind strongly to the ligand O3 hydroxyl group in the PNP6Pgal-BlBglC complex and to the ligand O4 hydroxyl group in the PNP6Pglc-BlBglC complex. Additionally, the BlBglC-His124 residue is a major contributor of hydrogen bonds to the PNP6Pgal O3 hydroxyl group but does not form any hydrogen bonds with PNP6Pglc. On the other hand, BlBglC residues Tyr173, Tyr301, Gln302, and Thr321 form hydrogen bonds with PNP6Pglc but not PNP6Pgal. These findings provide important details of the broad specificity of dual-phospho activity GH1 enzymes.

Capturing the Catalytic Proton of Dihydrofolate Reductase: Implications for General Acid-Base Catalysis

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H₃O⁺) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

Configurational Entropy of Folded Proteins and Its Importance for Intrinsically Disordered Proteins

Many pairwise additive force fields are in active use for intrinsically disordered proteins (IDPs) and regions (IDRs), some of which modify energetic terms to improve the description of IDPs/IDRs but are largely in disagreement with solution experiments for the disordered states. This work considers a new direction-the connection to configurational entropy-and how it might change the nature of our understanding of protein force field development to equally well encompass globular proteins, IDRs/IDPs, and disorder-to-order transitions. We have evaluated representative pairwise and many-body protein and water force fields against experimental data on representative IDPs and IDRs, a peptide that undergoes a disorder-to-order transition, for seven globular proteins ranging in size from 130 to 266 amino acids. We find that force fields with the largest statistical fluctuations consistent with the radius of gyration and universal Lindemann values for folded states simultaneously better describe IDPs and IDRs and disorder-to-order transitions. Hence, the crux of what a force field should exhibit to well describe IDRs/IDPs is not just the balance between protein and water energetics but the balance between energetic effects and configurational entropy of folded states of globular proteins.

Visualization of hydrogen atoms in a perdeuterated lectin-fucose complex reveals key details of protein-carbohydrate interactions

Carbohydrate-binding proteins from pathogenic bacteria and fungi have been shown to be implicated in various pathological processes, where they interact with glycans present on the surface of the host cells. These interactions are part of the initial processes of infection of the host and are very important to study at the atomic level. Here, we report the room temperature neutron structures of PLL lectin from Photorhabdus laumondii in its apo form and in complex with deuterated L-fucose, which is, to our knowledge, the first neutron structure of a carbohydrate-binding protein in complex with a fully deuterated carbohydrate ligand. A detailed structural analysis of the lectin-carbohydrate interactions provides information on the hydrogen bond network, the role of water molecules, and the extent of the CH-π stacking interactions between fucose and the aromatic amino acids in the binding site.

Active site dynamics and catalytic mechanism in arabinan hydrolysis catalyzed by GH43 endo-arabinanase from QM/MM molecular dynamics simulation and potential energy surface

The endo-1,5-α-L-arabinanases, belonging to glycoside hydrolase family 43 (GH43), catalyse the hydrolysis of α-1,5-arabinofuranosidic bonds in arabinose-containing polysaccharides. These enzymes are proposed targets for industrial and medical applications. Here, molecular dynamics (MD), potential energy surface and free energy (potential of mean force) simulations are undertaken using hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to understand the active site dynamics, catalytic mechanism and the electrostatic influence of active site residues of the GH43 endo-arabinanase from G. stearothermophilus. The calculated results give support to the single-displacement mechanism proposed for the inverting GH43 enzymes: first a proton is transferred from the general acid E201 to the substrate, followed by a nucleophilic attack by water, activated by the general base D27, on the anomer carbon. A conformational change (²E ↔E₃ ↔ ⁴E) in the -1 sugar ring is observed involving a transition state featuring an oxocarbenium ion character. Residues D87, K106, H271 are highlighted as potential targets for future mutation experiments in order to increase the efficiency of the reaction. To our knowledge, this is the first QM/MM study providing molecular insights into the glycosidic bond hydrolysis of a furanoside substrate by an inverting GH in solution.Communicated by Ramaswamy H. Sarma.

X-ray Structure, Bioinformatics Analysis, and Substrate Specificity of a 6-Phospho-β-glucosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis

In bacteria, mono- and disaccharides are phosphorylated during the uptake processes through the vastly spread transport system phosphoenolpyruvate-dependent phosphotransferase. As an initial step in the phosphorylated disaccharide metabolism pathway, 6-phospho-β-glucosidases and 6-phospho-β-galactosidases play a crucial role by releasing phosphorylated and nonphosphorylated monosaccharides. However, structural determinants for the specificity of these enzymes still need to be clarified. Here, an X-ray structure of a glycoside hydrolase family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglH, was determined at 2.2 Å resolution, and its substrate specificity was investigated. The sequence of BlBglH was compared to the sequences of 58 other GH1 enzymes using sequence alignments, sequence identity calculations, phylogenetic analysis, and motif discovery. Through these various analyses, BlBglH was found to have sequence features characteristic of the 6-phospho-β-glucosidase activity enzymes. Motif and structural observations highlighted the importance of loop L8 in 6-phospho-β-glucosidase activity enzymes. To further affirm enzyme specificity, molecular docking and molecular dynamics simulations were performed using the crystallographic structure of BlBglH. Docking was carried out with a 6-phospho-β-glucosidase enzyme activity positive and negative control ligand, followed by 400 ns of MD simulations. The positive and negative control ligands were PNP6Pglc and PNP6Pgal, respectively. PNP6Pglc maintained favorable interactions within the active site until the end of the MD simulation, while PNP6Pgal exhibited instability. The favorable binding of substrate stabilized the loops that surround the active site. Binding free energy calculations showed that the PNP6Pglc complex had a substantially lower binding energy compared to the PNP6Pgal complex. Altogether, the findings of this study suggest that BlBglH possesses 6-phospho-β-glucosidase enzymatic activity and revealed sequence and structural differences between bacterial GH1 enzymes of various activities.

Studying the Role of a Single Mutation of a Family 11 Glycoside Hydrolase Using High-Resolution X-ray Crystallography

XynII is a family 11 glycoside hydrolase that uses the retaining mechanism for catalysis. In the active site, E177 works as the acid/base and E86 works as the nucleophile. Mutating an uncharged residue (N44) to an acidic residue (D) near E177 decreases the enzyme's optimal pH by ~ 1.0 unit. D44 was previously suggested to be a second proton carrier for catalysis. To test this hypothesis, we abolished the activity of E177 by mutating it to be Q, and mutated N44 to be D or E. These double mutants have dramatically decreased activities. Our high-resolution crystallographic structures and the microscopic pK_a calculations show that D44 has similar position and pK_a value during catalysis, indicating that D44 changes electrostatics around E177, which makes it prone to rotate as the acid/base in acidic conditions, thus decreases the pH optimum. Our results could be helpful to design enzymes with different pH optimum.

Insight into the Catalytic Mechanism of GH11 Xylanase: Computational Analysis of Substrate Distortion Based on a Neutron Structure

The reaction mechanism of biomass decomposition by xylanases remains the subject of debate. To clarify the mechanism we investigated the glycosylation step of GH11 xylanase, an enzyme that catalyzes the hydrolysis of lignocellulosic hemicellulose (xylan). Making use of a recent neutron crystal structure, which revealed the protonation states of relevant residues, we used ab initio quantum mechanics/molecular mechanics (QM/MM) calculations to determine the detailed reaction mechanism of the glycosylation step. In particular, our focus is on the controversial question of whether or not an oxocarbenium ion intermediate is formed on the reaction pathway. The calculations support the validity of a basic retaining mechanism within a double-displacement scheme. The estimated free energy barrier of this reaction is ∼18 kcal/mol with QM/MM-CCSD(T)/6-31(+)G**//MP2/6-31+G**/AMBER calculations, and the rate-determining step of the glycosylation is scission of the glycosidic bond after proton transfer from the acidic Glu177. The estimated lifetime of the oxocarbenium ion intermediate (on the order of tens of ps) and the secondary kinetic isotope effect suggest that there is no accumulation of this intermediate on the reaction path, although the intermediate can be transiently formed. In the enzyme-substrate (ES) complex, the carbohydrate structure of the xylose residue at the -1 subsite has a rather distorted (skewed) geometry, and this xylose unit at the active site has an apparent half-chair conformation when the oxocarbenium ion intermediate is formed. The major catalytic role of the protein environment is to orient residues that take part in the initial proton transfer. Because of a fine alignment of catalytic residues, the enzyme can accelerate the glycosylation reaction without paying a reorganization energy penalty.

Maize WI5 encodes an endo-1,4-β-xylanase required for secondary cell wall synthesis and water transport in xylem

Water transport from roots to leaves through xylem is important for plant growth and development. Defects in water transport can cause drought stress, even when there is adequate water in the soil. Here, we identified the maize (Zea mays) wilty5 (wi5) mutant, which exhibits marked dwarfing and leaf wilting throughout most of its life cycle under normal growth conditions. wilty5 seedlings exhibited lower xylem conductivity and wilted more rapidly under drought, NaCl, and high temperature treatments than wild-type plants. Map-based cloning revealed that WI5 encodes an active endo-1,4-β-xylanase from glycosyl dehydration family 10, which mainly functions in degrading and reorganizing cell wall xylan. Reverse-transcription polymerase chain reaction and β-glucuronidase assays revealed that WI5 is highly expressed in stems, especially in internodes undergoing secondary wall assembly. RNA sequencing suggested that WI5 plays a unique role in internode growth. Immunohistochemistry and electron microscopy confirmed that wi5 is defective in xylan deposition and secondary cell wall thickening. Lignin deposition and xylan content were markedly reduced in wi5 compared to the wild-type plants. Our results suggest that WI5 functions in xylem cell wall thickening through its xylanase activity and thereby regulates xylem water transport, the drought stress response, and plant growth in maize.

Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D2O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

Titration of ionizable groups in proteins using multiple neutron data sets from a single crystal: application to the small GTPase Ras

Neutron protein crystallography (NPC) reveals the three-dimensional structures of proteins, including the positions of H atoms. The technique is particularly suited to elucidate ambiguous catalytic steps in complex biochemical reactions. While NPC uniquely complements biochemical assays and X-ray structural analyses by revealing the protonation states of ionizable groups at and around the active site of enzymes, the technique suffers from a major drawback: large single crystals must be grown to compensate for the relatively low flux of neutron beams. However, in addition to revealing the positions of hydrogens involved in enzyme catalysis, NPC has the advantage over X-ray crystallography that the crystals do not suffer radiation damage. The lack of radiation damage can be exploited to conduct in crystallo parametric studies. Here, the use of a single crystal of the small GTPase Ras to collect three neutron data sets at pD 8.4, 9.0 and 9.4 is reported, enabling an in crystallo titration study using NPC. In addition to revealing the behavior of titratable groups in the active site, the data sets will allow the analysis of allosteric water-mediated communication networks across the molecule, particularly regarding Cys118 and three tyrosine residues central to these networks, Tyr32, Tyr96 and Tyr137, with pKa values expected to be in the range sampled in our experiments.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Neutron and X-ray crystal structures of Lactobacillus brevis alcohol dehydrogenase reveal new insights into hydrogen-bonding pathways

Lactobacillus brevis alcohol dehydrogenase (LbADH) is a well studied homotetrameric enzyme which catalyzes the enantioselective reduction of prochiral ketones to the corresponding secondary alcohols. LbADH is stable and enzymatically active at elevated temperatures and accepts a broad range of substrates, making it a valuable tool in industrial biocatalysis. Here, the expression, purification and crystallization of LbADH to generate large, single crystals with a volume of up to 1 mm3 suitable for neutron diffraction studies are described. Neutron diffraction data were collected from an H/D-exchanged LbADH crystal using the BIODIFF instrument at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany to a resolution dmin of 2.15 Å in 16 days. This allowed the first neutron crystal structure of LbADH to be determined. The neutron structure revealed new details of the hydrogen-bonding network originating from the ion-binding site of LbADH and provided new insights into the reasons why divalent magnesium (Mg2+) or manganese (Mn2+) ions are necessary for its activity. X-ray diffraction data were obtained from the same crystal at the European Synchrotron Radiation Facility (ESRF), Grenoble, France to a resolution dmin of 1.48 Å. The high-resolution X-ray structure suggested partial occupancy of Mn2+ and Mg2+ at the ion-binding site. This is supported by the different binding affinity of Mn2+ and Mg2+ to the tetrameric structure calculated via free-energy molecular-dynamics simulations.

Improving the accuracy and resolution of neutron crystallographic data by three-dimensional profile fitting of Bragg peaks in reciprocal space

Neutron crystallography is a powerful technique for directly visualizing the locations of H atoms in biological macromolecules. This information has provided key new insights into enzyme mechanisms, ligand binding and hydration. However, despite the importance of this information, the application of neutron crystallography in biology has been limited by the relatively low flux of available neutron beams and the large incoherent neutron scattering from hydrogen, both of which contribute to weak diffraction data with relatively low signal-to-background ratios. A method has been developed to fit weak data based on three-dimensional profile fitting of Bragg peaks in reciprocal space by an Ikeda-Carpenter function with a bivariate Gaussian. When applied to data collected from three different proteins, three-dimensional profile fitting yields intensities with higher correlation coefficients (CC1/2) at high resolutions, decreased Rfree factors, extended resolutions and improved nuclear density maps. Importantly, additional features are revealed in nuclear density maps that may provide additional scientific information. These results suggest that three-dimensional profile fitting will help to extend the capabilities of neutron macromolecular crystallography.

The Neutron Macromolecular Crystallography Instruments at Oak Ridge National Laboratory: Advances, Challenges, and Opportunities

The IMAGINE and MaNDi instruments, located at Oak Ridge National Laboratory High Flux Isotope Reactor and Spallation Neutron Source, respectively, are powerful tools for determining the positions of hydrogen atoms in biological macromolecules and their ligands, orienting water molecules, and for differentiating chemical states in macromolecular structures. The possibility to model hydrogen and deuterium atoms in neutron structures arises from the strong interaction of neutrons with the nuclei of these isotopes. Positions can be unambiguously assigned from diffraction studies at the 1.5–2.5 Å resolutions, which are typical for protein crystals. Neutrons have the additional benefit for structural biology of not inducing radiation damage to protein crystals, which can be critical in the study of metalloproteins. Here we review the specifications of the IMAGINE and MaNDi beamlines and illustrate their complementarity. IMAGINE is suitable for crystals with unit cell edges up to 150 Å using a quasi-Laue technique, whereas MaNDi provides neutron crystallography resources for large unit cell samples with unit cell edges up to 300 Å using the time of flight (TOF) Laue technique. The microbial culture and crystal growth facilities which support the IMAGINE and MaNDi user programs are also described.

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.

Direct visualization of critical hydrogen atoms in a pyridoxal 5'-phosphate enzyme

Enzymes dependent on pyridoxal 5′-phosphate (PLP, the active form of vitamin B₆) perform a myriad of diverse chemical transformations. They promote various reactions by modulating the electronic states of PLP through weak interactions in the active site. Neutron crystallography has the unique ability of visualizing the nuclear positions of hydrogen atoms in macromolecules. Here we present a room-temperature neutron structure of a homodimeric PLP-dependent enzyme, aspartate aminotransferase, which was reacted in situ with α-methylaspartate. In one monomer, the PLP remained as an internal aldimine with a deprotonated Schiff base. In the second monomer, the external aldimine formed with the substrate analog. We observe a deuterium equidistant between the Schiff base and the C-terminal carboxylate of the substrate, a position indicative of a low-barrier hydrogen bond. Quantum chemical calculations and a low-pH room-temperature X-ray structure provide insight into the physical phenomena that control the electronic modulation in aspartate aminotransferase.

Pyridoxal 5’-phosphate (PLP) is a ubiquitous co factor for diverse enzymes, among them aspartate aminotransferase. Here the authors use neutron crystallography, which allows the visualization of the positions of hydrogen atoms, and computation to characterize the catalytic mechanism of the enzyme.

Protein crystallization and initial neutron diffraction studies of the photosystem II subunit PsbO

The PsbO protein of photosystem II stabilizes the active-site manganese cluster and is thought to act as a proton antenna. To enable neutron diffraction studies, crystals of the β-barrel core of PsbO were grown in capillaries. The crystals were optimized by screening additives in a counter-diffusion setup in which the protein and reservoir solutions were separated by a 1% agarose plug. Crystals were cross-linked with glutaraldehyde. Initial neutron diffraction data were collected from a 0.25 mm3 crystal at room temperature using the MaNDi single-crystal diffractometer at the Spallation Neutron Source, Oak Ridge National Laboratory.

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

Fifteen years of the Protein Crystallography Station: the coming of age of macromolecular neutron crystallography

The Protein Crystallography Station (PCS), located at the Los Alamos Neutron Scattering Center (LANSCE), was the first macromolecular crystallography beamline to be built at a spallation neutron source. Following testing and commissioning, the PCS user program was funded by the Biology and Environmental Research program of the Department of Energy Office of Science (DOE-OBER) for 13 years (2002-2014). The PCS remained the only dedicated macromolecular neutron crystallography station in North America until the construction and commissioning of the MaNDi and IMAGINE instruments at Oak Ridge National Laboratory, which started in 2012. The instrument produced a number of research and technical outcomes that have contributed to the field, clearly demonstrating the power of neutron crystallo-graphy in helping scientists to understand enzyme reaction mechanisms, hydrogen bonding and visualization of H-atom positions, which are critical to nearly all chemical reactions. During this period, neutron crystallography became a technique that increasingly gained traction, and became more integrated into macromolecular crystallography through software developments led by investigators at the PCS. This review highlights the contributions of the PCS to macromolecular neutron crystallography, and gives an overview of the history of neutron crystallography and the development of macromolecular neutron crystallography from the 1960s to the 1990s and onwards through the 2000s.

Structural and functional characterization of a highly stable endo-β-1,4-xylanase from Fusarium oxysporum and its development as an efficient immobilized biocatalyst

BACKGROUND

Replacing fossil fuel with renewable sources such as lignocellulosic biomass is currently a promising alternative for obtaining biofuel and for fighting against the consequences of climate change. However, the recalcitrant structure of lignocellulosic biomass residues constitutes a major limitation for its widespread use in industry. The efficient hydrolysis of lignocellulosic materials requires the complementary action of multiple enzymes including xylanases and β-xylosidases, which are responsible for cleaving exo- and endoxylan linkages, that release oligocarbohydrates that can be further processed by other enzymes.

RESULTS

We have identified the endo-β-1,4-xylanase Xyl2 from Fusarium oxysporum as a promising glycoside hydrolase family 11 enzyme for the industrial degradation of xylan. To characterize Xyl2, we have cloned the synthetic optimized gene and expressed and purified recombinant Xyl2 to homogeneity, finally obtaining 10 mg pure Xyl2 per liter of culture. The crystal structure of Xyl2 at 1.56 Å resolution and the structure of a methyl-xylopyranoside Xyl2 complex at 2.84 Å resolution cast a highly detailed view of the active site of the enzyme, revealing the molecular basis for the high catalytic efficiency of Xyl2. The kinetic analysis of Xyl2 demonstrates high xylanase activity and non-negligible β-xylosidase activity under a variety of experimental conditions including alkaline pH and elevated temperature. Immobilizing Xyl2 on a variety of solid supports enhances the enzymatic properties that render Xyl2 a promising industrial biocatalyst, which, together with the detailed structural data, may establish Xyl2 as a platform for future developments of industrially relevant xylanases.

CONCLUSIONS

F. oxysporum Xyl2 is a GH11 xylanase which is highly active in free form and immobilized onto a variety of solid supports in a wide pH range. Furthermore, immobilization of Xyl2 on certain supports significantly increases its thermal stability. A mechanistic rationale for Xyl2's remarkable catalytic efficiency at alkaline pH is proposed on the basis of two crystallographic structures. Together, these properties render Xyl2 an attractive biocatalyst for the sustainable industrial degradation of xylan.

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.