Investigations on the binding specificity of β-galactoside analogues with human galectin-1 using molecular dynamics simulations

Galectin-1 (Gal-1) is the first member of galectin family, which has a carbohydrate recognition domain, specifically binds towards b-galactoside containing oligosaccharides. Owing its association with carbohydrates, Gal-1 is involved in many biological processes such as cell signaling, adhesion and pathological pathways such as metastasis, apoptosis and increased tumour cell survival. The development of b-galactoside based inhibitors would help to control the Gal-1 expression. In the current study, we carried out molecular dynamics (MD) simulations to examine the structural and dynamic behaviour Gal-1thiodigalactoside (TDG), Gal-1-lactobionic acid (LBA) and Gal-1-beta-(1!6)-galactobiose (G16G) complexes. The analysis of glycosidic torsional angles revealed that b-galactoside analogues TDG and LBA have a single binding mode (BM1) whereas G16G has two binding modes (BM1 and BM2) for interacting with Gal-1 protein. We have computed the binding free energies for the complexes Gal-1-TDG, Gal-1-LBA and Gal-1-G16G using MM/PBSA and are -6.45, -6.22 and -3.08 kcal/mol, respectively. This trend agrees well with experiments that the binding of Gal-1 with TDG is stronger than LBA. Further analysis revealed that the interactions due to direct and water-mediated hydrogen bonds play a significant role to the structural stability of the complexes. The result obtained from this study is useful to formulate a set of rules and derive pharmacophore-based features for designing inhibitors against galectin-1