Molecular dynamics simulations of a monofucosylated biantennary glycan of the N-acetyllactosamine type: the human lactotransferrin glycan.

Molecular dynamics simulations were carried out to explore the conformational flexibility of the antennae of N-linked glycans. They were performed over 200 ps in vacuo on the complete disialylated monofucosylated biantennary glycan of the N-acetyllactosaminic type. Starting from a bird-conformation, the 3-D-structure evolved through 9 successive transitional states to a new, compact and energetically favorable conformation which had never been previously described. In this conformation, both antennae are organized in two coplanar loops rolled in a contrary direction and oriented perpendicularly to the plane of the di-N-acetyl chitobiose residue leading to a 'lobster conformation'. All the glycosidic linkages of the disialylated monofucosylated biantennary glycan, except the Fuc(alpha 1-6)GlcNAc(beta 1-), were modified by a phase transition. Particularly, the Man(beta 1-4) GlcNAc(beta 1-) linkage, which was previously described by NMR and X-ray diffraction as a rigid one, was involved in numerous conformational changes during 83 ps, even before the first transition phase. The freedom of mobility of the torsional angles of the Man(alpha 1-6)Man(beta 1-) linkage was limited, under these simulation conditions, to the angle psi which took three values: 30 degrees, 90 degrees and 180 degrees. Moreover, from 150 ps up to the end of the simulation, the value of the torsional angle omega of the NeuAc(alpha 2-6)Gal(beta 1-) linkage of the alpha-1,6-antenna continuously swung between 60 degrees and -60 degrees. Finally, we observed that the values of the torsional angles of the three linkages: NeuAc(alpha 2-6)Gal(beta 1-), Gal(beta 1-4)GlcNAc(beta 1-) and GlcNAc(beta 1-2)Man(beta 1-) of each of the two antennae were different, demonstrating their asymmetric conformation.     

Molecular dynamics simulations and the conformational mobility of blood group oligosaccharides

Molecular dynamics simulations were carried out without explicit consideration of solvent to explore the conformational mobility of blood group A and H oligosaccharides. The potential energy force field of Rasmussen and co-workers was used with the CHARMM program on a number of disaccharide and trisaccharide models composed of fucose, galactose, glucose, N-acetyl glucosamine, and N-acetyl galactosamine chosen to represent various fragments of blood group oligosaccharides. In agreement with results of earlier studies, stable chair conformations were found for each pyranoside from which no transitions were detected in simulations as long as 800 ps. Exocyclic dihedral angles, including that of C5-C6, generally show numerous transitions on a time scale of approximately 5-30 ps. The dihedral angles of some but not all glycosidic linkages of blood group oligosaccharides show transitions on the time scale of 30-50 ps, implying that the extent of internal motion in blood group oligosaccharides depends strongly on linkage stereochemistry. For certain blood group A and H oligosaccharides that show limited internal motion in these simulations, we argue that the calculations are consistent with our previous analysis of 1H nuclear Overhauser enhancement (NOE) data that imply single conformations over a wide range of temperature and solvent conditions. While the trajectories are consistent with 13C T1 values that have been interpreted as indicating rigid conformations, measurements of 13C-NOE and T1 as a function of magnetic field strength are proposed as an improved method for experimental detection of the internal motion that is suggested for certain oligosaccharides in these simulations. The results of these simulations differ substantially from those of peptides of a similar molecular weight in that the oligosaccharides show much less internal motion.     

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.     

Molecular dynamics simulation of Lewis blood groups and related oligosaccharides.

Molecular dynamics simulations without explicit inclusion of solvent molecules have been performed to study the motions of Lewisa and Lewisb blood group oligosaccharides, and two blood group A tetrasaccharides having type I and type II core chains. The blood group H trisaccharide has also been studied and compared with the blood group A type II core chain. The potential energy surface developed by Rasmussen and co-workers was used with the molecular mechanics code CHARMM. The lowest energy minima of the component disaccharide fragments were obtained from conformational energy mapping. The lowest energy minima of these disaccharide fragments were used to build the tri- and tetrasaccharides that were further minimized before the actual heating/equilibration and dynamics simulations. The trajectories of the disaccharide fragments, e.g., Fuc alpha- (1----4)GlcNAc, Gal beta-(1----4)GlcNAc, etc., show transitions among various minima. However, the oligosaccharides were found to be dynamically stable and no transitions to other minimum energy conformations were observed in the time series of the glycosidic dihedral angles even during trajectories as long as 300 ps. The stable conformations of the glycosidic linkages in the oligosaccharides are not necessarily the same as the minimum energy conformation of the corresponding isolated disaccharides. The average fluctuations of the glycosidic angles in the oligosaccharides were well within the range of +/- 15 degrees. The results of these trajectory calculations were consistent with the relatively rigid single-conformation models derived for these oligosaccharides from 1H-nmr data.     

Molecular dynamics simulations of cyclohenicosakis-[(1-->2)-beta-D-gluco-henicosapyranosyl], a cyclic (1-->2)-beta-D-glucan (a 'cyclosophoraose') of DP 21

We report molecular dynamics simulations of cyclohenicosakis-[(1-->2)-beta-D-gluco-henicosapyranosyl], termed 'cyclosophohenicosamer', a member of a class of cyclic (1-->2)-beta-D-glucans ('cyclosophoraoses'). Our goals were to provide insights into the conformational preferences of these cyclosophoraoses. Simulated annealing and constant-temperature molecular dynamics calculations were performed on the DP 21 cyclosophohenicosamer. The radius of gyration (R(G)) of the molecule and the conformation of glycosidic dihedral angles were used to analyze the result of computational studies. Most glycosidic linkages were concentrated in the lowest-energy region of the phi-psi energy map, and the values of radius of gyration from our simulations were consistent with the reported experimental value. The simulations produced various types of compact and asymmetric conformations within reasonable ranges of the glycosidic linkage conformation and radius of gyration. The results indicate the presence of a high degree of molecular flexibility of cyclosophohenicosamer and suggest the uniqueness of inclusion complexation with other molecules through this molecular flexibility.     

Solution conformation of a pectic acid fragment by 1H-nmr and molecular dynamics

¹H-NMR and molecular dynamics simulations in vacuo and in water of (1 → 4)-α-D -galacturono-disaccharide were performed. The results of the molecular dynamics simulations showed that the molecule fluctuates between two conformations characterized by different values of torsion angles around the glycosidic linkage and two different intramolecular hydrogen bonds. When these conformations are extrapolated to a regular polymeric structure, they generate pectic acid compatible with a 2₁- or a right-handed 3₁-helix. © 1994 John Wiley & Sons, Inc.     

The solution structure of the intramolecular photoproduct of d(TpA) derived with the use of NMR and a combination of distance geometry and molecular dynamics.

One and two dimensional NMR techniques have been used together with molecular modelling to obtain the solution structure for the photoproduct d(TpA)*. The NMR data confirm that the cyclobutane linkage is formed between the bonds thymine C6-C5 and adenine C5-C6. The 2D NOE data are used as constraints in a distance geometry calculation. The structures obtained show a trans-syn cyclobutane linkage and the glycosidic angles are SYN and ANTI for thymidine and deoxyadenosine, respectively. The coupling constant data are used to check the backbone torsion angles of the obtained structures. Typical torsion angles are a gamma+ and beta t for the deoxyadenosine residue. A free molecular dynamics simulation of a trans-syn d(TpA) photoproduct confirmed all these structural characteristics.     

Conformational and dynamical properties of disaccharides in water: a molecular dynamics study

Explicit-solvent molecular dynamics simulations (50 ns, 300 K) of the eight reducing glucose disaccharides (kojibiose, sophorose, nigerose, laminarabiose, maltose, cellobiose, isomaltose, and gentiobiose) have been carried out using the GROMOS 45A4 force field (including a recently reoptimized carbohydrate parameter set), to investigate and compare their conformational preferences, intramolecular hydrogen-bonding patterns, torsional dynamics, and configurational entropies. The calculated average values of the glycosidic torsional angles agree well with available experimental data, providing validation for the force field and simulation methodology employed in this study. These simulations show in particular that: 1) (1-->6)-linked disaccharides are characterized by an increased flexibility, the absence of any persistent intramolecular hydrogen bond and a significantly higher configurational entropy (compared to the other disaccharides); 2) cellobiose presents a highly persistent interresidue hydrogen bond and a significantly lower configurational entropy (compared to the other disaccharides); 3) persistent hydrogen bonds are observed for all disaccharides (except (1-->6)-linked) and typically involve a hydrogen donor in the reducing residue and an acceptor in the nonreducing one; 4) the probability distributions associated with the glycosidic dihedral angles and psi are essentially unimodal for all disaccharides, and full rotation around these angles occurs at most once or twice for (never for psi) on the 50-ns timescale; and 5) the timescales associated with torsional transitions (except around and psi) range from approximately 30 ps (rotation of hydroxyl groups) to the nanosecond range (rotation of the lactol and hydroxymethyl groups, and around the omega-glycosidic dihedral angle in (1-->6)-linked disaccharides).     

Depicting the MM3 potential energy surfaces of trisaccharides by single contour maps: application to beta-cellotriose and alpha-maltotriose

The adiabatic potential energy surfaces (PES) of two trisaccharides (beta-cellotriose and alpha-maltotriose) were obtained using the MM3 force field. Each PES can be described by a single 3D contour map for which the energy is plotted against the two psi glycosidic angles. Given the usually small variations of the phi glycosidic torsional angle in the low-energy regions of disaccharide maps (at least with MM3), it is valid to leave both phi glycosidic angles to relax in the process of building the conformational map of trisaccharides. The surfaces are those expected from the map of disaccharides containing the same linkages and monosaccharide units (i.e., beta-cellobiose and alpha-maltose), with second-order factors altering the 'symmetry' of both linkages. A large low-energy region appears for beta-cellotriose, comprising four minima in close proximity, with barriers between them below 0.6 kcal/mol. On the other hand, for alpha-maltotriose a main global minimum is observed, with several surrounding local minima. The surfaces obtained agree with single-crystal X-ray data on these trisaccharides and derivatives. A reduction of the linkage flexibilities is observed when passing from the disaccharides to the trisaccharides. Furthermore, the linkage closer to the reducing end appears to be less flexible than the linkage closer to the non-reducing end.     

A processive carbohydrate polymerase that mediates bifunctional catalysis using a single active site

Even in the absence of a template, glycosyltransferases can catalyze the synthesis of carbohydrate polymers of specific sequence. The paradigm has been that one enzyme catalyzes the formation of one type of glycosidic linkage, yet certain glycosyltransferases generate polysaccharide sequences composed of two distinct linkage types. In principle, bifunctional glycosyltransferases can possess separate active sites for each catalytic activity or one active site with dual activities. We encountered the fundamental question of one or two distinct active sites in our investigation of the galactosyltransferase GlfT2. GlfT2 catalyzes the formation of mycobacterial galactan, a critical cell-wall polymer composed of galactofuranose residues connected with alternating, regioisomeric linkages. We found that GlfT2 mediates galactan polymerization using only one active site that manifests dual regioselectivity. Structural modeling of the bifunctional glycosyltransferases hyaluronan synthase and cellulose synthase suggests that these enzymes also generate multiple glycosidic linkages using a single active site. These results highlight the versatility of glycosyltransferases for generating polysaccharides of specific sequence. We postulate that a hallmark of processive elongation of a carbohydrate polymer by a bifunctional enzyme is that one active site can give rise to two separate types of glycosidic bonds.     

Molecular dynamics simulation and nmr study of a blood group H trisaccharide

Molecular dynamics simulations in vacuum and solution have been carried out on 2′-α-L -fucosyllactitol, a model for blood group H in conjunction with two-dimensional nmr measurements on the same compound. Three independent starting conformations for the dynamics were chosen from low energy conformations obtained by a ?/ψ grid search. Nine 5 ns vacuum simulations of the trisaccharide were performed, employing three different ways to treat electrostatic interactions for each starting conformation: distance-dependent dielectric with ε = r, constant dielectric with ε = 1, or constant dielectric with ε = 80. In vacuum, transitions of ? and ψ for the α-L -Fuc-(1 → 2)-β-D -Gal element occur in a cooperative manner. The virtual distance obtained for H1 in fucose to H2 in galactose from nuclear Overhauser effect spectroscopy experiments agree with one of the conformations of the trisaccharide in one of the three 100 ps aqueous simulations (?/ψ ca. ?100°/150°), indicating this may be a dominant solution conformation. The rms fluctuations of the ?- and ψ-dihedral angles were ～ 10° for a conformational state, both in the vacuum and the aqueous simulations. For the simulations in vacuum, the agreement with experimental NOE data is reasonable when a constant dielectric of 1 is used (major conformers having ?/ψ ca. ?100°/150° and ?140°/100°), whereas the agreement was poor with a constant dielectric of 80. Translational diffusion coefficients calculated from the simulation of the oligosaccharides were 0.12–0.18 × 10^?5 cm²/s and from nmr measurements 0.27 × 10^?5 cm²/s. © 1994 John Wiley & Sons, Inc.     

Conformational analysis and molecular dynamics simulation of α-(1 → 2) and α-(1 → 3) linked rhamnose oligosaccharides: Reconciliation with optical rotation and NMR experiments

Molecular mechanics and dynamics calculations were carried out on the disaccharides α-L-Rhap-(1 → 2)-α-L-Rhap-(1 → OMe) (1) and α-L-Rhap-(1 → 3)-α-L-Rhap-(1 OMe) (2), and the trisaccharide α-L-Rhap-(1 → 2)-α-L-Rhap-(1 → 3)-α-L-Rhap-(1 → OMe) (3). The semiflexible conformational behavior of these molecules was characterized by the occupation of a combination of different glycosidic linkage and side-chain conformational positions whose relative occupations were sensitive to dielectric screening. Molecular dynamics simulations of the trisaccharide 3 showed little difference between the linkage conformations in the trisaccharide and the component disaccharides 1 and 2. Experimental optical rotation data of 1 and 2 were obtained as a function of temperature in varying solvents. The molecular models were combined with the semiempirical theory of Stevens and Sathyanarayana to yield calculated optical rotations. Interpretation of the data of both 1 and 2 implied that a combination of conformations, both in glycosidic and side-chain positions, could explain the experimental data. Solvents effects were important in influencing the conformational mix and averaged optical rotation. Three-bond heteronuclear coupling constants ³J_{C, H} were obtained for the glycosidic linkages of 1 and 2 in D₂O and DMSO. Analysis of the coupling constants with a Karplus curve showed that small reductions in the glycosidic torsion angles of the conformations of the models used here of ca. 10°–15° in ϕ and 5°–10° in ψ were required to give better agreement with experiment; a combination of conformations for both 1 and 2 was consistent with the data. There was a negligible influence on the coupling constants of 1 on changing the solvent from D₂O to DMSO. © 1997 John Wiley & Sons, Inc.     

Conformational analysis of an alpha-galactosyl trisaccharide epitope involved in hyperacute rejection upon xenotransplantation.

alpha-Galactosyl epitopes are carbohydrate structures bearing an alpha-Gal-(1-->3)-Gal terminus (alpha-Gal epitopes). The interaction of these epitopes on the surface of animal cells with anti alpha-Gal antibodies in human serum is believed to be the main cause in antibody-mediated hyperacute rejection in xenotransplantation. In this paper, conformational analysis of an N-linked alpha-D-Galp-(1-->3)-beta-D-Galp-(1-->4)-beta-D-Glcp trisaccharide epitope was conducted in terms of each monosaccharide residue conformation, primary hydroxymethyl group configuration, and interglycosidic conformations. Selective 2D J-delta INEPT experiments have been carried out at three different temperatures to evaluate three-bond, long-range 13C-1H coupling constants for the crucial alpha-(1-->3) linkage. The NMR experimental data were complemented by theoretical calculations. The flexibility and dynamics of the trisaccharide have been studied by Metropolis Monte Carlo simulations. Ensemble-averaged three-bond, long-range 13C-1H coupling constants and nuclear Overhauser effects were in good agreement with the experimental data. The alpha-(1-->3) glycosidic linkage has shown a restricted flexibility as indicated by NMR spectroscopy and molecular modeling.     

A statistical analysis of N- and O-glycan linkage conformations from crystallographic data

We have generated a database of 639 glycosidic linkage structures by an exhaustive survey of the available crystallographic data for isolated oligosaccharides, glycoproteins, and glycan-binding proteins. For isolated oligosaccharides there is relatively little crystallographic data available. A much larger number of glycoprotein and glycan-binding protein structures have now been solved in which two or more linked monosaccharides can be resolved. In the majority of these cases, only a few residues can be seen. Using the 639 glycosidic linkage structures, we have identified one or more distinct conformers for all the linkages. The O5-C1-O-C(x)' torsion angles for all these distinct conformers appear to be determined chiefly by the exo-anomeric effect. The Manalpha1-6Man linkage appears to be less restrained than the others, showing a wide degree of dispersion outside the ranges of the defined conformers. The identification of distinct conformers for glyco-sidic linkages allows "average" glycan structures to be modeled and also allows the easy identification of distorted glycosidic linkages. Such an analysis shows that the interactions between IgG Fc and its own N-linked glycan result in severe distortion of the terminal Galbeta1-4GlcNAc linkage only, indicating the strong interactions that must be present between the Gal residue and the protein surface. The applicability of this crystallographic based analysis to glycan structures in solution is discussed. This database of linkagestructures should be a very useful reference tool in three-dimensional structure determinations.     

Conformational properties of glucose-based disaccharides investigated using molecular dynamics simulations with local elevation umbrella sampling

Explicit-solvent molecular dynamics (MD) simulations of the 11 glucose-based disaccharides in water at 300 K and 1 bar are reported. The simulations were carried out with the GROMOS 45A4 force-field and the sampling along the glycosidic dihedral angles ? and ψ was artificially enhanced using the local elevation umbrella sampling (LEUS) method. The trajectories are analyzed in terms of free-energy maps, stable and metastable conformational states (relative free energies and estimated transition timescales), intramolecular H-bonds, single molecule configurational entropies, and agreement with experimental data. All disaccharides considered are found to be characterized either by a single stable (overwhelmingly populated) state ((1→n)-linked disaccharides with n = 1, 2, 3, or 4) or by two stable (comparably populated and differing in the third glycosidic dihedral angle ; gg or gt) states with a low interconversion barrier ((1→6)-linked disaccharides). Metastable (anti-? or anti-ψ) states are also identified with relative free energies in the range of 8-22 kJ mol⁻¹. The 11 compounds can be classified into four families: (i) the α(1→1)α-linked disaccharide trehalose (axial-axial linkage) presents no metastable state, the lowest configurational entropy, and no intramolecular H-bonds; (ii) the four α(1→n)-linked disaccharides (n = 1, 2, 3, or 4; axial-equatorial linkage) present one metastable (anti-ψ) state, an intermediate configurational entropy, and two alternative intramolecular H-bonds; (iii) the four β(1→n)-linked disaccharides (n = 1, 2, 3, or 4; equatorial-equatorial linkage) present two metastable (anti-? and anti-ψ) states, an intermediate configurational entropy, and one intramolecular H-bond; (iv) the two (1→6)-linked disaccharides (additional glycosidic dihedral angle) present no (isomaltose) or a pair of (gentiobiose) metastable (anti-?) states, the highest configurational entropy, and no intramolecular H-bonds. The observed conformational preferences appear to be dictated by four main driving forces (ring conformational preferences, exo-anomeric effect, steric constraints, and possible presence of a third glycosidic dihedral angle), leaving a secondary role to intramolecular H-bonding and specific solvation effects. In spite of the weak conformational driving force attributed to solvent-exposed H-bonds in water (highly polar protic solvent), intramolecular H-bonds may still have a significant influence on the physico-chemical properties of the disaccharide by decreasing its hydrophilicity. Along with previous work, the results also complete the suggestion of a spectrum of approximate transition timescales for carbohydrates up to the disaccharide level, namely: ∼30 ps (hydroxyl groups), ∼1 ns (free lactol group, free hydroxymethyl groups, glycosidic dihedral angle in (1→6)-linked disaccharides), ∼10 ns to 2 μs (ring conformation, glycosidic dihedral angles ? and ψ). The calculated average values of the glycosidic torsional angles agree well with the available experimental data, providing validation for the force-field and simulation methodology employed.     

Solution conformations of a trimannoside from nuclear magnetic resonance and molecular dynamics simulations

N-linked oligosaccharides often act as ligands for receptor proteins in a variety of cell recognition processes. Knowledge of the solution conformations, as well as protein-bound conformations, of these oligosaccharides is required to understand these important interactions. In this paper we present a model for the solution conformations sampled by a simple trimannoside, methyl 3, 6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside, which contains two of the most commonly found glycosidic linkages in N-linked oligosaccharides. This model was derived from simulated annealing protocols incorporating distance restraints extracted from NOESY spectra along with torsional restraints computed from three-bond (1)H-(13)C coupling constants measured across the glycosidic bonds. The model was refined in light of unrestrained molecular dynamics simulations conducted in the presence of solvent water. The resulting model depicts a molecule undergoing conformational averaging in solution, adopting four major and two minor conformations. The four major conformations arise from a pair of two-state transitions, one each at the alpha(1-->3) and alpha(1-->6) linkages, whereas the minor conformations result from an additional transition of the alpha(1-->6) linkage. Our data also suggest that the alpha(1-->3) transition is fast and changes the molecular shape slightly, whereas the alpha(1-->6) is much slower and alters the molecular shape dramatically.     

Forced detachment of the CD2-CD58 complex

The force-induced detachment of the adhesion protein complex CD2-CD58 was studied by steered molecular dynamics simulations. The forced detachment of CD2 and CD58 shows that the system can respond to an external force by two mechanisms, which depend on the loading rate. At the rapid loading rates of 70 and 35 pN/ps (pulling speeds of 1 and 0.5 A/ps) the two proteins unfold before they separate, whereas at slower loading rates of 7 and 3.5 pN/ps (pulling speeds of 0.1 and 0.05 A/ps), the proteins separate before the domains can unfold. When subjected to a constant force of 400 pN, the two proteins separated without significant structural distortion. These findings suggest that protein unfolding is not coupled to the adhesive function of CD2 and CD58. The simulations further confirm that salt bridges primarily determine the tensile strength of the protein-to-protein bond, and that the order of salt bridge rupture depends mainly on the position of the bond, relative to the line of action of the applied force. Salt bridges close to this line break first. The importance of each of the salt bridges for adhesion, determined from the simulations, correlates closely with their role in cell-to-cell adhesion and equilibrium binding determined by site-directed mutagenesis experiments.     

Solution conformation of the branch points of N-linked glycans: synthetic model compounds for tri'-antennary and tetraantennary glycans

The solution conformation of model compounds for the tri'-antennary and tetraantennary (six-arm) branch point of N-linked glycans has been determined through the use of chemical shift, relaxation, and nuclear Overhauser enhancement data. The object was to establish the conformation about the glycosidic linkages in the N-linked substructure GlcNAc(beta 1,6) [GlcNAc(beta 1,2)] Man(alpha)- by estimation of values for the appropriate glycosidic torsional angles. The GlcNAc(beta 1,6) linkage in a trisaccharide model compound was found to be constrained to a narrow rotameric subpopulation about the substituted Man C5-C6 bond (omega = -60 degrees) and a narrow range of possible phi - psi values. Free rotation about the Man C5-C6 bond was obstructed by unfavorable steric interactions between the GlcNAc(beta 1,6) and GlcNAc(beta 1,2) residues. A phi, psi value of 55 degrees, 190 degrees was found to be consistent with the NMR data for the GlcNAc(beta 1,6) linkage. However, the value of psi appears to be "virtual" in that the molecule is in equilibrium between two different values (90 degrees and 252 degrees). For the GlcNAc(beta 1,2) linkage, complete agreement between all the observed NMR parameters and all the calculated ensemble average values could only be obtained with a set of potential energy functions which included hydrogen bonding. Other choices of potentials yielded calculated values that disagreed with at least two of the observed quantities. As a result, we infer that an interresidue hydrogen bond is formed, and we find it to be between the GlcNAc(beta 1,2) ring oxygen and the Man C3 hydroxyl.(ABSTRACT TRUNCATED AT 250 WORDS)     

Uncertainties in structural determinations of oligosaccharide conformation, using measurements of nuclear Overhauser effects

A fundamental problem in the determination of molecular structure by n.m.r. spectroscopy is insufficient experimental constraints. This problem is particularly marked for oligosaccharides, where few constraints are available across glycosidic linkages. By calculating distances as a function of dihedral angle, it is shown that, in general, two n.O.e. constraints result in two possible conformations for each glycosidic linkage, one of which can usually be discarded on the basis of model building or energy calculations. Using these calculations, an estimate of the uncertainty in the structure can be obtained.     

500-picosecond molecular dynamics in water of the Man alpha 1----2Man alpha glycosidic linkage present in Asn-linked oligomannose-type structures on glycoproteins

Molecular dynamics simulations of the Man alpha 1----2Man alpha glycosidic linkage found in the N-linked glycans of glycoproteins were performed in vacuo and in the presence of water. In the latter case significant dampening of the molecular fluctuations was found when compared to the in vacuo simulation. A 500-ps dynamics simulation in water showed only occasional short-lived deviations from the minimum-energy conformation, more consistent with carbohydrate "breathing" than flexibility. These studies add further evidence that oligosaccharides can maintain "fixed" geometries with relatively long lifetimes and are in agreement with experimental NMR-derived parameters for the same linkage in oligomannose structures.