About this Journal Submit a Manuscript Table of Contents
International Journal of Carbohydrate Chemistry
Volume 2012 (2012), Article ID 725271, 7 pages
http://dx.doi.org/10.1155/2012/725271
Research Article

Conformational Search on the Lewis X Structure by Molecular Dynamic: Study of Tri- and Pentasaccharide

Laboratory of Natural Products and Bio Actives, Department of Chemistry, Faculty of Sciences, Aboubakr Belkaid University, Tlemcen 13000, Algeria

Received 4 July 2011; Revised 7 October 2011; Accepted 12 October 2011

Academic Editor: Yuriy A. Knirel

Copyright © 2012 N. Khebichat et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Carbohydrates play vital roles in many biological processes, such as recognition, adhesion, and signalling between cells. The Lewis X determinant is a trisaccharide fragment implicated as a specific differentiation antigen, tumor antigen, and key component of the ligand for the endothelial leukocyte adhesion molecule, so it is necessary or essential to determine and to know their conformational and structural properties. In this work, conformational analysis was performed using molecular dynamics (MD) simulation with the AMBER10 program package in order to study the dynamic behavior of of the Lewis X trisaccharide (β-D-Gal-(1,4)-[α-L-Fuc-(1,3)]-β-D-GlcNAc-OMe) and the Lewis X pentasaccharide (β-D-Gal-(1,4)-[α-L-Fuc-(1,3)]-β-D-GlcNAc-(1,3)-β-D-Gal-(1,4)-β-D-Glu-OMe) in explicit water model at 300 K for 10 ns using the GLYCAM 06 force field.

1. Introduction

In nature, carbohydrates form an important family of biomolecules, as simple or complex carbohydrates, either alone or covalently linked to proteins or lipids [1]. They play vital roles in many biological processes, such as recognition, adhesion, and communication between cells. Carbohydrates that are covalently linked to a nonsugar moiety (proteins, peptides, or lipids) are the most prominent cell-surface-exposed structures. The highly diverse structural variability of carbohydrates makes them good candidates for cell receptors and recognition molecules [2].

Determination of the three-dimensional structure of oligosaccharides and understanding the molecular basis of their recognition by receptors represent the main challenges of structural glycobiology. Elucidation of the three-dimensional structures and the dynamical properties of oligosaccharides is a prerequisite for a better understanding of the biochemistry of recognition processes and for the rational design of carbohydrate-derived drugs [1].

The human histo-blood group ABH (O) systems were the first major human alloantigens to be identified. The carbohydrate nature of the A, B, H, and Lewis antigens was established in 1950s [3, 4]. Figure 1 shows a schematic representation of current histo-blood group antigens although additional rare types also exist [5]. The oligosaccharidic epitope of ABH (O) and Lewis histo-blood groups have been the subject of many structural investigations. The energy maps of all of the disaccharides have been established by molecular mechanics methods [6, 7].

fig1
Figure 1: Schematic representation of histo-blood group oligosaccharides.

The determination of conformational preferences of oligosaccharides is best approached by describing their preferred conformations on potential energy surfaces as a function of the glycosidic linkage , torsional angles. In the early days of conformational analysis of oligosaccharides, the calculations of either rigid or, later, relaxed conformational maps predominated. Today, complex simulation techniques such as molecular dynamics are increasingly being used. Several examples can be found in the literature of “state-of-the-art” conformational search strategies in order to simplify the complexity of the multiple-minima problem [8].

In this work, we carried out extensive, explicit solvent MD simulation studies for tri- and pentasaccharide Lewis X for 10 ns. Each monosaccharide unit in the Lewis X molecules adopts the typical 4C1 chair conformation, with no significant deviation from the classical pyranose ring shape. The objective of this study was to explore the conformational behaviour of those oligosaccharides in explicit water model at 300 K for 10 ns using the GLYCAM 06 force field.

2. Theoretical Approach

The atomic labelling and the torsion angles of interest are given in Figure 2. The torsion angles are defined (05-Cl-0′2-C′2), (Cl-0′2-C′2-C′3), and (O5-C5-C6-O6). The conformations with specific torsion angles ω are referred to as gauche-trans (gt; ), gauche-gauche (gg; ) and trans-gauche (tg; ).

725271.fig.002
Figure 2: Atomic labeling and torsion angles of interest.

The MD simulations were carried out using AMBER10 [9] at a temperature of 300 K. This temperature was chosen because our primary objective was to study the dynamics of the oligosaccharides at physiologically relevant conditions. The starting Lewis X pentasaccharide and trisaccharide monomer coordinates and parameters were constructed with the “Glycam Biomolecule Builder” available online from the website of Woods group [10] (see Schemes 1 and 2).

725271.sch.001
Scheme 1: Lewis X trisaccharide.
725271.sch.002
Scheme 2: Lewis X pentasaccharide.

Since these coordinates are not optimized, a first minimization by molecular mechanics (1,000 steps of steepest descent followed by 1,000 steps of conjugate gradient minimization with a gradient tolerance of 0.05 kcal mol−1) was made with AMBER package version 10 [9]. Force field GLYCAM 06 was used for oligosaccharides (it used to simulate the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics and is currently suitable for all biologically relevant N- and O-linkages in oligo- and polysaccharides and glycoproteins [11]). The simulations were carried out using an explicit solvent environment with AMBER package version 10. The solvated system for the oligosaccharides molecules was prepared using an explicit water box of 8 Å in the xleap module of AMBER10. The number of water molecules was 569 for the trisaccharide and 1072 for the pentasaccharide. All simulations were run with the SANDER module of AMBER with SHAKE algorithm [12] (tolerance = 0.0005 Å) to constrain covalent bonds involving hydrogens, using periodic boundary conditions, a 2 fs time step, a temperature coupling [13], a 8  cutoff was applied to the Lennard-Jones interaction, and constant pressure of 1 atm. The nonbonded list was updated every 10 steps. This preliminary study was aimed at finding the most energetically favorable structure which was found at 4856.9 ps for the trisaccharide and at 2411.8 ps for the pentasaccharide. Those conformations were found with a “home-made” Perl script that merely checked all energies values, finding the lowest one and giving the associated structure [2].

The PTRAJ module of AMBER10 was used to analyze and process the trajectory and coordinate files from SANDER, where the analyses might include carrying out superimposition, clustering analysis of hydrogen bonds, calculating fluctuations in bonds, angles, or dihedrals, correlation functions, and so forth, and the root mean-square deviations (RMSDs).

Finally, isoenergy contour maps were generated with ORIGIN and SURFER 9 (Golden Software, Inc., Golden Co) [14]. The positions of the minima in these maps were subsequently calculated by minimizing the energy without restraints of those grid points indicating minima on the maps.

3. Results

The calculated adiabatic conformational energy surfaces (φ, ψ) for the tri- and pentasaccharide Lewis X are presented as potential energy contour maps in Figures 3 and 4, respectively. The contours are separated by 1 kcal/mol and range from 1 to 10 kcal/mol. The energies are given relative to the lowest minimum. Geometries and relative potential energies of the predicted minima are listed in Tables 1 and 2, and the lowest energy conformers for each compound are shown in Figures 3 and 4.

tab1
Table 1: Potential Energy, torsional angles for the conformers with low energies of the trisaccharide Lewis X.
tab2
Table 2: Potential energy, torsional angles for the conformers with low energies of the pentasaccharide Lewis X.
725271.fig.003
Figure 3: Relaxed conformational energy maps for the trisaccharide Lewis X calculated using the GLYCAM 06 force field.
725271.fig.004
Figure 4: Relaxed conformational energy maps for the pentasaccharide Lewis X calculated using the GLYCAM 06 force field.
3.1. System Root-Mean-Square Displacement

RMSD is an important parameter for describing the system stability. Figures 5 and 6 show the function curves between mass weighted RMSD fit and the time of the tri- and pentasaccharide Lewis X during dynamics simulation process, respectively.

725271.fig.005
Figure 5: RMSD of the trisaccharide Lewis X as a function of time during MD, compared to the lowest-energy conformation.
725271.fig.006
Figure 6: RMSD of the pentasaccharide Lewis X as a function of time during MD, compared to the lowest-energy conformation.
3.2. Computation of Phi/Psi for Modeled Oligosaccharide

In the present study, structure prediction of oligosaccharides has been made through the torsion angles phi (ϕ)/psi ( ). The combination of ϕ and angles fully determine the conformation of an oligosaccharide. Figures 3 and 4 show the distribution of sterically allowed and energy minimized conformations in ϕ/ψ glycosidic space obtained from our explicit solvent MD simulations for 10 ns in the case of the tri- and pentasaccharide Lewis X, respectively. These figures also show the best conformers (which have the low energies). These conformers and their corresponding ϕ/ angular values are reported in Tables 1 and 2 for the tri- and pentasaccharide Lewis X, respectively.

4. Discussion

The monitoring of kinetic, global, and potential energies along the trajectory, as well as the pressure and temperature, indicate that the global periodic system is stable and does not present simulation artefacts. Kinetic energy increased slowly, whereas potential energy decreased gradually during the simulations. The total energy of a system to be simulated should be conserved, whereas the fluctuation of a temperature is caused by the interconvention of the kinetic and potential energy components [15].

The RMS variations as a function of time are displayed on Figures 5 and 6 for the tri- and penta-saccharide Lewis X respectively, and indicate that both systems are equilibrated after 500 ps. The visual examination of the two 10 ns trajectories shows that both oligosaccharides are stable. An RMSD analysis was carried out to quantify the structural differences. The RMSD between the PDB starting structure of the simulations (obtained from Glycam biomolecule builder), and the low-energy conformers structure considering only the heavy atoms is 0.03 Å and 0.20 Å for the tri- and penta-saccharide respectively. This shows that both oligosaccharides are stable.

During the simulation of those oligosaccharides, we notice that for the pentasaccharide the most stable conformer is obtained at 2411.8 ps and that for the trisaccharide was obtained at 4856.9 ps. This means that the pentasaccharide is less flexible than the trisaccharide due to the steric effect.

5. Conclusion

A complete understanding of the role of carbohydrates in biological systems is to a large extent dependent on the information available about the equilibrium mixture and about the preferred conformation of the carbohydrate molecules in solution. The conformational analysis offers a tool which can determine all possible conformations which influence the solution behavior of carbohydrates. This work represents a theoretical study of tow oligosaccharides tri- and pentasaccharide Lewis X using molecular dynamic simulations in explicit water with the AMBER10 program package. RMSD measures of the accuracy of carbohydrates around their average conformations play an important indicator of many biological processes such as macromolecular recognition and complex formations [16].

Acknowledgments

The authors are grateful to Dr. R. Borsali, Dr. K. Mazeau, Dr. A. Rivet, and Dr. A. Nurisso (CERMAV) France, and they thank particularly Mr. Guendouzi A for his help.

References

  1. A. Imberty and S. Perez, “Structure, conformation, and dynamics of bioactive oligosaccharides: theoretical approaches and experimental validations,” Chemical Reviews, vol. 100, no. 12, pp. 4567–4588, 2000. View at Publisher · View at Google Scholar · View at Scopus
  2. Y. Luo, F. Barbault, C. Gourmala et al., “Cellular interaction through LewisX cluster: theoretical studies,” Journal of Molecular Modeling, vol. 14, no. 10, pp. 901–910, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  3. E. A. Kabat, “Immunochemical studies on the carbohydrate moiety of water-soluble blood group A, B, H, Lea, and Leb substances and their precursor I antigens,” in Carbohydrates in Solution, H. Isbell, Ed., vol. 117, pp. 334–355, American Chemical Society, Washington, DC, USA, 1973.
  4. W. M. Watkins, “Biochemistry and genetics of the ABO, Lewis, and P blood group systems,” in Advances in Human Genetics, H. Harris and K. Hirsch-Horn, Eds., vol. 10, pp. 1–136, Plenum Press, New York, NY, USA, 10th edition, 1980.
  5. S. I. Hakomori, “Antigen structure and genetic basis of histo-blood groups A, B and O: their changes associated with human cancer,” Biochimica et Biophysica Acta, vol. 1473, no. 1, pp. 247–266, 1999. View at Publisher · View at Google Scholar · View at Scopus
  6. F. Bizik and I. Tvaroska, “Conformational analysis of disaccharide fragments of blood group determinants in solution by molecular modelling,” Chemical Papers, vol. 49, no. 4, pp. 202–214, 1995.
  7. A. Imberty, E. Mikros, J. Koca, R. Mollicone, R. Oriol, and S. Perez, “Computer simulation of histo-blood group oligosaccharides: energy maps of all constituting disaccharides and potential energy surfaces of 14 ABH and Lewis carbohydrate antigens,” Glycoconjugate Journal, vol. 12, no. 3, pp. 331–349, 1995. View at Publisher · View at Google Scholar · View at Scopus
  8. T. Kozar and C.-W. von der Lieth, “Efficient modelling protocols for oligosaccharides: from vacuum to solvent,” Glycoconjugate Journal, vol. 14, no. 8, pp. 925–933, 1997. View at Publisher · View at Google Scholar · View at Scopus
  9. D. Case, T. Darden, T. Cheatham, et al., “AMBER,” University of California, CA, USA, 10th edition, 2006.
  10. R. J. Woods, R. A. Dwek, C. J. Edge, and B. Fraser-Reid, “Glycam biomolecule builder,” http://www.glycam.com/CCRC/biombuilder/biomb_index.jsp.
  11. C. Gourmala, Y. Luo, F. Barbault et al., “Elucidation of the LewisX-LewisX carbohydrate interaction with molecular dynamics simulations: a glycosynapse model,” Journal of Molecular Structure, vol. 821, no. 1–3, pp. 22–29, 2007. View at Publisher · View at Google Scholar · View at Scopus
  12. J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, “Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes,” Journal of Computational Physics, vol. 23, no. 3, pp. 327–341, 1977. View at Publisher · View at Google Scholar · View at Scopus
  13. H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. Dinola, and J. R. Haak, “Molecular dynamics with coupling to an external bath,” The Journal of Chemical Physics, vol. 81, no. 8, article 3684, 7 pages, 1984. View at Publisher · View at Google Scholar · View at Scopus
  14. E. Y. Yoon, “A conformational study of linkage positions in oligosaccharides investigated by 2-D NMR spectroscopy and molecular modeling,” Bulletin of the Korean Chemical Society, vol. 24, no. 3, pp. 339–344, 2003.
  15. J. A. McCammon and S. C. Harvey, Dynamics of Proteins and Nucleic Acids, Cambridge University Press, Cambridge, UK, 1987.
  16. E. B. Philip and W. Helge, Structural Boinformatics, 1st edition, 2003.