Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: The Multi-Targeting Drug DREAM Challenge

A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets (‘polypharmacology’). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.     

A network-based multi-target computational estimation scheme for anticoagulant activities of compounds

Background

Traditional virtual screening method pays more attention on predicted binding affinity between drug molecule and target related to a certain disease instead of phenotypic data of drug molecule against disease system, as is often less effective on discovery of the drug which is used to treat many types of complex diseases. Virtual screening against a complex disease by general network estimation has become feasible with the development of network biology and system biology. More effective methods of computational estimation for the whole efficacy of a compound in a complex disease system are needed, given the distinct weightiness of the different target in a biological process and the standpoint that partial inhibition of several targets can be more efficient than the complete inhibition of a single target.

Methodology

We developed a novel approach by integrating the affinity predictions from multi-target docking studies with biological network efficiency analysis to estimate the anticoagulant activities of compounds. From results of network efficiency calculation for human clotting cascade, factor Xa and thrombin were identified as the two most fragile enzymes, while the catalytic reaction mediated by complex IXa:VIIIa and the formation of the complex VIIIa:IXa were recognized as the two most fragile biological matter in the human clotting cascade system. Furthermore, the method which combined network efficiency with molecular docking scores was applied to estimate the anticoagulant activities of a serial of argatroban intermediates and eight natural products respectively. The better correlation (r = 0.671) between the experimental data and the decrease of the network deficiency suggests that the approach could be a promising computational systems biology tool to aid identification of anticoagulant activities of compounds in drug discovery.

Conclusions

This article proposes a network-based multi-target computational estimation method for anticoagulant activities of compounds by combining network efficiency analysis with scoring function from molecular docking.     

基于SVM 的药物靶点预测方法及其应用

目的:基于已知药物靶点和潜在药物靶点蛋白的一级结构相似性,结合SVM技术研究新的有效的药物靶点预测方法。方法:构造训练样本集,提取蛋白质序列的一级结构特征,进行数据预处理,选择最优核函数,优化参数并进行特征选择,训练最优预测模型,检验模型的预测效果。以G蛋白偶联受体家族的蛋白质为预测集,应用建立的最优分类模型对其进行潜在药物靶点挖掘。结果:基于SVM所建立的最优分类模型预测的平均准确率为81.03%。应用最优分类器对构造的G蛋白预测集进行预测,结果发现预测排位在前20的蛋白质中有多个与疾病相关。特别的,其中有两个G蛋白在治疗靶点数据库(TTD)中显示已作为临床试验的药物靶点。结论:基于SVM和蛋白质序列特征的药物靶点预测方法是有效的,应用该方法预测出的潜在药物靶点能够为发现新的药靶提供参考。     

Drug2ways: Reasoning over causal paths in biological networks for drug discovery

Elucidating the causal mechanisms responsible for disease can reveal potential therapeutic targets for pharmacological intervention and, accordingly, guide drug repositioning and discovery. In essence, the topology of a network can reveal the impact a drug candidate may have on a given biological state, leading the way for enhanced disease characterization and the design of advanced therapies. Network-based approaches, in particular, are highly suited for these purposes as they hold the capacity to identify the molecular mechanisms underlying disease. Here, we present drug2ways, a novel methodology that leverages multimodal causal networks for predicting drug candidates. Drug2ways implements an efficient algorithm which reasons over causal paths in large-scale biological networks to propose drug candidates for a given disease. We validate our approach using clinical trial information and demonstrate how drug2ways can be used for multiple applications to identify: i) single-target drug candidates, ii) candidates with polypharmacological properties that can optimize multiple targets, and iii) candidates for combination therapy. Finally, we make drug2ways available to the scientific community as a Python package that enables conducting these applications on multiple standard network formats.     

A Structure-Based Approach for Mapping Adverse Drug Reactions to the Perturbation of Underlying Biological Pathways

Adverse drug reactions (ADR), also known as side-effects, are complex undesired physiologic phenomena observed secondary to the administration of pharmaceuticals. Several phenomena underlie the emergence of each ADR; however, a dominant factor is the drug''s ability to modulate one or more biological pathways. Understanding the biological processes behind the occurrence of ADRs would lead to the development of safer and more effective drugs. At present, no method exists to discover these ADR-pathway associations. In this paper we introduce a computational framework for identifying a subset of these associations based on the assumption that drugs capable of modulating the same pathway may induce similar ADRs. Our model exploits multiple information resources. First, we utilize a publicly available dataset pairing drugs with their observed ADRs. Second, we identify putative protein targets for each drug using the protein structure database and in-silico virtual docking. Third, we label each protein target with its known involvement in one or more biological pathways. Finally, the relationships among these information sources are mined using multiple stages of logistic-regression while controlling for over-fitting and multiple-hypothesis testing. As proof-of-concept, we examined a dataset of 506 ADRs, 730 drugs, and 830 human protein targets. Our method yielded 185 ADR-pathway associations of which 45 were selected to undergo a manual literature review. We found 32 associations to be supported by the scientific literature.     

Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers

Common inflammatome gene signatures as well as disease‐specific signatures were identified by analyzing 12 expression profiling data sets derived from 9 different tissues isolated from 11 rodent inflammatory disease models. The inflammatome signature significantly overlaps with known drug targets and co‐expressed gene modules linked to metabolic disorders and cancer. A large proportion of genes in this signature are tightly connected in tissue‐specific Bayesian networks (BNs) built from multiple independent mouse and human cohorts. Both the inflammatome signature and the corresponding consensus BNs are highly enriched for immune response‐related genes supported as causal for adiposity, adipokine, diabetes, aortic lesion, bone, muscle, and cholesterol traits, suggesting the causal nature of the inflammatome for a variety of diseases. Integration of this inflammatome signature with the BNs uncovered 151 key drivers that appeared to be more biologically important than the non‐drivers in terms of their impact on disease phenotypes. The identification of this inflammatome signature, its network architecture, and key drivers not only highlights the shared etiology but also pinpoints potential targets for intervention of various common diseases.     

Target Inhibition Networks: Predicting Selective Combinations of Druggable Targets to Block Cancer Survival Pathways

A recent trend in drug development is to identify drug combinations or multi-target agents that effectively modify multiple nodes of disease-associated networks. Such polypharmacological effects may reduce the risk of emerging drug resistance by means of attacking the disease networks through synergistic and synthetic lethal interactions. However, due to the exponentially increasing number of potential drug and target combinations, systematic approaches are needed for prioritizing the most potent multi-target alternatives on a global network level. We took a functional systems pharmacology approach toward the identification of selective target combinations for specific cancer cells by combining large-scale screening data on drug treatment efficacies and drug-target binding affinities. Our model-based prediction approach, named TIMMA, takes advantage of the polypharmacological effects of drugs and infers combinatorial drug efficacies through system-level target inhibition networks. Case studies in MCF-7 and MDA-MB-231 breast cancer and BxPC-3 pancreatic cancer cells demonstrated how the target inhibition modeling allows systematic exploration of functional interactions between drugs and their targets to maximally inhibit multiple survival pathways in a given cancer type. The TIMMA prediction results were experimentally validated by means of systematic siRNA-mediated silencing of the selected targets and their pairwise combinations, showing increased ability to identify not only such druggable kinase targets that are essential for cancer survival either individually or in combination, but also synergistic interactions indicative of non-additive drug efficacies. These system-level analyses were enabled by a novel model construction method utilizing maximization and minimization rules, as well as a model selection algorithm based on sequential forward floating search. Compared with an existing computational solution, TIMMA showed both enhanced prediction accuracies in cross validation as well as significant reduction in computation times. Such cost-effective computational-experimental design strategies have the potential to greatly speed-up the drug testing efforts by prioritizing those interventions and interactions warranting further study in individual cancer cases.     

Network pharmacology: the next paradigm in drug discovery

The dominant paradigm in drug discovery is the concept of designing maximally selective ligands to act on individual drug targets. However, many effective drugs act via modulation of multiple proteins rather than single targets. Advances in systems biology are revealing a phenotypic robustness and a network structure that strongly suggests that exquisitely selective compounds, compared with multitarget drugs, may exhibit lower than desired clinical efficacy. This new appreciation of the role of polypharmacology has significant implications for tackling the two major sources of attrition in drug development--efficacy and toxicity. Integrating network biology and polypharmacology holds the promise of expanding the current opportunity space for druggable targets. However, the rational design of polypharmacology faces considerable challenges in the need for new methods to validate target combinations and optimize multiple structure-activity relationships while maintaining drug-like properties. Advances in these areas are creating the foundation of the next paradigm in drug discovery: network pharmacology.     

Approximating Optimal Controls for Networks when There Are Combinations of Population-Level and Targeted Measures Available: Chlamydia Infection as a Case-Study

Using a modified one-dimensional model for the spread of an SIS disease on a network, we show that the behaviour of complex network simulations can be replicated with a simpler model. This model is then used to design optimal controls for use on the network, which would otherwise be unfeasible to obtain, resulting in information about how best to combine a population-level random intervention with one that is more targeted. This technique is used to minimise intervention costs over a short time interval with a target prevalence, and also to minimise prevalence with a specified budget. When applied to chlamydia, we find results consistent with previous work; that is maximising targeted control (contact tracing) is important to using resources effectively, while high-intensity bursts of population control (screening) are more effective than maintaining a high level of coverage.     

Selecting targets for therapeutic validation through differential protein expression using chromatography-mass spectrometry

The identification of potential targets for therapeutic intervention can be accomplished on a systematic basis by a variety of techniques that include quantitative analysis of gene-specific mRNA levels and expressed proteins in normal and diseased cells. Differences in the expression levels of nucleic acid and protein gene products could suggest protein drug targets that are directly causative of disease, or reveal biochemical pathways that could be modulated by therapeutic molecules. Any effort based on mRNA or protein expression level comparisons could be confounded by a number of factors: level in steady-state may not be correlated with actual encoded protein levels; differentially expressed protein levels might be a result of disease process, and not causative of the process, and therapeutic intervention based on such a difference will be unproductive and the differential expression of mRNA or protein may be the result of biological variation unrelated to the disease process under study. In order to address these possibly confounding factors, it is necessary to validate potential targets by establishing their firm association with disease, and their minimal distribution in non-diseased tissues of any type. This requirement suggests that emphasis on true and reproducible quantitation of protein expression levels in a variety of samples will be an effective and highly efficient method of generating drug targets with a high degree of utility. To achieve this aim, we have established an industrial-scale proteomics-based discovery platform consisting of cell biology, protein chemistry, and mass spectrometry technical groups together with bioinformatics groups. The analytical method used for quantitation employs isotope labeling for differential analysis (ICATTM, Applied Biosystems, Inc.). With this technique, tryptic peptides are generated from labeled proteins that have been specifically captured from various subcellular locations or protein families. The resulting peptides are identified and quantified by mass spectrometry. To evaluate this approach on a large-scale, we have applied it to a study of continuous cell lines derived from human pancreatic adenocarcinomas. We have been able to establish processes for target discovery for small molecule drug targets as well as therapeutic antibody target identification for cell surface proteins. In addition, we have developed a process for identification of serum markers of this disease based upon standardized fractionation procedures. The results of these analyses will be presented together with the some of the issues from both the wet and dry (computational) lab that need to be addressed in such an undertaking.     

SMET: Systematic multiple enzyme targeting – a method to rationally design optimal strains for target chemical overproduction

Identifying multiple enzyme targets for metabolic engineering is very critical for redirecting cellular metabolism to achieve desirable phenotypes, e.g., overproduction of a target chemical. The challenge is to determine which enzymes and how much of these enzymes should be manipulated by adding, deleting, under-, and/or over-expressing associated genes. In this study, we report the development of a systematic multiple enzyme targeting method (SMET), to rationally design optimal strains for target chemical overproduction. The SMET method combines both elementary mode analysis and ensemble metabolic modeling to derive SMET metrics including l-values and c-values that can identify rate-limiting reaction steps and suggest which enzymes and how much of these enzymes to manipulate to enhance product yields, titers, and productivities. We illustrated, tested, and validated the SMET method by analyzing two networks, a simple network for concept demonstration and an Escherichia coli metabolic network for aromatic amino acid overproduction. The SMET method could systematically predict simultaneous multiple enzyme targets and their optimized expression levels, consistent with experimental data from the literature, without performing an iterative sequence of single-enzyme perturbation. The SMET method was much more efficient and effective than single-enzyme perturbation in terms of computation time and finding improved solutions.     

Intervention in gene regulatory networks via greedy control policies based on long-run behavior

Background  

A salient purpose for studying gene regulatory networks is to derive intervention strategies, the goals being to identify potential drug targets and design gene-based therapeutic intervention. Optimal stochastic control based on the transition probability matrix of the underlying Markov chain has been studied extensively for probabilistic Boolean networks. Optimization is based on minimization of a cost function and a key goal of control is to reduce the steady-state probability mass of undesirable network states. Owing to computational complexity, it is difficult to apply optimal control for large networks.     

A network pharmacology‐based approach to explore the effects of Chaihu Shugan powder on a non‐alcoholic fatty liver rat model through nuclear receptors

The pathogenesis of non‐alcoholic fatty liver disease (NAFLD) is still not fully understood, and currently, no effective pharmacotherapy is available. Nuclear receptors (NRs) are important biological participants in NAFLD that exhibit great therapeutic potential. Chaihu Shugan powder (CSP) is a traditional Chinese medicine (TCM) formula that has a wide therapeutic spectrum including NAFLD, but the effective components and functional mechanisms of CSP are unclear. We adopted a network pharmacology approach using multiple databases for Gene Ontology (GO) enrichment analysis and the molecular complex detection (MCODE) method for a protein‐protein interaction (PPI) analysis, and we used molecular docking method to screen the NR targets and determine the corresponding CSP components. The screening results were validated through a NAFLD rat model that was used to explain the possible relationship between CSP and NAFLD. Finally, we screened PPARγ, FXR, PPARα, RARα and PPARδ as target genes and quercetin, kaempferol, naringenin, isorhamnetin and nobiletin as target compounds. The five components were detected through high‐performance liquid chromatography‐mass spectrometry (HPLC‐MS), the results of which aligned with the docking experiments of PPARγ, PPARα and PPARδ. After CSP intervention, the NAFLD rat model showed ameliorated effects in terms of bodyweight, hepatic histopathology, and serum and liver lipids, and the mRNA levels of PPARγ, FXR, PPARα and RARα were significantly changed. The results from this study indicate that CSP exhibits healing effects in an NAFLD model and that the network pharmacology approach to screening NR targets and determining the corresponding CSP components is a practical strategy for explaining the mechanism by which CSP ameliorates NAFLD.     

Finding off‐targets,biological pathways,and target diseases for chymase inhibitors via structure‐based systems biology approach

Off‐target binding connotes the binding of a small molecule of therapeutic significance to a protein target in addition to the primary target for which it was proposed. Progressively such off‐targeting is emerging to be regular practice to reveal side effects. Chymase is an enzyme of hydrolase class that catalyzes hydrolysis of peptide bonds. A link between heart failure and chymase is ascribed, and a chymase inhibitor is in clinical phase II for treatment of heart failure. However, the underlying mechanisms of the off‐target effects of human chymase inhibitors are still unclear. Here, we develop a robust computational strategy that is applicable to any enzyme system and that allows the prediction of drug effects on biological processes. Putative off‐targets for chymase inhibitors were identified through various structural and functional similarity analyses along with molecular docking studies. Finally, literature survey was performed to incorporate these off‐targets into biological pathways and to establish links between pathways and particular adverse effects. Off‐targets of chymase inhibitors are linked to various biological pathways such as classical and lectin pathways of complement system, intrinsic and extrinsic pathways of coagulation cascade, and fibrinolytic system. Tissue kallikreins, granzyme M, neutrophil elastase, and mesotrypsin are also identified as off‐targets. These off‐targets and their associated pathways are elucidated for the effects of inflammation, cancer, hemorrhage, thrombosis, and central nervous system diseases (Alzheimer's disease). Prospectively, our approach is helpful not only to better understand the mechanisms of chymase inhibitors but also for drug repurposing exercises to find novel uses for these inhibitors. Proteins 2015; 83:1209–1224. © 2014 Wiley Periodicals, Inc.     

Dynamic eicosanoid responses upon different inhibitor and combination treatments on the arachidonic acid metabolic network

The arachidonic acid (AA) metabolic network produces key inflammatory mediators which have been considered as hallmark contributors in various inflammatory related diseases. Enzymes in this network, such as 5-lipoxygenase (5-LOX), cyclooxygenase (COX), leukotriene A(4) hydrolase (LTA4H) and prostaglandin E synthase (PGES), have been used as targets for anti-inflammatory drug discovery. Multi-target drugs and drug combinations have also been developed for this network. However, how the inhibitors alter the dynamics of metabolite production and which combinatorial target intervention solutions are better needs further exploration. We did a system based intervention analysis on the AA metabolic network. Using an LC-MS/MS method, we quantitatively studied the eicosanoid metabolites responses of AA metabolic network during stimulation of Sprague Dawley rat blood samples with the calcium ionophore. Our results indicate that inhibiting the upstream rather than the downstream target of 5-LOX pathway will simultaneously alter the AA metabolism to the COX pathway (and vice versa). Therefore, single-target inhibitors cannot control all the inflammatory mediators at the same time. We also suggest that in the case of multiple-target anti-inflammatory solutions, the combination of inhibitors of the downstream enzymes may have stronger inhibition efficiency and cause less side-effects compared to the other solutions. One therapeutic strategy, LTA4H/COX inhibition solution, was found promising for the intervention of inflammatory mediator biosynthesis and at the same time stimulating the production of anti-inflammatory agents.     

基于生物信息学方法发现潜在药物靶标

药物靶点通常是在代谢或信号通路中与特定疾病或病理状态有关的关键分子.通过绑定到特定活动区域抑制这个关键分子进行药物设计.确定特定疾病有关的靶标分子是现代新药开发的基础.在药物靶标发现的过程中,生物信息学方法发挥了不可替代的重要的作用,尤其适用于大规模多组学数据的分析.目前,已涌现了许多与疾病相关的数据库资源,基于生物网络特征、多基因芯片、蛋白质组、代谢组数据等建立了多种生物信息学方法发现潜在的药物靶标,并预测靶标可药性和药物副作用.     

Network medicine

To more effectively target complex diseases like cancer, diabetes and schizophrenia, we may need to rethink our strategies for drug development and the selection of molecular targets for pharmacological treatments. Here, we discuss the potential use of protein signaling networks as the targets for new therapeutic intervention. We argue that by targeting the architecture of aberrant signaling networks associated with cancer and other diseases new therapeutic strategies can be implemented. Transforming medicine into a network driven endeavour will require quantitative measurements of cell signaling processes; we will describe how this may be performed and combined with new algorithms to predict the trajectories taken by a cellular system either in time or through disease states. We term this approach, network medicine.     

PREDICT: a method for inferring novel drug indications with application to personalized medicine

Inferring potential drug indications, for either novel or approved drugs, is a key step in drug development. Previous computational methods in this domain have focused on either drug repositioning or matching drug and disease gene expression profiles. Here, we present a novel method for the large‐scale prediction of drug indications (PREDICT) that can handle both approved drugs and novel molecules. Our method is based on the observation that similar drugs are indicated for similar diseases, and utilizes multiple drug–drug and disease–disease similarity measures for the prediction task. On cross‐validation, it obtains high specificity and sensitivity (AUC=0.9) in predicting drug indications, surpassing existing methods. We validate our predictions by their overlap with drug indications that are currently under clinical trials, and by their agreement with tissue‐specific expression information on the drug targets. We further show that disease‐specific genetic signatures can be used to accurately predict drug indications for new diseases (AUC=0.92). This lays the computational foundation for future personalized drug treatments, where gene expression signatures from individual patients would replace the disease‐specific signatures.     

Identification of control targets in Boolean molecular network models via computational algebra

Background

Many problems in biomedicine and other areas of the life sciences can be characterized as control problems, with the goal of finding strategies to change a disease or otherwise undesirable state of a biological system into another, more desirable, state through an intervention, such as a drug or other therapeutic treatment. The identification of such strategies is typically based on a mathematical model of the process to be altered through targeted control inputs. This paper focuses on processes at the molecular level that determine the state of an individual cell, involving signaling or gene regulation. The mathematical model type considered is that of Boolean networks. The potential control targets can be represented by a set of nodes and edges that can be manipulated to produce a desired effect on the system.

Results

This paper presents a method for the identification of potential intervention targets in Boolean molecular network models using algebraic techniques. The approach exploits an algebraic representation of Boolean networks to encode the control candidates in the network wiring diagram as the solutions of a system of polynomials equations, and then uses computational algebra techniques to find such controllers. The control methods in this paper are validated through the identification of combinatorial interventions in the signaling pathways of previously reported control targets in two well studied systems, a p53-mdm2 network and a blood T cell lymphocyte granular leukemia survival signaling network. Supplementary data is available online and our code in Macaulay2 and Matlab are available via http://www.ms.uky.edu/~dmu228/ControlAlg.

Conclusions

This paper presents a novel method for the identification of intervention targets in Boolean network models. The results in this paper show that the proposed methods are useful and efficient for moderately large networks.