PharmaCore: The Automatic Generation of 3D Structure-Based Pharmacophore Models from Protein/Ligand Complexes.

In this work, we present PharmaCore: a new, completely automatic workflow aimed at generating three-dimensional (3D) structure-based pharmacophore models toward any target of interest. The proposed approach relies on using cocrystallized ligands to create the input files for generating the pharmacophore hypotheses, integrating not only the three-dimensional structural information on the ligand but also data concerning the binding mode of these molecules put in the protein cavity. We developed a Python library that, starting from the specific UniProt ID of the protein under investigation as the only element that requires user intervention, subsequently collects and aligns the corresponding structures bearing a known ligand in a fully automated fashion, bringing them all into the same coordinate system. The protocol includes a final phase in which the aligned small molecules are used to produce the pharmacophore hypotheses directly onto the protein structure using a specific software, e.g., Phase (Schrödinger LLC). To validate the entire procedure and highlight the possible applications in the field of drug discovery and repositioning, we first generated pharmacophores for soluble epoxide hydrolase (sEH) and compared with already-published ones. Then, we reproduced the binding profile of a reported selective binder of ATAD2 bromodomain (AM879), testing it against a panel of 1741 pharmacophores related to 16 epigenetic proteins and automatically generated with PharmaCore, finally disclosing putative unprecedented off-targets. The computational predictions were successfully validated with AlphaScreen assays, highlighting the applicability of the proposed workflow in drug discovery and repositioning. Finally, the process was also validated on tankyrase 2 and SARS-CoV-2 MPro, confirming the robustness of PharmaCore.

Chemoproteomics Reveals USP5 (Ubiquitin Carboxyl-Terminal Hydrolase 5) as Promising Target of the Marine Polyketide Gracilioether A.

Mass spectrometry-based chemical proteomic approaches using limited proteolysis have become a powerful tool for the identification and analysis of the interactions between a small molecule (SM) and its protein target(s). Gracilioether A (GeA) is a polyketide isolated from a marine sponge, for which we aimed to trace the interactome using this strategy. DARTS (Drug Affinity Responsive Target Stability) and t-LiP-MS (targeted-Limited Proteolysis-Mass Spectrometry) represented the main techniques used in this study. DARTS was applied on HeLa cell lysate for the identification of the GeA target proteins, and t-LiP-MS was employed to investigate the protein's regions involved in the binding with GeA. The results were complemented through the use of binding studies using Surface Plasmon Resonance (SPR) and in silico molecular docking experiments. Ubiquitin carboxyl-terminal hydrolase 5 (USP5) was identified as a promising target of GeA, and the interaction profile of the USP5-GeA complex was explained. USP5 is an enzyme involved in the pathway of protein metabolism through the disassembly of the polyubiquitin chains on degraded proteins into ubiquitin monomers. This activity is connected to different cellular functions concerning the maintenance of chromatin structure and receptors and the degradation of abnormal proteins and cancerogenic progression. On this basis, this structural information opens the way to following studies focused on the definition of the biological potential of Gracilioether A and the rational development of novel USP5 inhibitors based on a new structural skeleton.

Fatty Acid Synthase as Interacting Anticancer Target of the Terpenoid Myrianthic Acid Disclosed by MS-Based Proteomics Approaches.

This research focuses on the target deconvolution of the natural compound myrianthic acid, a triterpenoid characterized by an ursane skeleton isolated from the roots of Myrianthus arboreus and from Oenothera maritima Nutt. (Onagraceae), using MS-based chemical proteomic techniques. Application of drug affinity responsive target stability (DARTS) and targeted-limited proteolysis coupled to mass spectrometry (t-LiP-MS) led to the identification of the enzyme fatty acid synthase (FAS) as an interesting macromolecular counterpart of myrianthic acid. This result, confirmed by comparison with the natural ursolic acid, was thoroughly investigated and validated in silico by molecular docking, which gave a precise picture of the interactions in the MA/FAS complex. Moreover, biological assays showcased the inhibitory activity of myrianthic acid against the FAS enzyme, most likely related to its antiproliferative activity towards tumor cells. Given the significance of FAS in specific pathologies, especially cancer, the myrianthic acid structural moieties could serve as a promising reference point to start the potential development of innovative approaches in therapy.

Target identification by structure-based computational approaches: Recent advances and perspectives.

The use of computational techniques in the early stages of drug discovery has recently experienced a boost, especially in the target identification step. Finding the biological partner(s) for new or existing synthetic and/or natural compounds by "wet" approaches may be challenging; therefore, preliminary in silico screening is even more recommended. After a brief overview of some of the most known target identification techniques, recent advances in structure-based computational approaches for target identification are reported in this digest, focusing on Inverse Virtual Screening and its recent applications. Moreover, future perspectives concerning the use of such methodologies, coupled or not with other approaches, are analyzed.

2-Substituted 1,5-benzothiazepine-based HDAC inhibitors exert anticancer activities on human solid and acute myeloid leukemia cell lines.

Herein, we report the development of a new series of histone deacetylase inhibitors (HDACi) containing a 2-substituted 1,5-benzothiazepine scaffold. First, a virtual combinatorial library (∼1.6 × 103 items) was built according to a convenient synthetic route, and then it was submitted to molecular docking experiments on seven HDACs isoforms belonging to classes I and II. Integrated computational filters were used to select the most promising ones that were synthesized through an optimized approach, also amenable to generating both racemic and enantioenriched benzothiazepine-based derivatives. The obtained compounds showed potent HDAC inhibitory activity, especially those containing the sulphone moiety, endowed with IC50 in the nanomolar range. In addition, in vitro outcomes of our synthesized compounds demonstrated a cytotoxic effect on U937 and HCT116 cell lines and an arrest in the G2/M phase (13 ≤ IC50 ≤ 18 µM). Finally, Western blot analyses outlined the modulation of the histone acetyl markers such as H3K9/14, acetyl-tubulin, and the apoptotic indicator p21 in both cancer cell lines, disclosing a good HDAC inhibitor activity exerted by the designed items. Given the key role of HDACs in many cellular pathways, which makes these enzymes appealing and "hot" drug targets, our findings highlighted the importance of these 2-substituted 1,5-benzothiazepine-based compounds (both in the reduced and oxidized version) for the development of novel epidrugs.

Exploring the chemical space of functionalized [1,2,4]triazolo[4,3-a]quinoxaline-based compounds targeting the bromodomain of BRD9.

Here we report a detailed structure-activity relationship (SAR) study related to [1,2,4]triazolo[4,3-a]quinoxaline-based compounds targeting the reader module of bromodomain containing-protein 9 (BRD9). 3D structure-based pharmacophore models, previously introduced by us, were here employed to evaluate a second generation of compounds, exploring different substitution patterns on the heterocyclic core. Starting from the promising data obtained from our previously identified [1,2,4]triazolo[4,3-a]quinoxaline-based compounds 1-4, the combination of in silico studies, chemical synthesis, biophysical and in vitro assays led to the identification of a new set of derivatives, selected for thoroughly exploring the chemical space of the bromodomain binding site. In more details, the investigation of different linkers at C-4 position highlighted the amine spacer as mandatory for the binding with the protein counterpart and the crucial role of the alkyl substituents at C-1 for increasing the selectivity toward BRD9. Additionally, the importance of a hydrogen bond donor group, critical to anchor the ZA region and required for the interaction with Ile53 residue, was inferred from the analysis of our collected results. Herein we also propose an optimization and an update of our previously reported "pharm-druglike2" 3D structure-based pharmacophore model, introducing it as "pharm-druglike2.1". Compounds 24-26, 32, 34 and 36 were identified as new valuable BRD9 binders featuring IC50 values in the low micromolar range. Among them, 24 and 36 displayed an excellent selectivity towards BRD9 and a good antiproliferative effect on a panel of leukemia models, especially toward CCRF-CEM cell line, with no cytotoxicity on healthy cells. Notably, the interaction of 24 and 36 with the bromodomain and PHD finger-containing protein 1 (BRPF1) also emerged, disclosing them as new and unexplored dual inhibitors for these two proteins highly involved in leukemia. These findings highlight the potential for the identification of new attractive dual epidrugs as well as a promising starting point for the development of chemical degraders endowed with anticancer activities.

Diterpenoid Constituents of Psiadia punctulata and Evaluation of Their Antimicrobial Activity.

Sixteen diterpenes (1-16), along with 10 previously described compounds, including four flavonoids and six diterpenes, were isolated from the aerial parts of Psiadia punctulata growing in Saudi Arabia. The diterpene structures were elucidated using NMR spectroscopy and mass spectrometry data. Furthermore, a DFT/NMR procedure was used to suggest the relative configuration of several compounds. The labdane-derived skeletons, namely, ent-atisane, ent-beyerene, ent-trachylobane, and ent-kaurene, were identified. The extracts, fractions, and pure compounds were then tested against Staphylococcus aureus, Streptococcus mutans, Treponema denticola, and Lactobacillus plantarum. One diterpenoid, namely, psiadin, showed an additive effect with the antiseptic chlorhexidine, with a fractional inhibitory concentration index of less than 1. Additionally, psiadin showed a prospective inhibition activity for bacterial efflux pumps.

Cytotoxic Sesquiterpenoids from Ammoides atlantica Aerial Parts.

Seven new terpenoids, namely, guaiane (1-4), eudesmane (5), and bisabolane (6) sesquiterpenoids and a furanone (7), were isolated from the aerial parts of Ammoides atlantica, a herbaceous plant growing in Algeria, together with eight known compounds. All metabolites were characterized by their 1D and 2D NMR and HRESIMS data. A combined DFT/NMR method was applied to study the relative configurations of 1-4, 6, and 7. All compounds, except 2, were assayed against MCF-7, A375, A549, HaCaT, and Jurkat cell lines. Compounds 8, 10, and 11 induced a dose-dependent reduction in cell viability with different potency on almost all cell lines used. The most active compounds, 8 and 10, were studied to assess their potential apoptotic effects and cell cycle inhibition.

Introducing structure-based three-dimensional pharmacophore models for accelerating the discovery of selective BRD9 binders.

A well-structured in silico workflow is here reported for disclosing structure-based pharmacophore models against bromodomain-containing protein 9 (BRD9), accelerating virtual screening campaigns and facilitating the identification of novel binders. Specifically, starting from 23 known ligands co-crystallized with BRD9, three-dimensional pharmacophore models, namely placed in a reference protein structure, were developed. Specifically, we here introduce a fragment-related pharmacophore model, useful for the identification of new promising small chemical probes targeting the protein region responsible of the acetyllysine recognition, and two further pharmacophore models useful for the selection of compounds featuring drug-like properties. A pharmacophore-driven virtual screening campaign was then performed to facilitate the selection of new selective BRD9 ligands, starting from a large library of commercially available molecules. The identification of a promising BRD9 binder (7) prompted us to re-iterate this computational workflow on a second focused in-house built library of synthesizable compounds and, eventually, three further novel BRD9 binders were disclosed (8-10). Moreover, all these compounds were tested among a panel comprising other nine bromodomains, showing a high selectivity for BRD9. Preclinical bioscreens for potential anticancer activity highlighted compound 7 as that showing the most promising biological effects, proving the reliability of this in silico pipeline and confirming the applicability of the here introduced structure-based three-dimensional (3D) pharmacophore models as straightforward tools for the selection of new BRD9 ligands.

Repositioning of Quinazolinedione-Based Compounds on Soluble Epoxide Hydrolase (sEH) through 3D Structure-Based Pharmacophore Model-Driven Investigation.

The development of new bioactive compounds represents one of the main purposes of the drug discovery process. Various tools can be employed to identify new drug candidates against pharmacologically relevant biological targets, and the search for new approaches and methodologies often represents a critical issue. In this context, in silico drug repositioning procedures are required even more in order to re-evaluate compounds that already showed poor biological results against a specific biological target. 3D structure-based pharmacophoric models, usually built for specific targets to accelerate the identification of new promising compounds, can be employed for drug repositioning campaigns as well. In this work, an in-house library of 190 synthesized compounds was re-evaluated using a 3D structure-based pharmacophoric model developed on soluble epoxide hydrolase (sEH). Among the analyzed compounds, a small set of quinazolinedione-based molecules, originally selected from a virtual combinatorial library and showing poor results when preliminarily investigated against heat shock protein 90 (Hsp90), was successfully repositioned against sEH, accounting the related built 3D structure-based pharmacophoric model. The promising results here obtained highlight the reliability of this computational workflow for accelerating the drug discovery/repositioning processes.

Limonoids from Guarea guidonia and Cedrela odorata: Heat Shock Protein 90 (Hsp90) Modulator Properties of Chisomicine D.

Nine new limonoids (1-9) were isolated from the stem bark of Guarea guidonia (1-4) and Cedrela odorata (5-9). Their structures were elucidated using 1D and 2D NMR and MS data and chemical methods as three A2,B,D-seco-type limonoids (1-3), a mexicanolide (4), three nomilin-type (5-7) limonoids, and two limonol derivatives (8 and 9). A DFT/NMR procedure was used to define the relative configurations of 1 and 3. A surface plasmon resonance approach was used to screen the Hsp90 binding capability of the limonoids, and the A2,B,D-seco-type limonoid 8-hydro-(8S*,9S*)-dihydroxy-14,15-en-chisomicine A, named chisomicine D (1), demonstrated the highest affinity. By means of mass spectrometry data, biochemical and cellular assays, and molecular docking, 1 was found as a type of client-selective Hsp90 inhibitor binding to the C-terminus domain of the chaperone.

Insights into the Ligand Binding to Bromodomain-Containing Protein 9 (BRD9): A Guide to the Selection of Potential Binders by Computational Methods.

The estimation of the binding of a set of molecules against BRD9 protein was carried out through an in silico molecular dynamics-driven exhaustive analysis to guide the identification of potential novel ligands. Starting from eight crystal structures of this protein co-complexed with known binders and one apo form, we conducted an exhaustive molecular docking/molecular dynamics (MD) investigation. To balance accuracy and an affordable calculation time, the systems were simulated for 100 ns in explicit solvent. Moreover, one complex was simulated for 1 µs to assess the influence of simulation time on the results. A set of MD-derived parameters was computed and compared with molecular docking-derived and experimental data. MM-GBSA and the per-residue interaction energy emerged as the main indicators for the good interaction between the specific binder and the protein counterpart. To assess the performance of the proposed analysis workflow, we tested six molecules featuring different binding affinities for BRD9, obtaining promising outcomes. Further insights were reported to highlight the influence of the starting structure on the molecular dynamics simulations evolution. The data confirmed that a ranking of BRD9 binders using key parameters arising from molecular dynamics is advisable to discard poor ligands before moving on with the synthesis and the biological tests.

Elucidating heteroatom influence on homonuclear 4 J(H,H) coupling constants by DFT/NMR approach.

We report the structural dependency of long range scalar J-coupling constant across four bonds as function of the dihedral angles Φ1 and Φ3. The calculated homonuclear coupling constants 4 J(H,H ), obtained at a density functional theory level, were measured between C(1)─X(2) and X(2)─C(3) bonds in three-term models, where C, N, O, and S were systematically used as the second atom of the alkyl structures (1-4). The 4 J(H,H) calculated values, tabulated for variation of 30° for both Φ1 and Φ3, have disclosed an unexpected detectable coupling constant (4 J(H,H ) ≥ 1 Hz) across heteroatoms, useful to provide valuable structural information. A 2-methyl-1,3-dithiane sulfide (5) was used as a case study to prove the applicability and reliability of the calculated values to real issues. The 4 J(H,H ) values obtained at density functional theory for the system 4 have reproduced with good accuracy an unexpected experimental 4 J(H2ax-H4ax ) = 1.01 Hz of sulfide molecule (5), suggesting these calculated coupling constant values as a new powerful tool for the organic synthesis and stereochemical analysis.

Fusicoccane Diterpenes from Hypoestes forsskaolii as Heat Shock Protein 90 (Hsp90) Modulators.

Ten new (1-10) and six known (11-16) fusicoccane diterpenes were isolated from the roots of Hypoestes forsskaolii. The structural characterization of 1-10 was performed by spectroscopic analysis, including 1D and 2D NMR, ECD, and HRESIMS experiments. From a perspective of obtaining potential Hsp90α inhibitors, the isolates were screened by surface plasmon resonance measurements and their cytotoxic activity was assayed using Jurkat and HeLa cancer cells. Compound 6, 18-hydroxyhypoestenone, was shown to be the most active compound against Hsp90, and its interactions were studied also by biochemical and cellular assays and by molecular docking.

Guaianolides from Ormenis mixta: Structural Insights and Evaluation of Their Anti-inflammatory Profile.

In this paper, the isolation of five new guaianolides (1:  - 5: ) and four (6:  - 9: ) known sesquiterpenes from Ormenis mixta aerial parts is reported. The structural determination of the guaianolides was obtained by NMR spectroscopic data, as well as MS experiments. Their relative configurations were assigned by a combined quantum mechanical/NMR approach, comparing the experimental 13C/1H NMR chemical shift data and 1 J H-H homonuclear coupling constants with the related predicted values. The isolates were assayed for their anti-inflammatory potential evaluating nitric oxide release and cyclooxygenase-2 expression in J774A.1 macrophages treated with lipopolysaccharide from Escherichia coli. Our results indicated that, among the tested compounds, 1:  - 3: , and 7: were able to inhibit nitric oxide release, while all were able to inhibit cyclooxygenase-2 expression with different potencies.

Long-Lasting Anti-Inflammatory and Antinociceptive Effects of Acute Ammonium Glycyrrhizinate Administration: Pharmacological, Biochemical, and Docking Studies.

The object of the study was to estimate the long-lasting effects induced by ammonium glycyrrhizinate (AG) after a single administration in mice using animal models of pain and inflammation together with biochemical and docking studies. A single intraperitoneal injection of AG was able to produce anti-inflammatory effects in zymosan-induced paw edema and peritonitis. Moreover, in several animal models of pain, such as the writhing test, the formalin test, and hyperalgesia induced by zymosan, AG administered 24 h before the tests was able to induce a strong antinociceptive effect. Molecular docking studies revealed that AG possesses higher affinity for microsomal prostaglandin E synthase type-2 compared to type-1, whereas it seems to locate better in the binding pocket of cyclooxygenase (COX)-2 compared to COX-1. These results demonstrated that AG induced anti-inflammatory and antinociceptive effects until 24-48 h after a single administration thanks to its ability to bind the COX/mPGEs pathway. Taken together, all these findings highlight the potential use of AG for clinical treatment of pain and/or inflammatory-related diseases.

Down regulation of pro-inflammatory pathways by tanshinone IIA and cryptotanshinone in a non-genetic mouse model of Alzheimer's disease.

Alzheimer's disease (AD) is a common form of dementia mainly characterized by the deposition of neurofibrillary tangles and ß-amyloid (Aß) peptides in the brain. Additionally, increasing evidence demonstrates that a neuro-inflammatory state plays a key role in the development of this disease. Beside synthetic drugs, the use of natural compounds represents an alternative for the development of new potential drugs for the treatment of AD. Among these, the root of Salvia miltiorhiza Bunge (also known as Danshen) used for the treatment of cardiovascular, cerebrovascular disease and CNS functional decline in Chinese traditional medicine is one of the most representative examples. We therefore evaluated the effects of tanshinone IIA (TIIA) and cryptotanshinone (CRY) (the two major lipophilic compounds of Danshen) in a non-genetic mouse model of ß-amyloid (Aß)-induced AD, which is mainly characterized by reactive gliosis and neuro-inflammation in the brain. To this aim, mice were injected intracerebroventricularly (i.c.v.) with Aß1-42 peptide (3µg/3µl) and after with TIIA and CRY (1, 3, or 10mg/kg) intraperitoneally (i.p.) 3 times weekly for 21days following the induction of experimental AD. Spatial working memory was assessed as a measure of short-term memory in mice, whereas the level of GFAP, S100ß, COX-2, iNOS and NF-kBp65 monitored by western blot and ELISA assay, were selected as markers of reactive gliosis and neuro-inflammation. Finally, by docking studies, the modulation of key pro-inflammatory enzymes and pathways involved in the AD-related neuro-inflammation were also investigated. Results indicate that TIIA and CRY alleviate memory decline in Aß1-42-injected mice, in a dose dependent manner. Moreover, the analysis of gliosis-related and neuro-inflammatory markers in the hippocampal tissues reveal a remarkable reduction in the expression of GFAP, S100ß, COX-2, iNOS and NF-kBp65 after CRY (10mg/kg) treatment. These effects were less evident, but still significant, after TIIA (10mg/kg). Finally, in silico analysis also revealed that both compounds were able to interact with the binding sites of NF-kBp65 endorsing the data from biochemical analysis. We conclude that TIIA and CRY display anti-inflammatory and neuroprotective effect in a non-genetic mouse model of AD, thus playing a role in slowing down the course and onset of AD.

Identification and mechanism of action analysis of the new PARP-1 inhibitor 2″-hydroxygenkwanol A.

BACKGROUND: Poly(ADP-ribose) polymerase 1 (PARP-1) activity has been implicated in the pathogenesis of numerous diseases as cancer, inflammation, diabetes and neurodegenerative disorders, therefore the research for new PARP-1 inhibitors is still an active area. METHODS: To identify new potential PARP-1 inhibitors, we performed a screening of a small-molecule library consisting of polyphenols isolated from plants used in the traditional medicine, by Surface Plasmon Resonance (SPR). Biochemical and cellular assays were performed to confirm SPR results and select the promising candidate(s). Finally, limited proteolysis and ligand docking analyses allowed defining the protein region involved in the interaction with the putative inhibitor(s). RESULTS: The dimeric spiro-flavonoid 2″-hydroxygenkwanol A, member of a relatively recently discovered class of flavonoids containing a spirane C-atom, has been identified as possible PARP-1 inhibitor. This compound showed a high affinity for the polymerase (KD: 0.32±0.05µM); moreover PARP-1 activity in the presence of 2″-hydroxygenkwanol A was significantly affected both when using the recombinant protein and when measuring the cellular effects. Finally, our study suggests this compound to efficiently interact with the protein catalytic domain, into the nicotine binding pocket. CONCLUSION: 2″-hydroxygenkwanol A efficiently binds and inhibits PARP-1 at submicromolar concentrations, thus representing a promising lead for the design of a new class of PARP-1 modulators, useful as therapeutic agents and/or biochemical tools. GENERAL SIGNIFICANCE: Our study has identified an additional class of plant molecules, the spiro-biflavonoids, with known beneficial pharmacological properties but with an unknown mechanism of action, as a possible novel class of PARP-1 activity inhibitors.

Can Small Chemical Modifications of Natural Pan-inhibitors Modulate the Biological Selectivity? The Case of Curcumin Prenylated Derivatives Acting as HDAC or mPGES-1 Inhibitors.

Curcumin, or diferuloylmethane, a polyphenolic molecule isolated from the rhizome of Curcuma longa, is reported to modulate multiple molecular targets involved in cancer and inflammatory processes. On the basis of its pan-inhibitory characteristics, here we show that simple chemical modifications of the curcumin scaffold can regulate its biological selectivity. In particular, the curcumin scaffold was modified with three types of substituents at positions C-1, C-8, and/or C-8' [C5 (isopentenyl, 5-8), C10 (geranyl, 9-12), and C15 (farnesyl, 13, 14)] in order to make these molecules more selective than the parent compound toward two specific targets: histone deacetylase (HDAC) and microsomal prostaglandin E2 synthase-1 (mPGES-1). From combined in silico and in vitro analyses, three selective inhibitors by proper substitution at position 8 were revealed. Compound 13 has improved HDAC inhibitory activity and selectivity with respect to the parent compound, while 5 and 9 block the mPGES-1 enzyme. We hypothesize about the covalent interaction of curcumin, 5, and 9 with the mPGES-1 binding site.

Identification and Development of BRD9 Chemical Probes.

The development of BRD9 inhibitors involves the design and synthesis of molecules that can specifically bind the BRD9 protein, interfering with the function of the chromatin-remodeling complex ncBAF, with the main advantage of modulating gene expression and controlling cellular processes. Here, we summarize the work conducted over the past 10 years to find new BRD9 binders, with an emphasis on their structure-activity relationships, efficacies, and selectivities in preliminary studies. BRD9 is expressed in a variety of cancer forms, hence, its inhibition holds particular significance in cancer research. However, it is crucial to note that the expanding research in the field, particularly in the development of new degraders, may uncover new therapeutic potentials.