The topology of drug-target interaction networks: implicit dependence on drug properties and target families

Testing the predictive power of reverse screening to infer drug targets, with the help of machine learning

Estimating protein targets of compounds based on the similarity principle-similar molecules are likely to show comparable bioactivity-is a long-standing strategy in drug research. Having previously quantified this principle, we present here a large-scale evaluation of its predictive power for inferring macromolecular targets by reverse screening an unprecedented vast external test set of more than 300,000 active small molecules against another bioactivity set of more than 500,000 compounds. We show that machine-learning can predict the correct targets, with the highest probability among 2069 proteins, for more than 51% of the external molecules. The strong enrichment thus obtained demonstrates its usefulness in supporting phenotypic screens, polypharmacology, or repurposing. Moreover, we quantified the impact of the bioactivity knowledge available for proteins in terms of number and diversity of actives. Finally, we advise that developers of such approaches follow an application-oriented benchmarking strategy and use large, high-quality, non-overlapping datasets as provided here.

Is Target-Based Drug Discovery Efficient? Discovery and "Off-Target" Mechanisms of All Drugs

Target-based drug discovery is the dominant paradigm of drug discovery; however, a comprehensive evaluation of its real-world efficiency is lacking. Here, a manual systematic review of about 32000 articles and patents dating back to 150 years ago demonstrates its apparent inefficiency. Analyzing the origins of all approved drugs reveals that, despite several decades of dominance, only 9.4% of small-molecule drugs have been discovered through "target-based" assays. Moreover, the therapeutic effects of even this minimal share cannot be solely attributed and reduced to their purported targets, as they depend on numerous off-target mechanisms unconsciously incorporated by phenotypic observations. The data suggest that reductionist target-based drug discovery may be a cause of the productivity crisis in drug discovery. An evidence-based approach to enhance efficiency seems to be prioritizing, in selecting and optimizing molecules, higher-level phenotypic observations that are closer to the sought-after therapeutic effects using tools like artificial intelligence and machine learning.

Are fibrinaloid microclots a cause of autoimmunity in Long Covid and other post-infection diseases?

It is now well established that the blood-clotting protein fibrinogen can polymerise into an anomalous form of fibrin that is amyloid in character; the resultant clots and microclots entrap many other molecules, stain with fluorogenic amyloid stains, are rather resistant to fibrinolysis, can block up microcapillaries, are implicated in a variety of diseases including Long COVID, and have been referred to as fibrinaloids. A necessary corollary of this anomalous polymerisation is the generation of novel epitopes in proteins that would normally be seen as 'self', and otherwise immunologically silent. The precise conformation of the resulting fibrinaloid clots (that, as with prions and classical amyloid proteins, can adopt multiple, stable conformations) must depend on the existing small molecules and metal ions that the fibrinogen may (and is some cases is known to) have bound before polymerisation. Any such novel epitopes, however, are likely to lead to the generation of autoantibodies. A convergent phenomenology, including distinct conformations and seeding of the anomalous form for initiation and propagation, is emerging to link knowledge in prions, prionoids, amyloids and now fibrinaloids. We here summarise the evidence for the above reasoning, which has substantial implications for our understanding of the genesis of autoimmunity (and the possible prevention thereof) based on the primary process of fibrinaloid formation.

Activity Models of Key GPCR Families in the Central Nervous System: A Tool for Many Purposes

G protein-coupled receptors (GPCRs) are targets of many drugs, of which ∼25% are indicated for central nervous system (CNS) disorders. Drug promiscuity affects their efficacy and safety profiles. Predicting the polypharmacology profile of compounds against GPCRs can thus provide a basis for producing more precise therapeutics by considering the targets and the anti-targets in that family of closely related proteins. We provide a tool for predicting the polypharmacology of compounds within prominent GPCR families in the CNS: serotonin, dopamine, histamine, muscarinic, opioid, and cannabinoid receptors. Our in-house algorithm, "iterative stochastic elimination" (ISE), produces high-quality ligand-based models for agonism and antagonism at 31 GPCRs. The ISE models correctly predict 68% of CNS drug-GPCR interactions, while the "similarity ensemble approach" predicts only 33%. The activity models correctly predict 56% of reported activities of DrugBank molecules for these CNS receptors. We conclude that the combination of interactions and activity profiles generated by screening through our models form the basis for subsequent designing and discovering novel therapeutics, either single, multitargeting, or repurposed.

Compounds co-targeting kinases in axon regulatory pathways promote regeneration and behavioral recovery after spinal cord injury in mice

Recovery from spinal cord injury (SCI) and other central nervous system (CNS) trauma is hampered by limits on axonal regeneration in the CNS. Regeneration is restricted by the lack of neuron-intrinsic regenerative capacity and by the repressive microenvironment confronting damaged axons. To address this challenge, we have developed a therapeutic strategy that co-targets kinases involved in both extrinsic and intrinsic regulatory pathways. Prior work identified a kinase inhibitor (RO48) with advantageous polypharmacology (co-inhibition of targets including ROCK2 and S6K1), which promoted CNS axon growth in vitro and corticospinal tract (CST) sprouting in a mouse pyramidotomy model. We now show that RO48 promotes neurite growth from sensory neurons and a variety of CNS neurons in vitro, and promotes CST sprouting and/or regeneration in multiple mouse models of spinal cord injury. Notably, these in vivo effects of RO48 were seen in several independent experimental series performed in distinct laboratories at different times. Finally, in a cervical dorsal hemisection model, RO48 not only promoted growth of CST axons beyond the lesion, but also improved behavioral recovery in the rotarod, gridwalk, and pellet retrieval tasks. Our results provide strong evidence for RO48 as an effective compound to promote axon growth and regeneration. Further, they point to strategies for increasing robustness of interventions in pre-clinical models.

High-throughput metabolomics predicts drug-target relationships for eukaryotic proteins

Chemical probes are important tools for understanding biological systems. However, because of the huge combinatorial space of targets and potential compounds, traditional chemical screens cannot be applied systematically to find probes for all possible druggable targets. Here, we demonstrate a novel concept for overcoming this challenge by leveraging high‐throughput metabolomics and overexpression to predict drug–target interactions. The metabolome profiles of yeast treated with 1,280 compounds from a chemical library were collected and compared with those of inducible yeast membrane protein overexpression strains. By matching metabolome profiles, we predicted which small molecules targeted which signaling systems and recovered known interactions. Drug–target predictions were generated across the 86 genes studied, including for difficult to study membrane proteins. A subset of those predictions were tested and validated, including the novel targeting of GPR1 signaling by ibuprofen. These results demonstrate the feasibility of predicting drug–target relationships for eukaryotic proteins using high‐throughput metabolomics.

The Transporter-Mediated Cellular Uptake and Efflux of Pharmaceutical Drugs and Biotechnology Products: How and Why Phospholipid Bilayer Transport Is Negligible in Real Biomembranes

Over the years, my colleagues and I have come to realise that the likelihood of pharmaceutical drugs being able to diffuse through whatever unhindered phospholipid bilayer may exist in intact biological membranes in vivo is vanishingly low. This is because (i) most real biomembranes are mostly protein, not lipid, (ii) unlike purely lipid bilayers that can form transient aqueous channels, the high concentrations of proteins serve to stop such activity, (iii) natural evolution long ago selected against transport methods that just let any undesirable products enter a cell, (iv) transporters have now been identified for all kinds of molecules (even water) that were once thought not to require them, (v) many experiments show a massive variation in the uptake of drugs between different cells, tissues, and organisms, that cannot be explained if lipid bilayer transport is significant or if efflux were the only differentiator, and (vi) many experiments that manipulate the expression level of individual transporters as an independent variable demonstrate their role in drug and nutrient uptake (including in cytotoxicity or adverse drug reactions). This makes such transporters valuable both as a means of targeting drugs (not least anti-infectives) to selected cells or tissues and also as drug targets. The same considerations apply to the exploitation of substrate uptake and product efflux transporters in biotechnology. We are also beginning to recognise that transporters are more promiscuous, and antiporter activity is much more widespread, than had been realised, and that such processes are adaptive (i.e., were selected by natural evolution). The purpose of the present review is to summarise the above, and to rehearse and update readers on recent developments. These developments lead us to retain and indeed to strengthen our contention that for transmembrane pharmaceutical drug transport "phospholipid bilayer transport is negligible".

Dual Covalent Inhibition of PKM and IMPDH Targets Metabolism in Cutaneous Metastatic Melanoma

Overcoming acquired drug resistance is a primary challenge in cancer treatment. Notably, more than 50% of patients with BRAFV600E cutaneous metastatic melanoma (CMM) eventually develop resistance to BRAF inhibitors. Resistant cells undergo metabolic reprogramming that profoundly influences therapeutic response and promotes tumor progression. Uncovering metabolic vulnerabilities could help suppress CMM tumor growth and overcome drug resistance. Here we identified a drug, HA344, that concomitantly targets two distinct metabolic hubs in cancer cells. HA344 inhibited the final and rate-limiting step of glycolysis through its covalent binding to the pyruvate kinase M2 (PKM2) enzyme, and it concurrently blocked the activity of inosine monophosphate dehydrogenase, the rate-limiting enzyme of de novo guanylate synthesis. As a consequence, HA344 efficiently targeted vemurafenib-sensitive and vemurafenib-resistant CMM cells and impaired CMM xenograft tumor growth in mice. In addition, HA344 acted synergistically with BRAF inhibitors on CMM cell lines in vitro. Thus, the mechanism of action of HA344 provides potential therapeutic avenues for patients with CMM and a broad range of different cancers. SIGNIFICANCE: Glycolytic and purine synthesis pathways are often deregulated in therapy-resistant tumors and can be targeted by the covalent inhibitor described in this study, suggesting its broad application for overcoming resistance in cancer.

Potential Biological Targets Prediction, ADME Profiling, & Molecular Docking studies of Novel Steroidal Products from Cunninghamella Blakesleana

BACKGROUND

New potential biological targets prediction through inverse molecular docking technique is an another smart strategy to forecast the possibility of compounds being biologically active against various target receptors.

OBJECTIVES

In this case of designed study, we screened our recently obtained novel acetylinic steroidal biotransformed products [(1) 8-β-methyl-14-α-hydroxy∆4tibolone (2) 9-α-Hydroxy∆4 tibolone (3) 8-β-methyl-11-β-hydroxy∆4tibolone (4) 6-β-hydroxy∆4tibolone, (5) 6-β-9-α-dihydroxy∆4tibolone (6) 7-β-hydroxy∆4tibolone) ] from fungi Cunninghemella Blakesleana to predict their possible biological targets and profiling of ADME properties.

METHOD

The prediction of pharmacokinetics properties membrane permeability as well as bioavailability radar properties were carried out by using Swiss target prediction, and Swiss ADME tools, respectively these metabolites were also subjected to predict the possible mechanism of action along with associated biological network pathways by using Reactome data-base.

RESULTS

All the six screened compounds possess excellent drug ability criteria, and exhibited exceptionally excellent non inhibitory potential against all five isozymes of CYP450 enzyme complex, including (CYP1A2, CYP2C19, CYP2C9, CYP2D6, and CYP3A4) respectively. All the screened compounds are lying within the acceptable pink zone of bioavailability radar and showing excellent descriptive properties. Compounds [1-4 & 6] are showing high BBB (Blood Brain Barrier) permeation, while compound 5 is exhibiting high HIA (Human Intestinal Absorption) property of (Egan Egg).

CONCLUSION

In conclusion, the results of this study smartly reveals that in-silico based studies are considered to provide robustness towards a rational drug designing and development approach, therefore in this way it helps to avoid the possibility of failure of drug candidates in the later experimental stages of drug development phases.

Toward an Understanding of Pan-Assay Interference Compounds and Promiscuity: A Structural Perspective on Binding Modes

Pan-assay interference compounds (PAINS) are promiscuous compound classes that produce false positive hits in high-throughput screenings. Yet, the mechanisms of PAINS activity are poorly understood. Although PAINS are often associated with protein reactivity, several recent studies have shown that they also mediate noncovalent interactions. Aiming at a deep understanding of PAINS promiscuity, we performed an analysis of the Protein Data Bank to characterize the binding modes of PAINS. We explored the binding mode conservation of 34 PAINS classes present in 871 ligands and among 517 protein targets. The two major findings of this work are the following: First, different PAINS classes exhibit different levels of binding mode conservation. Our novel classification of PAINS based on binding mode similarity enables a rational assessment of PAINS from a structural perspective. Second, PAINS classes with variable binding modes can bind with high affinity. The evaluation of noncovalent binding modes of PAINS-like compounds sheds light on the mechanisms of promiscuous binding. Our findings could facilitate the decisions on how to deal with PAINS and help scientists to understand why PAINS produce hits in their screenings.

A protet-based, protonic charge transfer model of energy coupling in oxidative and photosynthetic phosphorylation

Textbooks of biochemistry will explain that the otherwise endergonic reactions of ATP synthesis can be driven by the exergonic reactions of respiratory electron transport, and that these two half-reactions are catalyzed by protein complexes embedded in the same, closed membrane. These views are correct. The textbooks also state that, according to the chemiosmotic coupling hypothesis, a (or the) kinetically and thermodynamically competent intermediate linking the two half-reactions is the electrochemical difference of protons that is in equilibrium with that between the two bulk phases that the coupling membrane serves to separate. This gradient consists of a membrane potential term Δψ and a pH gradient term ΔpH, and is known colloquially as the protonmotive force or pmf. Artificial imposition of a pmf can drive phosphorylation, but only if the pmf exceeds some 150-170mV; to achieve in vivo rates the imposed pmf must reach 200mV. The key question then is 'does the pmf generated by electron transport exceed 200mV, or even 170mV?' The possibly surprising answer, from a great many kinds of experiment and sources of evidence, including direct measurements with microelectrodes, indicates it that it does not. Observable pH changes driven by electron transport are real, and they control various processes; however, compensating ion movements restrict the Δψ component to low values. A protet-based model, that I outline here, can account for all the necessary observations, including all of those inconsistent with chemiosmotic coupling, and provides for a variety of testable hypotheses by which it might be refined.

Deep learning and generative methods in cheminformatics and chemical biology: navigating small molecule space intelligently

The number of 'small' molecules that may be of interest to chemical biologists - chemical space - is enormous, but the fraction that have ever been made is tiny. Most strategies are discriminative, i.e. have involved 'forward' problems (have molecule, establish properties). However, we normally wish to solve the much harder generative or inverse problem (describe desired properties, find molecule). 'Deep' (machine) learning based on large-scale neural networks underpins technologies such as computer vision, natural language processing, driverless cars, and world-leading performance in games such as Go; it can also be applied to the solution of inverse problems in chemical biology. In particular, recent developments in deep learning admit the in silico generation of candidate molecular structures and the prediction of their properties, thereby allowing one to navigate (bio)chemical space intelligently. These methods are revolutionary but require an understanding of both (bio)chemistry and computer science to be exploited to best advantage. We give a high-level (non-mathematical) background to the deep learning revolution, and set out the crucial issue for chemical biology and informatics as a two-way mapping from the discrete nature of individual molecules to the continuous but high-dimensional latent representation that may best reflect chemical space. A variety of architectures can do this; we focus on a particular type known as variational autoencoders. We then provide some examples of recent successes of these kinds of approach, and a look towards the future.

The biology of ergothioneine, an antioxidant nutraceutical

Ergothioneine (ERG) is an unusual thio-histidine betaine amino acid that has potent antioxidant activities. It is synthesised by a variety of microbes, especially fungi (including in mushroom fruiting bodies) and actinobacteria, but is not synthesised by plants and animals who acquire it via the soil and their diet, respectively. Animals have evolved a highly selective transporter for it, known as solute carrier family 22, member 4 (SLC22A4) in humans, signifying its importance, and ERG may even have the status of a vitamin. ERG accumulates differentially in various tissues, according to their expression of SLC22A4, favouring those such as erythrocytes that may be subject to oxidative stress. Mushroom or ERG consumption seems to provide significant prevention against oxidative stress in a large variety of systems. ERG seems to have strong cytoprotective status, and its concentration is lowered in a number of chronic inflammatory diseases. It has been passed as safe by regulatory agencies, and may have value as a nutraceutical and antioxidant more generally.

In silico molecular target prediction unveils mebendazole as a potent MAPK14 inhibitor

The concept of polypharmacology involves the interaction of drug molecules with multiple molecular targets. It provides a unique opportunity for the repurposing of already-approved drugs to target key factors involved in human diseases. Herein, we used an in silico target prediction algorithm to investigate the mechanism of action of mebendazole, an antihelminthic drug, currently repurposed in the treatment of brain tumors. First, we confirmed that mebendazole decreased the viability of glioblastoma cells in vitro (IC₅₀ values ranging from 288 nm to 2.1 µm). Our in silico approach unveiled 21 putative molecular targets for mebendazole, including 12 proteins significantly upregulated at the gene level in glioblastoma as compared to normal brain tissue (fold change > 1.5; P < 0.0001). Validation experiments were performed on three major kinases involved in cancer biology: ABL1, MAPK1/ERK2, and MAPK14/p38α. Mebendazole could inhibit the activity of these kinases in vitro in a dose-dependent manner, with a high potency against MAPK14 (IC₅₀  = 104 ± 46 nm). Its direct binding to MAPK14 was further validated in vitro, and inhibition of MAPK14 kinase activity was confirmed in live glioblastoma cells. Consistent with biophysical data, molecular modeling suggested that mebendazole was able to bind to the catalytic site of MAPK14. Finally, gene silencing demonstrated that MAPK14 is involved in glioblastoma tumor spheroid growth and response to mebendazole treatment. This study thus highlighted the role of MAPK14 in the anticancer mechanism of action of mebendazole and provides further rationale for the pharmacological targeting of MAPK14 in brain tumors. It also opens new avenues for the development of novel MAPK14/p38α inhibitors to treat human diseases.

Drugs, host proteins and viral proteins: how their promiscuities shape antiviral design

The reciprocal nature of drug specificity and target specificity implies that the same is true for their respective promiscuities. Protein promiscuity has two broadly different types of footprint in drug design. The first is relaxed specificity of binding sites for substrates, inhibitors, effectors or cofactors. The second involves protein-protein interactions of regulatory processes such as signal transduction and transcription, and here protein intrinsic disorder plays an important role. Both viruses and host cells exploit intrinsic disorder for their survival, as do the design and discovery programs for antivirals. Drug action, strictly speaking, always relies upon promiscuous activity, with drug promiscuity enlarging its scope. Drug repurposing searches for additional promiscuity on the part of both the drug and the target in the host. Understanding the subtle nuances of these promiscuities is critical in the design of novel and more effective antivirals.

VAE-Sim: A Novel Molecular Similarity Measure Based on a Variational Autoencoder

Molecular similarity is an elusive but core "unsupervised" cheminformatics concept, yet different "fingerprint" encodings of molecular structures return very different similarity values, even when using the same similarity metric. Each encoding may be of value when applied to other problems with objective or target functions, implying that a priori none are "better" than the others, nor than encoding-free metrics such as maximum common substructure (MCSS). We here introduce a novel approach to molecular similarity, in the form of a variational autoencoder (VAE). This learns the joint distribution p(z|x) where z is a latent vector and x are the (same) input/output data. It takes the form of a "bowtie"-shaped artificial neural network. In the middle is a "bottleneck layer" or latent vector in which inputs are transformed into, and represented as, a vector of numbers (encoding), with a reverse process (decoding) seeking to return the SMILES string that was the input. We train a VAE on over six million druglike molecules and natural products (including over one million in the final holdout set). The VAE vector distances provide a rapid and novel metric for molecular similarity that is both easily and rapidly calculated. We describe the method and its application to a typical similarity problem in cheminformatics.

How polypharmacologic is each chemogenomics library?

Aim:

High-throughput phenotypic screens have emerged as a promising avenue for small-molecule drug discovery. The challenge faced in high-throughput phenotypic screens is target deconvolution once a small molecule hit is identified. Chemogenomics libraries have emerged as an important tool for meeting this challenge. Here, we investigate their target-specificity by deriving a ‘polypharmacology index’ for broad chemogenomics screening libraries.

Methods:

All known targets of all the compounds in each library were plotted as a histogram and fitted to a Boltzmann distribution, whose linearized slope is indicative of the overall polypharmacology of the library.

Results & conclusion:

Comparison of libraries clearly distinguished the most target-specific library, which might be assumed to be more useful for target deconvolution in a phenotypic screen.

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.

Involvement of multiple influx and efflux transporters in the accumulation of cationic fluorescent dyes by Escherichia coli

Background

It is widely believed that most xenobiotics cross biomembranes by diffusing through the phospholipid bilayer, and that the use of protein transporters is an occasional adjunct. According to an alternative view, phospholipid bilayer transport is negligible, and several different transporters may be involved in the uptake of an individual molecular type. We recognise here that the availability of gene knockout collections allows one to assess the contributions of all potential transporters, and flow cytometry based on fluorescence provides a convenient high-throughput assay for xenobiotic uptake in individual cells.

Results

We used high-throughput flow cytometry to assess the ability of individual gene knockout strains of E coli to take up two membrane-permeable, cationic fluorescent dyes, namely the carbocyanine diS-C3(5) and the DNA dye SYBR Green. Individual strains showed a large range of distributions of uptake. The range of modal steady-state uptakes for the carbocyanine between the different strains was 36-fold. Knockouts of the ATP synthase α- and β-subunits greatly inhibited uptake, implying that most uptake was ATP-driven rather than being driven by a membrane potential. Dozens of transporters changed the steady-state uptake of the dye by more than 50% with respect to that of the wild type, in either direction (increased or decreased); knockouts of known influx and efflux transporters behaved as expected, giving credence to the general strategy. Many of the knockouts with the most reduced uptake were transporter genes of unknown function (‘y-genes’). Similarly, several overexpression variants in the ‘ASKA’ collection had the anticipated, opposite effects. Similar results were obtained with SYBR Green (the range being approximately 69-fold). Although it too contains a benzothiazole motif there was negligible correlation between its uptake and that of the carbocyanine when compared across the various strains (although the membrane potential is presumably the same in each case).

Conclusions

Overall, we conclude that the uptake of these dyes may be catalysed by a great many transporters of putatively broad and presently unknown specificity, and that the very large range between the ‘lowest’ and the ‘highest’ levels of uptake, even in knockouts of just single genes, implies strongly that phospholipid bilayer transport is indeed negligible. This work also casts serious doubt upon the use of such dyes as quantitative stains for representing either bioenergetic parameters or the amount of cellular DNA in unfixed cells (in vivo). By contrast, it opens up their potential use as transporter assay substrates in high-throughput screening.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1561-0) contains supplementary material, which is available to authorized users.

Reliable Target Prediction of Bioactive Molecules Based on Chemical Similarity Without Employing Statistical Methods

The prediction of biological targets of bioactive molecules from machine-readable materials can be routinely performed by computational target prediction tools (CTPTs). However, the prediction of biological targets of bioactive molecules from non-digital materials (e.g., printed or handwritten documents) has not been possible due to the complex nature of bioactive molecules and impossibility of employing computations. Improving the target prediction accuracy is the most important challenge for computational target prediction. A minimum structure is identified for each group of neighbor molecules in the proposed method. Each group of neighbor molecules represents a distinct structural class of molecules with the same function in relation to the target. The minimum structure is employed as a query to search for molecules that perfectly satisfy the minimum structure of what is guessed crucial for the targeted activity. The proposed method is based on chemical similarity, but only molecules that perfectly satisfy the minimum structure are considered. Structurally related bioactive molecules found with the same minimum structure were considered as neighbor molecules of the query molecule. The known target of the neighbor molecule is used as a reference for predicting the target of the neighbor molecule with an unknown target. A lot of information is needed to identify the minimum structure, because it is necessary to know which part(s) of the bioactive molecule determines the precise target or targets responsible for the observed phenotype. Therefore, the predicted target based on the minimum structure without employing the statistical significance is considered as a reliable prediction. Since only molecules that perfectly (and not partly) satisfy the minimum structure are considered, the minimum structure can be used without similarity calculations in non-digital materials and with similarity calculations (perfect similarity) in machine-readable materials. Nine tools (PASS online, PPB, SEA, TargetHunter, PharmMapper, ChemProt, HitPick, SuperPred, and SPiDER), which can be used for computational target prediction, are compared with the proposed method for 550 target predictions. The proposed method, SEA, PPB, and PASS online, showed the best quality and quantity for the accurate predictions.

High Impact: The Role of Promiscuous Binding Sites in Polypharmacology

The literature focuses on drug promiscuity, which is a drug’s ability to bind to several targets, because it plays an essential role in polypharmacology. However, little work has been completed regarding binding site promiscuity, even though its properties are now recognized among the key factors that impact drug promiscuity. Here, we quantified and characterized the promiscuity of druggable binding sites from protein-ligand complexes in the high quality Mother Of All Databases while using statistical methods. Most of the sites (80%) exhibited promiscuity, irrespective of the protein class. Nearly half were highly promiscuous and able to interact with various types of ligands. The corresponding pockets were rather large and hydrophobic, with high sulfur atom and aliphatic residue frequencies, but few side chain atoms. Consequently, their interacting ligands can be large, rigid, and weakly hydrophilic. The selective sites that interacted with one ligand type presented less favorable pocket properties for establishing ligand contacts. Thus, their ligands were highly adaptable, small, and hydrophilic. In the dataset, the promiscuity of the site rather than the drug mainly explains the multiple interactions between the drug and target, as most ligand types are dedicated to one site. This underlines the essential contribution of binding site promiscuity to drug promiscuity between different protein classes.

The human endogenous metabolome as a pharmacology baseline for drug discovery

We have limited understanding of the variation in in vitro affinities of drugs for their targets. An analysis of a highly curated set of 815 interactions between 566 drugs and 129 primary targets reveals that 71% of drug-target affinities have values above that of the corresponding endogenous ligand, 96% of them fitting within a range of two orders of magnitude. Our findings suggest that the evolutionary optimised affinity of endogenous ligands for their native proteins can serve as a baseline for the primary pharmacology of drugs. We show that the degree of off-target selectivity and safety risks of drugs derived from their secondary pharmacology depend very much on that baseline. Thus, we propose a new approach for estimating safety margins.

The Convergence of Stem Cell Technologies and Phenotypic Drug Discovery

Recent advances in induced pluripotent stem cell technologies and phenotypic screening shape the future of bioactive small-molecule discovery. In this review we analyze the impact of small-molecule phenotypic screens on drug discovery as well as on the investigation of human development and disease biology. We further examine the role of 3D spheroid/organoid structures, microfluidic systems, and miniaturized on-a-chip systems for future discovery strategies. In highlighting representative examples, we analyze how recent achievements can translate into future therapies. Finally, we discuss remaining challenges that need to be overcome for the adaptation of the next generation of screening approaches.

Information Loss in Network Pharmacology

With the advent of increasing computational power and large-scale data acquisition, network analysis has become an attractive tool to study the organisation of complex systems and the interrelation of their constituent entities in various scientific domains. In many cases, relations only occur between entities of two different subsets, thereby forming a bipartite network. Often, the analysis of such bipartite networks involves the consideration of its two monopartite projections in order to focus on each entity subset individually as a means to deduce properties of the underlying original network. Although it is broadly acknowledged that this type of projection is not lossless, the inherent limitations of their interpretability are rarely discussed. In this work, we introduce two approaches for measuring the information loss associated with bipartite network projection. Application to two structurally distinct cases in network pharmacology, namely, drug-target and disease-gene bipartite networks, confirms that the major determinant of information loss is the degree of vertices omitted during the monopartite projection.

Covalent binders in drug discovery

Covalent modulation of protein function can have multiple utilities including therapeutics, and probes to interrogate biology. While this field is still viewed with scepticism due to the potential for (idiosyncratic) toxicities, significant strides have been made in terms of understanding how to tune electrophilicity to selectively target specific residues. Progress has also been made in harnessing the potential of covalent binders to uncover novel biology and to provide an enhanced utility as payloads for Antibody Drug Conjugates. This perspective covers the tenets and applications of covalent binders.

Extensive Reliability Evaluation of Docking-Based Target-Fishing Strategies

The development of target-fishing approaches, aimed at identifying the possible protein targets of a small molecule, represents a hot topic in medicinal chemistry. A successful target-fishing approach would allow for the elucidation of the mechanism of action of all therapeutically interesting compounds for which the actual target is still unknown. Moreover, target-fishing would be essential for preventing adverse effects of drug candidates, by predicting their potential off-targets, and it would speed up drug repurposing campaigns. However, due to the huge number of possible protein targets that a small-molecule might interact with, experimental target-fishing approaches are out of reach. In silico target-fishing represents a valuable alternative, and examples of receptor-based approaches, exploiting the large number of crystallographic protein structures determined to date, have been reported in the literature. To the best of our knowledge, no proper evaluation of such approaches is, however, reported yet. In the present work, we extensively assessed the reliability of docking-based target-fishing strategies. For this purpose, a set of X-ray structures belonging to different targets was selected, and a dataset of compounds, including 10 experimentally active ligands for each target, was created. A target-fishing benchmark database was then obtained, and used to assess the performance of 13 different docking procedures, in identifying the correct target of the dataset ligands. Moreover, a consensus docking-based target-fishing strategy was developed and evaluated. The analysis highlighted that specific features of the target proteins could affect the reliability of the protocol, which however, proved to represent a valuable tool in the proper applicability domain. Our study represents the first extensive performance assessment of docking-based target-fishing approaches, paving the way for the development of novel efficient receptor-based target fishing strategies.

Web-Based Tools for Polypharmacology Prediction

Drug promiscuity or polypharmacology is the ability of small molecules to interact with multiple protein targets simultaneously. In drug discovery, understanding the polypharmacology of potential drug molecules is crucial to improve their efficacy and safety, and to discover the new therapeutic potentials of existing drugs. Over the past decade, several computational methods have been developed to study the polypharmacology of small molecules, many of which are available as Web services. In this chapter, we review some of these Web tools focusing on ligand based approaches. We highlight in particular our recently developed polypharmacology browser (PPB) and its application for finding the side targets of a new inhibitor of the TRPV6 calcium channel.

In Silico Target Prediction for Small Molecules

Drugs modulate disease states through their actions on targets in the body. Determining these targets aids the focused development of new treatments, and helps to better characterize those already employed. One means of accomplishing this is through the deployment of in silico methodologies, harnessing computational analytical and predictive power to produce educated hypotheses for experimental verification. Here, we provide an overview of the current state of the art, describe some of the well-established methods in detail, and reflect on how they, and emerging technologies promoting the incorporation of complex and heterogeneous data-sets, can be employed to improve our understanding of (poly)pharmacology.

Steroids-specific target library for steroids target prediction

Steroids exist universally and play critical roles in various biological processes. Identifying potential targets of steroids is of great significance in studying their physiological and biochemical activities, the side effects and for drug repurposing. Herein, aiming at more precise steroids targets prediction, a steroids-specific target library integrating 3325 PDB or homology modeling structures categorized into 196 proteins was built by considering chemical similarity from DrugBank and biological processes from KEGG. The main properties of this library include: (1) It was manually prepared and checked to eliminate mistakes. (2) The library enriched the possible steroids targets and could decrease the false positives of structure-based target screening for steroids. (3) The ranking by protein name instead of PDB ID could make the screening more efficiency and precise. (4) Protein flexibility was taken into account partially by the different active conformations through the structural redundancy of each category of protein, which leads to more accurate prediction. The case studies of glycocholic acid and 24-epibrassinolide proved its powerful predictive accuracy. In summary, our strategy to build the steroids-specific protein library for steroids target prediction is a promising approach and it provides a novel idea for the target prediction of small molecules.

A systematic and prospectively validated approach for identifying synergistic drug combinations against malaria

Background

Nearly half of the world’s population (3.2 billion people) were at risk of malaria in 2015, and resistance to current therapies is a major concern. While the standard of care includes drug combinations, there is a pressing need to identify new combinations that can bypass current resistance mechanisms. In the work presented here, a combined transcriptional drug repositioning/discovery and machine learning approach is proposed.

Methods

The integrated approach utilizes gene expression data from patient-derived samples, in combination with large-scale anti-malarial combination screening data, to predict synergistic compound combinations for three Plasmodium falciparum strains (3D7, DD2 and HB3). Both single compounds and combinations predicted to be active were prospectively tested in experiment.

Results

One of the predicted single agents, apicidin, was active with the AC50 values of 74.9, 84.1 and 74.9 nM in 3D7, DD2 and HB3 P. falciparum strains while its maximal safe plasma concentration in human is 547.6 ± 136.6 nM. Apicidin at the safe dose of 500 nM kills on average 97% of the parasite. The synergy prediction algorithm exhibited overall precision and recall of 83.5 and 65.1% for mild-to-strong, 48.8 and 75.5% for moderate-to-strong and 12.0 and 62.7% for strong synergies. Some of the prospectively predicted combinations, such as tacrolimus-hydroxyzine and raloxifene-thioridazine, exhibited significant synergy across the three P. falciparum strains included in the study.

Conclusions

Systematic approaches can play an important role in accelerating discovering novel combinational therapies for malaria as it enables selecting novel synergistic compound pairs in a more informed and cost-effective manner.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2294-5) contains supplementary material, which is available to authorized users.

Drug Target Commons 2.0: a community platform for systematic analysis of drug-target interaction profiles

Drug Target Commons (DTC) is a web platform (database with user interface) for community-driven bioactivity data integration and standardization for comprehensive mapping, reuse and analysis of compound-target interaction profiles. End users can search, upload, edit, annotate and export expert-curated bioactivity data for further analysis, using an application programmable interface, database dump or tab-delimited text download options. To guide chemical biology and drug-repurposing applications, DTC version 2.0 includes updated clinical development information for the compounds and target gene-disease associations, as well as cancer-type indications for mutant protein targets, which are critical for precision oncology developments.

Structure-Based Kinase Profiling To Understand the Polypharmacological Behavior of Therapeutic Molecules

Several drugs elicit their therapeutic efficacy by modulating multiple cellular targets and possess varied polypharmacological actions. The identification of the molecular targets of a potent bioactive molecule is essential in determining its overall polypharmacological profile. Experimental procedures are expensive and time-consuming. Therefore, computational approaches are actively implemented in rational drug discovery. Here, we demonstrate a computational pipeline, based on reverse virtual screening technique using several consensus scoring strategies, and perform structure-based kinase profiling of 12 FDA-approved drugs. This target prediction showed an overall good performance, with an average AU-ROC greater than 0.85 for most drugs, and identified the true targets even at the top 2% cutoff. In contrast, 10 non-kinase binder drugs exhibited lower binding efficiency and appeared in the bottom of ranking list. Subsequently, we validated this pipeline on a potent therapeutic molecule, mahanine, whose polypharmacological profile related to targeting kinases is unknown. Our target-prediction method identified different kinases. Furthermore, we have experimentally validated that mahanine is able to modulate multiple kinases that are involved in cross-talk with different signaling molecules, which thereby exhibits its polypharmacological action. More importantly, in vitro kinase assay exhibited the inhibitory effect of mahanine on two such predicted kinases' (mTOR and VEGFR2) activity, with IC₅₀ values being ∼12 and ∼22 μM, respectively. Next, we generated a comprehensive drug-protein interaction fingerprint that explained the basis of their target selectivity. We observed that it is controlled by variations in kinase conformations followed by significant differences in crucial hydrogen-bond and van der Waals interactions. Such structure-based kinase profiling could provide useful information in revealing the unknown targets of therapeutic molecules from their polypharmacological behavior and would assist in drug discovery.

Rare Diseases: Drug Discovery and Informatics Resource

A rare disease refers to any disease with very low prevalence individually. Although the impacted population is small for a single disease, more than 6000 rare diseases affect millions of people across the world. Due to the small market size, high cost and possibly low return on investment, only in recent years, the research and development of rare disease drugs have gradually risen globally, in several domains including gene therapy, enzyme replacement therapy, and drug repositioning. Due to the complex etiology and heterogeneous symptoms, there is a large gap between basic research and patient unmet needs for rare disease drug discovery. As computational biology increasingly arises researchers' awareness, the informatics database on rare disease have grown rapidly in the recent years, including drug targets, genetic variant and mutation, phenotype and ontology and patient registries. Along with the advances of informatics database and networks, new computational models will help accelerate the target identification and lead optimization process for rare disease pre-clinical drug development.

How Phenotypic Screening Influenced Drug Discovery: Lessons from Five Years of Practice

Since 2011, phenotypic screening has been a trend in the pharmaceutical industry as well as in academia. This renaissance was triggered by analyses that suggested that phenotypic screening is a superior strategy to discover first-in-class drugs. Despite these promises and considerable investments, pharmaceutical research organizations have encountered considerable challenges with the approach. Few success stories have emerged in the past 5 years and companies are questioning their investment in this area. In this contribution, we outline what we have learned about success factors and challenges of phenotypic screening. We then describe how our efforts in phenotypic screening have influenced our approach to drug discovery in general. We predict that concepts from phenotypic screening will be incorporated into target-based approaches and will thus remain influential beyond the current trend.

Network-based Approaches in Pharmacology

In drug discovery, network-based approaches are expected to spotlight our understanding of drug action across multiple layers of information. On one hand, network pharmacology considers the drug response in the context of a cellular or phenotypic network. On the other hand, a chemical-based network is a promising alternative for characterizing the chemical space. Both can provide complementary support for the development of rational drug design and better knowledge of the mechanisms underlying the multiple actions of drugs. Recent progress in both concepts is discussed here. In addition, a network-based approach using drug-target-therapy data is introduced as an example.

Covalent inhibitors design and discovery

In the history of therapeutics, covalent drugs occupy a very distinct category. While representing a significant fraction of the drugs on the market, very few have been deliberately designed to interact covalently with their biological target. In this review, the prevalence of covalent drugs will first be briefly covered, followed by an introduction to their mechanisms of action and more detailed discussions of their discovery and the development of safe and efficient covalent enzyme inhibitors. All stages of a drug discovery program will be covered, from target considerations to lead optimization, strategies to tune reactivity and computational methods. The goal of this article is to provide an overview of the field and to outline good practices that are needed for the proper assessment and development of covalent inhibitors as well as a good understanding of the potential and limitations of current computational methods for the design of covalent drugs.

Systematic Data Mining Reveals Synergistic H3R/MCHR1 Ligands

In this study, we report a ligand-centric data mining approach that guided the identification of suitable target profiles for treating obesity. The newly developed method is based on identifying target pairs for synergistic positive effects and also encompasses the exclusion of compounds showing a detrimental effect on obesity treatment (off-targets). Ligands with known activity against obesity-relevant targets were compared using fingerprint representations. Similar compounds with activities to different targets were evaluated for the mechanism of action since activation or deactivation of drug targets determines the pharmacological effect. In vitro validation of the modeling results revealed that three known modulators of melanin-concentrating hormone receptor 1 (MCHR1) show a previously unknown submicromolar affinity to the histamine H3 receptor (H₃R). This synergistic activity may present a novel therapeutic option against obesity.

PatchSearch: A Fast Computational Method for Off-Target Detection

Many therapeutic molecules are known to bind several proteins, which can be different from the initially targeted one. Such unexpected interactions with proteins called off-targets can lead to adverse effects. Potential off-target identification is important to predict to avoid drug side effects or to discover new targets for existing drugs. We propose a new program named PatchSearch that implements local nonsequential searching for similar binding sites on protein surfaces with a controlled amount of flexibility. It is based on detection of quasi-cliques in product graphs representing all the possible matchings between two compared structures. This method has been benchmarked on a large diversity of ligands and on five data sets ranging from 12 to more than 7000 protein structures. The experiments conducted in this study show that the PatchSearch method could be useful in the early identification of off-targets. The program and the benchmarks presented in this paper are available as an R package at https://github.com/MTiPatchSearch .

Entering the 'big data' era in medicinal chemistry: molecular promiscuity analysis revisited

The 'big data' concept plays an increasingly important role in many scientific fields. Big data involves more than unprecedentedly large volumes of data that become available. Different criteria characterizing big data must be carefully considered in computational data mining, as we discuss herein focusing on medicinal chemistry. This is a scientific discipline where big data is beginning to emerge and provide new opportunities. For example, the ability of many drugs to specifically interact with multiple targets, termed promiscuity, forms the molecular basis of polypharmacology, a hot topic in drug discovery. Compound promiscuity analysis is an area that is much influenced by big data phenomena. Different results are obtained depending on chosen data selection and confidence criteria, as we also demonstrate.

An improved approach for predicting drug-target interaction: proteochemometrics to molecular docking

Proteochemometric (PCM) methods, which use descriptors of both the interacting species, i.e. drug and the target, are being successfully employed for the prediction of drug-target interactions (DTI). However, unavailability of non-interacting dataset and determining the applicability domain (AD) of model are a main concern in PCM modeling. In the present study, traditional PCM modeling was improved by devising novel methodologies for reliable negative dataset generation and fingerprint based AD analysis. In addition, various types of descriptors and classifiers were evaluated for their performance. The Random Forest and Support Vector Machine models outperformed the other classifiers (accuracies >98% and >89% for 10-fold cross validation and external validation, respectively). The type of protein descriptors had negligible effect on the developed models, encouraging the use of sequence-based descriptors over the structure-based descriptors. To establish the practical utility of built models, targets were predicted for approved anticancer drugs of natural origin. The molecular recognition interactions between the predicted drug-target pair were quantified with the help of a reverse molecular docking approach. The majority of predicted targets are known for anticancer therapy. These results thus correlate well with anticancer potential of the selected drugs. Interestingly, out of all predicted DTIs, thirty were found to be reported in the ChEMBL database, further validating the adopted methodology. The outcome of this study suggests that the proposed approach, involving use of the improved PCM methodology and molecular docking, can be successfully employed to elucidate the intricate mode of action for drug molecules as well as repositioning them for new therapeutic applications.

A dual-signal strategy for the solid detection of both small molecules and proteins based on magnetic separation and highly fluorescent copper nanoclusters

Recently, a variety of analytical methods for the detection of small molecules or proteins based on small molecule-protein interaction have been developed. However, these methods often focus on either small molecules or proteins. Few efforts are made to detect both of them in the same system. In this work, a dual-signal strategy for the solid detection of both small molecules and proteins based on small molecule-protein interaction is proposed by using the streptavidin-biotin couple as a model. In our strategy, magnetic nanoparticles (MNPs) are adopted for target separation, and highly fluorescent copper nanoclusters (CuNCs) are synthesized in situ to give signals. In the absence of the targets, CuNCs are associated with the MNPs and present in the precipitate under magnetic field; whereas in the presence of either streptavidin or biotin, the CuNCs will present in the supernate. By monitoring the fluorescent intensity of each, dual-signal can be obtained for the solid detection of either the protein or the small molecule. Results show that sensitive and specific detection of both streptavidin (detection limit: 0.47nM) and biotin (detection limit: 3.1nM) can be achieved. This method can be extended for the detection of other small molecule-protein couples, and thereby has the potential for biomedical and clinical applications.

On the translocation of bacteria and their lipopolysaccharides between blood and peripheral locations in chronic, inflammatory diseases: the central roles of LPS and LPS-induced cell death

We have recently highlighted (and added to) the considerable evidence that blood can contain dormant bacteria. By definition, such bacteria may be resuscitated (and thus proliferate). This may occur under conditions that lead to or exacerbate chronic, inflammatory diseases that are normally considered to lack a microbial component. Bacterial cell wall components, such as the endotoxin lipopolysaccharide (LPS) of Gram-negative strains, are well known as potent inflammatory agents, but should normally be cleared. Thus, their continuing production and replenishment from dormant bacterial reservoirs provides an easy explanation for the continuing, low-grade inflammation (and inflammatory cytokine production) that is characteristic of many such diseases. Although experimental conditions and determinants have varied considerably between investigators, we summarise the evidence that in a great many circumstances LPS can play a central role in all of these processes, including in particular cell death processes that permit translocation between the gut, blood and other tissues. Such localised cell death processes might also contribute strongly to the specific diseases of interest. The bacterial requirement for free iron explains the strong co-existence in these diseases of iron dysregulation, LPS production, and inflammation. Overall this analysis provides an integrative picture, with significant predictive power, that is able to link these processes via the centrality of a dormant blood microbiome that can resuscitate and shed cell wall components.

Chemical Dimerizers in Three-Hybrid Systems for Small Molecule-Target Protein Profiling

The identification of the molecular targets and mechanisms underpinning the beneficial or detrimental effects of small-molecule leads and drugs constitutes a crucial aspect of current drug discovery. Over the last two decades, three-hybrid (3H) systems have progressively taken an important position in the armamentarium of small molecule-target protein profiling technologies. Yet, a prerequisite for successful 3H analysis is the availability of appropriate chemical inducers of dimerization. Herein, we present a comprehensive and critical overview of the chemical dimerizers specifically applied in both yeast and mammalian three-hybrid systems for small molecule-target protein profiling within the broader scope of target deconvolution and drug discovery. Furthermore, examples and alternative suggestions for typical components of chemical dimerizers for 3H systems are discussed. As illustrated, more tools have become available that increase the sensitivity and efficiency of 3H-based screening platforms. Hence, it is anticipated that the great potential of 3H systems will further materialize in important contributions to drug discovery.