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

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.

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.

Apigenin and its combination with Vorinostat induces apoptotic-mediated cell death in TNBC by modulating the epigenetic and apoptotic regulators and related miRNAs

Triple-negative breast cancer (TNBC) is a metastatic disease and a formidable treatment challenge as it does not respond to existing therapies. Epigenetic regulators play a crucial role in the progression and metastasis by modulating the expression of anti-apoptotic, pro-apoptotic markers and related miRNAs in TNBC cells. We have investigated the anti-TNBC potential of dietary flavonoid 'Apigenin' and its combination with Vorinostat on MDA-MB-231 cells. At Apigenin generated ROS, inhibited cell migration, arrested the cell cycle at subG0/G1 phases, and induced apoptotic-mediated cell death. Apigenin reduced the expression of the class-I HDACs at the transcriptomic and proteomic levels. In the immunoblotting study, Apigenin has upregulated pro-apoptotic markers and downregulated anti-apoptotic proteins. Apigenin inhibited the enzymatic activity of HDAC/DNMT and increased HAT activity. Apigenin has manifested its effect on miRNA expression by upregulating the tumor-suppressor miR-200b and downregulation oncomiR-21. Combination study reduced the growth of TNBC cells synergistically by modulating the expression of epigenetic and apoptotic regulators. Molecular docking and MD simulations explored the mechanism of catalytic inhibition of HDAC1 and HDAC3 and supported the in-vitro studies. The overall studies demonstrated an anti-TNBC potential of Apigenin and may help to design an effective strategy to treat metastatic phenotype of TNBC.

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.

Design, synthesis, and histone deacetylase inhibition study of novel 4-(2-aminoethyl) phenol derivatives

Histone deacetylases (HDACs) have been identified as promising targets for anticancer treatment. The study demonstrates virtual screening, molecular docking, and synthesis of 4-(2-aminoethyl) phenol derivatives as HDAC inhibitors. The virtual screening and molecular docking analysis led to the identification of 10 representative compounds, which were evaluated based on their drug-like properties. The results demonstrated that these compounds effectively interacted with the active site pocket of HDAC 3 through π-stacking, Zn2+ coordination, hydrogen bonding, and hydrophobic interactions with catalytic residues. Furthermore, a series of 4-(2-aminoethyl) phenol derivatives were synthesized, and their HDAC inhibitory activity was evaluated. Compounds 18 and 20 showed significant HDAC inhibitory activity of 64.94 ± 1.17% and 52.45 ± 1.45%, respectively, compared to the solvent control. The promising results of this study encourage further research on 4-(2-aminoethyl) phenol derivatives and may provide significant insight into the design of novel small molecule HDAC inhibitors to fight against target-specific malignancies of chronic obstructive pulmonary disease and nonsmall cell lung cancer in the future.

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 Review on Molecular Docking on HDAC Isoforms: Novel Tool for Designing Selective Inhibitors

Research into histone deacetylases (HDACs) has experienced a remarkable surge in recent years. These enzymes are key regulators of several fundamental biological processes, often associated with severe and potentially fatal diseases. Inhibition of their activity represents a promising therapeutic approach and a prospective strategy for the development of new therapeutic agents. A critical aspect of their inhibition is to achieve selectivity in terms of enzyme isoforms, which is essential to improve treatment efficacy while reducing undesirable pleiotropic effects. The development of computational chemistry tools, particularly molecular docking, is greatly enhancing the precision of designing molecules with inherent potential for specific activity. Therefore, it was considered necessary to review the molecular docking studies conducted on the major isozymes of the enzyme in order to identify the specific interactions associated with each selective HDAC inhibitor. In particular, the most critical isozymes of HDAC (1, 2, 3, 6, and 8) have been thoroughly investigated within the scope of this review.

Toxicity prediction and analysis of flavonoid apigenin as a histone deacetylase inhibitor: an in-silico approach

Occurrence of cancer is driving up on a global scale that exerts greater implications on the physical, psychological and economic stability of the human population. In the present context, numerous research studies are being conducted to explore and discover the drug molecule as an anticancer agent. Diverse scales of flavonoids entail the human diet, and they displayed prospective curative effects against an array of ailments. Among different categories of flavonoids, apigenin a trihydroxy flavone has been proven to have various pharmacological effects. Molecular docking is a key tool in structural molecular biology and computer assisted drug design. In this study, HDAC inhibitory action of apigenin and its probable toxicity was assessed by docking study using Auto dock platform. Molecular dynamics simulation was done by using iMODS server for elucidating the stability of the receptor-ligand complex. Toxicity predictions were evaluated by using tools such as CarcinoPred for carcinogenicity study, pkCSM for ADMET analysis, ProTox-II for rodent oral toxicity, lazar for estimating mutagenicity, BOILED Egg plot analysis for examining the gastrointestinal absorption and blood brain permeability, PASS prediction to identify the various biological functions and DruLiTo program to compute the drug likeness property.

Graphical abstract

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.

Epigenetic regulation exerted by Caliphruria subedentata and galantamine: an in vitro and in silico approach for mimic Alzheimer's disease

At the interface between genes and environment, epigenetic mechanisms, including DNA methylation and histone modification, regulate neurogenic processes such as differentiation, proliferation, and maturation of neural stem cells. However, these mechanisms are altered in Alzheimer's disease (AD), a neurodegenerative condition that mainly affects older adults. Since epigenetic mechanisms are known to be reversible, a number of molecules from natural sources are being studied as epigenetic regulators in AD. Recently, in vitro and in silico studies have shown that C. subedentata and its alkaloids modulated neurotoxicity. However, studies exploring the epigenetic activity of these alkaloids are limited. We conducted a set of bioassays to evaluate neuronal differentiation and the sensitivity of undifferentiated SH-SY5 cells against a neurotoxic stimulus. In addition, we analyzed the methylation profiles in genes such as APP, PSI, and BACE1 due to their role in amyloid processing. Docking and molecular dynamic analysis were used to explore the effect exerted by C. subedentata alkaloids on the regulation of histone deacetylases (HDAC2, HDAC3 and HDAC7). The results demonstrated that C. subedentata and galantamine induce neuronal differentiation and protect the undifferentiated SH-SY5Y cells against Aβ(1-42)-induced neurotoxicity. The methylation profiles of the studied genes show no statistically significant differences between C. subedentata, galantamine. However, these findings should be interpreted with caution, since small changes in methylation promoters in the brain could not be easily detected. Results from in silico approaches describe for the first time the potential promissing epigenetic effects of galantamine by regulating HDAC3 and HDAC7 modification.Communicated by Ramaswamy H. Sarma.

Repurposing of neprilysin inhibitor 'sacubitrilat' as an anti-cancer drug by modulating epigenetic and apoptotic regulators

Modifications in the epigenetic landscape have been considered a hallmark of cancer. Histone deacetylation is one of the crucial epigenetic modulations associated with the aggressive progression of various cancer subtypes. Herein, we have repurposed the neprilysin inhibitor sacubitrilat as a potent anticancer agent using in-silico protein-ligand interaction profiler (PLIP) analysis, molecular docking, and in vitro studies. The screening of PLIP profiles between vorinostat/panobinostat and HDACs/LTA4H followed by molecular docking resulted in five (Sacubitrilat, B65, BDS, BIR, and NPV) FDA-approved, experimental and investigational drugs. Sacubitrilat has demonstrated promising anticancer activity against colorectal cancer (SW-480) and triple-negative breast cancer (MDA-MB-231) cells, with IC₅₀ values of 14.07 μg/mL and 23.02 μg/mL, respectively. FACS analysis revealed that sacubitrilat arrests the cell cycle at the G0/G1 phase and induces apoptotic-mediated cell death in SW-480 cells. In addition, sacubitrilat inhibited HDAC isoforms at the transcriptomic level by 0.7-0.9 fold and at the proteomic level by 0.5-0.6 fold as compared to the control. Sacubitrilat increased the protein expression of tumor-suppressor (p53) and pro-apoptotic makers (Bax and Bid) by 0.2-2.5 fold while decreasing the expression of anti-apoptotic Bcl2 and Nrf2 proteins by 0.2-0.5 fold with respect to control. The observed cleaved PARP product indicates that sacubitrilat induces apoptotic-mediated cell death. This study may pave the way to identify the anticancer potential of sacubitrilat and can be explored in human clinical trials.

Structure based virtual screening, molecular dynamic simulation to identify the oxadiazole derivatives as inhibitors of Enterococcus D-Ala-D-Ser ligase for combating vancomycin resistance

Vancomycin resistance in enterococci mainly arises due to alteration in terminal peptidoglycan dipeptide. A comprehensive structural analysis for substrate specificity of dipeptide modifying d-Alanine: d-Serine ligase (Ddls) is essential to screen its inhibitors for combating vancomycin resistance. In this study modeled 3D structure of EgDdls from E. gallinarum was used for structure based virtual screening (SBVS) of oxadiazole derivatives. Initially, fifteen oxadiazole derivatives were identified as inhibitors at the active site of EgDdls from PubChem database. Further, four EgDdls inhibitors were evaluated using pharmacokinetic profile and molecular docking. The results of molecular docking showed that oxadiazole inhibitors could bind preferentially at ATP binding pocket with the lowest binding energy. Further, molecular dynamics simulation results showed stable behavior of EgDdls in complex with screened inhibitors. The residues Phe172, Lys174, Glu217, Phe292, and Asn302 of EgDdls were mainly involved in interactions with screened inhibitors. Furthermore, MM-PBSA calculation showed electrostatic and van der Waals interactions mainly contribute to overall binding energy. The PCA analysis showed motion of central domain and omega loop of EgDdls. This is involved in the formation of native dipeptide and stabilized after binding of 2-(1-(Ethylsulfonyl) piperidin-4-yl)-5-(furan-2-yl)-1,3,4-oxadiazole, which could be reason for the inhibition of EgDdls. Hence, in this study we have screened inhibitors of EgDdls which could be useful to alleviate the vancomycin resistance problem in enterococci, involved in hospital-acquired infections, especially urinary tract infections (UTI).

Identification of Activated Cdc42-Associated Kinase Inhibitors as Potential Anticancer Agents Using Pharmacoinformatic Approaches

BACKGROUND

Activated Cdc42-associated kinase (ACK1) is essential for numerous cellular functions, such as growth, proliferation, and migration. ACK1 signaling occurs through multiple receptor tyrosine kinases; therefore, its inhibition can provide effective antiproliferative effects against multiple human cancers. A number of ACK1-specific inhibitors were designed and discovered in the previous decade, but none have reached the clinic. Potent and selective ACK1 inhibitors are urgently needed.

METHODS

In the present investigation, the pharmacophore model (PM) was rationally built utilizing two distinct inhibitors coupled with ACK1 crystal structures. The generated PM was utilized to screen the drug-like database generated from the four chemical databases. The binding mode of pharmacophore-mapped compounds was predicted using a molecular docking (MD) study. The selected hit-protein complexes from MD were studied under all-atom molecular dynamics simulations (MDS) for 500 ns. The obtained trajectories were ranked using binding free energy calculations (ΔG kJ/mol) and Gibb's free energy landscape.

RESULTS

Our results indicate that the three hit compounds displayed higher binding affinity toward ACK1 when compared with the known multi-kinase inhibitor dasatinib. The inter-molecular interactions of Hit1 and Hit3 reveal that compounds form desirable hydrogen bond interactions with gatekeeper T205, hinge region A208, and DFG motif D270. As a result, we anticipate that the proposed scaffolds might help in the design of promising selective ACK1 inhibitors.

Ligands-induced open-close conformational change during DapE catalysis: Insights from molecular dynamics simulations

The microbial enzyme DapE plays a critical role in the lysine biosynthetic pathway and is considered as a potentially safe antibiotic target. In this study, atomistic simulations are employed to identify the modes of essential dynamics that define the conformational response of the enzyme to ligand binding and unbinding. The binding modes and the binding affinities of the products to the DapE enzyme are estimated from the MM-PBSA method, and the residues contributing to the ligand binding are identified. Various structural analyses and the principal component analysis of the molecular dynamics trajectories reveal that the removal of products from the active site causes a significant change in the overall enzyme structure. Both Cartesian and dihedral principal component analyses are used to characterize the structural changes in terms of domain unfolding and domain twisting motions. In the most dominant mode, that is, the domain unfolding motion, the two catalytic domains move away from the two dimerization domains of the dimeric enzyme, representing a closed-to-open conformational change. The conformational changes are initiated by the coordinated movement of three loops (Asp75-Pro82, Gly240-Asn244, and Thr347-Glu353) that trigger a domain-level movement. From multiple short trajectories, the time constant associated with the domain opening motion is estimated as 43.6 ns. Physiologically, this close-to-open conformational change is essential for the regeneration of the initial state of the enzyme for the subsequent cycle of catalytic action and provides the apo enzyme enough flexibility for efficient substrate binding.

Pharmacophore-Oriented Identification of Potential Leads as CCR5 Inhibitors to Block HIV Cellular Entry

Cysteine-cysteine chemokine receptor 5 (CCR5) has been discovered as a co-receptor for cellular entry of human immunodeficiency virus (HIV). Moreover, the role of CCR5 in a variety of cancers and various inflammatory responses was also discovered. Despite the fact that several CCR5 antagonists have been investigated in clinical trials, only Maraviroc has been licensed for use in the treatment of HIV patients. This indicates that there is a need for novel CCR5 antagonists. Keeping this in mind, the present study was designed. The active CCR5 inhibitors with known IC₅₀ value were selected from the literature and utilized to develop a ligand-based common feature pharmacophore model. The validated pharmacophore model was further used for virtual screening of drug-like databases obtained from the Asinex, Specs, InterBioScreen, and Eximed chemical libraries. Utilizing computational methods such as molecular docking studies, molecular dynamics simulations, and binding free energy calculation, the binding mechanism of selected inhibitors was established. The identified Hits not only showed better binding energy when compared to Maraviroc, but also formed stable interactions with the key residues and showed stable behavior throughout the 100 ns MD simulation. Our findings suggest that Hit1 and Hit2 may be potential candidates for CCR5 inhibition, and, therefore, can be considered for further CCR5 inhibition programs.

Ameliorative Sexual Behavior and Phosphodiesterase-5 Inhibitory Effects of Spondias mangifera Fruit Extract in Rodents: In Silico, In Vitro, and In Vivo Study

The ethanolic extracts of Spondias mangifera fruit (SMFE) were evaluated for aphrodisiac activity. The in-vitro phosphodiesterase-5 (PDE-5) inhibition was assessed based on in-silico molecular docking and simulation studies. In addition, the in-vivo sexual behavior was analyzed in the form of mount (MF, ML), intromission (IF, IL), and ejaculation (EF, EL) frequencies and latencies to validate the in-vitro results. Some biochemical parameters, including PDE-5, nitric oxide, and testosterone, were also observed. The above extract constituted β-amyrin, β-sitosterol, and oleanolic acid and showed tremendous binding with phosphodiesterase-5 and sildenafil. Both the sildenafil and ethanolic extracts (200 and 400 mg/kg/d bodyweight) significantly (p < 0.1, p < 0.05) increased MF, IF, and EF, respectively. In contrast, ML and IL significantly (p < 0.1) decreased, and EL significantly (p < 0.1) increased compared with a normal group of animals. The ethanolic extracts (200 and 400 mg/kg/d bodyweight) and sildenafil further significantly (p < 0.05, p < 0.1) diminished PDE-5 activity significantly (p < 0.05, p < 0.1) and enhanced nitric oxide and testosterone levels, as compared with normal rodents. Therefore, the S. mangifera ethanolic extract might be a valuable alternate aphrodisiac for erectile dysfunction.