Modelling of β-d-glucopyranose ring distortion in different force fields: a metadynamics study

High Yield of L-Sorbose via D-Glucose Isomerization in Ethanol over a Bifunctional Titanium-Boron-Beta Zeolite

D-Glucose-to-L-sorbose isomerization on Lewis acidic zeolite is a highly attractive avenue for sorbose production. But the L-sorbose yield is limited by the competing D-glucose-to-D-fructose isomerization and reaction equilibrium. In this work, it is suggested that ethanol directs the glucose conformation for selective D-glucose-to-L-sorbose isomerization. It also reacts with the produced L-sorbose to form ethyl-sorboside, which allows the D-glucose-to-L-sorbose isomerization to proceed beyond the thermodynamic equilibrium limit. It is shown that a bifunctional zeolite Beta containing framework titanium (Ti) and boron (B) is a selective catalyst for this tandem reaction: Lewis acidic framework Ti catalyzes the D-glucose-to-L-sorbose isomerization via an intramolecular 5,1-hydride shift process as confirmed by isotopic tracing experiments followed by 13C-NMR, while weak Brønsted acid framework B selectively promotes the sorbose ketalization with ethanol. A remarkably high yield of L-sorbose with a high fraction of sugar (>95 %: 27 % unreacted glucose, 11.4 % fructose, 57 % sorbose) was obtained after the mixture produced in ethanol was hydrolyzed.

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

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.

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.

Carbohydrate Force Fields: The Role of Small Partial Atomic Charges in Preventing Conformational Collapse

Although the quality of current additive all-atom force fields for carbohydrates has been demonstrated in many applications, occasional significant differences reported for the hydrodynamic behavior of specific polysaccharides modeled with different force fields is a cause for concern. In particular, irreversible conformational collapse has been noted for some polysaccharide simulations with the GLYCAM06j force field. Here, we investigate the cause of this phenomenon through comparative simulations of a range of saccharides with both the GLYCAM06j and the CHARMM36 carbohydrate force fields. We find that conformational collapse in GLYCAM06j occurs for saccharide chains containing the deoxy sugar α-l-rhamnose after relatively long simulation intervals. Further, we explore the mechanism of conformational collapse and show that this phenomenon arises because of the anomalous low energy in GLYCAM06j (as compared to quantum mechanical calculations) of a specific orientation of α-l-Rha to α-l-Rha glycosidic linkages, which are subsequently sustained by intramolecular interactions in the saccharide chain. We identify the lack of partial charges on aliphatic hydrogens in GLYCAM as the source of this anomaly, demonstrating that addition of small partial atomic charges on the aliphatic protons in rhamnose removes the conformational collapse phenomenon. This work reveals the large cumulative impact that small partial charges may have on the dynamic behavior of polysaccharides and indicates that future reparameterization of the GLYCAM06j force field should investigate the addition of partial charges on all aliphatic hydrogens.

Combining Computational Chemistry and Crystallography for a Better Understanding of the Structure of Cellulose

The approaches in this article seek to enhance understanding of cellulose at the molecular level, independent of the source and the particular crystalline form of cellulose. Four main areas of structure research are reviewed. Initially, the molecular shape is inferred from the crystal structures of many small molecules that have β-(1→4) linkages. Then, conformational analyses with potential energy calculations of cellobiose are covered, followed by the use of Atoms-In-Molecules theory to learn about interactions in experimental and theoretical structures. The last section covers models of cellulose nanoparticles. Controversies addressed include the stability of twofold screw-axis conformations, the influence of different computational methods, the predictability of crystalline conformations by studies of isolated molecules, and the twisting of model cellulose crystals.

Improvement of the Force Field for β-d-Glucose with Machine Learning

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β-d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β-d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly.

Structural Determinants of Substrate Recognition and Catalysis by Heparan Sulfate Sulfotransferases

Heparan sulfate (HS) and heparin contain imprinted "sulfation codes", which dictate their diverse physiological and pathological functions. A group of orchestrated biosynthetic enzymes cooperate in polymerizing and modifying HS chains. The biotechnological development of enzymes that can recreate this sulfation pattern on synthetic heparin is challenging, primarily due to the paucity of quantitative data for sulfotransferase enzymes. Herein, we identified critical structural characteristics that determine substrate specificity and shed light on the catalytic mechanism of sugar sulfation of two HS sulfotransferases, 2-O-sulfotransferase (HS2ST) and 6-O-sulfotransferase (HS6ST). Two sets of molecular clamps in HS2ST recognize appropriate substrates; these clamps flank the acceptor binding site on opposite sides. The hexuronic epimers, and not their puckers, have a critical influence on HS2ST selectivity. In contrast, HS6ST recognizes a broader range of substrates. This promiscuity is granted by a conserved tryptophan residue, W210, that positions the acceptor within the active site for catalysis by means of strong electrostatic interactions. Lysines K131 and K132 act in concert with a second tryptophan, W153, shedding water molecules from within the active site, thus providing HS6ST with a binding preference toward 2-O-sulfated substrates. QM/MM calculations provided valuable mechanistic insights into the catalytic process, identifying that the sulfation of both HS2ST and HS6ST follows a SN2-like mechanism. When they are taken together, our findings reveal the molecular basis of how these enzymes recognize different substrates and catalyze sugar sulfation, enabling the generation of enzymes that could create specific heparin epitopes.

Molecular Dynamics Simulation of Cellulose Synthase Subunit D Octamer with Cellulose Chains from Acetic Acid Bacteria: Insight into Dynamic Behaviors and Thermodynamics on Substrate Recognition

The present study reports the building of a computerized model and molecular dynamics (MD) simulation of cellulose synthase subunit D octamer (CesD) from Komagataeibacter hansenii. CesD was complexed with four cellulose chains having DP = 12 (G12) by model building, which revealed unexpected S-shaped pathways with bending regions. Combined conventional and accelerated MD simulations of CesD complex models were carried out, while the pyranose ring conformations of the glucose residues were restrained to avoid undesirable deviations of the ring conformation from the ⁴C₁ form. The N-terminal regions and parts of the secondary structures of CesD established appreciable contacts with the G12 chains. Hybrid quantum mechanical (QM) and molecular mechanical (MM) simulations of the CesD complex model were performed. Glucose residues located at the pathway bends exhibited reversible changes to the ring conformation into either skewed or boat forms, which might be related to the function of CesD in regulating microfibril production.

Thermal Response in Cellulose Iβ Based on Molecular Dynamics

The structural details of cellulose I β were discussed according to molecular dynamics simulations with the GLYCAM-06 force field. The simulation outcomes were in agreement with previous experimental data, including structural parameters and hydrogen bond pattern at 298 K. We found a new conformation of cellulose Iβ existed at the intermediate temperature that is between the low and high temperatures. Partial chain rotations along the backbone direction were found and conformations of hydroxymethyl groups that alternated from tg to either gt or gg were observed when the temperature increased from 298 K to 400 K. In addition, the gg conformation is preferred than gt. For the structure adopted at high temperature of 500 K, major chains were twisted and two chains detached from each plain. In contrast to the observation under intermediate temperature, the population of hydroxymethyl groups in gt exceeded that in gg conformation at high temperature. In addition, three patterns of hydrogen bonding were identified at low, intermediate and high temperatures in the simulations. The provided structural information indicated the transitions occurred around 350 K and 450 K, considered as the transitional temperatures of cellulose Iβ in this work.

Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations

The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 μs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on β-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of β-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the ¹C₄ ring conformation; this observation may have implications for understanding the structural basis of the biological function of β-GalNAc-containing glycosaminoglycans such as dermatan sulfate.

Perspective on computational simulations of glycosaminoglycans

Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.

Pyranose ring puckering in aldopentoses, ketohexoses and deoxyaldohexoses. A molecular dynamics study

Conformation of monosaccharides, including the ring shape, has for years been the subject of intensive research. Although d-aldohexopyranoses are the most extensively studied pyranoses, there also exist other groups of saccharides that contain analogous chemical system of the six-membered ring. Here we describe in details the results of the molecular dynamics-based conformational analysis concerning a series of pyranoses, namely: d-aldopentoses, d-ketohexoses as well as deoxy- (d-quinovose, l-fucose, l-rhamnose) and dideoxy- (abequose, paratose, tyvelose, digitoxose) derivatives of aldohexoses. By using the carbohydrate-dedicated GROMOS 56a6_CARBO force field, we determined the conformational properties of both the lactol and hydroxymethyl groups as well as the anomeric populations for all considered compounds. The orientation of the lactol group follows the trend expected on the basis of the exo-anomeric effect for all compounds whereas the conformation of the hydroxymethyl group in d-ketohexoses is represented by the two gauche (with respect to the ring oxygen atom) rotamers. The special emphasis is put on the ring-inversion properties studied in the context of both the full chair-chair inversion and the chair-boat/skew-boat rearrangement. The calculated ring-distortion energies, compared with those obtained for regular d-aldohexopyranoses allowed for estimating the influence of particular substituents on the ring flexibility. Overall, such influence is correlated with the dimension of the substituent and its orientation but is limited to the case of the chair-chair inversion whereas the chair-to-boat/skew-boat rearrangement exhibits roughly the same properties for all pyranoses. For all d-aldopyranoses the α anomers exhibit lower ring-inversion free energies in comparison to the β anomers whereas this trend is inverted in the case of d-ketohexopyranoses.

Molecular dynamics simulations of hexopyranose ring distortion in different force fields

The four classical, biomolecular force fields designed to study hexopyranose-based carbohydrates (GROMOS 56a6CARBO/56a6CARBO_R, GROMOS 53a6GLYC, CHARMM and GLYCAM06) have been tested in the context of ring-inversion properties. These properties were evaluated for both unfunctionalized monomers of all hexopyranoses of the d series and for residues in a chain composed of uniform units connected by α(1→4) and β(1→4) glycosidic linkages. The results indicate that the tested force fields differ in their predictions of the ring-inversion properties of both monomers and residues in a chain. The comparison with the available experimental data and with the semi-empirical Angyal scheme reveals that, at the level of monomers, GROMOS 56a6CARBO, GROMOS 53a6GLYC and CHARMM correctly reproduce the ring-inversion free energies. However, due to the lack of analogous reference data we cannot state which force field is more or less accurate in the context of ring distortion of residues in a chain. Therefore, the use of ab initio potentials is recommended in the prospective, quantitative studies on the related subject.

Effects of varying the 6-position oxidation state of hexopyranoses: a systematic comparative computational analysis of 48 monosaccharide stereoisomers

Knowledge of multi-dimensional carbohydrate structure is essential when delineating structure-function relationships in the development of analytical techniques such as ion mobility-mass spectrometry and of carbohydrate-based therapeutics, as well as in rationally modifying the chemical and physical properties of drugs and materials based on sugars. Although monosaccharides are conventionally presumed to adopt the canonical ⁴C₁ chair conformation, it is not well known how altering the substituent identity around the pyranose ring affects the favored conformational state. This work provides a comprehensive and systematic computational comparison of all eight aldohexose isomers in the gas phase with reduction and oxidation at the C-6 position using density functional theory (M05-2X/cc-pVTZ(-f)//B3LYP/6-31G**) to determine the conformational and anomeric preference for each sugar in the gas phase. All 6-deoxyhexose and aldohexose isomers favored the ⁴C₁ chair conformation, while oxidation at C-6 showed a shift in equilibrium to favor the ¹C₄ chair for β-alluronic acid, β-guluronic acid, and β-iduronic acid. The anomeric preference was found to be significantly affected by a remote change in oxidation state, with the alternate anomer favored for several isomers. These findings provide a fundamental platform to empirically test steric and electronic effects of pyranose substituents, with the goal of formulating straightforward rules that govern carbohydrate reactivity and drive quicker, more efficient syntheses. Graphical abstract A systematic comparative conformational analysis of all eight aldohexose isomers using DFT methods (M05-2X/cc-pVTZ(-f)) reveals changes in anomeric and ring conformational preference upon reduction or oxidation at the C-6 position for several sugars.

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.

Conformational sampling of oligosaccharides using Hamiltonian replica exchange with two-dimensional dihedral biasing potentials and the weighted histogram analysis method (WHAM)

Oligosaccharides and polysaccharides exert numerous functional roles in biology through their structural diversity and conformational properties. To investigate their conformational properties using computational methods, Hamiltonian replica exchange (H-REX) combined with two-dimensional grid-based correction maps as biasing potentials (bpCMAP) significantly improves the sampling efficiency about glycosidic linkages. In the current study, we extend the application of H-REX with bpCMAP to complex saccharides and establish systematic procedures for bpCMAP construction, determination of replica distribution, and data analysis. Our main findings are that (1) the bpCMAP for each type of glycosidic linkage can be constructed from the corresponding disaccharide using gas-phase umbrella sampling simulations, (2) the replica distribution can be conveniently determined following the exact definition of the average acceptance ratio based on the assigned distribution of biasing potentials, and (3) the extracted free energy surface (or potential of mean force (PMF)) can be improved using the weighted histogram analysis method (WHAM) allowing for the inclusion of data from the excited state replicas in the calculated probability distribution. The method is applied to a branched N-glycan found on the HIV gp120 protein, and a linear N-glycan. Considering the general importance of N-glycans and the wide appreciation of the sampling problem, the present method represents an efficient procedure for the conformational sampling of complex oligo- and polysaccharides under explicit solvent conditions. More generally, the use of WHAM is anticipated to be of general utility for the calculation of PMFs from H-REX simulations in a wide range of macromolecular systems.

Cationic Peptides and Peptidomimetics Bind Glycosaminoglycans as Potential Sema3A Pathway Inhibitors

Semaphorin3A (Sema3A) is a vertebrate-secreted protein that was initially characterized as a repulsive-guidance cue. Semaphorins have crucial roles in several diseases; therefore, the development of Sema3A inhibitors is of therapeutic interest. Sema3A interacts with glycosaminoglycans (GAGs), presumably through its C-terminal basic region. We used different biophysical techniques (i.e., NMR, surface plasmon resonance, isothermal titration calorimetry, fluorescence, and UV-visible spectroscopy) to characterize the binding of two Sema3A C-terminus-derived basic peptides (FS2 and NFS3) to heparin and chondroitin sulfate A. We found that these peptides bind to both GAGs with affinities in the low-micromolar range. On the other hand, a peptoid named SICHI (semaphorin-induced chemorepulsion inhibitor), which is positively charged at physiological pH, was first identified by our group as being able to block Sema3A chemorepulsion and growth-cone collapse in axons at the extracellular level. To elucidate the direct target for the reported SICHI inhibitory effect in the Sema3A signaling pathway, we looked first to the protein-protein interaction between secreted Sema3A and the Nrp1 receptor. However, our results show that SICHI does not bind directly to the Sema3A sema domain or to Nrp1 extracellular domains. We evaluated a new, to our knowledge, hypothesis, according to which SICHI binds to GAGs, thereby perturbing the Sema3A-GAG interaction. By using the above-mentioned techniques, we observed that SICHI binds to GAGs and competes with Sema3A C-terminus-derived basic peptides for binding to GAGs. These data support the ability of SICHI to block the biologically relevant interaction between Sema3A and GAGs, thus revealing SICHI as a new, to our knowledge, class of inhibitors that target the GAG-protein interaction.

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.

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.

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.

Computational study of synthetic agonist ligands of ionotropic glutamate receptors

Neurological glutamate receptors are among the most important and intensely studied protein ligand binding systems in humans. They are crucial for the functioning of the central nervous system and involved in a variety of pathologies. Apart from the neurotransmitter glutamate, several artificial, agonistic and antagonistic ligands are known. Of particular interest here are novel photoswitchable agonists that would open the field of optogenetics to glutamate receptors. The receptor proteins are complex, membrane-bound multidomain oligomers that undergo large scale functional conformational changes, making detailed studies of their atomic structure challenging. Therefore, a thorough understanding of the microscopic details of ligand binding and receptor activation remains elusive in many cases. This topic has been successfully addressed by theoretical studies in the past and in this paper, we present extensive molecular dynamics simulation and free energy calculation results on the binding of AMPA and an AMPA derivative, which is the basis for designing light-sensitive ligands. We provide a two-step model for ligand binding domain activation and predict binding free energies for novel compounds in good agreement to experimental observations.

Toward an accurate conformational modeling of iduronic acid

Iduronic acid (IdoA), unlike most other monosaccharides, can adopt different ring conformations, depending on the context of the molecular structure. Accurate modeling of this building block is essential for understanding the role of glycosaminoglycans and other glycoconjugates. Here, we use metadynamics to predict equilibria of (1)C(4), (4)C(1) and (2)S(O) conformations of α-L-IdoA-OMe and α-L-IdoA2S-OMe. Different schemes of scaling of atoms separated by three bonds (1-4 interaction) were tested. It was found that scaling (reduction) of 1-4 electrostatic interactions significantly changes conformational preferences toward the (4)C(1) conformation. More interestingly, scaling of 1-4 van der Waals interaction favors skew-boat conformations. This shows that a minor modification of noncovalent 1-4 interactions parameters can provide a good agreement between populations of conformers of iduronic acid in water from simulations and experiments.

Simulation of carbohydrates, from molecular docking to dynamics in water

Modeling of carbohydrates is particularly challenging because of the variety of structures resulting for the high number of monosaccharides and possible linkages and also because of their intrinsic flexibility. The development of carbohydrate parameters for molecular modeling is still an active field. Nowadays, main carbohydrates force fields are GLYCAM06, CHARMM36, and GROMOS 45A4. GLYCAM06 includes the largest choice of compounds and is compatible with the AMBER force fields and associated. Furthermore, AMBER includes tools for the implementation of new parameters. When looking at protein-carbohydrate interaction, the choice of the starting structure is of importance. Such complex can be sometimes obtained from the Protein Data Bank-although the stereochemistry of sugars may require some corrections. When no experimental data is available, molecular docking simulation is generally used to the obtain protein-carbohydrate complex coordinates. As molecular docking parameters are not specifically dedicated to carbohydrates, inaccuracies should be expected, especially for the docking of polysaccharides. This issue can be addressed at least partially by combining molecular docking with molecular dynamics simulation in water.

GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates

An improved parameter set for explicit-solvent simulations of carbohydrates (referred to as GROMOS 53A6GLYC) is presented, allowing proper description of the most stable conformation of all 16 possible aldohexopyranose-based monosaccharides. This set includes refinement of torsional potential parameters associated with the determination of hexopyranose rings conformation by fitting to their corresponding quantum-mechanical profiles. Other parameters, as the rules for third and excluded neighbors, are taken directly from the GROMOS 53A6 force field. Comparisons of the herein presented parameter set to our previous version (GROMOS 45A4), the GLYCAM06 force field, and available NMR data are presented in terms of ring puckering free energies, conformational distribution of the hydroxymethyl group, and glycosidic linkage geometries for 16 selected monosaccharides and eight disaccharides. The proposed parameter modifications have shown a significant improvement for the above-mentioned quantities over the two tested force fields, while retaining full compatibility with the GROMOS 53A6 and 54A7 parameter sets for other classes of biomolecules.

Dependence of pyranose ring puckering on anomeric configuration: methyl idopyranosides

In the aldohexopyranose idose, the unique presence of three axial ring hydroxyl groups causes considerable conformational flexibility, rendering it challenging to study experimentally and an excellent model for rationalizing the relationship between puckering and anomeric configuration. Puckering in methyl α- and β-L-idopyranosides was predicted from kinetically rigorous 10 μs simulations using GLYCAM11 and three explicit water models (TIP3P, TIP4P, and TIP4P-EW). In each case, computed pyranose ring three-bond (vicinal) (1)H-(1)H spin couplings ((3)J(H,H)) trended with NMR measurements. These values, calculated puckering exchange rates and free energies, were independent of the water model. The α- and β-anomers were (1)C(4) chairs for 85 and >99% of their respective trajectories and underwent (1)C(4)→(4)C(1) exchange at rates of 20 μs(-1) and 1 μs(-1). Computed α-anomer (1)C(4)↔(4)C(1) puckering rates depended on the exocyclic C6 substituent, comparing hydroxymethyl with carboxyl from previous work. The slower kinetics and restricted pseudorotational profile of the β-anomer were caused by water occupying a cavity bounded by the anomeric 1-O-methyl and the C6 hydroxymethyl groups. This finding rationalizes the different methyl α- and β-L-idopyranoside (3)J(H,H) values. Identifying a relationship between idopyranose anomeric configuration, microsecond puckering, and water structure facilitates engineering of biologically and commercially important derivatives and underpins deciphering presently elusive structure-function relationships in the glycome.

Synthesis and properties of septanose carbohydrates

Seven-atom ring sugars, called septanoses, are increasingly the focus of scientific inquiries because of their potential biological activities. This article details the synthesis, conformational analysis, and protein-binding properties of septanose carbohydrates. A distinction is drawn between septanoses that are substituted in the 6-position of the ring and those that are not. When a C-6 substituent is absent, the structure is essentially that of an aldohexose in its septanose, rather than furanose or pyranose, ring form; they may play as-of-yet unexplored roles in glycobiology. Septanoses having a hydroxymethyl group at C-6, on the other hand, are ring-expanded analogues of pyranoses. Syntheses have moved beyond the preparation of seven-membered ring monosaccharides to the development of septanosyl donors. These donors have been used in the synthesis of novel di- and trisaccharides that contain septanoses as well as a variety of glycoconjugates. Low-energy conformations adopted by septanoses have been organized based on ring substitution and stereochemistry. Instances where septanoses have been demonstrated to bind to natural proteins are presented and analyzed. The major conclusion drawn in the chapter is that advances in the synthesis of septanose carbohydrates now enable a detailed investigation of their activity in a number of biological contexts.

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.