Sugar Folding: A Novel Structural Prediction Tool for Oligosaccharides and Polysaccharides 2

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Comparison of Methods for Bulk Automated Simulation of Glycosidic Bond Conformations

Six empirical force fields were tested for applicability to calculations for automated carbohydrate database filling. They were probed on eleven disaccharide molecules containing representative structural features from widespread classes of carbohydrates. The accuracy of each method was queried by predictions of nuclear Overhauser effects (NOEs) from conformational ensembles obtained from 50 to 100 ns molecular dynamics (MD) trajectories and their comparison to the published experimental data. Using various ranking schemes, it was concluded that explicit solvent MM3 MD yielded non-inferior NOE accuracy with newer GLYCAM-06, and ultimately PBE0-D3/def2-TZVP (Triple-Zeta Valence Polarized) Density Functional Theory (DFT) simulations. For seven of eleven molecules, at least one empirical force field with explicit solvent outperformed DFT in NOE prediction. The aggregate of characteristics (accuracy, speed, and compatibility) made MM3 dynamics with explicit solvent at 300 K the most favorable method for bulk generation of disaccharide conformation maps for massive database filling.

Recognition of Dimeric Lewis X by Anti-Dimeric Lex Antibody SH2

The carbohydrate antigen dimeric Lewis X (DimLe^x), which accumulates in colonic and liver adenocarcinomas, is a valuable target to develop anti-cancer therapeutics. Using the native DimLe^x antigen as a vaccine would elicit an autoimmune response against the Le^x antigen found on normal, healthy cells. Thus, we aim to study the immunogenic potential of DimLe^x and search internal epitopes displayed by DimLe^x that remain to be recognized by anti-DimLe^x monoclonal antibodies (mAbs) but no longer possess epitopes recognized by anti-Le^x mAbs. In this context, we attempted to map the epitope recognized by anti-DimLe^x mAb SH2 by titrations and competitive inhibition experiments using oligosaccharide fragments of DimLe^x as well as Le^x analogues. We compare our results with that reported for anti-Le^x mAb SH1 and anti-polymeric Le^x mAbs 1G5F6 and 291-2G3-A. While SH1 recognizes an epitope localized to the non-reducing end Le^x trisaccharide, SH2, 1G5F6, and 291-2G3-A have greater affinity for DimLe^x conjugates than for Le^x conjugates. We show, however, that the Le^x trisaccharide is still an important recognition element for SH2, which (like 1G5F6 and 291-2G3-A) makes contacts with all three sugar units of Le^x. In contrast to mAb SH1, anti-polymeric Le^x mAbs make contact with the GlcNAc acetamido group, suggesting that epitopes extend further from the non-reducing end Le^x. Results with SH2 show that this epitope is only recognized when DimLe^x is presented by glycoconjugates. We have reported that DimLe^x adopts two conformations around the β-d-GlcNAc-(1→3)-d-Gal bond connecting the Le^x trisaccharides. We propose that only one of these conformations is recognized by SH2 and that this conformation is favored when the hexasaccharide is presented as part of a glycoconjugate such as DimLe^x-bovine serum albumin (DimLe^x-BSA). Proper presentation of the oligosaccharide candidate via conjugation to a protein or lipid is essential for the design of an anti-cancer vaccine or immunotherapeutic based on DimLe^x.

Interpreting atomic force microscopy nanoindentation of hierarchical biological materials using multi-regime analysis

We present a novel Multi-Regime Analysis (MRA) routine for interpreting force indentation measurements of soft materials using atomic force microscopy. The MRA approach combines both well established and semi-empirical theories of contact mechanics within a single framework to deconvolute highly complex and non-linear force-indentation curves. The fundamental assumption in the present form of the model is that each structural contribution to the mechanical response acts in series with other 'mechanical resistors'. This simplification enables interpretation of the micromechanical properties of materials with hierarchical structures and it allows automated processing of large data sets, which is particularly indispensable for biological systems. We validate the algorithm by demonstrating for the first time that the elastic modulus of polydimethylsiloxane (PDMS) films is accurately predicted from both approach and retraction branches of force-indentation curves. For biological systems with complex hierarchical structures, we show the unique capability of MRA to map the micromechanics of live plant cells, revealing an intricate sequence of mechanical deformations resolved with precision that is unattainable using conventional methods of analysis. We recommend the routine use of MRA to interpret AFM force-indentation measurements for other complex soft materials including mammalian cells, bacteria and nanomaterials.

Enhanced conformational sampling of carbohydrates by Hamiltonian replica-exchange simulation

Knowledge of the structure and conformational flexibility of carbohydrates in an aqueous solvent is important to improving our understanding of how carbohydrates function in biological systems. In this study, we extend a variant of the Hamiltonian replica-exchange molecular dynamics (MD) simulation to improve the conformational sampling of saccharides in an explicit solvent. During the simulations, a biasing potential along the glycosidic-dihedral linkage between the saccharide monomer units in an oligomer is applied at various levels along the replica runs to enable effective transitions between various conformations. One reference replica runs under the control of the original force field. The method was tested on disaccharide structures and further validated on biologically relevant blood group B, Lewis X and Lewis A trisaccharides. The biasing potential-based replica-exchange molecular dynamics (BP-REMD) method provided a significantly improved sampling of relevant conformational states compared with standard continuous MD simulations, with modest computational costs. Thus, the proposed BP-REMD approach adds a new dimension to existing carbohydrate conformational sampling approaches by enhancing conformational sampling in the presence of solvent molecules explicitly at relatively low computational cost.

Sucrose in aqueous solution revisited, Part 1: molecular dynamics simulations and direct and indirect dipolar coupling analysis

Although the crystal structure of the disaccharide sucrose was solved more than 30 years ago, its conformational distribution in aqueous solution is still a matter of debate. We report here a variety of molecular dynamics simulations (mostly of 100 ns) using the GLYCAM06 force field and various water models, paying particular attention to comparisons to NMR measurements of residual dipolar couplings and electron-mediated spin-spin couplings. We focus on the glycosidic linkage conformation, the puckering phase angle of the fructose ring, and intramolecular hydrogen bonds between the two sugars. Our results show that sucrose is indeed a dynamic molecule, but the crystal conformation is qualitatively the dominant one in dilute solution. A second conformational basin, populated in many force fields, is probably overstabilized in the calculations.

Sucrose in aqueous solution revisited, Part 2: adaptively biased molecular dynamics simulations and computational analysis of NMR relaxation

We report 100 ns molecular dynamics simulations, at various temperatures, of sucrose in water (with concentrations of sucrose ranging from 0.02 to 4M), and in a 7:3 water-DMSO mixture. Convergence of the resulting conformational ensembles was checked using adaptive-biased simulations along the glycosidic Φ and ψ torsion angles. NMR relaxation parameters, including longitudinal (R₁) and transverse (R₂) relaxation rates, nuclear Overhauser enhancements (NOE), and generalized order parameter (S²) were computed from the resulting time-correlation functions. The amplitude and time scales of molecular motions change with temperature and concentration in ways that track closely with experimental results, and are consistent with a model in which sucrose conformational fluctuations are limited (with 80-90% of the conformations having ϕ-ψ values within 20° of an average conformation), but with some important differences in conformation between pure water and DMSO-water mixtures.

Searching and optimizing structure ensembles for complex flexible sugars

NMR restrictions are suitable to specify the geometry of a molecule when a single well-defined global free energy minimum exists that is significantly lower than other local minima. Carbohydrates are quite flexible, and therefore, NMR observables do not always correlate with a single conformer but instead with an ensemble of low free energy conformers that can be accessed by thermal fluctuations. In this communication, we describe a novel procedure to identify and weight the contribution to the ensemble of local minima conformers based on comparison to residual dipolar couplings (RDCs) or other NMR observables, such as scalar couplings. A genetic algorithm is implemented to globally minimize the R factor comparing calculated RDCs to experiment. This is done by optimizing the weights of different conformers derived from the exhaustive local minima conformational search program, fast sugar structure prediction software (FSPS). We apply this framework to six human milk sugars, LND-1, LNF-1, LNF-2, LNF-3, LNnT, and LNT, and are able to determine corresponding population weights for the ensemble of conformers. Interestingly, our results indicate that in all cases the RDCs can be well represented by only a few most important conformers. This confirms that several, but not all of the glycosidic linkages in histo-blood group "epitopes" are quite rigid.

Molecular conformations in the pentasaccharide LNF-1 derived from NMR spectroscopy and molecular dynamics simulations

The conformational dynamics of the human milk oligosaccharide lacto-N-fucopentaose (LNF-1), α-L-Fucp-(1 → 2)-β-D-Galp-(1 → 3)-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, has been analyzed using NMR spectroscopy and molecular dynamics (MD) computer simulations. Employing the Hadamard (13)C-excitation technique and the J-HMBC experiment, (1)H,(13)C trans-glycosidic J coupling constants were obtained, and from one- and two-dimensional (1)H,(1)H T-ROESY experiments, proton-proton cross-relaxation rates were determined in isotropic D(2)O solution. In the lyotropic liquid-crystalline medium consisting of ditetradecylphosphatidylcholine, dihexylphosphatidylcholine, N-cetyl-N,N,N-trimethylammonium bromide, and D(2)O, (1)H, (1)H and one-bond (1)H, (13)C residual dipolar couplings (RDCs), as well as relative sign information on homonuclear RDCs, were determined for the pentasaccharide. Molecular dynamics simulations with explicit water were carried out from which the internal isomerization relaxation time constant, τ(N), was calculated for transitions at the ψ torsion angle of the β-(1 → 3) linkage to the lactosyl group in LNF-1. Compared to the global reorientation time, τ(M), of ∼0.6 ns determined experimentally in D(2)O solution, the time constant for the isomerization relaxation process, τ(N(scaled)), is about one-third as large. The NMR parameters derived from the isotropic solution show very good agreement with those calculated from the MD simulations. The only notable difference occurs at the reducing end, which should be more flexible than observed by the molecular simulation, a conclusion in complete agreement with previous (13)C NMR relaxation data. A hydrogen-bond analysis of the MD simulation revealed that inter-residue hydrogen bonds on the order of ∼30% were present across the glycosidic linkages to sugar ring oxygens. This finding highlights that intramolecular hydrogen bonds might be important in preserving well-defined structures in otherwise flexible molecules. An analysis including generalized order parameters obtained from nuclear spin relaxation experiments was performed and successfully shown to limit the conformational space accessible to the molecule when the number of experimental data are too scarce for a complete conformational analysis.

Plant O-hydroxyproline arabinogalactans are composed of repeating trigalactosyl subunits with short bifurcated side chains

Classical arabinogalactan proteins partially defined by type II O-Hyp-linked arabinogalactans (Hyp-AGs) are structural components of the plant extracellular matrix. Recently we described the structure of a small Hyp-AG putatively based on repetitive trigalactosyl subunits and suggested that AGs are less complex and varied than generally supposed. Here we describe three additional AGs with similar subunits. The Hyp-AGs were isolated from two different arabinogalactan protein fusion glycoproteins expressed in tobacco cells; that is, a 22-residue Hyp-AG and a 20-residue Hyp-AG, both isolated from interferon alpha2b-(Ser-Hyp)(20), and a 14-residue Hyp-AG isolated from (Ala-Hyp)(51)-green fluorescent protein. We used NMR spectroscopy to establish the molecular structure of these Hyp-AGs, which share common features: (i) a galactan main chain composed of two 1-->3 beta-linked trigalactosyl blocks linked by a beta-1-->6 bond; (ii) bifurcated side chains with Ara, Rha, GlcUA, and a Gal 6-linked to Gal-1 and Gal-2 of the main-chain trigalactosyl repeats; (iii) a common side chain structure composed of up to six residues, the largest consisting of an alpha-L-Araf-(1-->5)-alpha-L-Araf-(1-->3)-alpha-L-Araf-(1-->3- unit and an alpha-L-Rhap-(1-->4)-beta-D-GlcUAp-(1-->6)-unit, both linked to Gal. The conformational ensemble obtained by using nuclear Overhauser effect data in structure calculations revealed a galactan main chain with a reverse turn involving the beta-1-->6 link between the trigalactosyl blocks, yielding a moderately compact structure stabilized by H-bonds.

Molecular modeling and NMR studies of benzyl substituted mannosyl trisaccharide binding to two mannose-specific lectins: Allium sativam agglutinin I and Concanavalin A

The interaction of trimannoside, α-benzyl 3, 6-di-O-(α-D-mannopyranosyl)-α-D-mannopyranoside, 1 with ASAI (Allium sativam agglutinin I, garlic lectin) was studied to reveal the conformational preferences of this ligand in bound-state and detailed binding mode at atomic level. The binding phenomenon was then compared with another well-known mannose-binding lectin, ConA (Concanavalin A). Structural studies of the ligand in free state were done using NMR spectroscopy and Molecular Dynamics simulations. It is found that the substituted-trimannoside can undergo conformational transitions in solution, with one major and one minor conformation per glycosidic linkage (α 1→3 and α 1→6). On the other hand in the bound-state only one of the two major conformations was significantly populated. The role of phenyl ring in the binding process was explored. An extended binding site was observed for the trimannoside in ASAI utilizing the aromatic substituent, which is not seen in ConA. Binding data from difference absorption spectroscopy supported this fact that the binding of benzyl-substituted ligand is tighter with ASAI than ConA. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 952-967, 2010.

GROMOS96 43a1 performance in predicting oligosaccharide conformational ensembles within glycoproteins

In previous work [Pol-Fachin, L.; Fernandes, C. L.; Verli, H.; Carbohydr. Res.2009, 344, 491-500], we had demonstrated that GROMOS96 43a1 force field and Löwdin HF/6-31G * *-derived atomic charges, adequately represent a glycoprotein's conformational ensemble in aqueous solutions, taking as the starting geometries NMR-determined structures. Based on such data, the present work intends to evaluate the use of the main solution conformations of isolated disaccharides, to build the carbohydrate moiety of glycoproteins, for which no previous experimental information is available. The observed results suggested that the entire glycoprotein scaffold appears unable to promote major modifications in the conformational behavior of glycosidic linkages. Additionally, when compared to energy contour plots, the results support the use of solution ensembles, to refine vacuum conformations of carbohydrate databases in the assembling of glycoproteins 3D structures. Finally, such approach is applied to build a full glycosylated model for COX-1 and COX-2 enzymes.

When sugars get wet. A comprehensive study of the behavior of water on the surface of oligosaccharides

In this article, we characterize the behavior of water on the surface of a diverse group of carbohydrates and attempt to determine the role of saccharide size, linkage, and branching as well as secondary structure on the dynamics and structure of water at the surface. In order to better understand the similarities and differences in the behavior of the solvent on the carbohydrate surface, we explore residence times, rotational correlation functions, local solvent occupancy numbers, and diffusivities. We find that due to the differences in secondary structure water residence times are longer and translational and rotational dynamics are retarded when in contact with wide helices and branched sugars. In the case of extended helices and smaller oligosaccharides, water dynamics is faster and less hindered. This indicates that branching, the type of linkage between monomers, and the anomeric configuration all play a major role in determining the structure and dynamics of water on the surface of carbohydrates.

A tool for the prediction of structures of complex sugars

In two recent back to back articles(Xia et al., J Chem Theory Comput 3:1620-1628 and 1629-1643, 2007a, b) we have started to address the problem of complex oligosaccharide conformation and folding. The scheme previously presented was based on exhaustive searches in configuration space in conjunction with Nuclear Overhauser Effect (NOE) calculations and the use of a complex rotameric library that takes branching into account. NOEs are extremely useful for structural determination but only provide information about short range interactions and ordering. Instead, the measurement of residual dipolar couplings (RDC), yields information about molecular ordering or folding that is long range in nature. In this article we show the results obtained by incorporation RDC calculations into our prediction scheme. Using this new approach we are able to accurately predict the structure of six human milk sugars: LNF-1, LND-1, LNF-2, LNF-3, LNnT and LNT. Our exhaustive search in dihedral configuration space combined with RDC and NOE calculations allows for highly accurate structural predictions that, because of the non-ergodic nature of these molecules on a time scale compatible with molecular dynamics simulations, are extremely hard to obtain otherwise (Almond et al., Biochemistry 43:5853-5863, 2004). Molecular dynamics simulations in explicit solvent using as initial configurations the structures predicted by our algorithm show that the histo-blood group epitopes in these sugars are relatively rigid and that the whole family of oligosaccharides derives its conformational variability almost exclusively from their common linkage (beta-D: -GlcNAc-(1-->3)-beta-D: -Gal) which can exist in two distinct conformational states. A population analysis based on the conformational variability of this flexible glycosidic link indicates that the relative population of the two distinct states varies for different human milk oligosaccharides.