Off-target networks derived from ligand set similarity

Repurposing an epithelial sodium channel inhibitor as a therapy for murine and human skin inflammation

Inflammatory skin disease is characterized by a pathologic interplay between skin cells and immunocytes and can result in disfiguring cutaneous lesions and systemic inflammation. Immunosuppression is commonly used to target the inflammatory component; however, these drugs are often expensive and associated with side effects. To identify previously unidentified targets, we carried out a nonbiased informatics screen to identify drug compounds with an inverse transcriptional signature to keratinocyte inflammatory signals. Using psoriasis, a prototypic inflammatory skin disease, as a model, we used pharmacologic, transcriptomic, and proteomic characterization to find that benzamil, the benzyl derivative of the US Food and Drug Administration-approved diuretic amiloride, effectively reversed keratinocyte-driven inflammatory signaling. Through three independent mouse models of skin inflammation (Rac1G12V transgenic mice, topical imiquimod, and human skin xenografts from patients with psoriasis), we found that benzamil disrupted pathogenic interactions between the small GTPase Rac1 and its adaptor NCK1. This reduced STAT3 and NF-κB signaling and downstream cytokine production in keratinocytes. Genetic knockdown of sodium channels or pharmacological inhibition by benzamil prevented excess Rac1-NCK1 binding and limited proinflammatory signaling pathway activation in patient-derived keratinocytes without systemic immunosuppression. Both systemic and topical applications of benzamil were efficacious, suggesting that it may be a potential therapeutic avenue for treating skin inflammation.

Virtual Screening of Human Class-A GPCRs Using Ligand Profiles Built on Multiple Ligand-Receptor Interactions

G protein-coupled receptors (GPCRs) are a large family of integral membrane proteins responsible for cellular signal transductions. Identification of therapeutic compounds to regulate physiological processes is an important first step of drug discovery. We proposed MAGELLAN, a novel hierarchical virtual-screening (VS) pipeline, which starts with low-resolution protein structure prediction and structure-based binding-site identification, followed by homologous GPCR detections through structure and orthosteric binding-site comparisons. Ligand profiles constructed from the homologous ligand-GPCR complexes are then used to thread through compound databases for VS. The pipeline was first tested in a large-scale retrospective screening experiment against 224 human Class A GPCRs, where MAGELLAN achieved a median enrichment factor (EF) of 14.38, significantly higher than that using individual ligand profiles. Next, MAGELLAN was examined on 5 and 20 GPCRs from two public VS databases (DUD-E and GPCR-Bench) and resulted in an average EF of 9.75 and 13.70, respectively, which compare favorably with other state-of-the-art docking- and ligand-based methods, including AutoDock Vina (with EF = 1.48/3.16 in DUD-E and GPCR-Bench), DOCK 6 (2.12/3.47 in DUD-E and GPCR-Bench), PoLi (2.2 in DUD-E), and FINDSITECcomb2.0 (2.90 in DUD-E). Detailed data analyses show that the major advantage of MAGELLAN is attributed to the power of ligand profiling, which integrates complementary methods for ligand-GPCR interaction recognition and thus significantly improves the coverage and sensitivity of VS models. Finally, cases studies on opioid and motilin receptors show that new connections between functionally related GPCRs can be visualized in the minimum spanning tree built on the similarities of predicted ligand-binding ensembles, suggesting a novel use of MAGELLAN for GPCR deorphanization.

In Silico target fishing: addressing a "Big Data" problem by ligand-based similarity rankings with data fusion

Background

Ligand-based in silico target fishing can be used to identify the potential interacting target of bioactive ligands, which is useful for understanding the polypharmacology and safety profile of existing drugs. The underlying principle of the approach is that known bioactive ligands can be used as reference to predict the targets for a new compound.

Results

We tested a pipeline enabling large-scale target fishing and drug repositioning, based on simple fingerprint similarity rankings with data fusion. A large library containing 533 drug relevant targets with 179,807 active ligands was compiled, where each target was defined by its ligand set. For a given query molecule, its target profile is generated by similarity searching against the ligand sets assigned to each target, for which individual searches utilizing multiple reference structures are then fused into a single ranking list representing the potential target interaction profile of the query compound. The proposed approach was validated by 10-fold cross validation and two external tests using data from DrugBank and Therapeutic Target Database (TTD). The use of the approach was further demonstrated with some examples concerning the drug repositioning and drug side-effects prediction. The promising results suggest that the proposed method is useful for not only finding promiscuous drugs for their new usages, but also predicting some important toxic liabilities.

Conclusions

With the rapid increasing volume and diversity of data concerning drug related targets and their ligands, the simple ligand-based target fishing approach would play an important role in assisting future drug design and discovery.

Pharmacophore Alignment Search Tool (PhAST): Significance Assessment of Chemical Similarity

Previously, we proposed a ligand-based virtual screening technique (PhAST) based on global alignment of linearized interaction patterns. Here, we applied techniques developed for similarity assessment in local sequence alignments to our method resulting in p-values for chemical similarity. We compared two sampling strategies, a simple sampling strategy and a Markov Chain Monte Carlo (MCMC) method, and investigated the similarity of sampled distributions to Gaussian, Gumbel, modified Gumbel, and Gamma distributions. The Gumbel distribution with a Gaussian correction term was identified as the most similar to the observed empirical distributions. These techniques were applied in retrospective screenings on a drug-like dataset. Obtained p-values were adjusted to the size of the screening library with four different methods. Evaluation of E-value thresholds corroborated the Bonferroni correction as a preferred means to identify significant chemical similarity with PhAST. An online version of PhAST with significance estimation is available at http://modlab-cadd.ethz.ch/.

Finding novel pharmaceuticals in the systems biology era using multiple effective drug targets, phenotypic screening and knowledge of transporters: where drug discovery went wrong and how to fix it

Despite the sequencing of the human genome, the rate of innovative and successful drug discovery in the pharmaceutical industry has continued to decrease. Leaving aside regulatory matters, the fundamental and interlinked intellectual issues proposed to be largely responsible for this are: (a) the move from 'function-first' to 'target-first' methods of screening and drug discovery; (b) the belief that successful drugs should and do interact solely with single, individual targets, despite natural evolution's selection for biochemical networks that are robust to individual parameter changes; (c) an over-reliance on the rule-of-5 to constrain biophysical and chemical properties of drug libraries; (d) the general abandoning of natural products that do not obey the rule-of-5; (e) an incorrect belief that drugs diffuse passively into (and presumably out of) cells across the bilayers portions of membranes, according to their lipophilicity; (f) a widespread failure to recognize the overwhelmingly important role of proteinaceous transporters, as well as their expression profiles, in determining drug distribution in and between different tissues and individual patients; and (g) the general failure to use engineering principles to model biology in parallel with performing 'wet' experiments, such that 'what if?' experiments can be performed in silico to assess the likely success of any strategy. These facts/ideas are illustrated with a reasonably extensive literature review. Success in turning round drug discovery consequently requires: (a) decent systems biology models of human biochemical networks; (b) the use of these (iteratively with experiments) to model how drugs need to interact with multiple targets to have substantive effects on the phenotype; (c) the adoption of polypharmacology and/or cocktails of drugs as a desirable goal in itself; (d) the incorporation of drug transporters into systems biology models, en route to full and multiscale systems biology models that incorporate drug absorption, distribution, metabolism and excretion; (e) a return to 'function-first' or phenotypic screening; and (f) novel methods for inferring modes of action by measuring the properties on system variables at all levels of the 'omes. Such a strategy offers the opportunity of achieving a state where we can hope to predict biological processes and the effect of pharmaceutical agents upon them. Consequently, this should both lower attrition rates and raise the rates of discovery of effective drugs substantially.

Significance estimation for sequence-based chemical similarity searching (PhAST) and application to AuroraA kinase inhibitors

BACKGROUND

Chemical similarity searching allows the retrieval of preferred screening molecules from a compound database. Candidates are ranked according to their similarity to a reference compound (query). Assessing the statistical significance of chemical similarity scores helps prioritizing significant hits, and identifying cases where the database does not contain any promising compounds.

METHOD

Our text-based similarity measure, Pharmacophore Alignment Search Tool (PhAST), employs pair-wise sequence alignment. We adapted the concept of E-values as significance estimates and employed a sampling technique that incorporates the principle of importance sampling in a Markov chain Monte Carlo simulation to generate distributions of random alignment scores. These distributions were used to compute significance estimates for similarity scores in a preliminary prospective virtual screen for inhibitors of Aurora A kinase.

CONCLUSION

Assessing the significance of compound similarity computed with PhAST allows for a statistically motivated identification of candidate screening compounds. Inhibitors of Aurora A kinase were retrieved from a large compound library.

The promiscuous binding of pharmaceutical drugs and their transporter-mediated uptake into cells: what we (need to) know and how we can do so

A recent paper in this journal sought to counter evidence for the role of transport proteins in effecting drug uptake into cells, and questions that transporters can recognize drug molecules in addition to their endogenous substrates. However, there is abundant evidence that both drugs and proteins are highly promiscuous. Most proteins bind to many drugs and most drugs bind to multiple proteins (on average more than six), including transporters (mutations in these can determine resistance); most drugs are known to recognise at least one transporter. In this response, we alert readers to the relevant evidence that exists or is required. This needs to be acquired in cells that contain the relevant proteins, and we highlight an experimental system for simultaneous genome-wide assessment of carrier-mediated uptake in a eukaryotic cell (yeast).

A novel method of transcriptional response analysis to facilitate drug repositioning for cancer therapy

Little research has been done to address the huge opportunities that may exist to reposition existing approved or generic drugs for alternate uses in cancer therapy. In addition, there has been little work on strategies to reposition experimental cancer agents for testing in alternate settings that could shorten their clinical development time. Progress in each area has lagged, in part, because of the lack of systematic methods to define drug off-target effects (OTE) that might affect important cancer cell signaling pathways. In this study, we addressed this critical gap by developing an OTE-based method to repurpose drugs for cancer therapeutics, based on transcriptional responses made in cells before and after drug treatment. Specifically, we defined a new network component called cancer-signaling bridges (CSB) and integrated it with a Bayesian factor regression model (BFRM) to form a new hybrid method termed CSB-BFRM. Proof-of-concept studies were conducted in breast and prostate cancer cells and in promyelocytic leukemia cells. In each system, CSB-BFRM analysis could accurately predict clinical responses to more than 90% of drugs approved by the U.S. Food and Drug Administration and more than 75% of experimental clinical drugs that were tested. Mechanistic investigation of OTEs for several high-ranking drug-dose pairs suggested repositioning opportunities for cancer therapy, based on the ability to enforce retinoblastoma-dependent repression of important E2F-dependent cell-cycle genes. Together, our findings establish new methods to identify opportunities for drug repositioning or to elucidate the mechanisms of action of repositioned drugs.

The role of translational bioinformatics in drug discovery

The application of translational approaches (e.g. from bed to bench and back) is gaining momentum in the pharmaceutical industry. By utilizing the rapidly increasing volume of data at all phases of drug discovery, translational bioinformatics is poised to address some of the key challenges faced by the industry. Indeed, computational analysis of clinical data and patient records has informed decision-making in multiple aspects of drug discovery and development. Here, we review key examples of translational bioinformatics approaches to emphasize its potential to enhance the quality of drug discovery pipelines, reduce attrition rates and, ultimately, lead to more effective treatments.

The role of computational methods in the identification of bioactive compounds

Computational methods play an ever increasing role in lead finding. A vast repertoire of molecular design and virtual screening methods emerged in the past two decades and are today routinely used. There is increasing awareness that there is no single best computational protocol and correspondingly there is a shift recommending the combination of complementary methods. A promising trend for the application of computational methods in lead finding is to take advantage of the vast amounts of HTS (High Throughput Screening) data to allow lead assessment by detailed systems-based data analysis, especially for phenotypic screens where the identification of compound-target pairs is the primary goal. Herein, we review trends and provide examples of successful applications of computational methods in lead finding.

The chemical basis of pharmacology

Molecular biology now dominates pharmacology so thoroughly that it is difficult to recall that only a generation ago the field was very different. To understand drug action today, we characterize the targets through which they act and new drug leads are discovered on the basis of target structure and function. Until the mid-1980s the information often flowed in reverse: investigators began with organic molecules and sought targets, relating receptors not by sequence or structure but by their ligands. Recently, investigators have returned to this chemical view of biology, bringing to it systematic and quantitative methods of relating targets by their ligands. This has allowed the discovery of new targets for established drugs, suggested the bases for their side effects, and predicted the molecular targets underlying phenotypic screens. The bases for these new methods, some of their successes and liabilities, and new opportunities for their use are described.

Prediction and evaluation of protein farnesyltransferase inhibition by commercial drugs

The similarity ensemble approach (SEA) relates proteins based on the set-wise chemical similarity among their ligands. It can be used to rapidly search large compound databases and to build cross-target similarity maps. The emerging maps relate targets in ways that reveal relationships one might not recognize based on sequence or structural similarities alone. SEA has previously revealed cross talk between drugs acting primarily on G-protein coupled receptors (GPCRs). Here we used SEA to look for potential off-target inhibition of the enzyme protein farnesyltransferase (PFTase) by commercially available drugs. The inhibition of PFTase has profound consequences for oncogenesis, as well as a number of other diseases. In the present study, two commercial drugs, Loratadine and Miconazole, were identified as potential ligands for PFTase and subsequently confirmed as such experimentally. These results point toward the applicability of SEA for the prediction of not only GPCR-GPCR drug cross talk but also GPCR-enzyme and enzyme-enzyme drug cross talk.