Molecular docking studies of a group of hydroxamate inhibitors with gelatinase-A by molecular dynamics

We have performed docking and molecular dynamics simulations of hydroxamates complexed with human gelatinase-A (MMP-2) to gain insight into the structural and energetic preferences of these inhibitors. The study was conducted on a selected set of eleven compounds with variation in structure and activity. Molecular dynamics simulations were performed at 300 K for 100 ps with equilibration for 50 ps. The structural analyses of the trajectories indicate that the coordinate bond interactions, the hydrogen bond interactions, the van der Waals interactions as well as the hydrophobic interactions between ligand and receptor are responsible simultaneously for the preference of inhibition and potency. The ligand hydroxamate group is coordinated to the catalytic zinc ion and form stable hydrogen bonds with the carbonyl oxygen of Gly 162. The P1 group makes extensive van der Waals and hydrophobic contacts with the nonpolar side chains of several residues in the S1 subsite, including Leu 197, Val 198, Leu 218 and Tyr 223. Moreover, four to eight hydrogen bonds between hydroxamates and MMP-2 are formed to stabilize the inhibitors in the active site. Compared with the P2 and P3 groups, the P1 groups of inhibitors are oriented regularly, which is produced by the restrain of the S1 subsite. From the relationship between the length of the nonpolar P1 group and the biological activity, we confirm that MMP-2 has a pocket-like S1 subsite, not a channel-like S1 subsite proposed by Kiyama (Kiyama, R. et al., J. Med. Chem. 42 (1999), 1723). The energetic analyses show that the experimental binding free energies can be well correlated with the interactions between the inhibitors and their environments, which could be used as a simple score function to evaluate the binding affinities for other similar hydroxamates. The validity of the force field parameters and the MD simulations can be fully testified by the satisfactory agreements between the experimental structure-activity relationship and the information from the structural and energetic analyses. The information generated from the predicted complexes should be useful for further work in the area of structure-based design of new compounds.