Structure and dynamics of oligosaccharides: NMR and modeling studies

Probing altered receptor specificities of antigenically drifting human H3N2 viruses by chemoenzymatic synthesis, NMR, and modeling

Prototypic receptors for human influenza viruses are N-glycans carrying α2,6-linked sialosides. Due to immune pressure, A/H3N2 influenza viruses have emerged with altered receptor specificities that bind α2,6-linked sialosides presented on extended N-acetyl-lactosamine (LacNAc) chains. Here, binding modes of such drifted hemagglutinin's (HAs) are examined by chemoenzymatic synthesis of N-glycans having 13C-labeled monosaccharides at strategic positions. The labeled glycans are employed in 2D STD-1H by 13C-HSQC NMR experiments to pinpoint which monosaccharides of the extended LacNAc chain engage with evolutionarily distinct HAs. The NMR data in combination with computation and mutagenesis demonstrate that mutations distal to the receptor binding domain of recent HAs create an extended binding site that accommodates with the extended LacNAc chain. A fluorine containing sialoside is used as NMR probe to derive relative binding affinities and confirms the contribution of the extended LacNAc chain for binding.

Exploring multivalent carbohydrate-protein interactions by NMR

Nuclear Magnetic Resonance (NMR) has been widely employed to assess diverse features of glycan-protein molecular recognition events. Different types of qualitative and quantitative information at different degrees of resolution and complexity can be extracted from the proper application of the available NMR-techniques. In fact, affinity, structural, kinetic, conformational, and dynamic characteristics of the binding process are available. Nevertheless, except in particular cases, the affinity of lectin-sugar interactions is weak, mostly at the low mM range. This feature is overcome in biological processes by using multivalency, thus augmenting the strength of the binding. However, the application of NMR methods to monitor multivalent lectin-glycan interactions is intrinsically challenging. It is well known that when large macromolecular complexes are formed, the NMR signals disappear from the NMR spectrum, due to the existence of fast transverse relaxation, related to the large size and exchange features. Indeed, at the heart of the molecular recognition event, the associated free-bound chemical exchange process for both partners takes place in a particular timescale. Thus, these factors have to be considered and overcome. In this review article, we have distinguished, in a subjective manner, the existence of multivalent presentations in the glycan or in the lectin. From the glycan perspective, we have also considered whether multiple epitopes of a given ligand are presented in the same linear chain of a saccharide (i.e., poly-LacNAc oligosaccharides) or decorating different arms of a multiantennae scaffold, either natural (as in multiantennae N-glycans) or synthetic (of dendrimer or polymer nature). From the lectin perspective, the presence of an individual binding site at every monomer of a multimeric lectin may also have key consequences for the binding event at different levels of complexity.

Conformational and Structural characterization of carbohydrates and their interactions studied by NMR

Carbohydrates, either free or as glycans conjugated with other biomolecules, participate in many essential biological processes. Their apparent simplicity in terms of chemical functionality hides an extraordinary diversity and structural complexity. Deeply deciphering at the atomic level their structures is essential to understand their biological function and activities, but it is still a challenging task in need of complementary approaches and no generalized procedures are available to address the study of such complex, natural glycans. The versatility of Nuclear Magnetic Resonance spectroscopy (NMR) often makes it the preferred choice to study glycans and carbohydrates in solution media. The most basic NMR parameters, namely chemical shifts, coupling constants and nuclear Overhauser effects, allow defining short or repetitive chain sequences and characterize their structures and local geometries either in the free state or when interacting with other biomolecules, rendering additional information on the molecular recognition processes. The increased accessibility to carbohydrate molecules extensively or selectively labeled with 13C boosts the resolution and detail that analyzed glycan structures can reach. In turn, structural information derived from NMR, complemented with molecular modeling and theoretical calculations can also provide dynamic information on the conformational flexibility of carbohydrate structures. Furthermore, using partially oriented media or paramagnetic perturbations, it has been possible to introduce additional long-range observables rendering structural information on longer and branched glycan chains. In this review, we provide examples of these studies and an overview of the recent and most relevant NMR applications in the glycobiology field.

Structure of a protective epitope reveals the importance of acetylation of Neisseria meningitidis serogroup A capsular polysaccharide

Meningococcal meningitis remains a substantial cause of mortality and morbidity worldwide. Until recently, countries in the African meningitis belt were susceptible to devastating outbreaks, largely attributed to serogroup A Neisseria meningitidis (MenA). Vaccination with glycoconjugates of MenA capsular polysaccharide led to an almost complete elimination of MenA clinical cases. To understand the molecular basis of vaccine-induced protection, we generated a panel of oligosaccharide fragments of different lengths and tested them with polyclonal and monoclonal antibodies by inhibition enzyme-linked immunosorbent assay, surface plasmon resonance, and competitive human serum bactericidal assay, which is a surrogate for protection. The epitope was shown to optimize between three and six repeating units and to be O-acetylated. The molecular interactions between a protective monoclonal antibody and a MenA capsular polysaccharide fragment were further elucidated at the atomic level by saturation transfer difference NMR spectroscopy and X-ray crystallography. The epitope consists of a trisaccharide anchored to the antibody via the O- and N-acetyl moieties through either H-bonding or CH-π interactions. In silico docking showed that 3-O-acetylation of the upstream residue is essential for antibody binding, while O-acetate could be equally accommodated at three and four positions of the other two residues. These results shed light on the mechanism of action of current MenA vaccines and provide a foundation for the rational design of improved therapies.

Insight into Early-Stage Unfolding of GPI-Anchored Human Prion Protein

Prion diseases are fatal neurodegenerative disorders, which are characterized by the accumulation of misfolded prion protein (PrPSc) converted from a normal host cellular prion protein (PrPC). Experimental studies suggest that PrPC is enriched with α-helical structure, whereas PrPSc contains a high proportion of β-sheet. In this study, we report the impact of N-glycosylation and the membrane on the secondary structure stability utilizing extensive microsecond molecular dynamics simulations. Our results reveal that the HB (residues 173 to 194) C-terminal fragment undergoes conformational changes and helix unfolding in the absence of membrane environments because of the competition between protein backbone intramolecular and protein-water intermolecular hydrogen bonds as well as its intrinsic instability originated from the amino acid sequence. This initiation of the unfolding process of PrPC leads to a subsequent increase in the length of the HB-HC loop (residues 195 to 199) that may trigger larger rigid body motions or further unfolding around this region. Continuous interactions between prion protein and the membrane not only constrain the protein conformation but also decrease the solvent accessibility of the backbone atoms, thereby stabilizing the secondary structure, which is enhanced by N-glycosylation via additional interactions between the N-glycans and the membrane surface.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Paramagnetic NMR probes for characterization of the dynamic conformations and interactions of oligosaccharides

Paramagnetism-assisted nuclear magnetic resonance (NMR) techniques have recently been applied to a wide variety of biomolecular systems, using sophisticated immobilization methods to attach paramagnetic probes, such as spin labels and lanthanide-chelating groups, at specific sites of the target biomolecules. This is also true in the field of carbohydrate NMR spectroscopy. NMR analysis of oligosaccharides is often precluded by peak overlap resulting from the lack of variability of local chemical structures, by the insufficiency of conformational restraints from nuclear Overhauser effect (NOE) data due to low proton density, and moreover, by the inherently flexible nature of carbohydrate chains. Paramagnetic probes attached to the reducing ends of oligosaccharides cause paramagnetic relaxation enhancements (PREs) and/or pseudocontact shifts (PCSs) resolve the peak overlap problem. These spectral perturbations can be sources of long-range atomic distance information, which complements the local conformational information derived from J couplings and NOEs. Furthermore, paramagnetic NMR approaches, in conjunction with computational methods, have opened up possibilities for the description of dynamic conformational ensembles of oligosaccharides in solution. Several applications of paramagnetic NMR techniques are presented to demonstrate their utility for characterizing the conformational dynamics of oligosaccharides and for probing the carbohydrate-recognition modes of proteins. These techniques can be applied to the characterization of transient, non-stoichiometric interactions and will contribute to the visualization of dynamic biomolecular processes involving sugar chains.

Antibody recognition of carbohydrate epitopes†

Carbohydrate antigens are valuable as components of vaccines for bacterial infectious agents and human immunodeficiency virus (HIV), and for generating immunotherapeutics against cancer. The crystal structures of anti-carbohydrate antibodies in complex with antigen reveal the key features of antigen recognition and provide information that can guide the design of vaccines, particularly synthetic ones. This review summarizes structural features of anti-carbohydrate antibodies to over 20 antigens, based on six categories of glyco-antigen: (i) the glycan shield of HIV glycoproteins; (ii) tumor epitopes; (iii) glycolipids and blood group A antigen; (iv) internal epitopes of bacterial lipopolysaccharides; (v) terminal epitopes on polysaccharides and oligosaccharides, including a group of antibodies to Kdo-containing Chlamydia epitopes; and (vi) linear homopolysaccharides.

Synthesis of galactofuranosyl-(1 → 5)-thiodisaccharide glycomimetics as inhibitors of a β-d-galactofuranosidase

Description of the synthesis, molecular modeling and inhibitory properties of furanosyl thiodisaccharides that are mimetics of the motif β-d-Galf-(1 → 5)-d-Galf, found in glycoconjugates of pathogenic microorganisms.

Exploration of conformational spaces of high-mannose-type oligosaccharides by an NMR-validated simulation

Exploration of the conformational spaces of flexible biomacromolecules is essential for quantitatively understanding the energetics of their molecular recognition processes. We employed stable isotope- and lanthanide-assisted NMR approaches in conjunction with replica-exchange molecular dynamics (REMD) simulations to obtain atomic descriptions of the conformational dynamics of high-mannose-type oligosaccharides, which harbor intracellular glycoprotein-fate determinants in their triantennary structures. The experimentally validated REMD simulation provided quantitative views of the dynamic conformational ensembles of the complicated, branched oligosaccharides, and indicated significant expansion of the conformational space upon removal of a terminal mannose residue during the functional glycan-processing pathway.

Synthesis and conformational analysis of neoglycoconjugates derived from O- and S-glucose

Using olefin metathesis as a key step, four neoglycoconjugates incorporating α-O-glucose, α-S-glucose or β-S-glucose as a carbohydrate unit and L-serine or L-cysteine as an amino acid moiety have been synthesized. The four-atom carbon spacer allows the carbohydrate to explore a wide-ranging conformational space, which may have important implications for the molecular recognition of these molecules.

Solvation properties of N-acetyl-β-glucosamine: molecular dynamics study incorporating electrostatic polarization

N-Acetyl-β-glucosamine (NAG) is an important moiety of glycoproteins and is involved in many biological functions. However, conformational and dynamical properties of NAG molecules in aqueous solution, the most common biological environment, remain ambiguous due to limitations of experimental methods. Increasing efforts are made to probe structural properties of NAG and NAG-containing macromolecules, like peptidoglycans and polymeric chitin, at the atomic level using molecular dynamics simulations. In this work, we develop a polarizable carbohydrate force field for NAG and contrast simulation results of various properties using this novel force field and an analogous nonpolarizable (fixed charge) model. Aqueous solutions of NAG and its oligomers are investigated; we explore conformational properties (rotatable bond geometry), electrostatic properties (dipole moment distribution), dynamical properties (self-diffusion coefficient), hydrogen bonding (water bridge structure and dynamics), and free energy of hydration. The fixed-charge carbohydrate force field exhibits deviations from the gas phase relative rotation energy of exocyclic hydroxymethyl side chain and of chair/boat ring distortion. The polarizable force field predicts conformational properties in agreement with corresponding first-principles results. NAG-water hydrogen bonding pattern is studied through radial distribution functions (RDFs) and correlation functions. Intermolecular hydrogen bonding between solute and solvent is found to stabilize NAG solution structures while intramolecular hydrogen bonds define glycosidic linkage geometry of NAG oligomers. The electrostatic component of hydration free energy is highly dependent on force field atomic partial charges, influencing a more favorable free energy of hydration in the fixed-charge model compared to the polarizable model.

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.

N-Acetylated amino sugars: the dependence of NMR 3J(HNH2)-couplings on conformation, dynamics and solvent

N-Acetylated amino sugars are essential components of living organisms, but their dynamic conformational properties are poorly understood due to a lack of suitable experimental methodologies. Nuclear magnetic resonance (NMR) is ideally suited to these conformational studies, but accurate equations relating the conformation of key substituents (e.g., the acetamido group) to NMR observables are unavailable. To address this, density functional theory (DFT) methods have been used to calculate vicinal coupling constants in N-acetylated amino sugars and derive empirical Karplus equations for (3)J(H(N)H(2)) of N-acetyl-D-glucosamine (GlcNAc) and N-acetyl-D-galactosamine (GalNAc). The fitted Karplus parameters were found to be similar to those previously derived for peptide amide groups, but are consistently larger in magnitude. Local intramolecular interactions had a small effect on the calculated J-couplings and comparison with experimental data suggested that DFT slightly overestimated them. An implicit solvation model consistently lowered the magnitude of the calculated values, improving the agreement with the experimental data. However, an explicit solvent model, while having a small effect, worsened the agreement with experimental data. The largest contributor to experimentally-determined (3)J(H(N)H(2))-couplings is proposed to be librations of the amide group, which are well approximated by a Gaussian distribution about a mean dihedral angle. Exemplifying the usefulness of our derived Karplus equations, the libration of the amide group could be estimated in amino sugars from experimental data. The dynamical spread of the acetamido group in free alpha-GlcNAc, beta-GlcNAc and alpha-GalNAc was estimated to be 32 degrees , 42 degrees and 20 degrees , with corresponding mean dihedral angles of 160 degrees , 180 degrees and 146 degrees , respectively.

NMR analysis of carbohydrate-protein interactions

Carbohydrate-protein interactions are frequently characterized by dissociation constants in the microM to mM range. This is normally associated with fast dissociation rates of the corresponding complexes, in turn leading to fast exchange on the nuclear magnetic resonance (NMR) chemical shift time scale and on the NMR relaxation time scale. Therefore, NMR experiments that take advantage of fast exchange are well suited to study carbohydrate-protein interactions. In general, it is possible to analyze ligand binding by observing either protein signals or ligand resonances. Because most receptor proteins to which carbohydrates bind are rather large with molecular weights significantly exceeding 30 kDa, the analysis of the corresponding protein spectra is not trivial, and only very few studies have been addressing this issue so far. We, therefore, focus on NMR experiments that employ observation of free ligand, that is, carbohydrate signals to analyze the bound state. Two types of NMR experiments have been extremely valuable to analyze carbohydrate-protein interactions at atomic resolution. Whereas transferred nuclear Overhauser effect (NOE) experiments deliver bioactive conformations of carbohydrates binding to proteins, saturation transfer difference (STD) NMR spectra provide binding epitopes and valuable information about the binding thermodynamics and kinetics. We demonstrate the power of a combined transfer NOE/STD NMR approach for the analysis of carbohydrate-protein complexes using selected examples.

Structural view of glycosaminoglycan–protein interactions

The essential role of protein-glycosaminoglycan interactions in the regulation of various physiological processes has been recognized for several decades but it is only recently that the molecular basis underlying such interactions has emerged. The different methodologies to elucidate the three-dimensional features of glycosaminoglycans along with the interactions with proteins cover high resolution NMR spectroscopy, X-ray crystallography, molecular modeling, and hydrodynamic measurements. The structural results that have accumulated have been organized in databases that allow rapid searching with entries related either to the type of glycosaminoglycan or the type of protein. Finally, three selected examples enlightening the complexity of the nature of the interactions occurring between proteins and glycosaminoglycans are given. The example of interactions between heparin and antithrombin III illustrates how such a complex mechanism as the regulation of blood coagulation by a specific pentasaccharide can be dissected through the combined use of dedicated carbohydrate chemistry and structural glycobiology. The second example deals with the study of complexes between chemokines and heparin, and shows how multimolecular complexes of proteins can be organized in space throughout the action of glycosaminoglycans. Again, the synthesis of chemical mimetics offers an unexpected route to the development of novel glycotherapeutics. Finally, the area of enzymes/glycosaminoglycans complexes is briefly covered to realize the limited knowledge that we have for such an important class of biomacromolecular complexes.

Molecular dynamics and atomic charge calculations in the study of heparin conformation in aqueous solution

HF/6-31G** and molecular dynamics (MD) simulations were used to evaluate the performance of different atomic charge basis sets (i.e., Mulliken, Lowdin, and Electrostatic Potential Derived Charges--ESP) in heparin simulations. HF/3-21 G calculations were also used to study the NMR conformation of the IdoA residue. The results thus obtained indicated that ESP and Lowdin charges gave the better results in heparin simulations, followed by Mulliken charges, and that the minimum-energy conformation of IdoA can be different from that observed by NMR spectroscopy by less than 1 Angstrom. However, it was found that this small conformational modification is capable of inducing a change of almost 200 kJ/mol in the interactions of heparin with the surrounding environment, which is a meaningful amount of energy in the context of ligand-receptor interactions. This information can be potentially of great relevance in the design of heparin-derived antithrombotic compounds.

Studies of Disaccharide Solvation—Molecular Dynamics versus HPLC Retention

Molecular dynamics simulations have been used to assess the conformational behaviour of seven disaccharides in aqueous solution. Solvation decreased the overall conformational fluctuations of the sugars, compared to in vacuo simulations using a high dielectric constant. The most significant finding was a linear correlation between the experimental chromatographic retention parameter K´ and a molecular modelling parameter based on the next-nearest oxygen–oxygen distances in the disaccharides. The results support previous proposals for a stereospecific hydration model for carbohydrates and demonstrate the utility of a combined experimental/molecular modelling approach to its study.

Disialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling, Molecular Mechanics and Molecular Dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface disialogangliosides (GD1A, GD1B and GD3) in aqueous environment. The molecular mechanics calculation reveals that water mediated hydrogen bonding network plays a significant role in the structural stabilization of GD1A, GD1B and GD3. These water mediated hydrogen bonds not only exist between neighboring residues but also exist between residues that are separated by 2 to 3 residues in between. The conformational energy difference between different conformational states of gangliosides correlates very well with the number of water mediated and direct hydrogen bonds. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin is worked out. The NeuNAc has a limited allowed eulerian space at the binding site of Cholera Toxin (2.4%). The molecular modeling, molecular mechanics and molecular dynamics of disialoganglioside-cholera toxin complex reveal that cholera toxin can accommodate the disialoganglioside GD1A in three different modes. A single mode of binding is permissible for GD1B and GD3. Direct and water mediated hydrogen bonding interactions stabilizes these binding modes and play an essential role in defining the order of specificity for different disialogangliosides towards cholera toxin. This study not only provides models for the disialoganglioside-cholera toxin complexes but also identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.

Toward the understanding of the structure and dynamics of protein-carbohydrate interactions: molecular dynamics studies of the complexes between hevein and oligosaccharidic ligands

Herein we study, through all atom molecular dynamics simulations, the complex between hevein and two N-acetylated chitin oligomers, namely N,N(')-diacetylchitobiose and N,N('),N(")-triacetylchitotriose. The results of the simulations for two disaccharide complexes and one trisaccharide complex show that a carbohydrate oligomer is able to move on the surface of the relatively flat binding pocket of hevein, therefore occupying different binding subpockets. Statistical analysis methods were also applied in order to define the principal overall motions in the complexes, showing how the different ligands in the simulations modulate the protein motions. The oligosaccharide binding can be considered as defined by a subtle balance between enthalpic (formation of intermolecular interactions between the ligand and the receptor) and entropic (due mainly to the possibility for the sugar to move on the surface of the protein domain) effects, determining multiple binding conformations. This structural and dynamical view could parallel the results obtained by regularly used restrained MD simulations based on NOE NMR data that provide a well defined structure for both the disaccharide and trisaccharide complexes, and agrees with the observations for longer oligosaccharide chains.

Monosialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling and Molecular Mechanics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface monosialogangliosides (GM3, GM2 and GM1) in aqueous environment. Water mediated hydrogen bonding network plays a significant role in the structural stabilization of GM3, GM2 and GM1. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin reveals a limited allowed eulerian space of 2.4% with a much less allowed eulerian space (1.4%) for external galactose of GM1. The molecular mechanics of monosialoganglioside-cholera toxin complex reveals that cholera toxin can accommodate the monosialogangliosides in three different modes. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different monosialogangliosides towards cholera toxin. This study identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.

Hetaryleneaminopolyols and hetarylenecarbopeptoids: a new type of glyco- and peptidomimetics. Syntheses and studies on solution conformation and dynamics

Ready access to a new class of oligomers has been demonstrated by the synthesis of hetaryleneaminopolyols and hetarylenecarbopeptoids using 3-hydroxymethyl-5-(4-amino-4-deoxy-d-arabinotetritol-1-yl)-2-methylfuran and 5-(4-amino-4-deoxy-d-arabinotetritol-1-yl)-2-methyl-3-furoic acid as novel scaffolds. The conformational behavior of peptidomimetics 22, 23, 25, 26, and 36 have been analyzed by NMR spectroscopy and extensive molecular dynamics simulations. MD simulations using the GB/SA continuum solvent model for water and the MM3 force field provide a population distribution of conformers which satisfactorily agrees with the experimental NMR data for the torsional degrees of freedom of the molecule.

Conformation of glycomimetics in the free and protein-bound state: structural and binding features of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man

The conformational properties of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man in solution have been carefully analyzed by a combination of NMR spectroscopy and time-averaged restrained molecular dynamics. It has been found that both the alpha-1,3- and the alpha-1,6-glycosidic linkages show a major conformational averaging. Unusual Phi ca. 60 degrees orientations for both Phi torsion angles are found. Moreover, a major conformational distinction between the natural compound and the glycomimetic affects to the behavior of the omega(16) torsion angle around the alpha-1 --> 6-linkage. Despite this increased flexibility, the C-glycosyl analogue is recognized by three mannose binding lectins, as shown by NMR (line broadening, TR-NOE, and STD) and surface plasmon resonance (SPR) methods. Moreover, a process of conformational selection takes place, so that these lectins probably bind the glycomimetic similarly to the way they recognize the natural analogue. Depending upon the architecture and extension of the binding site of the lectin, loss or gain of binding affinity with respect to the natural analogue is found.

Synthesis and conformational analysis of a pentasaccharide corresponding to the cell-wall polysaccharide of the Group A Streptococcus

The synthesis and conformational analysis of a pentasaccharide corresponding to a fragment of the cell-wall polysaccharide (CWPS) of the bacteria Streptococcus Group A are described. The polysaccharide consists of alternating alpha-(1 --> 2)- and alpha-(1 --> 3)-linked L-rhamnopyranose (Rhap) residues with branching 2-acetamido-2-deoxy-D-glucopyranose (GlcpNAc) residues linked beta-(1 --> 3) to alternate rhamnose rings. The pentasaccharide is of interest as a possible terminal unit on the CWPS, for use in a vaccine. The syntheses employed a trichloroacetimidate glycosyl donor. Molecular dynamics (MD) calculations of the pentasaccharide with the force fields CVFF and PARM22, both in gas phase and with explicit water present, gave different predictions for the flexibility and preferred conformational space. Metropolis Monte Carlo (MMC) calculations with the HSEA force field were also performed. Experimental data were obtained from 1D transient NOE measurements. Complete build-up curves were compared to those obtained by full relaxation matrix calculations in order to derive a model of the conformation. Overall, the best fit between experimental and calculated data was obtained with MMC simulations using the HSEA force field. Molecular dynamics and MMC simulations of a tetrasaccharide corresponding to the Group A-variant polysaccharide, which differs in structure from Group A in lacking the GlcpNAc residues, were also performed for purposes of comparison.

Experimental evidence for the existence of non-exo-anomeric conformations in branched oligosaccharides: the neomycin-B case

For the first time in natural O-glycosides, a large amount of non-exo-anomeric conformation is experimentally detected, in solution.

Solution-state conformational study of the hevamine inhibitor allosamidin and six potential inhibitor analogues by NMR spectroscopy and molecular modeling

The solution-state conformations of the hevamine inhibitor allosamidin and six potential inhibitor analogues were studied by various NMR spectroscopic techniques and molecular modeling using force field calculations. Determination solely of the global energy minimum conformation was found to be insufficient for consensus with the NMR results, and agreement between the NMR experimental data and the theoretical calculations was only reached by assessing the structures as population-weighted average conformers on the basis of Boltzmann distributions derived from the calculated relative energies. The conformations of the glycosidic linkages in the compounds were found to be similar when the sugar residues were the same, but differences were markedly evident otherwise and also for the various heterocyclic group linkages. The binding of the compounds to hevamine, which may also complex to chitinases in general, was assessed using HMQC, transfer-NOESY, and both 1-D and 2-D saturation transfer difference NMR experiments. Under the conditions employed, only allosamidin was implicated to be bound to hevamine, and then only by HMQC with the dipolar coupling-based experiments failing to substantiate the formation of the complex. However, the results are consistent with the biochemical activities of the compounds whereby only allosamidin has been shown to act as a competitive inhibitor.

Application of the anomeric samarium route for the convergent synthesis of the C-linked trisaccharide alpha-D-Man-(1-->3)-[alpha-D-Man-(1-->6)]-D-Man and the disaccharides alpha-D-Man-(1-->3)-D-Man and alpha-D-Man-(1-->6)-D-Man

Studies are reported on the assembly of the branched C-trisaccharide, alpha-D-Man-(1-->3)-[alpha-D-Man-(1-->6)]-D-Man, representing the core region of the asparagine-linked oligosaccharides. The key step in this synthesis uses a SmI(2)-mediated coupling of two mannosylpyridyl sulfones to a C3,C6-diformyl branched monosaccharide unit, thereby assembling all three sugar units in one reaction and with complete stereocontrol at the two anomeric carbon centers. Subsequent tin hydride-based deoxygenation followed by a deprotection step produces the target C-trimer. In contrast to many of the other C-glycosylation methods, this approach employes intact carbohydrate units as C-glycosyl donors and acceptors, which in many instances parallels the well-studied O-glycosylation reactions. The synthesis of the C-disaccharides alpha-D-Man-(1-->3)-D-Man and alpha-D-Man-(1-->6)-D-Man is also described, they being necessary for the following conformational studies of all three carbohydrate analogues both in solution and bound to several mannose-binding proteins.

An effective strategy for structural elucidation of oligosaccharides through NMR spectroscopy combined with peracetylation using doubly 13C-labeled acetyl groups

The use of NMR spectroscopy for the elucidation of larger carbohydrate structures isolated from natural sources is principally limited by severe overlap of 1H signals, poor sensitivity when experiments involve 13C nuclei, and difficulties in conclusively establishing linkage positions. Peracetylation of oligosaccharides with doubly 13C-labeled acetyl groups provides several major advantages for their structural elucidation when combined with specifically tailored NMR pulse sequences. The 2.5–4.7 Hz J-coupling constants between acetyl carbonyl-13C nuclei and protons of the sugar ring at the sites of acetylation enables these sites to be readily assigned. By inference, glycosidic linkage positions on monosaccharides can be unambiguously determined. This can be used in lieu of permethylation analysis, yet does not require degradation of oligosaccharides. Spectral dispersion in the directly detected (1H) dimension is increased ~2.6–2.7-fold due to the downfield shifting of sugar-ring protons at the positions of acetylation. Peracetylation also introduces three new frequency dimensions for NMR studies, namely the 13CO, 13CMe, and 1HMe frequencies of the acetyl groups. These frequencies can be correlated to sugar protons, either independently or in combination, in alternative 2-, 3-, or 4-D experiments. The use of Hartmann–Hahn coherence transfer combined with zero-quantum dephasing periods permits purely absorptive in-phase multiplets to be extracted and enables accurate scalar couplings between ring protons to be measured, even in multidimensional experiments. Results are illustrated on a nonasaccharide-alditol derived from N-linked glycoproteins and on some smaller structures containing sialic acids and N-acetylhexosamines. Methods for small-scale sample acetylation using the superacylation catalyst, 4-dimethylamino pyridine, are described. A brief historical perspective pertinent to the fundamental contributions of Dr. R.U. Lemieux to the field of carbohydrate NMR is also presented.Key words: NMR, oligosaccharides, peracetylation, doubly 13C-labeled acetyl groups, tailored pulse sequences, heteronuclear Hartmann–Hahn.

How to detect internal motion by homonuclear NMR

NOESY and ROESY cross-peak intensities depend on internuclear distances and internal motion. Internal motion is usually ignored, and NOESY cross-peak intensities are interpreted in terms of internuclear distances only. Off-resonance ROESY experiments measure a weighted average of NOE and ROE. The weight can be described and experimentally set by an angle theta;. For large enough molecules, NOE and ROE have opposite signs. Therefore, each cross-peak intensity becomes zero for an angle theta;(0). For any sample, the maximum angle theta;(0) is determined by the overall motion of the molecule. Smaller theta;(0) values reflect the angular component of internal motions. Because individual cross-peaks are analyzed, the method evaluates internal motions of individual H,H vectors. The reduction of theta;(0) is largest for internal motions on a time scale of 100-300 ps. The sensitivity of theta;(0) for internal motions decreases with increasing molecular weight. We estimate that detecting internal motions will be practicable for molecules up to about 15 kDa. We describe a protocol to measure theta;(0) from a series of off-resonance ROESY spectra. For such a series, we describe the choice of experimental parameters, a procedure to extract theta;(0) from the raw data, and the interpretation of theta;(0) in terms of internal motions. In the small protein BPTI, we analyzed 75 cross-peaks. The precision of theta;(0) was 0.25 degrees, as compared to typical reductions of theta;(0) of 3 degrees. We found a well-defined maximum theta;(0) for cross-peaks in rigid parts of the molecule, which reflects the overall motion of the molecule. For BPTI, also many structurally important long-range cross-peaks appear rigid. The lower theta;(0) values of long-range contacts involving methyl groups are consistent with methyl rotation on the 25-ps time scale. The lower theta;(0) values of the flexible C-terminus and of flexible side chains translate into upper limits for the angular order parameter of 0.4 and 0.5-0.8, respectively. Off-resonance ROESY can monitor internal motions of H,H contacts that are used in a structure calculation. Because no isotope labeling is needed, off-resonance ROESY can be used to detect internal motions in a wide range of natural products.

Conformational selection of glycomimetics at enzyme catalytic sites: experimental demonstration of the binding of distinct high-energy distorted conformations of C-, S-, and O-glycosides by E. Coli beta-galactosidases

We show that the conformational features of the molecular complexes of E. coli beta-galactosidase and O-glycosides may differ from those formed with closely related compounds in their chemical nature, such as C- and S-glycosyl analogues. In the particular case presented here, NMR and ab initio quantum mechanical results show that the 3D-shapes of the ligand/inhibitor within the enzyme binding site depend on the chemical nature of the compounds. In fact, they depend on the relative size of the stereoelectronic barriers for chair deformation or for rotation around Phi glycosidic linkage.

Molecular dynamics simulations of alpha2 --> 8-linked disialoside: conformational analysis and implications for binding to proteins

Computational methods have played a key role in elucidating the various three-dimensional structures of oligosaccharides. Such structural information, together with other experimental data, leads to a better understanding of the role of oligosaccharide in various biological processes. The disialoside Neu5Ac-alpha2-->8-Neu5Ac appears as the terminal glycan in glycoproteins and glycolipids, and is known to play an important role in various events of cellular communication. Neurotoxins such as botulinum and tetanus require Neu5Ac-alpha2 --> 8-Neu5Ac for infecting the host. Glycoconjugates containing this disialoside and the enzymes catalyzing their biosynthesis are also regulated during cell growth, development, and differentiation. Unlike other biologically relevant disaccharides that have only two linkage bonds, the alpha2-->8-linked disialoside has four: C2-O, O-C8', C8'-C7', and C7'-C6'. The present report describes the results from nine 1 ns MD simulations of alpha2-->8-linked disialoside (Neu5Ac-alpha2-->8-Neu5Ac); simulations were run using GROMOS96 by explicitly considering the solvent molecules. Conformations around the O-C8' bond are restricted to the +sc/+ap regions due to stereochemical reasons. In contrast, conformations around the C2-O and C8'-C7' bonds were found to be largely unrestricted and all the three staggered regions are accessible. The conformations around the C7'-C6' bond were found to be in either the -sc or the anti region. These results are in excellent agreement with the available NMR and potential energy calculation studies. Overall, the disaccharide is flexible and adopts mainly two ensembles of conformations differing in the conformation around the C7'-C6' bond. The flexibility associated with this disaccharide allows for better optimization of intermolecular contacts while binding to proteins and this may partially compensate for the loss of conformational entropy that may be incurred due to disaccharide's flexibility.