On the influence of solvent on the stereoselectivity of glycosylation reactions.

Methodology development in carbohydrate chemistry entails the stereoselective formation of C-O bonds as a key step in the synthesis of oligo- and polysaccharides. The anomeric selectivity of a glycosylation reaction is affected by a multitude of parameters, such as the nature of the donor and acceptor, activator/promotor system, temperature and solvent. The influence of different solvents on the stereoselective outcome of glycosylation reactions employing thioglucopyranosides as glycosyl donors with a non-participating protecting group at position 2 has been studied. A large change in selectivity as a function of solvent was observed and a correlation between selectivity and the Kamlet-Taft solvent parameter π* was found. Furthermore, molecular modeling using density functional theory methodology was conducted to decipher the role of the solvent and possible reaction pathways were investigated.

Conformational flexibility of the disaccharide ß-L-Fucp-(1→4)-α-D-Glcp-OMe as deduced from NMR spectroscopy experiments and computer simulations.

Carbohydrates in biological systems are referred to as glycans and modification of their structures is a hallmark indicator of disease. Analysis of the three-dimensional structure forms the basis for further insight into how they function and comparison of crystal structure with solution-state conformation(s) is particularly relevant, which has been performed for the disaccharide ß-L-Fucp-(1→4)-α-D-Glcp-OMe. In water solution the conformational space at the glycosidic linkage between the two sugar residues is identified from molecular dynamics (MD) simulations as having a low-energy exo-syn conformation, deviating somewhat from the solid-state conformation, and two anti-conformational states, i.e., anti-ϕ and anti-ψ, indicating flexibility at the glycosidic linkage. NMR data were obtained from 1D 1H,1H-NOESY and STEP-NOESY experiments, measurement of transglycosidic 3JCH coupling constants and NMR spin-simulation. The free energy profile of the ω torsion angle computed from MD simulation was in excellent agreement with the rotamer distribution from NMR experiment being for gt:gg:tg 38 : 53 : 9, respectively, with a proposed inter-residue O5'⋯HO6 hydrogen bond being predominant in the gg rotamer. Quantum mechanics methodology was used to calculate transglycosidic NMR 3JCH coupling constants, averaged over a conformational ensemble of structures representing various rotamers of exocyclic groups, in good to excellent agreement with Karplus-type relationships previously developed. Furthermore, 1H and 13C NMR chemical shifts were calculated using the same methodology and were found to be in excellent agreement with experimental data.

Glycosidic α-linked mannopyranose disaccharides: an NMR spectroscopy and molecular dynamics simulation study employing additive and Drude polarizable force fields.

D-Mannose is a structural component in N-linked glycoproteins from viruses and mammals as well as in polysaccharides from fungi and bacteria. Structural components often consist of D-Manp residues joined via α-(1→2)-, α-(1→3)-, α-(1→4)- or α-(1→6)-linkages. As models for these oligo- and polysaccharides, a series of mannose-containing disaccharides have been investigated with respect to conformation and dynamics. Translational diffusion NMR experiments were performed to deduce rotational correlation times for the molecules, 1D 1H,1H-NOESY and 1D 1H,1H-T-ROESY NMR experiments were carried out to obtain inter-residue proton-proton distances and one-dimensional long-range and 2D J-HMBC experiments were acquired to gain information about conformationally dependent heteronuclear coupling constants across glycosidic linkages. To attain further spectroscopic data, the doubly 13C-isotope labeled α-D-[1,2-13C2]Manp-(1→4)-α-D-Manp-OMe was synthesized thereby facilitating conformational analysis based on 13C,13C coupling constants as interpreted by Karplus-type relationships. Molecular dynamics simulations were carried out for the disaccharides with explicit water as solvent using the additive CHARMM36 and Drude polarizable force fields for carbohydrates, where the latter showed broader population distributions. Both simulations sampled conformational space in such a way that inter-glycosidic proton-proton distances were very well described whereas in some cases deviations were observed between calculated conformationally dependent NMR scalar coupling constants and those determined from experiment, with closely similar root-mean-square differences for the two force fields. However, analyses of dipole moments and radial distribution functions with water of the hydroxyl groups indicate differences in the underlying physical forces dictating the wider conformational sampling with the Drude polarizable versus additive C36 force field and indicate the improved utility of the Drude polarizable model in investigating complex carbohydrates.

Complete 1H and 13C NMR chemical shift assignments of mono-to tetrasaccharides as basis for NMR chemical shift predictions of oligo- and polysaccharides using the computer program CASPER.

Carbohydrate structure can be elucidated or confirmed by using NMR spectroscopy as the prime technique. Prediction of 1H and 13C NMR chemical shifts by computational approaches makes this assignment process more efficient and the program CASPER can perform this task rapidly. It does so by relying on chemical shift data of mono-, di-, and trisaccharides. In order to improve accuracy and quality of these predictions we have assigned 1H and 13C NMR chemical shifts of 30 monosaccharides, 17 disaccharides, 10 trisaccharides and one tetrasaccharide; in total 58 compounds. Due to different rotamers, ring forms, α- and ß-anomeric forms and pD conditions this resulted in 74 1H and 13C NMR chemical shift data sets, all of which were refined using total line-shape analysis for the 1H resonances in order to obtain accurate chemical shifts. Subsequent NMR chemical shift predictions for three sialic acid-containing oligosaccharides, viz., GD1a, a disialyl-LNnT hexasaccharide and a polysialic acid-lactose decasaccharide, and NMR-based structural elucidations of two O-antigen polysaccharides from E. coli O174 were performed by the CASPER program (http://www.casper.organ.su.se/casper/) resulting in very good to excellent agreement between experimental and predicted data thereby demonstrating its utility for carbohydrate compounds that have been chemically or enzymatically synthesized, structurally modified or isolated from nature.

Differential Dynamics at Glycosidic Linkages of an Oligosaccharide as Revealed by 13C NMR Spin Relaxation and Stochastic Modeling.

Among biomolecules, carbohydrates are unique in that not only can linkages be formed through different positions, but the structures may also be branched. The trisaccharide ß-d-Glcp-(1→3)[ß-d-Glcp-(1→2)]-α-d-Manp-OMe represents a model of a branched vicinally disubstituted structure. A 13C site-specific isotopologue, with labeling in each of the two terminal glucosyl residues, enabled the acquisition of high-quality 13C NMR relaxation parameters, T1 and T2, and heteronuclear NOE, with standard deviations of ≤0.5%. For interpretation of the experimental NMR data, a diffusive chain model was used, in which the dynamics of the glycosidic linkages is coupled to the global reorientation motion of the trisaccharide. Brownian dynamics simulations relying on the potential of mean force at the glycosidic linkages were employed to evaluate spectral densities of the spin probes. Calculated NMR relaxation parameters showed a very good agreement with experimental data, deviating <3%. The resulting dynamics are described by correlation times of 196 and 174 ps for the ß-(1→2)- and ß-(1→3)-linked glucosyl residues, respectively, i.e., different and linkage dependent. Notably, the devised computational protocol was performed without any fitting of parameters.

Structure-Reactivity Relationships of Conformationally Armed Disaccharide Donors and Their Use in the Synthesis of a Hexasaccharide Related to the Capsular Polysaccharide from Streptococcus pneumoniae Type 37.

To advance the field of glycobiology, efficient synthesis methods of oligosaccharides and glycoconjugates are a requisite. In glycosylation reactions using superarmed donors, both selectivity and reactivity issues must be considered, and we herein investigate these aspects for differently protected ß-linked 2-O-glycosylated glucosyl donors carrying bulky tert-butyldimethylsilyl groups to different extents. The acceptors in reactions being secondary alcohols presents a challenging situation with respect to steric crowding. Conformational pyranose ring equilibria of the superarmed disaccharide donors with axial-rich substituents contained skew and boat conformations, and three-state models were generally assumed. With NIS/TfOH as the promotor, 2,6-di-tert-butyl-4-methylpyridine as the base, and a dichloromethane/toluene solvent mixture, ethyl 1-thio-ß-d-glucosyl disaccharide donors having 6-O-benzyl group(s) besides tert-butyldimethylsilyl groups were efficiently coupled at -40 °C to the hydroxyl group at position 3 of glucopyranosyl acceptors to form ß-(1 → 2),ß-(1 → 3)-linked trisaccharides, isolated in excellent 95% yield. The more axial-rich donors in skew and boat conformations are thus preorganized closer to the assumed transition state in these glycosylation reactions. The developed methodology was subsequently applied in the synthesis of a multibranched hexasaccharide related to the capsular polysaccharide from Streptococcus pneumoniae type 37, which consists of a ß-(1 → 3)-linked backbone and a ß-(1 → 2)-linked side chain of d-glucosyl residues in disaccharide repeating units.

A molecular dynamics study of the effect of glycosidic linkage type in the hemicellulose backbone on the molecular chain flexibility.

The macromolecular conformation of the constituent polysaccharides in lignocellulosic biomass influences their supramolecular interactions, and therefore their function in plants and their performance in technical products. The flexibility of glycosidic linkages from the backbone of hemicelluloses was studied by evaluating the conformational freedom of the φ and ψ dihedral angles using molecular dynamic simulations, additionally selected molecules were correlated with experimental data by nuclear magnetic resonance spectroscopy. Three types of ß-(1→4) glycosidic linkages involving the monosaccharides (Glcp, Xylp and Manp) present in the backbone of hemicelluloses were defined. Different di- and tetrasaccharides with combinations of such sugar monomers from hemicelluloses were simulated, and free energy maps of the φ - ψ space and hydrogen-bonding patterns were obtained. The glycosidic linkage between Glc-Glc or Glc-Man (C-type) was the stiffest with mainly one probable conformation; the linkage from Man-Man or Man-Glc (M-type) was similar but with an increased probability for an alternative conformation making it more flexible, and the linkage between two Xyl-units (X-type) was the most flexible with two almost equally populated conformations. Glycosidic linkages of the same type showed essentially the same conformational space in both disaccharides and in the central region of tetrasaccharides. Different probabilities of glycosidic linkage conformations in the backbone of hemicelluloses can be directly estimated from the free energy maps, which to a large degree affect the overall macromolecular conformations of these polymers. The information gained contributes to an increased understanding of the function of hemicelluloses both in the cell wall and in technical products.

Delineating the conformational flexibility of trisaccharides from NMR spectroscopy experiments and computer simulations.

The conformation of saccharides in solution is challenging to characterize in the context of a single well-defined three-dimensional structure. Instead, they are better represented by an ensemble of conformations associated with their structural diversity and flexibility. In this study, we delineate the conformational heterogeneity of five trisaccharides via a combination of experimental and computational techniques. Experimental NMR measurements target conformationally sensitive parameters, including J couplings and effective distances around the glycosidic linkages, while the computational simulations apply the well-calibrated additive CHARMM carbohydrate force field in combination with efficient enhanced sampling molecular dynamics simulation methods. Analysis of conformational heterogeneity is performed based on sampling of discreet states as defined by dihedral angles, on root-mean-square differences of Cartesian coordinates and on the extent of volume sampled. Conformational clustering, based on the glycosidic linkage dihedral angles, shows that accounting for the full range of sampled conformations is required to reproduce the experimental data, emphasizing the utility of the molecular simulations in obtaining an atomic detailed description of the conformational properties of the saccharides. Results show the presence of differential conformational preferences as a function of primary sequence and glycosidic linkage types. Significant differences in conformational ensembles associated with the anomeric configuration of a single glycosidic linkage reinforce the impact of such changes on the conformational properties of carbohydrates. The present structural insights of the studied trisaccharides represent a foundation for understanding the range of conformations adopted in larger oligosaccharides and how these molecules encode their conformational heterogeneity into the monosaccharide sequence.

Temperature Dependence of Hydroxymethyl Group Rotamer Populations in Cellooligomers.

Empirical force fields for computer simulations of carbohydrates are often implicitly assumed to be valid also at temperatures different from room temperature for which they were optimized. Herein, the temperature dependence of the hydroxymethyl group rotamer populations in short oligosaccharides is investigated using molecular dynamics simulations and NMR spectroscopy. Two oligosaccharides, viz., methyl ß-cellobioside and ß-cellotetraose were simulated using three different carbohydrate force fields (CHARMM C35, GLYCAM06, and GROMOS 56Acarbo) in combination with different water models (SPC, SPC/E, and TIP3P) using replica exchange molecular dynamics simulations. For comparison, hydroxymethyl group rotamer populations were investigated for methyl ß-cellobioside and cellopentaose based on measured NMR (3)JH5,H6 coupling constants, in the latter case by using a chemical shift selective NMR-filter. Molecular dynamics simulations in combination with NMR spectroscopy show that the temperature dependence of the hydroxymethyl rotamer population in these short cellooligomers, in the range 263-344 K, generally becomes exaggerated in simulations when compared to experimental data, but also that it is dependent on simulation conditions, and most notably properties of the water model.