BioGPS: Navigating biological space to predict polypharmacology, off-targeting, and selectivity

Design, synthesis, and biological evaluation of first-in-class indomethacin-based PROTACs degrading SARS-CoV-2 main protease and with broad-spectrum antiviral activity

To date, Proteolysis Targeting Chimera (PROTAC) technology has been successfully applied to mediate proteasomal-induced degradation of several pharmaceutical targets mainly related to oncology, immune disorders, and neurodegenerative diseases. On the other hand, its exploitation in the field of antiviral drug discovery is still in its infancy. Recently, we described two indomethacin (INM)-based PROTACs displaying broad-spectrum antiviral activity against coronaviruses. Here, we report the design, synthesis, and characterization of a novel series of INM-based PROTACs that recruit either Von-Hippel Lindau (VHL) or cereblon (CRBN) E3 ligases. The panel of INM-based PROTACs was also enlarged by varying the linker moiety. The antiviral activity resulted very susceptible to this modification, particularly for PROTACs hijacking VHL as E3 ligase, with one piperazine-based compound (PROTAC 6) showing potent anti-SARS-CoV-2 activity in infected human lung cells. Interestingly, degradation assays in both uninfected and virus-infected cells with the most promising PROTACs emerged so far (PROTACs 5 and 6) demonstrated that INM-PROTACs do not degrade human PGES-2 protein, as initially hypothesized, but induce the concentration-dependent degradation of SARS-CoV-2 main protease (Mpro) both in Mpro-transfected and in SARS-CoV-2-infected cells. Importantly, thanks to the target degradation, INM-PROTACs exhibited a considerable enhancement in antiviral activity with respect to indomethacin, with EC50 values in the low-micromolar/nanomolar range. Finally, kinetic solubility as well as metabolic and chemical stability were measured for PROTACs 5 and 6. Altogether, the identification of INM-based PROTACs as the first class of SARS-CoV-2 Mpro degraders demonstrating activity also in SARS-CoV-2-infected cells represents a significant advance in the development of effective, broad-spectrum anti-coronavirus strategies.

DeepGRID: Deep Learning Using GRID Descriptors for BBB Prediction

Deep Learning approaches are able to automatically extract relevant features from the input data and capture nonlinear relationships between the input and output. In this work, we present the GRID-derived AI (GrAId) descriptors, a simple modification to GRID MIFs that facilitate their use in combination with Convolutional Neural Networks (CNNs) to build Deep Learning models in a rotationally, conformationally, and alignment-independent approach we are calling DeepGRID. To our knowledge, this is the first time that GRID MIFs have been combined with CNNs in a Deep Learning approach. We applied the approach to build regression and classification models for blood-brain barrier permeation, an important factor when designing CNS drugs and conversely when designing to avoid off-target effects for CNS-inactive drugs. The VolSurf approach was one of the first to successfully model this property from three-dimensional structures, using descriptors derived from their GRID Molecular Interaction Fields (MIFs) in combination with PLS. We compared the DeepGRID models with others built using the hand-crafted VolSurf descriptors in combination with both PLS and Random Forest (RF). Both the DeepGRID and RF regression models performed best according to the % of compounds with a Geometric Mean Fold Error (GMFE) within 2-fold of the experimental data. Applying these regression models as classifiers, for the smaller 332 and 416 compound data sets all models performed well with ROC AUC values of ∼0.9 on the external test set. For the larger 2105 compound data set, the DeepGRID classifier performed the best with an AUC of 0.87 on the external test set with the RF model having an AUC of 0.84 and the original VolSurf lgBB model having an AUC of 0.83.

Growth hormone secretagogues modulate inflammation and fibrosis in mdx mouse model of Duchenne muscular dystrophy

Introduction

Growth hormone secretagogues (GHSs) exert multiple actions, being able to activate GHS-receptor 1a, control inflammation and metabolism, to enhance GH/insulin-like growth factor-1 (IGF-1)-mediated myogenesis, and to inhibit angiotensin-converting enzyme. These mechanisms are of interest for potentially targeting multiple steps of pathogenic cascade in Duchenne muscular dystrophy (DMD).

Methods

Here, we aimed to provide preclinical evidence for potential benefits of GHSs in DMD, via a multidisciplinary in vivo and ex vivo comparison in mdx mice, of two ad hoc synthesized compounds (EP80317 and JMV2894), with a wide but different profile. 4-week-old mdx mice were treated for 8 weeks with EP80317 or JMV2894 (320 µg/kg/d, s.c.).

Results

In vivo, both GHSs increased mice forelimb force (recovery score, RS towards WT: 20% for EP80317 and 32% for JMV2894 at week 8). In parallel, GHSs also reduced diaphragm (DIA) and gastrocnemius (GC) ultrasound echodensity, a fibrosis-related parameter (RS: ranging between 26% and 75%). Ex vivo, both drugs ameliorated DIA isometric force and calcium-related indices (e.g., RS: 40% for tetanic force). Histological analysis highlighted a relevant reduction of fibrosis in GC and DIA muscles of treated mice, paralleled by a decrease in gene expression of TGF-β1 and Col1a1. Also, decreased levels of pro-inflammatory genes (IL-6, CD68), accompanied by an increment in Sirt-1, PGC-1α and MEF2c expression, were observed in response to treatments, suggesting an overall improvement of myofiber metabolism. No detectable transcript levels of GHS receptor-1a, nor an increase of circulating IGF-1 were found, suggesting the presence of a novel receptor-independent mechanism in skeletal muscle. Preliminary docking studies revealed a potential binding capability of JMV2894 on metalloproteases involved in extracellular matrix remodeling and cytokine production, such as ADAMTS-5 and MMP-9, overactivated in DMD.

Discussion

Our results support the interest of GHSs as modulators of pathology progression in mdx mice, disclosing a direct anti-fibrotic action that may prove beneficial to contrast pathological remodeling.

ELIOT: A platform to navigate the E3 pocketome and aid the design of new PROTACs

Proteolysis-targeting chimeras (PROTACs) are novel therapeutics for the treatment of human disease. They exploit the enormous potential of the E3 ligases, a class of proteins that mark a target protein for degradation via the ubiquitin-proteasome system. Despite the existence of several E3 ligase-related databases, the choice of the functioning ligase is limited to only 1.6% of those available, probably due to the fragmentary understanding of their structures and their known ligands; in fact, none of the existing databases report detailed studies covering their 3D structure or their pockets. Here, we report ELIOT (E3 LIgase pocketOme navigaTor), an accurate and complete platform containing the E3 ligase pocketome to enable navigation and selection of new E3 ligases and new ligands for the design of new PROTACs. All E3 ligase pockets were characterized with innovative 3D descriptors including their PROTAC-ability score, and similarity analyses between E3 pockets are presented. Tissue specificity and their degree of involvement in patients with specific cancer types are also annotated for each E3 ligase, enabling appropriate selection for the design of a PROTAC with improved specificity. All data are available at https://eliot.moldiscovery.com.

An Integrated Machine Learning Model To Spot Peptide Binding Pockets in 3D Protein Screening

The prediction of peptide-protein binding sites is of utmost importance to tackle the onset of severe neurodegenerative diseases and cancer. In this work, we detail a novel machine learning model based on Linear Discriminant Analysis (LDA) demonstrating to be highly predictive in detecting the putative protein binding regions of small peptides. Starting from 439 high-quality pockets derived from peptide-protein crystallographic complexes, three sets of well-established peptide-binding regions were first selected through a Partitioning Around Medoids (PAM) clustering algorithm based on morphological and energetic 3D GRID-MIF molecular descriptors. Next, the best combination between all the putative interacting peptide pockets and related GRID-MIF scores was automatically explored by using the LDA-based protocol implemented in BioGPS. This approach proved successful to recognize the actual interacting peptide regions (that is, AUC = 0.86 and partial ROC enrichment at 5% of 0.48) from all the other pockets of the protein. Validated on two external collections sets, including 445 and 347 crystallographic peptide-protein complexes, our LDA-based model could be effective to further run peptide-protein virtual screening campaigns.

Discovery of a new class of triazole based inhibitors of acetyl transferase KAT2A

We have recently developed a new synthetic methodology that provided both N-aryl-5-hydroxytriazoles and N-pyridine-4-alkyl triazoles. A selection of these products was carried through virtual screening towards targets that are contemporary and validated for drug discovery and development. This study determined a number of potential structure target dyads of which N-pyridinium-4-carboxylic-5-alkyl triazole displayed the highest score specificity towards KAT2A. Binding affinity tests of abovementioned triazole and related analogs towards KAT2A confirmed the predictions of the in-silico assay. Finally, we have run in vitro inhibition assays of selected triazoles towards KAT2A; the ensemble of binding and inhibition assays delivered pyridyl-triazoles carboxylates as the prototype of a new class of inhibitors of KAT2A.

Exploiting ELIOT for Scaffold-Repurposing Opportunities: TRIM33 a Possible Novel E3 Ligase to Expand the Toolbox for PROTAC Design

The field of targeted protein degradation, through the control of the ubiquitin-proteasome system (UPS), is progressing considerably; to exploit this new therapeutic modality, the proteolysis targeting chimera (PROTAC) technology was born. The opportunity to use PROTACs engaging of new E3 ligases that can hijack and control the UPS system could greatly extend the applicability of degrading molecules. To this end, here we show a potential application of the ELIOT (E3 LIgase pocketOme navigaTor) platform, previously published by this group, for a scaffold-repurposing strategy to identify new ligands for a novel E3 ligase, such as TRIM33. Starting from ELIOT, a case study of the cross-relationship using GRID Molecular Interaction Field (MIF) similarities between TRIM24 and TRIM33 binding sites was selected. Based on the assumption that similar pockets could bind similar ligands and considering that TRIM24 has 12 known co-crystalised ligands, we applied a scaffold-repurposing strategy for the identification of TRIM33 ligands exploiting the scaffold of TRIM24 ligands. We performed a deeper computational analysis to identify pocket similarities and differences, followed by docking and water analysis; selected ligands were synthesised and subsequently tested against TRIM33 via HTRF binding assay, and we obtained the first-ever X-ray crystallographic complexes of TRIM33α with three of the selected compounds.

Exploring Ligand Binding Domain Dynamics in the NRs Superfamily

Nuclear receptors (NRs) are transcription factors that play an important role in multiple diseases, such as cancer, inflammation, and metabolic disorders. They share a common structural organization composed of five domains, of which the ligand-binding domain (LBD) can adopt different conformations in response to substrate, agonist, and antagonist binding, leading to distinct transcription effects. A key feature of NRs is, indeed, their intrinsic dynamics that make them a challenging target in drug discovery. This work aims to provide a meaningful investigation of NR structural variability to outline a dynamic profile for each of them. To do that, we propose a methodology based on the computation and comparison of protein cavities among the crystallographic structures of NR LBDs. First, pockets were detected with the FLAPsite algorithm and then an "all against all" approach was applied by comparing each pair of pockets within the same sub-family on the basis of their similarity score. The analysis concerned all the detectable cavities in NRs, with particular attention paid to the active site pockets. This approach can guide the investigation of NR intrinsic dynamics, the selection of reference structures to be used in drug design and the easy identification of alternative binding sites.

CROMATIC: Cross-Relationship Map of Cavities from Coronaviruses

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiologic agent of COVID-19 disease, has rapidly imposed an urgent need to identify effective drug candidates. In this context, the high resolution and non-redundant beta-Coronavirus protein cavities database is pivotal to help virtual screening protocols. Furthermore, the cross-relationship among cavities can lead to highlighting multitarget therapy chances. Here, we first collect all protein cavities on SARS-CoV-2, SARS-CoV, and MERS-CoV X-ray structures, and then, we compute a similarity map by using molecular interaction fields (MIFs). All the results come together in CROMATIC (CROss-relationship MAp of CaviTIes from Coronaviruses). CROMATIC encloses both a comprehensive and a non-redundant version of the cavities collection and a similarity map revealing, on the one hand, cavities that are conserved among the three Coronaviruses and, on the other hand, unexpected similarities among cavities that can represent a key starting point for multitarget therapy strategies. Similarity analysis was also performed for the available structures of SARS-CoV-2 spike variants, linking sequence mutations to three-dimensional interaction alterations. The CROMATIC repository is freely available to the scientific community at https://github.com/moldiscovery/sars-cromatic.

Discovery of novel SARS-CoV-2 inhibitors targeting the main protease Mpro by virtual screenings and hit optimization

Two years after its emergence, SARS-CoV-2 still represents a serious and global threat to human health. Antiviral drug development usually takes a long time and, to increase the chances of success, chemical variability of hit compounds represents a valuable source for the discovery of new antivirals. In this work, we applied a platform of variably oriented virtual screening campaigns to seek for novel chemical scaffolds for SARS-CoV-2 main protease (M^pro) inhibitors. The study on the resulting 30 best hits led to the identification of a series of structurally unrelated M^pro inhibitors. Some of them exhibited antiviral activity in the low micromolar range against SARS-CoV-2 and other human coronaviruses (HCoVs) in different cell lines. Time-of-addition experiments demonstrated an antiviral effect during the viral replication cycle at a time frame consistent with the inhibition of SARS-CoV-2 M^pro activity. As a proof-of-concept, to validate the pharmaceutical potential of the selected hits against SARS-CoV-2, we rationally optimized one of the hit compounds and obtained two potent SARS-CoV-2 inhibitors with increased activity against M^pro both in vitro and in a cellular context, as well as against SARS-CoV-2 replication in infected cells. This study significantly contributes to the expansion of the chemical variability of SARS-CoV-2 M^pro inhibitors and provides new scaffolds to be exploited for pan-coronavirus antiviral drug development.

PLATO: A Predictive Drug Discovery Web Platform for Efficient Target Fishing and Bioactivity Profiling of Small Molecules

PLATO (Polypharmacology pLATform predictiOn) is an easy-to-use drug discovery web platform, which has been designed with a two-fold objective: to fish putative protein drug targets and to compute bioactivity values of small molecules. Predictions are based on the similarity principle, through a reverse ligand-based screening, based on a collection of 632,119 compounds known to be experimentally active on 6004 protein targets. An efficient backend implementation allows to speed-up the process that returns results for query in less than 20 s. The graphical user interface is intuitive to give practitioners easy input and transparent output, which is available as a standard report in portable document format. PLATO has been validated on thousands of external data, with performances better than those of other parallel approaches. PLATO is available free of charge (http://plato.uniba.it/ accessed on 13 April 2022).

Getting Insights into Structural and Energetic Properties of Reciprocal Peptide-Protein Interactions

Peptide-protein interactions play a key role for many cellular and metabolic processes involved in the onset of largely spread diseases such as cancer and neurodegenerative pathologies. Despite the progress in the structural characterization of peptide-protein interfaces, the in-depth knowledge of the molecular details behind their interactions is still a daunting task. Here, we present the first comprehensive in silico morphological and energetic study of peptide binding sites by focusing on both peptide and protein standpoints. Starting from the PixelDB database, a nonredundant benchmark collection of high-quality 3D crystallographic structures of peptide-protein complexes, a classification analysis of the most representative categories based on the nature of each cocrystallized peptide has been carried out. Several interpretable geometrical and energetic descriptors have been computed both from peptide and target protein sides in the attempt to unveil physicochemical and structural causative correlations. Finally, we investigated the most frequent peptide-protein residue pairs at the binding interface and made extensive energetic analyses, based on GRID MIFs, with the aim to study the peptide affinity-enhancing interactions to be further exploited in rational drug design strategies.

FragExplorer: GRID-Based Fragment Growing and Replacement

Understanding which chemical modifications can be made to known ligands is a key aspect of structure-based drug design and one that was pioneered by the software GRID. We developed FragExplorer with the explicit aim of showing GRID users which fragments would best match the GRID molecular interaction fields in a protein binding site, given a bound ligand as a starting point. Users can grow ligands or replace existing moieties; the R-Group Exploration mode identifies all potential R-Groups and searches for replacements automatically; the Scaffold Exploration mode does the same for all potential scaffolds. For a ligand with three points of variation, R-Group Exploration will typically explore a chemical space of 10¹⁶ potential molecules; including Scaffold Exploration increases this to 10²². FragExplorer was designed to be integrated within an interactive 3D Editor/Designer; therefore, the speed of computation was an important consideration; a typical fragment search takes 20 seconds. In a fragment reprediction test, FragExplorer demonstrates an overall fragment retrieval rate of 55%, increasing to 69% for smaller fragments. At a 90% substructural match, the retrieval rate increases to ∼80%. We also show how the approach could have been used to hop from olmesartan to azilsartan or to optimize a p38 MAP kinase lead to a compound that bears similarity to a known nanomolar inhibitor.

Polypharmacology: The science of multi-targeting molecules

Polypharmacology is a concept where a molecule can interact with two or more targets simultaneously. It offers many advantages as compared to the conventional single-targeting molecules. A multi-targeting drug is much more efficacious due to its cumulative efficacy at all of its individual targets making it much more effective in complex and multifactorial diseases like cancer, where multiple proteins and pathways are involved in the onset and development of the disease. For a molecule to be polypharmacologic in nature, it needs to possess promiscuity which is the ability to interact with multiple targets; and at the same time avoid binding to antitargets which would otherwise result in off-target adverse effects. There are certain structural features and physicochemical properties which when present would help researchers to predict if the designed molecule would possess promiscuity or not. Promiscuity can also be identified via advanced state-of-the-art computational methods. In this review, we also elaborate on the methods by which one can intentionally incorporate promiscuity in their molecules and make them polypharmacologic. The polypharmacology paradigm of "one drug-multiple targets" has numerous applications especially in drug repurposing where an already established drug is redeveloped for a new indication. Though designing a polypharmacological drug is much more difficult than designing a single-targeting drug, with the current technologies and information regarding different diseases and chemical functional groups, it is plausible for researchers to intentionally design a polypharmacological drug and unlock its advantages.

Connecting chemistry and biology through molecular descriptors

Through the representation of small molecule structures as numerical descriptors and the exploitation of the similarity principle, chemoinformatics has made paramount contributions to drug discovery, from unveiling mechanisms of action and repurposing approved drugs to de novo crafting of molecules with desired properties and tailored targets. Yet, the inherent complexity of biological systems has fostered the implementation of large-scale experimental screenings seeking a deeper understanding of the targeted proteins, the disrupted biological processes and the systemic responses of cells to chemical perturbations. After this wealth of data, a new generation of data-driven descriptors has arisen providing a rich portrait of small molecule characteristics that goes beyond chemical properties. Here, we give an overview of biologically relevant descriptors, covering chemical compounds, proteins and other biological entities, such as diseases and cell lines, while aligning them to the major contributions in the field from disciplines, such as natural language processing or computer vision. We now envision a new scenario for chemical and biological entities where they both are translated into a common numerical format. In this computational framework, complex connections between entities can be unveiled by means of simple arithmetic operations, such as distance measures, additions, and subtractions.

First-in-Class Isonipecotamide-Based Thrombin and Cholinesterase Dual Inhibitors with Potential for Alzheimer Disease

Recently, the direct thrombin (thr) inhibitor dabigatran has proven to be beneficial in animal models of Alzheimer’s disease (AD). Aiming at discovering novel multimodal agents addressing thr and AD-related targets, a selection of previously and newly synthesized potent thr and factor Xa (fXa) inhibitors were virtually screened by the Multi-fingerprint Similarity Searching aLgorithm (MuSSeL) web server. The N-phenyl-1-(pyridin-4-yl)piperidine-4-carboxamide derivative 1, which has already been experimentally shown to inhibit thr with a K_i value of 6 nM, has been flagged by a new, upcoming release of MuSSeL as a binder of cholinesterase (ChE) isoforms (acetyl- and butyrylcholinesterase, AChE and BChE), as well as thr, fXa, and other enzymes and receptors. Interestingly, the inhibition potency of 1 was predicted by the MuSSeL platform to fall within the low-to-submicromolar range and this was confirmed by experimental K_i values, which were found equal to 0.058 and 6.95 μM for eeAChE and eqBChE, respectively. Thirty analogs of 1 were then assayed as inhibitors of thr, fXa, AChE, and BChE to increase our knowledge of their structure-activity relationships, while the molecular determinants responsible for the multiple activities towards the target enzymes were rationally investigated by molecular cross-docking screening.

Chemical Modulation of the 1-(Piperidin-4-yl)-1,3-dihydro-2H-benzo[d]imidazole-2-one Scaffold as a Novel NLRP3 Inhibitor

In the search for new chemical scaffolds able to afford NLRP3 inflammasome inhibitors, we used a pharmacophore-hybridization strategy by combining the structure of the acrylic acid derivative INF39 with the 1-(piperidin-4-yl)1,3-dihydro-2H-benzo[d]imidazole-2-one substructure present in HS203873, a recently identified NLRP3 binder. A series of differently modulated benzo[d]imidazole-2-one derivatives were designed and synthesised. The obtained compounds were screened in vitro to test their ability to inhibit NLRP3-dependent pyroptosis and IL-1β release in PMA-differentiated THP-1 cells stimulated with LPS/ATP. The selected compounds were evaluated for their ability to reduce the ATPase activity of human recombinant NLRP3 using a newly developed assay. From this screening, compounds 9, 13 and 18, able to concentration-dependently inhibit IL-1β release in LPS/ATP-stimulated human macrophages, emerged as the most promising NLRP3 inhibitors of the series. Computational simulations were applied for building the first complete model of the NLRP3 inactive state and for identifying possible binding sites available to the tested compounds. The analyses led us to suggest a mechanism of protein–ligand binding that might explain the activity of the compounds.

Inverse Virtual Screening for the rapid re-evaluation of the presumed biological safe profile of natural products. The case of steviol from Stevia rebaudiana glycosides on farnesoid X receptor (FXR)

Nonnutritive sweeteners (NNSs) are widely employed as dietary substitutes for classical sugars thanks to their safety profile and low toxicity. In this study, a re-evaluation of the biological effects of steviol (1), the main metabolite from Stevia rebaudiana glycosides, was performed using the Inverse Virtual Screening (IVS) target fishing computational approach. Starting from well-known pharmacological properties of Stevia rebaudiana glycosides, this computational tool was employed for predicting the putative interacting targets of 1 and, afterwards, of its five synthetic ester derivatives 2-6, accounting a large panel of proteins involved in cancer and inflammation events. Applying this methodology, the farnesoid X receptor (FXR) was identified as the putative target partner of 1-6. The predicted ligand-protein interactions were corroborated by transactivation assays, specifically disclosing the agonistic activity of 1 and the antagonistic activities of 2-6 on FXR. The reported results highlight the feasibility of IVS as a fast and potent tool for predicting the interacting targets of query compounds, addressing the re-evaluation of their bioactivity. In light of the obtained results, the presumably safe profile of known compounds, such as the case of steviol (1), is critically discussed.

Away from Flatness: Unprecedented Nitrogen-Bridged Cyclopenta[a]indene Derivatives as Novel Anti-Alzheimer Multitarget Agents

Nature-inspired, bridged polycyclic molecules share low similarity with currently available drugs, containing preferentially planar and/or achiral moieties. This "Escape from Flatland" scenario, aimed at exploring pharmacological properties of atypical molecular scaffolds, finds interest in synthetic routes leading to tridimensional-shaped molecules. Herein we report on the synthesis of N-bridged cyclopenta[a]indene derivatives, achieved through microwave-assisted thermal rearrangement of allene 3-benzazecines with high diastereoselectivity. The biological evaluation disclosed selective inhibition of human acetylcholinesterase or butyrylcholinesterase, depending on the substitution around the molecular core, which was rationalized by means of docking simulations. The most potent BChE inhibitor 31 was effective in neuroprotection from glutamatergic excitotoxicity and displayed low intrinsic cytotoxicity and good brain penetration. Overall, compound 31 and its close congeners 34 and 35 acted as multitarget agents addressing different biological events involved in neurodegeneration, particularly in the progression of Alzheimer's disease.

Targeting Nsp9 as an anti-SARS-CoV-2 strategy

Non-structural protein 9 (Nsp9) plays a key role in viral replication of coronavirus and represents a promising target for anti-SARS-CoV-2 strategies.

A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs

We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.

Bcr-Abl Allosteric Inhibitors: Where We Are and Where We Are Going to

The fusion oncoprotein Bcr-Abl is an aberrant tyrosine kinase responsible for chronic myeloid leukemia and acute lymphoblastic leukemia. The auto-inhibition regulatory module observed in the progenitor kinase c-Abl is lost in the aberrant Bcr-Abl, because of the lack of the N-myristoylated cap able to bind the myristoyl binding pocket also conserved in the Bcr-Abl kinase domain. A way to overcome the occurrence of resistance phenomena frequently observed for Bcr-Abl orthosteric drugs is the rational design of allosteric ligands approaching the so-called myristoyl binding pocket. The discovery of these allosteric inhibitors although very difficult and extremely challenging, represents a valuable option to minimize drug resistance, mostly due to the occurrence of mutations more frequently affecting orthosteric pockets, and to enhance target selectivity with lower off-target effects. In this perspective, we will elucidate at a molecular level the structural bases behind the Bcr-Abl allosteric control and will show how artificial intelligence can be effective to drive the automated de novo design towards off-patent regions of the chemical space.

Artificial Neural Networks in Computer-Aided Drug Design: An Overview of Recent Advances

Computer-aided drug design (CADD) is the framework in which the huge amount of data accumulated by high-throughput experimental methods used in drug design is quantitatively studied. Its objectives include pattern recognition, biomarker identification and/or classification, etc. In order to achieve these objectives, machine learning algorithms and especially artificial neural networks (ANNs) have been used over ADMET factor testing and QSAR modeling evaluation. This paper provides an overview of the current trends in CADD-applied ANNs, since their use was re-boosted over a decade ago.

An up-to-date overview of computational polypharmacology in modern drug discovery

INTRODUCTION

In recent years, computational polypharmacology has gained significant attention to study the promiscuous nature of drugs. Despite tremendous challenges, community-wide efforts have led to a variety of novel approaches for predicting drug polypharmacology. In particular, some rapid advances using machine learning and artificial intelligence have been reported with great success.

AREAS COVERED

In this article, the authors provide a comprehensive update on the current state-of-the-art polypharmacology approaches and their applications, focusing on those reports published after our 2017 review article. The authors particularly discuss some novel, groundbreaking concepts, and methods that have been developed recently and applied to drug polypharmacology studies.

EXPERT OPINION

Polypharmacology is evolving and novel concepts are being introduced to counter the current challenges in the field. However, major hurdles remain including incompleteness of high-quality experimental data, lack of in vitro and in vivo assays to characterize multi-targeting agents, shortage of robust computational methods, and challenges to identify the best target combinations and design effective multi-targeting agents. Fortunately, numerous national/international efforts including multi-omics and artificial intelligence initiatives as well as most recent collaborations on addressing the COVID-19 pandemic have shown significant promise to propel the field of polypharmacology forward.

Protein electrostatics: From computational and structural analysis to discovery of functional fingerprints and biotechnological design

Computationally driven engineering of proteins aims to allow them to withstand an extended range of conditions and to mediate modified or novel functions. Therefore, it is crucial to the biotechnological industry, to biomedicine and to afford new challenges in environmental sciences, such as biocatalysis for green chemistry and bioremediation. In order to achieve these goals, it is important to clarify molecular mechanisms underlying proteins stability and modulating their interactions. So far, much attention has been given to hydrophobic and polar packing interactions and stability of the protein core. In contrast, the role of electrostatics and, in particular, of surface interactions has received less attention. However, electrostatics plays a pivotal role along the whole life cycle of a protein, since early folding steps to maturation, and it is involved in the regulation of protein localization and interactions with other cellular or artificial molecules. Short- and long-range electrostatic interactions, together with other forces, provide essential guidance cues in molecular and macromolecular assembly. We report here on methods for computing protein electrostatics and for individual or comparative analysis able to sort proteins by electrostatic similarity. Then, we provide examples of electrostatic analysis and fingerprints in natural protein evolution and in biotechnological design, in fields as diverse as biocatalysis, antibody and nanobody engineering, drug design and delivery, molecular virology, nanotechnology and regenerative medicine.

Multiple catalytic activities of human 17β-hydroxysteroid dehydrogenase type 7 respond differently to inhibitors

Cholesterol biosynthesis is a multistep process in mammals that includes the aerobic removal of three methyl groups from the intermediate lanosterol, one from position 14 and two from position 4. During the demethylations at position 4, a 3-ketosteroid reductase catalyses the conversion of both 4-methylzymosterone and zymosterone to 4-methylzymosterol and zymosterol, respectively, restoring the alcoholic function of lanosterol, which is also maintained in cholesterol. Unlike other eukaryotes, mammals also use the same enzyme as an estrone reductase that can transform estrone (E1) into estradiol (E2). This enzyme, named 17β-hydroxysteroid dehydrogenase type 7 (HSD17B7), is therefore a multifunctional protein in mammals, and one that belongs to both the HSD17B family, which is involved in steroid-hormone metabolism, and to the family of post-squalene cholesterol biosynthesis enzymes. In the present study, a series of known inhibitors of human HSD17B7's E1-reductase activity have been assayed for potential inhibition against 3-ketosteroid reductase activity. Surprisingly, the assayed compounds lost their inhibition activity when tested in HepG2 cells that were incubated with radiolabelled acetate and against the recombinant overexpressed human enzyme incubated with 4-methylzymosterone (both radiolabelled and not). Preliminary kinetic analyses suggest a mixed or non-competitive inhibition on the E1-reductase activity, which is in agreement with Molecular Dynamics simulations. These results raise questions about the mechanism(s) of action of these possible inhibitors, the enzyme dynamic regulation and the interplay between the two activities.

Acute pancreatitis promotes the generation of two different exosome populations

Exosomes are small extracellular vesicles that act as intercellular messengers. Previous studies revealed that, during acute pancreatitis, circulating exosomes could reach the alveolar compartment and activate macrophages. However, proteomic analysis suggested that the most likely origin of these exosomes could be the liver instead of the pancreas. The present study aimed to characterize the exosomes released by pancreas to pancreatitis-associated ascitic fluid (PAAF) as well as those circulating in plasma in an experimental model of taurocholate-induced acute pancreatitis in rats. We provide evidence that during acute pancreatitis two different populations of exosomes are generated with relevant differences in cell distribution, protein and microRNA content as well as different implications in their physiological effects. During pancreatitis plasma exosomes, but not PAAF exosomes, are enriched in the inflammatory miR-155 and show low levels of miR-21 and miR-122. Mass spectrometry-based proteomic analysis showed that PAAF exosomes contains 10–30 fold higher loading of histones and ribosomal proteins compared to plasma exosomes. Finally, plasma exosomes have higher pro-inflammatory activity on macrophages than PAAF exosomes. These results confirm the generation of two different populations of exosomes during acute pancreatitis. Deep understanding of their specific functions will be necessary to use them as therapeutic targets at different stages of the disease.

A platform for target prediction of phenotypic screening hit molecules

Many drug discovery programmes, particularly for infectious diseases, are conducted phenotypically. Identifying the targets of phenotypic screening hits experimentally can be complex, time-consuming, and expensive. However, it would be valuable to know what the molecular target(s) is, as knowledge of the binding pose of the hit molecule in the binding site can facilitate the compound optimisation. Furthermore, knowing the target would allow de-prioritisation of less attractive chemical series or molecular targets. To generate target-hypotheses for phenotypic active compounds, an in silico platform was developed that utilises both ligand and protein-structure information to generate a ranked set of predicted molecular targets. As a result of the web-based workflow the user obtains a set of 3D structures of the predicted targets with the active molecule bound. The platform was exemplified using Mycobacterium tuberculosis, the causative organism of tuberculosis. In a test that we performed, the platform was able to predict the targets of 60% of compounds investigated, where there was some similarity to a ligand in the protein database.

•

An algorithm to predict the molecular target(s) of phenotypic hits against TB.

•

Uses information based on the ligand and protein structure.

•

Allow visualisation of proposed binding pose.

•

Web interface developed.

Methods for Discovering and Targeting Druggable Protein-Protein Interfaces and Their Application to Repurposing

Drug repurposing is a creative and resourceful approach to increase the number of therapies by exploiting available and approved drugs. However, identifying new protein targets for previously approved drugs is challenging. Although new strategies have been developed for drug repurposing, there is broad agreement that there is room for further improvements. In this chapter, we review protein-protein interaction (PPI) interface-targeting strategies for drug repurposing applications. We discuss certain features, such as hot spot residue and hot region prediction and their importance in drug repurposing, and illustrate common methods used in PPI networks to identify drug off-targets. We also collect available online resources for hot spot prediction, binding pocket identification, and interface clustering which are effective resources in polypharmacology. Finally, we provide case studies showing the significance of protein interfaces and hot spots in drug repurposing.

Selection of protein conformations for structure-based polypharmacology studies

Several drugs exert their therapeutic effect through the modulation of multiple targets. Structure-based approaches hold great promise for identifying compounds with the desired polypharmacological profiles. These methods use knowledge of the protein binding sites to identify stereoelectronically complementary ligands. The selection of the most suitable protein conformations to be used in the design process is vital, especially for multitarget drug design in which the same ligand has to be accommodated in multiple binding pockets. Herein, we focus on currently available techniques for the selection of the most suitable protein conformations for multitarget drug design, compare the potential advantages and limitations of each method, and comment on how their combination could help in polypharmacology drug design.

Truly Target-Focused Pharmacophore Modeling: A Novel Tool for Mapping Intermolecular Surfaces

Pharmacophore models are an accurate and minimal tridimensional abstraction of intermolecular interactions between chemical structures, usually derived from a group of molecules or from a ligand-target complex. Only a limited amount of solutions exists to model comprehensive pharmacophores using the information of a particular target structure without knowledge of any binding ligand. In this work, an automated and customable tool for truly target-focused (T²F) pharmacophore modeling is introduced. Key molecular interaction fields of a macromolecular structure are calculated using the AutoGRID energy functions. The most relevant points are selected by a newly developed filtering cascade and clustered to pharmacophore features with a density-based algorithm. Using five different protein classes, the ability of this method to identify essential pharmacophore features was compared to structure-based pharmacophores derived from ligand-target interactions. This method represents an extremely valuable instrument for drug design in a situation of scarce ligand information available, but also in the case of underexplored therapeutic targets, as well as to investigate protein allosteric pockets and protein-protein interactions.

Electrostatic recognition in substrate binding to serine proteases

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2‐step mechanism of protein‐protein recognition suggested in the literature.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Identifying DNase I hypersensitive sites as driver distal regulatory elements in breast cancer

Efforts to identify driver mutations in cancer have largely focused on genes, whereas non-coding sequences remain relatively unexplored. Here we develop a statistical method based on characteristics known to influence local mutation rate and a series of enrichment filters in order to identify distal regulatory elements harboring putative driver mutations in breast cancer. We identify ten DNase I hypersensitive sites that are significantly mutated in breast cancers and associated with the aberrant expression of neighboring genes. A pan-cancer analysis shows that three of these elements are significantly mutated across multiple cancer types and have mutation densities similar to protein-coding driver genes. Functional characterization of the most highly mutated DNase I hypersensitive sites in breast cancer (using in silico and experimental approaches) confirms that they are regulatory elements and affect the expression of cancer genes. Our study suggests that mutations of regulatory elements in tumors likely play an important role in cancer development.

Cancer driver mutations can occur within noncoding genomic sequences. Here, the authors develop a statistical approach to identify candidate noncoding driver mutations in DNase I hypersensitive sites in breast cancer and experimentally demonstrate they are regulatory elements of known cancer genes.

Detecting similar binding pockets to enable systems polypharmacology

In the era of systems biology, multi-target pharmacological strategies hold promise for tackling disease-related networks. In this regard, drug promiscuity may be leveraged to interfere with multiple receptors: the so-called polypharmacology of drugs can be anticipated by analyzing the similarity of binding sites across the proteome. Here, we perform a pairwise comparison of 90,000 putative binding pockets detected in 3,700 proteins, and find that 23,000 pairs of proteins have at least one similar cavity that could, in principle, accommodate similar ligands. By inspecting these pairs, we demonstrate how the detection of similar binding sites expands the space of opportunities for the rational design of drug polypharmacology. Finally, we illustrate how to leverage these opportunities in protein-protein interaction networks related to several therapeutic classes and tumor types, and in a genome-scale metabolic model of leukemia.

Traditionally, the fact that most drugs are promiscuous binders has been a major concern in pharmacology, due to the occurrence of undesired off-target clinical events. In the recent years, however, the realization that many diseases are the result of complex biological processes has encouraged rethinking of drug promiscuity as a promising feature, since it is sometimes necessary to interfere with multiple receptors in order to overcome the robustness of disease-related networks. One way to identify groups of proteins that could be targeted simultaneously is to look for similar binding sites. We have massively done so for all human proteins with a known high-resolution three-dimensional structure, unveiling a vast space of ‘polypharmacology’ opportunities. Of these, we know, a great majority is not of therapeutic interest. To pinpoint promising multi-target combinations, we advocate for the use of computational tools that are able to rapidly simulate the effect of drug-target interactions on biological networks.

Involvement of exosomes in lung inflammation associated with experimental acute pancreatitis

A frequent complication of acute pancreatitis is the lung damage associated with the systemic inflammatory response. Although various pro-inflammatory mediators generated at both local and systemic levels have been identified, the pathogenic mechanisms of the disease are still poorly understood. In recent years, exosomes have emerged as a new intercellular communication system able to transfer encapsulated proteins and small RNAs and protect them from degradation. Using an experimental model of taurocholate-induced acute pancreatitis in rats, we aimed to evaluate the role of exosomes in the extent of the systemic inflammatory response. Induction of pancreatitis increased the concentration of circulating exosomes, which showed a different proteomic profile to those obtained from control animals. A series of tracking experiments using PKH26-stained exosomes revealed that circulating exosomes effectively reached the alveolar compartment and were internalized by macrophages. In vitro experiments revealed that exosomes obtained under inflammatory conditions activate and polarize these alveolar macrophages towards a pro-inflammatory phenotype. Interestingly, the proteomic analysis of circulating exosomes during acute pancreatitis suggested a multi-organ origin with a relevant role for the liver as a source of these vesicles. Tracking experiments also revealed that the liver retains the majority of exosomes from the peritoneal cavity. We conclude that exosomes are involved in the lung damage associated with experimental acute pancreatitis and could be relevant mediators in the systemic effects of pancreatitis. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Computational polypharmacology: a new paradigm for drug discovery

INTRODUCTION

Over the past couple of years, the cost of drug development has sharply increased along with the high rate of clinical trial failures. Such increase in expenses is partially due to the inability of the "one drug - one target" approach to predict drug side effects and toxicities. To tackle this issue, an alternative approach, known as polypharmacology, is being adopted to study small molecule interactions with multiple targets. Apart from developing more potent and effective drugs, this approach allows for studies of off-target activities and the facilitation of drug repositioning. Although exhaustive polypharmacology studies in-vitro or in-vivo are not practical, computational methods of predicting unknown targets or side effects are being developed. Areas covered: This article describes various computational approaches that have been developed to study polypharmacology profiles of small molecules. It also provides a brief description of the algorithms used in these state-of-the-art methods. Expert opinion: Recent success in computational prediction of multi-targeting drugs has established polypharmacology as a promising alternative approach to tackle some of the daunting complications in drug discovery. This will not only help discover more effective agents, but also present tremendous opportunities to study novel target pharmacology and facilitate drug repositioning efforts in the pharmaceutical industry.

Improved Potency of Indole-Based NorA Efflux Pump Inhibitors: From Serendipity toward Rational Design and Development

The NorA efflux pump is a potential drug target for reversal of resistance to selected antibacterial agents, and recently we described indole-based inhibitor candidates. Herein we report a second class of inhibitors derived from them but with significant differences in shape and size. In particular, compounds 13 and 14 are very potent inhibitors in that they demonstrated the lowest IC₅₀ values (2 μM) ever observed among all indole-based compounds we have evaluated.

Identification of KasA as the cellular target of an anti-tubercular scaffold

Phenotypic screens for bactericidal compounds are starting to yield promising hits against tuberculosis. In this regard, whole-genome sequencing of spontaneous resistant mutants generated against an indazole sulfonamide (GSK3011724A) identifies several specific single-nucleotide polymorphisms in the essential Mycobacterium tuberculosis β-ketoacyl synthase (kas) A gene. Here, this genomic-based target assignment is confirmed by biochemical assays, chemical proteomics and structural resolution of a KasA-GSK3011724A complex by X-ray crystallography. Finally, M. tuberculosis GSK3011724A-resistant mutants increase the in vitro minimum inhibitory concentration and the in vivo 99% effective dose in mice, establishing in vitro and in vivo target engagement. Surprisingly, the lack of target engagement of the related β-ketoacyl synthases (FabH and KasB) suggests a different mode of inhibition when compared with other Kas inhibitors of fatty acid biosynthesis in bacteria. These results clearly identify KasA as the biological target of GSK3011724A and validate this enzyme for further drug discovery efforts against tuberculosis.

Comparing Drug Images and Repurposing Drugs with BioGPS and FLAPdock: The Thymidylate Synthase Case

Repurposing and repositioning drugs has become a frequently pursued and successful strategy in the current era, as new chemical entities are increasingly difficult to find and get approved. Herein we report an integrated BioGPS/FLAPdock pipeline for rapid and effective off-target identification and drug repurposing. Our method is based on the structural and chemical properties of protein binding sites, that is, the ligand image, encoded in the GRID molecular interaction fields (MIFs). Protein similarity is disclosed through the BioGPS algorithm by measuring the pockets' overlap according to which pockets are clustered. Co-crystallized and known ligands can be cross-docked among similar targets, selected for subsequent in vitro binding experiments, and possibly improved for inhibitory potency. We used human thymidylate synthase (TS) as a test case and searched the entire RCSB Protein Data Bank (PDB) for similar target pockets. We chose casein kinase IIα as a control and tested a series of its inhibitors against the TS template. Ellagic acid and apigenin were identified as TS inhibitors, and various flavonoids were selected and synthesized in a second-round selection. The compounds were demonstrated to be active in the low-micromolar range.

Research Resources for Nuclear Receptor Signaling Pathways

Nuclear receptor (NR) signaling pathways impact cellular function in a broad variety of tissues in both normal physiology and disease states. The complex tissue-specific biology of these pathways is an enduring impediment to the development of clinical NR small-molecule modulators that combine therapeutically desirable effects in specific target tissues with suppression of off-target effects in other tissues. Supporting the important primary research in this area is a variety of web-based resources that assist researchers in gaining an appreciation of the molecular determinants of the pharmacology of a NR pathway in a given tissue. In this study, selected representative examples of these tools are reviewed, along with discussions on how current and future generations of tools might optimally adapt to the future of NR signaling research.

Identification and Visualization of Kinase-Specific Subpockets

The identification and design of selective compounds is important for the reduction of unwanted side effects as well as for the development of tool compounds for target validation studies. This is, in particular, true for therapeutically important protein families that possess conserved folds and have numerous members such as kinases. To support the design of selective kinase inhibitors, we developed a novel approach that allows identification of specificity determining subpockets between closely related kinases solely based on their three-dimensional structures. To account for the intrinsic flexibility of the proteins, multiple X-ray structures of the target protein of interest as well as of unwanted off-target(s) are taken into account. The binding pockets of these protein structures are calculated and fused to a combined target and off-target pocket, respectively. Subsequently, shape differences between these two combined pockets are identified via fusion rules. The approach provides a user-friendly visualization of target-specific areas in a binding pocket which should be explored when designing selective compounds. Furthermore, the approach can be easily combined with in silico alanine mutation studies to identify selectivity determining residues. The potential impact of the approach is demonstrated in four retrospective experiments on closely related kinases, i.e., p38α vs Erk2, PAK1 vs PAK4, ITK vs AurA, and BRAF vs VEGFR2. Overall, the presented approach does not require any profiling data for training purposes, provides an intuitive visualization of a large number of protein structures at once, and could also be applied to other target classes.

In silico methods to address polypharmacology: current status, applications and future perspectives

Polypharmacology, a new paradigm in drug discovery that focuses on multi-target drugs (MTDs), has potential application for drug repurposing, the process of finding new uses for existing approved drugs, prediction of off-target toxicities and rational design of MTDs. In this scenario, computational strategies have demonstrated great potential in predicting polypharmacology and in facilitating drug repurposing. Here, we provide a comprehensive overview of various computational approaches that enable the prediction and analysis of in vitro and in vivo drug-response phenotypes and outline their potential for drug discovery. We give an outlook on the latest advances in rational design of MTDs and discuss possible future directions of algorithm development in this field.