The discovery of novel HDAC3 inhibitors via virtual screening and in vitro bioassay

Stack-HDAC3i: A high-precision identification of HDAC3 inhibitors by exploiting a stacked ensemble-learning framework

Epigenetics involves reversible modifications in gene expression without altering the genetic code itself. Among these modifications, histone deacetylases (HDACs) play a key role by removing acetyl groups from lysine residues on histones. Overexpression of HDACs is linked to the proliferation and survival of tumor cells. To combat this, HDAC inhibitors (HDACi) are commonly used in cancer treatments. However, pan-HDAC inhibition can lead to numerous side effects. Therefore, isoform-selective HDAC inhibitors, such as HDAC3i, could be advantageous for treating various medical conditions while minimizing off-target effects. To date, computational approaches that use only the SMILES notation without any experimental evidence have become increasingly popular and necessary for the initial discovery of novel potential therapeutic drugs. In this study, we develop an innovative and high-precision stacked-ensemble framework, called Stack-HDAC3i, which can directly identify HDAC3i using only the SMILES notation. Using an up-to-date benchmark dataset, we first employed both molecular descriptors and Mol2Vec embeddings to generate feature representations that cover multi-view information embedded in HDAC3i, such as structural and contextual information. Subsequently, these feature representations were used to train baseline models using nine popular ML algorithms. Finally, the probabilistic features derived from the selected baseline models were fused to construct the final stacked model. Both cross-validation and independent tests showed that Stack-HDAC3i is a high-accuracy prediction model with great generalization ability for identifying HDAC3i. Furthermore, in the independent test, Stack-HDAC3i achieved an accuracy of 0.926 and Matthew's correlation coefficient of 0.850, which are 0.44-6.11% and 0.83-11.90% higher than its constituent baseline models, respectively.

Exploring different classification-dependent QSAR modelling strategies for HDAC3 inhibitors in search of meaningful structural contributors

Histone deacetylase 3 (HDAC3), a Zn2+-dependent class I HDACs, contributes to numerous disorders such as neurodegenerative disorders, diabetes, cardiovascular disease, kidney disease and several types of cancers. Therefore, the development of novel and selective HDAC3 inhibitors might be promising to combat such diseases. Here, different classification-based molecular modelling studies such as Bayesian classification, recursive partitioning (RP), SARpy and linear discriminant analysis (LDA) were conducted on a set of HDAC3 inhibitors to pinpoint essential structural requirements contributing to HDAC3 inhibition followed by molecular docking study and molecular dynamics (MD) simulation analyses. The current study revealed the importance of hydroxamate function for Zn2+ chelation as well as hydrogen bonding interaction with Tyr298 residue. The importance of hydroxamate function for higher HDAC3 inhibition was noticed in the case of Bayesian classification, recursive partitioning and SARpy models. Also, the importance of substituted thiazole ring was revealed, whereas the presence of linear alkyl groups with carboxylic acid function, any type of ester function, benzodiazepine moiety and methoxy group in the molecular structure can be detrimental to HDAC3 inhibition. Therefore, this study can aid in the design and discovery of effective novel HDAC3 inhibitors in the future.

HDAC3 inhibitors: a patent review of their broad-spectrum applications as therapeutic agents

INTRODUCTION

Histone deacetylases (HDACs) are a class of zinc-dependent enzymes. They maintain acetylation homeostasis, with numerous biological functions and are associated with many diseases. HDAC3 strictly requires multi-subunit complex formation for activity. It is associated with the progression of numerous non-communicable diseases. Its widespread involvement in diseases makes it an epigenetic drug target. Preexisting HDAC3 inhibitors have many uses, highlighting the need for continued research in the discovery of HDAC3-selective inhibitors.

AREA COVERED

This review provides an overview of 24 patents published from 2010 to 2023, focusing on compounds that inhibit the HDAC3 isoenzyme.

EXPERT OPINION

HDAC3-selective inhibitors - pivotal for pharmacological applications, as single or combination therapies - are gaining traction as a strategy to move away from complications laden pan-HDAC inhibitors. Moreover, there is an unmet need for HDAC3 inhibitors with alternative zinc-binding groups (ZBGs) because some preexisting ZBGs have limitations related to toxicity and side effects. Difficulties in achieving HDAC3 selectivity may be due to isoform selectivity. However, advancements in computer-aided drug design and experimental data of HDAC3 3D co-crystallized models could lead to the discovery of novel HDAC3-selective inhibitors, which bear alternative ZBGs with balanced selectivity for HDAC3 and potency.

Discovery of Novel Metalloenzyme Inhibitors Based on Property Characterization: Strategy and Application for HDAC1 Inhibitors

Metalloenzymes are ubiquitously present in the human body and are relevant to a variety of diseases. However, the development of metalloenzyme inhibitors is limited by low specificity and poor drug-likeness associated with metal-binding fragments (MBFs). A generalized drug discovery strategy was established, which is characterized by the property characterization of zinc-dependent metalloenzyme inhibitors (ZnMIs). Fifteen potential Zn2+-binding fragments (ZnBFs) were identified, and a customized pharmacophore feature was defined based on these ZnBFs. The customized feature was set as a required feature and applied to a search for novel inhibitors for histone deacetylase 1 (HDAC1). Ten potential HDAC1 inhibitors were recognized, and one of them (compound 9) was a known potent HDAC1 inhibitor. The results demonstrated the effectiveness of our strategy to identify novel inhibitors for zinc-dependent metalloenzymes.

Hydroxamic acid derivatives as selective HDAC3 inhibitors: computer-aided drug design strategies

Histone deacetylases (HDACs) are critical epigenetic drug targets that have gained significant attention in the scientific community for the treatment of cancer. The currently marketed HDAC inhibitors lack selectivity for the various HDAC isoenzymes. Here, we describe our protocol for the discovery of novel potential hydroxamic acid based HDAC3 inhibitors through pharmacophore modeling, virtual screening, docking, molecular dynamics (MD) simulation and toxicity studies. The ten pharmacophore hypotheses were established, and their reliability was validated by different ROC (receiving operator curve) analysis. Among them, the best model (Hypothesis 9 or RRRA) was employed for searching SCHEMBL, ZINC and MolPort database to screen out hit molecules as selective HDAC3 inhibitors, followed by different docking stages. MD simulation (50 ns) and MMGBSA study were performed to study the stability of ligand binding modes and with the help of trajectory analysis, to calculate the ligand-receptor complex RMSD (root-mean-square deviation), RMSF (root-mean-square fluctuation) and H-bond distance, etc. Finally, in-silico toxicity studies were performed on top screened molecules and compared with reference drug SAHA and established structure-activity relationship (SAR). The results indicated that compound 31, with high inhibitory potency and less toxicity (probability value 0.418), is suitable for further experimental analysis.Communicated by Ramaswamy H. Sarma.

Pharmacophore-based virtual screening, 3D QSAR, Docking, ADMET, and MD simulation studies: An in silico perspective for the identification of new potential HDAC3 inhibitors

Histone deacetylase 3 (HDAC3) is an epigenetic regulator that involves gene expression, apoptosis, and cell cycle progression, and the overexpression of HDAC3 is accountable for several cancers, neurodegeneracy, and many other diseases. Therefore, HDAC3 emerged as a promising drug target for the novel drug design. Here, we carried out the pharmacophore modeling using 50 benzamide-based HDAC3 selective inhibitors and utilized it for PHASE ligand screening to retrieve the hits with similar pharmacophore features. The dataset inhibitors of best hypotheses used to build the 3D QSAR model and the generated 3D QSAR model resulted in good PLS statistics with a regression coefficient (R2) of 0.89, predictive coefficient (Q2) of 0.88, and Pearson-R factor of 0.94 indicating its excellent predictive ability. The hits retrieved from pharmacophore-based virtual screening were subjected to docking against HDAC3 for the identification of potential inhibitors. A total of 10 hitsM1 to M10 were ranked using their scoring functions and further subject to lead optimization. The Prime MM/GBSA, AutoDock binding free energies, and ADMET studies were implemented for the selection of lead candidates. The four ligand molecules M1, M2, M3, and M4 were identified as potential leads against HDAC3 after lead optimization. The top two leads M1 and M2 were subjected to MD simulations for their stability evaluation with HDAC3. The newly designed leads M11 and M12 were identified as HDAC3 potential inhibitors from MD simulations studies. Therefore, the outcomes of the present study could provide insights into the discovery of new potential HDAC3 inhibitors with improved selectivity and activity against a variety of cancers and neurodegenerative diseases.

Medicinal chemistry advances in targeting class I histone deacetylases

Histone deacetylases (HDACs) are a class of zinc (Zn)-dependent metalloenzymes that are responsible for epigenetic modifications. HDACs are largely associated with histone proteins that regulate gene expression at the DNA level. This tight regulation is controlled by acetylation [via histone acetyl transferases (HATs)] and deacetylation (via HDACs) of histone and non-histone proteins that alter the coiling state of DNA, thus impacting gene expression as a downstream effect. For the last two decades, HDACs have been studied extensively and indicated in a range of diseases where HDAC dysregulation has been strongly correlated with disease emergence and progression-most prominently, cancer, neurodegenerative diseases, HIV, and inflammatory diseases. The involvement of HDACs as regulators in these biochemical pathways established them as an attractive therapeutic target. This review summarizes the drug development efforts exerted to create HDAC inhibitors (HDACis), specifically class I HDACs, with a focus on the medicinal chemistry, structural design, and pharmacology aspects of these inhibitors.

Nutritional Provision of Iron Complexes by the Major Allergen Alt a 1 to Human Immune Cells Decreases Its Presentation

Alternaria alternata is a common fungus strongly related with severe allergic asthma, with 80% of affected individuals being sensitized solely to its major allergen Alt a 1. Here, we assessed the function of Alt a 1 as an innate defense protein binding to micronutrients, such as iron-quercetin complexes (FeQ2), and its impact on antigen presentation in vitro. Binding of Alt a 1 to FeQ2 was determined in docking calculations. Recombinant Alt a 1 was generated, and binding ability, as well as secondary and quaternary structure, assessed by UV-VIS, CD, and DLS spectroscopy. Proteolytic functions were determined by casein and gelatine zymography. Uptake of empty apo- or ligand-filled holoAlt a 1 were assessed in human monocytic THP1 cells under the presence of dynamin and clathrin-inhibitors, activation of the Arylhydrocarbon receptor (AhR) using the human reporter cellline AZ-AHR. Human PBMCs were stimulated and assessed for phenotypic changes in monocytes by flow cytometry. Alt a 1 bound strongly to FeQ2 as a tetramer with calculated Kd values reaching pico-molar levels and surpassing affinities to quercetin alone by a factor of 5000 for the tetramer. apoAlt a 1 but not holoAlta 1 showed low enzymatic activity against casein as a hexamer and gelatin as a trimer. Uptake of apo- and holo-Alt a 1 occurred partly clathrin-dependent, with apoAlt a 1 decreasing labile iron in THP1 cells and holoAlt a 1 facilitating quercetin-dependent AhR activation. In human PBMCs uptake of holoAlt a 1 but not apoAlt a 1 significantly decreased the surface expression of the costimulatory CD86, but also of HLADR, thereby reducing effective antigen presentation. We show here for the first time that the presence of nutritional iron complexes, such as FeQ2, significantly alters the function of Alt a 1 and dampens the human immune response, thereby supporting the notion that Alt a 1 only becomes immunogenic under nutritional deprivation.

Identification of novel leads as potent inhibitors of HDAC3 using ligand-based pharmacophore modeling and MD simulation

In the landscape of epigenetic regulation, histone deacetylase 3 (HDAC3) has emerged as a prominent therapeutic target for the design and development of candidate drugs against various types of cancers and other human disorders. Herein, we have performed ligand-based pharmacophore modeling, virtual screening, molecular docking, and MD simulations to design potent and selective inhibitors against HDAC3. The predicted best pharmacophore model 'Hypo 1' showed excellent correlation (R2 = 0.994), lowest RMSD (0.373), lowest total cost value (102.519), and highest cost difference (124.08). Hypo 1 consists of four salient pharmacophore features viz. one hydrogen bond acceptor (HBA), one ring aromatic (RA), and two hydrophobic (HYP). Hypo 1 was validated by Fischer's randomization with a 95% of confidence level and the external test set of 60 compounds with a good correlation coefficient (R2 = 0.970). The virtual screening of chemical databases, drug-like properties calculations followed by molecular docking resulted in identifying 22 representative hit compounds. Performed 50 ns of MD simulations on top three hits were retained the salient π-stacking, Zn2+ coordination, hydrogen bonding, and hydrophobic interactions with catalytic residues from the active site pocket of HDAC3. Total binding energy calculated by MM-PBSA showed that the Hit 1 and Hit 2 formed stable complexes with HDAC3 as compared to reference TSA. Further, the PLIP analysis showed a close resemblance between the salient pharmacophore features of Hypo 1 and the presence of molecular interactions in co-crystallized FDA-approved drugs. We conclude that the screened hit compounds may act as potent inhibitors of HDAC3 and further preclinical and clinical studies may pave the way for developing them as effective therapeutic agents for the treatment of different cancers and neurodegenerative disorders.

Identification of novel topoisomerase II alpha inhibitors by virtual screening, molecular docking, and bioassay

Breast cancer is one of the most common tumors, and its treatment still leaves room for improvement. Topoisomerase II alpha is a potential target for the treatment of human diseases such as breast cancer. In this article, we attempted to discover a novel anticancer drug. We have used the topoisomerase II alpha protein-Homo sapiens (Human) to hierarchically screen the Maybridge database. Based on their docking score, the top hit compounds have been assayed for inhibition in a topoisomerase II pBR322 DNA relaxation assay in vitro. Candidate compound 6 (CP6) was found to have the best inhibitory effect for topoisomerase II among the 20 tested compounds. In addition, CP6 had potent cytotoxicity against eight tested tumor cell lines. At the same time, CP6 was shown to have potential anti-multidrug resistance capabilities. This study identifies CP6, which can contribute to the development of new topoisomerase II inhibitors as anticancer agents.

Prospect of Anterior Gradient 2 homodimer inhibition via repurposing FDA-approved drugs using structure-based virtual screening

Anterior Gradient 2 (AGR2) has recently been reported as a tumor biomarker in various cancers, i.e., breast, prostate and lung cancer. Predominantly, AGR2 exists as a homodimer via a dimerization domain (E60-K64); after it is self-dimerized, it helps FGF2 and VEGF to homo-dimerize and promotes the angiogenesis and the invasion of vascular endothelial cells and fibroblasts. Up till now, no small molecule has been discovered to inhibit the AGR2-AGR2 homodimer. Therefore, the present study was performed to prepare a validated 3D structure of AGR2 by homology modeling and discover a small molecule by screening the FDA-approved drugs library on AGR2 homodimer as a target protein. Thirteen different homology models of AGR2 were generated based on different templates which were narrowed down to 5 quality models sorted by their overall Z-scores. The top homology model based on PDB ID = 3PH9 was selected having the best Z-score and was further assessed by Verify-3D, ERRAT and RAMPAGE analysis. Structure-based virtual screening narrowed down the large library of FDA-approved drugs to ten potential AGR2-AGR2 homodimer inhibitors having FRED score lower than - 7.8 kcal/mol in which the top 5 drugs' binding stability was counter-validated by molecular dynamic simulation. To sum up, the present study prepared a validated 3D structure of AGR2 and, for the first time reported the discovery of 5 FDA-approved drugs to inhibit AGR2-AGR2 homodimer by using structure-based virtual screening. Moreover, the binding of the top 5 hits with AGR2 was also validated by molecular dynamic simulation. A validated 3D structure of Anterior Gradient 2 (AGR2) was prepared by homology modeling, which was used in virtual screening of FDA-approved drugs library for the discovery of prospective inhibitors of AGR2-AGR2 homodimer.

Computational representations of protein-ligand interfaces for structure-based virtual screening

Introduction: Structure-based virtual screening (SBVS) is an essential strategy for hit identification. SBVS primarily uses molecular docking, which exploits the protein-ligand binding mode and associated affinity score for compound ranking. Previous studies have shown that computational representation of protein-ligand interfaces and the later establishment of machine learning models are efficacious in improving the accuracy of SBVS.Areas covered: The authors review the computational methods for representing protein-ligand interfaces, which include the traditional ones that use deliberately designed fingerprints and descriptors and the more recent methods that automatically extract features with deep learning. The effects of these methods on the performance of machine learning models are briefly discussed. Additionally, case studies that applied various computational representations to machine learning are cited with remarks.Expert opinion: It has become a trend to extract binding features automatically by deep learning, which uses a completely end-to-end representation. However, there is still plenty of scope for improvement . The interpretability of deep-learning models, the organization of data management, the quantity and quality of available data, and the optimization of hyperparameters could impact the accuracy of feature extraction. In addition, other important structural factors such as water molecules and protein flexibility should be considered.

Hybrid inhibitors of DNA and HDACs remarkably enhance cytotoxicity in leukaemia cells

Chlorambucil is a nitrogen mustard-based DNA alkylating drug, which is widely used as a front-line treatment of chronic lymphocytic leukaemia (CLL). Despite its widespread application and success for the initial treatment of leukaemia, a majority of patients eventually develop acquired resistance to chlorambucil. In this regard, we have designed and synthesised a novel hybrid molecule, chloram-HDi that simultaneously impairs DNA and HDAC enzymes. Chloram-HDi efficiently inhibits the proliferation of HL-60 and U937 leukaemia cells with GI₅₀ values of 1.24 µM and 1.75 µM, whereas chlorambucil exhibits GI₅₀ values of 21.1 µM and 37.7 µM against HL-60 and U937 leukaemia cells, respectively. The mechanism behind its remarkably enhanced cytotoxicity is that chloram-HDi not only causes a significant DNA damage of leukaemia cells but also downregulates DNA repair protein, Rad52, resulting in the escalation of its DNA-damaging effect. Furthermore, chloram-HDi inhibits HDAC enzymes to induce the acetylation of α-tubulin and histone H3.

HDAC3i-Finder: A Machine Learning-based Computational Tool to Screen for HDAC3 Inhibitors

Histone deacetylase 3 (HDAC3) is a potential drug target for treatment of human diseases such as cancer, chronic inflammation, neurodegenerative diseases and diabetes. Machine learning (ML) as an essential cheminformatics approach has been widely used for QSAR modeling. However, none of them has been applied to HDAC3. To this end, we carefully compiled a set of 1098 compounds from the ChEMBL database that have been assayed against HDAC3 and calculated three different sets of molecular features for each compound, i. e. two-dimensional Mordred descriptors, MACCS keys (166 bits) and Morgan2 fingerprints (1024 bits). Five ML classifiers, i. e. k-Nearest Neighbour (KNN), Support Vector Machine (SVM), Random forest (RF), eXtreme Gradient Boosting (XGBoost) and Deep Neural Network (DNN) were trained on each feature set and optimized for classification. A total of 15 models were generated and carefully compared, among which the best-performing one was the XGBoost model based on the Morgan2 fingerprints, i. e. XGBoost_morgan2. Evaluated on a well-curated benchmarking set named MUBD-HDAC3, this model achieved a high early ROC enrichment (ROCE0.5 %: 41.02). A further retrospective screening of an annotated chemical library in PubChem demonstrated that the best model could identify 8 novel-scaffold HDAC3 inhibitors while assaying only 1 % of the compounds. To make this model accessible for the scientific community, we developed a python GUI application named HDAC3i-Finder to facilitate prospective screening for HDAC3 inhibitors. The source code of HDAC3i-Finder is available at https://github.com/jwxia2014/HDAC3i-Finder.

Two New 1,3,4-Oxadiazoles With Effective Antifungal Activity Against Candida albicans

Candida infections have become a serious public health problem with high mortality rates, especially in immunocompromised patients, since Candida albicans is the major opportunistic pathogen responsible for systemic or invasive candidiasis. Commercially available antifungal agents are restricted and fungal resistance to such drugs has increased; therefore, the development of a more specific antifungal agent is necessary. Using assays for antifungal activity, here we report that two new compounds of 1,3,4-oxadiazoles class (LMM5 and LMM11), which were discovered by in silico methodologies as possible thioredoxin reductase inhibitors, were effective against C. albicans. Both compounds had in vitro antifungal activity with MIC 32 μg/ml. Cytotoxicity in vitro demonstrated that LMM5 and LMM11 were non-toxic in the cell lines evaluated. The kinetic of the time-kill curve suggested a fungistatic profile and showed an inhibitory effect of LMM5 and LMM11 in 12 h that remained for 24 and 36 h, which is better than fluconazole. In the murine systemic candidiasis model by C. albicans, the two compounds significantly reduced the renal and spleen fungal burden. According to the SEM and TEM images, we hypothesize that the mechanism of action of LMM5 and LMM11 is directly related to the inhibition of the enzyme thioredoxin reductase and internally affect the fungal cell. In view of all in vitro and in vivo results, LMM5 and LMM11 are effective therapeutic candidates for the development of new antifungal drugs addressing the treatment of human infections caused by C. albicans.

Cheminformatics techniques in antimalarial drug discovery and development from natural products 1: basic concepts

A large number of natural products, especially those used in ethnomedicine of malaria, have shown varying in vitro antiplasmodial activities. Facilitating antimalarial drug development from this wealth of natural products is an imperative and laudable mission to pursue. However, limited manpower, high research cost coupled with high failure rate during preclinical and clinical studies might militate against the pursuit of this mission. These limitations may be overcome with cheminformatic techniques. Cheminformatics involves the organization, integration, curation, standardization, simulation, mining and transformation of pharmacology data (compounds and bioactivity) into knowledge that can drive rational and viable drug development decisions. This chapter will review the application of cheminformatics techniques (including molecular diversity analysis, quantitative-structure activity/property relationships and Machine learning) to natural products with in vitro and in vivo antiplasmodial activities in order to facilitate their development into antimalarial drug candidates and design of new potential antimalarial compounds.

Rescoring and Linearly Combining: A Highly Effective Consensus Strategy for Virtual Screening Campaigns

The study proposes a novel consensus strategy based on linear combinations of different docking scores to be used in the evaluation of virtual screening campaigns. The consensus models are generated by applying the recently proposed Enrichment Factor Optimization (EFO) method, which develops the linear equations by exhaustively combining the available docking scores and by optimizing the resulting enrichment factors. The performances of such a consensus strategy were evaluated by simulating the entire Directory of Useful Decoys (DUD datasets). In detail, the poses were initially generated by the PLANTS docking program and then rescored by ReScore+ with and without the minimization of the complexes. The so calculated scores were then used to generate the mentioned consensus models including two or three different scoring functions. The reliability of the generated models was assessed by a per target validation as performed by default by the EFO approach. The encouraging performances of the here proposed consensus strategy are emphasized by the average increase of the 17% in the Top 1% enrichment factor (EF) values when comparing the single best score with the linear combination of three scores. Specifically, kinases offer a truly convincing demonstration of the efficacy of the here proposed consensus strategy since their Top 1% EF average ranges from 6.4 when using the single best performing primary score to 23.5 when linearly combining scoring functions. The beneficial effects of this consensus approach are clearly noticeable even when considering the entire DUD datasets as evidenced by the area under the curve (AUC) averages revealing a 14% increase when combining three scores. The reached AUC values compare very well with those reported in literature by an extended set of recent benchmarking studies and the three-variable models afford the highest AUC average.