Motional Properties of Two Vicinally Disubstituted Trisaccharides As Studied by Multiple-Field Carbon-13 NMR Relaxation

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.

Delineating the conformational flexibility of trisaccharides from NMR spectroscopy experiments and computer simulations

The conformation of saccharides in solution is challenging to characterize in the context of a single well-defined three-dimensional structure. Instead, they are better represented by an ensemble of conformations associated with their structural diversity and flexibility. In this study, we delineate the conformational heterogeneity of five trisaccharides via a combination of experimental and computational techniques. Experimental NMR measurements target conformationally sensitive parameters, including J couplings and effective distances around the glycosidic linkages, while the computational simulations apply the well-calibrated additive CHARMM carbohydrate force field in combination with efficient enhanced sampling molecular dynamics simulation methods. Analysis of conformational heterogeneity is performed based on sampling of discreet states as defined by dihedral angles, on root-mean-square differences of Cartesian coordinates and on the extent of volume sampled. Conformational clustering, based on the glycosidic linkage dihedral angles, shows that accounting for the full range of sampled conformations is required to reproduce the experimental data, emphasizing the utility of the molecular simulations in obtaining an atomic detailed description of the conformational properties of the saccharides. Results show the presence of differential conformational preferences as a function of primary sequence and glycosidic linkage types. Significant differences in conformational ensembles associated with the anomeric configuration of a single glycosidic linkage reinforce the impact of such changes on the conformational properties of carbohydrates. The present structural insights of the studied trisaccharides represent a foundation for understanding the range of conformations adopted in larger oligosaccharides and how these molecules encode their conformational heterogeneity into the monosaccharide sequence.

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

The conformational space available to the flexible molecule α-D-Manp-(1-->2)-α-D-Manp-OMe, a model for the α-(1-->2)-linked mannose disaccharide in N- or O-linked glycoproteins, is determined using experimental data and molecular simulation combined with a maximum entropy approach that leads to a converged population distribution utilizing different input information. A database survey of the Protein Data Bank where structures having the constituent disaccharide were retrieved resulted in an ensemble with >200 structures. Subsequent filtering removed erroneous structures and gave the database (DB) ensemble having three classes of mannose-containing compounds, viz., N- and O-linked structures, and ligands to proteins. A molecular dynamics (MD) simulation of the disaccharide revealed a two-state equilibrium with a major and a minor conformational state, i.e., the MD ensemble. These two different conformation ensembles of the disaccharide were compared to measured experimental spectroscopic data for the molecule in water solution. However, neither of the two populations were compatible with experimental data from optical rotation, NMR (1)H,(1)H cross-relaxation rates as well as homo- and heteronuclear (3)J couplings. The conformational distributions were subsequently used as background information to generate priors that were used in a maximum entropy analysis. The resulting posteriors, i.e., the population distributions after the application of the maximum entropy analysis, still showed notable deviations that were not anticipated based on the prior information. Therefore, reparameterization of homo- and heteronuclear Karplus relationships for the glycosidic torsion angles Φ and Ψ were carried out in which the importance of electronegative substituents on the coupling pathway was deemed essential resulting in four derived equations, two (3)J(COCC) and two (3)J(COCH) being different for the Φ and Ψ torsions, respectively. These Karplus relationships are denoted JCX/SU09. Reapplication of the maximum entropy analysis gave excellent agreement between the MD- and DB-posteriors. The information entropies show that the current reparametrization of the Karplus relationships constitutes a significant improvement. The Φ(H) torsion angle of the disaccharide is governed by the exo-anomeric effect and for the dominating conformation Φ(H) = -40 degrees and Ψ(H) = 33 degrees. The minor conformational state has a negative Ψ(H) torsion angle; the relative populations of the major and the minor states are approximately 3 : 1. It is anticipated that application of the methodology will be useful to flexible molecules ranging from small organic molecules to large biomolecules.

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.

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.

Dynamics of Hyaluronan Oligosaccharides Revealed by 15N Relaxation

Complex carbohydrates are considered to be flexible biomolecules, yet few experimental techniques are available to characterize their dynamics. In this study, we investigate the potential of 15N relaxation to probe the dynamics of hyaluronan oligosaccharides by adapting approaches previously applied to proteins. Unlike the 13C nucleus, 15N provides considerably enhanced spectral resolution, allowing position-specific information to be measured even in the middle of oligomers as large as decasaccharides. Moreover, isotopic incorporation maintains the 1H-15N group as an isolated spin-pair, allowing relaxation experiments to be performed and interpreted at low concentrations. A methodology is described for calculating the Lipari and Szabo model-free parameters at specific positions in hyaluronan oligomers and is used to produce a dynamic representation for the hexasaccharide. In this model, the glycosidic linkages and acetamido rotamer were determined to deviate by 18 degrees and 24 degrees from their mean positions, respectively. This approach allows the dynamic structural characterization of hyaluronan and other nitrogen-containing carbohydrates. The resultant models provide crucial insights into the physical properties and biology of these flexible molecules, which are at present poorly understood.

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.

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.

A Conformational Study of the Trisaccharideβ-D-Glcp-(1→2)[β-D-Glcp-(1→3)]α-D-Glcp-OMe by NMR NOESY and TROESY Experiments, Computer Simulations, and X-Ray Crystal Structure Analysis

Proton-proton cross-relaxation rates have been measured for the trisaccharide beta-D-Glcp-(l --> 2)[beta-D-Glcp-(1 --> 3)]alpha-D-Glcp-OMe in D2O as well as in D2O/[D6]DMSO 7:3 solution at 30 degrees C by means of one-dimensional NMR pulsed field gradient 1H,1H NOESY and TROESY experiments. Interatomic distances for the trisaccharide in D2O were calculated from the cross-relaxation rates for two intraresidue and three interglycosidic proton pairs, using the isolated spin-pair approximation. In the solvent mixture one intraresidue and three interglycosidic distances were derived without the use of a specific molecular model. In this case the distances were calculated from the cross-relaxation rates in combination with "model-free" motional parameters previously derived from 13C relaxation measurements. The proton-proton distances for interglycosidic pairs were compared with those averaged from Metropolis Monte Carlo and Langevin Dynamics simulations with the HSEA, PARM22, and CHEAT95 force fields. The crystal structure of the trisaccharide was solved by analysis of X-ray data. Interresidue proton pairs from the crystal structure and those observed by NMR experiments were similar. However, the corresponding proton-proton distances generated by computer simulations were longer. For the (1 --> 2) linkage the glycosidic torsion angles of the crystal structure were found in a region of conformational space populated by all three force fields, whereas for the (1 --> 3) linkage they occupied a region of low population density, as seen from the simulations.

Comparison of the conformation and dynamics of a polysaccharide and of its isolated heptasaccharide repeating unit on the basis of nuclear Overhauser effect, long-range C-C and C-H coupling constants, and NMR relaxation data

A comparison of the conformation and dynamics of the cell wall polysaccharide of S. mitis J22 and the heptasaccharide repeating unit made from this polysaccharide was performed on the basis on nmr data. We have previously reported a model for this highly flexible polysaccharide in which four residues of the antigenic epitope adopt a defined conformation as do the two residues of the lectin-binding epitope. These domains are connected by a 6-substituted galactofuranoside residue that acts as a flexible hinge and the repeating subunits are joined by phosphodiester linkages that provide further flexibility. Homonuclear nuclear Overhauser effect (NOE) and long-range C-C and C-H scalar coupling constants measured in uniform (13)C-labeled samples of the polysaccharide and heptasaccharide were very similar, indicating a similar conformational average in solution. Significant differences in the solution dynamics were found from the heteronuclear relaxation data, T(1), T(1 rho), and NOE, which reflect the faster molecular tumbling of the heptasaccharide. Internal motions occurring on a picosecond time scale are relatively uniform along the polymer while dynamics on the time scale longer than a few nanoseconds is characteristic of hinge motion.

13C high-resolution solid state NMR studies of the structure and dynamics of 1,6:3,4-dianhydro-2-O-tosyl-beta-D-galactopyranose

13C CP/MAS, dipolar dephasing MAS and theoretical GIAO calculations were employed to assign 13C resonances to the molecular structure of 1,6:3,4-dianhydro-2-O-tosyl-beta-D-galactopyranose 1. From spinning sideband intensities, employing the graphical method of Herzfeld and Berger the 13C delta(ii) parameters for aromatic residue were calculated. The experimental data were compared with computed results obtained by means of the B3PW91 hybrid method and 6-311G (df, p) basis set. The X-ray geometry of 1 with the correlated position of hydrogen atoms was taken as input data for theoretical calculations. As concluded from Cambridge Crystallographic Database (CSD) search, there are two reports describing the X-ray studies of 1 that show the slightly different geometry of the compound under investigation. This work shows that such discrepancies in geometry can generate differences between computed 13C delta(ii) parameters up to 6 ppm. 13C T1 and 1H T1rho relaxation times reveal that 1 is very rigid in crystal lattice. This structure is characterized by extremely long 1H T1rho, found to be in range ca. 200 ms.

A conformational study of the vicinally branched trisaccharide beta-D-glcp-(1 --> 2)[beta-D-glcp-(1 --> 3)]alpha-D-Manp-OMe by nuclear Overhauser effect spectroscopy (NOESY) and transverse rotating-frame Overhauser effect spectroscopy (TROESY) experiments: comparison to Monte Carlo and Langevin dynamics simulations

The trisaccharide beta-D-Glcp-(1 --> 2)[beta-D-Glcp-(1 --> 3)]alpha-D-Manp-OMe, a model for branching regions in oligosaccharides, has been investigated by one-dimensional DPFGSE (1)H, (1)H nuclear Overhauser effect spectroscopy (NOESY) and transverse rotating-frame Overhauser effect spectroscopy (TROESY) experiments at 30 degrees C in water and in the solvent mixture water : dimethyl sulfoxide (7 : 3). Cross-relaxation rates were obtained from the nmr experiments and interpreted as proton-proton distances. From Metropolis Monte Carlo and Langevin dynamics simulations, distances were calculated and compared to those obtained from experiment. Using the previously determined dynamics from carbon-13 nmr relaxation measurements of the trisaccharide in the solvent mixture, intraresidue proton distances could be obtained that were in excellent to reasonable agreement with those calculated from simulations. In water, the isolated spin-pair approximation was used for comparison of interproton distances. The experimentally derived distances in both solvents showed that the trans-glycosidic distances were shorter between the anomeric proton of the glucosyl group and the proton at the linkage position, respectively, than to the proton on the adjacent carbon on the mannosyl residue. The interresidue distances calculated from the computer simulations, performed with three different force fields, namely HSEA, PARM22, and CHEAT95, resulted in the reverse order in all cases but one.