A lazy learning-based QSAR classification study for screening potential histone deacetylase 8 (HDAC8) inhibitors

Unveiling critical structural features for effective HDAC8 inhibition: a comprehensive study using quantitative read-across structure-activity relationship (q-RASAR) and pharmacophore modeling

Histone deacetylases constitute a group of enzymes that participate in several biological processes. Notably, inhibiting HDAC8 has become a therapeutic strategy for various diseases. The current inhibitors for HDAC8 lack selectivity and target multiple HDACs. Consequently, there is a growing recognition of the need for selective HDAC8 inhibitors to enhance the effectiveness of therapeutic interventions. In our current study, we have utilized a multi-faceted approach, including Quantitative Structure-Activity Relationship (QSAR) combined with Quantitative Read-Across Structure-Activity Relationship (q-RASAR) modeling, pharmacophore mapping, molecular docking, and molecular dynamics (MD) simulations. The developed q-RASAR model has a high statistical significance and predictive ability (Q2F1:0.778, Q2F2:0.775). The contributions of important descriptors are discussed in detail to gain insight into the crucial structural features in HDAC8 inhibition. The best pharmacophore hypothesis exhibits a high regression coefficient (0.969) and a low root mean square deviation (0.944), highlighting the importance of correctly orienting hydrogen bond acceptor (HBA), ring aromatic (RA), and zinc-binding group (ZBG) features in designing potent HDAC8 inhibitors. To confirm the results of q-RASAR and pharmacophore mapping, molecular docking analysis of the five potent compounds (44, 54, 82, 102, and 118) was performed to gain further insights into these structural features crucial for interaction with the HDAC8 enzyme. Lastly, MD simulation studies of the most active compound (54, mapped correctly with the pharmacophore hypothesis) and the least active compound (34, mapped poorly with the pharmacophore hypothesis) were carried out to validate the observations of the studies above. This study not only refines our understanding of essential structural features for HDAC8 inhibition but also provides a robust framework for the rational design of novel selective HDAC8 inhibitors which may offer insights to medicinal chemists and researchers engaged in the development of HDAC8-targeted therapeutics.

The Histone Deacetylase Family: Structural Features and Application of Combined Computational Methods

Histone deacetylases (HDACs) are crucial in gene transcription, removing acetyl groups from histones. They also influence the deacetylation of non-histone proteins, contributing to the regulation of various biological processes. Thus, HDACs play pivotal roles in various diseases, including cancer, neurodegenerative disorders, and inflammatory conditions, highlighting their potential as therapeutic targets. This paper reviews the structure and function of the four classes of human HDACs. While four HDAC inhibitors are currently available for treating hematological malignancies, numerous others are undergoing clinical trials. However, their non-selective toxicity necessitates ongoing research into safer and more efficient 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, and structure-based virtual screening (molecular docking). Moreover, recent developments in the field of molecular dynamics simulations, combined with Poisson-Boltzmann/molecular mechanics generalized Born surface area techniques, have improved the prediction of ligand binding affinity. In this review, we delve into the ways in which these methods have contributed to designing and identifying HDAC inhibitors.

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.

Improving binary crow search algorithm for feature selection

The feature selection (FS) process has an essential effect in solving many problems such as prediction, regression, and classification to get the optimal solution. For solving classification problems, selecting the most relevant features of a dataset leads to better classification accuracy with low training time. In this work, a hybrid binary crow search algorithm (BCSA) based quasi-oppositional (QO) method is proposed as an FS method based on wrapper mode to solve a classification problem. The QO method was employed in tuning the value of flight length in the BCSA which is controlling the ability of the crows to find the optimal solution. To evaluate the performance of the proposed method, four benchmark datasets have been used which are human intestinal absorption, HDAC8 inhibitory activity (IC50), P-glycoproteins, and antimicrobial. Accordingly, the experimental results are discussed and compared against other standard algorithms based on the accuracy rate, the average number of selected features, and running time. The results have proven the robustness of the proposed method relied on the high obtained value of accuracy (84.93–95.92%), G-mean (0.853–0.971%), and average selected features (4.36–11.8) with a relatively low computational time. Moreover, to investigate the effectiveness of the proposed method, Friedman test was used which declared that the performance supremacy of the proposed BCSA-QO with four datasets was very evident against BCSA and CSA by selecting the minimum relevant features and producing the highest accuracy classification rate. The obtained results verify the usefulness of the proposed method (BCSA-QO) in the FS with classification in terms of high classification accuracy, a small number of selected features, and low computational time.

Comparative assessment of favipiravir and remdesivir against human coronavirus NL63 in molecular docking and cell culture models

Human coronavirus NL63 (HCoV-NL63) mainly affects young children and immunocompromised patients, causing morbidity and mortality in a subset of patients. Since no specific treatment is available, this study aims to explore the anti-SARS-CoV-2 agents including favipiravir and remdesivir for treating HCoV-NL63 infection. We first successfully modelled the 3D structure of HCoV-NL63 RNA-dependent RNA polymerase (RdRp) based on the experimentally solved SARS-CoV-2 RdRp structure. Molecular docking indicated that favipiravir has similar binding affinities to SARS-CoV-2 and HCoV-NL63 RdRp with LibDock scores of 75 and 74, respectively. The LibDock scores of remdesivir to SARS-CoV-2 and HCoV-NL63 were 135 and 151, suggesting that remdesivir may have a higher affinity to HCoV-NL63 compared to SARS-CoV-2 RdRp. In cell culture models infected with HCoV-NL63, both favipiravir and remdesivir significantly inhibited viral replication and production of infectious viruses. Overall, remdesivir compared to favipiravir is more potent in inhibiting HCoV-NL63 in cell culture. Importantly, there is no evidence of resistance development upon long-term exposure to remdesivir. Furthermore, combining favipiravir or remdesivir with the clinically used antiviral cytokine interferon-alpha resulted in synergistic effects. These findings provided a proof-of-concept that anti-SARS-CoV-2 drugs, in particular remdesivir, have the potential to be repurposed for treating HCoV-NL63 infection.

QSAR Models for Active Substances against Pseudomonas aeruginosa Using Disk-Diffusion Test Data

Pseudomonas aeruginosa is a Gram-negative bacillus included among the six “ESKAPE” microbial species with an outstanding ability to “escape” currently used antibiotics and developing new antibiotics against it is of the highest priority. Whereas minimum inhibitory concentration (MIC) values against Pseudomonas aeruginosa have been used previously for QSAR model development, disk diffusion results (inhibition zones) have not been apparently used for this purpose in the literature and we decided to explore their use in this sense. We developed multiple QSAR methods using several machine learning algorithms (support vector classifier, K nearest neighbors, random forest classifier, decision tree classifier, AdaBoost classifier, logistic regression and naïve Bayes classifier). We used four sets of molecular descriptors and fingerprints and three different methods of data balancing, together with the “native” data set. In total, 32 models were built for each set of descriptors or fingerprint and balancing method, of which 28 were selected and stacked to create meta-models. In terms of balanced accuracy, the best performance was provided by KNN, logistic regression and decision tree classifier, but the ensemble method had slightly superior results in nested cross-validation.

Discriminations of active from inactive HDAC8 inhibitors Part II: Bayesian classification study to find molecular fingerprints

In continuation of our earlier work (Doi: 10.1080/07391102.2019.1661876), a statistically validated and robust Bayesian model was developed on a large diverse set of HDAC8 inhibitors. The training set comprised of 676 small molecules and 293 compounds were considered as test set molecules. The findings of this analysis will help to explore some major directions regarding the HDAC8 inhibitor designing approach. Acrylamide (G1-G3, G9), N-substituted 2-phenylimidazole (G4-G8, G9, G12-G13, G16-G19), benzimidazole (G10-G11), piperidine substituted pyrrole (G13-G14) groups, alkyl/aryl amide (G15) and aryloxy carboxamide (G20) fingerprints were found to play a crucial role in HDAC8 inhibitory activity whereas -CH-N=CH- (B1, B4-B6, B14) motif, benzamide (B2-B3, B9-B13, B16-B17) groups and heptazepine (B7-B8, B15, B18-B20) group were found to influence negatively the HDAC8 inhibitory activity. The importance of such fingerprints was further validated by the HDAC8 enzyme and related inhibitor interactions at the receptor level. These results are in close agreement with those of our previous work that validate each other. Moreover, this comparative learning may enrich future endeavours regarding the designing strategy of HDAC8 inhibitors.

Virtual Screening Meets Deep Learning

BACKGROUND

Automated compound testing is currently the de facto standard method for drug screening, but it has not brought the great increase in the number of new drugs that was expected. Computer- aided compounds search, known as Virtual Screening, has shown the benefits to this field as a complement or even alternative to the robotic drug discovery. There are different methods and approaches to address this problem and most of them are often included in one of the main screening strategies. Machine learning, however, has established itself as a virtual screening methodology in its own right and it may grow in popularity with the new trends on artificial intelligence.

OBJECTIVE

This paper will attempt to provide a comprehensive and structured review that collects the most important proposals made so far in this area of research. Particular attention is given to some recent developments carried out in the machine learning field: the deep learning approach, which is pointed out as a future key player in the virtual screening landscape.

Histone deacetylase 8 (HDAC8) and its inhibitors with selectivity to other isoforms: An overview

The histone deacetylases (HDACs) enzymes provided crucial role in transcriptional regulation of cells through deacetylation of nuclear histone proteins. Discoveries related to the HDAC8 enzyme activity signified the importance of HDAC8 isoform in cell proliferation, tumorigenesis, cancer, neuronal disorders, parasitic/viral infections and other epigenetic regulations. The pan-HDAC inhibitors can confront these conditions but have chances to affect epigenetic functions of other HDAC isoforms. Designing of selective HDAC8 inhibitors is a key feature to combat the pathophysiological and diseased conditions involving the HDAC8 activity. This review is concerned about the structural and positional aspects of HDAC8 in the HDAC family. It also covers the contributions of HDAC8 in the pathophysiological conditions, a preliminary discussion about the recent scenario of HDAC8 inhibitors. This review might help to deliver the structural, functional and computational information in order to identify and design potent and selective HDAC8 inhibitors for target specific treatment of diseases involving HDAC8 enzymatic activity.

Combiphore (Structure and Ligand Based Pharmacophore) - Approach for the Design of GPR40 Modulators in the Management of Diabetes

BACKGROUND

Pharmacophore mapping and molecular docking can be synergistically integrated to improve the drug design and discovery process. A rational strategy, combiphore approach, derived from the combined study of Structure and Ligand based pharmacophore has been described to identify novel GPR40 modulators.

METHODS

DISCOtech module from Discovery studio was used for the generation of the Structure and Ligand based pharmacophore models which gave hydrophobic aromatic, ring aromatic and negative ionizable as essential pharmacophoric features. The generated models were validated by screening active and inactive datasets, GH scoring and ROC curve analysis. The best model was exposed as a 3D query to screen the hits from databases like GLASS (GPCR-Ligand Association), GPCR SARfari and Mini-Maybridge. Various filters were applied to retrieve the hit molecules having good drug-like properties. A known protein structure of hGPR40 (pdb: 4PHU) having TAK-875 as ligand complex was used to perform the molecular docking studies; using SYBYL-X 1.2 software.

RESULTS AND CONCLUSION

Clustering both the models gave RMSD of 0.89. Therefore, the present approach explored the maximum features by combining both ligand and structure based pharmacophore models. A common structural motif as identified in combiphore for GPR40 modulation consists of the para-substituted phenyl propionic acid scaffold. Therefore, the combiphore approach, whereby maximum structural information (from both ligand and biological protein) is explored, gives maximum insights into the plausible protein-ligand interactions and provides potential lead candidates as exemplified in this study.

Novel HDAC8 inhibitors: A multi-computational approach

Abnormal HDAC function triggers irregular gene transcription that hampers the essential cellular activities leading to tumour activation and progression. HDAC inhibition has, therefore, been reported as a potential target for cancer treatment. In the present study, a sequential computational framework was carried out to discover newer lead compounds, namely HDAC8 inhibitors for cancer therapy. Pharmacophoric hypotheses were generated based on hydroxamic acid derivatives reported earlier for HDAC inhibition. The model AAADR.122, demonstrated statistical significance (r² = 0.93, Q² = 0.81) and proved robust on validation with a cross-validated correlation coefficient of 0.89. It was utilized to arrive at novel hits through a virtual screening workflow. The specificity of the process was enhanced further by analysing the crucial interactions of the ligands with key catalytic residues, achieved by induced fit docking (PDB ID: 1T64). On assessment, the filtered leads displayed optimal drug like features. Investigations using density functional theory (DFT) also facilitated the recognition of molecular spots in the leads beneficial for HDAC8 interaction. Overall, two leads were proposed for HDAC8 inhibition with potential anti-cancer activity.

QSAR models for predicting the toxicity of piperidine derivatives against Aedes aegypti

QSAR models are proposed for predicting the toxicity of 33 piperidine derivatives against Aedes aegypti. From 2D topological descriptors, calculated with the PaDEL software, ordinary least squares multilinear regression (OLS-MLR) treatment from the QSARINS software and machine learning and related approaches including linear and radial support vector machine (SVM), projection pursuit regression (PPR), radial basis function neural network (RBFNN), general regression neural network (GRNN) and k-nearest neighbours (k-NN), led to four-variable models. Their robustness and predictive ability were evaluated through both internal and external validation. Determination coefficients (r²) greater than 0.85 on the training sets and 0.8 on the test sets were obtained with OLS-MLR and linear SVM. They slightly outperform PPR, radial SVM and RBFNN, whereas GRNN and k-NN showed lower performance. The easy availability of the involved structural descriptors and the simplicity of the MLR model make the corresponding model attractive at an exploratory level for proposing, from this limited dataset, guidelines in the design of new potentially active molecules.

Quantitative structure-activity relationship analysis and virtual screening studies for identifying HDAC2 inhibitors from known HDAC bioactive chemical libraries

Histone deacetylases (HDAC) are emerging as promising targets in cancer, neuronal diseases and immune disorders. Computational modelling approaches have been widely applied for the virtual screening and rational design of novel HDAC inhibitors. In this study, different machine learning (ML) techniques were applied for the development of models that accurately discriminate HDAC2 inhibitors form non-inhibitors. The obtained models showed encouraging results, with the global accuracy in the external set ranging from 0.83 to 0.90. Various aspects related to the comparison of modelling techniques, applicability domain and descriptor interpretations were discussed. Finally, consensus predictions of these models were used for screening HDAC2 inhibitors from four chemical libraries whose bioactivities against HDAC1, HDAC3, HDAC6 and HDAC8 have been known. According to the results of virtual screening assays, structures of some hits with pair-isoform-selective activity (between HDAC2 and other HDACs) were revealed. This study illustrates the power of ML-based QSAR approaches for the screening and discovery of potent, isoform-selective HDACIs.

Targeting histone deacetylase 8 as a therapeutic approach to cancer and neurodegenerative diseases

Histone deacetylase 8 (HDAC8), a unique class I zinc-dependent HDAC, is an emerging target in cancer and other diseases. Its substrate repertoire extends beyond histones to many nonhistone proteins. Besides being a deacetylase, HDAC8 also mediates signaling via scaffolding functions. Aberrant expression or deregulated interactions with transcription factors are critical in HDAC8-dependent cancers. Many potent HDAC8-selective inhibitors with cellular activity and anticancer effects have been reported. We present HDAC8 as a druggable target and discuss inhibitors of different chemical scaffolds with cellular effects. Furthermore, we review HDAC8 activators that revert activity of mutant enzymes. Isotype-selective HDAC8 targeting in patients with HDAC8-relevant cancers is challenging, however, is promising to avoid adverse side effects as observed with pan-HDAC inhibitors.