Distinguishing Anisotropy and Flexibility of the Pentasaccharide LNF-1 in Solution by Carbon-13 NMR Relaxation and Hydrodynamic Modeling

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.

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.

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

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

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.

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.

A perspective on the primary and three-dimensional structures of carbohydrates

Carbohydrates, in more biologically oriented areas referred to as glycans, constitute one of the four groups of biomolecules. The glycans, often present as glycoproteins or glycolipids, form highly complex structures. In mammals ten monosaccharides are utilized in building glycoconjugates in the form of oligo- (up to about a dozen monomers) and polysaccharides. Subsequent modifications and additions create a large number of different compounds. In bacteria, more than a hundred monosaccharides have been reported to be constituents of lipopolysaccharides, capsular polysaccharides, and exopolysaccharides. Thus, the number of polysaccharide structures possible to create is huge. NMR spectroscopy plays an essential part in elucidating the primary structure, that is, monosaccharide identity and ring size, anomeric configuration, linkage position, and sequence, of the sugar residues. The structural studies may also employ computational approaches for NMR chemical shift predictions (CASPER program). Once the components and sequence of sugar residues have been unraveled, the three-dimensional arrangement of the sugar residues relative to each other (conformation), their flexibility (transitions between and populations of conformational states), together with the dynamics (timescales) should be addressed. To shed light on these aspects we have utilized a combination of experimental liquid state NMR techniques together with molecular dynamics simulations. For the latter a molecular mechanics force field such as our CHARMM-based PARM22/SU01 has been used. The experimental NMR parameters acquired are typically (1)H,(1)H cross-relaxation rates (related to NOEs), (3)JCH and (3)JCCtrans-glycosidic coupling constants and (1)H,(13)C- and (1)H,(1)H-residual dipolar couplings. At a glycosidic linkage two torsion angles ϕ and ψ are defined and for 6-substituted residues also the ω torsion angle is required. Major conformers can be identified for which highly populated states are present. Thus, in many cases a well-defined albeit not rigid structure can be identified. However, on longer timescales, oligosaccharides must be considered as highly flexible molecules since also anti-conformations have been shown to exist with H-C-O-C torsion angles of ∼180°, compared to syn-conformations in which the protons at the carbon atoms forming the glycosidic linkage are in close proximity. The accessible conformational space governs possible interactions with proteins and both minor changes and significant alterations occur for the oligosaccharides in these interaction processes. Transferred NOE NMR experiments give information on the conformation of the glycan ligand when bound to the proteins whereas saturation transfer difference NMR experiments report on the carbohydrate part in contact with the protein. It is anticipated that the subtle differences in conformational preferences for glycan structures facilitate a means to regulate biochemical processes in different environments. Further developments in the analysis of glycan structure and in particular its role in interactions with other molecules, will lead to clarifications of the importance of structure in biochemical regulation processes essential to health and disease.

Stochastic modeling of flexible biomolecules applied to NMR relaxation. 2. Interpretation of complex dynamics in linear oligosaccharides

A computational stochastic approach is applied to the description of flexible molecules. By combining (i) molecular dynamics simulations, (ii) hydrodynamics approaches, and (iii) a multidimensional diffusive description for internal and global dynamics, it is possible to build an efficient integrated approach to the interpretation of relaxation processes in flexible systems. In particular, the model is applied to the interpretation of nuclear magnetic relaxation measurements of linear oligosaccharides, namely a mannose-containing trisaccharide and the pentasaccharide LNF-1. Experimental data are reproduced with sufficient accuracy without free model parameters.

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.

Sucrose in aqueous solution revisited, Part 2: adaptively biased molecular dynamics simulations and computational analysis of NMR relaxation

We report 100 ns molecular dynamics simulations, at various temperatures, of sucrose in water (with concentrations of sucrose ranging from 0.02 to 4M), and in a 7:3 water-DMSO mixture. Convergence of the resulting conformational ensembles was checked using adaptive-biased simulations along the glycosidic Φ and ψ torsion angles. NMR relaxation parameters, including longitudinal (R₁) and transverse (R₂) relaxation rates, nuclear Overhauser enhancements (NOE), and generalized order parameter (S²) were computed from the resulting time-correlation functions. The amplitude and time scales of molecular motions change with temperature and concentration in ways that track closely with experimental results, and are consistent with a model in which sucrose conformational fluctuations are limited (with 80-90% of the conformations having ϕ-ψ values within 20° of an average conformation), but with some important differences in conformation between pure water and DMSO-water mixtures.

Searching and optimizing structure ensembles for complex flexible sugars

NMR restrictions are suitable to specify the geometry of a molecule when a single well-defined global free energy minimum exists that is significantly lower than other local minima. Carbohydrates are quite flexible, and therefore, NMR observables do not always correlate with a single conformer but instead with an ensemble of low free energy conformers that can be accessed by thermal fluctuations. In this communication, we describe a novel procedure to identify and weight the contribution to the ensemble of local minima conformers based on comparison to residual dipolar couplings (RDCs) or other NMR observables, such as scalar couplings. A genetic algorithm is implemented to globally minimize the R factor comparing calculated RDCs to experiment. This is done by optimizing the weights of different conformers derived from the exhaustive local minima conformational search program, fast sugar structure prediction software (FSPS). We apply this framework to six human milk sugars, LND-1, LNF-1, LNF-2, LNF-3, LNnT, and LNT, and are able to determine corresponding population weights for the ensemble of conformers. Interestingly, our results indicate that in all cases the RDCs can be well represented by only a few most important conformers. This confirms that several, but not all of the glycosidic linkages in histo-blood group "epitopes" are quite rigid.

Molecular conformations in the pentasaccharide LNF-1 derived from NMR spectroscopy and molecular dynamics simulations

The conformational dynamics of the human milk oligosaccharide lacto-N-fucopentaose (LNF-1), α-L-Fucp-(1 → 2)-β-D-Galp-(1 → 3)-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, has been analyzed using NMR spectroscopy and molecular dynamics (MD) computer simulations. Employing the Hadamard (13)C-excitation technique and the J-HMBC experiment, (1)H,(13)C trans-glycosidic J coupling constants were obtained, and from one- and two-dimensional (1)H,(1)H T-ROESY experiments, proton-proton cross-relaxation rates were determined in isotropic D(2)O solution. In the lyotropic liquid-crystalline medium consisting of ditetradecylphosphatidylcholine, dihexylphosphatidylcholine, N-cetyl-N,N,N-trimethylammonium bromide, and D(2)O, (1)H, (1)H and one-bond (1)H, (13)C residual dipolar couplings (RDCs), as well as relative sign information on homonuclear RDCs, were determined for the pentasaccharide. Molecular dynamics simulations with explicit water were carried out from which the internal isomerization relaxation time constant, τ(N), was calculated for transitions at the ψ torsion angle of the β-(1 → 3) linkage to the lactosyl group in LNF-1. Compared to the global reorientation time, τ(M), of ∼0.6 ns determined experimentally in D(2)O solution, the time constant for the isomerization relaxation process, τ(N(scaled)), is about one-third as large. The NMR parameters derived from the isotropic solution show very good agreement with those calculated from the MD simulations. The only notable difference occurs at the reducing end, which should be more flexible than observed by the molecular simulation, a conclusion in complete agreement with previous (13)C NMR relaxation data. A hydrogen-bond analysis of the MD simulation revealed that inter-residue hydrogen bonds on the order of ∼30% were present across the glycosidic linkages to sugar ring oxygens. This finding highlights that intramolecular hydrogen bonds might be important in preserving well-defined structures in otherwise flexible molecules. An analysis including generalized order parameters obtained from nuclear spin relaxation experiments was performed and successfully shown to limit the conformational space accessible to the molecule when the number of experimental data are too scarce for a complete conformational analysis.

Comparing simulated and experimental translation and rotation constants: range of validity for viscosity scaling

Proper simulation of dynamic properties, including molecular diffusion, is an important goal of empirical force fields. However, the widely used TIP3P water model does not reproduce the experimental viscosity of water. Consequently, scaling of simulated diffusion constants of solutes in aqueous solutions is required to effectively compare them with experiment. It is proposed that scaling by the ratio of viscosities of model and real water is appropriate in the regime where the concentration dependence of simulated and experimental solution viscosities is parallel. With this ansatz, viscosity scaling can be carried out for glucose and trehalose up to 20 wt % for simulations carried out with the CHARMM additive carbohydrate force field C35 and TIP3P water; above this value, the concentration dependence of simulated viscosities lags that of experiment, and scaling is not advised. Scaled translational diffusion constants for glucose and the disaccharides trehalose, maltose, and melibiose at low concentration agree nearly quantitatively with experiment, as do NMR (13)C T(1)'s for glucose, trehalose, and maltose; these results support the use of C35 for simulations of sugar transport properties at low concentration. At high concentrations the scaled diffusion constants for glucose and trehalose underestimate and overestimate experiment, respectively. Hydrodynamic bead model calculations indicate a hydration level of approximately 1 water/hydroxyl for glucose. Patterns for the disaccharides are more complicated, though trehalose binds 0.5 to 1 more water than does maltose depending on the analysis.

Pathways and populations: stereoelectronic insights into the exocyclic torsion of 5-(hydroxymethyl)tetrahydropyran

High level ab initio computations in vacuum and with the IEFPCM implicit solvent model are carried out on 5-(hydroxymethyl)tetrahydropyran to investigate the effects of water on the exocyclic torsional surface. Rotamer populations evaluated from the omega(C-C-C-O), theta(C-C-C-O) solvent surface agree almost quantitatively with experimental values for the closely related methyl 4-deoxy-alpha-D-xylohexopyranoside. Potentials of mean force obtained from the two surfaces show substantial solvent stabilization of the TG (omega = 180 +/- 60 degrees) rotamer and the barriers at omega= 120 and 240 degrees but solvent destabilization at the cis barrier (omega = 0 degrees). Natural bond orbital analyses indicate that energetics of these effects are largely explained by overstabilization of the vacuum GT (omega= 60 +/- 60 degrees) and GG (omega = 300 +/- 60 degrees) rotamers. Solvent stabilization of theta conformations provides entropic stabilization.

Molecular dynamics studies of polyethylene oxide and polyethylene glycol: hydrodynamic radius and shape anisotropy

A revision (C35r) to the CHARMM ether force field is shown to reproduce experimentally observed conformational populations of dimethoxyethane. Molecular dynamics simulations of 9, 18, 27, and 36-mers of polyethylene oxide (PEO) and 27-mers of polyethylene glycol (PEG) in water based on C35r yield a persistence length lambda = 3.7 A, in quantitative agreement with experimentally obtained values of 3.7 A for PEO and 3.8 A for PEG; agreement with experimental values for hydrodynamic radii of comparably sized PEG is also excellent. The exponent upsilon relating the radius of gyration and molecular weight (R(g) proportional, variantM(w)(upsilon)) of PEO from the simulations equals 0.515 +/- 0.023, consistent with experimental observations that low molecular weight PEG behaves as an ideal chain. The shape anisotropy of hydrated PEO is 2.59:1.44:1.00. The dimension of the middle length for each of the polymers nearly equals the hydrodynamic radius R(h)obtained from diffusion measurements in solution. This explains the correspondence of R(h) and R(p), the pore radius of membrane channels: a polymer such as PEG diffuses with its long axis parallel to the membrane channel, and passes through the channel without substantial distortion.

The structural plasticity of heparan sulfate NA-domains and hence their role in mediating multivalent interactions is confirmed by high-accuracy (15)N-NMR relaxation studies

Considering the biological importance of heparan sulfate (HS) and the significant activity of its highly-sulfated regions (S-domains), the paucity of known functions for the non-sulfated NA-domains is somewhat puzzling. It has been suggested that chain dynamics within the NA-domains are the key to their functional role in HS. In this study, we investigate this hypothesis using state-of-the-art nuclear magnetic resonance (NMR) experiments at multiple frequencies. To resolve the problem of severe overlap in (1)H-NMR spectra of repetitive polysaccharides from proteoglycans, we have prepared oligosaccharides with the chemical structure of HS NA-domains containing the (15)N nucleus, which has enough chemical shift dispersion to probe the central residues of octasaccharides at atomic resolution using 600 MHz NMR. By performing NMR relaxation experiments at three magnetic-field strengths, high quality data on internal dynamics and rotational diffusion was obtained. Furthermore, translational diffusion could also be measured by NMR using pulse field gradients. These experimental data were used, in concert with molecular dynamics simulations, to provide information on local molecular shape, greatly aiding our relaxation analyses. Our results, which are more accurate than those presented previously, confirm the higher flexibility of the NA-domains as compared with reported data on S-domains. It is proposed that this flexibility has two functional roles. First, it confers a greater area of interaction from the anchoring point on the core protein for the bioactive S-domains. Secondly, it allows multiple interactions along the same HS chain that are dynamically independent of each other.

A conformational dynamics study of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution by NMR experiments and molecular simulations

The conformational preference of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution has been studied by NMR spectroscopy using one-dimensional (1)H,(1)H T-ROESY experiments and measurement of trans-glycosidic (3)J(C,H) coupling constants. Molecular dynamics (MD) simulations with a CHARMM22 type of force field modified for carbohydrates were performed with water as the explicit solvent. The homonuclear cross-relaxation rates, interpreted as effective proton-proton distances, were compared to those obtained from simulation. Via a Karplus torsional relationship, (3)J(C,H) values were calculated from simulation and compared to experimental data. Good agreement was observed between experimental data and the MD simulation, except for one inter-residue T-ROE between protons in the terminal sugar residues. The results show that the trisaccharide exhibits substantial conformational flexibility, in particular along the psi glycosidic torsion angles. Notably, for these torsions, a high degree of correlation (77%) was observed in the MD simulation revealing either psi(2)(+) psi(3)(+) or psi(2)(-)psi(3)(-) states. The simulations also showed that non-exoanomeric conformations were present at the phi torsion angles, but to a limited extent, with the phi(3) state populated to a larger extent than the phi(2) state. Further NMR analysis of the trisaccharide by translational diffusion measurements and (13)C T(1) relaxation experiments quantified global reorientation using an anisotropic model together with interpretation of the internal dynamics via the "model-free" approach. Fitting of the dynamically averaged states to experimental data showed that the psi(2)(+)psi(3)(+) state is present to approximately 49%, psi(2)(-) psi(3)(-) to approximately 39%, and phi(3) (non-exo) to approximately 12%. Finally, using a dynamic and population-averaged model, (1)H,(1)H T-ROE buildup curves were calculated using a full relaxation matrix approach and were found to be in excellent agreement with experimental data, in particular for the above inter-residue proton-proton interaction between the terminal residues.

Conformational analysis of a tetrasaccharide based on NMR spectroscopy and molecular dynamics simulations

The conformational preference of the human milk oligosaccharide lacto-N-neotetraose, beta-d-Galp-(1 --> 4)-beta-d-GlcpNAc-(1 --> 3)-beta-d-Galp-(1 --> 4)-d-Glcp, has been analyzed using (1)H,(1)H T-ROESY and (1)H,(13)C trans-glycosidic J coupling experiments in isotropic solution and (1)H,(13)C residual dipolar couplings (RDCs) obtained in lyotropic liquid crystalline media. Molecular dynamics simulations of the tetrasaccharide with explicit water as the solvent revealed that two conformational states are significantly populated at the psi glycosidic torsion angle, defined by C(anomeric)-O-C-H, of the (1 --> 3)-linkage. Calculation of order parameters, related to the molecular shape, were based on the inertia tensor and fitting of experimental RDCs to different conformational states showed that psi(+) > 0 degrees is the major and psi(-) < 0 degrees is the minor conformation in solution, in complete agreement with a two-state analysis based on the T-ROESY data. Attention was also given to the effect of salt (200 mM NaCl) in the anisotropic medium, which was a ternary mixture of n-octyl-penta(ethylene glycol), n-octanol, and D(2)O.

Dynamic properties of biologically active synthetic heparin-like hexasaccharides

A complete study of the dynamics of two synthetic heparin-like hexasaccharides, D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-1-->iPr (1) and -->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHAc-6-SO4-alpha-(1-->4)-L-IdoA-alpha-(1-->4)-D-GlcNHSO3-alpha-(1-->4)-L-IdoA-2-SO4-alpha-1-->iPr (2), has been performed using 13C-nuclear magnetic resonance (NMR) relaxation parameters, T1, T2, and heteronuclear nuclear Overhauser effect (NOEs). Compound 1 is constituted from sequences corresponding to the major polysaccharide heparin region, while compound 2 contains a sequence never found in natural heparin. They differ from each other only in sulphation patterns, and are capable of stimulating fibroblast growth factors (FGFs)-1 induced mitogenesis. Both oligosaccharides exhibit a remarkable anisotropic overall motion in solution as revealed by their anisotropic ratios (tau /tau||), 4.0 and 3.0 respectively. This is a characteristic behaviour of natural glycosaminoglycans (GAG) which has also been observed for the antithrombin (AT) binding pentasaccharide D-GlcNHSO3-6-SO4-alpha-(1-->4)-D-GlcA-beta-(1-->4)-D-GlcNHSO3-(3,6-SO4)-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-1-->Me (3) (Hricovíni, M., Guerrini, M., Torri, G., Piani, S., and Ungarelli, F. (1995) Conformational analysis of heparin epoxide in aqueous solution. An NMR relaxation study. Carbohydr. Res., 277, 11-23). The motional properties observed for 1 and 2 provide additional support to the suitability of these compounds as heparin models in agreement with previous structural (de Paz, J.L., Angulo, J., Lassaletta, J.M., Nieto, P.M., Redondo-Horcajo, M., Lozano, R.M., Jiménez-Gallego, G., and Martín-Lomas, M. (2001) The activation of fibroblast growth factors by heparin: synthesis, structure and biological activity of heparin-like oligosaccharides. Chembiochem, 2, 673-685; Ojeda, R., Angulo, J., Nieto, P.M., and Martin-Lomas. M. (2002) The activation of fibroblast growth factors by heparin: synthesis and structural study of rationally modified heparin-like oligosaccharides. Can. J. Chem,. 80, 917-936; Lucas, R., Angulo, J., Nieto, P.M., and Martin-Lomas, M. (2003) Synthesis and structural studies of two new heparin-like hexasaccharides. Org. Biomol. Chem., 1, 2253-2266) and biological data (Angulo, J., Ojeda, R., de Paz, J.L., Lucas, R., Nieto, P.M., Lozano, R.M., Redondo-Horcajo, M., Giménez-Gallego, G., and Martín-Lomas, M. (2004) The activation of fibroblast growth factors (FGFs) by glycosaminoglycans: influence of the sulphation pattern on the biological activity of FGF-1. Chembiochem, 5, 55-61). Fast internal motions observed for the less sulphated compound 2, as compared with 1, may be related to their different behavior in stimulating FGF1-induced mitogenic activity.

Dynamics in the cyclic Enterobacterial common antigen as studied by 13C NMR relaxation

The motional properties of the cyclic enterobacterial common antigen (cECA), consisting of four trisaccharide repeat units, have been investigated by carbon-13 spin relaxation. R(1), R(2) and NOE relaxation parameters have been determined at three magnetic field strengths. The data were interpreted within the model-free framework to include the possibility of motional anisotropy, and overall as well as local dynamical parameters were fitted separately for each ring carbon. The motional anisotropy was addressed by assuming an axially symmetric diffusion tensor, which was fitted from the overall correlation times for each site in the sugar residues using the previously determined crystal structure. The data were found to be in agreement with an oblate shape of the molecule, and the values for D(iso) and D(||)/D(perpendicular sign) were in good agreement with translational diffusion data and an estimate based on calculation of the moment of inertia tensor, respectively. The local dynamics in cECA were found to be residue-dependent. Somewhat lower values for the order parameters, as well as longer local correlation times, were observed for the beta-linked ManNAcA residue compared to the two alpha-linked residues in the trisaccharide repeat unit.

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.

Oligosaccharides Implicated in Recognition Are Predicted to Have Relatively Ordered Structures

Fucosylated O- and N-linked glycans are essential recognition molecules in plants and animals. To understand how they impart their functions, through interactions with proteins, requires a detailed analysis of structure and dynamics, but this is presently lacking. In this study, the three-dimensional structure and dynamics of three fucosylated oligosaccharides are investigated using a combination of high field (800 MHz) nuclear magnetic resonance and long (50 ns) molecular dynamics simulations in explicit water. Predictions from dynamics simulations were in agreement with nuclear Overhauser cross-peak intensities. Similarly, a theory of weak alignment in neutral media resulted in reasonable predictions of residual dipolar couplings for the trisaccharide fucosyllactose. However, for larger penta- and hexasaccharides (LNF-1 and LND-1), the anisotropic component of the alignment was underestimated, attributed to shape irregularities of the fucosyl branches on an otherwise linear core, being more pronounced in a singly branched than a doubly branched oligosaccharide. Simulations, confirmed by experiment, predicted fucosylated molecules that are restricted to librations about a single average conformation. This restriction is partly due to microscopic water interactions, which act to stabilize intramolecular hydrogen bonds and maintain tight and ordered conformations; a view not forthcoming from simpler, nonaqueous simulations. Such a conclusion is crucial for understanding how these molecules interact with proteins and impart their recognition properties.

Dynamics of the Escherichia coli O91 O-Antigen Polysaccharide in Solution as Studied by Carbon-13 NMR Relaxation

The dynamics of the O-antigen part of the lipopolysaccharide from the enterohemorrhagic Escherichia coli O91 has been determined in solution using (13)C NMR relaxation measurements at two magnetic field strengths, 9.4 and 14.1 T, thereby facilitating the testing of several dynamical models. The biological repeating unit, consisting of five sugar residues and substituents, could be determined by spectral analysis of different (1)H,(13)C correlations and corroborated by the relaxation data. The site specifically (13)C-labeled material was shown to have approximately 10 repeating units with a narrow distribution. A model-free analysis of the relaxation data revealed a complex dynamical behavior where the sugar residues could be described by a global correlation time (tau(m) = 5.4 ns), generalized order parameters (S(2) approximately 0.63), and different correlation times for internal motions related to their position in the repeating unit along the polymer (tau(e) approximately 360-520 ps). One of the sugar residues showed, in addition, a chemical exchange contribution. Furthermore, a substituent on another sugar residue was described by two order parameters (S(f)(2) = 0.51 and S(s)(2) = 0.21). The solution dynamics of the polysaccharide are thus described by highly intricate motions, both in amplitude and time scales. These results are of significance in the general description of polysaccharides surrounding bacterial cell surfaces and in the presentation of antigenic epitopes to the immune system of an invaded host.

Application of NMR, molecular simulation, and hydrodynamics to conformational analysis of trisaccharides

The preferred conformations and conformational flexibilities of the trisaccharides alpha-D-Glcp-(1-->2)-beta-D-Glcp-(1-->3)-alpha-D-Glcp-OMe (I) and alpha-D-Glcp-(1-->3)[beta-D-Glcp-(1-->4)]-alpha-D-Glcp-OMe (II) in aqueous solution were determined using nuclear magnetic resonance (NMR) spectroscopy, molecular dynamics (MD) and Langevin dynamics (LD) simulations, and hydrodynamics calculations. Both trisaccharides have a vicinal substitution pattern in which long range (nonsequential) interactions may play an important role. LD simulation at 600 K indicated that the all-syn conformation predominated, though other conformations were apparent. NOE data and MD and LD simulations at 298 K all indicated that trisaccharide I is a single all-syn conformer in solution. Given that previous studies showed evidence of anti-conformers in beta-D-Glcp-(1-->2)-beta-D-Glcp-(1-->3)-alpha-D-Glcp-OMe, this result provides an example of how changing the anomeric configuration of one residue from beta to alpha can make an oligosaccharide more rigid. Discrepancies in inter-ring distances obtained by experiment and by simulation of the all-syn conformer suggest the presence of an anti-psi conformation at the beta-(1-->4)-linkage for II. A combined analysis of measured and calculated translational diffusion constants and (13)C T(1) relaxation times yield order parameters of 0.9 for each trisaccharide. This implies that any interconversion among conformations is significantly slower than tumbling. Anisotropies of approximately 1.6 and 1.3 calculated for I and II, respectively, are consistent with the observed relatively flat T(1) profiles because the tumbling is not in the motional narrowing regime.

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.

Micelle-bound conformation of a hairpin-forming peptide: combined NMR and molecular dynamics study

A peptide fragment from a protein hairpin turn region was modified by addition of isoleucine residues to both ends to enhance binding to lipid micelles; the resulting peptide (I(1)-I(2)-C(3)-N(4)-N(5)-P(6)-H(7)-I(8)-I(9)) contains the core sequence I-C-N-N-P-H from an antibody-binding region of hemagglutinin A. Nuclear magnetic resonance (NMR) diffusion measurements indicated partial binding (43-65%) of the peptide to micelles of n-octylglucoside and significantly stronger binding (85%) to dodecylphosphocholine (DPC) micelles. Simulated annealing and conformational analysis using nuclear Overhauser enhancement restraints revealed a type I or III hairpin turn between residues N(5) and I(8) of the DPC-bound peptide. Amide exchange experiments support the possibility that a hydrogen bond forms between N(5) and I(8), stabilizing the turn. In contrast, no discernable structure was observed for the peptide in aqueous solution by either NMR or circular dichroism. Molecular dynamics simulations of DPC micelles and peptide-micelle complexes suggested that the peptide lies flat on the micelle surface and showed rapid rearrangement of the lipids to accommodate the bound peptide. According to a search performed using the basic local alignment search tool (BLAST), the sequences N-P-H-I and N-P-H-V are present as hairpin turns in eight of the nine proteins whose crystal structures were available. The addition of isoleucine residues and the use of lipid micelles to stabilize hairpin conformations equivalent to those found in proteins generates new possibilities for reproducing biologically important hairpin turns from short, linear peptides.

Dipole-dipole coupling constant for a directly bonded CH pair--a carbon-13 relaxation study

Multiple-field (4.7, 9.4, 14.1 T) carbon-13 relaxation data are reported for hexamethylenetetramine (HMTA) in the cryosolvent D(2)O/DMSO at 243 K. Under these conditions, the reorientational motion of HMTA is outside of the extreme narrowing range and the relaxation data can be subjected to a quantitative interpretation. Because of the high symmetry of the HMTA molecule, the reorientation must be isotropic. Treating the reorientation as a small-step rotational diffusion of a rigid body, we obtain a rotational correlation time of 1.0 ns and a carbon-proton dipole-dipole coupling constant corresponding to an effective internuclear distance of 114. 2 pm. The harmonic vibrational correction to the dipole-dipole coupling constant, based on a known force field, yields an NMR estimate of the r(alpha) distance of 110.8 +/- 0.3 pm.

An alternative method for pucker determination in carbohydrates from residual dipolar couplings: a solution NMR study of the fructofuranosyl ring of sucrose

A simple solution NMR method is presented for pucker determination of five-membered rings using only residual dipolar couplings obtained in a single liquid crystalline medium, DMPC/DHPC bicelles (DMPC = dimyristoylphosphatidylcholine; DHPC = dihexanoylphosphatidylcholine). The method was applied to determine the pucker of the fructofuranosyl ring of sucrose. The results indicate a fructofuranosyl pucker phase in the 20 degrees - 70 degrees range. The pucker phases are in agreement with those from previous NMR and optical spectroscopic studies and, importantly, do not rely on empirically parametrized Karplus curves. Furthermore, the analysis implies more than one stable pucker phase and rapid ring interconversion in this range. The present results suggest that using residual dipolar couplings alone can reveal multiple conformations present in solution. Hence, when a sufficient number of residual dipolar coupling constants is measured, the outcome is a robust, reliable, and independent route for determining carbohydrate and nucleic acid structure by NMR spectroscopy.

The diluted aqueous solvation of carbohydrates as inferred from molecular dynamics simulations and NMR spectroscopy

The purpose of this paper is to review our understanding of the dilute hydration (aqueous solvation) behaviour of disaccharide compounds. To this end we discuss and scrutinize the results that have been obtained for the three model disaccharides: maltose, sucrose and trehalose from experimental NMR studies and from theoretical molecular dynamics studies in explicit aqueous solutions. The focus is on the description of molecular hydration features that will influence macroscopic entities such as diffusion and relaxation: residence times of hydration waters, hydration numbers and hydration densities. The principles of molecular dynamics simulation are briefly outlined while a detailed presentation is given of the key features that characterise hydration: the solvation of the glycosidic linkage, the radial hydration of the solute, the water density anisotropy around the solute, the residential behaviour of water molecules in the periphery of the solute, and rotational and translational diffusion coefficients. With respect to the use of NMR in characterising the structure and dynamics of the hydration, the hydrodynamic theory of rotational and translational diffusion of biomolecules as well as the use of pulse field gradient spin echo experiments are briefly presented. The NMR-defined rotational diffusion coefficients (D(r)) and the experimentally determined translational diffusion (D(t)) coefficients are reported for 4% (w/w) solutions of sucrose, trehalose and maltose. These results are compared with theoretical data obtained from molecular dynamics simulations of sucrose, trehalose and maltose under identical conditions (concentration, temperature, etc.). With our present level of knowledge we can propose that although carbohydrates share a number of hydration characteristics, evidence is accumulating in support of the notion that it is not the amount or overall hydration but rather the detailed individual carbohydrate-water interaction that is likely to determine carbohydrate structure and functionality.

Oligosaccharides display both rigidity and high flexibility in water as determined by 13C NMR relaxation and 1H,1H NOE spectroscopy: evidence of anti-phi and anti-psi torsions in the same glycosidic linkage

The trisaccharide beta-D-Glcp-(1-->2)-beta-D-Glcp-(1-->3)-alpha-D-Glcp-OMe has been investigated by molecular dynamics (MD) simulations and NMR experiments in water. 13C spin-lattice (T1) and spin-spin (T2) relaxation times, together with 1H,13C NOE data were measured at two magnetic field strengths (9.4 and 14.1 T) in a 277 K D2O solution. Relaxation data interpreted by means of the model-free formalism revealed a rigid (S2 approximately 0.9) oligosaccharide tumbling in solution. 1H,1H Cross-relaxation rates were determined at 600 MHz by 1D DPFGSE NOESY and T-ROESY experiments, which provided high quality data and subsequently proton-proton distances within the trisaccharide. The presence of anti conformers at both torsions of a glycosidic linkage is demonstrated for the first time. MD simulations were carried out to facilitate analysis of the NOE data. In total, 15 simulations-starting from five different conformational states--were performed, with production runs of up to 10 ns, resulting in 83 ns of oligosaccharide dynamics in water. anti Conformers were populated to different degrees in the simulations, especially at the phi2 torsion angle. By combining the results from the NOE experiments and the MD simulations, the anti conformers at the (1-->2)-linkage were quantified as 7% anti-phi2 and 2% anti-psi2, revealing a highly flexible trisaccharide in which large conformational changes occur. From the MD simulations, interresidue hydrogen bonding, from HO2" to O2 or O3, was significantly populated (approximately 40%) in both of the anti conformational states. The contentious issue over rigidity versus flexibility in oligosaccharides has thus been thoroughly examined, showing that the dynamics should be taken into account for a relevant description of the molecular system.

Comparison of aqueous molecular dynamics with NMR relaxation and residual dipolar couplings favors internal motion in a mannose oligosaccharide

An investigation has been performed to assess how aqueous dynamical simulations of flexible molecules can be compared against NMR data. The methodology compares state-of-the-art NMR data (residual dipolar coupling, NOESY, and (13)C relaxation) to molecular dynamics simulations in water over several nanoseconds. In contrast to many previous applications of residual dipolar coupling in structure investigations of biomolecules, the approach described here uses molecular dynamics simulations to provide a dynamic representation of the molecule. A mannose pentasaccharide, alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->2)-D-Manp, was chosen as the model compound for this study. The presence of alpha-linked mannan is common to many glycopeptides, and therefore an understanding of the structure and the dynamics of this molecule is of both chemical and biological importance. This paper sets out to address the following questions. (1) Are the single structures which have been used to interpret residual dipolar couplings a useful representation of this molecule? (2) If dynamic flexibility is included in a representation of the molecule, can relaxation and residual dipolar coupling data then be simultaneously satisfied? (3) Do aqueous molecular dynamics simulations provide a reasonable representation of the dynamics present in the molecule and its interaction with water? In summary, two aqueous molecular dynamics simulations, each of 20 ns, were computed. They were started from two distant conformations and both converged to one flexible ensemble. The measured residual dipolar couplings were in agreement with predictions made by averaging the whole ensemble and from a specific single structure selected from the ensemble. However, the inclusion of internal motion was necessary to rationalize the relaxation data. Therefore, it is proposed that although residual dipolar couplings can be interpreted as a single-structure, this may not be a correct interpretation of molecular conformation in light of other experimental data. Second, the methodology described here shows that the ensembles from aqueous molecular dynamics can be effectively tested against experimental data sets. In the simulation, significant conformational motion was observed at each of the linkages, and no evidence for intramolecular hydrogen bonds at either alpha(1-->2) or alpha(1-->3) linkages was found. This is in contrast to simulations of other linkages, such as beta(1-->4), which are often predicted to maintain intramolecular hydrogen bonds and are coincidentally predicted to have less conformational freedom in solution.

Conformational analysis of a flexible oligosaccharide using residual dipolar couplings

We present a new approach to the analysis of the conformational and the motional properties of an oligosaccharide, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside. The approach relies on an order matrix analysis of residual dipolar couplings in the solution state. By combining a number of different types of couplings, (1)D(CH), (2)D(CH), and D(HH), an order matrix is solved for each ring of the trimannoside. The resulting order parameters indicate the internal motion at the alpha (1,3) linkage to be limited, while significant motion is suggested at the alpha (1,6) linkage. Two structures for the trimannoside were determined by aligning the order tensor principal axes obtained from two different orienting media, bicelles and phage. The very similar conformations at the alpha (1,3) linkage of these two structures confirm that the internal motion at the alpha (1,3) linkage is small and the conformation is a good representation of a single preferred structure. The different conformations at the alpha (1,6) linkage suggest that the motional amplitudes are large and the conformations must be viewed as virtual conformers. Compared with traditional NMR methods, data acquisition is easy and data analysis is straightforward.