Molecular dynamics simulations of alpha2 --> 8-linked disialoside: conformational analysis and implications for binding to proteins.
Computational methods have played a key role in elucidating the various three-dimensional structures of oligosaccharides. Such structural information, together with other experimental data, leads to a better understanding of the role of oligosaccharide in various biological processes. The disialoside Neu5Ac-alpha2-->8-Neu5Ac appears as the terminal glycan in glycoproteins and glycolipids, and is known to play an important role in various events of cellular communication. Neurotoxins such as botulinum and tetanus require Neu5Ac-alpha2 --> 8-Neu5Ac for infecting the host. Glycoconjugates containing this disialoside and the enzymes catalyzing their biosynthesis are also regulated during cell growth, development, and differentiation. Unlike other biologically relevant disaccharides that have only two linkage bonds, the alpha2-->8-linked disialoside has four: C2-O, O-C8', C8'-C7', and C7'-C6'. The present report describes the results from nine 1 ns MD simulations of alpha2-->8-linked disialoside (Neu5Ac-alpha2-->8-Neu5Ac); simulations were run using GROMOS96 by explicitly considering the solvent molecules. Conformations around the O-C8' bond are restricted to the +sc/+ap regions due to stereochemical reasons. In contrast, conformations around the C2-O and C8'-C7' bonds were found to be largely unrestricted and all the three staggered regions are accessible. The conformations around the C7'-C6' bond were found to be in either the -sc or the anti region. These results are in excellent agreement with the available NMR and potential energy calculation studies. Overall, the disaccharide is flexible and adopts mainly two ensembles of conformations differing in the conformation around the C7'-C6' bond. The flexibility associated with this disaccharide allows for better optimization of intermolecular contacts while binding to proteins and this may partially compensate for the loss of conformational entropy that may be incurred due to disaccharide's flexibility.
Vasudevan
Sheeja V
SV
Biotechnology Center, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India.
Balaji
Petety V
PV
D000465:Algorithms; D002236:Carbohydrate Conformation; D002240:Carbohydrate Sequence; D003198:Computer Simulation; D004187:Disaccharides / Q000737:chemistry* / Q000378:metabolism*; D005732:Gangliosides / Q000737:chemistry; D006860:Hydrogen Bonding; D008958:Models, Molecular; D008969:Molecular Sequence Data; D019158:N-Acetylneuraminic Acid / Q000737:chemistry* / Q000378:metabolism*; D009438:Neuraminic Acids / Q000737:chemistry* / Q000378:metabolism*; D011485:Protein Binding; D012984:Software; D012996:Solutions; D012997:Solvents; D013816:Thermodynamics; D014102:Torsion Abnormality
D004187:Disaccharides; D005732:Gangliosides; C453731:N-acetylneuraminyl-2-8-N-acetylneuraminic acid; D009438:Neuraminic Acids; D012996:Solutions; D012997:Solvents; D019158:N-Acetylneuraminic Acid