Examining the stability of binding modes of the co-crystallized inhibitors of human HDAC8 by molecular dynamics simulation

In the quest for histone deacetylase inhibitors: current trends in the application of multilayered computational methods

Histone deacetylase (HDAC) inhibitors have gained attention over the past three decades because of their potential in the treatment of different diseases including various forms of cancers, neurodegenerative disorders, autoimmune, inflammatory diseases, and other metabolic disorders. To date, 5 HDAC inhibitor drugs are marketed for the treatment of hematological malignancies and several drug-candidate HDAC inhibitors are at different stages of clinical trials. However, due to the toxic side effects of these drugs resulting from the lack of target selectivity, active studies are ongoing to design and develop either class-selective or isoform-selective inhibitors. Computational methods have aided the discovery of HDAC inhibitors with the desired potency and/or selectivity. These methods include ligand-based approaches such as scaffold hopping, pharmacophore modeling, three-dimensional quantitative structure-activity relationships (3D-QSAR); and structure-based virtual screening (molecular docking). The current trends involve the application of the combination of these methods and incorporating molecular dynamics simulations coupled with Poisson-Boltzmann/molecular mechanics generalized Born surface area (MM-PBSA/MM-GBSA) to improve the prediction of ligand binding affinity. This review aimed at understanding the current trends in applying these multilayered strategies and their contribution to the design/identification of HDAC inhibitors.

Synthesis of N-aminophalimides derived from α-amino acids: Theoretical study to find them as HDAC8 inhibitors by docking simulations and in vitro assays

Phthalimides are valuable for synthesis and biological properties. New acetamides 3(a-c) and 4(a-c) were synthesized and characterized as precursors for novel N-aminophalimides 5(a-c) and 6(a-c). Structures of 4a, 5(a-b), and 6(a-b) were confirmed by single crystal X-ray. Docking studies identified compounds with favorable Gibbs free energy values for binding to histone deacetylase 8 (HDAC8), an enzyme targeted for anticancer drug development. These compounds bound to both the orthosteric and allosteric pockets of HDAC8, similar to Trichostatin A (TSA), an HDAC8 inhibitor. 6(a-c) contain hydroxyacetamide moiety as a zinc-binding group, a phthalimide moiety as a capping group, and aminoacetamide moiety as a linker group, which are important for ligand-receptor binding. ΔG values indicated that compounds 5b, 6b, and 6c had higher affinity for HDAC8 in the allosteric pocket compared to TSA. In vitro evaluation of inhibitory activities on HDAC8 revealed that compounds 3(a-c) and 5(a-c) showed similar inhibitory effects (IC50 ) ranging from 0.445 to 0.751 μM. Compounds 6(a-c) showed better affinity, with 6a (IC50  = 28 nM) and 6b (IC50  = 0.18 μM) showing potent inhibitory effects slightly lower than TSA (IC50  = 26 nM). These findings suggest that the studied compounds hold promise as potential candidates for further biological investigations.

HDAC8 as an emerging target in drug discovery with special emphasis on medicinal chemistry

HDAC8 catalyzes the deacetylation of both histones and nonhistone proteins. The abnormal expression of HDAC8 is associated with various pathological conditions causing cancer and other diseases like myopathies, Cornelia de Lange syndrome, renal fibrosis, and viral and parasitic infections. The substrates of HDAC8 are involved in diverse molecular mechanisms of cancer such as cell proliferation, invasion, metastasis and drug resistance. Based on the crystal structures and the key residues at the active site, HDAC8 inhibitors have been designed along the canonical pharmacophore. This article details the importance, recent advancements, and the structural and functional aspects of HDAC8 with special emphasis on the medicinal chemistry aspect of HDAC8 inhibitors that will help in developing novel epigenetic therapeutics.

Discovery of novel VEGFR2-TK inhibitors by phthalimide pharmacophore based virtual screening, molecular docking, MD simulation and DFT

Currently, numerous potent chemotherapeutic agents are available in the market but most of them show poor pharmacokinetics, lethal effects and drug resistance during their enduring use. The increased cancer cases, deaths and need of better treatment stimulates us to give newer lifesaving anticancer drugs. The phthalimide derivatives are structurally diverse and exert potential anticancer activity. In this regard, the 3D QSAR Pharmacophore model was developed and validated using fifty-eight phthalimide derivatives. The validation parameters corroborated the reliability and statistical robustness of CEASER Hypo 1. Three databases-NCI Open, Drug Bank, and Asinex were submitted to ADMET and drug-like filtering; 117893 drug-like compounds were mapped on CEASER Hypo 1; and 362 hits with IC50 <1 µM were discovered. These hits were docked on VEGFR2-TK, and in the form of results fifteen hits exhibited greater affinity than sorafenib. The top lead ASN 03206926 was subjected for MD simulation (100 ns) and RMSD, Rg, RMSF, number of hydrogen bonds, and SASA verified that the complex was stable, rigid and highly compact. Results demonstrated GLU885, PHE918, CYS919, LYS920, HIS1026, CYS1045, ASP1046 are the essential residues for favourable interactions. The binding free energy calculations support the affinity and stability revealed by docking and MD simulation. The DFT calculations, negative binding energy and lower HOMO-LUMO band gap revealed that the process is spontaneous and ASN 03206926 is very reactive. Following extensive analysis we suggest that the ASN 03206926 might be employed as a new VEGFR2-TK inhibitor for the treatment of breast and VEGFR2-TK associated cancers.Communicated by Ramaswamy H. Sarma.

A Therapeutic Perspective of HDAC8 in Different Diseases: An Overview of Selective Inhibitors

Histone deacetylases (HDACs) are epigenetic enzymes which participate in transcriptional repression and chromatin condensation mechanisms by removing the acetyl moiety from acetylated ε-amino group of histone lysines and other non-histone proteins. In recent years, HDAC8, a class I HDAC, has emerged as a promising target for different disorders, including X-linked intellectual disability, fibrotic diseases, cancer, and various neuropathological conditions. Selective HDAC8 targeting is required to limit side effects deriving from the treatment with pan-HDAC inhibitors (HDACis); thus, many endeavours have focused on the development of selective HDAC8is. In addition, polypharmacological approaches have been explored to achieve a synergistic action on multi-factorial diseases or to enhance the drug efficacy. In this frame, proteolysis-targeting chimeras (PROTACs) might be regarded as a dual-targeting approach for attaining HDAC8 proteasomal degradation. This review highlights the most relevant and recent advances relative to HDAC8 validation in various diseases, providing a snapshot of the current selective HDAC8is, with a focus on polyfunctional modulators.

Potential inhibitors of methionine aminopeptidase type II identified via structure-based pharmacophore modeling

Methionine aminopeptidase (MetAP2) is a metal-containing enzyme that removes initiator methionine from the N-terminus of a newly synthesized protein. Inhibition of the enzyme is crucial in diminishing cancer growth and metastasis. Fumagillin-a natural irreversible inhibitor of MetAP2-and its derivatives are used as potent MetAP2 inhibitors. However, because of their adverse effects, none of them has progressed to clinical studies. In search for potential reversible inhibitors, we built structure-based pharmacophore models using the crystal structure of MetAP2 complexed with fumagillin (PDB ID: 1BOA). The pharmacophore models were validated using Gunner-Henry scoring method. The best pharmacophore consisting of 1 H-bond donor, 1 H-bond acceptor, and 3 hydrophobic features was used to conduct pharmacophore-based virtual screening of ZINC15 database against MetAP2. The top 10 compounds with pharmacophore fit values > 3.00 were selected for further analysis. These compounds were subjected to absorption, distribution, metabolism, elimination, and toxicity (ADMET) prediction and found to have druglike properties. Furthermore, molecular docking calculations was performed on these hits using AutoDock4 to predict their binding mode and binding energy. Three diverse compounds: ZINC000014903160, ZINC000040174591, and ZINC000409110720 with respective binding energy/docking scores of - 9.22, - 9.21, and -817 kcal/mol, were submitted to 100 ns (MD) simulations using Nanoscale MD (NAMD) software. The compounds showed stable binding mode over time. Therefore, they may serve as a scaffold for further computational and experimental optimization toward the design of more potent and safer MetAP2 inhibitors.

Matairesinol, an active constituent of HC9 polyherbal formulation, exhibits HDAC8 inhibitory and anticancer activity

Histone deacetylase 8 (HDAC8) has emerged as a promising drug target for cancer therapeutics development. HDAC8 has been reported to regulate cancer cell proliferation, invasion and promote metastasis through modulation of cell cycle associated proteins. Of late, phytocompounds have been demonstrated to exhibit anticancer and anti-HDAC8 activity. Here, we have shown the HDAC8 inhibitory potential of an active phytocompound from HC9 (herbal composition-9), a polyherbal anticancer formulation based on the traditional Ayurvedic drug, Stanya Shodhan Kashaya. HC9 was recently reported to exhibit anticancer activity against breast cancer cells through induction of cell cycle arrest, decrease in migration and invasion as well as regulation of inflammation and chromatin modulators. In silico studies such as molecular docking, molecular dynamics (MD) simulation and binding free energy analyses showed greater binding energy values and interaction stability of MA with HDAC8 compared to other phytocompounds of HC9. Interestingly, in vitro validation confirmed the anti-HDAC8 activity of MA. Further, in vitro studies showed that MA significantly decreased the viability of breast and prostate cancer cell lines, thereby confirming its anticancer potential.

Discovery of novel N-substituted thiazolidinediones (TZDs) as HDAC8 inhibitors: in-silico studies, synthesis, and biological evaluation

Epigenetics plays a fundamental role in cancer progression, and developing agents that regulate epigenetics is crucial for cancer management. Among Class I and Class II HDACs, HDAC8 is one of the essential epigenetic players in cancer progression. Therefore, we designed, synthesized, purified, and structurally characterized novel compounds containing N-substituted TZD (P1-P25). Cell viability assay of all compounds on leukemic cell lines (CEM, K562, and KCL22) showed the cytotoxic potential of P8, P9, P10, P12, P19, and P25. In-vitro screening of different HDACs isoforms revealed that P19 was the most potent and selective inhibitor for HDAC8 (IC50 - 9.3 μM). Thermal shift analysis (TSA) confirmed the binding of P19 to HDAC8. In-vitro screening of all compounds on the transport activity of GLUT1, GLUT4, and GLUT5 indicated that P19 inhibited GLUT1 (IC50 - 28.2 μM). P10 and P19 induced apoptotic cell death in CEM cells (55.19% and 60.97% respectively) and P19 was less cytotoxic on normal WBCs (CC50 - 104.2 μM) and human fibroblasts (HS27) (CC50 - 105.0 μM). Thus, among this novel series of TZD derivatives, compound P19 was most promising HDAC8 inhibitor and cytotoxic on leukemic cells. Thus, P19 could serve as a lead for further development of optimized molecules with enhanced selectivity and potency.

Identification of HDAC6 selective inhibitors: pharmacophore based virtual screening, molecular docking and molecular dynamics simulation

HDAC6 regulates the expression and activity of various tumor-related proteins, but currently there is no selective inhibitor targeting HDAC6 for clinical application. In order to discover novel HDAC6 inhibitors, virtual screening methods comprised of pharmacophore based virtual screening, molecular docking and molecular dynamics (MD) simulations were employed. 15 molecules were obtained after virtual screening. After in vitro bioassays, two of the hits showed inhibition activity against HDAC6, among which the inhibition activity of G1 to HDAC6 reached 81% at concentration of 20 μM. In addition, the inhibitory activity against HDAC1 and HDAC10 demonstrated that G1 and G10 were highly selective to HDAC6. The analysis of the binding modes of G1 and G10 provides a reference for further development of highly active HDAC6 inhibitors. Communicated by Ramaswamy H. Sarma.

Crystallographic structure versus homology model: a case study of molecular dynamics simulation of human and zebrafish histone deacetylase 10

Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100 ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design.Communicated by Ramaswamy H. Sarma.