Exploration of Novel Inhibitors for Class I Histone Deacetylase Isoforms by QSAR Modeling and Molecular Dynamics Simulation Assays

<div><p>Histone deacetylases (HDAC) are metal-dependent enzymes and considered as important targets for cell functioning. Particularly, higher expression of class I HDACs is common in the onset of multiple malignancies which results in deregulation of many target genes involved in cell growth, differentiation and survival. Although substantial attempts have been made to control the irregular functioning of HDACs by employing various inhibitors with high sensitivity towards transformed cells, limited success has been achieved in epigenetic cancer therapy. Here in this study, we used ligand-based pharmacophore and 2-dimensional quantitative structure activity relationship (QSAR) modeling approaches for targeting class I HDAC isoforms. Pharmacophore models were generated by taking into account the known IC₅₀ values and experimental energy scores with extensive validations. The QSAR model having an external R² value of 0.93 was employed for virtual screening of compound libraries. 10 potential lead compounds (C1-C10) were short-listed having strong binding affinities for HDACs, out of which 2 compounds (C8 and C9) were able to interact with all members of class I HDACs. The potential binding modes of HDAC2 and HDAC8 to C8 were explored through molecular dynamics simulations. Overall, bioactivity and ligand efficiency (binding energy/non-hydrogen atoms) profiles suggested that proposed hits may be more effective inhibitors for cancer therapy.</p></div

Classification and biological roles of HDACs.

<p>Classification and biological roles of HDACs.</p

Ligand efficiency (LE), Lipophilic Efficiency (LipE) and Fit Quality (FQ) values of C1-C10 hits.

<p>Ligand efficiency (LE), Lipophilic Efficiency (LipE) and Fit Quality (FQ) values of C1-C10 hits.</p

Statistical validation of designed pharmacophore and QSAR models.

<p><b>(a)</b> Top five pharmacophore models labeled as Model (1–5). <b>(b)</b> Correlation analysis of QSAR regression model. Green dots depict the observed vs predicted IC₅₀ values and correlation best fit dotted line is shown in blue.</p

Hydrogen bonds plot for 12 ns MD simulation.

<p>(a) Hydrogen bonds for HDAC2-C8 complex (yellow). (b) Hydrogen bonds for HDAC8-C8 complex (green).</p

2D structures and binding energies of compounds C1-C10.

<p><b>(a)</b> 2D structures <b>(b)</b> Binding energy overview.</p

Binding pattern of compound C8 with class 1 HDACs.

<p>C8 (7-methoxy-N-((4-sulfamoylphenyl)methyl)-1-benzofuran-2-carboxamide) was shown in blue, whereas hydrogen bonding and hydrophobic residues were shown in pink and yellow, respectively. <b>(a)</b> Hit C8 forms hydrogen bonds with HIS140, HIS148, HIS178, ASP264 and TYR303 of HDAC1; <b>(b)</b> HIS145, HIS146, HIS183 and TYR308 of HDAC2; <b>(c)</b> HIS134, HIS135, HIS172, ASP259 and <b>(d)</b>TYR298 of HDAC3 and HIS142, HIS143, HIS180 and TYR306 of HDAC8. C8 bonding with Zn+2 is shown in green color with an average distance of 2Å. Hydrophobic residues involved in interaction are <b>(a)</b> ASP99, GLY149, PHE150 and GLY301 in HDAC1; <b>(b)</b> PHE155, PHE210, ASP269 and GLY306 in HDAC2; <b>(c)</b> ASP93, PHE144, ASP170, PHE209 and GLY296 in HDAC3 and <b>(d)</b> PHE152, PHE208, MET274 and GLY304 in HDAC8.</p

Radius of gyration (Rg) analysis for 12 ns MD simulation.

<p>(a) Rg/RMSD plot for HDAC2-C8 complex (yellow) and HDAC2 without ligand (red). (b) Rg/RMSD for HDAC8-C8 complex (green) and HDAC8 without ligand (blue).</p

Energy plot for 12 ns MD simulation.

<p>(a) Energy plot for HDAC2-C8 complex (yellow). (b) Energy plot for HDAC8-C8 complex (green).</p

Ligand efficiency analysis of 10 selected compounds (C1-C10). (a) ClogP and (b) Fit Quality graph.

<p>Ligand efficiency analysis of 10 selected compounds (C1-C10). (a) ClogP and (b) Fit Quality graph.</p