Molecular modeling study on tunnel behavior in different histone deacetylase isoforms

Histone deacetylases (HDACs) have emerged as effective therapeutic targets in the treatment of various diseases including cancers as these enzymes directly involved in the epigenetic regulation of genes. However the development of isoform-selective HDAC inhibitors has been a challenge till date since all HDAC enzymes possess conserved tunnel-like active site. In this study, using molecular dynamics simulation we have analyzed the behavior of tunnels present in HDAC8, 10, and 11 enzymes of class I, II, and IV, respectively. We have identified the equivalent tunnel forming amino acids in these three isoforms and found that they are very much conserved with subtle differences to be utilized in selective inhibitor development. One amino acid, methionine of HDAC8, among six tunnel forming residues is different in isoforms of other classes (glutamic acid (E) in HDAC10 and leucine (L) in HDAC 11) based on which mutations were introduced in HDAC11, the less studied HDAC isoform, to observe the effects of this change. The HDAC8-like (L268M) mutation in the tunnel forming residues has almost maintained the deep and narrow tunnel as present in HDAC8 whereas HDAC10-like (L268E) mutation has changed the tunnel wider and shallow as observed in HDAC10. These results explained the importance of the single change in the tunnel formation in different isoforms. The observations from this study can be utilized in the development of isoform-selective HDAC inhibitors.