Virtual screening and experimental validation of novel histone deacetylase inhibitors

STS1 and STS2 Phosphatase Inhibitor Baicalein Enhances the Expansion of Hematopoietic and Progenitor Stem Cells and Alleviates 5-Fluorouracil-Induced Myelosuppression

STS1 and STS2, as the protein phosphatases that dephosphorylate FLT3 and cKIT, negatively regulate the self-renewal and differentiation of hematopoietic stem and progenitor cells (HSPCs). To obtain the small molecule inhibitors of STS1/STS2 phosphatase activity used to expand HSPCs both in vitro and in vivo, we establish an in vitro phosphatase assay using the recombinant proteins of the STS1/STS2 histidine phosphatase (HP) domain, by which we screened out baicalein (BC) as one of the effective inhibitors targeting STS1 and STS2. Then, we further demonstrate the direct binding of BC with STS1/STS2 using molecular docking and capillary electrophoresis and verify that BC can restore the phosphorylation of FLT3 and cKIT from STS1/STS2 inhibition. In a short-term in vitro culture, BC promotes profound expansion and enhances the colony-forming capacity of both human and mouse HSPCs along with the elevation of phospho-FLT3 and phospho-cKIT levels. Likewise, in vivo administration with BC significantly increases the proportions of short-term hematopoietic stem cells (ST-HSCs), multipotent progenitors (MPPs) and especially long-term HSCs (LT-HSCs) in healthy mouse bone marrow and increases the numbers of colony-forming units (CFU) formed by HSPCs as well. More importantly, pre-administration of BC significantly enhances the survival of mice with lethal 5-fluorouracil (5-FU) injection due to the alleviation of 5-FU-induced myelosuppression, as evidenced by the recovery of bone marrow histologic injury, the increased proportions of LT-HSCs, ST-HSCs and MPPs, and enhanced colony-forming capacity. Collectively, our study not only suggests BC as one of the small molecule candidates to stimulate HSPC expansion both in vitro and in vivo when needed in either physiologic or pathologic conditions, but also supports STS1/STS2 as potential therapeutic drug targets for HSPC expansion and hematopoietic injury recovery.

Discovery of a novel HDACi structure that inhibits the proliferation of ovarian cancer cells in vivo and in vitro

Histone deacetylases (HDACs) exhibit increased expression in cancer and promote oncogenesis via the acetylation of or interactions with key transcriptional regulators. HDAC inhibitors (HDACis) decrease HDAC activity to selectively inhibit the occurrence and development of tumors. Our study screened and obtained a new HDACi structure. In vitro experiments have showed that among the leads, Z31216525 significantly inhibited the proliferation and induced the apoptosis of epithelial ovarian cancer (EOC) cells. In vivo experiments demonstrated that compared to the control, Z31216525 significantly inhibited tumor growth and showed very low toxicity. Further mechanistic studies revealed that Z31216525 may exert an antitumor effect by inhibiting the expression of the c-Myc gene. Collectively, our studies identified a novel HDACi that is expected to become a new potential therapeutic drug for EOC and has important value for the design of new HDACi structures.

Discovery of novel isoform-selective histone deacetylases 5 and 9 inhibitors through combined ligand-based pharmacophore modeling, molecular mocking, and molecular dynamics simulations for cancer treatment

Class IIa histone deacetylases (HDACs) 5 and 9 play crucial roles in several human disorders such as cancer, making them important targets for drug design. Continuous research is pursed to overcome the cytotoxicity side effect that comes with the currently available broad-spectrum HDACs inhibitors. Herein, common features of active HDACs inhibitors in clinical trials and use have been calculated to generate the best pharmacophore hypothesis. Guner-Henry scoring system was used to validate the generated hypotheses. Hypo1 of HDAC5 and Hypo2 of HDAC9 exhibited the most statistically significance hypotheses. Compounds with fit value of 3 and more were examined by QuickVina 2 docking tool to calculate their binding affinity toward all class IIa HDACs. A total of 6 potential selective compounds were subjected to 100 molecular dynamics (MD) simulation to examine their binding modes. The free binding energy calculations were computed according to the MM-PBSA method. Proposed selective compounds displayed good stability with their targets and thus they may offer potent leads for the designing of HDAC5 and HDAC9 isoform selective inhibitors.

A unique ligand-steered strategy for CC chemokine receptor 2 homology modeling to facilitate structure-based virtual screening

CC chemokine receptor 2 (CCR2) antagonists that disrupt CCR2/MCP-1 interaction are expected to treat a variety of inflammatory and autoimmune diseases. The lack of CCR2 crystal structure limits the application of structure-based drug design (SBDD) to this target. Although a few three-dimensional theoretical models have been reported, their accuracy remains to be improved in terms of templates and modeling approaches. In this study, we developed a unique ligand-steered strategy for CCR2 homology modeling. It starts with an initial model based on the X-ray structure of the closest homolog so far, that is, CXCR4. Then, it uses Elastic Network Normal Mode Analysis (EN-NMA) and flexible docking (FD) by AutoDock Vina software to generate ligand-induced fit models. It selects optimal model(s) as well as scoring function(s) via extensive evaluation of model performance based on a unique benchmarking set constructed by our in-house tool, that is, MUBD-DecoyMaker. The model of 81_04 presents the optimal enrichment when combined with the scoring function of PMF04, and the proposed binding mode between CCR2 and Teijin lead by this model complies with the reported mutagenesis data. To highlight the advantage of our strategy, we compared it with the only reported ligand-steered strategy for CCR2 homology modeling, that is, Discovery Studio/Ligand Minimization. Lastly, we performed prospective virtual screening based on 81_04 and CCR2 antagonist bioassay. The identification of two hit compounds, that is, E859-1281 and MolPort-007-767-945, validated the efficacy of our model and the ligand-steered strategy.

Cannabinoid Receptor Subtype 2 (CB2R) in a Multitarget Approach: Perspective of an Innovative Strategy in Cancer and Neurodegeneration

The cannabinoid receptor subtype 2 (CB2R) represents an interesting and new therapeutic target for its involvement in the first steps of neurodegeneration as well as in cancer onset and progression. Several studies, focused on different types of tumors, report a promising anticancer activity induced by CB2R agonists due to their ability to reduce inflammation and cell proliferation. Moreover, in neuroinflammation, the stimulation of CB2R, overexpressed in microglial cells, exerts beneficial effects in neurodegenerative disorders. With the aim to overcome current treatment limitations, new drugs can be developed by specifically modulating, together with CB2R, other targets involved in such multifactorial disorders. Building on successful case studies of already developed multitarget strategies involving CB2R, in this Perspective we aim at prompting the scientific community to consider new promising target associations involving HDACs (histone deacetylases) and σ receptors by employing modern approaches based on molecular hybridization, computational polypharmacology, and machine learning algorithms.

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.

Hypericin attenuates nonalcoholic fatty liver disease and abnormal lipid metabolism via the PKA-mediated AMPK signaling pathway in vitro and in vivo

Nonalcoholic fatty liver disease (NAFLD) is the most common liver disease worldwide and constitutes a major risk factor for progression to cirrhosis, liver failure and hepatocellular carcinoma (HCC). The occurrence of NAFLD is closely associated with abnormal lipid metabolism and implies a high risk of type 2 diabetes and cardiovascular disease. Therefore, specific and effective drugs for the prevention and treatment of NAFLD are necessary. Hypericin (HP) is one of the main active ingredients of Hypericum perforatum L., and we previously revealed its protective role in islet β-cells and its effects against type 2 diabetes. In this study, we aimed to explore the preventive and therapeutic effects of HP against NAFLD and the underlying mechanisms in vitro and in vivo. Here, we demonstrated that HP improved cell viability by reducing apoptosis and attenuated lipid accumulation in hepatocytes both in vitro and in vivovia attenuating oxidative stress, inhibiting lipogenesis and enhancing lipid oxidization. Thus, HP exhibited significant preventive and therapeutic effects against HFHS-induced NAFLD and dyslipidemia in mice. Furthermore, we demonstrated that HP directly bound to PKACα and activated PKA/AMPK signaling to elicit its effects against NAFLD, suggesting that PKACα is one of the drug targets of HP. In addition, the enhancing effect of HP on lipolysis in adipocytes through the activation of PKACα was also elucidated. Together, the conclusions indicated that HP, of which one of the targets is PKACα, has the potential to be used as a preventive or therapeutic drug against NAFLD or abnormal lipid metabolism in the future.

An Improved Protocol for the Virtual Screening Discovery of Novel Histone Deacetylase Inhibitors

Histone deacetylases (HDACs) are enzymes that play a key role in structural modification and gene expression. The overexpression of HDAC is associated with cancer, and thus inhibiting the enzyme could be an efficient cancer therapy. To discover new HDAC inhibitors (HDACis), we proposed an improved protocol combining a hierarchical pharmacophore search, molecular docking, and molecular dynamic simulations. The test results showed that the improved screening protocol effectively reduced the false-positive rates of drug-like chemicals. Based on the protocol, we obtained 16 hit compounds as potential HDACis from the Life Chemicals database. Enzyme inhibition experiments showed that two of the hit chemical compounds had HDAC-inhibitory effects. In vitro assays showed that Z165155756 could selectively inhibit the proliferation of cancer cells and specifically promoted apoptosis and induced G1/S phase arrest in A2780 cells. It may have potential therapeutic effects in ovarian cancer and is worthy of further investigation.

Turning to Computer-aided Drug Design in the Treatment of Diffuse Large B-cell Lymphoma: Has it been Helpful?

INTRODUCTION

Amidst the numerous effective therapeutic options available for the treatment of Diffuse Large B-cell Lymphoma (DLBCL), about 30-40% of patients treated with first-line chemoimmunotherapy still experience a relapse or refractory DLBCL. This has necessitated a continuous search for new therapeutic agents to augment the existing therapeutic arsenal.

METHODS

The dawn of Computer-Aided Drug Design (CADD) in the drug discovery process has accounted for persistency in the application of computational approaches either alone or in combinatorial strategies with experimental methods towards the identification of potential hit compounds with high therapeutic efficacy in abrogating DLBCL.

RESULTS

This review showcases the interventions of structure-based and ligand-based computational approaches which have led to the identification of numerous small molecule inhibitors against implicated targets in DLBCL therapy, even though many of these potential inhibitors are piled-up awaiting further experimental validation and exploration.

CONCLUSION

We conclude that a successful and a conscious amalgamation of CADD and experimental approaches could pave the way for the discovery of the next generation potential leads in DLBCL therapy with improved activities and minimal toxicities.

Discovery of novel acylhydrazone neuraminidase inhibitors

Neuraminidase (NA) plays a crucial role in the replication and transmission of influenza virus. NA inhibitors have been developed as effective treatments for influenza A and B infections. In this paper, a new lead neuraminidase inhibitor 6a (IC₅₀ = 7.10 ± 0.2 μM) was discovered by ligand-based virtual screening, receptor-based virtual screening, molecular dynamics simulation (MD), and bioassay validation. MD simulation indicates that the morpholinyl group of 6a could be embedded in 430-loop of NA. To exploit the 430-loop in the active site, a series of novel acylhydrazone NA inhibitors 6b-6g were designed and synthesized based on the lead compound 6a. Compound 6e exerts the most potency, with IC₅₀ value of 2.37 ± 0.5 μM against NA, which is lower than that of oseltamivir carboxylate (OC) (IC₅₀ = 3.84 μM). Overall, this work provided unique insights in the discovery of potent inhibitors against NA.

A novel benzamine lead compound of histone deacetylase inhibitor ZINC24469384 can suppresses HepG2 cells proliferation by upregulating NR1H4

Histone deacetylases (HDACs) can enzymatically transferred acetyl functional group from protein or lysine residues of histone, so they can regulate the expression of lots of genes. Now HDACs are used as drug targets and many HDAC inhibitors (HDACis) were approved for cancer therapy or in clinical trials. However, the physiological mechanisms and regulatory processes of HDACi anti-cancer effects are largely unexplored and uncompleted. Here we use the virtual screening workflow obtained 25 hit compounds and ZINC24469384 can significantly inhibit HDAC activity while arrest cell cycle at G1/S phase and significantly induced HepG2 cell apoptosis, time-course RNA-seq demonstrate that HepG2 cells transcriptionally respond to ZINC24469384. Pathway analysis of DEGs and DASGs reveal that NR1H4 may play an important role in ZINC24469384-induced anti-proliferation effect and is dramatically alleviated by down-regulating the SOCS2 expression and promoting STAT3 phosphorylation in knockdown NR1H4 HepG2 cells. Analysis based on TCGA database indicated that NR1H4 and SOCS2 were downregulated in liver cancer, this suggest NR1H4 and SOCS2 may play an important role in tumorigenesis. These results indicated that ZINC24469384 is a novel benzamine lead compound of HDACi and provides a novel mechanism for HDACi to inhibit cancer.

Pharmacophore-based virtual screening for identification of potential selective inhibitors of human histone deacetylase 6

Histone deacetylase (HDAC) 6 plays a role in oncogenic transformation and cancer metastasis via tubulin deacetylation, making it a critical target for anticancer drug design. However, lack of selectivity shown by many of the current HDAC6 inhibitors in clinical use and trials prompts the continuous search for selective inhibitors. Here, 10 pharmacophore hypotheses were developed based on the 3D common features of training set of 20 HDAC inhibitors in clinical use and trials. The hypotheses were validated using a test set of another 20 HDAC inhibitors along with 400 inactive (decoys) molecules based on Güner-Henry pharmacophore scoring method. Hypothesis 1 consisting of 1 H-bond donor, 1 H-bond acceptor and 2 hydrophobic features, was used to screen "DruglikeDiverse" database using Biovia Discovery Studio 4.5. The top 10 hit compounds were selected based on the pharmacophore fit values (>3.00). Their binding affinity against HDAC6 compared to class I HDACs (1, 2, 3 & 8) and a class IIa member (HDAC7), was calculated by molecular docking using AutoDock4. The stability of binding modes of 2 potential HDAC6-selective inhibitors (ENA501965 and IBS399024) was examined by 30 ns-molecular dynamics (MD) simulation using nanoscale MD (NAMD) software. Both ligands showed potential stability in HDAC6 active site over time. Therefore, these may provide additional scaffolds for further optimization towards the design of safe, potent and selective HDAC6 inhibitors.

Prediction Methods of Herbal Compounds in Chinese Medicinal Herbs

Chinese herbal medicine has recently gained worldwide attention. The curative mechanism of Chinese herbal medicine is compared with that of western medicine at the molecular level. The treatment mechanism of most Chinese herbal medicines is still not clear. How do we integrate Chinese herbal medicine compounds with modern medicine? Chinese herbal medicine drug-like prediction method is particularly important. A growing number of Chinese herbal source compounds are now widely used as drug-like compound candidates. An important way for pharmaceutical companies to develop drugs is to discover potentially active compounds from related herbs in Chinese herbs. The methods for predicting the drug-like properties of Chinese herbal compounds include the virtual screening method, pharmacophore model method and machine learning method. In this paper, we focus on the prediction methods for the medicinal properties of Chinese herbal medicines. We analyze the advantages and disadvantages of the above three methods, and then introduce the specific steps of the virtual screening method. Finally, we present the prospect of the joint application of various methods.

Genome-wide analysis of circular RNAs in bovine cumulus cells treated with BMP15 and GDF9

Circular RNAs (circRNAs) are important members of the non-coding RNA family, and those relating to animal physiologies have been widely studied in recent years. This study aimed to explore the roles of circRNAs in the regulation of follicular development. We constructed four bovine cumulus cell cDNA libraries, including a negative control group (NC) and groups treated with BMP15, GDF9 and BMP15 + GDF9, and we sequenced the libraries on the Illumina HiSeq Xten platform. We identified 1706 circRNAs and screened for differential circRNA expression. We conducted a bioinformatics analysis of these circRNAs and screened for differential circRNAs. Functional annotation and enrichment analysis of the host genes showed that the differential circRNAs were related to locomotion, reproduction, biological adhesion, growth, rhythmic processes, biological phases and hormone secretion. According to the differential expression of circRNA between groups, there were 3 up-regulated and 6 down-regulated circRNAs in the BMP15 group as well as 12 up-regulated and 24 down-regulated circRNAs in the GDF9 group. Co-addition of both BMP15 and GDF9 resulted in 15 up-regulated and 13 down-regulated circRNAs. circ_n/a_75,circ_12691_1 and circ_n/a_303 were altered in both the BMP15 and GDF9 groups as well as in the BMP15 + GDF9 combination group. We focused on these three circRNAs because they were potentially associated with the additive effect of BMP15 and GDF9. Quantitative PCR analysis showed that the expression levels of these three circRNAs were consistent with the sequencing results. In addition, the target miRNAs of circ_n/a_75 and circ_n/a_303, miR-339a, miR-2400 and miR-30c, were down-regulated in the experimental group, which was in contrast to the circRNAs trend. These findings demonstrated that BMP15 and GDF9 may regulate the target gene through circRNA, as a miRNA sponge, in order to regulate the status of bovine cumulus cells and affect follicular development.

Carboxylic acid derivatives display potential selectivity for human histone deacetylase 6: Structure-based virtual screening, molecular docking and dynamics simulation studies

Human histone deacetylase 6 (HDAC6) has been shown to play a major role in oncogenic cell transformation via deacetylation of α-tubulin, making it a viable target of anticancer drug design and development. The crystal structure of HDAC6 catalytic domain 2 has been recently made available, providing avenues for structure-based drug design campaign. Here, in our continuous effort to identify potentially selective HDAC6 inhibitors, structure-based virtual screening of ∼72 461 compounds was carried out using Autodock Vina. The top 100 compounds with calculated ΔG < -10 kcal/mol were manually inspected for binding mode orientation. Furthermore, the top 20 compounds with reasonable binding modes were evaluated for selectivity by further docking against HDAC6 and HDAC7 using Autodock4. Four compounds with a carboxylic fragment, displayed potential selectivity for HDAC6 over HDAC7, and were found to have good druglike and ADMET properties. Their docking complexes were then submitted to 10 ns-molecular dynamics (MD) simulation using nanoscale MD (NAMD) software, to examine the stability of ligand binding modes. These predicted inhibitors remained bound to HDAC6 in the presence of water and ions, and the root-mean-square deviation (RMSD), radius of gyration (Rg) and nonbond distance (protein-ligand) profiles suggested that they might be stable over time of the simulation. This study may provide scaffolds for further lead optimization towards the design of HDAC6 inhibitors with improved selectivity.

Exploration of the binding pocket of histone deacetylases: the design of potent and isoform-selective inhibitors

Histone deacetylases (HDACs) are enzymes that act on histone proteins to remove the acetyl group and thereby regulate the chromatin state. HDACs act not only on histone protein but also nonhistone proteins that are key players in cellular processes such as the cell cycle, signal transduction, apoptosis, and more. "Classical" HDACs have been shown to be promising targets for anticancer drug design and development. However, the selectivity of HDAC inhibitors for HDAC isoforms remains the motivation of current research in this field. Here, we explored Class I HDACs and HDAC6 by sequence alignment and structural superimposition, catalytic channel extraction, and identification of critical residues involved in HDAC catalysis. Based on the general pharmacophore features of known HDAC inhibitors, we developed a library of compounds by scaffold hopping on a fragment hit identified via structurebased virtual screening of the molecular fragment library retrieved from the Otava database. Molecular docking assay revealed five of these compounds to have increased potency and selectivity for HDACs 1 and 2. Furthermore, their predicted absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties were consistent with those of drug-like compounds. With further modelingbased and experimental investigations, we believe that these findings may offer additional potential HDAC inhibitors with improved selectivity.

Attenuated SUMOylation of sirtuin 1 in premature neonates with bronchopulmonary dysplasia

A prospective study was performed to investigate the effects of hyperoxia on the expression of small ubiquitin‑related modifier (SUMO) and sirtuin 1 (SIRT1) proteins, and to examine interactions between these proteins in premature neonates with bronchopulmonary dysplasia (BPD). Peripheral blood mononuclear cells (PBMCs) were isolated from residual venous blood samples of 20 premature infants with BPD and 20 gender‑matched premature infants without BPD (non‑BPD group). Expression levels of SUMO and SIRT1 proteins in PBMCs were assessed by western blot analysis, and their interactions in PBMCs were detected using the immunoprecipitation assay. Based on the fraction of inspired oxygen (FiO2) administered, neonates were divided into normoxia, low‑(21%<FiO2<30%), medium‑(30%≤FiO2<40%) and high‑oxygen (FiO2≥40%) groups. Expression levels of SUMO1 and SUMO2/3 proteins in the normoxia group were significantly lower than those in the medium‑ or high‑oxygen groups (P<0.01), but were comparable to those in the low‑oxygen group. SIRT1 expression levels in both the medium‑ and high‑oxygen groups were significantly lower than those in the normoxia group (P<0.01). In the BPD group, the expression of SIRT1 protein was significantly lower (P<0.01), and its interaction with SUMO1 and SUMO2/3 was significantly attenuated compared with that in the non‑BPD group (P<0.01). Supplemental oxygen with FiO2≥30% was associated with upregulation of SUMO1 and SUMO2/3 expression and downregulation of SIRT1 expression. The present findings suggest that decreased SIRT1 expression and its SUMOylation by SUMO1 and SUMO2/3 may be associated with the development of BPD.

The Search for Potent, Small-Molecule HDACIs in Cancer Treatment: A Decade After Vorinostat

Histone deacetylases (HDACs) play a crucial role in the remodeling of chromatin, and are involved in the epigenetic regulation of gene expression. In the last decade, inhibition of HDACs came out as a target for specific epigenetic changes associated with cancer and other diseases. Until now, more than 20 HDAC inhibitors (HDACIs) have entered clinical studies, and some of them (e.g., vorinostat, romidepsin) have been approved for the treatment of cutaneous T-cell lymphoma. This review provides an overview of current knowledge, progress, and molecular mechanisms of HDACIs, covering a period from 2011 until 2015.

A Thoroughly Validated Virtual Screening Strategy for Discovery of Novel HDAC3 Inhibitors

Histone deacetylase 3 (HDAC3) has been recently identified as a potential target for the treatment of cancer and other diseases, such as chronic inflammation, neurodegenerative diseases, and diabetes. Virtual screening (VS) is currently a routine technique for hit identification, but its success depends on rational development of VS strategies. To facilitate this process, we applied our previously released benchmarking dataset, i.e., MUBD-HDAC3 to the evaluation of structure-based VS (SBVS) and ligand-based VS (LBVS) combinatorial approaches. We have identified FRED (Chemgauss4) docking against a structural model of HDAC3, i.e., SAHA-3 generated by a computationally inexpensive “flexible docking”, as the best SBVS approach and a common feature pharmacophore model, i.e., Hypo1 generated by Catalyst/HipHop as the optimal model for LBVS. We then developed a pipeline that was composed of Hypo1, FRED (Chemgauss4), and SAHA-3 sequentially, and demonstrated that it was superior to other combinations in terms of ligand enrichment. In summary, we present the first highly-validated, rationally-designed VS strategy specific to HDAC3 inhibitor discovery. The constructed pipeline is publicly accessible for the scientific community to identify novel HDAC3 inhibitors in a time-efficient and cost-effective way.