Discovery of novel hit compounds as potential HDAC1 inhibitors: The case of ligand- and structure-based virtual screening

Histone deacetylases (HDACs) as an important family of epigenetic regulatory enzymes are implicated in the onset and progression of carcinomas. As a result, HDAC inhibition has been proven as a compelling strategy for reversing the aberrant epigenetic changes associated with cancer. However, non-selective profile of most developed HDAC inhibitors (HDACIs) leads to the occurrence of various side effects, limiting their clinical utility. This evidence provides a solid ground for ongoing research aimed at identifying isoform-selective inhibitors. Among the isoforms, HDAC1 have particularly gained increased attention as a preferred target for the design of selective HDACIs. Accordingly, in this paper, we have developed a reliable virtual screening process, combining different ligand- and structure–based methods, to identify novel benzamide-based analogs with potential HDAC1 inhibitory activity. For this purpose, a focused library of 736,160 compounds from PubChem database was first compiled based on 80% structural similarity with four known benzamide-based HDAC1 inhibitors, Mocetinostat, Entinostat, Tacedinaline, and Chidamide. Our inclusive in-house 3D-QSAR model, derived from pharmacophore-based alignment, was then employed as a 3D-query to discriminate hits with the highest predicted HDAC1 inhibitory activity. The selected hits were subjected to subsequent structure-based approaches (induced-fit docking (IFD), MM-GBSA calculations and molecular dynamics (MD) simulation) to retrieve potential compounds with the highest binding affinity for HDAC1 active site. Additionally, in silico ADMET properties and PAINS filtration were also considered for selecting an enriched set of the best drug-like molecules. Finally, six top-ranked hit molecules, CID_38265326, CID_56064109, CID_8136932, CID_55802151, CID_133901641 and CID_18150975 were identified to expose the best stability profiles and binding mode in the HDAC1 active site. The IFD and MD results cooperatively confirmed the interactions of the promising selected hits with critical residues within HDAC1 active site. In summary, the presented computational approach can provide a set of guidelines for the further development of improved benzamide-based derivatives targeting HDAC1 isoform.