In silico identification of potential Hsp90 inhibitors via ensemble docking, DFT and molecular dynamics simulations

The molecular chaperone heat shock protein 90 (Hsp90) has emerged as one of the most exciting targets for anticancer drug development and Hsp90 inhibitors are potentially useful chemotherapeutic agents in cancer. Within the current study, Hsp90 inhibitors that entered different phases of clinical trials were subjected to Zinc15 structure query to find similar compounds (>= 78%). Obtained small molecules (1-29) with defined similarity cut-off were docked into ensemble of Hsp90-a NTDs. Docked complexes were ranked on the basis of binding modes and Gibbs free energies as Hsp90 binders (cutoff point; Delta G(b) <= -12 kcal/mol). Top-ranked compounds were subjected to energy decomposition analysis per residue of binding pocket via density functional theory (DFT) calculations in B3LYP level of theory. Subsequent MD simulations of the top-ranked complexes were performed for 100 ns to explore the stable binding modes during a reasonable period in explicit water. Results of molecular docking and intermolecular binding analysis indicated that H-bond, hydrophobic and salt bridge interactions were determinant forces in complex formation. Compounds 19 and 20 were well accommodated in binding pocket of Hsp90 via relatively varied conformations. It was revealed that Asn51 and Phe138 were key residues that interacted stably to 19 and 20. Although primary mechanism of action for proposed molecules are unknown and yet to be explored, results of the present study revealed key structural features for future structure-guided optimization toward potent inhibitors of Hsp90-alpha NTD.