Quinazolinone–hydroxamate derivatives as potential HDAC8 inhibitors: An integrated docking, ADME, DFT, and molecular dynamics study

Abstract

Cancer remains a major global health challenge, highlighting the need for novel therapeutic agents with improved selectivity and efficacy. Histone deacetylase (HDAC) inhibition represents an important epigenetic strategy in anticancer drug discovery because it can influence gene expression, cell-cycle regulation, and tumor cell survival. In this computational study, fifteen 4-oxoquinazolin-3(4 H )-yl- N -hydroxybenzamide derivatives were evaluated as potential HDAC8-targeting compounds using an integrated in silico workflow. Molecular docking was performed against HDAC8 using the crystal structure deposited under Protein Data Bank (PDB) ID 1T69, with suberoylanilide hydroxamic acid (SAHA) used as the reference inhibitor. The investigated derivatives showed favorable accommodation within the HDAC8 catalytic pocket, with predicted interaction patterns involving zinc coordination, hydrogen bonding, and π-related interactions with key active-site residues. Several derivatives showed favorable predicted docking scores and interaction profiles relative to SAHA, suggesting that the designed quinazolinone–hydroxamate scaffold may be suitable for further prioritization. Prime molecular mechanics generalized Born surface area (MM-GBSA) analysis identified selected derivatives with comparatively favorable calculated binding free-energy contributions, particularly through van der Waals, electrostatic, and lipophilic interactions. In silico absorption, distribution, metabolism, and excretion (ADME) prediction indicated that most derivatives possessed acceptable drug-like and pharmacokinetic descriptors, including physicochemical properties, predicted oral absorption, and limited central nervous system (CNS) activity. Density functional theory (DFT) analysis provided supportive electronic-level information by identifying molecular regions that may contribute to hydrogen bonding and zinc coordination. A 100 ns molecular dynamics simulation of the HDAC8–derivative 15 complex suggested the short-timescale persistence of one representative predicted binding pose under the selected simulation conditions. Overall, this study identifies selected 4-oxoquinazoline hydroxamic acid derivatives as preliminary computational leads for further experimental validation as potential HDAC8-targeting compounds.