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.

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.

Prodrug Approach to Exploit (S) Alanine as Arginine Mimic Moiety in the Development of Protein Arginine Methyltransferase 4 Inhibitors.

Protein arginine methyltransferase (PRMT) 4 (also known as coactivator-associated arginine methyltransferase 1; CARM1) is involved in a variety of biological processes and is considered as an emerging target class in oncology and other diseases. A successful strategy to identify PRMT substrate-competitive inhibitors has been to exploit chemical scaffolds able to mimic the arginine substrate. (S)-Alanine amide moiety is a valuable arginine mimic for the development of potent and selective PRMT4 inhibitors; however, its high hydrophilicity led to derivatives with poor cellular outcomes. Here, we describe the development of PRMT4 inhibitors featuring a central pyrrole core and an alanine amide moiety. Rounds of optimization, aimed to increase lipophilicity and simultaneously preserve the inhibitory activity, produced derivatives that, despite good potency and physicochemical properties, did not achieve on-target effects in cells. On the other hand, masking the amino group with a NAD(P)H:quinone oxidoreductase 1 (NQO1)-responsive trigger group, led to prodrugs able to reduce arginine dimethylation of the PRMT4 substrates BRG1-associated factor 155 (BAF155). These results indicate that prodrug strategies can be successfully applied to alanine-amide containing PRMT4 inhibitors and provide an option to enable such compounds to achieve sufficiently high exposures in vivo.

Repositioning of Small Molecules through the Inverse Virtual Screening in silico Tool: Case of Benzothiazole-Based Inhibitors of Soluble Epoxide Hydrolase (sEH).

Computational techniques accelerate drug discovery by identifying bioactive compounds for specific targets, optimizing molecules with moderate activity, or facilitating the repositioning of inactive items onto new targets. Among them, the Inverse Virtual Screening (IVS) approach is aimed at the evaluation of one or a small set of molecules against a panel of targets for addressing target identification. In this work, a focused library of benzothiazole-based compounds was re-investigated by IVS. Four items, originally synthesized and tested on bromodomain-containing protein 9 (BRD9) but yielding poor binding, were critically re-analyzed, disclosing only a partial fit with 3D structure-based pharmacophore models, which, in the meanwhile, were developed for this target. Afterwards, these compounds were re-evaluated through IVS on a panel of proteins involved in inflammation and cancer, identifying soluble epoxide hydrolase (sEH) as a putative interacting target. Three items were subsequently confirmed as able to interfere with sEH activity, leading to inhibition percentages spanning from 70 % up to 30 % when tested at 10 µM. Finally, one benzothiazole-based compound emerged as the most promising inhibitor featuring an IC50 in the low micromolar range (IC50=6.62±0.13 µM). Our data confirm IVS as a predictive tool for accelerating the target identification and repositioning processes.

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.

NMR Metabolomics and Chemometrics of Commercial Varieties of Phaseolus vulgaris L. Seeds from Italy and In Vitro Antioxidant and Antifungal Activity.

The metabolite fingerprinting of four Italian commercial bean seed cultivars, i.e., Phaseolus Cannellino (PCANN), Controne (PCON), Vellutina (PVEL), and Occhio Nero (PON), were investigated by Nuclear Magnetic Resonance (NMR) spectroscopy and multivariate data analysis. The hydroalcoholic and organic extract analysis disclosed more than 32 metabolites from various classes, i.e., carbohydrates, amino acids, organic acids, nucleosides, alkaloids, and fatty acids. PVEL, PCON, and PCANN varieties displayed similar chemical profiles, albeit with somewhat different quantitative results. The PON metabolite composition was slightly different from the others; it lacked GABA and pipecolic acid, featured a higher percentage of malic acid than the other samples, and showed quantitative variations of several metabolites. The lipophilic extracts from all four cultivars demonstrated the presence of omega-3 and omega-6 unsaturated fatty acids. After the determination of the total phenolic, flavonoids, and condensed tannins content, in vitro antioxidant activity was then assessed using the DPPH scavenging activity, the ABTS scavenging assay, and ferric-reducing antioxidant power (FRAP). Compared to non-dark seeds (PCON, PCANN), brown seeds (PVEL, PON) featured a higher antioxidant capacity. Lastly, only PON extract showed in vitro antifungal activity against the sclerotia growth of S. rolfsii, by inhibiting halo growth by 75%.

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.

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 Anticancer Potential of Premna resinosa (Hochst.) Leaf Surface Extract: Discovering New Diterpenes as Heat Shock Protein 70 (Hsp70) Binding Agents.

Premna, a genus consisting of approximately 200 species, predominantly thrives in tropical and subtropical areas. Many of these species have been utilized in ethnopharmacology for diverse medicinal applications. In Saudi Arabia, Premna resinosa (Hochst.) Schauer (Lamiaceae) grows wildly, and its slightly viscid leaves are attributed to the production of leaf accession. In this study, we aimed to extract the surface accession from fresh leaves using dichloromethane to evaluate the anticancer potential. The plant exudate yielded two previously unknown labdane diterpenes, Premnaresone A and B, in addition to three already described congeners and four known flavonoids. The isolation process was accomplished using a combination of silica gel column chromatography and semi-preparative HPLC, the structures of which were identified by NMR and HRESIMS analyses and a comparison with the literature data of associated compounds. Furthermore, we employed a density functional theory (DFT)/NMR approach to suggest the relative configuration of different compounds. Consequently, we investigated the possibility of developing new chaperone inhibitors by subjecting diterpenes 1-5 to a Surface Plasmon Resonance-screening, based on the knowledge that oridonin, a diterpene, interacts with Heat Shock Protein 70 (Hsp70) 1A in cancer cells. Additionally, we studied the anti-proliferative activity of compounds 1-5 on human Jurkat (human T-cell lymphoma) and HeLa (epithelial carcinoma) cell lines, where diterpene 3 exhibited activity in Jurkat cell lines after 48 h, with an IC50 of 15.21 ± 1.0 µM. Molecular docking and dynamic simulations revealed a robust interaction between compound 3 and Hsp70 key residues.

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.

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.

Identification of 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based small molecules as selective BRD9 binders.

Targeting bromodomain-containing protein 9 (BRD9) represents a promising strategy for the development of new agents endowed with anticancer properties. With this aim, a set of 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based compounds was investigated following a combined approach that relied on in silico studies, chemical synthesis, biophysical and biological evaluation of the most promising items. The protocol was initially based on molecular docking experiments, accounting a library of 1896 potentially synthesizable items tested in silico against the bromodomain of BRD9. A first set of 21 compounds (1-21) was selected and the binding on BDR9 was assessed through AlphaScreen assays. The obtained results disclosed compounds 17 and 20 able to bind BRD9 in the submicromolar range (IC50 = 0.35 ± 0.18 µM and IC50 = 0.14 ± 0.03 µM, respectively) showing a promising selectivity profile when tested against further nine bromodomains. Taking advantage of 3D structure-based pharmacophore models, additional 10 derivatives were selected in silico for the synthetic step and binding assessment, highlighting seven compounds (22, 23, 25, 26, 28, 29, 31) able to selectively bind BRD9 among different bromodomains. The ability of the identified BRD9 binders to cross artificial membranes in vitro was also assessed, revealing a very good passive permeability profile. Preliminary studies were carried out on a panel of healthy and cancer human cell lines to explore the biological behavior of the selected compounds, disclosing a moderate activity and significant selectivity profile towards leukaemia cells. These results highlighted the applicability of the reported multidisciplinary approach for accelerating the selection of promising items and for driving the chemical synthesis of novel selective BRD9 binders. Moreover, the low molecular weight of the reported 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based BRD9 binders suggests the possibility for further exploring the chemical space in order to obtain new analogues with improved potency.

Phytochemical Analysis of the Methanolic Extract and Essential Oil from Leaves of Industrial Hemp Futura 75 Cultivar: Isolation of a New Cannabinoid Derivative and Biological Profile Using Computational Approaches.

Cannabis sativa L. is a plant belonging to the Cannabaceae family, cultivated for its psychoactive cannabinoid (Δ9-THC) concentration or for its fiber and nutrient content in industrial use. Industrial hemp shows a low Δ9-THC level and is a valuable source of phytochemicals, mainly represented by cannabinoids, flavones, terpenes, and alkaloids, with health-promoting effects. In the present study, we investigated the phytochemical composition of leaves of the industrial hemp cultivar Futura 75, a monoecious cultivar commercially used for food preparations or cosmetic purposes. Leaves are generally discarded, and represent waste products. We analyzed the methanol extract of Futura 75 leaves by HPLC and NMR spectroscopy and the essential oil by GC-MS. In addition, in order to compare the chemical constituents, we prepared the water infusion. One new cannabinoid derivative (1) and seven known components, namely, cannabidiol (2), cannabidiolic acid (3), ß-cannabispirol (4), ß-cannabispirol (5), canniprene (6), cannabiripsol (7), and cannflavin B (8) were identified. The content of CBD was highest in all preparations. In addition, we present the outcomes of a computational study focused on elucidating the role of 2α-hydroxy-Δ3,7-cannabitriol (1), CBD (2), and CBDA (3) in inflammation and thrombogenesis.

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.

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.

Carnosol Attenuates LPS-Induced Inflammation of Cardiomyoblasts by Inhibiting NF-κB: A Mechanistic in Vitro and in Silico Study.

Carnosol possesses several beneficial pharmacological properties. However, its role in lipopolysaccharide (LPS) induced inflammation and cardiomyocyte cell line (H9C2) has never been investigated. Therefore, the effect of carnosol and an NF-κB inhibitor BAY 11-7082 was examined, and the underlying role of the NF-κB-dependent inflammatory pathway was analyzed as the target enzyme. Cell viability, inflammatory cytokines levels (tumor necrosis factor-α (TNF-α), interleukin (IL)-1ß, IL-6, and prostaglandin E 2 (PGE2)), and related gene expression (TNF-α, IL-1ß, IL-6, and cyclooxygenase-2 (COX-2)) were analyzed by ELISA and real-time PCR. In addition, docking studies analyzed carnosol's molecular interactions and binding modes to NF-κB and IKK. We report that LPS caused the reduction of cell viability while enhancing both cytokines protein and mRNA levels (P < 0.001, for all cases). However, the BAY 11-7082 pretreatment of the cells and carnosol increased cell viability and reduced cytokine protein and mRNA levels (P < 0.001 vs. LPS, for all cases). Furthermore, our in silico analyses also supported the modulation of NF-κB and IKK by carnosol. This evidence highlights the defensive effects of carnosol against sepsis-induced myocardial dysfunction and, contextually, paved the rationale for the next in vitro and in vivo studies aimed to precisely describe its mechanism(s) of action.

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.

In Silico, In Vitro, and In Vivo Analysis of Tanshinone IIA and Cryptotanshinone from Salvia miltiorrhiza as Modulators of Cyclooxygenase-2/mPGES-1/Endothelial Prostaglandin EP3 Pathway.

Tanshinone IIA (TIIA) and cryptotanshinone (CRY) from Salvia miltiorrhiza Bunge were investigated for their inhibitory activity against the cyclooxygenase-2 (COX-2)/microsomal prostaglandin E synthase-1 (mPGES-1)/endothelial prostaglandin 3 (EP3) pathway using in silico, in vitro, in vivo, and ex vivo assays. From the analysis of the docking poses, both diterpenoids were able to interact significantly with COX-2, 5-lipoxygenase (5-LO), platelet-activating factor receptor (PAFR), and mPGES-1. This evidence was further corroborated by data obtained from a cell-free assay, where CRY displayed a significant inhibitory potency against mPGES-1 (IC50 = 1.9 ± 0.4 µM) and 5-LO (IC50 = 7.1 µM), while TIIA showed no relevant inhibition of these targets. This was consistent with their activity to increase mice bleeding time (CRY: 2.44 ± 0.13 min, p ≤ 0.001; TIIA: 2.07 ± 0.17 min p ≤ 0.01) and with the capability to modulate mouse clot retraction (CRY: 0.048 ± 0.011 g, p ≤ 0.01; TIIA: 0.068 ± 0.009 g, p ≤ 0.05). For the first time, our results show that TIIA and, in particular, CRY are able to interact significantly with the key proteins involved not only in the onset of inflammation but also in platelet activity (and hyper-reactivity). Future preclinical and clinical investigations, together with this evidence, could provide the scientific basis to consider these compounds as an alternative therapeutic approach for thrombotic- and thromboembolic-based diseases.

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.

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.