Molecular dynamics simulations of asialoglycoprotein receptor ligands

Several recent studies have implicated carbohydrates in cell adhesion, inflammation, clearance of glycoproteins from blood circulation, embryonic development, and metastasis among others. Understanding the conformation of these carbohydrate recognition elements and their interaction at the molecular level is essential for the design of oligosaccharide inhibitors/drugs. Given the difficulty in solving carbohydrate structures by X-ray crystallography and since NMR experiments give only time-averaged conformation, molecular dynamics simulations are well suited to determine all the accessible conformations of oligosaccharides. Present communication reports the simulation of some of the oligosaccharide ligands of asialoglycoprotein receptor for 1 ns using Biosym's Insight II molecular modeling package on NCI-FCRDC's Y-MP 8D/8128 supercomputer. Results obtained from these simulations, in addition to explaining the observed differences in the binding affinities of these ligands to the asialoglycoprotein receptor, have led to a modified model for the recognition of the oligosaccharides by the receptor. Accordingly, only the two terminal galactose residues on the 1,3-arm of the triantennary oligosaccharide (GlcNAc2Man3 core of the N-linked oligosaccharides with N-acetyllactosamine in beta 1,2- and beta 1,4-linkages on the 1,3-linked core mannose) are primarily required for recognition, and the terminal galactose on the 1,6-arm (N-acetyllactosamine in beta 1,2-linkage on the 1,6-linked core mannose) provides additional binding energy. It has been shown that the oligosaccharides studied here have significant flexibility and the flexibility is more around the 1,3-linkage than the 1,6-linkage. The need for simulation for longer periods and with multiple initial conformations is also discussed in the present report.