Identification of Drug Combinations Containing Imatinib for Treatment of BCR-ABL+ Leukemias

The BCR-ABL translocation is found in chronic myeloid leukemia (CML) and in Ph+ acute lymphoblastic leukemia (ALL) patients. Although imatinib and its analogues have been used as front-line therapy to target this mutation and control the disease for over a decade, resistance to the therapy is still observed and most patients are not cured but need to continue the therapy indefinitely. It is therefore of great importance to find new therapies, possibly as drug combinations, which can overcome drug resistance. In this study, we identified eleven candidate anti-leukemic drugs that might be combined with imatinib, using three approaches: a kinase inhibitor library screen, a gene expression correlation analysis, and literature analysis. We then used an experimental search algorithm to efficiently explore the large space of possible drug and dose combinations and identified drug combinations that selectively kill a BCR-ABL+ leukemic cell line (K562) over a normal fibroblast cell line (IMR-90). Only six iterations of the algorithm were needed to identify very selective drug combinations. The efficacy of the top forty-nine combinations was further confirmed using Ph+ and Ph- ALL patient cells, including imatinib-resistant cells. Collectively, the drug combinations and methods we describe might be a first step towards more effective interventions for leukemia patients, especially those with the BCR-ABL translocation.     

Prediction of drug combinations by integrating molecular and pharmacological data

Combinatorial therapy is a promising strategy for combating complex disorders due to improved efficacy and reduced side effects. However, screening new drug combinations exhaustively is impractical considering all possible combinations between drugs. Here, we present a novel computational approach to predict drug combinations by integrating molecular and pharmacological data. Specifically, drugs are represented by a set of their properties, such as their targets or indications. By integrating several of these features, we show that feature patterns enriched in approved drug combinations are not only predictive for new drug combinations but also provide insights into mechanisms underlying combinatorial therapy. Further analysis confirmed that among our top ranked predictions of effective combinations, 69% are supported by literature, while the others represent novel potential drug combinations. We believe that our proposed approach can help to limit the search space of drug combinations and provide a new way to effectively utilize existing drugs for new purposes.     

Multiobjective optimization identifies cancer-selective combination therapies

Combinatorial therapies are required to treat patients with advanced cancers that have become resistant to monotherapies through rewiring of redundant pathways. Due to a massive number of potential drug combinations, there is a need for systematic approaches to identify safe and effective combinations for each patient, using cost-effective methods. Here, we developed an exact multiobjective optimization method for identifying pairwise or higher-order combinations that show maximal cancer-selectivity. The prioritization of patient-specific combinations is based on Pareto-optimization in the search space spanned by the therapeutic and nonselective effects of combinations. We demonstrate the performance of the method in the context of BRAF-V600E melanoma treatment, where the optimal solutions predicted a number of co-inhibition partners for vemurafenib, a selective BRAF-V600E inhibitor, approved for advanced melanoma. We experimentally validated many of the predictions in BRAF-V600E melanoma cell line, and the results suggest that one can improve selective inhibition of BRAF-V600E melanoma cells by combinatorial targeting of MAPK/ERK and other compensatory pathways using pairwise and third-order drug combinations. Our mechanism-agnostic optimization method is widely applicable to various cancer types, and it takes as input only measurements of a subset of pairwise drug combinations, without requiring target information or genomic profiles. Such data-driven approaches may become useful for functional precision oncology applications that go beyond the cancer genetic dependency paradigm to optimize cancer-selective combinatorial treatments.     

Combination therapies for combating antimicrobial resistance

New drug development strategies are needed to combat antimicrobial resistance. The object of this perspective is to highlight one such strategy: treating infections with sets of drugs rather than individual drugs. We will highlight three categories of combination therapy: those that inhibit targets in different pathways; those that inhibit distinct nodes in the same pathway; and those that inhibit the very same target in different ways. We will then consider examples of naturally occurring combination therapies produced by micro-organisms, and conclude by discussing key opportunities and challenges for making more widespread use of drug combinations.     

In silico drug combination discovery for personalized cancer therapy

Background

Drug combination therapy, which is considered as an alternative to single drug therapy, can potentially reduce resistance and toxicity, and have synergistic efficacy. As drug combination therapies are widely used in the clinic for hypertension, asthma, and AIDS, they have also been proposed for the treatment of cancer. However, it is difficult to select and experimentally evaluate effective combinations because not only is the number of cancer drug combinations extremely large but also the effectiveness of drug combinations varies depending on the genetic variation of cancer patients. A computational approach that prioritizes the best drug combinations considering the genetic information of a cancer patient is necessary to reduce the search space.

Results

We propose an in-silico method for personalized drug combination therapy discovery. We predict the synergy between two drugs and a cell line using genomic information, targets of drugs, and pharmacological information. We calculate and predict the synergy scores of 583 drug combinations for 31 cancer cell lines. For feature dimension reduction, we select the mutations or expression levels of the genes in cancer-related pathways. We also used various machine learning models. Extremely Randomized Trees (ERT), a tree-based ensemble model, achieved the best performance in the synergy score prediction regression task. The correlation coefficient between the synergy scores predicted by ERT and the actual observations is 0.738. To compare with an existing drug combination synergy classification model, we reformulate the problem as a binary classification problem by thresholding the synergy scores. ERT achieved an F1 score of 0.954 when synergy scores of 20 and -20 were used as the threshold, which is 8.7% higher than that obtained by the state-of-the-art baseline model. Moreover, the model correctly predicts the most synergistic combination, from approximately 100 candidate drug combinations, as the top choice for 15 out of the 31 cell lines. For 28 out of the 31 cell lines, the model predicts the most synergistic combination in the top 10 of approximately 100 candidate drug combinations. Finally, we analyze the results, generate synergistic rules using the features, and validate the rules through the literature survey.

Conclusion

Using various types of genomic information of cancer cell lines, targets of drugs, and pharmacological information, a drug combination synergy prediction pipeline is proposed. The pipeline regresses the synergy level between two drugs and a cell line as well as classifies if there exists synergy or antagonism between them. Discovering new drug combinations by our pipeline may improve personalized cancer therapy.

     

Magnetic resonance spectroscopy detectable metabolomic fingerprint of response to antineoplastic treatment

Targeted therapeutic approaches are increasingly being implemented in the clinic, but early detection of response frequently presents a challenge as many new therapies lead to inhibition of tumor growth rather than tumor shrinkage. Development of novel non-invasive methods to monitor response to treatment is therefore needed. Magnetic resonance spectroscopy (MRS) and magnetic resonance spectroscopic imaging are non-invasive imaging methods that can be employed to monitor metabolism, and previous studies indicate that these methods can be useful for monitoring the metabolic consequences of treatment that are associated with early drug target modulation. However, single-metabolite biomarkers are often not specific to a particular therapy. Here we used an unbiased 1H MRS-based metabolomics approach to investigate the overall metabolic consequences of treatment with the phosphoinositide 3-kinase inhibitor LY294002 and the heat shock protein 90 inhibitor 17AAG in prostate and breast cancer cell lines. LY294002 treatment resulted in decreased intracellular lactate, alanine fumarate, phosphocholine and glutathione. Following 17AAG treatment, decreased intracellular lactate, alanine, fumarate and glutamine were also observed but phosphocholine accumulated in every case. Furthermore, citrate, which is typically observed in normal prostate tissue but not in tumors, increased following 17AAG treatment in prostate cells. This approach is likely to provide further information about the complex interactions between signaling and metabolic pathways. It also highlights the potential of MRS-based metabolomics to identify metabolic signatures that can specifically inform on molecular drug action.     

Prediction in the face of uncertainty: a Monte Carlo-based approach for systems biology of cancer treatment

Cancer is known to be a complex disease and its therapy is difficult. Much information is available on molecules and pathways involved in cancer onset and progression and this data provides a valuable resource for the development of predictive computer models that can help to identify new potential drug targets or to improve therapies. Modeling cancer treatment has to take into account many cellular pathways usually leading to the construction of large mathematical models. The development of such models is complicated by the fact that relevant parameters are either completely unknown, or can at best be measured under highly artificial conditions. Here we propose an approach for constructing predictive models of such complex biological networks in the absence of accurate knowledge on parameter values, and apply this strategy to predict the effects of perturbations induced by anti-cancer drug target inhibitions on an epidermal growth factor (EGF) signaling network. The strategy is based on a Monte Carlo approach, in which the kinetic parameters are repeatedly sampled from specific probability distributions and used for multiple parallel simulations. Simulation results from different forms of the model (e.g., a model that expresses a certain mutation or mutation pattern or the treatment by a certain drug or drug combination) can be compared with the unperturbed control model and used for the prediction of the perturbation effects. This framework opens the way to experiment with complex biological networks in the computer, likely to save costs in drug development and to improve patient therapy.     

A quantitative basis for antiretroviral therapy for HIV-1 infection

Highly active antiretroviral therapy (HAART) has dramatically decreased mortality from HIV-1 infection and is a major achievement of modern medicine. However, there is no fundamental theory of HAART. Elegant models describe the dynamics of viral replication, but a metric for the antiviral activity of drug combinations relative to a target value needed for control of replication is lacking. Treatment guidelines are based on empirical results of clinical trials in which other factors such as regimen tolerability also affect outcome. Why only certain drug combinations control viral replication remains unclear. Here we quantify the intrinsic antiviral activity of antiretroviral drug combinations. We show that most single antiretroviral drugs show previously unappreciated complex nonlinear pharmacodynamics that determine their inhibitory potential at clinical concentrations. We demonstrate that neither of the major theories for drug combinations accurately predicts the combined effects of multiple antiretrovirals. However, the combined effects can be understood with a new approach that considers the degree of independence of drug effects. This analysis allows a direct comparison of the inhibitory potential of different drug combinations under clinical concentrations, reconciles the results of clinical trials, defines a target level of inhibition associated with treatment success and provides a rational basis for treatment simplification and optimization.     

Mathematical modeling of multi-drugs therapy: a challenge for determining the optimal combinations of antiviral drugs

In the current era of antiviral drug therapy, combining multiple drugs is a primary approach for improving antiviral effects, reducing the doses of individual drugs, relieving the side effects of strong antiviral drugs, and preventing the emergence of drug-resistant viruses. Although a variety of new drugs have been developed for HIV, HCV and influenza virus, the optimal combinations of multiple drugs are incompletely understood. To optimize the benefits of multi-drugs combinations, we must investigate the interactions between the combined drugs and their target viruses. Mathematical models of viral infection dynamics provide an ideal tool for this purpose. Additionally, whether drug combinations computed by these models are synergistic can be assessed by two prominent drug combination theories, Loewe additivity and Bliss independence. By combining the mathematical modeling of virus dynamics with drug combination theories, we could show the principles by which drug combinations yield a synergistic effect. Here, we describe the theoretical aspects of multi-drugs therapy and discuss their application to antiviral research.     

Optimizing combination chemotherapy by controlling drug ratios

Cancer chemotherapy treatments typically employ drug combinations in which the dose of each agent is pushed to the brink of unacceptable toxicity; however, emerging evidence indicates that this approach may not be providing optimal efficacy due to the manner in which drugs interact. Specifically, whereas certain ratios of combined drugs can be synergistic, other ratios of the same agents may be antagonistic, implying that the most efficacious combinations may be those that utilize certain agents at reduced doses. Advances in nano-scale drug delivery vehicles now enable the translation of in vitro information on synergistic drug ratios into improved anticancer combination therapies in which the desired drug ratio can be controlled and maintained following administration in vivo, so that synergistic effects can be exploited. This "ratiometric" approach to combination chemotherapy opens new opportunities to enhance the effectiveness of existing and future treatment regimens across a spectrum of human diseases.     

p70S6 kinase is a critical node that integrates HER-family and PI3 kinase signaling networks

Therapies targeting oncogenic drivers rapidly induce compensatory adaptive responses that blunt drug effectiveness, contributing to therapeutic resistance. Adaptive responses are characteristic of robust cell signaling networks, and thus there is increasing interest in drug combinations that co-target the driver and the adaptive response. An alternative approach to co-inhibiting oncogenic and adaptive targets is to identify a critical node where the activities of these targets converge. Nodes of convergence between signaling modules represent potential therapeutic vulnerabilities because their inhibition could result in the collapse of the network, leading to enhanced cytotoxicity. In this report we demonstrate that p70S6 kinase (p70S6K) can function as a critical node linking HER-family and phosphoinositide-3-kinase (PI3K) pathway signaling. We used high-throughput combinatorial drug screening to identify adaptive survival responses to targeted therapies, and found that HER-family and PI3K represented compensatory signaling pathways. Co-targeting these pathways with drug combinations caused synergistic cytotoxicity in cases where inhibition of neither target was effective as a monotherapy. We utilized Reverse Phase Protein Arrays and determined that phosphorylation of ribosomal protein S6 was synergistically down-regulated upon HER-family and PI3K/mammalian target of rapamycin (mTOR) co-inhibition. Expression of constitutively active p70S6K protected against apoptosis induced by combined HER-family and PI3K/mTOR inhibition. Direct inhibition of p70S6K with small molecule inhibitors phenocopied HER-family and PI3K/mTOR co-inhibition. These data implicate p70S6K as a critical node in the HER-family/PI3K signaling network. The ability of direct inhibitors of p70S6K to phenocopy co-inhibition of two upstream signaling targets indicates that identification and targeting of critical nodes can overcome adaptive resistance to targeted therapies.     

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.     

Modeling the short-, middle- and long-term viral load responses for comparing estimated dynamic parameters

A virologic marker, the number of HIV RNA copies or viral load, is currently used to evaluate antiviral therapies in AIDS clinical trials. This marker can be used to assess the antiviral potency of therapies, but is easily affected by drug exposures, drug resistance and other factors during the long-term treatment evaluation process. The study of HIV dynamics is one of the most important development in recent AIDS research for understanding the pathogenesis of HIV-1 infection and antiviral treatment strategies. Although many HIV dynamic models have been proposed by AIDS researchers in the last decade, they have only been used to quantify short-term viral dynamics and do not correctly describe long-term virologic responses to antiretroviral treatment. In other words, these simple viral dynamic models can only be used to fit short-term viral load data for estimating dynamic parameters. In this paper, a mechanism-based differential equation models is introduced for characterizing the long-term viral dynamics with antiretroviral therapy. We applied this model to fit different segments of the viral load trajectory data from a simulation experiment and an AIDS clinical trial study, and found that the estimates of dynamic parameters from our modeling approach are very consistent. We may conclude that our model can not only characterize long-term viral dynamics, but can also quantify short- and middle-term viral dynamics. It suggests that if there are enough data in the early stage of the treatment, the results from our modeling based on short-term information can be used to capture the performance of long-term care with HIV-1 infected patients.     

基于深度学习与多层次信息融合的药物靶标亲和力预测*

药物研发是非常重要但也十分耗费人力物力的过程。利用计算机辅助预测药物与蛋白质亲和力的方法可以极大地加快药物研发过程。药物靶标亲和力预测的关键在于对药物和蛋白质进行准确详细地信息表征。提出一种基于深度学习与多层次信息融合的药物靶标亲和力的预测模型,试图通过综合药物与蛋白质的多层次信息,来获得更好的预测表现。首先将药物表述成分子图和扩展连接指纹两种形式,分别利用图卷积神经网络模块和全连接层进行学习;其次将蛋白质序列和蛋白质K-mer特征分别输入卷积神经网络模块和全连接层来学习蛋白质潜在特征;随后将4个通道学习到的特征进行融合,再利用全连接层进行预测。在两个基准药物靶标亲和力数据集上验证了所提方法的有效性,并与其他已有模型作对比研究。结果说明提出的模型相比基准模型能得到更好的预测性能,表明提出的综合药物与蛋白质多层次信息的药物靶标亲和力预测策略是有效的。     

Simulating Serial-Target Antibacterial Drug Synergies Using Flux Balance Analysis

Flux balance analysis (FBA) is an increasingly useful approach for modeling the behavior of metabolic systems. However, standard FBA modeling of genetic knockouts cannot predict drug combination synergies observed between serial metabolic targets, even though such synergies give rise to some of the most widely used antibiotic treatments. Here we extend FBA modeling to simulate responses to chemical inhibitors at varying concentrations, by diverting enzymatic flux to a waste reaction. This flux diversion yields very similar qualitative predictions to prior methods for single target activity. However, we find very different predictions for combinations, where flux diversion, which mimics the kinetics of competitive metabolic inhibitors, can explain serial target synergies between metabolic enzyme inhibitors that we confirmed in Escherichia coli cultures. FBA flux diversion opens the possibility for more accurate genome-scale predictions of drug synergies, which can be used to suggest treatments for infections and other diseases.     

基于深度学习的蛋白质建模与设计

随着蛋白质序列及结构数据的大量累积,在获得了大量描述性信息之后如何有效利用海量数据,从已有数据中高效提取信息并且应用到下游任务当中就成为了研究者亟待解决的问题。蛋白质的设计可使新蛋白的研发不再受限于实验条件,这对药物靶点预测、新药研发和材料设计等领域具有重要意义。深度学习作为一种高效的数据特征提取方法,可以通过它对蛋白质数据进行建模,进而加入先验信息对蛋白质进行设计。故此基于深度学习的蛋白质设计就成为一个具有广阔前景的研究领域。文中主要阐述基于深度学习的蛋白质序列与结构数据的建模和设计方法。详述该方法的策略、原理、适用范围、应用实例。讨论了深度学习方法在本领域的应用前景及局限性,以期为相关研究提供参考。     

Rational deployment of antimalarial drugs in Africa: should first-line combination drugs be reserved for paediatric malaria cases?

Artemisinin-based combination therapy is exerting novel selective pressure upon populations of Plasmodium falciparum across Africa. Levels of resistance to non-artemisinin partner drugs differ among parasite populations, and so the artemisinins are not uniformly protected from developing resistance, already present in South East Asia. Here, we consider strategies for prolonging the period of high level efficacy of combination therapy for two particular endemicities common in Africa. Under high intensity transmission, two alternating first-line combinations, ideally with antagonistic selective effects on the parasite genome, are advocated for paediatric malaria cases. This leaves second-line and other therapies for adult cases, and for intermittent preventive therapy. The drug portfolio would be selected to protect the 'premier' combination regimen from selection for resistance, while maximising impact on severe disease and mortality in children. In endemic areas subject to low, seasonal transmission of Plasmodium falciparum, such a strategy may deliver little benefit, as children represent a minority of cases. Nevertheless, the deployment of other drug-based interventions in low transmission and highly seasonal areas, such as mass drug administration aimed to interrupt malaria transmission, or intermittent preventive therapy, does provide an opportunity to diversify drug pressure. We thus propose an integrated approach to drug deployment, which minimises direct selective pressure on parasite populations from any one drug component. This approach is suitable for qualitatively and quantitatively different burdens of malaria, and should be supported by a programme of routine surveillance for emerging resistance.     

Efficient Modeling and Active Learning Discovery of Biological Responses

High throughput and high content screening involve determination of the effect of many compounds on a given target. As currently practiced, screening for each new target typically makes little use of information from screens of prior targets. Further, choices of compounds to advance to drug development are made without significant screening against off-target effects. The overall drug development process could be made more effective, as well as less expensive and time consuming, if potential effects of all compounds on all possible targets could be considered, yet the cost of such full experimentation would be prohibitive. In this paper, we describe a potential solution: probabilistic models that can be used to predict results for unmeasured combinations, and active learning algorithms for efficiently selecting which experiments to perform in order to build those models and determining when to stop. Using simulated and experimental data, we show that our approaches can produce powerful predictive models without exhaustive experimentation and can learn them much faster than by selecting experiments at random.