3D-QSAR (CoMFA, CoMSIA) and Molecular Docking Studies on Histone Deacetylase 1 Selective Inhibitors

Identification of structural fingerprints among natural inhibitors of HDAC1 to accelerate nature-inspired drug discovery in cancer epigenetics

Histone deacetylase 1 (HDAC1), a class I HDAC enzyme, is crucial for histone modification. Currently, it is emerged as one of the important biological targets for designing small molecule drugs through cancer epigenetics. Along with synthetic inhibitors different natural inhibitors are showing potential HDAC1 inhibitions. In order to gain insights into the relationship between the molecular structures of the natural inhibitors and HDAC1, different molecular modelling techniques (Bayesian classification, recursive partitioning, molecular docking and molecular dynamics simulations) have been applied on a dataset of 155 HDAC1 nature-inspired inhibitors with diverse scaffolds. The Bayesian study showed acceptable ROC values for both the training set and test sets. The Recursive partitioning study produced decision tree 1 with 6 leaves. Further, molecular docking study was processed for generating the protein ligand complex which identified some potential amino acid residues such as F205, H28, L271, P29, F150, Y204 for the binding interactions in case of natural inhibitors. Stability of these HDAC1-natutal inhibitors complexes has been also evaluated by molecular dynamics simulation study. The current modelling study is an attempt to get a deep insight into the different important structural fingerprints among different natural compounds modulating HDAC1 inhibition.Communicated by Ramaswamy H. Sarma.

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.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule's binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R2 = 0.8521, Q2LOO = 0.8438, and Q2LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: -8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies.Communicated by Ramaswamy H. Sarma.

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.

A quantitative structural analysis of AR-42 derivatives as HDAC1 inhibitors for the identification of promising structural contributors

Alteration and abnormal epigenetic mechanisms can lead to the aberration of normal biological functions and the occurrence of several diseases. The histone deacetylase (HDAC) family of enzymes is one of the prime regulators of epigenetic functions modifying the histone proteins, and thus, regulating epigenetics directly. HDAC1 is one of those HDACs which have important contributions to cellular epigenetics. The abnormality of HDAC is correlated to the occurrence, progression, and poor prognosis in several disease conditions namely neurodegenerative disorders, cancer cell proliferation, metastasis, chemotherapy resistance, and survival in various cancers. Therefore, the progress of potent and effective HDAC1 inhibitors is one of the prime approaches to combat such diseases. In this study, both regression and classification-based molecular modelling studies were conducted on some AR-42 derivatives as HDAC1 inhibitors to elucidate the crucial structural aspects that are responsible for regulating their biological responses. This study revealed that the molecular polarizability, van der Waals volume, the presence of aromatic rings as well as the higher number of hydrogen bond acceptors might affect prominently their inhibitory activity and might be responsible for proper fitting and interactions at the HDAC1 active site to pertain effective inhibition.

Computational modelling studies of some 1,3-thiazine derivatives as anti-influenza inhibitors targeting H1N1 neuraminidase via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Background

Influenza virus disease remains one of the most contagious diseases that aided the deaths of many patients, especially in this COVID-19 pandemic era. Recent discoveries have shown that the high prevalence of influenza and SARS-CoV-2 coinfection can rapidly increase the death rate of patients. Hence, it became necessary to search for more potent inhibitors for influenza disease therapy. The present study utilized some computational modeling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions of some 1,3-thiazine derivatives as inhibitors of influenza neuraminidase (NA).

Results

The 2D-QSAR modeling results showed GFA-MLR ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9192, Q2 = 0.8767, R2adj = 0.8991, RMSE = 0.0959, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8943, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7745) and GFA-ANN ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9227, Q2 = 0.9212, RMSE = 0.0940, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8831, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7763) models with the computed descriptors as ATS7s, SpMax5_Bhv, nHBint6, and TDB9m for predicting the NA inhibitory activities of compounds which have passed the global criteria of accepting QSAR model. The 3D-QSAR modeling was carried out based on the comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA). The CoMFA_ES ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9620, Q2 = 0.643) and CoMSIA_SED ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.8770, Q2 = 0.702) models were found to also have good and reliable predicting ability. The compounds were also virtually screened based on their binding scores via molecular docking simulations with the active site of the NA (H1N1) target receptor which also confirms their resilient potency. Four potential lead compounds (4, 7, 14, and 15) with the relatively high inhibitory rate (> 50%) and docking (> − 6.3 kcal/mol) scores were identified as the possible lead candidates for in silico exploration of improved anti-influenza agents.

Conclusion

The drug-likeness and ADMET predictions of the lead compounds revealed non-violation of Lipinski’s rule and good pharmacokinetic profiles as important guidelines for rational drug design. Hence, the outcome of this research set a course for the in silico design and exploration of novel NA inhibitors with improved potency.

In-silico modelling studies of 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase inhibitors via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Influenza virus disease is one of the most infectious diseases responsible for many human deaths, and the high mutability of the virus causes drug resistance effects in recent times. As such, it became necessary to explore more inhibitors that could avert future influenza pandemics. The present research utilized some in-silico modelling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions on some 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase (NA) inhibitors. The 2D-QSAR modelling results revealed GFA-MLR (R train 2 =0.8414, Q2 = 0.7680) and GFA-ANN (R train ​ 2 =0.8754, Q2 = 0.8753) models with the most relevant descriptors (MATS3i, SpMax5_Bhe, minsOH and VE3_D) for predicting the inhibitory activities of the molecules which has passed the global criteria of accepting QSAR models. The results of the 3D-QSAR modelling results showed that CoMFA_ES (R train ​ 2 =0.9030, Q2 = 0.5390) and CoMSIA_EA (R train 2 =0.880, Q2 = 0.547) models are having good predicting ability among other developed models. The molecules were virtually screened via molecular docking simulation with the active site of NA protein receptor (pH1N1) which confirms their resilient potency when compared with zanamivir standard drug. Molecule 11 as the most potent molecule formed more H-bond interactions with the key residues such as TRP178, ARG152, ARG292, ARG371, and TYR406 that triggered the catalytic reactions for NA inhibition. Furthermore, six (6) molecules (9, 10, 11, 17, 22, and 31) with relatively high inhibitory activities and docking scores were identified as the possible leads for in-silico exploration of novel NA inhibitors. The drug-likeness and ADMET predictions of the lead molecules revealed non-violation of Lipinski's rule and good pharmacokinetic profiles respectively, which are important guidelines for rational drug design. Hence, the outcome of this study overlaid a solid foundation for the in-silico design and exploration of novel NA inhibitors with improved potency.

Classification models and SAR analysis on HDAC1 inhibitors using machine learning methods

Histone deacetylase (HDAC) 1, a member of the histone deacetylases family, plays a pivotal role in various tumors. In this study, we collected 7313 human HDAC1 inhibitors with bioactivities to form a dataset. Then, the dataset was divided into a training set and a test set using two splitting methods: (1) Kohonen's self-organizing map and (2) random splitting. The molecular structures were represented by MACCS fingerprints, RDKit fingerprints, topological torsions fingerprints and ECFP4 fingerprints. A total of 80 classification models were built by using five machine learning methods, including decision tree (DT), random forest, support vector machine, eXtreme Gradient Boosting and deep neural network. Model 15A_2 built by the XGBoost algorithm based on ECFP4 fingerprints showed the best performance, with an accuracy of 88.08% and an MCC value of 0.76 on the test set. Finally, we clustered the 7313 HDAC1 inhibitors into 31 subsets, and the substructural features in each subset were investigated. Moreover, using DT algorithm we analyzed the structure-activity relationship of HDAC1 inhibitors. It may conclude that some substructures have a significant effect on high activity, such as N-(2-amino-phenyl)-benzamide, benzimidazole, AR-42 analogues, hydroxamic acid with a middle chain alkyl and 4-aryl imidazole with a midchain of alkyl whose α carbon is chiral.

3D-QSAR-Based Pharmacophore Modeling, Virtual Screening, and Molecular Docking Studies for Identification of Tubulin Inhibitors with Potential Anticancer Activity

In this study, we aimed to develop a pharmacophore-based three-dimensional quantitative structure activity relationship (3D-QSAR) for a set including sixty-two cytotoxic quinolines (1-62) as anticancer agents with tubulin inhibitory activity. A total of 279 pharmacophore hypotheses were generated based on the survival score to build QSAR models. A six-point pharmacophore model (AAARRR.1061) was identified as the best model which consisted of three hydrogen bond acceptors (A) and three aromatic ring (R) features. The model showed a high correlation coefficient (R 2 = 0.865), cross-validation coefficient (Q 2 = 0.718), and F value (72.3). The best pharmacophore model was then validated by the Y-Randomization test and ROC-AUC analysis. The generated 3D contour maps were used to reveal the structure activity relationship of the compounds. The IBScreen database was screened against AAARRR.1061, and after calculating ADMET properties, 10 compounds were selected for further docking study. Molecular docking analysis showed that compound STOCK2S-23597 with the highest docking score (-10.948 kcal/mol) had hydrophobic interactions and can form four hydrogen bonds with active site residues.

QSAR Modeling, Molecular Docking and Molecular Dynamics Simulations Studies of Lysine-Specific Demethylase 1 (LSD1) Inhibitors as Anticancer Agents

BACKGROUND

Histone Lysine Demetylases1 (LSD1) is a promising medication to treat cancer, which plays a crucial role in epigenetic modulation of gene expression. Inhibition of LSD1with small molecules has emerged as a vital mechanism to treat cancer.

OBJECTIVE

In the present research, molecular modeling investigations, such as CoMFA, CoMFA-RF, CoMSIA and HQSAR, molecular docking and Molecular Dynamics (MD) simulations were carried out on some tranylcypromine derivatives as LSD1 inhibitors.

METHODS

The QSAR models were carried out on a series of Tranylcypromine derivatives as data set via the SYBYL-X2.1.1 program. Molecular docking and MD simulations were carried out by the MOE software and the SYBYL program, respectively. The internal and external predictability performances related to the generated models for these LSD1 inhibitors were justified by evaluating cross-validated correlation coefficient (q²), noncross- validated correlation coefficient (r²_ncv) and predicted correlation coefficient (r²_pred) of the training and test set molecules, respectively.

RESULTS

The CoMFA (q², 0.670; r²_ncv, 0.930; r²_pred, 0.968), CoMFA-RF (q², 0.694; r²_ncr, 0.926; r²_pred, 0.927), CoMSIA (q², 0.834; r²_ncv, 0.956; r²_pred, 0.958) and HQSAR models (q², 0.854; r²_ncv, 0.900; r²_pred, 0.728) for training as well as the test set of LSD1 inhibition resulted in significant findings.

CONCLUSION

These QSAR models were found to be perfect and strong with better predictability. Contour maps of all models were generated and it was proven by molecular docking studies and molecular dynamics simulation that the hydrophobic, electrostatic and hydrogen bonding fields are crucial in these models for improving the binding affinity and determining the structure-activity relationship. These theoretical results are possibly beneficial to design new strong LSD1 inhibitors with enhanced activity to treat cancer.