Combining Different Docking Engines and Consensus Strategies to Design and Validate Optimized Virtual Screening Protocols for the SARS-CoV-2 3CL Protease.

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

Binding Mode Exploration of B1 Receptor Antagonists' by the Use of Molecular Dynamics and Docking Simulation-How Different Target Engagement Can Determine Different Biological Effects.

The kinin B1 receptor plays a critical role in the chronic phase of pain and inflammation. The development of B1 antagonists peaked in recent years but almost all promising molecules failed in clinical trials. Little is known about these molecules' mechanisms of action and additional information will be necessary to exploit the potential of the B1 receptor. With the aim of contributing to the available knowledge of the pharmacology of B1 receptors, we designed and characterized a novel class of allosteric non-peptidic inhibitors with peculiar binding characteristics. Here, we report the binding mode analysis and pharmacological characterization of a new allosteric B1 antagonist, DFL20656. We analyzed the binding of DFL20656 by single point mutagenesis and radioligand binding assays and we further characterized its pharmacology in terms of IC50, B1 receptor internalization and in vivo activity in comparison with different known B1 antagonists. We highlighted how different binding modes of DFL20656 and a Merck compound (compound 14) within the same molecular pocket can affect the biological and pharmacological properties of B1 inhibitors. DFL20656, by its peculiar binding mode, involving tight interactions with N114, efficiently induced B1 receptor internalization and evoked a long-lasting effect in an in vivo model of neuropathic pain. The pharmacological characterization of different B1 antagonists highlighted the effects of their binding modes on activity, receptor occupancy and internalization. Our results suggest that part of the failure of most B1 inhibitors could be ascribed to a lack of knowledge about target function and engagement.

ADPredict: ADP-ribosylation site prediction based on physicochemical and structural descriptors.

Motivation: ADP-ribosylation is a post-translational modification (PTM) implicated in several crucial cellular processes, ranging from regulation of DNA repair and chromatin structure to cell metabolism and stress responses. To date, a complete understanding of ADP-ribosylation targets and their modification sites in different tissues and disease states is still lacking. Identification of ADP-ribosylation sites is required to discern the molecular mechanisms regulated by this modification. This motivated us to develop a computational tool for the prediction of ADP-ribosylated sites. Results: Here, we present ADPredict, the first dedicated computational tool for the prediction of ADP-ribosylated aspartic and glutamic acids. This predictive algorithm is based on (i) physicochemical properties, (ii) in-house designed secondary structure-related descriptors and (iii) three-dimensional features of a set of human ADP-ribosylated proteins that have been reported in the literature. ADPredict was developed using principal component analysis and machine learning techniques; its performance was evaluated both internally via intensive bootstrapping and in predicting two external experimental datasets. It outperformed the only other available ADP-ribosylation prediction tool, ModPred. Moreover, a novel secondary structure descriptor, HM-ratio, was introduced and successfully contributed to the model development, thus representing a promising tool for bioinformatics studies, such as PTM prediction. Availability and implementation: ADPredict is freely available at www.ADPredict.net. Supplementary information: Supplementary data are available at Bioinformatics online.

"DompeKeys": a set of novel substructure-based descriptors for efficient chemical space mapping, development and structural interpretation of machine learning models, and indexing of large databases.

The conversion of chemical structures into computer-readable descriptors, able to capture key structural aspects, is of pivotal importance in the field of cheminformatics and computer-aided drug design. Molecular fingerprints represent a widely employed class of descriptors; however, their generation process is time-consuming for large databases and is susceptible to bias. Therefore, descriptors able to accurately detect predefined structural fragments and devoid of lengthy generation procedures would be highly desirable. To meet additional needs, such descriptors should also be interpretable by medicinal chemists, and suitable for indexing databases with trillions of compounds. To this end, we developed-as integral part of EXSCALATE, Dompé's end-to-end drug discovery platform-the DompeKeys (DK), a new substructure-based descriptor set, which encodes the chemical features that characterize compounds of pharmaceutical interest. DK represent an exhaustive collection of curated SMARTS strings, defining chemical features at different levels of complexity, from specific functional groups and structural patterns to simpler pharmacophoric points, corresponding to a network of hierarchically interconnected substructures. Because of their extended and hierarchical structure, DK can be used, with good performance, in different kinds of applications. In particular, we demonstrate how they are very well suited for effective mapping of chemical space, as well as substructure search and virtual screening. Notably, the incorporation of DK yields highly performing machine learning models for the prediction of both compounds' activity and metabolic reaction occurrence. The protocol to generate the DK is freely available at https://dompekeys.exscalate.eu and is fully integrated with the Molecular Anatomy protocol for the generation and analysis of hierarchically interconnected molecular scaffolds and frameworks, thus providing a comprehensive and flexible tool for drug design applications.

Identification and characterization of a new potent inhibitor targeting CtBP1/BARS in melanoma cells.

BACKGROUND: The C-terminal-binding protein 1/brefeldin A ADP-ribosylation substrate (CtBP1/BARS) acts both as an oncogenic transcriptional co-repressor and as a fission inducing protein required for membrane trafficking and Golgi complex partitioning during mitosis, hence for mitotic entry. CtBP1/BARS overexpression, in multiple cancers, has pro-tumorigenic functions regulating gene networks associated with "cancer hallmarks" and malignant behavior including: increased cell survival, proliferation, migration/invasion, epithelial-mesenchymal transition (EMT). Structurally, CtBP1/BARS belongs to the hydroxyacid-dehydrogenase family and possesses a NAD(H)-binding Rossmann fold, which, depending on ligands bound, controls the oligomerization of CtBP1/BARS and, in turn, its cellular functions. Here, we proposed to target the CtBP1/BARS Rossmann fold with small molecules as selective inhibitors of mitotic entry and pro-tumoral transcriptional activities. METHODS: Structured-based screening of drug databases at different development stages was applied to discover novel ligands targeting the Rossmann fold. Among these identified ligands, N-(3,4-dichlorophenyl)-4-{[(4-nitrophenyl)carbamoyl]amino}benzenesulfonamide, called Comp.11, was selected for further analysis. Fluorescence spectroscopy, isothermal calorimetry, computational modelling and site-directed mutagenesis were employed to define the binding of Comp.11 to the Rossmann fold. Effects of Comp.11 on the oligomerization state, protein partners binding and pro-tumoral activities were evaluated by size-exclusion chromatography, pull-down, membrane transport and mitotic entry assays, Flow cytometry, quantitative real-time PCR, motility/invasion, and colony assays in A375MM and B16F10 melanoma cell lines. Effects of Comp.11 on tumor growth in vivo were analyzed in mouse tumor model. RESULTS: We identify Comp.11 as a new, potent and selective inhibitor of CtBP1/BARS (but not CtBP2). Comp.11 directly binds to the CtBP1/BARS Rossmann fold affecting the oligomerization state of the protein (unlike other known CtBPs inhibitors), which, in turn, hinders interactions with relevant partners, resulting in the inhibition of both CtBP1/BARS cellular functions: i) membrane fission, with block of mitotic entry and cellular secretion; and ii) transcriptional pro-tumoral effects with significantly hampered proliferation, EMT, migration/invasion, and colony-forming capabilities. The combination of these effects impairs melanoma tumor growth in mouse models.  CONCLUSIONS: This study identifies a potent and selective inhibitor of CtBP1/BARS active in cellular and melanoma animal models revealing new opportunities to study the role of CtBP1/BARS in tumor biology and to develop novel melanoma treatments.

ProfhEX: AI-based platform for small molecules liability profiling.

Off-target drug interactions are a major reason for candidate failure in the drug discovery process. Anticipating potential drug's adverse effects in the early stages is necessary to minimize health risks to patients, animal testing, and economical costs. With the constantly increasing size of virtual screening libraries, AI-driven methods can be exploited as first-tier screening tools to provide liability estimation for drug candidates. In this work we present ProfhEX, an AI-driven suite of 46 OECD-compliant machine learning models that can profile small molecules on 7 relevant liability groups: cardiovascular, central nervous system, gastrointestinal, endocrine, renal, pulmonary and immune system toxicities. Experimental affinity data was collected from public and commercial data sources. The entire chemical space comprised 289'202 activity data for a total of 210'116 unique compounds, spanning over 46 targets with dataset sizes ranging from 819 to 18896. Gradient boosting and random forest algorithms were initially employed and ensembled for the selection of a champion model. Models were validated according to the OECD principles, including robust internal (cross validation, bootstrap, y-scrambling) and external validation. Champion models achieved an average Pearson correlation coefficient of 0.84 (SD of 0.05), an R2 determination coefficient of 0.68 (SD = 0.1) and a root mean squared error of 0.69 (SD of 0.08). All liability groups showed good hit-detection power with an average enrichment factor at 5% of 13.1 (SD of 4.5) and AUC of 0.92 (SD of 0.05). Benchmarking against already existing tools demonstrated the predictive power of ProfhEX models for large-scale liability profiling. This platform will be further expanded with the inclusion of new targets and through complementary modelling approaches, such as structure and pharmacophore-based models. ProfhEX is freely accessible at the following address: https://profhex.exscalate.eu/ .

DHFR Inhibitors Display a Pleiotropic Anti-Viral Activity against SARS-CoV-2: Insights into the Mechanisms of Action.

During the COVID-19 pandemic, drug repurposing represented an effective strategy to obtain quick answers to medical emergencies. Based on previous data on methotrexate (MTX), we evaluated the anti-viral activity of several DHFR inhibitors in two cell lines. We observed that this class of compounds showed a significant influence on the virus-induced cytopathic effect (CPE) partly attributed to the intrinsic anti-metabolic activity of these drugs, but also to a specific anti-viral function. To elucidate the molecular mechanisms, we took advantage of our EXSCALATE platform for in-silico molecular modelling and further validated the influence of these inhibitors on nsp13 and viral entry. Interestingly, pralatrexate and trimetrexate showed superior effects in counteracting the viral infection compared to other DHFR inhibitors. Our results indicate that their higher activity is due to their polypharmacological and pleiotropic profile. These compounds can thus potentially give a clinical advantage in the management of SARS-CoV-2 infection in patients already treated with this class of drugs.

"Molecular Anatomy": a new multi-dimensional hierarchical scaffold analysis tool.

The scaffold representation is widely employed to classify bioactive compounds on the basis of common core structures or correlate compound classes with specific biological activities. In this paper, we present a novel approach called "Molecular Anatomy" as a flexible and unbiased molecular scaffold-based metrics to cluster large set of compounds. We introduce a set of nine molecular representations at different abstraction levels, combined with fragmentation rules, to define a multi-dimensional network of hierarchically interconnected molecular frameworks. We demonstrate that the introduction of a flexible scaffold definition and multiple pruning rules is an effective method to identify relevant chemical moieties. This approach allows to cluster together active molecules belonging to different molecular classes, capturing most of the structure activity information, in particular when libraries containing a huge number of singletons are analyzed. We also propose a procedure to derive a network visualization that allows a full graphical representation of compounds dataset, permitting an efficient navigation in the scaffold's space and significantly contributing to perform high quality SAR analysis. The protocol is freely available as a web interface at https://ma.exscalate.eu .

ELR+ CXC chemokines and their receptors (CXC chemokine receptor 1 and CXC chemokine receptor 2) as new therapeutic targets.

ELR+ CXC chemokines, by direct interaction with their cell surface receptors CXC chemokine receptor 1 (CXCR1) and CXC chemokine receptor 2 (CXCR2), are believed to be crucially involved in the direct migration and activation of leukocytes. ELR+ CXC chemokines are supposed to play a key role in several inflammatory diseases and this makes ELR+ CXC chemokines and their receptors attractive therapeutic targets. The first aim of this review is to discuss the potential pathological role of ELR+ CXC chemokines in different pathologies, including ulcerative colitis (UC), ischaemia/reperfusion injury (RI), bronchiolitis obliterans syndrome (BOS) and tumor progression. Moreover, the most recently described inhibitors of ELR+ CXC chemokines and their therapeutic indications will be reviewed. Finally, the mode of action and the potential therapeutical use of reparixin, a new potent and selective inhibitor of CXCR1/2 activity, and its chemical derivatives are also discussed.

Design of noncompetitive interleukin-8 inhibitors acting on CXCR1 and CXCR2.

Chemokines CXCL8 and CXCL1 play a key role in the recruitment of neutrophils at the site of inflammation. CXCL8 binds two membrane receptors, CXCR1 and CXCR2, whereas CXCL1 is a selective agonist for CXCR2. In the past decade, the physiopathological role of CXCL8 and CXCL1 has been investigated. A novel class of small molecular weight allosteric CXCR1 inhibitors was identified, and reparixin, the first drug candidate, is currently under clinical investigation in the prevention of ischemia/reperfusion injury in organ transplantation. Reparixin binding mode to CXCR1 has been studied and used for a computer-assisted design program of dual allosteric CXCR1 and CXCR2 inhibitors. In this paper, the results of modeling-driven SAR studies for the identification of potent dual inhibitors are discussed, and three new compounds (56, 67, and 79) sharing a common triflate moiety have been selected as potential leads with optimized pharmacokinetic characteristics.

PARP1-produced poly-ADP-ribose causes the PARP12 translocation to stress granules and impairment of Golgi complex functions.

Poly-ADP-ribose-polymerases (PARPs) 1 and 2 are nuclear enzymes that catalyze the poly-ADP-ribosylation of nuclear proteins transferring poly-ADP-ribose (PAR) polymers to specific residues. PARPs and PAR intervene in diverse functions, including DNA repair in the nucleus and stress granule assembly in the cytoplasm. Stress granules contribute to the regulation of translation by clustering and stabilizing mRNAs as well as several cytosolic PARPs and signaling proteins to modulate cell metabolism and survival. Our study is focused on one of these PARPs, PARP12, a Golgi-localized mono-ADP-ribosyltransferase that under stress challenge reversibly translocates from the Golgi complex to stress granules. PARP1 activation and release of nuclear PAR drive this translocation by direct PAR binding to the PARP12-WWE domain. Thus, PAR formation functionally links the activity of the nuclear and cytosolic PARPs during stress response, determining the release of PARP12 from the Golgi complex and the disassembly of the Golgi membranes, followed by a block in anterograde-membrane traffic. Notably, these functions can be rescued by reverting the stress condition (by drug wash-out). Altogether these data point at a novel, reversible nuclear signaling that senses stress to then act on cytosolic PARP12, which in turn converts the stress response into a reversible block in intracellular-membrane traffic.