Probing the mechanisms of hydrazide-based HDAC inhibitors binding to HDAC3 using Gaussian accelerated molecular dynamics (GaMD) simulations

Abstract

Histone deacetylases (HDACs) have emerged as promising targets for anticancer drug development. They regulate gene expression by removing acetyl groups from lysine residues on histone tails, leading to chromatin condensation. A hydrazide-based HDAC inhibitor, N-(4-(2-Propylhydrazine-1-carbonyl)benzyl)-1H-indole-2-carboxamide (11h), has been reported to exhibit significant in vivo antitumor activity. In comparison to the lead compound N-(4-(2-Propylhydrazine-1-carbonyl)benzyl)cinnamamide (17), compound 11h demonstrates 2- to 5-fold higher HDAC inhibition and cell-based antitumor activity. However, the inhibitory mechanism of 11h remains insufficiently explored. In this study, we conducted 500 ns Gaussian Accelerated Molecular Dynamics (GaMD) simulations on Histone deacetylase 3 (HDAC3) and two complex systems (HDAC3-17 and HDAC3-11h). Our findings revealed that upon inhibitor binding, the active pocket volume of HDAC3 undergone alterations, and the movement of the L6-loop toward the active site impeded substrate entry. Moreover, we observed a destabilization of the α-helix in the aa75-89 region of HDAC3 compared to the two complex systems, indicating partial unwinding. Notably, 11h exhibited a closer proximity of its carbonyl oxygen to the active pocket’s Zn2+ metal compared to 17, increasing the likelihood of coordination with the Zn2+ metal. The analysis of protein-ligand interactions highlighted a greater number of hydrogen bonds and other interactions between 11h and the receptor protein when compared to 17, underscoring the stronger binding of 11h to HDAC3. In conclusion, our study provided theoretical insights into the inhibitory mechanism of hydrazide-based HDAC inhibitors on HDAC3, thereby contributing to the development of improved drug targets for cancer therapy. Communicated by Ramaswamy H. Sarma