Population distribution of flexible molecules from maximum entropy analysis using different priors as background information: application to the Φ, Ψ-conformational space of the α-(1-->2)-linked mannose disaccharide present in N- and O-linked glycoproteins

Nonreducing Sugar Scaffold Enables the Development of Immunomodulatory TLR4-specific LPS Mimetics with Picomolar Potency

Innate immune defense mechanisms against infection and cancer encompass the modulation of pattern recognition receptor (PRR)-mediated inflammation, including upregulation of various transcription factors and the activation of pro-inflammatory pathways important for immune surveillance. Dysfunction of PRRs-mediated signaling has been implicated in cancer and autoimmune diseases, while the overactivation of PRRs-driven responses during infection can lead to devastating consequences such as acute lung injury or sepsis. We used crystal structure-based design to develop immunomodulatory lipopolysaccharide (LPS) mimetics targeting one of the ubiquitous PRRs, Toll-like Receptor 4 (TLR4). Taking advantage of an exo-anomeric conformation and specific molecular shape of synthetic nonreducing β,β-diglucosamine, which was investigated by NMR, we developed two sets of lipid A mimicking glycolipids capable of either potently activating innate immune responses or inhibiting pro-inflammatory signaling. Stereoselective 1,1'-glycosylation towards fully orthogonally protected nonreducing GlcNβ(1↔1')βGlcN followed by stepwise assembly of differently functionalised phosphorylated glycolipids provided biologically active molecules that were evaluated for their ability to trigger or to inhibit cellular innate immune responses. Two LPS mimetics, identified as potent TLR4-specific inducers of the intracellular signaling pathways, serve as vaccine adjuvant- and immunotherapy candidates, while anionic glycolipids with TLR4-inhibitory potential hold therapeutic promise for the management of acute or chronic inflammation.

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.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

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.

Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Capsular Polysaccharides

Capsule is one of the common virulence factors in Gram-negative bacteria protecting pathogens from host defenses and consists of long-chain capsular polysaccharides (CPS) anchored in the outer membrane (OM). Elucidating structural properties of CPS is important to understand its biological functions as well as the OM properties. However, the outer leaflet of the OM in current simulation studies is represented exclusively by LPS due to the complexity and diversity of CPS. In this work, representative Escherichia coli CPS, K_LPS (a lipid A-linked form) and K_PG (a phosphatidylglycerol-linked form), are modeled and incorporated into various symmetric bilayers with co-existing LPS in different ratios. All-atom molecular dynamics simulations of these systems have been conducted to characterize various bilayer properties. Incorporation of K_LPS makes the acyl chains of LPS more rigid and ordered, while incorporation of K_PG makes them less ordered and flexible. These results are consistent with the calculated area per lipid (APL) of LPS, in which the APL of LPS becomes smaller when K_LPS is incorporated, whereas it gets larger when K_PG is included. Torsional analysis reveals that the influence of the CPS presence on the conformational distributions of the glycosidic linkages of LPS is small, and minor differences are also detected for the inner and outer regions of the CPS. Combined with previously modeled enterobacterial common antigens (ECAs) in the form of mixed bilayers, this work provides more realistic OM models as well as the basis for characterization of interactions between the OM and OM proteins.

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.

Elucidating the enhanced binding affinity of a double mutant SP-D with trimannose on the influenza A virus using molecular dynamics

•

The Asp325Ala mutation in SP-D promotes a trimannose conformational change to a more stable state.

•

The Arg343Val mutation in SP-D reduces its interaction with Glu333 to increase the binding affinity with trimannose.

•

The Arg343Val mutation contributes more to the increase of SP-D’s binding affinity with trimannose than Asp325Ala.

Surfactant protein D (SP-D) is an essential component of the human pulmonary surfactant system, which is crucial in the innate immune response against glycan-containing pathogens, including Influenza A viruses (IAV) and SARS-CoV-2. Previous studies have shown that wild-type (WT) SP-D can bind IAV but exhibits poor antiviral activities. However, a double mutant (DM) SP-D consisting of two point mutations (Asp325Ala and Arg343Val) inhibits IAV more potently. Presently, the structural mechanisms behind the point mutations’ effects on SP-D’s binding affinity with viral surface glycans are not fully understood. Here we use microsecond-scale, full-atomistic molecular dynamics (MD) simulations to understand the molecular mechanism of mutation-induced SP-D’s higher antiviral activity. We find that the Asp325Ala mutation promotes a trimannose conformational change to a more stable state. Arg343Val increases the binding with trimannose by increasing the hydrogen bonding interaction with Glu333. Free energy perturbation (FEP) binding free energy calculations indicate that the Arg343Val mutation contributes more to the increase of SP-D’s binding affinity with trimannose than Asp325Ala. This study provides a molecular-level exploration of how the two mutations increase SP-D binding affinity with trimannose, which is vital for further developing preventative strategies for related diseases.

Use of Raman and Raman optical activity to extract atomistic details of saccharides in aqueous solution

Sugars are crucial components in biosystems and industrial applications. In aqueous environments, the natural state of short saccharides or charged glycosaminoglycans is floating and wiggling in solution. Therefore, tools to characterize their structure in a native aqueous environment are crucial but not always available. Here, we show that a combination of Raman/ROA and, on occasions, NMR experiments with Molecular Dynamics (MD) and Quantum Mechanics (QM) is a viable method to gain insights into structural features of sugars in solutions. Combining these methods provides information about accessible ring puckering conformers and their proportions. It also provides information about the conformation of the linkage between the sugar monomers, i.e., glycosidic bonds, allowing for identifying significantly accessible conformers and their relative abundance. For mixtures of sugar moieties, this method enables the deconvolution of the Raman/ROA spectra to find the actual amounts of its molecular constituents, serving as an effective analytical technique. For example, it allows calculating anomeric ratios for reducing sugars and analyzing more complex sugar mixtures to elucidate their real content. Altogether, we show that combining Raman/ROA spectroscopies with simulations is a versatile method applicable to saccharides. It allows for accessing many features with precision comparable to other methods routinely used for this task, making it a viable alternative. Furthermore, we prove that the proposed technique can scale up by studying the complicated raffinose trisaccharide, and therefore, we expect its wide adoption to characterize sugar structural features in solution.

Preferred conformations of lipooligosaccharides and oligosaccharides of Moraxella catarrhalis

Moraxella catarrhalis (M. catarrhalis) is a pathogenic gram-negative bacterium that causes otitis media and sinusitis in children. Three major serotypes A, B and C are identified to account for approximately 95% of the clinical isolates. Understanding the conformational properties of different serotypes of M. catarrhalis provides insights into antigenic determinants. In this work, all-atom molecular dynamics simulations were conducted for M. catarrhalis lipooligosaccharide (LOS) bilayer systems and oligosaccharides (OS) in water solution to investigate the conformational similarities and differences of three serotypes. For up to 10 neutral monosaccharides in the core part, the conformational ensembles described by the pair-wise root mean square deviation distributions are similar among the three serotypes of either the LOS or OS. At the central β-($1\to4$)-linkage, anti-$\psi$ conformation in conjunction with the gauche-gauche (g-) conformation of the central trisubstituted glucosyl residue is observed as the dominant conformation to sustain the structural characteristics of M. catarrhalis three types, which is further supported by calculated transglycosidic ${}^3{J}_{C,H}\Big({\psi}_H\Big)$ of serotype A in comparison to experimental data. Interestingly, the conformational variability of three serotypes is more restricted for the OS in water solution than that in the LOS bilayer systems. The LOS-LOS interactions in the bilayer systems are responsible for the increased conformational diversity despite of tight packing. Solvent-accessible surface area analysis suggests that a trisaccharide attached to the β-($1\to 6$)-linked sugar in all three serotypes of LOS could be the common epitope and have the possibility to interact with antibodies.

A detailed picture of a protein-carbohydrate hydrogen-bonding network revealed by NMR and MD simulations

Cyanovirin-N (CV-N) is a cyanobacterial lectin with antiviral activity towards HIV and several other viruses. Here, we identify mannoside hydroxyl protons that are hydrogen bonded to the protein backbone of the CV-N domain B binding site, using NMR spectroscopy. For the two carbohydrate ligands Manα(1→2)ManαOMe and Manα(1→2) Manα(1→6)ManαOMe five hydroxyl protons are involved in hydrogen-bonding networks. Comparison with previous crystallographic results revealed that four of these hydroxyl protons donate hydrogen bonds to protein backbone carbonyl oxygens in solution and in the crystal. Hydrogen bonds were not detected between the side chains of Glu41 and Arg76 with sugar hydroxyls, as previously proposed for CV-N binding of mannosides. Molecular dynamics simulations of the CV-N/Manα(1→2)Manα(1→6)ManαOMe complex confirmed the NMR-determined hydrogen-bonding network. Detailed characterization of CV-N/mannoside complexes provides a better understanding of lectin-carbohydrate interactions and opens up to the use of CV-N and similar lectins as antiviral agents.

Inter-residual Hydrogen Bonding in Carbohydrates Unraveled by NMR Spectroscopy and Molecular Dynamics Simulations

Carbohydrates, also known as glycans in biological systems, are omnipresent in nature where they as glycoconjugates occur as oligo- and polysaccharides linked to lipids and proteins. Their three-dimensional structure is defined by two or three torsion angles at each glycosidic linkage. In addition, transglycosidic hydrogen bonding between sugar residues may be important. Herein we investigate the presence of these inter-residue interactions by NMR spectroscopy in D₂ O/[D₆ ]DMSO (70:30) or D₂ O and by molecular dynamics (MD) simulations with explicit water as solvent for disaccharides with structural elements α-d-Manp-(1→2)-d-Manp, β-d-GlcpNAc-(1→2)-d-Manp, and α-d-Glcp-(1→4)-β-d-Glcp, all of which have been suggested to exhibit inter-residue hydrogen bonding. For the disaccharide β-d-GlcpNAc-(1→2)-β-d-Manp-OMe, the large extent of O5'⋅⋅⋅HO3 hydrogen bonding as seen from the MD simulation is implicitly supported by the ¹ H NMR chemical shift and ³ J_HO3,H3 value of the hydroxy proton. In the case of α-d-Glcp-(1→4)-β-d-Glcp-OMe, the existence of a transglycosidic hydrogen bond O2'⋅⋅⋅HO3 was proven by the presence of a cross-peak in ¹ H,¹³ C HSQC-TOCSY experiments as a result of direct TOCSY transfer between HO3 of the reducing end residue and H2' (detected at C2') of the terminal residue. The occurrence of inter-residue hydrogen bonding, albeit transient, is judged important for the stabilization of three-dimensional structures, which may be essential in maintaining a conformational state for carbohydrate-protein interactions of glycans to take place in biologically important environments.

Elucidation of the O-antigen structure of Escherichia coli O63

The structure of the O-antigen polysaccharide (PS) from the Shiga-toxin producing Escherichia coli O63 has been elucidated using a combination of bioinformatics, component analyses and NMR spectroscopy. The O-antigen is comprised of tetrasaccharide repeating units with the following structure: →2)-β-d-Quip3N(d-allo-ThrAc)-(1→2)-β-d-Ribf-(1→4)-β-d-Galp-(1→3)-α-d-GlcpNAc-(1→ in which the N-acetylated d-allo-threonine is amide-linked to position 3 of the 3-amino-3-deoxy-d-Quip sugar residue. The presence of a predicted flippase and polymerase encoded in the O63 gene cluster is consistent with the Wzx/Wzy biosynthetic pathway and consequently the biological repeating unit has likely an N-acetyl-d-glucosamine residue at its reducing end. A bioinformatics approach based on predictive glycosyltransferase function present in ECODAB (E. coli O-antigen database) suggested the structural element β-d-Galp-(1→3)-d-GlcpNAc in the O-antigen. Notably, multiple gene sequence alignment of fdtA and qdtA from E. coli to that in E. coli O63 resulted in discrimination between the two, confirmation of the latter in E. coli O63, and consequently, together with qdtB, biosynthesis of dTDP-d-Quip3N. The E. coli O63 O-antigen polysaccharide differs in two aspects from that of E. coli O114 where the latter carries instead an l-serine residue, and the glycosidic linkage positions to and from the Quip3N residue are both changed. The structural characterization of the O63 antigen repeat supports the predicted functional assignment of the O-antigen cluster genes.

Solution Structure of Mannobioses Unravelled by Means of Raman Optical Activity

Structural analysis of carbohydrates is a complicated endeavour, due to the complexity and diversity of the samples at hand. Herein, we apply a combined computational and experimental approach, employing molecular dynamics (MD) and density functional theory (DFT) calculations together with NMR and Raman optical activity (ROA) measurements, in the structural study of three mannobiose disaccharides, consisting of two mannoses with varying glycosidic linkages. The disaccharide structures make up the scaffold of high mannose glycans and are therefore important targets for structural analysis. Based on the MD population analysis and NMR, the major conformers of each mannobiose were identified and used as input for DFT analysis. By systematically varying the solvent models used to describe water interacting with the molecules and applying overlap integral analysis to the resulting calculational ROA spectra, we found that a full quantum mechanical/molecular mechanical approach is required for an optimal calculation of the ROA parameters. Subsequent normal mode analysis of the predicted vibrational modes was attempted in order to identify possible marker bands for glycosidic linkages. However, the normal mode vibrations of the mannobioses are completely delocalised, presumably due to conformational flexibility in these compounds, rendering the identification of isolated marker bands unfeasible.

Synthesis and O-Glycosidic Linkage Conformational Analysis of 13C-Labeled Oligosaccharide Fragments of an Antifreeze Glycolipid

NMR studies of two ¹³C-labeled disaccharides and a tetrasaccharide were undertaken that comprise the backbone of a novel thermal hysteresis glycolipid containing a linear glycan sequence of alternating [βXyl p-(1→4)-βMan p-(1→4)] _n dimers. Experimental trans-glycoside NMR J-couplings, parameterized equations obtained from density functional theory (DFT) calculations, and an in-house circular statistics package ( MA'AT) were used to derive conformational models of linkage torsion angles ϕ and ψ in solution, which were compared to those obtained from molecular dynamics simulations. Modeling using different probability distribution functions showed that MA'AT models of ϕ in βMan(1→4)βXyl and βXyl(1→4)βMan linkages are very similar in the disaccharide building blocks, whereas MA'AT models of ψ differ. This pattern is conserved in the tetrasaccharide, showing that linkage context does not influence linkage geometry in this linear system. Good agreement was observed between the MA'AT and MD models of ψ with respect to mean values and circular standard deviations. Significant differences were observed for ϕ, indicating that revision of the force-field employed by GLYCAM is probably needed. Incorporation of the experimental models of ϕ and ψ into the backbone of an octasaccharide fragment leads to a helical amphipathic topography that may affect the thermal hysteresis properties of the glycolipid.

Use of Circular Statistics To Model αMan-(1→2)-αMan and αMan-(1→3)-α/βMan O-Glycosidic Linkage Conformation in 13C-Labeled Disaccharides and High-Mannose Oligosaccharides

A new experimental method, MA' AT analysis, has been applied to investigate the conformational properties of O-glycosidic linkages in several biologically important mannose-containing di- and oligosaccharides. Methyl α-d-mannopyranosyl-(1→2)-α-d-mannopyranoside (2), methyl α-d-mannopyranosyl-(1→3)-α-d-mannopyranoside (3), and methyl α-d-mannopyranosyl-(1→3)-β-d-mannopyranoside (4) were prepared with selective ¹³C-enrichment to enable the measurement of NMR scalar couplings across their internal O-glycosidic linkages. Density functional theory (DFT) was used to parameterize equations for J_CH and J_CC values in 2-4 that are sensitive to phi (ϕ) and psi (ψ). The experimental J-couplings and parameterized equations were treated using a circular statistics algorithm encoded in the MA' AT program. Conformations about ϕ and ψ treated using single-state von Mises models gave excellent fits to the ensembles of redundant J-couplings. Mean values and circular standard deviations (CSDs) for each linkage torsion angle ϕ (CSD) and ψ (CSD) in 2, -29° (25°) and 20° (22°); in 3, -36° (36°) and 8° (27°); in 4, -37° (34°) and 10° (26°); ϕ = H1'-C1'-O1'-CX and ψ = C1'-O1'-CX-HX (CX = aglycone carbon) were compared to histograms obtained from 1 μs aqueous molecular dynamics (MD) simulations and X-ray database statistical analysis. MA' AT-derived models of ψ were in very good agreement with the MD and X-ray data, but not those of ϕ, suggesting a need for force field revision. The effect of structural context on linkage conformation was also investigated in four selectively ¹³C-labeled homomannose tri- and tetrasaccharides using the MA' AT method. In the cases examined, context effects were found to be small.

Bisecting GlcNAc restricts conformations of branches in model N-glycans with GlcNAc termini

Bisected N-glycans play significant roles in tumor migration and Alzheimer's disease through modulating the action and localization of their carrier proteins. Such biological functions are often discussed in terms of the conformation of the attached N-glycans with or without bisecting GlcNAc. To obtain insights into the effects of bisecting GlcNAc on glycan conformation, a systematic NMR structural analysis was performed on two pairs of synthetic N-glycans, with and without bisecting GlcNAc. The analysis reveals that terminal GlcNAcs and bisecting GlcNAc cooperate to restrict the conformations of both the α1-3 and α1-6 branches of N-glycans. ¹H and ¹³C chemical shift comparisons suggest that bisecting GlcNAc directly modulates local conformation. Unique NOE correlations between core-mannose and the α1-3 branch mannose as well as the ³J_C-H constant of the glycoside linkage indicate that bisecting GlcNAc restricts the conformation of the 1-3 branch. The angles of the glycosidic bonds between core-mannose and α1-6 branch mannose derived from ³J_C-H and ³J_H-H coupling constants show that terminal GlcNAcs restrict the distribution of the ψ angle to 180° and the bisecting GlcNAc increases the distribution of the ω angle +60° in the presence of terminal GlcNAcs. It is feasible that restriction of branch conformations by bisecting GlcNAc has important consequences for protein-glycan interplay and following biological events.

Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution

Molecular dynamics simulations of native α-, β-, and γ-cyclodextrin in aqueous solution have been conducted with the goal to investigate the performance of the CHARMM36 force field, the AMBER-compatible q4md-CD force field, and five variants of the GROMOS force field. The properties analyzed are structural parameters derived from X-ray diffraction and NMR experiments as well as hydrogen bonds and hydration patterns, including hydration free enthalpies. Recent revisions of the torsional-angle parameters for carbohydrate systems within the GROMOS family of force fields lead to a significant improvement of the agreement between simulated and experimental NMR data. Therefore, we recommend using the variant 53A6_GLYC instead of 53A6 and 56A6_{CARBO_R} or 2016H66 instead of 56A6_CARBO to simulate cyclodextrins in solution. The CHARMM36 and q4md-CD force fields show a similar performance as the three recommended GROMOS parameter sets. A significant difference is the more flexible nature of the cyclodextrins modeled with the CHARMM36 and q4md-CD force fields compared to the three recommended GROMOS parameter sets.

Modeling of Oligosaccharides within Glycoproteins from Free-Energy Landscapes

In spite of the abundance of glycoproteins in biological processes, relatively little three-dimensional structural data is available for glycan structures. Here, we study the structure and flexibility of the vast majority of mammalian oligosaccharides appearing in N- and O-glycosylated proteins using a bottom up approach. We report the conformational free-energy landscapes of all relevant glycosidic linkages as obtained from local elevation simulations and subsequent umbrella sampling. To the best of our knowledge, this represents the first complete conformational library for the construction of N- and O-glycan structures. Next, we systematically study the effect of neighboring residues, by extensively simulating all relevant trisaccharides and one tetrasaccharide. This allows for an unprecedented comparison of disaccharide linkages in large oligosaccharides. With a small number of exceptions, the conformational preferences in the larger structures are very similar as in the disaccharides. This, finally, allows us to suggest several efficient approaches to construct complete N- and O-glycans on glycoproteins, as exemplified on two relevant examples.

Exploiting Uniformly 13C-Labeled Carbohydrates for Probing Carbohydrate-Protein Interactions by NMR Spectroscopy

NMR of a uniformly ¹³C-labeled carbohydrate was used to elucidate the atomic details of a sugar-protein complex. The structure of the ¹³C-labeled Manα(1-2)Manα(1-2)ManαOMe trisaccharide ligand, when bound to cyanovirin-N (CV-N), was characterized and revealed that in the complex the glycosidic linkage torsion angles between the two reducing-end mannoses are different from the free trisaccharide. Distances within the carbohydrate were employed for conformational analysis, and NOE-based distance mapping between sugar and protein revealed that Manα(1-2)Manα(1-2)ManαOMe is bound more intimately with its two reducing-end mannoses into the domain A binding site of CV-N than with the nonreducing end unit. Taking advantage of the ¹³C spectral dispersion of ¹³C-labeled carbohydrates in isotope-filtered experiments is a versatile means for a simultaneous mapping of the binding interactions on both, the carbohydrate and the protein.

Flexibility at a glycosidic linkage revealed by molecular dynamics, stochastic modeling, and (13)C NMR spin relaxation: conformational preferences of α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe in water and dimethyl sulfoxide solutions

The monosaccharide L-rhamnose is common in bacterial polysaccharides and the disaccharide α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe represents a structural model for a part of Shigella flexneri O-antigen polysaccharides. Utilization of [1'-(13)C]-site-specific labeling in the anomeric position at the glycosidic linkage between the two sugar residues facilitated the determination of transglycosidic NMR (3)JCH and (3)JCC coupling constants. Based on these spin-spin couplings the major state and the conformational distribution could be determined with respect to the ψ torsion angle, which changed between water and dimethyl sulfoxide (DMSO) as solvents, a finding mirrored by molecular dynamics (MD) simulations with explicit solvent molecules. The (13)C NMR spin relaxation parameters T1, T2, and heteronuclear NOE of the probe were measured for the disaccharide in DMSO-d6 at two magnetic field strengths, with standard deviations ≤1%. The combination of MD simulation and a stochastic description based on the diffusive chain model resulted in excellent agreement between calculated and experimentally observed (13)C relaxation parameters, with an average error of <2%. The coupling between the global reorientation of the molecule and the local motion of the spin probe is deemed essential if reproduction of NMR relaxation parameters should succeed, since decoupling of the two modes of motion results in significantly worse agreement. Calculation of (13)C relaxation parameters based on the correlation functions obtained directly from the MD simulation of the solute molecule in DMSO as solvent showed satisfactory agreement with errors on the order of 10% or less.

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.

Molecular dynamics simulations and NMR spectroscopy studies of trehalose-lipid bilayer systems

The disaccharide trehalose (TRH) strongly affects the physical properties of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH using NMR spectroscopy and molecular dynamics (MD) computer simulations. We compare dipolar couplings derived from DMPC/TRH trajectories with those determined (i) experimentally in TRH using conventional high-resolution NMR in a weakly ordered solvent (bicelles), and (ii) by solid-state NMR in multilamellar vesicles (MLV) formed by DMPC. Analysis of the experimental and MD-derived couplings in DMPC indicated that the force field used in the simulations reasonably well describes the experimental results with the exception for the glycerol fragment that exhibits significant deviations. The signs of dipolar couplings, not available from the experiments on highly ordered systems, were determined from the trajectory analysis. The crucial step in the analysis of residual dipolar couplings (RDCs) in TRH determined in a bicelle-environment was access to the conformational distributions derived from the MD trajectory. Furthermore, the conformational behavior of TRH, investigated by J-couplings, in the ordered and isotropic phases is essentially identical, indicating that the general assumptions in the analyses of RDCs are well founded.

NMR of glycans: shedding new light on old problems

The diversity in molecular arrangements and dynamics displayed by glycans renders traditional NMR strategies, employed for proteins and nucleic acids, insufficient. Because of the unique properties of glycans, structural studies often require the adoption of a different repertoire of tailor-made experiments and protocols. We present an account of recent developments in NMR techniques that will deepen our understanding of structure-function relations in glycans. We open with a survey and comparison of methods utilized to determine the structure of proteins, nucleic acids and carbohydrates. Next, we discuss the structural information obtained from traditional NMR techniques like chemical shifts, NOEs/ROEs, and coupling-constants, along with the limitations imposed by the unique intrinsic characteristics of glycan structure on these approaches: flexibility, range of conformers, signal overlap, and non-first-order scalar (strong) coupling. Novel experiments taking advantage of isotopic labeling are presented as an option for overcoming spectral overlap and raising sensitivity. Computational tools used to explore conformational averaging in conjunction with NMR parameters are described. In addition, recent developments in hydroxyl detection and hydrogen bond detection in protonated solvents, in contrast to traditional sample preparations in D2O for carbohydrates, further increase the tools available for both structure information and chemical shift assignments. We also include previously unpublished data in this context. Accurate determination of couplings in carbohydrates has been historically challenging due to the common presence of strong-couplings. We present new strategies proposed for dealing with their influence on NMR signals. We close with a discussion of residual dipolar couplings (RDCs) and the advantages of using (13)C isotope labeling that allows gathering one-bond (13)C-(13)C couplings with a recently improved constant-time COSY technique, in addition to the commonly measured (1)H-(13)C RDCs.

Conformational properties of α- or β-(1→6)-linked oligosaccharides: Hamiltonian replica exchange MD simulations and NMR experiments

Conformational sampling for a set of 10 α- or β-(1→6)-linked oligosaccharides has been studied using explicit solvent Hamiltonian replica exchange (HREX) simulations and NMR spectroscopy techniques. Validation of the force field and simulation methodology is done by comparing calculated transglycosidic J coupling constants and proton-proton distances with the corresponding NMR data. Initial calculations showed poor agreement, for example, with >3 Hz deviation of the calculated (3)J(H5,H6R) values from the experimental data, prompting optimization of the ω torsion angle parameters associated with (1→6)-linkages. The resulting force field is in overall good agreement (i.e., within ∼0.5 Hz deviation) from experimental (3)J(H5,H6R) values, although some small limitations are evident. Detailed hydrogen bonding analysis indicates that most of the compounds lack direct intramolecular H-bonds between the two monosaccharides; however, minor sampling of the O6···HO2' hydrogen bond is present in three compounds. The results verify the role of the gauche effect between O5 and O6 atoms in gluco- and manno-configured pyranosides causing the ω torsion angle to sample an equilibrium between the gt and gg rotamers. Conversely, galacto-configured pyranosides sample a population distribution in equilibrium between gt and tg rotamers, while the gg rotamer populations are minor. Water radial distribution functions suggest decreased accessibility to the O6 atom in the (1→6)-linkage as compared to the O6' atom in the nonreducing sugar. The role of bridging water molecules between two sugar moieties on the distributions of ω torsion angles in oligosaccharides is also explored.

Conformation and dynamics at a flexible glycosidic linkage revealed by NMR spectroscopy and molecular dynamics simulations: analysis of β-L-Fucp-(1→6)-α-D-Glcp-OMe in water solution

The intrinsic flexibility of carbohydrates facilitates different 3D structures in response to altered environments. At glycosidic (1→6)-linkages, three torsion angles are variable, and herein the conformation and dynamics of β-L-Fucp-(1→6)-α-D-Glcp-OMe are investigated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. The disaccharide shows evidence of conformational averaging for the ψ and ω torsion angles, best explained by a four-state conformational distribution. Notably, there is a significant population of conformations having ψ = 85° (clinal) in addition to those having ψ = 180° (antiperiplanar). Moderate differences in (13)C R1 relaxation rates are found to be best explained by axially symmetric tumbling in combination with minor differences in librational motion for the two residues, whereas the isomerization motions are occurring too slowly to be contributing significantly to the observed relaxation rates. The MD simulation was found to give a reasonably good agreement with experiment, especially with respect to diffusive properties, among which the rotational anisotropy, D∥/D⊥, is found to be 2.35. The force field employed showed too narrow ω torsion angles in the gauche-trans and gauche-gauche states as well as overestimating the population of the gauche-trans conformer. This information can subsequently be used in directing parameter developments and emphasizes the need for refinement of force fields for (1→6)-linked carbohydrates.

Glycan synthesis, structure, and dynamics: A selection

Glycan structural information is a prerequisite for elucidation of carbohydrate function in biological systems. To this end we employ a tripod approach for investigation of carbohydrate 3D structure and dynamics based on organic synthesis; different experimental spectroscopy techniques, NMR being of prime importance; and molecular simulations using, in particular, molecular dynamics (MD) simulations. The synthesis of oligosaccharides in the form of glucosyl fluorides is described, and their use as substrates for the Lam16A E115S glucosyl synthase is exemplified as well as a conformational analysis of a cyclic &beta;-(1&rarr;3)-heptaglucan based on molecular simulations. The flexibility of the N-acetyl group of aminosugars is by MD simulations indicated to function as a gatekeeper for transitions of glycosidic torsion angles to other regions of conformational space. A novel approach to visualize glycoprotein (GP) structures is presented in which the protein is shown by, for example, ribbons, but instead of stick or space-filling models for the carbohydrate portion it is visualized by the colored geometrical figures known as CFG representation in a 3D way, which we denote 3D-CFG, thereby effectively highlighting the sugar residues of the glycan part of the GP and the position(s) on the protein.

Glycan fragment database: a database of PDB-based glycan 3D structures

The glycan fragment database (GFDB), freely available at http://www.glycanstructure.org, is a database of the glycosidic torsion angles derived from the glycan structures in the Protein Data Bank (PDB). Analogous to protein structure, the structure of an oligosaccharide chain in a glycoprotein, referred to as a glycan, can be characterized by the torsion angles of glycosidic linkages between relatively rigid carbohydrate monomeric units. Knowledge of accessible conformations of biologically relevant glycans is essential in understanding their biological roles. The GFDB provides an intuitive glycan sequence search tool that allows the user to search complex glycan structures. After a glycan search is complete, each glycosidic torsion angle distribution is displayed in terms of the exact match and the fragment match. The exact match results are from the PDB entries that contain the glycan sequence identical to the query sequence. The fragment match results are from the entries with the glycan sequence whose substructure (fragment) or entire sequence is matched to the query sequence, such that the fragment results implicitly include the influences from the nearby carbohydrate residues. In addition, clustering analysis based on the torsion angle distribution can be performed to obtain the representative structures among the searched glycan structures.

Methyl α-l-rhamnosyl-(1→2)[α-l-rhamnosyl-(1→3)]-α-l-rhamnoside penta-hydrate: synchrotron study

The title hydrate, C₁₉H₃₄O₁₃·5H₂O, contains a vicinally disubstituted tris­accharide in which the two terminal rhamnosyl sugar groups are positioned adjacent to each other. The conformation of the tris­accharide is described by the glycosidic torsion angles ϕ2 = 48 (1)°, ψ2 = −29 (1)°, ϕ3 = 44 (1)° and ψ3 = 4 (1)°, whereas the ψ2 torsion angle represents a conformation from the major state in solution, the ψ3 torsion angle conformation may have been caught near a potential energy saddle-point when compared to its solution structure, in which at least two but probably three conformational states are populated. Extensive inter­molecular O—H⋯O hydrogen bonding is present in the crystal and a water-containing channel is formed along the b-axis direction.

Conformational flexibility of the pentasaccharide LNF-2 deduced from NMR spectroscopy and molecular dynamics simulations

Human milk oligosaccharides (HMOs) are important as prebiotics since they stimulate the growth of beneficial bacteria in the intestine and act as receptor analogues that can inhibit the binding of pathogens. The conformation and dynamics of the HMO Lacto-N-fucopentaose 2 (LNF-2), α-L-Fucp-(1 → 4)[β-D-Galp-(1 → 3)]-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, having a Lewis A epitope, has been investigated employing NMR spectroscopy and molecular dynamics (MD) computer simulations. 1D (1)H,(1)H-NOESY experiments were used to obtain proton-proton cross-relaxation rates from which effective distances were deduced and 2D J-HMBC and 1D long-range experiments were utilized to measure trans-glycosidic (3)J(CH) coupling constants. The MD simulations using the PARM22/SU01 force field for carbohydrates were carried out for 600 ns with explicit water as solvent which resulted in excellent sampling for flexible glycosidic torsion angles. In addition, in vacuo MD simulations were performed using an MM3-2000 force field, but the agreement was less satisfactory based on an analysis of heteronuclear trans-glycosidic coupling constants. LNF-2 has a conformationally well-defined region consisting of the terminal branched part of the pentasaccharide, i.e., the Lewis A epitope, and a flexible β-D-GlcpNAc-(1 → 3)-β-D-Galp-linkage towards the lactose unit, which is situated at the reducing end. For this β-(1 → 3)-linkage a negative ψ torsion angle is favored, when experimental NMR data is combined with the MD simulation in the analysis. In addition, flexibility on a similar time scale, i.e., on the order of the global overall molecular reorientation, may also be present for the ϕ torsion angle of the β-D-Galp-(1 → 4)-D-Glcp-linkage as suggested by the simulation. It was further observed from a temperature variation study that some (1)H NMR chemical shifts of LNF-2 were highly sensitive and this study indicates that Δδ/ΔT may be an additional tool for revealing conformational dynamics of oligosaccharides.

Overexpression of a homogeneous oligosaccharide with 13C labeling by genetically engineered yeast strain

This report describes a novel method for overexpression of (13)C-labeled oligosaccharides using genetically engineered Saccharomyces cerevisiae cells, in which a homogeneous high-mannose-type oligosaccharide accumulates because of deletions of genes encoding three enzymes involved in the processing pathway of asparagine-linked oligosaccharides in the Golgi complex. Using uniformly (13)C-labeled glucose as the sole carbon source in the culture medium of these engineered yeast cells, high yields of the isotopically labeled Man(8)GlcNAc(2) oligosaccharide could be successfully harvested from glycoprotein extracts of the cells. Furthermore, (13)C labeling at selected positions of the sugar residues in the oligosaccharide could be achieved using a site-specific (13)C-enriched glucose as the metabolic precursor, facilitating NMR spectral assignments. The (13)C-labeling method presented provides the technical basis for NMR analyses of structures, dynamics, and interactions of larger, branched oligosaccharides.

NMR analysis of conformationally dependent (n)J(C, H) and (n)J(C, C) in the trisaccharide α-L-Rhap-(1 → 2)[α-L-Rhap-(1 → 3)]-α-L-Rhap-OMe and a site-specifically labeled isotopologue thereof

An array of NMR spectroscopy experiments have been carried out to obtain conformationally dependent (1)H,(13)C- and (13)C,(13)C-spin-spin coupling constants in the trisaccharide α-L-Rhap-(1 → 2)[α-L-Rhap-(1 → 3)]-α-L-Rhap-OMe. The trisaccharide was synthesized with (13)C site-specific labeling at C2' and C2″, i.e. in the rhamnosyl groups in order to alleviate (1)H spectral overlap. This facilitated the measurement of a key trans-glycosidic proton-proton cross-relaxation rate using 1D (1)H,(1)H-T-ROESY experiments as well as a (3)J(C, H) coupling employing 1D (1)H,(13)C-long-range experiments, devoid of potential interference from additional J coupling. By means of both the natural abundance compound and the (13)C-labeled sample 2D (1)H,(13)C-J-HMBC and (1)H,(13)C-HSQC-HECADE NMR experiments, total line-shape analysis of (1)H NMR spectra and 1D (13)C NMR experiments were employed to extract (3)J(C, H) , (2)J(C, H), (3)J(C, C), and (1)J(C, C) coupling constants. The (13)C site-specific labeling facilitates straightforward determination of (n)J(C, C) as the splitting of the (13)C natural abundance resonances. This study resulted in eight conformationally dependent coupling constants for the trisaccharide and illustrates the use of (13)C site-specific labeling as a valuable approach that extends the 1D and 2D NMR methods in current use to attain both hetero- and homonuclear spin-spin coupling constants that subsequently can be utilized for conformational analysis.

Conformational properties of methyl β-maltoside and methyl α- and β-cellobioside disaccharides

An investigation of the conformational properties of methyl β-maltoside, methyl α-cellobioside, and methyl β-cellobioside disaccharides using NMR spectroscopy and molecular dynamics (MD) techniques, is presented. Emphasis is placed on validation of a recently presented force field for hexopyranose disaccharides followed by elucidation of the conformational properties of two different types of glycosidic linkages, α-(1 → 4) and β-(1 → 4). Both gas-phase and aqueous-phase simulations are performed to gain insight into the effect of solvent on the conformational properties. A number of transglycosidic J-coupling constants and proton-proton distances are calculated from the simulations and are used to identify the percent sampling of the three glycosidic conformations (syn, anti-φ, and anti-ψ) and, in turn, describe the flexibility around the glycosidic linkage. The results show the force field to be in overall good agreement with experiment, although some very small limitations are evident. Subsequently, a thorough hydrogen bonding analysis is performed to obtain insights into the conformational properties of the disaccharides. In methyl β-maltoside, competition between HO2'-O3 intramolecular hydrogen bonding and intermolecular hydrogen bonding of those groups with solvent leads to increased sampling of syn, anti-φ, and anti-ψ conformations and better agreement with NMR J-coupling constants. In methyl α- and β-cellobioside, O5'-HO6 and HO2'-O3 hydrogen bonding interactions are in competition with intermolecular hydrogen bonding involving the solvent molecules. This competition leads to retention of the O5'-HO3 hydrogen bond and increased sampling of the syn region of the φ/ψ map. Moreover, glycosidic torsions are correlated to the intramolecular hydrogen bonding occurring in the molecules. The present results verify that in the β-(1 → 4)-linkage intramolecular hydrogen bonding in the aqueous phase is due to the decreased ability of water to successfully compete for the O5' and HO3 hydrogen bonding moieties, in contrast to that occurring between the O5' and HO6 atoms in this α-(1 → 4)-linkage.