Virtual Screening and Scaffold Hopping Based on GRID Molecular Interaction Fields

Abstract

In this study, a set of strategies for structure-based design using GRID molecular interaction fields (MIFs) to derive a pharmacophoric representation of a protein is reported. Thrombin, one of the key enzymes involved in the blood coagulation cascade, was chosen as the model system since abundant published experimental data are available related to both crystal structures and structurally diverse sets of inhibitors. First, a virtual screening methodology was developed either using a pharmacophore representation of the protein based on GRID MIFs or using GRID MIFs from the 3D structure of a set of chosen thrombin inhibitors. The search was done in a 3D multiconformation version of the Available Chemical Directory (ACD) database, which had been spiked with 262 known thrombin inhibitors (multiple conformers available per compound). The model managed to find 80% of the known thrombin inhibitors among the 74,291 conformers in the ACD by only searching 5% of the database; hence, a 15-fold enrichment of the library was achieved. Second, a scaffold hopping methodology was developed using GRID MIFs, giving the scaffold interaction pattern and the shape of the scaffold, together with the distance between the anchor points. The scaffolds reported by Dolle in the Journal of Combinatorial Chemistry summaries (2000 and 2001) and scaffolds built or derived from ligands cocomplexed with the thrombin enzyme were parameterized using a new set of descriptors and saved into a searchable database. The scaffold representation from the database was then compared to a template scaffold (from a thrombin crystal structure), and the thrombin-derived scaffolds included in the database were found among the top solutions. To validate the usefulness of the methodology to replace the template scaffold, the entire molecule was built (scaffold and side chains) and the resulting compounds were docked into the active site of thrombin. The docking solutions showed the same binding pattern as the cocomplexed compound, hence, showing that this method can be a valuable tool for medicinal chemists to select interchangeable core structures (scaffolds) in an easy manner and retaining the binding properties from the original ligand.