Molecular modelling of protein-carbohydrate interactions. Docking of monosaccharides in the binding site of concanavalin A

A general procedure is described for addressing the computer simulation of protein-carbohydrate interactions. First, a molecular mechanical force field capable of performing conformational analysis of oligosaccharides has been derived by the addition of new parameters to the Tripos force field; it is also compatible with the simulation of protein. Second, a docking procedure which allows for a systematic exploration of the orientations and positions of a ligand into a protein cavity has been designed. This so-called ‘crankshaft’ method uses rotations and variations about/of virtual bonds connecting, via dummy atoms, the ligand to the protein binding site. Third, calculation of the relative stability of protein ligand complexes is performed. This strategy has been applied to search for all favourable interactions occurring between a lectin [concanavalin A (ConA)] and methyl a-D-mannopyranoside or methyl α-D-glucopyranoside. For each monosaccharide, different stable orientations and positions within the binding site can be distinguished. Among them, one corresponds to very favourable interactions, not only in terms of hydrogen bonding, but also in terms of van der Waals interactions. It corresponds precisely to the binding mode of methyl α-D-mannopyranoside into ConA as revealed by the 2.9 Ä resolution of the crystalline complex (Derewenda et al., 1989). Some implications of the present modelling study with respect to the molecular basis of the specificity of the interaction of lectins with various monosaccharides are presented.