Abstract
BACKGROUND: Computational drug design approaches are important for shortening the time and reducing the cost for drug discovery and development. Among these methods, molecular docking and quantitative structure activity relationship (QSAR) play key roles for lead discovery and optimization. Here, we propose an integrated approach with core strategies to identify the protein-ligand hot spots for QSAR models and lead optimization. These core strategies are: 1) to generate both residue-based and atom-based interactions as the features; 2) to identify compound common and specific skeletons; and 3) to infer consensus features for QSAR models. RESULTS: We evaluated our methods and new strategies on building QSAR models of human acetylcholinesterase (huAChE). The leave-one-out cross validation values q (2) and r (2) of our huAChE QSAR model are 0.82 and 0.78, respectively. The experimental results show that the selected features (resides/atoms) are important for enzymatic functions and stabling the protein structure by forming key interactions (e.g., stack forces and hydrogen bonds) between huAChE and its inhibitors. Finally, we applied our methods to arthrobacter globiformis histamine oxidase (AGHO) which is correlated to heart failure and diabetic. CONCLUSIONS: Based on our AGHO QSAR model, we identified a new substrate verified by bioassay experiments for AGHO. These results show that our methods and new strategies can yield stable and high accuracy QSAR models. We believe that our methods and strategies are useful for discovering new leads and guiding lead optimization in drug discovery.