Identification of novel HDAC8 selective inhibitors through ligand and structure based studies: Exploiting the acetate release channel differences among class I isoforms

Histone deacetylases (HDACs) are key regulators of gene expression and have emerged as crucial therapeutic targets for cancer. Among the HDACs, inhibition of HDAC8 enzyme has been reported to be a novel strategy in the treatment of female-specific cancers. Most of the HDAC inhibitors discovered so far inhibit multiple HDAC isoforms causing toxicities in the clinic thus limiting their potential. Therefore, the discovery of isoform-selective HDAC8 inhibitors is highly desirable. In the present study, a combination of ligand and structure based drug design tools were utilized to build a statistically significant pharmacophore based 3D QSAR model with statistical parameters R²: 0.9964, and Q²: 0.7154, from a series of 31 known HDAC8 inhibitors. Top 1000 hits obtained from Virtual screening of Phase database were subjected to docking studies against HDAC8. Top 100 hits obtained were redocked into HDAC Class I (HDAC 1,2,3) and Class II isoforms (HDAC 4, 6) and rescored with XP Glide Score. Based on fitness score, XP glide score and interacting amino acid residues, five HDAC8 inhibitors (1–5) were selected for in vitro studies. The HDAC8 activity assay followed by enzyme kinetics clearly indicated Compounds 1, 2 and 3 to be potent HDAC8 selective inhibitors with IC₅₀ of 126 pM, 112 nM, and 442 nM respectively. These compounds were cytotoxic to HeLa cells where HDAC8 is overexpressed but not to normal cells, HEK293. Also, they were able to induce apoptosis by modulating Bax/Bcl2, cleavage of PARP and release of Cytochrome C. Molecular Dynamics simulations observed most favorable interaction patterns and presented a rationale for the activities of the identified compounds. Selectivity against HDAC8 was due to exploitation of the architectural difference in the acetate release channel among class I HDAC isoforms.