Discovery of selective HDAC6 inhibitors driven by artificial intelligence and molecular dynamics simulation approaches

Increasing evidence showed that histone deacetylase 6 (HDAC6) dysfunction is directly associated with the onset and progression of various diseases, especially cancers, making the development of HDAC6-targeted anti-tumor agents a research hotspot. In this study, artificial intelligence (AI) technology and molecular simulation strategies were fully integrated to construct an efficient and precise drug screening pipeline, which combined Voting strategy based on compound-protein interaction (CPI) prediction models, cascade molecular docking, and molecular dynamic (MD) simulations. The biological potential of the screened compounds was further evaluated through enzymatic and cellular activity assays. Among the identified compounds, Cmpd.18 exhibited more potent HDAC6 enzyme inhibitory activity (IC₅₀ = 5.41 nM) than that of tubastatin A (TubA) (IC₅₀ = 15.11 nM), along with a favorable subtype selectivity profile (selectivity index ≈ 117.23 for HDAC1), which was further verified by the Western blot analysis. Additionally, Cmpd.18 induced G2/M phase arrest and promoted apoptosis in HCT-116 cells, exerting desirable antiproliferative activity (IC₅₀ = 2.59 μM). Furthermore, based on long-term MD simulation trajectory, the key residues facilitating Cmpd.18's binding were identified by decomposition free energy analysis, thereby elucidating its binding mechanism. Moreover, the representative conformation analysis also indicated that Cmpd.18 could stably bind to the active pocket in an effective conformation, thus demonstrating the potential for in-depth research of the 2-(2-phenoxyethyl)pyridazin-3(2H)-one scaffold.