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

Escherichia coli O176 LPS structure and dynamics: A NMR spectroscopy and MD simulation study

A lipopolysaccharide (LPS) molecule is a key component of the bacterial outer membrane used to protect the bacterium and to interact with the environment. To gain insight into its function, the study of the LPS conformation and dynamics at the molecular and cellular levels is necessary, but these highly diverse and dynamic membrane-LPS systems are difficult to study. In this work, by using NMR spectroscopy and molecular dynamics (MD) simulations, we determined the conformational preferences of an E. coli O176 O-antigen polysaccharide at the atomic level. Moreover, we analyzed the use of non-uniform sampling (NUS) for the acquisition of high dynamic range spectra, like ¹H,¹H-NOESY NMR experiments. A comparison of the effective transglycosidic distances derived from conventional uniformly sampled and NUS ¹H,¹H-NOESY data showed high similarity under equal measuring time conditions. Furthermore, the experimentally derived internuclear distances of the O-antigen polysaccharide with ten repeating units (RUs) showed very good agreement to those calculated from the MD simulations of the same O-antigen polysaccharide in solution. Analysis of the LPS bilayer simulations with five and with ten RUs revealed that, although similar with respect to populated states in solution, the O-antigen in LPS bilayers had more extended chains as a result of spatial limitations due to close packing. Additional MD simulations of O-antigen polysaccharides from E. coli O6 (branched repeating unit) and O91 (negatively charged linear repeating unit) in solution and LPS bilayers were performed and compared to those of O176 (linear polymer). For all three O-antigens, the ensemble of structures present for the polysaccharides in solution were consistent with the results from their ¹H,¹H-NOESY experiments. In addition, the similarities between the O-antigen on its own and as a constituent of the full LPS in bilayer environment makes it possible to realistically describe the LPS conformation and dynamics from the MD simulations.

•

Uniform and non-uniform sampled NOESY NMR data yield similar internuclear distances.

•

O-antigen internuclear distances from NMR and MD show excellent agreement.

•

O-antigen ensemble structures from MD are consistent with NMR observations.

•

O-antigen structures are more extended in LPS bilayers than in solution.

•

MD simulations can describe realistic LPS conformation and dynamics.

Review: Use of residual dipolar couplings to determine the structure of carbohydrates

Solution nuclear magnetic resonance spectroscopy is especially useful in the carbohydrate field. The measurement of residual dipolar couplings provides long-range structural information, a valuable complement for the structural study of carbohydrates either in its free form or in the bound state to proteins. They permit to deduce the geometry and the flexibility of the glycosidic linkages, which have a major influence on the conformation of carbohydrates and their overall shape. This article reviews the current application of the residual dipolar couplings methodology to carbohydrates.

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 1: molecular dynamics simulations and direct and indirect dipolar coupling analysis

Although the crystal structure of the disaccharide sucrose was solved more than 30 years ago, its conformational distribution in aqueous solution is still a matter of debate. We report here a variety of molecular dynamics simulations (mostly of 100 ns) using the GLYCAM06 force field and various water models, paying particular attention to comparisons to NMR measurements of residual dipolar couplings and electron-mediated spin-spin couplings. We focus on the glycosidic linkage conformation, the puckering phase angle of the fructose ring, and intramolecular hydrogen bonds between the two sugars. Our results show that sucrose is indeed a dynamic molecule, but the crystal conformation is qualitatively the dominant one in dilute solution. A second conformational basin, populated in many force fields, is probably overstabilized in the calculations.

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.

Structure and orientation of small molecules dissolved in the liquid crystalline phases of CsPFO/water system by multinuclear NMR

Pyridine, L-alanine and L-phenylalanine dissolved in the liquid crystalline phases of the lyotropic system CsPFO/water are studied by means of (2)H, (13)C and (1)H NMR. The orientational order parameters of the solutes are determined in a wide temperature range, together with some relevant geometrical parameters. In particular, a prevailing conformation for L-phenylalanine interacting with the micelle is suggested and, for all solutes, convincing representations of their specific interactions with the micelle surface are inferred. From (19)F and (13)C NMR spectra in the nematic phase, the orientational order parameters for the perfluorooctanoate chain inside the micelles are estimated.

A tool for the prediction of structures of complex sugars

In two recent back to back articles(Xia et al., J Chem Theory Comput 3:1620-1628 and 1629-1643, 2007a, b) we have started to address the problem of complex oligosaccharide conformation and folding. The scheme previously presented was based on exhaustive searches in configuration space in conjunction with Nuclear Overhauser Effect (NOE) calculations and the use of a complex rotameric library that takes branching into account. NOEs are extremely useful for structural determination but only provide information about short range interactions and ordering. Instead, the measurement of residual dipolar couplings (RDC), yields information about molecular ordering or folding that is long range in nature. In this article we show the results obtained by incorporation RDC calculations into our prediction scheme. Using this new approach we are able to accurately predict the structure of six human milk sugars: LNF-1, LND-1, LNF-2, LNF-3, LNnT and LNT. Our exhaustive search in dihedral configuration space combined with RDC and NOE calculations allows for highly accurate structural predictions that, because of the non-ergodic nature of these molecules on a time scale compatible with molecular dynamics simulations, are extremely hard to obtain otherwise (Almond et al., Biochemistry 43:5853-5863, 2004). Molecular dynamics simulations in explicit solvent using as initial configurations the structures predicted by our algorithm show that the histo-blood group epitopes in these sugars are relatively rigid and that the whole family of oligosaccharides derives its conformational variability almost exclusively from their common linkage (beta-D: -GlcNAc-(1-->3)-beta-D: -Gal) which can exist in two distinct conformational states. A population analysis based on the conformational variability of this flexible glycosidic link indicates that the relative population of the two distinct states varies for different human milk oligosaccharides.

Structural glycobiology: a game of snakes and ladders

Oligo- and polysaccharides are infamous for being extremely flexible molecules, populating a series of well-defined rotational isomeric states under physiological conditions. Characterization of this heterogeneous conformational ensemble has been a major obstacle impeding high-resolution structure determination of carbohydrates and acting as a bottleneck in the effort to understand the relationship between the carbohydrate structure and function. This challenge has compelled the field to develop and apply theoretical and experimental methods that can explore conformational ensembles by both capturing and deconvoluting the structural and dynamic properties of carbohydrates. This review focuses on computational approaches that have been successfully used in combination with experiment to detail the three-dimensional structure of carbohydrates in a solution and in a complex with proteins. In addition, emerging experimental techniques for three-dimensional structural characterization of carbohydrate-protein complexes and future challenges in the field of structural glycobiology are discussed. The review is divided into five sections: (1) The complexity and plasticity of carbohydrates, (2) Predicting carbohydrate-protein interactions, (3) Calculating relative and absolute binding free energies for carbohydrate-protein complexes, (4) Emerging and evolving techniques for experimental characterization of carbohydrate-protein structures, and (5) Current challenges in structural glycoscience.