Energy-optimised pharmacophore approach to identify potential hotspots during inhibition of Class II HDAC isoforms

Histone deacetylases (HDACs) are conjugated enzymes that modulate chromatin architecture by deacetylating lysine residues on the histone tails leading to transcriptional repression. Pharmacological interventions of these enzymes with small molecule inhibitors called Histone deacetylase inhibitors (HDACi) have shown enhanced acetylation of the genome and are hence emerging as potential targets at the clinic. Type-specific inhibition of Class II HDACs has shown enhanced therapeutic benefits against developmental and neurodegenerative disorders. However, the structural identity of class-specific isoforms limits the potential of their inhibitors in precise targeting of their enzymes. Diverse strategies have been implemented to recognise the features in HDAC enzymes which may help in identifying isoform specificity factors. This work attempts a computational approach that combines in silico docking and energy-optimised pharmacophore (E-pharmacophore) mapping of 18 known HDAC inhibitors and has identified structural variations that regulate their interactions against the six Class II HDAC enzymes considered for the study. This combined approach establishes that inhibitors possessing higher number of aromatic rings in different structural regions might function as potent inhibitors, while inhibitors with scarce ring structures might point to compromised potency. This would aid the rationale for chemical optimisation and design of isoform selective HDAC inhibitors with enhanced affinity and therapeutic efficiency.