Energy based pharmacophore mapping of HDAC inhibitors against class I HDAC enzymes

Histone deacetylases (HDACs) are important class of enzymes that deacetylate the ε-amino group of the lysine residues in the histone tails to form a closed chromatin configuration resulting in the regulation of gene expression. Inhibition of these HDACs enzymes have been identified as one of the promising approaches for cancer treatment. The type-specific inhibition of class I HDAC enzymes is known to elicit improved therapeutic effects and thus, the search for promising type-specific HDAC inhibitors compounds remains an ongoing research interest in cancer drug discovery. Several different strategies are employed to identify the features that could identify the isoform specificity factors in these HDAC enzymes. This study combines the insilico docking and energy-optimized pharmacophore (e-pharmacophore) mapping of several known HDACi's to identify the structural variants that are significant for the interactions against each of the four class I HDAC enzymes. Our hybrid approach shows that all the inhibitors with at least one aromatic ring in their linker regions hold higher affinities against the target enzymes, while those without any aromatic rings remain as poor binders. We hypothesize the e-pharmacophore models for the HDACi's against all the four Class I HDAC enzymes which are not reported elsewhere. The results from this work will be useful in the rational design and virtual screening of more isoform specific HDACi's against the class I HDAC family of proteins.