Structure-Activity Studies of 1,2,4-Oxadiazoles for the Inhibition of the NAD+-Dependent Lysine Deacylase Sirtuin 2.

The NAD+-dependent lysine deacylase sirtuin 2 (Sirt2) is involved in multiple pathological conditions such as cancer. Targeting Sirt2 has thus received an increased interest for therapeutic purposes. Furthermore, the orthologue from Schistosoma mansoni (SmSirt2) has been considered for the potential treatment of the neglected tropical disease schistosomiasis. We previously identified a 1,2,4-oxadiazole-based scaffold from the screening of the "Kinetobox" library as a dual inhibitor of human Sirt2 (hSirt2) and SmSirt2. Herein, we describe the structure-activity studies on 1,2,4-oxadiazole-based analogues, which are potent inhibitors of human Sirt2 deacetylation. As proposed by docking studies, a substrate-competitive and cofactor-noncompetitive binding mode of inhibition could be determined in vitro via binding assays and kinetic analysis and further confirmed by a crystal structure of an oxadiazole inhibitor in complex with hSirt2. Optimized analogues reduced cell viability and inhibited prostate cancer cell migration, in correlation with Sirt2 deacetylase inhibition both in vitro and in cells.

Direct GPCR-EGFR interaction enables synergistic membrane-to-nucleus information transfer.

We addressed the heteromerization of the epidermal growth factor receptor (EGFR) with G-protein coupled receptors (GPCR) on the basis of angiotensin-II-receptor-subtype-1(AT1R)-EGFR interaction as proof-of-concept and show its functional relevance during synergistic nuclear information transfer, beyond ligand-dependent EGFR transactivation. Following in silico modelling, we generated EGFR-interaction deficient AT1R-mutants and compared them to AT1R-wildtype. Receptor interaction was assessed by co-immunoprecipitation (CoIP), Förster resonance energy transfer (FRET) and fluorescence-lifetime imaging microscopy (FLIM). Changes in cell morphology, ERK1/2-phosphorylation (ppERK1/2), serum response factor (SRF)-activation and cFOS protein expression were determined by digital high content microscopy at the single cell level. FRET, FLIM and CoIP confirmed the physical interaction of AT1R-wildtype with EGFR that was strongly reduced for the AT1R-mutants. Responsiveness of cells transfected with AT1R-WT or -mutants to angiotensin II or EGF was similar regarding changes in cell circularity, ppERK1/2 (direct and by ligand-dependent EGFR-transactivation), cFOS-expression and SRF-activity. By contrast, the EGFR-AT1R-synergism regarding these parameters was completely absent for in the interaction-deficient AT1R mutants. The results show that AT1R-EGFR heteromerisation enables AT1R-EGFR-synergism on downstream gene expression regulation, modulating the intensity and the temporal pattern of nuclear AT1R/EGFR-information transfer. Furthermore, remote EGFR transactivation, via ligand release or cytosolic tyrosine kinases, is not sufficient for the complete synergistic control of gene expression.

Fatty acid binding to the human transport proteins FABP3, FABP4, and FABP5 from a Ligand's perspective.

Fatty acid binding proteins (FABPs) are a family of amphiphilic transport proteins with high diversity in terms of their amino acid sequences and binding preferences. Beyond their main biological role as cytosolic fatty acid transporters, many aspects regarding their binding mechanism and functional specializations in human cells remain unclear. In this work, the binding properties and thermodynamics of FABP3, FABP4, and FABP5 were analyzed under various physical conditions. For this purpose, the FABPs were loaded with fatty acids bearing fluorescence or spin probes as model ligands, comparing their binding affinities via microscale thermophoresis (MST) and continuous-wave electron paramagnetic resonance (CW EPR) spectroscopy. The CW EPR spectra of non-covalently bound 5- and 16-DOXYL stearic acid (5/16-DSA) deliver in-depth information about the dynamics and chemical environments of ligands inside the binding pockets of the FABPs. EPR spectral simulations allow the construction of binding curves, revealing two different binding states ('intermediately' and 'strongly' bound). The proportion of bound 5/16-DSA depends strongly on the FABP concentration and the temperature but with remarkable differences between the three isoforms. Additionally, the more dynamic state ('intermediately bound') seems to dominate at body temperature with thermodynamic preference. The ligand binding studies were supplemented by aggregation studies via dynamic light scattering and bioinformatic analyses. Beyond the remarkably fine-tuned binding properties exhibited by each FABP, which were discernible with our EPR-centered approach, the results of this work attest to the power of simple spectroscopic experiments to provide new insights into the ligand binding mechanisms of proteins in general on a molecular level.

Exploring Aromatic Cage Flexibility Using Cosolvent Molecular Dynamics Simulations─An In-Silico Case Study of Tudor Domains.

Cosolvent molecular dynamics (MD) simulations have proven to be powerful in silico tools to predict hotspots for binding regions on protein surfaces. In the current study, the method was adapted and applied to two Tudor domain-containing proteins, namely Spindlin1 (SPIN1) and survival motor neuron protein (SMN). Tudor domains are characterized by so-called aromatic cages that recognize methylated lysine residues of protein targets. In the study, the conformational transitions from closed to open aromatic cage conformations were investigated by performing MD simulations with cosolvents using six different probe molecules. It is shown that a trajectory clustering approach in combination with volume and atomic distance tracking allows a reasonable discrimination between open and closed aromatic cage conformations and the docking of inhibitors yields very good reproducibility with crystal structures. Cosolvent MDs are suitable to capture the flexibility of aromatic cages and thus represent a promising tool for the optimization of inhibitors.

Comparative Structure-Based Virtual Screening Utilizing Optimized AlphaFold Model Identifies Selective HDAC11 Inhibitor.

HDAC11 is a class IV histone deacylase with no crystal structure reported so far. The catalytic domain of HDAC11 shares low sequence identity with other HDAC isoforms, which makes conventional homology modeling less reliable. AlphaFold is a machine learning approach that can predict the 3D structure of proteins with high accuracy even in absence of similar structures. However, the fact that AlphaFold models are predicted in the absence of small molecules and ions/cofactors complicates their utilization for drug design. Previously, we optimized an HDAC11 AlphaFold model by adding the catalytic zinc ion and minimization in the presence of reported HDAC11 inhibitors. In the current study, we implement a comparative structure-based virtual screening approach utilizing the previously optimized HDAC11 AlphaFold model to identify novel and selective HDAC11 inhibitors. The stepwise virtual screening approach was successful in identifying a hit that was subsequently tested using an in vitro enzymatic assay. The hit compound showed an IC50 value of 3.5 µM for HDAC11 and could selectively inhibit HDAC11 over other HDAC subtypes at 10 µM concentration. In addition, we carried out molecular dynamics simulations to further confirm the binding hypothesis obtained by the docking study. These results reinforce the previously presented AlphaFold optimization approach and confirm the applicability of AlphaFold models in the search for novel inhibitors for drug discovery.

Specific Inhibitors of Mitochondrial Deacylase Sirtuin 4 Endowed with Cellular Activity.

Sirtuins are NAD+-dependent protein lysine deacylases implicated in aging-related diseases. Mammalian Sirtuin 4 (Sirt4) is located in mitochondria and a potential therapeutic target for cancer and metabolic diseases, but no potent and selective Sirt4 inhibitors have been reported. Here, we describe the identification of potent Sirt4-specific small-molecule inhibitors. Testing hits from a target-based virtual screen revealed 12 active compounds. A focused screen based on two top compounds, followed by structure-assisted design of derivatives, yielded four first-in-class potent Sirt4 inhibitors. Kinetic analyses indicate compound competition with the acyl peptide substrate, consistent with the docking models and implicating Sirt4's unique acyl binding site. The compounds indeed show preference for Sirt4 over other isoforms, with one of them (69) being highly isoform selective, and they are active in cells. Our results provide first lead compounds and mechanistic insights for optimization toward Sirt4-specific inhibitors useful as experimental tools and potential therapeutics.

Utilization of AlphaFold models for drug discovery: Feasibility and challenges. Histone deacetylase 11 as a case study.

Histone deacetylase 11 (HDAC11), an enzyme that cleaves acyl groups from acylated lysine residues, is the sole member of class IV of HDAC family with no reported crystal structure so far. The catalytic domain of HDAC11 shares low sequence identity with other HDAC isoforms which complicates the conventional template-based homology modeling. AlphaFold is a neural network machine learning approach for predicting the 3D structures of proteins with atomic accuracy even in absence of similar structures. However, the structures predicted by AlphaFold are missing small molecules as ligands and cofactors. In our study, we first optimized the HDAC11 AlphaFold model by adding the catalytic zinc ion followed by assessment of the usability of the model by docking of the selective inhibitor FT895. Minimization of the optimized model in presence of transplanted inhibitors, which have been described as HDAC11 inhibitors, was performed. Four complexes were generated and proved to be stable using three replicas of 50 ns MD simulations and were successfully utilized for docking of the selective inhibitors FT895, MIR002 and SIS17. For SIS17, The most reasonable pose was selected based on structural comparison between HDAC6, HDAC8 and the HDAC11 optimized AlphaFold model. The manually optimized HDAC11 model is thus able to explain the binding behavior of known HDAC11 inhibitors and can be used for further structure-based optimization.

Application of Ligand- and Structure-Based Prediction Models for the Design of Alkylhydrazide-Based HDAC3 Inhibitors as Novel Anti-Cancer Compounds.

Histone deacetylases (HDAC) represent promising epigenetic targets for several diseases including different cancer types. The HDAC inhibitors approved to date are pan-HDAC inhibitors and most show a poor selectivity profile, side effects, and in particular hydroxamic-acid-based inhibitors lack good pharmacokinetic profiles. Therefore, the development of isoform-selective non-hydroxamic acid HDAC inhibitors is a highly regarded field in medicinal chemistry. In this study, we analyzed different ligand-based and structure-based drug design techniques to predict the binding mode and inhibitory activity of recently developed alkylhydrazide HDAC inhibitors. Alkylhydrazides have recently attracted more attention as they have shown promising effects in various cancer cell lines. In this work, pharmacophore models and atom-based quantitative structure-activity relationship (QSAR) models were generated and evaluated. The binding mode of the studied compounds was determined using molecular docking as well as molecular dynamics simulations and compared with known crystal structures. Calculated free energies of binding were also considered to generate QSAR models. The created models show a good explanation of in vitro data and were used to develop novel HDAC3 inhibitors.

Chemically Diverse S. mansoni HDAC8 Inhibitors Reduce Viability in Worm Larval and Adult Stages.

Schistosoma mansoni HDAC8 is a reliable target to fight schistosomiasis, and several inhibitors have been reported in the literature up to now. Nevertheless, only a few displayed selectivity over the human deacetylases and some exhibited very low or no activity against parasite larvae and/or adult worms. We report here the in vitro enzyme and biological activity of a small library of HDAC inhibitors from our lab, in many cases exhibiting submicromolar/nanomolar potency against smHDAC8 and diverse degrees of selectivity over hHDAC1 and/or hHDAC6. Such compounds were tested against schistosomula, and a selection of them against the adult forms of S. mansoni, to detect their effect on viability. Some of them showed the highest viability reduction for the larval stage with IC50 values around 1 µM and/or displayed ∼40-50 % activity in adult worms at 10 µM, joined to moderate to no toxicity in human fibroblast MRC-5 cells.

Development of Alkylated Hydrazides as Highly Potent and Selective Class I Histone Deacetylase Inhibitors with T cell Modulatory Properties.

Histone deacetylases (HDACs) are epigenetic regulators and additionally control the activity of non-histone substrates. We recently demonstrated that inhibition of HDAC8 overexpressed in various of cancers reduces hepatocellular carcinoma tumorigenicity in a T cell-dependent manner. Here, we present alkylated hydrazide-based class I HDAC inhibitors in which the n-hexyl side chain attached to the hydrazide moiety shows HDAC8 selectivity in vitro. Analysis of the mode of inhibition of the most promising compound 7d against HDAC8 revealed a substrate-competitive binding mode. 7d marked induced acetylation of the HDAC8 substrates H3K27 and SMC3 but not tubulin in CD4+ T lymphocytes, and significantly upregulated gene expressions for memory and effector functions. Furthermore, intraperitoneal injection of 7d (10 mg/kg) in C57BL/6 mice increased interleukin-2 expression in CD4+ T cells and CD8+ T cell proportion with no apparent toxicity. This study expands a novel chemotype of HDAC8 inhibitors with T cell modulatory properties for future therapeutic applications.

Design, Synthesis and Biological Characterization of Histone Deacetylase 8 (HDAC8) Proteolysis Targeting Chimeras (PROTACs) with Anti-Neuroblastoma Activity.

In addition to involvement in epigenetic gene regulation, histone deacetylases (HDACs) regulate multiple cellular processes through mediating the activity of non-histone protein substrates. The knockdown of HDAC8 isozyme is associated with the inhibition of cell proliferation and apoptosis enhancement in several cancer cell lines. As shown in several studies, HDAC8 can be considered a potential target in the treatment of cancer forms such as childhood neuroblastoma. The present work describes the development of proteolysis targeting chimeras (PROTACs) of HDAC8 based on substituted benzhydroxamic acids previously reported as potent and selective HDAC8 inhibitors. Within this study, we investigated the HDAC8-degrading profiles of the synthesized PROTACs and their effect on the proliferation of neuroblastoma cells. The combination of in vitro screening and cellular testing demonstrated selective HDAC8 PROTACs that show anti-neuroblastoma activity in cells.

First Fluorescent Acetylspermidine Deacetylation Assay for HDAC10 Identifies Selective Inhibitors with Cellular Target Engagement.

Histone deacetylases (HDACs) are important epigenetic regulators involved in many diseases, especially cancer. Five HDAC inhibitors have been approved for anticancer therapy and many are in clinical trials. Among the 11 zinc-dependent HDACs, HDAC10 has received relatively little attention by drug discovery campaigns, despite its involvement, e. g., in the pathogenesis of neuroblastoma. This is due in part to a lack of robust enzymatic conversion assays. In contrast to the protein lysine deacetylase and deacylase activity of most other HDAC subtypes, it has recently been shown that HDAC10 has strong preferences for deacetylation of oligoamine substrates like acetyl-putrescine or -spermidine. Hence, it is also termed a polyamine deacetylase (PDAC). Here, we present the first fluorescent enzymatic conversion assay for HDAC10 using an aminocoumarin-labelled acetyl-spermidine derivative to measure its PDAC activity, which is suitable for high-throughput screening. Using this assay, we identified potent inhibitors of HDAC10-mediated spermidine deacetylation in vitro. Based on the oligoamine preference of HDAC10, we also designed inhibitors with a basic moiety in appropriate distance to the zinc binding hydroxamate that showed potent inhibition of HDAC10 with high selectivity, and we solved a HDAC10-inhibitor structure using X-ray crystallography. We could demonstrate selective cellular target engagement for HDAC10 but a lysosomal phenotype in neuroblastoma cells that was previously associated with HDAC10 inhibition was not observed. Thus, we have developed new chemical probes for HDAC10 that allow further clarification of the biological role of this enzyme.

Docking, Binding Free Energy Calculations and In Vitro Characterization of Pyrazine Linked 2-Aminobenzamides as Novel Class I Histone Deacetylase (HDAC) Inhibitors.

Class I histone deacetylases, HDAC1, HDAC2, and HDAC3, represent potential targets for cancer treatment. However, the development of isoform-selective drugs for these enzymes remains challenging due to their high sequence and structural similarity. In the current study, we applied a computational approach to predict the selectivity profile of developed inhibitors. Molecular docking followed by MD simulation and calculation of binding free energy was performed for a dataset of 2-aminobenzamides comprising 30 previously developed inhibitors. For each HDAC isoform, a significant correlation was found between the binding free energy values and in vitro inhibitory activities. The predictive accuracy and reliability of the best preforming models were assessed on an external test set of newly designed and synthesized inhibitors. The developed binding free-energy models are cost-effective methods and help to reduce the time required to prioritize compounds for further studies.

Identification of histone deacetylase 10 (HDAC10) inhibitors that modulate autophagy in transformed cells.

Histone deacetylases (HDACs) are a family of 18 epigenetic modifiers that fall into 4 classes. Histone deacetylase inhibitors (HDACi) are valid tools to assess HDAC functions. HDAC6 and HDAC10 belong to the class IIb subgroup of the HDAC family. The targets and biological functions of HDAC10 are ill-defined. This lack of knowledge is due to a lack of specific and potent HDAC10 inhibitors with cellular activity. Here, we have synthesized and characterized piperidine-4-acrylhydroxamates as potent and highly selective inhibitors of HDAC10. This was achieved by targeting the acidic gatekeeper residue Glu274 of HDAC10 with a basic piperidine moiety that mimics the interaction of the polyamine substrate of HDAC10. We have confirmed the binding modes of selected inhibitors using X-ray crystallography. Promising candidates were selected based on their specificity by in vitro profiling using recombinant HDACs. The most promising HDAC10 inhibitors 10c and 13b were tested for specificity in acute myeloid leukemia (AML) cells with the FLT3-ITD oncogene. By immunoblot experiments we assessed the hyperacetylation of histones and tubulin-α, which are class I and HDAC6 substrates, respectively. As validated test for HDAC10 inhibition we used flow cytometry assessing autolysosome formation in neuroblastoma and AML cells. We demonstrate that 10c and 13b inhibit HDAC10 with high specificity over HDAC6 and with no significant impact on class I HDACs. The accumulation of autolysosomes is not a consequence of apoptosis and 10c and 13b are not toxic for normal human kidney cells. These data show that 10c and 13b are nanomolar inhibitors of HDAC10 with high specificity. Thus, our new HDAC10 inhibitors are tools to identify the downstream targets and functions of HDAC10 in cells.

Structure-Activity Relationship of Anti-Mycobacterium abscessus Piperidine-4-carboxamides, a New Class of NBTI DNA Gyrase Inhibitors.

Mycobacterium abscessus causes difficult-to-cure pulmonary infections. The bacterium is resistant to most anti-infective agents, including first line antituberculosis (anti-TB) drugs. MMV688844 (844) is a piperidine-4-carboxamide (P4C) with bactericidal properties against M. abscessus. We recently identified DNA gyrase as the molecular target of 844. Here, we present in silico docking and genetic evidence suggesting that P4Cs display a similar binding mode to DNA gyrase as gepotidacin. Gepotidacin is a member of the Novel Bacterial Topoisomerase Inhibitors (NBTIs), a new class of nonfluoroquinolone DNA gyrase poisons. Thus, our work suggests that P4Cs present a novel structural subclass of NBTI. We describe structure-activity relationship studies of 844 leading to analogues showing increased antibacterial activity. Selected derivatives were tested for their inhibitory activity against recombinant M. abscessus DNA gyrase. Further optimization of the lead structures led to improved stability in mouse plasma and increased oral bioavailability.

Histone Deacetylase (HDAC) Inhibitors for the Treatment of Schistosomiasis.

Schistosomiasis is a major neglected parasitic disease that affects more than 240 million people worldwide and for which the control strategy consists of mass treatment with the only available drug, praziquantel. Schistosomes display morphologically distinct stages during their life cycle and the transformations between stages are controlled by epigenetic mechanisms. The targeting of epigenetic actors might therefore represent the parasites' Achilles' heel. Specifically, histone deacetylases have been recently characterized as drug targets for the treatment of schistosomiasis. This review focuses on the recent development of inhibitors for schistosome histone deacetylases. In particular, advances in the development of inhibitors of Schistosoma mansoni histone deacetylase 8 have indicated that targeting this enzyme is a promising approach for the treatment of this infection.

Synthesis, structure-activity relationships, cocrystallization and cellular characterization of novel smHDAC8 inhibitors for the treatment of schistosomiasis.

Schistosomiasis is a major neglected parasitic disease that affects more than 265 million people worldwide and for which the control strategy consists of mass treatment with the only available drug, praziquantel. In this study, we chemically optimized our previously reported benzhydroxamate-based inhibitors of Schistosoma mansoni histone deacetylase 8 (smHDAC8). Crystallographic analysis provided insights into the inhibition mode of smHDAC8 activity by the highly potent inhibitor 5o. Structure-based optimization of the novel inhibitors was carried out using the available crystal structures as well as docking studies on smHDAC8. The compounds were evaluated in screens for inhibitory activity against schistosome and human HDACs (hHDAC). The in vitro and docking results were used for detailed structure activity relationships. The synthesized compounds were further investigated for their lethality against the schistosome larval stage using a fluorescence-based assay. The most promising inhibitor 5o showed significant dose-dependent killing of the schistosome larvae and markedly impaired egg laying of adult worm pairs maintained in culture.

Evaluation of γ-carboline-phenothiazine conjugates as simultaneous NMDA receptor blockers and cholinesterase inhibitors.

Alzheimer's disease (AD) is the most common form of dementia. It is associated with the impairment of memory and other cognitive functions that are mainly caused by progressive defects in cholinergic and glutamatergic signaling in the central nervous system. Inhibitors of acetylcholinesterase (AChE) and ionotropic glutamate receptors of the N-methyl-d-aspartate (NMDA) receptor family are currently approved as AD therapeutics. We previously showed using a cell-based assay of NMDA receptor-mediated glutamate-induced excitotoxicity that bis-γ-carbolinium conjugates are useful NMDA receptor blockers. However, these compounds also act as subnanomolar AChE inhibitors, which may cause serious anticholinergic side effects when applied in vivo. Here, we evaluated new structures containing γ-carbolines linked to phenothiazine via a propionyl spacer. These compounds were superior to the previously characterized bis-γ-carbolinium conjugates because they blocked NMDA receptors without requiring a quaternary pyridine N-atom and inhibited AChE with moderate IC50 values of 0.54-5.3 µM. In addition, these new compounds displayed considerable selectivity for the inhibition of butyrylcholinesterase (BChE; IC50 = 0.008-0.041 µM), which may be favorable for AD treatment. Inhibitory activities towards the NMDA receptors and AChE were in the same micromolar range, which may be beneficial for equal dosing against multiple targets in AD patients.

Structure-Based Design, Docking and Binding Free Energy Calculations of A366 Derivatives as Spindlin1 Inhibitors.

The chromatin reader protein Spindlin1 plays an important role in epigenetic regulation, through which it has been linked to several types of malignant tumors. In the current work, we report on the development of novel analogs of the previously published lead inhibitor A366. In an effort to improve the activity and explore the structure-activity relationship (SAR), a series of 21 derivatives was synthesized, tested in vitro, and investigated by means of molecular modeling tools. Docking studies and molecular dynamics (MD) simulations were performed to analyze and rationalize the structural differences responsible for the Spindlin1 activity. The analysis of MD simulations shed light on the important interactions. Our study highlighted the main structural features that are required for Spindlin1 inhibitory activity, which include a positively charged pyrrolidine moiety embedded into the aromatic cage connected via a propyloxy linker to the 2-aminoindole core. Of the latter, the amidine group anchor the compounds into the pocket through salt bridge interactions with Asp184. Different protocols were tested to identify a fast in silico method that could help to discriminate between active and inactive compounds within the A366 series. Rescoring the docking poses with MM-GBSA calculations was successful in this regard. Because A366 is known to be a G9a inhibitor, the most active developed Spindlin1 inhibitors were also tested over G9a and GLP to verify the selectivity profile of the A366 analogs. This resulted in the discovery of diverse selective compounds, among which 1s and 1t showed Spindlin1 activity in the nanomolar range and selectivity over G9a and GLP. Finally, future design hypotheses were suggested based on our findings.

Exploring aromatic cage flexibility of the histone methyllysine reader protein Spindlin1 and its impact on binding mode prediction: an in silico study.

Some of the main challenges faced in drug discovery are pocket flexibility and binding mode prediction. In this work, we explored the aromatic cage flexibility of the histone methyllysine reader protein Spindlin1 and its impact on binding mode prediction by means of in silico approaches. We first investigated the Spindlin1 aromatic cage plasticity by analyzing the available crystal structures and through molecular dynamic simulations. Then we assessed the ability of rigid docking and flexible docking to rightly reproduce the binding mode of a known ligand into Spindlin1, as an example of a reader protein displaying flexibility in the binding pocket. The ability of induced fit docking was further probed to test if the right ligand binding mode could be obtained through flexible docking regardless of the initial protein conformation. Finally, the stability of generated docking poses was verified by molecular dynamic simulations. Accurate binding mode prediction was obtained showing that the herein reported approach is a highly promising combination of in silico methods able to rightly predict the binding mode of small molecule ligands in flexible binding pockets, such as those observed in some reader proteins.