Novel inhibitors of human histone deacetylase (HDAC) identified by QSAR modeling of known inhibitors, virtual screening, and experimental validation

Inhibitors of histone deacetylases (HDACIs) have emerged as a new class of drugs for the treatment of human cancers and other diseases because of their effects on cell growth, differentiation, and apoptosis. In this study we have developed several quantitative structure - activity relationship (QSAR) models for 59 chemically diverse histone deacetylase class 1 (HDAC1) inhibitors. The variable selection k nearest neighbor (kNN) and support vector machines (SVM) QSAR modeling approaches using both MolconnZ and MOE chemical descriptors generated from two-dimensional rendering of compounds as chemical graphs have been employed. We have relied on a rigorous model development workflow including the division of the data set into training, test, and external sets and extensive internal and external validation. Highly predictive QSAR models were generated with leave-one-out cross-validated (LOO-CV) q² and external R² values as high as 0.80 and 0.87, respectively, using the kNN/MolconnZ approach and 0.93 and 0.87, respectively, using the SVM/MolconnZ approach. All validated QSAR models were employed concurrently for virtual screening (VS) of an in-house compound collection including 9.5 million molecules compiled from the ZINC7.0 database, the World Drug Index (WDI) database, the ASINEX Synergy libraries, and other commercial databases. VS resulted in 45 structurally unique consensus hits that were considered novel putative HDAC1 inhibitors. These computational hits had several novel structural features that were not present in the original data set. Four computational hits with novel scaffolds were tested experimentally, and three of them were confirmed active against HDAC 1, with IC₅₀ values for the most active compound of 1.00 μM. The fourth compound was later identified to be a selective inhibitor of HDAC6, a Class II HDAC. Moreover, two of the confirmed hits are marketed drugs, which could potentially facilitate their further development as anticancer agents. This study illustrates the power of the combined QS AR-VS method as a general approach for the effective identification of structurally novel bioactive compounds.