Zinc chelation with hydroxamate in histone deacetylases modulated by water access to the linker binding channel

Density functional theory (DFT) studies in HDAC-based chemotherapeutics: Current findings, case studies and future perspectives

Density Functional Theory (DFT) is a quantum chemical computational method used to predict and analyze the electronic properties of atoms, molecules, and solids based on the density of electrons rather than wavefunctions. It provides insights into the structure, bonding, and behavior of different molecules, including those involved in the development of chemotherapeutic agents, such as histone deacetylase inhibitors (HDACis). HDACs are a wide group of metalloenzymes that facilitate the removal of acetyl groups from acetyl-lysine residues situated in the N-terminal tail of histones. Abnormal HDAC recruitment has been linked to several human diseases, especially cancer. Therefore, it has been recognized as a prospective target for accelerating the development of anticancer therapies. Researchers have studied HDACs and its inhibitors extensively using a combination of experimental methods and diverse in-silico approaches such as machine learning and quantitative structure-activity relationship (QSAR) methods, molecular docking, molecular dynamics, pharmacophore mapping, and more. In this context, DFT studies can make significant contribution by shedding light on the molecular properties, interactions, reaction pathways, transition states, reactivity and mechanisms involved in the development of HDACis. This review attempted to elucidate the scope in which DFT methodologies may be used to enhance our comprehension of the molecular aspects of HDAC inhibitors, aiding in the rational design and optimization of these compounds for therapeutic applications in cancer and other ailments. The insights gained can guide experimental efforts toward developing more potent and selective HDAC inhibitors.

Metalloenzyme Inhibitors against Zoonotic Infections: Focus on Leishmania and Schistosoma

The term "zoonosis" denotes diseases transmissible among vertebrate animals and humans. These diseases constitute a significant public health challenge, comprising 61% of human pathogens and causing an estimated 2.7 million deaths annually. Zoonoses not only affect human health but also impact animal welfare and economic stability, particularly in low- and middle-income nations. Leishmaniasis and schistosomiasis are two important neglected tropical diseases with a high prevalence in tropical and subtropical areas, imposing significant burdens on affected regions. Schistosomiasis, particularly rampant in sub-Saharan Africa, lacks alternative treatments to praziquantel, prompting concerns regarding parasite resistance. Similarly, leishmaniasis poses challenges with unsatisfactory treatments, urging the development of novel therapeutic strategies. Effective prevention demands a One Health approach, integrating diverse disciplines to enhance diagnostics and develop safer drugs. Metalloenzymes, involved in parasite biology and critical in different biological pathways, emerged in the last few years as useful drug targets for the treatment of human diseases. Herein we have reviewed recent reports on the discovery of inhibitors of metalloenzymes associated with zoonotic diseases like histone deacetylases (HDACs), carbonic anhydrase (CA), arginase, and heme-dependent enzymes.

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.

Exploration of novel hydroxamate zinc binding group inhibitors against HDAC-1-3 enzymes by AI-based virtual screening: atomistic insights from steered molecular dynamics

Overexpression of histone deacetylase (HDAC) enzymes is linked to a wide variety of illnesses, including malignancies and neurological disorders, which makes HDAC inhibitors potentially therapeutic. However, most HDAC inhibitors lack subclass or isoform selectivity, which can be dangerous. Featuring both enhanced selectivity and toxicity profiles, slow-binding HDAC inhibitors offer promising treatment options for a variety of disorders. Diseases like cardiac, neurodegenerative disorders and diabetes are mainly associated with the HDAC1, HDAC2 and HDAC3 enzymes. The AI-based virtual screening tool PyRMD is implemented to identify the potential inhibitors from ∼2 million compounds. Based on the IC50 values, the top 10 compounds were selected for molecular docking. From the docking and ADMET study, the top-ranked three compounds were selected for molecular dynamics (MD) simulations. Further, to get more insights into the binding/unbinding mechanism of the ligand, we have employed the steered molecular dynamics (SMD) simulations. This study assists in developing Amber force field parameters for the HDAC1, HDAC2 and HDAC3 proteins and sheds light on the discovery of a potent drug. Our study suggests that hydroxamic acid derivative (i.e. referred to as Comp-1, CHEMBL600072) is the potential inhibitor for the series of HDAC-related diseases.Communicated by Ramaswamy H. Sarma.

Discussion on structure classification and regulation function of histone deacetylase and their inhibitor

Epigenetic regulation of genes through posttranslational regulation of proteins is a well-explored approach for disease treatment, particularly in cancer chemotherapy. Histone deacetylases have shown significant potential as effective drug targets in therapeutic studies aiming to restore epigenetic normality in oncology. Besides their role in modifying histones, histone deacetylases can also catalyze the deacetylation of various nonhistone proteins and participate in the regulation of multiple biological processes. This paper provides a review of the classification, structure, and functional characteristics of the four classes of human histone deacetylases. The increasing abundance of structural information on HDACs has led to the gradual elucidation of structural differences among subgroups and subtypes. This has provided a reasonable explanation for the selectivity of certain HDAC inhibitors. Currently, the US FDA has approved a total of six HDAC inhibitors for marketing, primarily for the treatment of various hematological tumors and a few solid tumors. These inhibitors all have a common pharmacodynamic moiety consisting of three parts: CAP, ZBG, and Linker. In this paper, the structure-effect relationship of HDAC inhibitors is explored by classifying the six HDAC inhibitors into three main groups: isohydroxamic acids, benzamides, and cyclic peptides, based on the type of inhibitor ZBG. However, there are still many questions that need to be answered in this field. In this paper, the structure-functional characteristics of HDACs and the structural information of the pharmacophore model and enzyme active region of HDAC is are considered, which can help to understand the inhibition mechanism of the compounds as well as the rational design of HDACs. This paper integrates the structural-functional characteristics of HDACs as well as the pharmacophore model of HDAC is and the structural information of the enzymatic active region, which not only contributes to the understanding of the inhibition mechanism of the compounds, but also provides a basis for the rational design 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.

A Nuclear-Targeted AIE Photosensitizer for Enzyme Inhibition and Photosensitization in Cancer Cell Ablation

The nucleus is considered the ideal target for anti-tumor therapy because DNA and some enzymes in the nucleus are the main causes of cell canceration and malignant proliferation. However, nuclear target drugs with good biosafety and high efficiency in cancer treatment are rare. Herein, a nuclear-targeted material MeTPAE with aggregation-induced emission (AIE) characteristics was developed based on a triphenylamine structure skeleton. MeTPAE can not only interact with histone deacetylases (HDACs) to inhibit cell proliferation but also damage telomere and nucleic acids precisely through photodynamic treatment (PDT). The cocktail strategy of MeTPAE caused obvious cell cycle arrest and showed excellent PDT anti-tumor activity, which offered new opportunities for the effective treatment of malignant tumors.

Histone deacetylase in neuropathology

Neuroepigenetics, a new branch of epigenetics, plays an important role in the regulation of gene expression. Neuroepigenetics is associated with holistic neuronal function and helps in formation and maintenance of memory and learning processes. This includes neurodevelopment and neurodegenerative defects in which histone modification enzymes appear to play a crucial role. These modifications, carried out by acetyltransferases and deacetylases, regulate biologic and cellular processes such as apoptosis and autophagy, inflammatory response, mitochondrial dysfunction, cell-cycle progression and oxidative stress. Alterations in acetylation status of histone as well as non-histone substrates lead to transcriptional deregulation. Histone deacetylase decreases acetylation status and causes transcriptional repression of regulatory genes involved in neural plasticity, synaptogenesis, synaptic and neural plasticity, cognition and memory, and neural differentiation. Transcriptional deactivation in the brain results in development of neurodevelopmental and neurodegenerative disorders. Mounting evidence implicates histone deacetylase inhibitors as potential therapeutic targets to combat neurologic disorders. Recent studies have targeted naturally-occurring biomolecules and micro-RNAs to improve cognitive defects and memory. Multi-target drug ligands targeting HDAC have been developed and used in cell-culture and animal-models of neurologic disorders to ameliorate synaptic and cognitive dysfunction. Herein, we focus on the implications of histone deacetylase enzymes in neuropathology, their regulation of brain function and plausible involvement in the pathogenesis of neurologic defects.

Mechanism of zinc ejection by disulfiram in nonstructural protein 5A

Hepatitis C virus (HCV) is a notorious member of the Flaviviridae family of enveloped, positive-strand RNA viruses. Non-structural protein 5A (NS5A) plays a key role in HCV replication and assembly. NS5A is a multi-domain protein which includes an N-terminal amphipathic membrane anchoring alpha helix, a highly structured domain-1, and two intrinsically disordered domains 2-3. The highly structured domain-1 contains a zinc finger (Zf)-site, and binding of zinc stabilizes the overall structure, while ejection of this zinc from the Zf-site destabilizes the overall structure. Therefore, NS5A is an attractive target for anti-HCV therapy by disulfiram, through ejection of zinc from the Zf-site. However, the zinc ejection mechanism is poorly understood. To disclose this mechanism based on three different states, A-state (NS5A protein), B-state (NS5A + Zn), and C-state (NS5A + Zn + disulfiram), we have performed molecular dynamics (MD) simulation in tandem with DFT calculations in the current study. The MD results indicate that disulfiram triggers Zn ejection from the Zf-site predominantly through altering the overall conformation ensemble. On the other hand, the DFT assessment demonstrates that the Zn adopts a tetrahedral configuration at the Zf-site with four Cys residues, which indicates a stable protein structure morphology. Disulfiram binding induces major conformational changes at the Zf-site, introduces new interactions of Cys39 with disulfiram, and further weakens the interaction of this residue with Zn, causing ejection of zinc from the Zf-site. The proposed mechanism elucidates the therapeutic potential of disulfiram and offers theoretical guidance for the advancement of drug candidates.

Harnessing the Role of HDAC6 in Idiopathic Pulmonary Fibrosis: Design, Synthesis, Structural Analysis, and Biological Evaluation of Potent Inhibitors

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease characterized by a progressive-fibrosing phenotype. IPF has been associated with aberrant HDAC activities confirmed by our immunohistochemistry studies on HDAC6 overexpression in IPF lung tissues. We herein developed a series of novel hHDAC6 inhibitors, having low inhibitory potency over hHDAC1 and hHDAC8, as potential pharmacological tools for IPF treatment. Their inhibitory potency was combined with low in vitro and in vivo toxicity. Structural analysis of 6h and structure-activity relationship studies contributed to the optimization of the binding mode of the new molecules. The best-performing analogues were tested for their efficacy in inhibiting fibrotic sphere formation and cell viability, proving their capability in reverting the IPF phenotype. The efficacy of analogue 6h was also determined in a validated human lung model of TGF-β1-dependent fibrogenesis. The results highlighted in this manuscript may pave the way for the identification of first-in-class molecules for the treatment of IPF.

Histone deacetylase 10, a potential epigenetic target for therapy

Histone deacetylase (HDAC) 10, a class II family, has been implicated in various tumors and non-tumor diseases, which makes the discovery of biological functions and novel inhibitors a fundamental endeavor. In cancers, HDAC10 plays crucial roles in regulating various cellular processes through its epigenetic functions or targeting some decisive molecular or signaling pathways. It also has potential clinical utility for targeting tumors and non-tumor diseases, such as renal cell carcinoma, prostate cancer, immunoglobulin A nephropathy (IgAN), intracerebral hemorrhage, human immunodeficiency virus (HIV) infection and schizophrenia. To date, relatively few studies have investigated HDAC10-specific inhibitors. Therefore, it is important to study the biological functions of HDAC10 for the future development of specific HDAC10 inhibitors. In this review, we analyzed the biological functions, mechanisms and inhibitors of HDAC10, which makes HDAC10 an appealing therapeutic target.

Fragment-Based Drug Design of Selective HDAC6 Inhibitors

Medicinal chemistry society has enough arguments to justify the usage of fragment-based drug design (FBDD) methodologies for the identification of lead compounds. Since the FDA approval of three kinase inhibitors - vemurafenib, venetoclax, and erdafitinib, FBDD has become a challenging alternative to high-throughput screening methods in drug discovery. The following protocol presents in silico drug design of selective histone deacetylase 6 (HDAC6) inhibitors through a fragment-based approach. To date, structural motifs that are important for HDAC inhibitory activity and selectivity are described as: surface recognition group (CAP group), aliphatic or aromatic linker, and zinc-binding group (ZBG). The main idea of this FBDD method is to identify novel and target-selective CAP groups by virtual scanning of publicly available fragment databases. Template structure used to search for novel heterocyclic and carbocyclic fragments is 1,8-naphthalimide (CAP group of scriptaid, a potent HDAC inhibitor). Herein, the design of HDAC6 inhibitors is based on linking the identified fragments with the aliphatic or aromatic linker and hydroxamic acid (ZBG) moiety. Final selection of potential selective HDAC6 inhibitors is based on combined structure-based (molecular docking) and ligand-based (three-dimensional quantitative structure-activity relationships, 3D-QSAR) techniques. Designed compounds are docked in the active site pockets of human HDAC1 and HDAC6 isoforms, and their docking conformations used to predict their HDAC inhibitory and selectivity profiles through two developed 3D-QSAR models (describing HDAC1 and HDAC6 inhibitory activities).

Novel 4-Oxoquinazoline-Based N-Hydroxypropenamides as Histone Deacetylase Inhibitors: Design, Synthesis, and Biological Evaluation

Two series of novel 4-oxoquinazoline-based N-hydroxypropenamides (9a-m and 10a-m) were designed, synthesized, and evaluated for their inhibitory and cytotoxicity activities against histone deacetylase (HDAC). The compounds showed good to potent HDAC inhibitory activity and cytotoxicity against three human cancer cell lines (SW620, colon; PC-3, prostate; NCI-H23, lung cancer). In this series, compounds with the N-hydroxypropenamide functionality impeded at position 7 on the 4-oxoquinazoline skeleton (10a-m) were generally more potent than compounds with the N-hydroxypropenamide moiety at position 6 (9a-m). Also, the N ³-benzyl-substituted derivatives (9h-m, 10h-m) exhibited stronger bioactivity than the N ³-alkyl-substituted ones (9a-e, 10a-e). Two compounds 10l and 10m were the most potent ones. Their HDAC inhibitory activity (IC₅₀ values, 0.041-0.044 μM) and cytotoxicity (IC₅₀ values, 0.671-1.211 μM) were approximately 2- to 3-fold more potent than suberoylanilide hydroxamic acid (SAHA). Some compounds showed up to 10-fold more potent HDAC6 inhibition compared to their inhibitory activity in total HDAC extract assay. Analysis of selected compounds 10l and 10m revealed that these compounds strongly induced both early and late apoptosis and arrested SW620 cells at the G2/M phase. Docking studies were carried out on the HDAC6 isoform for series 10a-m and revealed some important features contributing to the inhibitory activity of synthesized compounds.

Novel Conjugated Quinazolinone-Based Hydroxamic Acids: Design, Synthesis and Biological Evaluation

BACKGROUND

The target-based approach to drug discovery currently attracts a great deal of interest from medicinal chemists in anticancer drug discovery and development. Histone deacetylase (HDAC) inhibitors represent an extensive class of targeted anti-cancer agents. Among the most explored structure moieties, hydroxybenzamides and hydroxypropenamides have been demonstrated to have potential HDAC inhibitory effects. Several compounds of these structural classes have been approved for clinical uses to treat different types of cancer, such as vorinostat and belinostat.

AIMS

This study aims at developing novel HDAC inhibitors bearing conjugated quinazolinone scaffolds with potential cytotoxicity against different cancer cell lines.

METHODS

A series of novel N-hydroxyheptanamides incorporating conjugated 6-hydroxy-2 methylquinazolin- 4(3H)-ones (15a-l) was designed, synthesized and evaluated for HDAC inhibitory potency as well as cytotoxicity against three human cancer cell lines, including HepG-2, MCF-7 and SKLu-1. Molecular simulations were finally performed to gain more insight into the structureactivity relationships.

RESULTS

It was found that among novel conjugated quinazolinone-based hydroxamic acids synthesized, compounds 15a, 15c and 15f were the most potent, both in terms of HDAC inhibition and cytotoxicity. Especially, compound 15f displayed up to nearly 4-fold more potent than SAHA (vorinostat) in terms of cytotoxicity against MCF-7 cell line with IC50 value of 1.86 μM, and HDAC inhibition with IC50 value of 6.36 μM. Docking experiments on HDAC2 isozyme showed that these compounds bound to HDAC2 with binding affinities ranging from -10.08 to -14.93 kcal/mol compared to SAHA (-15.84 kcal/mol). It was also found in this research that most of the target compounds seemed to be more cytotoxic toward SKLu-1than MCF-7 and HepG-2.

CONCLUSION

The resesrch results suggest that some hydroxamic acids could emerge for further evaluation and the results are well served as basics for further design of more potent HDAC inhibitors and antitumor agents.

Molecular Docking Simulations on Histone Deacetylases (HDAC)-1 and -2 to Investigate the Flavone Binding

Histone modifications through acetylation are fundamental for remodelling chromatin and consequently activating gene expression. The imbalance between acetylation and deacetylation activity causes transcriptional dysregulation associated with several disorders. Flavones, small molecules of plant origin, are known to interfere with class I histone deacetylase (HDAC) enzymes and to enhance acetylation, restoring cell homeostasis. To investigate the possible physical interactions of flavones on human HDAC1 and 2, we carried out in silico molecular docking simulations. Our data have revealed how flavone, and other two flavones previously investigated, i.e., apigenin and luteolin, can interact as ligands with HDAC1 and 2 at the active site binding pocket. Regulation of HDAC activity by dietary flavones could have important implications in developing epigenetic therapy to regulate the cell gene expression.

OleD Loki as a Catalyst for Hydroxamate Glycosylation

Herein we describe the ability of the permissive glycosyltransferase (GT) OleD Loki to convert a diverse set of >15 histone deacetylase (HDAC) inhibitors (HDACis) into their corresponding hydroxamate glycosyl esters. Representative glycosyl esters were subsequently evaluated in assays for cancer cell line cytotoxicity, chemical and enzymatic stability, and axolotl embryo tail regeneration. Computational substrate docking models were predictive of enzyme-catalyzed turnover and suggest certain HDACis may form unproductive, potentially inhibitory, complexes with GTs.

A short guide to histone deacetylases including recent progress on class II enzymes

The interaction between histones and DNA is important for eukaryotic gene expression. A loose interaction caused, for example, by the neutralization of a positive charge on the histone surface by acetylation, induces a less compact chromatin structure, resulting in feasible accessibility of RNA polymerase and increased gene expression. In contrast, the formation of a tight chromatin structure due to the deacetylation of histone lysine residues on the surface by histone deacetylases enforces the interaction between the histones and DNA, which minimizes the chance of RNA polymerases contacting DNA, resulting in decreased gene expression. Therefore, the balance of the acetylation of histones mediated by histone acetylases (HATs) and histone deacetylases (HDACs) is an issue of transcription that has long been studied in relation to posttranslational modification. In this review, current knowledge of HDACs is briefly described with an emphasis on recent progress in research on HDACs, especially on class IIa HDACs.

Targeting specific structural and functional features of enzymes involved in regulating the interactions between DNA and the histone proteins associated with it could lead to the development of more effective cancer therapeutics. Histone deacetylases (HDACs), enzymes which remove acetyl groups from histones, make the histones wrap more tightly around the DNA so that it becomes inaccessible to the initial steps in gene expression. Drugs that target these enzymes have shown limited efficacy due to lack of specificity and off-target toxicity. Jeong-Sun Kim at Chonnam National University, Gwangju, and Suk-Youl Park at Pohang Accelerator Laboratory, Pohang University of Science and Technology, South Korea, review the latest knowledge about class II HDACs. They suggest that their unique structural features and low enzymatic activity are important features to consider when designing new, more selective HDAC inhibitors.

Catalytic Roles of Coenzyme Pyridoxal-5'-phosphate (PLP) in PLP-dependent Enzymes: Reaction Pathway for Methionine-γ-lyase-catalyzed L-methionine Depletion

Pyridoxal-5'-phosphate (PLP), the active form of vitamin B₆, is an important and versatile coenzyme involved in a variety of enzymatic reactions, accounting for about 4% of all classified activities. However, the detailed catalytic reaction pathways for PLP-dependent enzymes remain to be explored. Methionine-γ-lyase (MGL), a promising alternative anti-tumor agent to conventional chemotherapies whose catalytic mechanism is highly desired for guiding further development of re-engineered enzymes, was used as a representative PLP-dependent enzyme, and the catalytic mechanism for L-Met elimination by MGL was explored at the first-principles quantum mechanical/molecular mechanical (QM/MM) level with umbrella sampling. The QM/MM calculations revealed that the enzymatic reaction pathway consists of 4 stages for a total of 19 reaction steps with five intermediates captured in available crystal structures. Furthermore, the more comprehensive role of PLP was revealed. Besides the commonly known role of "electron sink", coenzyme PLP can also assist proton transfer and temporarily store the excess proton generated in some intermediate states by using its hydroxyl group and phosphate group. Thus, PLP is participated in most of the 19 steps. This study not only provided a theoretical basis for further development and re-engineering MGL as a potential anti-tumor agent, but also revealed the comprehensive role of PLP which could be used to explore the mechanisms of other PLP-dependent enzymes.

Conformational dynamics and allosteric effect modulated by the unique zinc-binding motif in class IIa HDACs

Class IIa histone deacetylases (HDACs) have been considered as potential targets for the treatment of several diseases. Compared to other HDACs, class IIa HDACs have an additional second zinc binding motif. So far, the function of the unique zinc-binding motif is still not very clear. In this work, extensive classical molecular dynamics (MD) simulations were employed to illuminate the conformational change modulated by the unique zinc-binding motif. It has been revealed that the unique zinc-binding motif is a crucial structural stabilization factor in retaining the catalytic activity of the enzyme and the stability of the multi-protein complex, by remotely modulating the active site pocket in a "closed" conformation. Moreover, it is also revealed that the Loop2 motion in HDAC4 is less flexible than that in HDAC7, which opens a new avenue to design selective inhibitors by utilizing the local conformational dynamics difference between the structurally highly similar HDAC4 and HDAC7. Finally, by comparative studies with class I HDACs (HDAC1-3), it is found that the reversible "in-out" conformational transformation of the binding rail (highly conserved both in class I and IIa HDACs) occurs spontaneously in HDAC1-3, whereas the binding rail is trapped in an "in" conformation owing to the strong metal coordination interaction of the unique CCHC zinc-binding motif in class IIa HDACs. Thus, the CCHC zinc-binding motif may be a feasible allosteric site for the development of class IIa-selective inhibitors.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Quinazolin-4(3H)-one-Based Hydroxamic Acids: Design, Synthesis and Evaluation of Histone Deacetylase Inhibitory Effects and Cytotoxicity

The present article describes the synthesis and biological activity of various series of novel hydroxamic acids incorporating quinazolin-4(3H)-ones as novel small molecules targeting histone deacetylases. Biological evaluation showed that these hydroxamic acids were potently cytotoxic against three human cancer cell lines (SW620, colon; PC-3, prostate; NCI-H23, lung). Most compounds displayed superior cytotoxicity than SAHA (suberoylanilide hydroxamic acid, Vorinostat) in term of cytotoxicity. Especially, N-hydroxy-7-(7-methyl-4-oxoquinazolin-3(4H)-yl)heptanamide (5b) and N-hydroxy-7-(6-methyl-4-oxoquinazolin-3(4H)-yl)heptanamide (5c) (IC₅₀ values, 0.10-0.16 μm) were found to be approximately 30-fold more cytotoxic than SAHA (IC₅₀ values of 3.29-3.67 μm). N-Hydroxy-7-(4-oxoquinazolin-3(4H)-yl)heptanamide (5a; IC₅₀ values of 0.21-0.38 μm) was approximately 10- to 15-fold more potent than SAHA in cytotoxicity assay. These compounds also showed comparable HDAC inhibition potency with IC₅₀ values in sub-micromolar ranges. Molecular docking experiments indicated that most compounds, as represented by 5b and 5c, strictly bound to HDAC2 at the active binding site with binding affinities much higher than that of SAHA.

Purification and enzymatic assay of class I histone deacetylase enzymes

The reversible acetylation of histones has a profound influence on transcriptional status. Histone acetyltransferases catalyze the addition of these chemical modifications to histone lysine residues. Conversely, histone deacetylases (HDACs) catalyze the removal of these acetyl groups from histone lysine residues. As modulators of transcription, HDACs have found themselves as targets of several FDA-approved chemotherapeutic compounds which aim to inhibit enzyme activity. The ongoing efforts to develop targeted and isoform-specific HDAC inhibitors necessitates tools to study these modifications and the enzymes that maintain an equilibrium of these modifications. In this chapter, we present an optimized workflow for the isolation of recombinant protein and subsequent assay of class I HDAC activity. We demonstrate the application of this assay by assessing the activities of recombinant HDAC1, HDAC2, and SIN3B. This assay system utilizes readily available reagents and can be used to assess the activity and responsiveness of class I HDAC complexes to HDAC inhibitors.

Perfluorinated HDAC inhibitors as selective anticancer agents

A series of potent histone deacetylase inhibitors is presented that incorporate alkyl or perfluorinated alkyl chains. Several new compounds show greater in vitro antiproliferative activity than the clinically approved inhibitor, SAHA. Furthermore, the new compounds show up to 5-fold greater activity against cancer cells than healthy cells. This selectivity is in contrast to SAHA, which is more active against the healthy cell line than the cancer cell line tested. Finally, we report an increase in activity for SAHA under mild hyperthermia, indicating that it could be an interesting candidate to use in combination with thermal therapy.

Computational insights for the hydride transfer and distinctive roles of key residues in cholesterol oxidase

Cholesterol oxidase (ChOx), a member of the glucose-methanol-choline (GMC) family, catalyzes the oxidation of the substrate via a hydride transfer mechanism and concomitant reduction of the FAD cofactor. Unlike other GMC enzymes, the conserved His447 is not the catalytic base that deprotonates the substrate in ChOx. Our QM/MM MD simulations indicate that the Glu361 residue acts as a catalytic base facilitating the hydride transfer from the substrate to the cofactor. We find that two rationally chosen point mutations (His447Gln and His447Asn) cause notable decreases in the catalytic activity. The binding free energy calculations show that the Glu361 and His447 residues are important in substrate binding. We also performed high-level double-hybrid density functional theory simulations using small model systems, which support the QM/MM MD results. Our work provides a basis for unraveling the substrate oxidation mechanism in GMC enzymes in which the conserved histidine does not act as a base.

Intrinsic Dynamics of the Binding Rail and Its Allosteric Effect in the Class I Histone Deacetylases

The development of novel isoform/class-selective inhibitors is still of great biological and medical significance to conquer the continuously reported side effects for the histone deacetylase (HDAC) drugs. The first potent HDAC allosteric inhibitor was discovered last year, and this allosteric inhibitor design is thought to be a promising strategy to overcome the current challenges in HDAC inhibitor design. However, the detailed allosteric mechanism and its remote regulatory effects on the catalytic/inhibitor activity of HDAC are still unclear. In this work, on the basis of microsecond-time-scale all-atom molecular dynamics (MD) simulations and picosecond-time-scale density functional theory/molecular mechanics MD simulations on HDAC8, we propose that the allostery is achieved by the intrinsic conformational flexibility of the binding rail (constituted by a highly conserved X-D residue dyad), which steers the loop-loop motion and creates the diverse shapes of the allosteric sites in different HDAC isoforms. Additionally, the rotatability of the binding rail is an inherent structural feature that regulates the hydrophobicity of the linker binding channel and thus further affects the HDAC enzyme inhibitory/catalytic activity by utilizing the promiscuity of X-D dyad. Since the plastic X residue is different among class I HDACs, these new findings provide a deeper understanding of the allostery, which is guidable for the design of new allosteric inhibitors toward the allosteric site and structure modifications on the conventional inhibitors binding into the active pocket by exploiting the intrinsic dynamic features of the conserved X-D dyad.

Mechanisms for the deamination reaction of 8-oxoguanine catalyzed by 8-oxoguanine deaminase: A combined QM/MM molecular dynamics study

The deamination reaction of 8-oxoguanine (8-oxoG) catalyzed by 8-oxoguanine deaminase (8-oxoGD) plays a critically important role in the DNA repair activity for oxidative damage. In order to elucidate the complete enzymatic catalysis mechanism at the stages of 8-oxoguanine binding, departure of 2-hydroxy-1H-purine-6,8(7H,9H)-dione from the active site, and formation of 8-oxoxanthine, extensive combined QM(PM3)/MM molecular dynamics simulations have been performed. Computations show that the rate-limiting step corresponds to the nucleophilic attack from zinc-coordinate hydroxide group to free 8-oxoguanine. Through conformational analyses, we demonstrate that Trp115, Trp123 and Leu119 connect to O8@8-oxoguanine with hydrogen bonds, and we suggest that mutations of tryptophan (115 and 123) to histidine or phenylalanine and mutation of leucine (119) to alanine could potentially lead to a mutant with enhanced activity. On this ground, a proton transfer mechanism for the formation of 8-oxoxanthine was further discussed. Both Glu218 and water molecule could be used as proton shuttles, and water molecule plays a major role in proton transfer in substrate. On the other hand, comparative simulations on the deamination of guanine and isocytosine reveal that, for the helping of hydrogen bonds between O8@8-oxoguanine and enzyme, O8@8-oxoguanine is the fastest to be deaminated among the three substrates which are also supported by the experimental kinetic constants.

Targeting prostate cancer with compounds possessing dual activity as androgen receptor antagonists and HDAC6 inhibitors

While enzalutamide and abiraterone are approved for treatment of metastatic castration-resistant prostate cancer (mCRPC), approximately 20-40% of patients have no response to these agents. It has been stipulated that the lack of response and the development of secondary resistance to these drugs may be due to the presence of AR splice variants. HDAC6 has a role in regulating the androgen receptor (AR) by modulating heat shock protein 90 (Hsp90) acetylation, which controls the nuclear localization and activation of the AR in androgen-dependent and independent scenarios. With dual-acting AR-HDAC6 inhibitors it should be possible to target patients who don't respond to enzalutamide. Herein, we describe the design, synthesis and biological evaluation of dual-acting compounds which target AR and are also specific towards HDAC6. Our efforts led to compound 10 which was found to have potent dual activity (HDAC6 IC₅₀=0.0356μM and AR binding IC₅₀=<0.03μM). Compound 10 was further evaluated for antagonist and other cell-based activities, in vitro stability and pharmacokinetics.

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.

Inside HDACs with more selective HDAC inhibitors

Inhibitors of histone deacetylases (HDACs) are nowadays part of the therapeutic arsenal mainly against cancers, with four compounds approved by the Food and Drug Administration. During the last five years, several groups have made continuous efforts to improve this class of compounds, designing more selective compounds or compounds with multiple capacities. After a survey of the HDAC biology and structures, this review summarizes the results of the chemists working in this field, and highlights when possible the behavior of the molecules inside their targets.

An Automated Strategy for Binding-Pose Selection and Docking Assessment in Structure-Based Drug Design

Molecular docking is a widely used technique in drug design to predict the binding pose of a candidate compound in a defined therapeutic target. Numerous docking protocols are available, each characterized by different search methods and scoring functions, thus providing variable predictive capability on a same ligand-protein system. To validate a docking protocol, it is necessary to determine a priori the ability to reproduce the experimental binding pose (i.e., by determining the docking accuracy (DA)) in order to select the most appropriate docking procedure and thus estimate the rate of success in docking novel compounds. As common docking programs use generally different root-mean-square deviation (RMSD) formulas, scoring functions, and format results, it is both difficult and time-consuming to consistently determine and compare their predictive capabilities in order to identify the best protocol to use for the target of interest and to extrapolate the binding poses (i.e., best-docked (BD), best-cluster (BC), and best-fit (BF) poses) when applying a given docking program over thousands/millions of molecules during virtual screening. To reduce this difficulty, two new procedures called Clusterizer and DockAccessor have been developed and implemented for use with some common and "free-for-academics" programs such as AutoDock4, AutoDock4(Zn), AutoDock Vina, DOCK, MpSDockZn, PLANTS, and Surflex-Dock to automatically extrapolate BD, BC, and BF poses as well as to perform consistent cluster and DA analyses. Clusterizer and DockAccessor (code available over the Internet) represent two novel tools to collect computationally determined poses and detect the most predictive docking approach. Herein an application to human lysine deacetylase (hKDAC) inhibitors is illustrated.

A salt bridge turns off the foot-pocket in class-II HDACs

It is for the first time revealed that a conserved R–E salt bridge turns off the foot-pocket in class-II HDACs.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Inhibition mechanism of SAHA in HDAC: a revisit

SAHA (vorinostat, Merck) is a famous clinical drug for zinc-containing histone deacetylase (HDAC) targets against cancer and several other human disorders, whose inhibition mechanism (namely the protonation mechanism) upon binding to HDAC has been debated for more than ten years. It is very challenging to verify experimentally and is still controversial theoretically. The popular "Class-dependent" (namely "Tyr-dependent") hypothesis is that the deprotonation of SAHA is mostly regulated by the conserved Tyr308 in class I HDAC while it is replaced by the His843 in class IIa HDAC. Herein, by elaborate QM(DFT)/MM MD simulations, we exclude the prevalent "Class-dependent" mechanism and advance a novel "Metal-dependent" mechanism, where the remote second metal site (K(+) in most HDAC and Ca(2+) in HDAC2) determines the protonation of SAHA. This proof-of-principle "Metal-dependent" mechanism opens up a new avenue to utilize the second metal site for isoform-selective inhibitor design.

Thiol versus hydroxamate as zinc binding group in HDAC inhibition: An ab initio QM/MM molecular dynamics study

Zinc-dependent histone deacetylases (HDACs) play a critical role in transcriptional repression and gene silencing, and are among the most attractive targets for the development of new therapeutics against cancer and various other diseases. Two HDAC inhibitors have been approved by FDA as anti-cancer drugs: one is SAHA whose hydroxamate is directly bound to zinc, the other is FK228 whose active form may use thiol as the zinc binding group. In spite of extensive studies, it remains to be ambiguous regarding how thiol and hydroxamate are bound to the zinc active site of HDACs. In this work, our computational approaches center on Born-Oppenheimer ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics with umbrella sampling, which allow for modeling of the zinc active site with reasonable accuracy while properly including dynamics and effects of protein environment. Meanwhile, an improved short-long effective function (SLEF2) to describe non-bonded interactions between zinc and other atoms has been employed in initial MM equilibrations. Our ab initio QM/MM MD simulations have confirmed that hydroxamate is neutral when it is bound to HDAC8, and found that thiol is deprotonated when directly bound to zinc in the HDAC active site. By comparing thiol and hydroxamate, our results elucidated the differences in their binding environment in the HDAC active sites, and emphasized the importance of the linker design to achieve more specific binding toward class IIa HDACs.

Computational design of a time-dependent histone deacetylase 2 selective inhibitor

Development of isoform-selective histone deacetylase (HDAC) inhibitors is of great biological and medical interest. Among 11 zinc-dependent HDAC isoforms, it is particularly challenging to achieve isoform inhibition selectivity between HDAC1 and HDAC2 due to their very high structural similarities. In this work, by developing and applying a novel de novo reaction-mechanism-based inhibitor design strategy to exploit the reactivity difference, we have discovered the first HDAC2-selective inhibitor, β-hydroxymethyl chalcone. Our bioassay experiments show that this new compound has a unique time-dependent selective inhibition on HDAC2, leading to about 20-fold isoform-selectivity against HDAC1. Furthermore, our ab initio QM/MM molecular dynamics simulations, a state-of-the-art approach to study reactions in biological systems, have elucidated how the β-hydroxymethyl chalcone can achieve the distinct time-dependent inhibition toward HDAC2.