O-Methylation in Carbohydrates: An NMR and MD Simulation Study with Application to Methylcellulose

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

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 ¹H,¹H-NOESY and 1D ¹H,¹H-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 ¹³C-isotope labeled α-D-[1,2-¹³C₂]Manp-(1→4)-α-D-Manp-OMe was synthesized thereby facilitating conformational analysis based on ¹³C,¹³C 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.

A Molecular Dynamic Study of the Effects of Surface Partitioning on the OH Radical Interactions with Solutes in Multicomponent Aqueous Aerosols

The surface-bulk partitioning of small saccharide and amide molecules in aqueous droplets was investigated using molecular dynamics. The air-particle interface was modeled using a 80 Å cubic water box containing a series of organic molecules and surrounded by gaseous OH radicals. The properties of the organic solutes within the interface and the water bulk were examined at a molecular level using density profiles and radial pair distribution functions. Molecules containing only polar functional groups such as urea and glucose are found predominantly in the water bulk, forming an exclusion layer near the water surface. Substitution of a single polar group by an alkyl group in sugars and amides leads to the migration of the molecule toward the interface. Within the first 2 nm from the water surface, surface-active solutes lose their rotational freedom and adopt a preferred orientation with the alkyl group pointing toward the surface. The different packing within the interface leads to different solvation shell structures and enhanced interaction between the organic molecules and absorbed OH radicals. The simulations provide quantitative information about the dimension, composition, and organization of the air-water interface as well as about the nonreactive interaction of the OH radicals with the organic solutes. It suggests that increased concentrations, preferred orientations, and decreased solvation near the air-water surface may lead to differences in reactivities between surface-active and surface-inactive molecules. The results are important to explain how heterogeneous oxidation mechanisms and kinetics within interfaces may differ from those of the bulk.

Primary Structure of Glycans by NMR Spectroscopy

Glycans, carbohydrate molecules in the realm of biology, are present as biomedically important glycoconjugates and a characteristic aspect is that their structures in many instances are branched. In determining the primary structure of a glycan, the sugar components including the absolute configuration and ring form, anomeric configuration, linkage(s), sequence, and substituents should be elucidated. Solution state NMR spectroscopy offers a unique opportunity to resolve all these aspects at atomic resolution. During the last two decades, advancement of both NMR experiments and spectrometer hardware have made it possible to unravel carbohydrate structure more efficiently. These developments applicable to glycans include, inter alia, NMR experiments that reduce spectral overlap, use selective excitations, record tilted projections of multidimensional spectra, acquire spectra by multiple receivers, utilize polarization by fast-pulsing techniques, concatenate pulse-sequence modules to acquire several spectra in a single measurement, acquire pure shift correlated spectra devoid of scalar couplings, employ stable isotope labeling to efficiently obtain homo- and/or heteronuclear correlations, as well as those that rely on dipolar cross-correlated interactions for sequential information. Refined computer programs for NMR spin simulation and chemical shift prediction aid the structural elucidation of glycans, which are notorious for their limited spectral dispersion. Hardware developments include cryogenically cold probes and dynamic nuclear polarization techniques, both resulting in enhanced sensitivity as well as ultrahigh field NMR spectrometers with a ¹H NMR resonance frequency higher than 1 GHz, thus improving resolution of resonances. Taken together, the developments have made and will in the future make it possible to elucidate carbohydrate structure in great detail, thereby forming the basis for understanding of how glycans interact with other molecules.

Rationalizing the Stereoelectronic Influence of Interglycosidic Bond Conformations on the Reactivity of 1,4-O-Linked Disaccharide Donors

Disaccharide donors are key precursors in convergent glycan synthesis strategies. Unexpectedly, we observed that disaccharide thioglycosyl donors containing 1,4-O-linked α-glycosidic bonds are much more reactive than their β-analogues with the same protecting group pattern. Herein, we rationalized that such a difference in their reactivity is attributed to the conformation of the 1,4-O-interglycosidic bond which is controlled by anomeric and exo-anomeric effects. Moreover, the conformational preferences of these donors are dictated by the dihedral angles ϕ and ψ of their interglycosidic linkages and the torsional angle ω of their side chain along the C5-C6 bond. This fundamental research clarifies how the long-range stereoelectronic effects from the nonreducing end sugar can influence the reactivity of the leaving group at the reducing end and the behavior of disaccharide donors thereof.