Combining in silico tools and NMR data to validate protein-ligand structural models: Application to matrix metalloproteinases

A combination of in silico tools and experimental NMR data is proposed for relatively fast determination of protein-ligand structural models and demonstrated from known inhibitors of matrix metalloproteinases (MMP). The 15N 1H heteronuclear single quantum coherence (HSQC) spectral assignment and the 3D structure, either X-ray or NMR, are needed. In this method, the HSQC spectrum with or without the ligand is used to determine the interaction region of the ligand. Docking calculations are then performed to obtain a set of structural models. From the latter, the nuclear Overhauser effects (NOEs) between the ligand and the protein can be predicted. Guided by these predictions, a number of NOEs can be detected and assigned through a HSQC NOESY experiment. These data are used as structural restraints to reject/refine the initial structural models through further in silico work. For a test protein (MMP-12, human macrophage metalloelastase), a final structure of a protein-ligand adduct was obtained that matches well with the full structural determination. A number of structural predictions were then made for adducts of a similar protein (MMP-1, human fibroblast collagenase) with the same and different ligands. The quality of the final results depended on the type and number of experimental NOEs, but in all cases, a well-defined ligand conformation in the protein binding site was obtained. This protocol is proposed as a viable alternative to the many approaches described in the literature.