Pharmmaker: Pharmacophore modeling and hit identification based on druggability simulations

Recent years have seen progress in druggability simulations, that is, molecular dynamics simulations of target proteins in solutions containing drug‐like probe molecules to characterize their drug‐binding abilities, if any. An important consecutive step is to analyze the trajectories to construct pharmacophore models (PMs) to use for virtual screening of libraries of small molecules. While considerable success has been observed in this type of computer‐aided drug discovery, a systematic tool encompassing multiple steps from druggability simulations to pharmacophore modeling, to identifying hits by virtual screening of libraries of compounds, has been lacking. We address this need here by developing a new tool, Pharmmaker, building on the DruGUI module of our ProDy application programming interface. Pharmmaker is composed of a suite of steps: (Step 1) identification of high affinity residues for each probe molecule type; (Step 2) selecting high affinity residues and hot spots in the vicinity of sites identified by DruGUI; (Step 3) ranking of the interactions between high affinity residues and specific probes; (Step 4) obtaining probe binding poses and corresponding protein conformations by collecting top‐ranked snapshots; and (Step 5) using those snapshots for constructing PMs. The PMs are then used as filters for identifying hits in structure‐based virtual screening. Pharmmaker, accessible online at http://prody.csb.pitt.edu/pharmmaker/, can be used in conjunction with other tools available in ProDy.