The combination of molecular dynamics with crystallography for elucidating protein-ligand interactions: a case study involving peanut lectin complexes with T-antigen and lactose.

Peanut lectin binds T-antigen [Galbeta(1-3)GalNAc] with an order of magnitude higher affinity than it binds the disaccharide lactose. The crystal structures of the two complexes indicate that the higher affinity for T-antigen is generated by two water bridges involving the acetamido group. Fresh calorimetric measurements on the two complexes have been carried out in the temperature range 280-313 K. Four sets of nanosecond molecular-dynamics (MD) simulations, two at 293 K and the other two at 313 K, were performed on each of the two complexes. At each temperature, two somewhat different protocols were used to hydrate the complex in the two runs. Two MD runs under slightly different conditions for each complex served to assess the reliability of the approach for exploring protein-ligand interactions. Enthalpies based on static calculations and on MD simulations favour complexation involving T-antigen. The simulations also brought to light ensembles of direct and water-mediated protein-sugar interactions in both the cases. These ensembles provide a qualitative explanation for the temperature dependence of the thermodynamic parameters of peanut lectin-T-antigen interaction and for the results of one of the two mutational studies on the lectin. They also support the earlier conclusion that the increased affinity of peanut lectin for T-antigen compared with that for lactose is primarily caused by additional water bridges involving the acetamido group. The calculations provide a rationale for the observed sugar-binding affinity of one of the two available mutants. Detailed examination of the calculations point to the need for exercising caution in interpreting results of MD simulations: while long simulations are not possible owing to computational reasons, it is desirable to carry out several short simulations with somewhat different initial conditions.

X-ray structure determination and modeling of the cyclic tetrasaccharide cyclo.

The cyclic tetrasaccharide cyclo-[-->6)-alpha-D-Glcp-(1-->3)-alpha-D-Glcp-(1-->6)-alpha-D-Glcp-(1-->3)-alpha-D-Glcp-(1-->] is the major compound obtained by the action of endo-alternases on the alternan polysaccharide. Crystals of this cyclo-tetra-glucose belong to the orthorhombic space group P2(1)2(1)2(1) with a = 7.620(5), b = 12.450(5) and c = 34.800(5) A. The asymmetric unit contains one tetrasaccharide together with five water molecules. The tetrasaccharide adopts a plate-like overall shape with a very shallow depression on one side. The shape is not fully symmetrical and this is clearly apparent on comparing the (phi, psi) torsion angles of the two alpha-(1-->6) linkages. There is almost 10 degrees differences in phi and more than 20 degrees differences in psi. The hydrogen bond network is asymmetric, with a single intramolecular hydrogen bond: O-2 of glucose ring 1 being the donor to O-2 of glucose ring 3. These two hydroxyl groups are located below the ring and their orientation, dictated by this hydrogen bond, makes the floor of the plate. Among the five water molecules, one located above the center of the plate occupies perfectly the shallow depression in the plate shape formed by the tetrasaccharide. Molecular dynamics simulation of the tetrasaccharide in explicit water allows rationalization of the discrepancies observed between the X-ray structures and data obtained previously by NMR.

Structure/thermodynamics relationships of lectin-saccharide complexes: the Erythrina corallodendron case.

Molecular dynamics (MD) simulations of Erythrina corallodendron lectin binding to a monosaccharide, alpha-galactose, and a disaccharide, N-acetyl lactosamine, have been performed in order to investigate the relationship between structure and thermodynamics. A simulated annealing protocol has been used to generate ensembles of structures for the two complexes, from which both qualitative and quantitative information on binding dynamics have been extracted. The ensembled averaged lectin-saccharide interaction enthalpy is equivalent for both sugars, whereas the calculation based on the X-ray structures does show a difference. Within large statistical errors, the calculated 'binding enthalpy' is also the same for the two systems. These errors arise largely from terms involving solvent and are a typical limitation of current MD simulations. Significant qualitative differences in binding between the two complexes are, however, observed over the ensembles. These could be important for unraveling the structure/thermodynamic relationship. Stated simply, there are a greater number of binding options available to the disaccharide compared to the monosaccharide. The implications of alternative binding states on thermodynamic parameters and the 'breaking of enthalpy-entropy compensation' are discussed. The role of solvent in lectin-saccharide complex formation is suggested to be significant.