Crystallographic observation of a covalent catalytic intermediate in a beta-glycosidase

Electrostatic fingerprints of catalytically active amino acids in enzymes

The computed electrostatic and proton transfer properties are studied for 20 enzymes that represent all six major enzyme commission classes and a variety of different folds. The properties of aspartate, glutamate, and lysine residues that have been previously experimentally determined to be catalytically active are reported. The catalytic aspartate and glutamate residues studied here are strongly coupled to at least one other aspartate or glutamate residue and often to multiple other carboxylate residues with intrinsic pKa differences less than 1 pH unit. Sometimes these catalytic acidic residues are also coupled to a histidine residue, such that the intrinsic pKa of the acidic residue is higher than that of the histidine. All catalytic lysine residues studied here are strongly coupled to tyrosine or cysteine residues, wherein the intrinsic pKa of the anion-forming residue is higher than that of the lysine. Some catalytic lysines are also coupled to other lysines with intrinsic pKa differences within 1 pH unit. Some evidence of the possible types of interactions that facilitate nucleophilicity is discussed. The interactions reported here provide important clues about how side chain functional groups that are weak Brønsted acids or bases for the free amino acid in solution can achieve catalytic potency and become strong acids, bases or nucleophiles in the enzymatic environment.

Tuning the Transglycosylation Reaction of a GH11 Xylanase by a Delicate Enhancement of its Thumb Flexibility

Glycoside hydrolases (GHs) are attractive tools for multiple biotechnological applications. In conjunction with their hydrolytic function, GHs can perform transglycosylation under specific conditions. In nature, oligosaccharide synthesis is performed by glycosyltransferases (GTs); however, the industrial use of GTs is limited by their instability in solution. A key difference between GTs and GHs is the flexibility of their binding site architecture. We have used the xylanase from Bacillus circulans (BCX) to study the interplay between active-site flexibility and transglycosylation. Residues of the BCX "thumb" were substituted to increase the flexibility of the enzyme binding site. Replacement of the highly conserved residue P116 with glycine shifted the balance of the BCX enzymatic reaction toward transglycosylation. The effects of this point mutation on the structure and dynamics of BCX were investigated by NMR spectroscopy. The P116G mutation induces subtle changes in the configuration of the thumb and enhances the millisecond dynamics of the active site. Based on our findings, we propose the remodelling of the GH enzymes glycon site flexibility as a strategy to improve the transglycosylation efficiency of these biotechnologically important catalysts.

QM/MM Study of the Catalytic Reaction of Myrosinase; Importance of Assigning Proper Protonation States of Active-Site Residues

Myrosinase from Sinapis alba hydrolyzes glycosidic bonds of β-d-S-glucosides. The enzyme shows an enhanced activity in the presence of l-ascorbic acid. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics simulations to study the catalytic reaction of wild-type myrosinase and its E464A, Q187A, and Q187E mutants. Test calculations show that a proper QM region to study the myrosinase reaction must contain the whole substrate, models of Gln-187, Glu-409, Gln-39, His-141, Asn-186, Tyr-330, Glu-464, Arg-259, and a water molecule. Furthermore, to make the deglycosylation step possible, Arg-259 must be charged, Glu-464 must be protonated on OE2, and His-141 must be protonated on the NE2 atom. The results indicate that assigning proper protonation states of the residues is more important than the size of the model QM system. Our model reproduces the anomeric retaining characteristic of myrosinase and also reproduces the experimental fact that ascorbate increases the rate of the reaction. A water molecule in the active site, positioned by Gln-187, helps the aglycon moiety of the substrate to stabilize the buildup of negative charge during the glycosylation reaction and this in turn makes the moiety a better leaving group. The water molecule also lowers the glycosylation barrier by ∼9 kcal/mol. The results indicate that the Q187E and E464A mutants but not the Q187A mutant can perform the glycosylation step. However, the energy profiles for the deglycosylation step of the mutants are not similar to that of the wild-type enzyme. The Glu-464 residue lowers the barriers of the glycosylation and deglycosylation steps. The ascorbate ion can act as a general base in the reaction of the wild-type enzyme only if the Glu-464 and His-141 residues are properly protonated.

Dynamics of Ligand Binding to a Rigid Glycosidase*

The single-domain GH11 glycosidase from Bacillus circulans (BCX) is involved in the degradation of hemicellulose, which is one of the most abundant renewable biomaterials in nature. We demonstrate that BCX in solution undergoes minimal structural changes during turnover. NMR spectroscopy results show that the rigid protein matrix provides a frame for fast substrate binding in multiple conformations, accompanied by slow conversion, which is attributed to an enzyme-induced substrate distortion. A model is proposed in which the rigid enzyme takes advantage of substrate flexibility to induce a conformation that facilitates the acyl formation step of the hydrolysis reaction.

Modeling catalytic reaction mechanisms in glycoside hydrolases

Modeling catalysis in carbohydrate-active enzymes is a daunting challenge because of the high flexibility and diversity of both enzymes and carbohydrates. Glycoside hydrolases (GHs) are an illustrative example, where conformational changes and subtle interactions have been shown to be critical for catalysis. GHs have pivotal roles in industry (e.g. biofuel or detergent production) and biomedicine (e.g. targets for cancer and diabetes), and thus, a huge effort is devoted to unveil their molecular mechanisms. Besides experimental techniques, computational methods have served to provide an in-depth understanding of GH mechanisms, capturing complex reaction coordinates and the conformational itineraries that substrates follow during the whole catalytic pathway, providing a framework that ultimately may assist the engineering of these enzymes and the design of new inhibitors.

New paths of cyanogenesis from enzymatic-promoted cleavage of β-cyanoglucosides are suggested by a mixed DFT/QTAIM approach

Cyanogenesis is an enzyme-promoted cleavage of β-cyanoglucosides; the release of hydrogen cyanide is believed to produce food poisoning by consumption of certain crops as Cassava (Manihot esculenta Crantz). The production of hydrogen cyanide by some disruption of the plant wall is related to the content of two β-cyanoglucosides (linamarin and lotaustralin) which are stored within the tuber. Some features about the mechanistic bases of these transformations have been published; nevertheless, there are still questions about the exact mechanism, such as the feasibility of a difference in the kinetics of cyanogenesis between both cyanoglucosides. In this work, we have performed a theoretical analysis using DFT and QTAIM theoretical frameworks to propose a feasible mechanism of the observed first step of the enzyme-catalyzed rupture of these glucosides; our results led us to explain the observed difference between linamarin and lotaustralin. Meanwhile, DFT studies suggest that there are no differences between local reactivity indexes of both glucosides; QTAIM topological analysis suggests two important intramolecular interactions which we found to fix the glucoside in such a way that suggests the linamarin as a more reactive system towards a nucleophilic attack, thus explaining the readiness to liberate hydrogen cyanide.

Enhanced Azidolysis by the Formation of Stable Ser-His Catalytic Dyad in a Glycoside Hydrolase Family 10 Xylanase Mutant

Glycoside hydrolases require carboxyl groups as catalysts for their activity. A retaining xylanase from Streptomyces olivaceoviridis E-86 belonging to glycoside hydrolase family 10 possesses Glu128 and Glu236 that respectively function as acid/base and nucleophile. We previously developed a unique mutant of the retaining xylanase, N127S/E128H, whose deglycosylation is triggered by azide. A crystallographic study reported that the transient formation of a Ser–His catalytic dyad in the reaction cycle possibly reduced the azidolysis reaction. In the present study, we engineered a catalytic dyad with enhanced stability by site-directed mutagenesis and crystallographic study of N127S/E128H. Comparison of the Michaelis complexes of N127S/E128H with pNP-X₂ and with xylopentaose showed that Ser127 could form an alternative hydrogen bond with Thr82, which disrupts the formation of the Ser–His catalytic dyad. The introduction of T82A mutation in N127S/E128H produces an enhanced first-order rate constant (6 times that of N127S/E128H). We confirmed the presence of a stable Ser–His hydrogen bond in the Michaelis complex of the triple mutant, which forms the productive tautomer of His128 that acts as an acid catalyst. Because the glycosyl azide is applicable in the bioconjugation of glycans by using click chemistry, the enzyme-assisted production of the glycosyl azide may contribute to the field of glycobiology.

Insight into the functional roles of Glu175 in the hyperthermostable xylanase XYL10C-ΔN through structural analysis and site-saturation mutagenesis

BACKGROUND

Improving the hydrolytic performance of hemicellulases to degrade lignocellulosic biomass is of considerable importance for second-generation biorefinery. Xylanase, as the crucial hemicellulase, must be thermostable and have high activity for its potential use in the bioethanol industry. To obtain excellent xylanase candidates, it is necessary to understand the structure-function relationships to provide a meaningful reference to improve the enzyme properties. This study aimed to investigate the catalytic mechanism of a highly active hyperthermophilic xylanase variant, XYL10C-ΔN, for hemicellulose degradation.

RESULTS

By removing the N-terminal 66 amino acids, the variant XYL10C-ΔN showed a 1.8-fold improvement in catalytic efficiency and could hydrolyze corn stover more efficiently in hydrolysis of corn stover; however, it retained similar thermostability to the wild-type XYL10C. Based on the crystal structures of XYL10C-ΔN and its complex with xylobiose, Glu175 located on loop 3 was found to be specific to GH10 xylanases and probably accounts for the excellent enzyme properties by interacting with Lys135 and Met137 on loop 2. Site-saturation mutagenesis confirmed that XYL10C-ΔN with glutamate acid at position 175 had the highest catalytic efficiency, specific activity, and the broadest pH-activity profile. The functional roles of Glu175 were also verified in the mutants of another two GH10 xylanases, XylE and XynE2, which showed increased catalytic efficiencies and wider pH-activity profiles.

CONCLUSIONS

XYL10C-ΔN, with excellent thermostability, high catalytic efficiency, and great lignocellulose-degrading capability, is a valuable candidate xylanase for the biofuel industry. The mechanism underlying improved activity of XYN10C-ΔN was thus investigated through structural analysis and functional verification, and Glu175 was identified to play the key role in the improved catalytic efficiency. This study revealed the importance of a key residue (Glu175) in XYN10C-ΔN and provides a reference to modify GH10 xylanases for improved catalytic performance.

Structural Considerations on the Use of Endo-Xylanases for the Production of prebiotic Xylooligosaccharides from Biomass

Xylooligosaccharides (XOS) have gained increased interest as prebiotics during the last years. XOS and arabinoxylooligosaccharides (AXOS) can be produced from major fractions of biomass including agricultural by-products and other low cost raw materials. Endo-xylanases are key enzymes for the production of (A)XOS from xylan. As the xylan structure is broadly diverse due to different substitutions, diverse endo-xylanases have evolved for its degradation. In this review structural and functional aspects are discussed, focusing on the potential applications of endo-xylanases in the production of differently substituted (A)XOS as emerging prebiotics, as well as their implication in the processing of the raw materials. Endo-xylanases are found in at least eight different glycoside hydrolase families (GH), and can either have a retaining or an inverting catalytic mechanism. To date, it is mainly retaining endo-xylanases that are used in applications to produce (A)XOS. Enzymes from these GH-families (mainly GH10 and GH11, and the more recently investigated GH30) are taken as prototypes to discuss substrate preferences and main products obtained. Finally, the need of new and accessory enzymes (new specificities from new families or sources) to increase the yield of different types of (A)XOS is discussed, along with in vitro tests of produced oligosaccharides and production of enzymes in GRAS organisms to facilitate use in functional food manufacturing.

Enzymatic Cleavage of Glycosidic Bonds: Strategies on How to Set Up and Control a QM/MM Metadynamics Simulation

Carbohydrates play crucial roles in many biological processes, from cell-cell adhesion to chemical signaling. Their complexity and diversity, related to α/β anomeric configuration, ring substituents, and conformational variations, require a diverse set of enzymes for their processing. Among them, glycoside hydrolases (GHs) are responsible for the hydrolysis of one of the strongest bonds in nature: the glycosidic bond. These highly specialized biological catalysts select particular conformations their carbohydrate substrates to enhance catalysis. The evolution of this conformation during the reaction of glycosidic bond cleavage, known as the conformational catalytic itinerary, is of fundamental interest in glycobiology, with impact on inhibitor and drug design. Here we review some of the aspects and the main strategies one needs to take into account when simulating a reaction in a GH enzyme using QM/MM metadynamics. Several specific aspects are highlighted, from the importance of the distortion of the substrate at the Michaelis complex to the variable control during the metadynamics simulation or the analysis of the reaction mechanism and conformational itinerary. The increasing speed of computer power and methodological advances have added a vital tool to the study of GH mechanisms, as shown here and recent reviews. It is hoped that this chapter will serve as a first guide for those attempting to perform a metadynamics simulation of these relevant and fascinating enzymes.

Mechanism-based inhibitors of glycosidases: design and applications

This article covers recent developments in the design and application of activity-based probes (ABPs) for glycosidases, with emphasis on the different enzymes involved in metabolism of glucosylceramide in humans. Described are the various catalytic reaction mechanisms employed by inverting and retaining glycosidases. An understanding of catalysis at the molecular level has stimulated the design of different types of ABPs for glycosidases. Such compounds range from (1) transition-state mimics tagged with reactive moieties, which associate with the target active site—forming covalent bonds in a relatively nonspecific manner in or near the catalytic pocket—to (2) enzyme substrates that exploit the catalytic mechanism of retaining glycosidase targets to release a highly reactive species within the active site of the enzyme, to (3) probes based on mechanism-based, covalent, and irreversible glycosidase inhibitors. Some applications in biochemical and biological research of the activity-based glycosidase probes are discussed, including specific quantitative visualization of active enzyme molecules in vitro and in vivo, and as strategies for unambiguously identifying catalytic residues in glycosidases in vitro.

Biochemical properties and atomic resolution structure of a proteolytically processed β-mannanase from cellulolytic Streptomyces sp. SirexAA-E

β-Mannanase SACTE_2347 from cellulolytic Streptomyces sp. SirexAA-E is abundantly secreted into the culture medium during growth on cellulosic materials. The enzyme is composed of domains from the glycoside hydrolase family 5 (GH5), fibronectin type-III (Fn3), and carbohydrate binding module family 2 (CBM2). After secretion, the enzyme is proteolyzed into three different, catalytically active variants with masses of 53, 42 and 34 kDa corresponding to the intact protein, loss of the CBM2 domain, or loss of both the Fn3 and CBM2 domains. The three variants had identical N-termini starting with Ala51, and the positions of specific proteolytic reactions in the linker sequences separating the three domains were identified. To conduct biochemical and structural characterizations, the natural proteolytic variants were reproduced by cloning and heterologously expressed in Escherichia coli. Each SACTE_2347 variant hydrolyzed only β-1,4 mannosidic linkages, and also reacted with pure mannans containing partial galactosyl- and/or glucosyl substitutions. Examination of the X-ray crystal structure of the GH5 domain of SACTE_2347 suggests that two loops adjacent to the active site channel, which have differences in position and length relative to other closely related mannanases, play a role in producing the observed substrate selectivity.

Structure, dynamics, and specificity of endoglucanase D from Clostridium cellulovorans

The enzymatic degradation of cellulose is a critical step in the biological conversion of plant biomass into an abundant renewable energy source. An understanding of the structural and dynamic features that cellulases utilize to bind a single strand of crystalline cellulose and hydrolyze the β-1,4-glycosidic bonds of cellulose to produce fermentable sugars would greatly facilitate the engineering of improved cellulases for the large-scale conversion of plant biomass. Endoglucanase D (EngD) from Clostridium cellulovorans is a modular enzyme comprising an N-terminal catalytic domain and a C-terminal carbohydrate-binding module, which is attached via a flexible linker. Here, we present the 2.1-Å-resolution crystal structures of full-length EngD with and without cellotriose bound, solution small-angle X-ray scattering (SAXS) studies of the full-length enzyme, the characterization of the active cleft glucose binding subsites, and substrate specificity of EngD on soluble and insoluble polymeric carbohydrates. SAXS data support a model in which the linker is flexible, allowing EngD to adopt an extended conformation in solution. The cellotriose-bound EngD structure revealed an extended active-site cleft that contains seven glucose-binding subsites, but unlike the majority of structurally determined endocellulases, the active-site cleft of EngD is partially enclosed by Trp162 and Tyr232. EngD variants, which lack Trp162, showed a significant reduction in activity and an alteration in the distribution of cellohexaose degradation products, suggesting that Trp162 plays a direct role in substrate binding.

Protein dielectric constants determined from NMR chemical shift perturbations

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

The crystallization and structural analysis of cellulases (and other glycoside hydrolases): strategies and tactics

The three-dimensional (3-D) structures of cellulases, and other glycoside hydrolases, are a central feature of research in carbohydrate chemistry and biochemistry. 3-D structure is used to inform protein engineering campaigns, both academic and industrial, which are typically used to improve the stability or activity of an enzyme. Examples of classical protein engineering goals include higher thermal stability, reduced metal-ion dependency, detergent and protease resistance, decreased product inhibition, and altered specificity. 3-D structure may also be used to interpret the behavior of enzyme variants that are derived from screening or random mutagenesis approaches, with a view to establishing an iterative design process. In other areas, 3-D structure is used as one of the many tools to probe enzymatic catalysis, typically dovetailing with physical organic chemistry approaches to provide complete reaction mechanisms for enzymes by visualizing catalytic site interactions at different stages of the reaction. Such mechanistic insight is not only fundamentally important, impacting on inhibitor and drug design approaches with ramifications way beyond cellulose hydrolysis, but also provides the framework for the design of enzyme variants to use as biocatalysts for the synthesis of bespoke oligosaccharides. Here we review some of the strategies and tactics that may be applied to the X-ray structure solution of cellulases (and other carbohydrate-active enzymes). The general approach is first to decide why you are doing the work, then to establish correct domain boundaries for truncated constructs (typically the catalytic domain only), and finally to pursue crystallization of pure, homogeneous, and monodisperse protein with appropriate ligand and additive combinations. Cellulase-specific strategies are important for the delineation of domain boundaries, while glycoside hydrolases generally also present challenges and opportunities for the selection and optimization of ligands to both aid crystallization, and also provide structural and mechanistic insight. As the many roles for plant cell wall degrading enzymes increase, so does the need for rapid high-quality structure determination to provide a sound structural foundation for understanding mechanism and specificity, and for future protein engineering strategies.

Cloning and molecular characterization of a new fungal xylanase gene from Sclerotinia sclerotiorum S2

Sclerotinia sclerotiorum fungus has three endoxylanases induced by wheat bran. In the first part, a partial xylanase sequence gene (90 bp) was isolated by PCR corresponding to catalytic domains (β 5 and β 6 strands of this protein). The high homology of this sequence with xylanase of Botryotinia fuckeliana has permitted in the second part to amplify the XYN1 gene. Sequence analysis of DNA and cDNA revealed an ORF of 746 bp interrupted by a 65 bp intron, thus encoding a predicted protein of 226 amino acids. The mature enzyme (20.06 kDa), is coded by 188 amino acid (pI 9.26). XYN1 belongs to G/11 glycosyl hydrolases family with a conserved catalytic domain containing E(86) and E(178) residues. Bioinformatics analysis revealed that there was no Asn-X-Ser/Thr motif required for N-linked glycosylation in the deduced sequence however, five O-glycosylation sites could intervene in the different folding of xylanses isoforms and in their secretary pathway.

Mechanistic insights into the 1,3-xylanases: useful enzymes for manipulation of algal biomass

Xylanases capable of degrading the crystalline microfibrils of 1,3-xylan that reinforce the cell walls of some red and siphonous green algae have not been well studied, yet they could prove to be of great utility in algaculture for the production of food and renewable chemical feedstocks. To gain a better mechanistic understanding of these enzymes, a suite of reagents was synthesized and evaluated as substrates and inhibitors of an endo-1,3-xylanase. With these reagents, a retaining mechanism was confirmed for the xylanase, its catalytic nucleophile identified, and the existence of -3 to +2 substrate-binding subsites demonstrated. Protein crystal X-ray diffraction methods provided a high resolution structure of a trapped covalent glycosyl-enzyme intermediate, indicating that the 1,3-xylanases likely utilize the (1)S(3) → (4)H(3) → (4)C(1) conformational itinerary to effect catalysis.

Solution and gas-phase H/D exchange of protein-small-molecule complexes: Cex and its inhibitors

The properties of noncovalent complexes of the enzyme exo-1,4-β-D-glycanase ("Cex") with three aza-sugar inhibitors, deoxynojirimycin (X(2)DNJ), isofagomine lactam (X(2)IL), and isofagomine (X(2)IF), have been studied with solution and gas-phase hydrogen deuterium exchange (H/Dx) and measurements of collision cross sections of gas-phase ions. In solution, complexes have lower H/Dx levels than free Cex because binding the inhibitors blocks some sites from H/Dx and reduces fluctuations of the protein. In mass spectra of complexes, abundant Cex ions are seen, which mostly are formed by dissociation of complexes in the ion sampling interface. Both complex ions and Cex ions formed from a solution containing complexes have lower cross sections than Cex ions from a solution of Cex alone. This suggests the Cex ions formed by dissociation "remember" their solution conformations. For a given charge, ions of the complexes have greater gas-phase H/Dx levels than ions of Cex. Unlike cross sections, H/Dx levels of the complexes do not correlate with the relative gas-phase binding strengths measured by MS/MS. Cex ions from solutions with or without inhibitors, which have different cross sections, show the same H/Dx level after 15 s, indicating the ions may fold or unfold on the seconds time scale of the H/Dx experiment. Thus, cross sections show that complexes have more compact conformations than free protein ions on the time scale of ca. 1 ms. The gas-phase H/Dx measurements show that at least some complexes retain different conformations from the Cex ions on a time scale of seconds.

Functional analysis of hyperthermophilic endocellulase from Pyrococcus horikoshii by crystallographic snapshots

A hyperthermophilic membrane-related β-1,4-endoglucanase (family 5, cellulase) of the archaeon Pyrococcus horikoshii was found to be capable of hydrolysing cellulose at high temperatures. The hyperthermophilic cellulase has promise for applications in biomass utilization. To clarify its detailed function, we determined the crystal structures of mutants of the enzyme in complex with either the substrate or product ligands. We were able to resolve different kinds of complex structures at 1.65–2.01 Å (1 Å=0.1 nm). The structural analysis of various mutant enzymes yielded a sequence of crystallographic snapshots, which could be used to explain the catalytic process of the enzyme. The substrate position is fixed by the alignment of one cellobiose unit between the two aromatic amino acid residues at subsites +1 and +2. During the enzyme reaction, the glucose structure of cellulose substrates is distorted at subsite −1, and the β-1,4-glucoside bond between glucose moieties is twisted between subsites −1 and +1. Subsite −2 specifically recognizes the glucose residue, but recognition by subsites +1 and +2 is loose during the enzyme reaction. This type of recognition is important for creation of the distorted boat form of the substrate at subsite −1. A rare enzyme–substrate complex was observed within the low-activity mutant Y299F, which suggested the existence of a trapped ligand structure before the formation by covalent bonding of the proposed intermediate structure. Analysis of the enzyme–substrate structure suggested that an incoming water molecule, essential for hydrolysis during the retention process, might be introduced to the cleavage position after the cellobiose product at subsites +1 and +2 was released from the active site.

Co-crystals of 3-deoxy-3-fluoro-α-D-glucopyranose and 3-deoxy-3-fluoro-β-D-glucopyranose

3-Deoxy-3-fluoro-D-glucopyranose crystallizes from acetone to give a unit cell containing two crystallographically independent molecules. One of these molecules (at site A) is structurally homogeneous and corresponds to 3-deoxy-3-fluoro-β-D-glucopyranose, C(6)H(11)FO(5), (I). The second molecule (at site B) is structurally heterogeneous and corresponds to a mixture of (I) and 3-deoxy-3-fluoro-α-D-glucopyranose, (II); treatment of the diffraction data using partial-occupancy oxygen at the anomeric center gave a high-quality packing model with an occupancy ratio of 0.84:0.16 for (II):(I) at site B. The mixture of α- and β-anomers at site B appears to be accommodated in the lattice because hydrogen-bonding partners are present to hydrogen bond to the anomeric OH group in either an axial or equatorial orientation. Cremer-Pople analysis of (I) and (II) shows the pyranosyl ring of (II) to be slightly more distorted than that of (I) [θ((I)) = 3.85 (15)° and θ((II)) = 6.35 (16)°], but the general direction of distortion is similar in both structures [ϕ((I)) = 67 (2)° (B(C1,C4)) and ϕ((II)) = 26.0 (15)° ((C3)TB(C1)); B = boat conformation and TB = twist-boat conformation]. The exocyclic hydroxymethyl (-CH(2)OH) conformation is gg (gauche-gauche) (H5 anti to O6) in both (I) and (II). Structural comparisons of (I) and (II) to related unsubstituted, deoxy and fluorine-substituted monosaccharides show that the gluco ring can assume a wide range of distorted chair structures in the crystalline state depending on ring substitution patterns.

4-Deoxy-4-fluoro-β-D-glucopyranose

4-Deoxy-4-fluoro-β-D-glucopyranose, C(6)H(11)FO(5), (I), crystallizes from water at room temperature in a slightly distorted (4)C(1) chair conformation. The observed chair distortion differs from that observed in β-D-glucopyranose [Kouwijzer, van Eijck, Kooijman & Kroon (1995). Acta Cryst. B51, 209-220], (II), with the former skewed toward a B(C3,O5) (boat) conformer and the latter toward an (O5)TB(C2) (twist-boat) conformer, based on Cremer-Pople analysis. The exocyclic hydroxymethyl group conformations in (I) and (II) are similar; in both cases, the O-C-C-O torsion angle is ∼-60° (gg conformer). Intermolecular hydrogen bonding in the crystal structures of (I) and (II) is conserved in that identical patterns of donors and acceptors are observed for the exocyclic substituents and the ring O atom of each monosaccharide. Inspection of the crystal packing structures of (I) and (II) reveals an essentially identical packing configuration.

Crystallographic snapshots of an entire reaction cycle for a retaining xylanase from Streptomyces olivaceoviridis E-86

Retaining glycosyl hydrolases, which catalyse both glycosylation and deglycosylation in a concerted manner, are the most abundant hydrolases. To date, their visualization has tended to be focused on glycosylation because glycosylation reactions can be visualized by inactivating deglycosylation step and/or using substrate analogues to isolate covalent intermediates. Furthermore, during structural analyses of glycosyl hydrolases with hydrolytic reaction products by the conventional soaking method, mutarotation of an anomeric carbon in the reaction products promptly and certainly occurs. This undesirable structural alteration hinders visualization of the second step in the reaction. Here, we investigated X-ray crystallographic visualization as a possible method for visualizing the conformational itinerary of a retaining xylanase from Streptomyces olivaceoviridis E-86. To clearly define the stereochemistry at the anomeric carbon during the deglycosylation step, extraneous nucleophiles, such as azide, were adopted to substitute for the missing base catalyst in an appropriate mutant. The X-ray crystallographic visualization provided snapshots of the components of the entire reaction, including the E*S complex, the covalent intermediate, breakdown of the intermediate and the enzyme-product (E*P)complex.

Use of Fluorinated Functionality in Enzyme Inhibitor Development: Mechanistic and Analytical Advantages

On the one hand, owing to its electronegativity, relatively small size, and notable leaving group ability from anionic intermediates, fluorine offers unique opportunities for mechanism-based enzyme inhibitor design. On the other, the "bio-orthogonal" and NMR-active 19-fluorine nucleus allows the bioorganic chemist to follow the mechanistic fate of fluorinated substrate analogues or inhibitors as they are enzymatically processed. This article takes an overview of the field, highlighting key developments along these lines. It begins by highlighting new screening methodologies for drug discovery that involve appropriate tagging of either substrate or the target protein itself with (19)F-markers, that then report back on turnover and binding, respectively, via an the NMR screen. Taking this one step further, substrate-tagging with fluorine can be done is such a manner as to provide stereochemical information on enzyme mechanism. For example, substitution of one of the terminal hydrogens in phosphoenolpyruvate, provides insight into the, otherwise latent, facial selectivity of C-C bond formation in KDO synthase. Perhaps, most importantly, from the point of view of this discussion, appropriately tailored fluorinated functionality can be used to form to stabilized "transition state analogue" complexes with a target enzymes. Thus, 5-fluorinated pyrimidines, alpha-fluorinated ketones, and 2-fluoro-2-deoxysugars each lead to covalent adduction of catalytic active site residues in thymidylate synthase, serine protease and glycosidase enzymes, respectively. In all such cases, (19)F NMR allows the bioorganic chemist to spectrally follow "transition state analogue" formation. Finally, the use of specific fluorinated functionality to engineer "suicide substrates" is highlighted in a discussion of the development of the alpha-(2'Z-fluoro)vinyl trigger for amino acid decarboxylase inactivation. Here (19)F NMR allows the bioorganic chemist to glean useful partition ratio data directly out of the NMR tube.

Engineering the thermostability of a TIM-barrel enzyme by rational family shuffling

A possible approach to generate enzymes with an engineered temperature optimum is to create chimeras of homologous enzymes with different temperature optima. We tested this approach using two family-10 xylanases from Thermotoga maritima: the thermophilic xylanase A catalytic domain (TmxAcat, T(opt)=68 degrees C), and the hyperthermophilic xylanase B (TmxB, T(opt)=102 degrees C). Twenty-one different chimeric constructs were created by mimicking family shuffling in a rational manner. The measured temperature optima of the 16 enzymatically active chimeras do not monotonically increase with the percentage of residues coming from TmxB. Only four chimeras had a higher temperature optimum than TmxAcat, the most stable variant (T(opt)=80 degrees C) being the one in which both terminal segments came from TmxB. Further analysis suggests that the interaction between the N- and C-terminal segments has a disproportionately high contribution to the overall thermostability. The results may be generalizable to other enzymes where the N- and C-termini are in contact.

Synthesis and testing of 2-deoxy-2,2-dihaloglycosides as mechanism-based inhibitors of alpha-glycosidases

The synthesis of a series of 2-deoxy-2,2-dihaloglycosyl halides as potential alpha-glycosidase inactivators has been achieved via the halogenation of protected 2-fluoroglycal precursors. Direct chlorination of per-O-acetylated 2-fluoro-d-glucal and 2-fluoromaltal followed by basic deprotection yielded the corresponding 2-chloro-2-deoxy-2-fluoroglycosyl chlorides. Reaction of the per-O-acetylated 2-fluoroglycals with acetyl hypofluorite or Selectfluor yielded the 2-deoxy-2,2-difluoroglycosyl derivatives, which were converted to their alpha-chlorides using thionyl chloride and deprotected under basic conditions. Trinitrophenyl glycosides of the 2-deoxy-2,2-difluoro mono- and disaccharides were synthesized by arylation of the hemiacetals with picryl fluoride, then deprotected with HCl in methanol. All three monosaccharide derivatives caused active site-directed, time-dependent inactivation of yeast alpha-glucosidase via the trapping of covalent glycosyl-enzyme intermediates, and kinetic parameters for inactivation by each compound were determined. Surprisingly neither of the 2-deoxy-2,2-dihalomaltosyl chlorides caused time-dependent inactivation of human pancreatic alpha-amylase, despite the fact that the trinitrophenyl 2-deoxy-2,2-difluoromaltoside functioned in that mode. The trinitrophenyl glycosides appear to be approximately 1000-fold more reactive than the corresponding chlorides in the enzyme active sites.

Predicting and annotating catalytic residues: an information theoretic approach

We introduce a computational method to predict and annotate the catalytic residues of a protein using only its sequence information, so that we describe both the residues' sequence locations (prediction) and their specific biochemical roles in the catalyzed reaction (annotation). While knowing the chemistry of an enzyme's catalytic residues is essential to understanding its function, the challenges of prediction and annotation have remained difficult, especially when only the enzyme's sequence and no homologous structures are available. Our sequence-based approach follows the guiding principle that catalytic residues performing the same biochemical function should have similar chemical environments; it detects specific conservation patterns near in sequence to known catalytic residues and accordingly constrains what combination of amino acids can be present near a predicted catalytic residue. We associate with each catalytic residue a short sequence profile and define a Kullback-Leibler (KL) distance measure between these profiles, which, as we show, effectively captures even subtle biochemical variations. We apply the method to the class of glycohydrolase enzymes. This class includes proteins from 96 families with very different sequences and folds, many of which perform important functions. In a cross-validation test, our approach correctly predicts the location of the enzymes' catalytic residues with a sensitivity of 80% at a specificity of 99.4%, and in a separate cross-validation we also correctly annotate the biochemical role of 80% of the catalytic residues. Our results compare favorably to existing methods. Moreover, our method is more broadly applicable because it relies on sequence and not structure information; it may, furthermore, be used in conjunction with structure-based methods.

Peptide insertion, positioning, and stabilization in a membrane: insight from an all-atom molecular dynamics simulation

Peptide insertion, positioning, and stabilization in a model membrane are probed via an all-atom molecular dynamics (MD) simulation. One peptide (WL5) is simulated in each leaflet of a solvated dimyristoylglycero-3-phosphate (DMPC) membrane. Within the first 5 ns, the peptides spontaneously insert into the membrane and then stabilize during the remaining 70 ns of simulation time. In both leaflets, the peptides localize to the membrane interface, and this localization is attributed to the formation of peptide-lipid hydrogen bonds. We show that the single tryptophan residue in each peptide contributes significantly to these hydrogen bonds; specifically, the nitrogen heteroatom of the indole ring plays a critical role. The tilt angles of the indole rings relative to the membrane normal in the upper and lower leaflets are approximately 26 degrees and 54 degrees , respectively. The tilt angles of the entire peptide chain are 62 degrees and 74 degrees . The membrane induces conformations of the peptide that are characteristic of beta-sheets, and the peptide enhances the lipid ordering in the membrane. Finally, the diffusion rate of the peptides in the membrane plane is calculated (based on experimental peptide concentrations) to be approximately 6 A(2)/ns, thus suggesting a 500 ns time scale for intermolecular interactions.

The Crystal Structure of Human Cytosolic β-Glucosidase Unravels the Substrate Aglycone Specificity of a Family 1 Glycoside Hydrolase

Human cytosolic beta-glucosidase (hCBG) is a xenobiotic-metabolizing enzyme that hydrolyses certain flavonoid glucosides, with specificity depending on the aglycone moiety, the type of sugar and the linkage between them. In this study, the substrate preference of this enzyme was investigated by mutational analysis, X-ray crystallography and homology modelling. The crystal structure of hCBG was solved by the molecular replacement method and refined at 2.7 A resolution. The main-chain fold of the enzyme belongs to the (beta/alpha)(8) barrel structure, which is common to family 1 glycoside hydrolases. The active site is located at the bottom of a pocket (about 16 A deep) formed by large surface loops, surrounding the C termini of the barrel of beta-strands. As for all the clan of GH-A enzymes, the two catalytic glutamate residues are located on strand 4 (the acid/base Glu165) and on strand 7 (the nucleophile Glu373). Although many features of hCBG were shown to be very similar to previously described enzymes from this family, crucial differences were observed in the surface loops surrounding the aglycone binding site, and these are likely to strongly influence the substrate specificity. The positioning of a substrate molecule (quercetin-4'-glucoside) by homology modelling revealed that hydrophobic interactions dominate the binding of the aglycone moiety. In particular, Val168, Trp345, Phe225, Phe179, Phe334 and Phe433 were identified as likely to be important in determining substrate specificity in hCBG, and site-directed mutagenesis supported a key role for some of these residues.

Crystal structures of Paenibacillus polymyxa beta-glucosidase B complexes reveal the molecular basis of substrate specificity and give new insights into the catalytic machinery of family I glycosidases

Bacteria species involved in degradation of cellulosic substrates produce a variety of enzymes for processing related compounds along the hydrolytic pathway. Paenibacillus polymyxa encodes two homologous beta-glucosidases, BglA and BglB, presenting different quaternary structures and substrate specificities. We previously reported the 3D-structure of BglA, which is highly specific against cellobiose. Here, we present structural analysis of BglB, a monomeric enzyme that acts as an exo-beta-glucosidase hydrolyzing cellobiose and cellodextrins of higher degree of polymerization. The crystal structure of BglB shows that several polar residues narrow the active site pocket and contour additional subsites. The structure of the BglB-cellotetraose complex confirms these subsites, revealing the substrate-binding mode, and shows the oligosaccharide-enzyme recognition pattern in detail. Comparison between BglA and BglB crystal structures suggests that oligomerization in BglA can assist in fine-tuning the specificity of the active centre by modulating the loops surrounding the cavity. We have solved the crystal structure of BglB with bound thiocellobiose, a competitive inhibitor, which together with the BglB-cellotetraose complex delineate the general features of the aglycon site. The detailed characterization of the atomic interactions at the aglycon site show a recognition pattern common to all bacterial beta-glucosidases, and presents some differences with the aglycon site in plant beta-glycosidases essentially by means of a different orientation of the basal Trp. The crystal structures of of BglB with a covalently bound inhibitor (derived from 2-fluoroglucoside) and glucose (produced by hydrolysis of the substrate in the crystal), provide additional pictures of the binding events and the intermediates formed during the reaction. Altogether, this information can assist in the understanding of subtle differences of the enzyme mechanism and substrate recognition within this family of enzymes, and consequently it can help in the development of new enzymes with improved activity or specificity.