Bimodal substrate binding in the active site of the glycosidase BcX

Bacillus circulans xylanase (BcX) from the glycoside hydrolase family 11 degrades xylan through a retaining, double-displacement mechanism. The enzyme is thought to hydrolyze glycosidic bonds in a processive manner and has a large, active site cleft, with six subsites allowing the binding of six xylose units. Such an active site architecture suggests that oligomeric xylose substrates can bind in multiple ways. In the crystal structure of the catalytically inactive variant BcX E78Q, the substrate xylotriose is observed in the active site, as well as bound to the known secondary binding site and a third site on the protein surface. Nuclear magnetic resonance (NMR) titrations with xylose oligomers of different lengths yield nonlinear chemical shift trajectories for active site nuclei resonances, indicative of multiple binding orientations for these substrates for which binding and dissociation are in fast exchange on the NMR timescale, exchanging on the micro- to millisecond timescale. Active site binding can be modeled with a 2 : 1 model with dissociation constants in the low and high millimolar range. Extensive mutagenesis of active site residues indicates that tight binding occurs in the glycon binding site and is stabilized by Trp9 and the thumb region. Mutations F125A and W71A lead to large structural rearrangements. Binding at the glycon site is sensed throughout the active site, whereas the weak binding mostly affects the aglycon site. The interactions with the two active site locations are largely independent of each other and of binding at the secondary binding site.

Chemometric optimisation of enzymatic hydrolysis of beechwood xylan to target desired xylooligosaccharides

Generation of specific xylooligosaccharides (XOS) is attractive to the pharmaceutical and food industries due to the importance of their structure upon their application. This study used chemometrics to develop a comprehensive computational modelling set to predict the parameters maximising the generation of the desired XOS during enzymatic hydrolysis. The evaluated parameters included pH, temperature, substrate concentration, enzyme dosage and reaction time. A Box-Behnken design was combined with response surface methodology to develop the models. High-performance anion-exchange chromatography coupled with triple-quadrupole mass spectrometry (HPAEC-QqQ-MS) allowed the identification of 22 XOS within beechwood xylan hydrolysates. These data were used to validate the developed models and demonstrated their accuracy in predicting the parameters maximising the generation of the desired XOS. The maximum yields for X2-X6 were 314.2 ± 1.2, 76.6 ± 4.5, 38.4 ± 0.4, 17.8 ± 0.7, and 5.3 ± 0.2 mg/g xylan, respectively. These values map closely to the model predicted values 311.7, 92.6, 43.0, 16.3, and 4.9 mg/g xylan, respectively.

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.

Current Status of Carbohydrates Information in the Protein Data Bank

Carbohydrates are well known for their physicochemical, biological, functional, and therapeutic characteristics. Unfortunately, their chemical nature imposes severe challenges for the structural elucidation of these phenomena, impairing not only the depth of our understanding of carbohydrates but also the development of new biotechnological and therapeutic applications based on these molecules. In the recent past, the amount of structural information, obtained mainly from X-ray crystallography, has increased progressively, as well as its quality. In this context, the current work presents a global analysis of the carbohydrate information available in the Protein Data Bank (PDB). From high quality structures, it is clear that most of the data are highly concentrated on a few sets of residue types, on their monosaccharidic forms, and connected by a small diversity of glycosidic linkages. The geometries of these linkages can be mostly associated with the types of linkages instead of residues, while the level of puckering distortion was characterized, quantified, and located in a pseudorotational equilibrium landscape, not only to local minima but also to transitional states. These qualitative and quantitative analyses offer a global picture of the carbohydrate structural content in the PDB, potentially supporting the building of new models for carbohydrate-related biological phenomena at the atomistic level, including new developments on force field parameters.

Rationalising pKa shifts in Bacillus circulans xylanase with computational studies

Bacillus circulans xylanase (BcX), a family 11 glycoside hydrolase, catalyses the hydrolysis of xylose polymers with a net retention of stereochemistry. Glu172 in BcX is believed to act as a general acid by protonating the aglycone during glycosylation, and then as a general base to facilitate the deglycosylation step. The key to the dual role of this general acid/base lies in its protonation states, which depend on its intrinsic pK_a value and the specific environment which it resides within. To fully understand the detailed molecular features in BcX to establish the dual role of Glu172, we present a combined study based on both atomistic simulations and empirical models to calculate pK_a shifts for the general acid/base Glu172 in BcX at different functional states. Its pK_a values and those of nearby residues, obtained based on QM/MM free energy calculations, MCCE and PROPKA, show a good agreement with available experimental data. Additionally, our study provides additional insights into the effects of structural and electrostatic perturbations caused by mutations and chemical modifications, suggesting that the local solvation environment and mutagenesis of the residues adjacent to Glu172 establish its dual role during hydrolysis. The strengths and limitations of various methods for calculating pK_as and pK_a shifts have also been discussed.

Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains

This article describes a revised version 56A6(CARBO_R) of the GROMOS 56A6(CARBO) force field for hexopyranose-based carbohydrates. The simulated properties of unfunctionalized hexopyranoses are unaltered with respect to 56A6CARBO . In the context of both O1 -alkylated hexopyranoses and oligosaccharides, the revision stabilizes the regular (4) C1 chair for α-anomers, with the opposite effect for β-anomers. As a result, spurious ring inversions observed in α(1→4)-linked chains when using the original 56A6(CARBO) force field are alleviated. The (4) C1 chair is now the most stable conformation for all d-hexopyranose residues, irrespective of the linkage type and anomery, and of the position of the residue along the chain. The methylation of a d-hexopyranose leads to a systematic shift in the ring-inversion free energy ((4) C1 to (1) C4 ) by 7-8 kJ mol(-1), positive for the α-anomers and negative for the β-anomers, which is qualitatively compatible with the expected enhancement of the anomeric effect upon methylation at O1. The ring-inversion free energies for residues within chains are typically smaller in magnitude compared to those of the monomers, and correlate rather poorly with the latter. This suggests that the crowding of ring substituents upon chain formation alters the ring flexibility in a nonsystematic fashion. In general, the description of carbohydrate chains afforded by 56A6(CARBO_R) suggests a significant extent of ring flexibility, i.e., small but often non-negligible equilibrium populations of inverted chairs, and challenges the "textbook" picture of conformationally locked carbohydrate rings.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

The complete conformational free energy landscape of β-xylose reveals a two-fold catalytic itinerary for β-xylanases

Unraveling the conformational catalytic itinerary of glycoside hydrolases (GHs) is a growing topic of interest in glycobiology, with major impact in the design of GH inhibitors. β-xylanases are responsible for the hydrolysis of glycosidic bonds in β-xylans, a group of hemicelluloses of high biotechnological interest that are found in plant cell walls. The precise conformations followed by the substrate during catalysis in β-xylanases have not been unambiguously resolved, with three different pathways being proposed from structural analyses. In this work, we compute the conformational free energy landscape (FEL) of β-xylose to predict the most likely catalytic itineraries followed by β-xylanases. The calculations are performed by means of ab initio metadynamics, using the Cremer-Pople puckering coordinates as collective variables. The computed FEL supports only two of the previously proposed itineraries, 2SO → [2,5B]ǂ → 5S1 and 1S3 → [4H3]ǂ → 4C1, which clearly appear in low energy regions of the FEL. Consistently, 2SO and 1S3 are conformations preactivated for catalysis in terms of free energy/anomeric charge and bond distances. The results however exclude the OE → [OS2]ǂ → B2,5 itinerary that has been recently proposed for a family 11 xylanase. Classical and ab initio QM/MM molecular dynamics simulations reveal that, in this case, the observed OE conformation has been enforced by enzyme mutation. These results add a word of caution on using modified enzymes to inform on catalytic conformational itineraries of glycoside hydrolases.

Dynamics of the thumb-finger regions in a GH11 xylanase Bacillus circulans: comparison between the Michaelis and covalent intermediate

The structure and dynamics of B. circulans β-1,4-xylanase (BCX) were comparatively studied utilizing molecular dynamics.

Exploration of conformational flexibility and hydrogen bonding of xylosides in different solvents, as a model system for enzyme active site interactions

The predominantly populated conformation of carbohydrates in solution does not necessarily represent the biologically active species; rather, any conformer accessible without too large an energy penalty may be present in a biological pathway. Thus, the conformational preferences of a naphthyl xyloside, which initiates in vivo synthesis of antiproliferative glycosaminoglycans, have been studied by using NMR spectroscopy in a variety of solvents. Equilibria comprising the conformations (4)C1, (2)SO and (1)C4 were found, with a strong dependence on the hydrogen bonding ability of the solvent. Studies of fluorinated analogues revealed a direct hydrogen bond from the hydroxyl group at C2 to the fluorine atom at C4 by a (1h)JF4,HO2 coupling. Hydrogen bond directionality was further established via comparisons of fluorinated levoglucosan molecules.

The dynamics of the conformational changes in the hexopyranose ring: a transition path sampling approach

The transition paths corresponding to the conformational rearrangements in the ring of hexapyranose (α-d- and β-d-glucose) molecules were described by applying the transition path sampling method.

Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

Molecular approaches for ameliorating microbial xylanases

In industrial processes, chemical catalysis is being replaced by enzyme catalysis, since the latter is environmentally benign, non-persistent and cost effective. Microbial xylanases have significant applications in textile, baking, food and feed industries, and in paper and pulp industries for reducing the chlorine requirement. The hazardous chlorine required for bleaching can be reduced up to 25-30% by including an enzymatic step in the pulp bleaching process. The paper pulp bleaching requires xylanases that are active at alkaline pH and elevated temperatures. The enzymes from the cultured microbes do not perform optimally in the paper industry due to their inadequate stability under the process conditions of high temperature and alkaline pH. This review, therefore, deals with the rationale of molecular approaches such as protein engineering for designing xylanases with improved characteristics to suit the process conditions in industries, and prospects and problems.

Conformational analyses of the reaction coordinate of glycosidases

The enzymatic hydrolysis of the glycosidic bond is catalyzed by diverse enzymes generically termed glycoside hydrolases (hereafter GHs) or glycosidases. The many sequence-based families of glycosidases have served as a rich hunting ground for enzymologists for years. Not only are these enzymes of fundamental interest, providing paradigms for enzymatic catalysis that extend beyond the bounds of carbohydrate chemistry, but the enzymes themselves play myriad essential roles in diverse biological processes. The wide utility of glycosidases, from their industrial harnessing in the hydrolysis of plant biomass to their roles in human physiology and disease, has engendered a large scientific constituency with an interest in glycosidase chemistry. A fascinating thread of this research, and one with major impact on the design of enzyme inhibitors, is the conformational analysis of reaction pathways within the diverse families. These GH families provide a large pallet of enzymes with which chemists have attempted to depict the conformational landscape of glycosidase action. In this Account, we review three-dimensional insight into the conformational changes directed by glycosidases, primarily from structural observations of the stable enzyme-ligand species adjacent to the transition state (or states) and of enzyme-inhibitor complexes. We further show how recent computational advances dovetail with structural insight to provide a quantum mechanical basis for glycosidase action. The glycosidase-mediated hydrolysis of the acetal or ketal bond in a glycoside may occur with either inversion or retention of the configuration of the anomeric carbon. Inversion involves a single step and transition state, whereas retention, often referred to as the double displacement, is a two-step process with two transition states. The single transition state for the inverting enzymes and the two transition states (those flanking the covalent intermediate) in the double displacement have been shown to have substantial oxocarbenium ion character. The dissociative nature of these transition states results in significant relative positive charge accumulation on the pyranose ring. The delocalization of lone-pair electrons from the ring oxygen that stabilizes the cationic transition state implies that at, or close to, the transition states the pyranose will be distorted away from its lowest energy conformation to one that favors orbital overlap. Over the preceding decade, research has highlighted the harnessing of noncovalent interactions to aid this distortion of the sugar substrates from their lowest energy chair conformation to a variety of different boat, skew boat, and half-chair forms, each of which favors catalysis with a given enzyme and substrate. Crystallographic observation of stable species that flank the transition state (or states), of both retaining and inverting glycosidases, has allowed a description of their conformational itineraries, illustrating how enzymes facilitate the "electrophilic migration" of the anomeric center along the reaction coordinate. The blossoming of computational approaches, such as ab initio metadynamics, has underscored the quantum mechanical basis for glycoside hydrolysis. Conformational analyses highlight not only the itineraries used by enzymes, enabling their inhibition, but are also reflected in the nonenzymatic synthesis of glycosides, wherein chemists mimic strategies found in nature.

Three-dimensional structure of a thermophilic family GH11 xylanase from Thermobifida fusca

Thermostable enzymes employ various structural features dictated at the amino-acid sequence level that allow them to maintain their integrity at higher temperatures. Many hypotheses as to the nature of thermal stability have been proposed, including optimized core hydrophobicity and an increase in charged surface residues to enhance polar solvent interactions for solubility. Here, the three-dimensional structure of the family GH11 xylanase from the moderate thermophile Thermobifida fusca in its trapped covalent glycosyl-enzyme intermediate complex is presented. Interactions with the bound ligand show fewer direct hydrogen bonds from ligand to protein than observed in previous complexes from other species and imply that binding of the xylan substrate involves several water-mediated hydrogen bonds.

A reoptimized GROMOS force field for hexopyranose-based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers

This article presents a reoptimization of the GROMOS 53A6 force field for hexopyranose-based carbohydrates (nearly equivalent to 45A4 for pure carbohydrate systems) into a new version 56A(CARBO) (nearly equivalent to 53A6 for non-carbohydrate systems). This reoptimization was found necessary to repair a number of shortcomings of the 53A6 (45A4) parameter set and to extend the scope of the force field to properties that had not been included previously into the parameterization procedure. The new 56A(CARBO) force field is characterized by: (i) the formulation of systematic build-up rules for the automatic generation of force-field topologies over a large class of compounds including (but not restricted to) unfunctionalized polyhexopyranoses with arbritrary connectivities; (ii) the systematic use of enhanced sampling methods for inclusion of experimental thermodynamic data concerning slow or unphysical processes into the parameterization procedure; and (iii) an extensive validation against available experimental data in solution and, to a limited extent, theoretical (quantum-mechanical) data in the gas phase. At present, the 56A(CARBO) force field is restricted to compounds of the elements C, O, and H presenting single bonds only, no oxygen functions other than alcohol, ether, hemiacetal, or acetal, and no cyclic segments other than six-membered rings (separated by at least one intermediate atom). After calibration, this force field is shown to reproduce well the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers. As a result, the 56A(CARBO) force field should be suitable for: (i) the characterization of the dynamics of pyranose ring conformational transitions (in simulations on the microsecond timescale); (ii) the investigation of systems where alternative ring conformations become significantly populated; (iii) the investigation of anomerization or epimerization in terms of free-energy differences; and (iv) the design of simulation approaches accelerating the anomerization process along an unphysical pathway.

Catalysis: transition-state molecular recognition?

THE KEY TO UNDERSTANDING THE FUNDAMENTAL PROCESSES OF CATALYSIS IS THE TRANSITION STATE (TS): indeed, catalysis is a transition-state molecular recognition event. Practical objectives, such as the design of TS analogues as potential drugs, or the design of synthetic catalysts (including catalytic antibodies), require prior knowledge of the TS structure to be mimicked. Examples, both old and new, of computational modelling studies are discussed, which illustrate this fundamental concept. It is shown that reactant binding is intrinsically inhibitory, and that attempts to design catalysts that focus simply upon attractive interactions in a binding site may fail. Free-energy changes along the reaction coordinate for S(N)2 methyl transfer catalysed by the enzyme catechol-O-methyl transferase are described and compared with those for a model reaction in water, as computed by hybrid quantum-mechanical/molecular-mechanical molecular dynamics simulations. The case is discussed of molecular recognition in a xylanase enzyme that stabilises its sugar substrate in a (normally unfavourable) boat conformation and in which a single-atom mutation affects the free-energy of activation dramatically.

Analysis of the reaction coordinate of alpha-L-fucosidases: a combined structural and quantum mechanical approach

The enzymatic hydrolysis of alpha-L-fucosides is of importance in cancer, bacterial infections, and fucosidosis, a neurodegenerative lysosomal storage disorder. Here we show a series of snapshots along the reaction coordinate of a glycoside hydrolase family GH29 alpha-L-fucosidase unveiling a Michaelis (ES) complex in a (1)C(4) (chair) conformation and a covalent glycosyl-enzyme intermediate in (3)S(1) (skew-boat). First principles metadynamics simulations on isolated alpha-L-fucose strongly support a (1)C(4)<-->(3)H(4)<-->(3)S(1) conformational itinerary for the glycosylation step of the reaction mechanism and indicate a strong "preactivation" of the (1)C(4) complex to nucleophilic attack as reflected by free energy, C1-O1/O5-C1 bond length elongation/reduction, C1-O1 bond orientation, and positive charge development around the anomeric carbon. Analysis of an imino sugar inhibitor is consistent with tight binding of a chair-conformed charged species.

Characterization of temperature dependent and substrate-binding cleft movements in Bacillus circulans family 11 xylanase: a molecular dynamics investigation

BACKGROUND

Xylanases (EC 3.2.1.8) hydrolyze xylan, one of the most abundant plant polysaccharides found in nature, and have many potential applications in biotechnology.

METHODS

Molecular dynamics simulations were used to investigate the effects of temperature between 298 to 338 K and xylobiose binding on residues located in the substrate-binding cleft of the family 11 xylanase from Bacillus circulans (BcX).

RESULTS

In the absence of xylobiose the BcX exhibits temperature dependent movement of the thumb region which adopts an open conformation exposing the active site at the optimum catalytic temperature (328 K). In the presence of substrate, the thumb region restricts access to the active site at all temperatures, and this conformation is maintained by substrate/protein hydrogen bonds involving active site residues, including hydrogen bonds between Tyr69 and the 2' hydroxyl group of the substrate. Substrate access to the active site is regulated by temperature dependent motions that are restricted to the thumb region, and the BcX/substrate complex is stabilized by extensive intermolecular hydrogen bonding with residues in the active site.

GENERAL SIGNIFICANCE

These results call for a revision of both the "hinge-bending" model for the activity of group 11 xylanases, and the role of Tyr69 in the catalytic mechanism.