The Conformational Free-Energy Landscape of β-d-Mannopyranose: Evidence for a 1S5 → B2,5 → OS2 Catalytic Itinerary in β-Mannosidases

Accelerated Lignocellulosic Molecule Adsorption Structure Determination

Here, we present a study combining Bayesian optimization structural inference with the machine learning interatomic potential Neural Equivariant Interatomic Potential (NequIP) to accelerate and enable the study of the adsorption of the conformationally flexible lignocellulosic molecules β-d-xylose and 1,4-β-d-xylotetraose on a copper surface. The number of structure evaluations needed to map out the relevant potential energy surfaces are reduced by Bayesian optimization, while NequIP minimizes the time spent on each evaluation, ultimately resulting in cost-efficient and reliable sampling of large systems and configurational spaces. Although the applicability of Bayesian optimization for the conformational analysis of the more flexible xylotetraose molecule is restricted by the sample complexity bottleneck, the latter can be effectively bypassed with external conformer search tools, such as the Conformer-Rotamer Ensemble Sampling Tool, facilitating the subsequent lower-dimensional global minimum adsorption structure determination. Finally, we demonstrate the applicability of the described approach to find adsorption structures practically equivalent to the density functional theory counterparts at a fraction of the computational cost.

Control of Substrate Conformation by Hydrogen Bonding in a Retaining β-Endoglycosidase

Bacterial β-glycosidases are hydrolytic enzymes that depolymerize polysaccharides such as β-cellulose, β-glucans and β-xylans from different sources, offering diverse biomedical and industrial uses. It has been shown that a conformational change of the substrate, from a relaxed 4 C1 conformation to a distorted 1 S3 /1,4 B conformation of the reactive sugar, is necessary for catalysis. However, the molecular determinants that stabilize the substrate's distortion are poorly understood. Here we use quantum mechanics/molecular mechanics (QM/MM)-based molecular dynamics methods to assess the impact of the interaction between the reactive sugar, i. e. the one at subsite -1, and the catalytic nucleophile (a glutamate) on substrate conformation. We show that the hydrogen bond involving the C2 exocyclic group and the nucleophile controls substrate conformation: its presence preserves sugar distortion, whereas its absence (e.g. in an enzyme mutant) knocks it out. We also show that 2-deoxy-2-fluoro derivatives, widely used to trap the reaction intermediates by X-ray crystallography, reproduce the conformation of the hydrolysable substrate at the experimental conditions. These results highlight the importance of the 2-OH⋅⋅⋅nucleophile interaction in substrate recognition and catalysis in endo-glycosidases and can inform mutational campaigns aimed to search for more efficient enzymes.

Glycosidase mechanisms: Sugar conformations and reactivity in endo- and exo-acting enzymes

The enzymatic breakdown of carbohydrates plays a critical role in several biological events and enables the development of sustainable processes to obtain bioproducts and biofuels. In this scenario, the design of efficient inhibitors for glycosidases that can act as drug targets and the engineering of carbohydrate-active enzymes with tailored catalytic properties is of remarkable importance. To guide rational approaches, it is necessary to elucidate enzyme molecular mechanisms, in particular understanding how the microenvironment modulates the conformational space explored by the substrate. Computer simulations, especially those based on ab initio methods, have provided a suitable atomic description of carbohydrate conformations and catalytic reactions in several glycosidase families. In this review, we will focus on how the active-site topology (pocket or cleft) and mode of cleavage (endo or exo) can affect the catalytic mechanisms adopted by glycosidases, in particular the substrate conformations along the reaction coordinate.

Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides

Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable ⁴C₁ and ¹C₄ chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the ⁴C₁ and ¹C₄ conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the ⁴C₁ conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.

Cellulose Fast Pyrolysis Activated by Intramolecular Hydrogen Bonds

The conversion of inedible biomass by fast pyrolysis is a promising route for sustainable production of renewable fuels and value-added chemicals, but low selectivity toward desired products hampers its economic viability. Understanding the molecular-level reaction pathways of biomass fast pyrolysis could be the key to overcoming this challenge. However, the effects of intramolecular and interchain hydrogen bonds near the reaction center have not been thoroughly explored. In this work, the reaction pathways and kinetics of fast pyrolysis of cellulose, a major component of biomass, were investigated using the density functional theory. A new intramolecular hydroxyl-activated mechanism is presented for cellulose activation. Our calculations incorporating noncovalent interactions accurately captured the activation energy of 50.8 kcal mol^-1, agreeable with the apparent activation energy measured experimentally. The findings of cellulose pyrolysis provide insights into the investigation of interactions during real-life biomass pyrolysis.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Characterization of the low energy conformations and differential reactivities of D-glucose and D-mannose based oxepines

Carbohydrate-based oxepines are valuable intermediates for the synthesis of septanose carbohydrates. Here we report the characterization of the preferred conformations of D-glucose and D-mannose based oxepines 1 and 2 using computational chemistry and NMR spectroscopy. Monte Carlo conformational searches on 1 and 2 were performed, followed by DFT optimization and single-point energy calculations on the low energy conformations of each oxepine. Coupling constants were computed for all unique conformations at a B3LYP/6-31G(d,p)u+1s level of theory and weighted based on a Boltzmann distribution. These values were then compared to the experimental values collected using ³J_H,H values collected from ¹H NMR spectra. Information from the MC/DFT approach was then used in a least squares method that correlated DFT calculated and observed ³J_H,H coupling constants. The conformations of 1 and 2 are largely governed by a combination of the rigidifying enol ether element in combination with the reduction of unfavorable interactions. The vinylogous anomeric effect (VAE) emerged as a consequence, rather than a driver of conformations. Oxepine 1 showed greater reactivity in Ferrier rearrangement reactions relative to oxepine 2, in line with its greater %VAE.

Glycosyl Oxocarbenium Ions: Structure, Conformation, Reactivity, and Interactions

Carbohydrates (glycans, saccharides, and sugars) are essential molecules in all domains of life. Research on glycoscience spans from chemistry to biomedicine, including material science and biotechnology. Access to pure and well-defined complex glycans using synthetic methods depends on the success of the employed glycosylation reaction. In most cases, the mechanism of the glycosylation reaction is believed to involve the oxocarbenium ion. Understanding the structure, conformation, reactivity, and interactions of this glycosyl cation is essential to predict the outcome of the reaction. In this Account, building on our contributions on this topic, we discuss the theoretical and experimental approaches that have been employed to decipher the key features of glycosyl cations, from their structures to their interactions and reactivity.We also highlight that, from a chemical perspective, the glycosylation reaction can be described as a continuum, from unimolecular SN1 with naked oxocarbenium cations as intermediates to bimolecular SN2-type mechanisms, which involve the key role of counterions and donors. All these factors should be considered and are discussed herein. The importance of dissociative mechanisms (involving contact ion pairs, solvent-separated ion pairs, solvent-equilibrated ion pairs) with bimolecular features in most reactions is also highlighted.The role of theoretical calculations to predict the conformation, dynamics, and reactivity of the oxocarbenium ion is also discussed, highlighting the advances in this field that now allow access to the conformational preferences of a variety of oxocarbenium ions and their reactivities under SN1-like conditions.Specifically, the ground-breaking use of superacids to generate these cations is emphasized, since it has permitted characterization of the structure and conformation of a variety of glycosyl oxocarbenium ions in superacid solution by NMR spectroscopy.We also pay special attention to the reactivity of these glycosyl ions, which depends on the conditions, including the counterions, the possible intra- or intermolecular participation of functional groups that may stabilize the cation and the chemical nature of the acceptor, either weak or strong nucleophile. We discuss recent investigations from different experimental perspectives, which identified the involved ionic intermediates, estimating their lifetimes and reactivities and studying their interactions with other molecules. In this context, we also emphasize the relationship between the chemical methods that can be employed to modulate the sensitivity of glycosyl cations and the way in which glycosyl modifying enzymes (glycosyl hydrolases and transferases) build and cleave glycosidic linkages in nature. This comparison provides inspiration on the use of molecules that regulate the stability and reactivity of glycosyl cations.

Rational Design of Mechanism-Based Inhibitors and Activity-Based Probes for the Identification of Retaining α-l-Arabinofuranosidases

Identifying and characterizing the enzymes responsible for an observed activity within a complex eukaryotic catabolic system remains one of the most significant challenges in the study of biomass-degrading systems. The debranching of both complex hemicellulosic and pectinaceous polysaccharides requires the production of α-l-arabinofuranosidases among a wide variety of coexpressed carbohydrate-active enzymes. To selectively detect and identify α-l-arabinofuranosidases produced by fungi grown on complex biomass, potential covalent inhibitors and probes which mimic α-l-arabinofuranosides were sought. The conformational free energy landscapes of free α-l-arabinofuranose and several rationally designed covalent α-l-arabinofuranosidase inhibitors were analyzed. A synthetic route to these inhibitors was subsequently developed based on a key Wittig–Still rearrangement. Through a combination of kinetic measurements, intact mass spectrometry, and structural experiments, the designed inhibitors were shown to efficiently label the catalytic nucleophiles of retaining GH51 and GH54 α-l-arabinofuranosidases. Activity-based probes elaborated from an inhibitor with an aziridine warhead were applied to the identification and characterization of α-l-arabinofuranosidases within the secretome of A. niger grown on arabinan. This method was extended to the detection and identification of α-l-arabinofuranosidases produced by eight biomass-degrading basidiomycete fungi grown on complex biomass. The broad applicability of the cyclophellitol-derived activity-based probes and inhibitors presented here make them a valuable new tool in the characterization of complex eukaryotic carbohydrate-degrading systems and in the high-throughput discovery of α-l-arabinofuranosidases.

Strategies for carbohydrate model building, refinement and validation

Sugars are the most stereochemically intricate family of biomolecules and present substantial challenges to anyone trying to understand their nomenclature, reactions or branched structures. Current crystallographic programs provide an abstraction layer allowing inexpert structural biologists to build complete protein or nucleic acid model components automatically either from scratch or with little manual intervention. This is, however, still not generally true for sugars. The need for carbohydrate-specific building and validation tools has been highlighted a number of times in the past, concomitantly with the introduction of a new generation of experimental methods that have been ramping up the production of protein-sugar complexes and glycoproteins for the past decade. While some incipient advances have been made to address these demands, correctly modelling and refining carbohydrates remains a challenge. This article will address many of the typical difficulties that a structural biologist may face when dealing with carbohydrates, with an emphasis on problem solving in the resolution range where X-ray crystallography and cryo-electron microscopy are expected to overlap in the next decade.

DFT/PCM theoretical study of the conversion of methyl 4-O-methyl-α-d-galactopyranoside 6-sulfate and its 2-sulfated derivative into their 3,6-anhydro counterparts

Modeling of the conversion of methyl 4-O-methyl-α-d-galactopyranoside 6-sulfate (2) and 2,6-disulfate (1) into methyl 3,6-anhydro-4-O-methyl-α-d-galactopyranoside (4) and its 2-sulfate (3), respectively (Scheme 1) has been carried out using DFT at the M06-2X/6-311 + G(d,p)//M06-2X/6-31 + G(d,p) level with the polarizable continuum model (PCM) in water. The three steps necessary for the alkaline transformation of 6-sulfated (and 2,6-disulfated) galactose units into 3,6-anhydro derivatives were evaluated. The final substitution step appears to be the rate limiting, involving an activation energy of ca. 23 kcal/mol. The other two steps (deprotonation and chair inversion) combined involve lower activation energies (9-12 kcal/mol). Comparison of the thermodynamics and kinetics of the reactions suggest that if the deprotonation step precedes the chair inversion, the reaction should be faster for both compounds. No major differences in reaction rate can be theoretically predicted to be caused by the presence of sulfate on O-2, although one experimental result suggested that O-2 sulfation should increase the reaction rate. The conformational pathways are complex, given the large number of rotamers available for each compound, and the way that some of these rotamers combine into some of the pathways. In any case, the conformation (O)S2 appears as a common intermediate for the chair inversion processes.

The loop structure of Actinomycete glycoside hydrolase family 5 mannanases governs substrate recognition

Endo-β-1,4-mannanases from Streptomyces thermolilacinus (StMan) and Thermobifida fusca (TfMan) demonstrated different substrate specificities. StMan hydrolyzed galactosylmannooligosaccharide (GGM5; 6(III) ,6(IV) -α-d-galactosyl mannopentaose) to GGM3 and M2, whereas TfMan hydrolyzed GGM5 to GGM4 and M1. To determine the region involved in the substrate specificity, we constructed chimeric enzymes of StMan and TfMan and evaluated their substrate specificities. Moreover, the crystal structure of the catalytic domain of StMan (StMandC) and the complex structure of the inactive mutant StE273AdC with M6 were solved at resolutions of 1.60 and 1.50 Å, respectively. Structural comparisons of StMandC and the catalytic domain of TfMan lead to the identification of a subsite around -1 in StMandC that could accommodate a galactose branch. These findings demonstrate that the two loops (loop7 and loop8) are responsible for substrate recognition in GH5 actinomycete mannanases. In particular, Trp281 in loop7 of StMan, which is located in a narrow and deep cleft, plays an important role in its affinity toward linear substrates. Asp310 in loop8 of StMan specifically bound to the galactosyl unit in the -1 subsite.

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.

Evaluation of Selected Classical Force Fields for Alchemical Binding Free Energy Calculations of Protein-Carbohydrate Complexes

Protein-carbohydrate recognition is crucial in many vital biological processes including host-pathogen recognition, cell-signaling, and catalysis. Accordingly, computational prediction of protein-carbohydrate binding free energies is of enormous interest for drug design. However, the accuracy of current force fields (FFs) for predicting binding free energies of protein-carbohydrate complexes is not well understood owing to technical challenges such as the highly polar nature of the complexes, anomerization, and conformational flexibility of carbohydrates. The present study evaluated the performance of alchemical predictions of binding free energies with the GAFF1.7/AM1-BCC and GLYCAM06j force fields for modeling protein-carbohydrate complexes. Mean unsigned errors of 1.1 ± 0.06 (GLYCAM06j) and 2.6 ± 0.08 (GAFF1.7/AM1-BCC) kcal·mol(-1) are achieved for a large data set of monosaccharide ligands for Ralstonia solanacearum lectin (RSL). The level of accuracy provided by GLYCAM06j is sufficient to discriminate potent, moderate, and weak binders, a goal that has been difficult to achieve through other scoring approaches. Accordingly, the protocols presented here could find useful applications in carbohydrate-based drug and vaccine developments.

Conformational Plasticity in Glycomimetics: Fluorocarbamethyl-L-idopyranosides Mimic the Intrinsic Dynamic Behaviour of Natural Idose Rings

Sugar function, structure and dynamics are intricately correlated. Ring flexibility is intrinsically related to biological activity; actually plasticity in L-iduronic rings modulates their interactions with biological receptors. However, the access to the experimental values of the energy barriers and free-energy difference for conformer interconversion in water solution has been elusive. Here, a new generation of fluorine-containing glycomimetics is presented. We have applied a combination of organic synthesis, NMR spectroscopy and computational methods to investigate the conformational behaviour of idose- and glucose-like rings. We have used low-temperature NMR spectroscopic experiments to slow down the conformational exchange of the idose-like rings. Under these conditions, the exchange rate becomes slow in the (19) F NMR spectroscopic chemical shift timescale and allows shedding light on the thermodynamic and kinetic features of the equilibrium. Despite the minimal structural differences between these compounds, a remarkable difference in their dynamic behaviour indeed occurs. The importance of introducing fluorine atoms in these sugars mimics is also highlighted. Only the use of (19) F NMR spectroscopic experiments has permitted the unveiling of key features of the conformational equilibrium that would have otherwise remained unobserved.

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.

Small endohedral metallofullerenes: exploration of the structure and growth mechanism in the Ti@C2n (2n = 26-50) family

Analysis of the structure and the bottom-up growth mechanism in the family of small endohedral metallofullerenes Ti@C_2n (2n = 26–50).

The formation of the smallest fullerene, C₂₈, was recently reported using gas phase experiments combined with high-resolution FT-ICR mass spectrometry. An internally located group IV metal stabilizes the highly strained non-IPR C₂₈ cage by charge transfer (IPR = isolated pentagon rule). Ti@C₄₄ also appeared as a prominent peak in the mass spectra, and U@C₂₈ was demonstrated to form by a bottom-up growth mechanism. We report here a computational analysis using standard DFT calculations and Car–Parrinello MD simulations for the family of the titled compounds, aiming to identify the optimal cage for each endohedral fullerene and to unravel key aspects of the intriguing growth mechanisms of fullerenes. We show that all the optimal isomers from C₂₆ to C₅₀ are linked by a simple C₂ insertion, with the exception of a few carbon cages that require an additional C₂ rearrangement. The ingestion of a C₂ unit is always an exergonic/exothermic process that can occur through a rather simple mechanism, with the most energetically demanding step corresponding to the closure of the carbon cage. The large formation abundance observed in mass spectra for Ti@C₂₈ and Ti@C₄₄ can be explained by the special electronic properties of these cages and their higher relative stabilities with respect to C₂ reactivity. We further verify that extrusion of C atoms from an already closed fullerene is much more energetically demanding than forming the fullerene by a bottom-up mechanism. Independent of the formation mechanism, the present investigations strongly support that, among all the possible isomers, the most stable, smaller non-IPR carbon cages are formed, a conclusion that is also valid for medium and large cages.

Protein engineering study of β-mannosidase to set up a potential chemically efficient biocatalyst

This study is focused on the analysis and mutagenesis of β-mannosidase from Bacteroides thetaiotaomicron with the aim of broadening its substrate specificity to 2-acetamido-2-deoxy-β-d-mannopyranosyl (β-ManNAc) derivatives. Various conformations ((4)C1, (4)H5 and (1)S5) of native and modified ligands were docked to the binding site of the protein to determine the most suitable conformation of sugars for further hydrolysis. Key amino acid residues were mutated in silico focusing on stabilizing the acetamido group of β-ManNAc as well as forming the oxazoline intermediate needed for hydrolysis. The results of large set of 5 ns molecular dynamic simulations showed that the majority of the active site residues are involved in substrate interaction and do not exhibit a higher flexibility except for Asn178. Mutations of Asn178 to alanine and Asp199 to serine could lead to a stabilization of the acetamido group in the binding site. So far, in vitro mutagenesis and the screen of a large variety of biological sources were unable to extend β-mannosidase's activity to include β-ManNAc derivatives.

Molecular mechanism of a hotdog-fold acyl-CoA thioesterase

Thioesterases are enzymes that hydrolyze thioester bonds between a carbonyl group and a sulfur atom. They catalyze key steps in fatty acid biosynthesis and metabolism, as well as polyketide biosynthesis. The reaction molecular mechanism of most hotdog-fold acyl-CoA thioesterases remains unknown, but several hypotheses have been put forward in structural and biochemical investigations. The reaction of a human thioesterase (hTHEM2), representing a thioesterase family with a hotdog fold where a coenzyme A moiety is cleaved, was simulated by quantum mechanics/molecular mechanics metadynamics techniques to elucidate atomic and electronic details of its mechanism, its transition-state conformation, and the free energy landscape of the process. A single-displacement acid-base-like mechanism, in which a nucleophilic water molecule is activated by an aspartate residue acting as a base, was found, confirming previous experimental proposals. The results provide unambiguous evidence of the formation of a tetrahedral-like transition state. They also explain the roles of other conserved active-site residues during the reaction, especially that of a nearby histidine/serine pair that protonates the thioester sulfur atom, the participation of which could not be elucidated from mutation analyses alone.

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.

Combined inhibitor free-energy landscape and structural analysis reports on the mannosidase conformational coordinate

Mannosidases catalyze the hydrolysis of a diverse range of polysaccharides and glycoconjugates, and the various sequence-based mannosidase families have evolved ingenious strategies to overcome the stereoelectronic challenges of mannoside chemistry. Using a combination of computational chemistry, inhibitor design and synthesis, and X-ray crystallography of inhibitor/enzyme complexes, it is demonstrated that mannoimidazole-type inhibitors are energetically poised to report faithfully on mannosidase transition-state conformation, and provide direct evidence for the conformational itinerary used by diverse mannosidases, including β-mannanases from families GH26 and GH113. Isofagomine-type inhibitors are poor mimics of transition-state conformation, owing to the high energy barriers that must be crossed to attain mechanistically relevant conformations, however, these sugar-shaped heterocycles allow the acquisition of ternary complexes that span the active site, thus providing valuable insight into active-site residues involved in substrate recognition.

Structural analysis and insights into the glycon specificity of the rice GH1 Os7BGlu26 β-D-mannosidase

Rice Os7BGlu26 is a GH1 family glycoside hydrolase with a threefold higherkcat/Kmvalue for 4-nitrophenyl β-D-mannoside (4NPMan) compared with 4-nitrophenyl β-D-glucoside (4NPGlc). To investigate its selectivity for β-D-mannoside and β-D-glucoside substrates, the structures of apo Os7BGlu26 at a resolution of 2.20 Å and of Os7BGlu26 with mannose at a resolution of 2.45 Å were elucidated from isomorphous crystals in space groupP212121. The (β/α)8-barrel structure is similar to other GH1 family structures, but with a narrower active-site cleft. The Os7BGlu26 structure with D-mannose corresponds to a product complex, with β-D-mannose in the1S5skew-boat conformation. Docking of the1S3,1S5,2SOand3S1pyranose-ring conformations of 4NPMan and 4NPGlc substrates into the active site of Os7BGlu26 indicated that the lowest energies were in the1S5and1S3skew-boat conformations. Comparison of these docked conformers with other rice GH1 structures revealed differences in the residues interacting with the catalytic acid/base between enzymes with and without β-D-mannosidase activity. The mutation of Tyr134 to Trp in Os7BGlu26 resulted in similarkcat/Kmvalues for 4NPMan and 4NPGlc, while mutation of Tyr134 to Phe resulted in a 37-fold higherkcat/Kmfor 4NPMan than 4NPGlc. Mutation of Cys182 to Thr decreased both the activity and the selectivity for β-D-mannoside. It was concluded that interactions with the catalytic acid/base play a significant role in glycon selection.

String method for calculation of minimum free-energy paths in Cartesian space in freely-tumbling systems

The string method is a molecular-simulation technique that aims to calculate the minimum free-energy path of a chemical reaction or conformational transition, in the space of a pre-defined set of reaction coordinates that is typically highly dimensional. Any descriptor may be used as a reaction coordinate, but arguably the Cartesian coordinates of the atoms involved are the most unprejudiced and intuitive choice. Cartesian coordinates, however, present a non-trivial problem, in that they are not invariant to rigid-body molecular rotations and translations, which ideally ought to be unrestricted in the simulations. To overcome this difficulty, we reformulate the framework of the string method to integrate an on-the-fly structural-alignment algorithm. This approach, referred to as SOMA (String method with Optimal Molecular Alignment), enables the use of Cartesian reaction coordinates in freely tumbling molecular systems. In addition, this scheme permits the dissection of the free-energy change along the most probable path into individual atomic contributions, thus revealing the dominant mechanism of the simulated process. This detailed analysis also provides a physically-meaningful criterion to coarse-grain the representation of the path. To demonstrate the accuracy of the method we analyze the isomerization of the alanine dipeptide in vacuum and the chair-to-inverted-chair transition of β-D mannose in explicit water. Notwithstanding the simplicity of these systems, the SOMA approach reveals novel insights into the atomic mechanism of these isomerizations. In both cases, we find that the dynamics and the energetics of these processes are controlled by interactions involving only a handful of atoms in each molecule. Consistent with this result, we show that a coarse-grained SOMA calculation defined in terms of these subsets of atoms yields nearidentical minimum free-energy paths and committor distributions to those obtained via a highly-dimensional string.

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.

Is N-acetyl-D-glucosamine a rigid 4C1 chair?

Understanding microsecond-timescale dynamics is crucial to establish three-dimensional (3D) structure-activity relationships in sugars but has been intractable to experiments and simulations. As a consequence, whether arguably the most important chemical scaffold in glycobiology, N-acetyl-d-glucosamine (GlcNAc), deviates from a rigid (4)C(1) chair is unknown. Here, conformer populations and exchange kinetics were quantified from the longest aqueous carbohydrate simulations to date (0.2 ms total) of GlcNAc, four derivatives from heparan sulfate and their methylglycosides. Unmodified GlcNAc took 3-5 μs to reach a conformational equilibrium, which comprised a metastable (4)C(1) chair that underwent (4)C(1) ↔ (1)C(4) transitions at a predicted forward rate of 0.8 μs(-1) with an average (1)C(4)-chair lifetime of 3 ns. These predictions agree with high-resolution crystallography and nuclear magnetic resonance but not with the hypothesis that GlcNAc is a rigid (4)C(1) chair, concluded from previous experimental analyses and non-aqueous modeling. The methylglycoside was calculated to have a slower forward rate (0.3 μs(-1)) and a more stable (4)C(1) conformer (0.2 kcal mol(-1)), suggesting that pivotal 3D intermediates (particularly (2)S(O), (1)S(5) and B(2,5)) increased in energy, and water was implicated as a major cause. Sulfonation (N-, 3-O and 6-O) significantly augmented this effect by blocking pseudorotation, but did not alter the rotational preferences of hydroyxl or hydroxymethyl groups. We therefore propose that GlcNAc undergoes puckering exchange that is dependent on polymerization and sulfo substituents. Our analyses, and 3D model of the equilibrium GlcNAc conformer in water, can be used as dictionary data and present new opportunities to rationally modify puckering and carbohydrate bioactivity, with diverse applications from improving crop yields to disease amelioration.