Systems biology approaches for advancing the discovery of effective drug combinations

Systems biology: perspectives on multiscale modeling in research on endocrine-related cancers

Drawing on concepts from experimental biology, computer science, informatics, mathematics and statistics, systems biologists integrate data across diverse platforms and scales of time and space to create computational and mathematical models of the integrative, holistic functions of living systems. Endocrine-related cancers are well suited to study from a systems perspective because of the signaling complexities arising from the roles of growth factors, hormones and their receptors as critical regulators of cancer cell biology and from the interactions among cancer cells, normal cells and signaling molecules in the tumor microenvironment. Moreover, growth factors, hormones and their receptors are often effective targets for therapeutic intervention, such as estrogen biosynthesis, estrogen receptors or HER2 in breast cancer and androgen receptors in prostate cancer. Given the complexity underlying the molecular control networks in these cancers, a simple, intuitive understanding of how endocrine-related cancers respond to therapeutic protocols has proved incomplete and unsatisfactory. Systems biology offers an alternative paradigm for understanding these cancers and their treatment. To correctly interpret the results of systems-based studies requires some knowledge of how in silico models are built, and how they are used to describe a system and to predict the effects of perturbations on system function. In this review, we provide a general perspective on the field of cancer systems biology, and we explore some of the advantages, limitations and pitfalls associated with using predictive multiscale modeling to study endocrine-related cancers.

Perturbational Gene-Expression Signatures for Combinatorial Drug Discovery

Cancer is a complex disease that relies on both oncogenic mutations and non-mutated genes for survival, and therefore coined as oncogene and non-oncogene addictions. The need for more effective combination therapies to overcome drug resistance in oncology has been increasingly recognized, but the identification of potentially synergistic drugs at scale remains challenging. Here we propose a gene-expression-based approach, which uses the recurrent perturbation-transcript regulatory relationships inferred from a large compendium of chemical and genetic perturbation experiments across multiple cell lines, to engender a testable hypothesis for combination therapies. These transcript-level recurrences were distinct from known compound-protein target counterparts, were reproducible in external datasets, and correlated with small-molecule sensitivity. We applied these recurrent relationships to predict synergistic drug pairs for cancer and experimentally confirmed two unexpected drug combinations in vitro. Our results corroborate a gene-expression-based strategy for combinatorial drug screening as a way to target non-mutated genes in complex diseases.

•

Compound signatures targeting non-oncogene addiction for combinatorial drug discovery

•

These signatures are reproducible and linked to cancer hallmarks and drug sensitivity

•

Two synergistic drug combinations are experimentally confirmed in vitro

Systems biology in inflammatory bowel diseases: on the way to precision medicine

Inflammatory bowel diseases (IBD) are chronic and recurrent inflammatory disorders of the gastrointestinal tract. The elucidation of their etiopathology requires complex and multiple approaches. Systems biology has come to fulfill this need in approaching the pathogenetic mechanisms of IBD and its etiopathology, in a comprehensive way, by combining data from different scientific sources. In combination with bioinformatics and network medicine, it uses principles from computer science, mathematics, physics, chemistry, biology, medicine and computational tools to achieve its purposes. Systems biology utilizes scientific sources that provide data from omics studies (e.g., genomics, transcriptomics, etc.) and clinical observations, whose combined analysis leads to network formation and ultimately to a more integrative image of disease etiopathogenesis. In this review, we analyze the current literature on the methods and the tools utilized by systems biology in order to cover an innovative and exciting field: IBD-omics.

Recent perspectives on the virulent factors and treatment options for multidrug-resistant Acinetobacter baumannii

Acinetobacter baumannii (AB) is one of the most notorious and opportunistic pathogens, which caused high morbidity and mortality rate and World Health Organization (WHO) declared this bacterium as priority-1 pathogen in 2017. The current antibacterial agents, such as colistins, carbapenems, and tigecyclines have limited applications, which necessitate novel and alternative therapeutic remedies. Thus, the understanding of recent perspectives on the virulent factors and antibiotic resistance mechanism exhibited by the bacteria are extremely important. In addition to many combinatorial therapies of antibacterial, there is several natural compounds demonstrated significant antibacterial potential towards these bacteria. The computational systems biology and high throughput screening approaches provide crucial insights in identifying novel drug targets and lead molecules with therapeutics potential. Hence, this review provides profound insight on the recent aspects of the virulent factors associated with AB, role of biofilm formation in drug resistance and the mechanisms of multidrug resistance. This review further illustrates the status of current therapeutic agents, scope, and applications of natural therapeutics, such as herbal medicines and role of computational biology, immunoinformatics and virtual screening in novel lead developments. Thus, this review provides novel insight on latest developments in drug-resistance mechanism of multidrug-resistant A. baumannii (MDRAB) and discovery of probable therapeutic interventions.

Computational Approaches in Theranostics: Mining and Predicting Cancer Data

The ability to understand the complexity of cancer-related data has been prompted by the applications of (1) computer and data sciences, including data mining, predictive analytics, machine learning, and artificial intelligence, and (2) advances in imaging technology and probe development. Computational modelling and simulation are systematic and cost-effective tools able to identify important temporal/spatial patterns (and relationships), characterize distinct molecular features of cancer states, and address other relevant aspects, including tumor detection and heterogeneity, progression and metastasis, and drug resistance. These approaches have provided invaluable insights for improving the experimental design of therapeutic delivery systems and for increasing the translational value of the results obtained from early and preclinical studies. The big question is: Could cancer theranostics be determined and controlled in silico? This review describes the recent progress in the development of computational models and methods used to facilitate research on the molecular basis of cancer and on the respective diagnosis and optimized treatment, with particular emphasis on the design and optimization of theranostic systems. The current role of computational approaches is providing innovative, incremental, and complementary data-driven solutions for the prediction, simplification, and characterization of cancer and intrinsic mechanisms, and to promote new data-intensive, accurate diagnostics and therapeutics.

Freshwater neotropical oligochaetes as native test species for the toxicity evaluation of cadmium, mercury and their mixtures

The toxicity of metals, whether isolated or in mixtures, involves changes in biochemical processes as well as in cell membranes, which may lead to deleterious short- and long-term effects on the affected organisms. Among metals, cadmium and mercury stand out due to their abundance in nature, frequent use for industrial activities and biological accumulation, with high levels of residence in trophic chains. Benthic communities are particularly prone to metal pollution since metals usually accumulate in sediments. The aim of this study was to evaluate the acute toxicity of mercury and cadmium, single and in mixture, to two native species of epibenthic oligochaetes: Allonais inaequalis and Dero furcatus. In order to assess the potential of these species as bioindicators, we compared their sensitivity with those of other internationally used species by applying the species sensitivity distribution approach. The 96h-LC₅₀ of cadmium chloride was 627 and 364 μg L^-1 for A. inaequalis and D. furcatus, respectively, evidencing that the latter species is almost twice as sensitive to this metal than A. inaequalis. For mercury chloride, the 96h-LC₅₀ was 129 μg L^-1 for A. inaequalis and 92 μg L^-1 for D. furcatus. The sensitivities of these oligochaetes were superior or similar to that of other frequently used oligochaete test species such as Tubifex tubifex and Lumbriculus variegatus. The metal mixtures had synergism in general (D. furcatus) or at high doses only (A. inaequalis), implying a potentiation of their toxic effects when both metals co-occur in the environment. By comparing the derived toxicity values with concentrations of cadmium and mercury measured in the field, it can be concluded that aquatic organisms are likely to be at risk when exposed to the environmental relevant concentrations of cadmium and mercury here tested, especially when they are both present.

Predictive approaches for drug combination discovery in cancer

Drug combinations have been proposed as a promising therapeutic strategy to overcome drug resistance and improve efficacy of monotherapy regimens in cancer. This strategy aims at targeting multiple components of this complex disease. Despite the increasing number of drug combinations in use, many of them were empirically found in the clinic, and the molecular mechanisms underlying these drug combinations are often unclear. These challenges call for rational, systematic approaches for drug combination discovery. Although high-throughput screening of single-agent therapeutics has been successfully implemented, it is not feasible to test all possible drug combinations, even for a reduced subset of anticancer drugs. Hence, in vitro and in vivo screening of a large number of drug combinations are not practical. Therefore, devising computational methods to efficiently explore the space of drug combinations and to discover efficacious combinations has attracted a lot of attention from the scientific community in the past few years. Nevertheless, in the absence of consensus regarding the computational approaches used to predict efficacious drug combinations, a plethora of methods, techniques and hypotheses have been developed to date, while the research field lacks an elaborate categorization of the existing computational methods and the available data sources. In this manuscript, we review and categorize the state-of-the-art computational approaches for drug combination prediction, and elaborate on the limitations of these methods and the existing challenges. We also discuss about the recent pan-cancer drug combination data sets and their importance in revising the available methods or developing more performant approaches.

Synergy from gene expression and network mining (SynGeNet) method predicts synergistic drug combinations for diverse melanoma genomic subtypes

Systems biology perspectives are crucial for understanding the pathophysiology of complex diseases, and therefore hold great promise for the discovery of novel treatment strategies. Drug combinations have been shown to improve durability and reduce resistance to available first-line therapies in a variety of cancers; however, traditional drug discovery approaches are prohibitively cost and labor-intensive to evaluate large-scale matrices of potential drug combinations. Computational methods are needed to efficiently model complex interactions of drug target pathways and identify mechanisms underlying drug combination synergy. In this study, we employ a computational approach, SynGeNet (Synergy from Gene expression and Network mining), which integrates transcriptomics-based connectivity mapping and network centrality analysis to analyze disease networks and predict drug combinations. As an exemplar of a disease in which combination therapies demonstrate efficacy in genomic-specific contexts, we investigate malignant melanoma. We employed SynGeNet to generate drug combination predictions for each of the four major genomic subtypes of melanoma (BRAF, NRAS, NF1, and triple wild type) using publicly available gene expression and mutation data. We validated synergistic drug combinations predicted by our method across all genomic subtypes using results from a high-throughput drug screening study across. Finally, we prospectively validated the drug combination for BRAF-mutant melanoma that was top ranked by our approach, vemurafenib (BRAF inhibitor) + tretinoin (retinoic acid receptor agonist), using both in vitro and in vivo models of BRAF-mutant melanoma and RNA-sequencing analysis of drug-treated melanoma cells to validate the predicted mechanisms. Our approach is applicable to a wide range of disease domains, and, importantly, can model disease-relevant protein subnetworks in precision medicine contexts.

Dynamic Targeting in Cancer Treatment

With the advent of personalized medicine, design and development of anti-cancer drugs that are specifically targeted to individual or sets of genes or proteins has been an active research area in both academia and industry. The underlying motivation for this approach is to interfere with several pathological crosstalk pathways in order to inhibit or at the very least control the proliferation of cancer cells. However, after initially conferring beneficial effects, if sub-lethal, these artificial perturbations in cell function pathways can inadvertently activate drug-induced up- and down-regulation of feedback loops, resulting in dynamic changes over time in the molecular network structure and potentially causing drug resistance as seen in clinics. Hence, the targets or their combined signatures should also change in accordance with the evolution of the network (reflected by changes to the structure and/or functional output of the network) over the course of treatment. This suggests the need for a "dynamic targeting" strategy aimed at optimizing tumor control by interfering with different molecular targets, at varying stages. Understanding the dynamic changes of this complex network under various perturbed conditions due to drug treatment is extremely challenging under experimental conditions let alone in clinical settings. However, mathematical modeling can facilitate studying these effects at the network level and beyond, and also accelerate comparison of the impact of different dosage regimens and therapeutic modalities prior to sizeable investment in risky and expensive clinical trials. A dynamic targeting strategy based on the use of mathematical modeling can be a new, exciting research avenue in the discovery and development of therapeutic drugs.

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.

Bioinformatics Analysis of Key Differentially Expressed Genes in Nonalcoholic Fatty Liver Disease Mice Models

Nonalcoholic fatty liver disease (NAFLD) is a global health problem characterized by excessive accumulation of fat in the liver without effect of other pathological factors including hepatitis infection and alcohol abuse. Current studies indicate that gene factors play important roles in the development of NAFLD. However, the molecular characteristics of differentially expressed genes (DEGs) and associated mechanisms with NAFLD have not been well elucidated. Using two microarray data associated with the gene expression profiling in liver tissues of NAFLD mice models, we identified and selected several common key DEGs that contributed to NAFLD. Based on bioinformatics analysis, we discovered that the DEGs were associated with a variety of biological processes, cellular components, and molecular functions and were also related to several significant pathways. Via pathway crosstalk analysis based on overlapping DEGs, we observed that the identified pathways could form large and complex crosstalk networks. Besides, large and complex protein interaction networks of DEGs were further constructed. In addition, many hub host factors with a high degree of connectivity were identified based on interaction networks. Furthermore, significant modules in interaction networks were found, and the DEGs in the identified modules were found to be enriched with distinct pathways. Taken together, these results suggest that the key DEGs, associated pathways, and modules contribute to the development of NAFLD and might be used as novel molecular targets for the treatment of NAFLD.

Sensitivities of three tropical indigenous freshwater invertebrates to single and mixture exposures of diuron and carbofuran and their commercial formulations

As compared to their temperate counterparts, few toxicity tests have been conducted so far into the evaluation of the sensitivity of indigenous tropical species to pesticides. Especially mixture toxicity assessments appear to be scarce. To contribute to increase our knowledge in this arena, we evaluated the acute toxicity of diuron and carbofuran and their mixtures to the neotropical oligochaetes Allonais inaequalis and Dero furcatus, and the ostracod Strandesia trispinosa. Tests were performed with both the pure active ingredients, as well as their formulated products. The toxicity of the latter to the three test organisms was generally greater than that of the pure active ingredients, although absolute differences were rather small. The sensitivity of the indigenous species was slightly greater than temperate test species from the same taxonomic groups. The concentration addition conceptual model best described the results of the mixture toxicity data. Derived deviations of this model appeared to be dependent on the test organism and as to whether the pesticides were applied as active ingredients or their commercial products. Reported field concentrations of the two pesticides indicate risks to freshwater biota, especially if they are both present. The test species used in the present study are concluded to be suitable candidates as surrogate test organisms in local pesticide risk evaluations.

Mechanism-Based Disease Progression Model Describing Host-Pathogen Interactions During the Pathogenesis of Acinetobacter baumannii Pneumonia

The emergence of highly resistant bacteria is a serious threat to global public health. The host immune response is vital for clearing bacteria from the infected host; however, the current drug development paradigm does not take host‐pathogen interactions into consideration. Here, we used a systems‐based approach to develop a quantitative, mechanism‐based disease progression model to describe bacterial dynamics, host immune response, and lung injury in an immunocompetent rat pneumonia model. Previously, Long‐Evans rats were infected with Acinetobacter baumannii (A. baumannii) strain 307‐0294 at five different inocula and total lung bacteria, interleukin‐1beta (IL‐1β), tumor necrosis factor‐α (TNF‐α), cytokine‐induced neutrophil chemoattractant 1 (CINC‐1), neutrophil counts, and albumin were quantified. Model development was conducted in ADAPT5 version 5.0.54 using a pooled approach with maximum likelihood estimation; all data were co‐modeled. The final model characterized host‐pathogen interactions during the natural time course of bacterial pneumonia. Parameters were estimated with good precision. Our expandable model will integrate drug effects to aid in the design of optimized antibiotic regimens.

Immuno-Oncology: Emerging Targets and Combination Therapies

Host immunity recognizes and eliminates most early tumor cells, yet immunological checkpoints, exemplified by CTLA-4, PD-1, and PD-L1, pose a significant obstacle to effective antitumor immune responses. T-lymphocyte co-inhibitory pathways influence intensity, inflammation and duration of antitumor immunity. However, tumors and their immunosuppressive microenvironments exploit them to evade immune destruction. Recent PD-1 checkpoint inhibitors yielded unprecedented efficacies and durable responses across advanced-stage melanoma, showcasing potential to replace conventional radiotherapy regimens. Neverthless, many clinical problems remain in terms of efficacy, patient-to-patient variability, and undesirable outcomes and side effects. In this review, we evaluate recent advances in the immuno-oncology field and discuss ways forward. First, we give an overview of current immunotherapy modalities, involving mainy single agents, including inhibitor monoclonal antibodies (mAbs) targeting T-cell checkpoints of PD-1 and CTLA-4. However, neoantigen recognition alone cannot eliminate tumors effectively in vivo given their inherent complex micro-environment, heterogeneous nature and stemness. Then, based mainly upon CTLA-4 and PD-1 checkpoint inhibitors as a "backbone," we cover a range of emerging ("second-generation") therapies incorporating other immunotherapies or non-immune based strategies in synergistic combination. These include targeted therapies such as tyrosine kinase inhibitors, co-stimulatory mAbs, bifunctional agents, epigenetic modulators (such as inhibitors of histone deacetylases or DNA methyltransferase), vaccines, adoptive-T-cell therapy, nanoparticles, oncolytic viruses, and even synthetic "gene circuits." A number of novel immunotherapy co-targets in pre-clinical development are also introduced. The latter include metabolic components, exosomes and ion channels. We discuss in some detail of the personalization of immunotherapy essential for ultimate maximization of clinical outcomes. Finally, we outline possible future technical and conceptual developments including realistic in vitro and in vivo models and inputs from physics, engineering, and artificial intelligence. We conclude that the breadth and quality of immunotherapeutic approaches and the types of cancers that can be treated will increase significantly in the foreseeable future.

Ensemble Prediction of Synergistic Drug Combinations Incorporating Biological, Chemical, Pharmacological, and Network Knowledge

Combinatorial therapy may reduce drug side effects and improve drug efficacy, making combination therapy a promising strategy to treat complex diseases. However, in the existing computational methods, the natural properties and network knowledge of drugs have not been adequately and simultaneously considered, making it difficult to identify effective drug combinations. Computational methods that incorporate multiple sources of information (biological, chemical, pharmacological, and network knowledge) offer more opportunities to screen synergistic drug combinations. Therefore, we developed a novel Ensemble Prediction framework of Synergistic Drug Combinations (EPSDC) to accurately and efficiently predict drug combinations by integrating information from multiple-sources. EPSDC constructs feature vector of drug pair by concatenating different types of drug similarities, and then uses these groups in a feature-based base predictor. Next, transductive learning is applied on heterogeneous drug-target networks to achieve a network-based score for the drug pair. Finally, two types of ensemble rules are introduced to combine the feature-based score and the network-based score, and then potential drug combinations are prioritized. To demonstrate the effect of the ensemble rule, comprehensive experiments were conducted to compare single models and ensemble models. The experimental results indicated that our method outperformed the state-of-the-art method in five-fold cross validation and de novo prediction tests on the two benchmark datasets. We further analyzed the effect of maximum length of the meta-path and the impacts of different types of features. Moreover, the practical usefulness of our method was confirmed in the predicted novel drug combinations. The source code of EPSDC is available at https://github.com/KDDing/EPSDC.

From gene networks to drugs: systems pharmacology approaches for AUD

The alcohol research field has amassed an impressive number of gene expression datasets spanning key brain areas for addiction, species (humans as well as multiple animal models), and stages in the addiction cycle (binge/intoxication, withdrawal/negative effect, and preoccupation/anticipation). These data have improved our understanding of the molecular adaptations that eventually lead to dysregulation of brain function and the chronic, relapsing disorder of addiction. Identification of new medications to treat alcohol use disorder (AUD) will likely benefit from the integration of genetic, genomic, and behavioral information included in these important datasets. Systems pharmacology considers drug effects as the outcome of the complex network of interactions a drug has rather than a single drug-molecule interaction. Computational strategies based on this principle that integrate gene expression signatures of pharmaceuticals and disease states have shown promise for identifying treatments that ameliorate disease symptoms (called in silico gene mapping or connectivity mapping). In this review, we suggest that gene expression profiling for in silico mapping is critical to improve drug repurposing and discovery for AUD and other psychiatric illnesses. We highlight studies that successfully apply gene mapping computational approaches to identify or repurpose pharmaceutical treatments for psychiatric illnesses. Furthermore, we address important challenges that must be overcome to maximize the potential of these strategies to translate to the clinic and improve healthcare outcomes.

DRUG-NEM: Optimizing drug combinations using single-cell perturbation response to account for intratumoral heterogeneity

Single-cell high-throughput technologies enable the ability to identify combination cancer therapies that account for intratumoral heterogeneity, a phenomenon that has been shown to influence the effectiveness of cancer treatment. We developed and applied an approach that identifies top-ranking drug combinations based on the single-cell perturbation response when an individual tumor sample is screened against a panel of single drugs. This approach optimizes drug combinations by choosing the minimum number of drugs that produce the maximal intracellular desired effects for an individual sample.

An individual malignant tumor is composed of a heterogeneous collection of single cells with distinct molecular and phenotypic features, a phenomenon termed intratumoral heterogeneity. Intratumoral heterogeneity poses challenges for cancer treatment, motivating the need for combination therapies. Single-cell technologies are now available to guide effective drug combinations by accounting for intratumoral heterogeneity through the analysis of the signaling perturbations of an individual tumor sample screened by a drug panel. In particular, Mass Cytometry Time-of-Flight (CyTOF) is a high-throughput single-cell technology that enables the simultaneous measurements of multiple (>40) intracellular and surface markers at the level of single cells for hundreds of thousands of cells in a sample. We developed a computational framework, entitled Drug Nested Effects Models (DRUG-NEM), to analyze CyTOF single-drug perturbation data for the purpose of individualizing drug combinations. DRUG-NEM optimizes drug combinations by choosing the minimum number of drugs that produce the maximal desired intracellular effects based on nested effects modeling. We demonstrate the performance of DRUG-NEM using single-cell drug perturbation data from tumor cell lines and primary leukemia samples.

Developing an Agent-Based Drug Model to Investigate the Synergistic Effects of Drug Combinations

The growth and survival of cancer cells are greatly related to their surrounding microenvironment. To understand the regulation under the impact of anti-cancer drugs and their synergistic effects, we have developed a multiscale agent-based model that can investigate the synergistic effects of drug combinations with three innovations. First, it explores the synergistic effects of drug combinations in a huge dose combinational space at the cell line level. Second, it can simulate the interaction between cells and their microenvironment. Third, it employs both local and global optimization algorithms to train the key parameters and validate the predictive power of the model by using experimental data. The research results indicate that our multicellular system can not only describe the interactions between the microenvironment and cells in detail, but also predict the synergistic effects of drug combinations.

In silico-based screen synergistic drug combinations from herb medicines: a case using Cistanche tubulosa

Neuroinflammation is characterized by the elaborated inflammatory response repertoire of central nervous system tissue. The limitations of the current treatments for neuroinflammation are well-known side effects in the clinical trials of monotherapy. Drug combination therapies are promising strategies to overcome the compensatory mechanisms and off-target effects. However, discovery of synergistic drug combinations from herb medicines is rare. Encouraged by the successfully applied cases we move on to investigate the effective drug combinations based on system pharmacology among compounds from Cistanche tubulosa (SCHENK) R. WIGHT. Firstly, 63 potential bioactive compounds, the related 133 direct and indirect targets are screened out by Drug-likeness evaluation combined with drug targeting process. Secondly, Compound-Target network is built to acquire the data set for predicting drug combinations. We list the top 10 drug combinations which are employed by the algorithm Probability Ensemble Approach (PEA), and Compound-Target-Pathway network is then constructed by the 12 compounds of the combinations, targets, and pathways to unearth the corresponding pharmacological actions. Finally, an integrating pathway approach is developed to elucidate the therapeutic effects of the herb in different pathological features-relevant biological processes. Overall, the method may provide a productive avenue for developing drug combination therapeutics.

Integrative network and transcriptomics-based approach predicts genotype- specific drug combinations for melanoma

Computational methods for drug combination predictions are needed to identify effective therapies that improve durability and prevent drug resistance in an efficient manner. In this paper, we present SynGeNet, a computational method that integrates transcriptomics data characterizing disease and drug z-score profiles with network mining algorithms in order to predict synergistic drug combinations. We compare SynGeNet to other available transcriptomics-based tools to predict drug combinations validated across melanoma cell lines in three genotype groups: BRAF-mutant, NRAS-mutant and combined. We showed that SynGeNet outperforms other available tools in predicting validated drug combinations and single agents tested as part of additional drug pairs. Interestingly, we observed that the performance of SynGeNet decreased when the network construction step was removed and improved when the proportion of matched-genotype validation cell lines increased. These results suggest that delineating functional information from transcriptomics data via network mining and genomic features can improve drug combination predictions.

Identification of optimal strategies for state transition of complex biological networks

Complex biological networks typically contain numerous parameters, and determining feasible strategies for state transition by parameter perturbation is not a trivial task. In the present study, based on dynamical and structural analyses of the biological network, we optimized strategies for controlling variables in a two-node gene regulatory network and a T-cell large granular lymphocyte signaling network associated with blood cancer by using an efficient dynamic optimization method. Optimization revealed the critical value for each decision variable to steer the system from an undesired state into a desired attractor. In addition, the minimum time for the state transition was determined by defining and solving a time-optimal control problem. Moreover, time-dependent variable profiles for state transitions were achieved rather than constant values commonly adopted in previous studies. Furthermore, the optimization method allows multiple controls to be simultaneously adjusted to drive the system out of an undesired attractor. Optimization improved the results of the parameter perturbation method, thus providing a valuable guidance for experimental design.

Toxicity of abamectin and difenoconazole mixtures to a Neotropical cladoceran after simulated run-off and spray drift exposure

Aquatic risk assessments of pesticides in tropical countries have often been disputed for being largely based on risk evaluations conducted in temperate regions. Although pesticide sensitivity comparisons between temperate and tropical freshwater organisms have indeed not revealed consistent differences, risk assessments are currently still based on a relatively small tropical toxicity dataset. In addition, greater levels of runoff and spray drift may be expected in tropical than in temperate agroecosystems, indicating that aquatic life in edge-of-field water bodies is likely to be subjected to higher concentrations of pesticides and their mixtures. The aim of the present study was to evaluate the toxicity of Kraft^® 36 EC (a.i. abamectin), Score^® 250 EC (a.i. difenoconazole) and their mixture to the Neotropical cladoceran Macrothrix flabelligera. Laboratory toxicity tests with the individual formulated products indicated EC50-48h values of 3.1 and 659μg a.i./L given as nominal test concentrations, respectively. Mixtures of the two pesticides revealed a concentration-dependent deviation of the independent action model, with antagonism at low and synergism at high pesticide mixture concentrations. Laboratory toxicity tests were also conducted with microcosm water that was treated with the individual or mixtures through runoff or direct overspray. Microcosm tanks receiving runoff water from experimental soil plots applied with recommended doses of the individual pesticides did not show toxicity to the test organism. Microcosms that received runoff water containing the pesticide mixture, however, did cause a short-term effect on immobility. The microcosms that were treated by direct overspray of both pesticide formulations showed the most pronounced toxic effects. Study findings suggest a potential risk of these pesticides at environmentally relevant concentrations, especially when they are both present.

Informatics Approaches for Predicting, Understanding, and Testing Cancer Drug Combinations

Making cancer treatment more effective is one of the grand challenges in our health care system. However, many drugs have entered clinical trials but so far showed limited efficacy or induced rapid development of resistance. We urgently need multi-targeted drug combinations, which shall selectively inhibit the cancer cells and block the emergence of drug resistance. The book chapter focuses on mathematical and computational tools to facilitate the discovery of the most promising drug combinations to improve efficacy and prevent resistance. Data integration approaches that leverage drug-target interactions, cancer molecular features, and signaling pathways for predicting, understanding, and testing drug combinations are critically reviewed.

Combine and conquer: challenges for targeted therapy combinations in early phase trials

Our increasing understanding of cancer biology has led to the development of molecularly targeted anticancer drugs. The full potential of these agents has not, however, been realised, owing to the presence of de novo (intrinsic) resistance, often resulting from compensatory signalling pathways, or the development of acquired resistance in cancer cells via clonal evolution under the selective pressures of treatment. Combinations of targeted treatments can circumvent some mechanisms of resistance to yield a clinical benefit. We explore the challenges in identifying the best drug combinations and the best combination strategies, as well as the complexities of delivering these treatments to patients. Recognizing treatment-induced toxicity and the inability to use continuous pharmacodynamically effective doses of many targeted treatments necessitates creative intermittent scheduling. Serial tumour profiling and the use of parallel co-clinical trials can contribute to understanding mechanisms of resistance, and will guide the development of adaptive clinical trial designs that can accommodate hypothesis testing, in order to realize the full potential of combination therapies.

ChemEngine: harvesting 3D chemical structures of supplementary data from PDF files

Digital access to chemical journals resulted in a vast array of molecular information that is now available in the supplementary material files in PDF format. However, extracting this molecular information, generally from a PDF document format is a daunting task. Here we present an approach to harvest 3D molecular data from the supporting information of scientific research articles that are normally available from publisher's resources. In order to demonstrate the feasibility of extracting truly computable molecules from PDF file formats in a fast and efficient manner, we have developed a Java based application, namely ChemEngine. This program recognizes textual patterns from the supplementary data and generates standard molecular structure data (bond matrix, atomic coordinates) that can be subjected to a multitude of computational processes automatically. The methodology has been demonstrated via several case studies on different formats of coordinates data stored in supplementary information files, wherein ChemEngine selectively harvested the atomic coordinates and interpreted them as molecules with high accuracy. The reusability of extracted molecular coordinate data was demonstrated by computing Single Point Energies that were in close agreement with the original computed data provided with the articles. It is envisaged that the methodology will enable large scale conversion of molecular information from supplementary files available in the PDF format into a collection of ready- to- compute molecular data to create an automated workflow for advanced computational processes. Software along with source codes and instructions available at https://sourceforge.net/projects/chemengine/files/?source=navbar.Graphical abstract.

Current Trends in Multidrug Optimization: An Alley of Future Successful Treatment of Complex Disorders

The identification of effective and long-lasting cancer therapies still remains elusive, partially due to patient and tumor heterogeneity, acquired drug resistance, and single-drug dose-limiting toxicities. The use of drug combinations may help to overcome some limitations of current cancer therapies by challenging the robustness and redundancy of biological processes. However, effective drug combination optimization requires the careful consideration of numerous parameters. The complexity of this optimization problem is clearly nontrivial and likely requires the assistance of advanced heuristic optimization techniques. In the current review, we discuss the application of optimization techniques for the identification of optimal drug combinations. More specifically, we focus on the application of phenotype-based screening approaches in the field of cancer therapy. These methods are divided into three categories: (1) modeling methods, (2) model-free approaches based on biological search algorithms, and (3) merged approaches, particularly phenotypically driven network biology methods and computation network models relying on phenotypic data. In addition to a brief description of each approach, we include a critical discussion of the advantages and disadvantages of each method, with a strong focus on the limitations and considerations needed to successfully apply such methods in biological research.

Neuroblastoma, a Paradigm for Big Data Science in Pediatric Oncology

Pediatric cancers rarely exhibit recurrent mutational events when compared to most adult cancers. This poses a challenge in understanding how cancers initiate, progress, and metastasize in early childhood. Also, due to limited detected driver mutations, it is difficult to benchmark key genes for drug development. In this review, we use neuroblastoma, a pediatric solid tumor of neural crest origin, as a paradigm for exploring "big data" applications in pediatric oncology. Computational strategies derived from big data science-network- and machine learning-based modeling and drug repositioning-hold the promise of shedding new light on the molecular mechanisms driving neuroblastoma pathogenesis and identifying potential therapeutics to combat this devastating disease. These strategies integrate robust data input, from genomic and transcriptomic studies, clinical data, and in vivo and in vitro experimental models specific to neuroblastoma and other types of cancers that closely mimic its biological characteristics. We discuss contexts in which "big data" and computational approaches, especially network-based modeling, may advance neuroblastoma research, describe currently available data and resources, and propose future models of strategic data collection and analyses for neuroblastoma and other related diseases.

Current Trends in Multidrug Optimization

The identification of effective and long-lasting cancer therapies still remains elusive, partially due to patient and tumor heterogeneity, acquired drug resistance, and single-drug dose-limiting toxicities. The use of drug combinations may help to overcome some limitations of current cancer therapies by challenging the robustness and redundancy of biological processes. However, effective drug combination optimization requires the careful consideration of numerous parameters. The complexity of this optimization problem is clearly nontrivial and likely requires the assistance of advanced heuristic optimization techniques. In the current review, we discuss the application of optimization techniques for the identification of optimal drug combinations. More specifically, we focus on the application of phenotype-based screening approaches in the field of cancer therapy. These methods are divided into three categories: (1) modeling methods, (2) model-free approaches based on biological search algorithms, and (3) merged approaches, particularly phenotypically driven network biology methods and computation network models relying on phenotypic data. In addition to a brief description of each approach, we include a critical discussion of the advantages and disadvantages of each method, with a strong focus on the limitations and considerations needed to successfully apply such methods in biological research.

Looking beyond the cancer cell for effective drug combinations

Combinations of therapies are being actively pursued to expand therapeutic options and deal with cancer’s pervasive resistance to treatment. Research efforts to discover effective combination treatments have focused on drugs targeting intracellular processes of the cancer cells and in particular on small molecules that target aberrant kinases. Accordingly, most of the computational methods used to study, predict, and develop drug combinations concentrate on these modes of action and signaling processes within the cancer cell. This focus on the cancer cell overlooks significant opportunities to tackle other components of tumor biology that may offer greater potential for improving patient survival. Many alternative strategies have been developed to combat cancer; for example, targeting different cancer cellular processes such as epigenetic control; modulating stromal cells that interact with the tumor; strengthening physical barriers that confine tumor growth; boosting the immune system to attack tumor cells; and even regulating the microbiome to support antitumor responses. We suggest that to fully exploit these treatment modalities using effective drug combinations it is necessary to develop multiscale computational approaches that take into account the full complexity underlying the biology of a tumor, its microenvironment, and a patient’s response to the drugs. In this Opinion article, we discuss preliminary work in this area and the needs—in terms of both computational and data requirements—that will truly empower such combinations.

A network-based phenotype mapping approach to identify genes that modulate drug response phenotypes

To better address the problem of drug resistance during cancer chemotherapy and explore the possibility of manipulating drug response phenotypes, we developed a network-based phenotype mapping approach (P-Map) to identify gene candidates that upon perturbed can alter sensitivity to drugs. We used basal transcriptomics data from a panel of human lymphoblastoid cell lines (LCL) to infer drug response networks (DRNs) that are responsible for conferring response phenotypes for anthracycline and taxane, two common anticancer agents use in clinics. We further tested selected gene candidates that interact with phenotypic differentially expressed genes (PDEGs), which are up-regulated genes in LCL for a given class of drug response phenotype in triple-negative breast cancer (TNBC) cells. Our results indicate that it is possible to manipulate a drug response phenotype, from resistant to sensitive or vice versa, by perturbing gene candidates in DRNs and suggest plausible mechanisms regulating directionality of drug response sensitivity. More important, the current work highlights a new way to formulate systems-based therapeutic design: supplementing therapeutics that aim to target disease culprits with phenotypic modulators capable of altering DRN properties with the goal to re-sensitize resistant phenotypes.

Systematic Analysis of Quantitative Logic Model Ensembles Predicts Drug Combination Effects on Cell Signaling Networks

A major challenge in developing anticancer therapies is determining the efficacies of drugs and their combinations in physiologically relevant microenvironments. We describe here our application of "constrained fuzzy logic" (CFL) ensemble modeling of the intracellular signaling network for predicting inhibitor treatments that reduce the phospho-levels of key transcription factors downstream of growth factors and inflammatory cytokines representative of hepatocellular carcinoma (HCC) microenvironments. We observed that the CFL models successfully predicted the effects of several kinase inhibitor combinations. Furthermore, the ensemble predictions revealed ambiguous predictions that could be traced to a specific structural feature of these models, which we resolved with dedicated experiments, finding that IL-1α activates downstream signals through TAK1 and not MEKK1 in HepG2 cells. We conclude that CFL-Q2LM (Querying Quantitative Logic Models) is a promising approach for predicting effective anticancer drug combinations in cancer-relevant microenvironments.

A simple gene set-based method accurately predicts the synergy of drug pairs

Background

The advance in targeted therapy has greatly increased the effectiveness of clinical cancer therapy and reduced the cytotoxicity of treatments to normal cells. However, patients still suffer from cancer relapse due to the occurrence of drug resistance. It is of great need to explore potential combinatorial drug therapy since individual drug alone may not be sufficient to inhibit continuous activation of cancer-addicted genes or pathways. The DREAM challenge has confirmed the potentiality of computational methods for predicting synergistic drug combinations, while the prediction accuracy can be further improved.

Methods

Based on previous reports, we hypothesized the similarity in biological functions or genes perturbed by two drugs can determine their synergistic effects. To test the feasibility of the hypothesis, we proposed three scoring systems: co-gene score, co-GS score, and co-gene/GS score, measuring the similarities in genes with significant expressional changes, enriched gene sets, and significantly changed genes within an enriched gene sets between a pair of drugs, respectively. Performances of these scoring systems were evaluated by the probabilistic c-index (PC-index) devised by the DREAM consortium. We also applied the proposed method to the Connectivity Map dataset to explore more potential synergistic drug combinations.

Results

Using a gold standard derived by the DREAM consortium, we confirmed the prediction power of the three scoring systems (all P-values < 0.05). The co-gene/GS score achieved the best prediction of drug synergy (PC-index = 0.663, P-value < 0.0001), outperforming all methods proposed during DREAM challenge. Furthermore, a binary classification test showed that co-gene/GS scoring was highly accurate and specific. Since our method is constructed on a gene set-based analysis, in addition to synergy prediction, it provides insights into the functional relevance of drug combinations and the underlying mechanisms by which drugs achieve synergy.

Conclusions

Here we proposed a novel and simple method to predict and investigate drug synergy, and validated its efficacy to accurately predict synergistic drug combinations and to comprehensively explore their underlying mechanisms. The method is widely applicable to expression profiles of other drug treatments and is expected to accelerate the realization of precision cancer treatment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0310-3) contains supplementary material, which is available to authorized users.

Novel proteases from the genome of the carnivorous plant Drosera capensis: Structural prediction and comparative analysis

In his 1875 monograph on insectivorous plants, Darwin described the feeding reactions of Drosera flypaper traps and predicted that their secretions contained a "ferment" similar to mammalian pepsin, an aspartic protease. Here we report a high-quality draft genome sequence for the cape sundew, Drosera capensis, the first genome of a carnivorous plant from order Caryophyllales, which also includes the Venus flytrap (Dionaea) and the tropical pitcher plants (Nepenthes). This species was selected in part for its hardiness and ease of cultivation, making it an excellent model organism for further investigations of plant carnivory. Analysis of predicted protein sequences yields genes encoding proteases homologous to those found in other plants, some of which display sequence and structural features that suggest novel functionalities. Because the sequence similarity to proteins of known structure is in most cases too low for traditional homology modeling, 3D structures of representative proteases are predicted using comparative modeling with all-atom refinement. Although the overall folds and active residues for these proteins are conserved, we find structural and sequence differences consistent with a diversity of substrate recognition patterns. Finally, we predict differences in substrate specificities using in silico experiments, providing targets for structure/function studies of novel enzymes with biological and technological significance. Proteins 2016; 84:1517-1533. © 2016 Wiley Periodicals, Inc.

Computational modeling in melanoma for novel drug discovery

INTRODUCTION

There is a growing body of evidence highlighting the applications of computational modeling in the field of biomedicine. It has recently been applied to the in silico analysis of cancer dynamics. In the era of precision medicine, this analysis may allow the discovery of new molecular targets useful for the design of novel therapies and for overcoming resistance to anticancer drugs. According to its molecular behavior, melanoma represents an interesting tumor model in which computational modeling can be applied. Melanoma is an aggressive tumor of the skin with a poor prognosis for patients with advanced disease as it is resistant to current therapeutic approaches.

AREAS COVERED

This review discusses the basics of computational modeling in melanoma drug discovery and development. Discussion includes the in silico discovery of novel molecular drug targets, the optimization of immunotherapies and personalized medicine trials.

EXPERT OPINION

Mathematical and computational models are gradually being used to help understand biomedical data produced by high-throughput analysis. The use of advanced computer models allowing the simulation of complex biological processes provides hypotheses and supports experimental design. The research in fighting aggressive cancers, such as melanoma, is making great strides. Computational models represent the key component to complement these efforts. Due to the combinatorial complexity of new drug discovery, a systematic approach based only on experimentation is not possible. Computational and mathematical models are necessary for bringing cancer drug discovery into the era of omics, big data and personalized medicine.

Transcriptional response networks for elucidating mechanisms of action of multitargeted agents

Drug discovery is moving away from the single target-based approach towards harnessing the potential of polypharmacological agents that modulate the activity of multiple nodes in the complex networks of deregulations underlying disease phenotypes. Computational network pharmacology methods that use systems-level drug-response phenotypes, such as those originating from genome-wide transcriptomic profiles, have proved particularly effective for elucidating the mechanisms of action of multitargeted compounds. Here, we show, via the case study of the natural product pinosylvin, how the combination of two complementary network-based methods can provide novel, unexpected mechanistic insights. This case study also illustrates that elucidating the mechanism of action of multitargeted natural products through transcriptional response-based approaches is a challenging endeavor, often requiring multiple computational-experimental iterations.

Synthetic biology for pharmaceutical drug discovery

Synthetic biology (SB) is an emerging discipline, which is slowly reorienting the field of drug discovery. For thousands of years, living organisms such as plants were the major source of human medicines. The difficulty in resynthesizing natural products, however, often turned pharmaceutical industries away from this rich source for human medicine. More recently, progress on transformation through genetic manipulation of biosynthetic units in microorganisms has opened the possibility of in-depth exploration of the large chemical space of natural products derivatives. Success of SB in drug synthesis culminated with the bioproduction of artemisinin by microorganisms, a tour de force in protein and metabolic engineering. Today, synthetic cells are not only used as biofactories but also used as cell-based screening platforms for both target-based and phenotypic-based approaches. Engineered genetic circuits in synthetic cells are also used to decipher disease mechanisms or drug mechanism of actions and to study cell-cell communication within bacteria consortia. This review presents latest developments of SB in the field of drug discovery, including some challenging issues such as drug resistance and drug toxicity.

Modelling of compound combination effects and applications to efficacy and toxicity: state-of-the-art, challenges and perspectives

The development of treatments involving combinations of drugs is a promising approach towards combating complex or multifactorial disorders. However, the large number of compound combinations that can be generated, even from small compound collections, means that exhaustive experimental testing is infeasible. The ability to predict the behaviour of compound combinations in biological systems, whittling down the number of combinations to be tested, is therefore crucial. Here, we review the current state-of-the-art in the field of compound combination modelling, with the aim to support the development of approaches that, as we hope, will finally lead to an integration of chemical with systems-level biological information for predicting the effect of chemical mixtures.

Using quantitative systems pharmacology for novel drug discovery

INTRODUCTION

Over the past three decades, the predominant paradigm in drug discovery was designing selective ligands for a specific target to avoid unwanted side effects. However, in the last 5 years, the aim has shifted to take into account the biological network in which they interact. Quantitative and Systems Pharmacology (QSP) is a new paradigm that aims to understand how drugs modulate cellular networks in space and time, in order to predict drug targets and their role in human pathophysiology.

AREAS COVERED

This review discusses existing computational and experimental QSP approaches such as polypharmacology techniques combined with systems biology information and considers the use of new tools and ideas in a wider 'systems-level' context in order to design new drugs with improved efficacy and fewer unwanted off-target effects.

EXPERT OPINION

The use of network biology produces valuable information such as new indications for approved drugs, drug-drug interactions, proteins-drug side effects and pathways-gene associations. However, we are still far from the aim of QSP, both because of the huge effort needed to model precisely biological network models and the limited accuracy that we are able to reach with those. Hence, moving from 'one molecule for one target to give one therapeutic effect' to the 'big systems-based picture' seems obvious moving forward although whether our current tools are sufficient for such a step is still under debate.

Synergy Maps: exploring compound combinations using network-based visualization

Background

The phenomenon of super-additivity of biological response to compounds applied jointly, termed synergy, has the potential to provide many therapeutic benefits. Therefore, high throughput screening of compound combinations has recently received a great deal of attention. Large compound libraries and the feasibility of all-pairs screening can easily generate large, information-rich datasets. Previously, these datasets have been visualized using either a heat-map or a network approach—however these visualizations only partially represent the information encoded in the dataset.

Results

A new visualization technique for pairwise combination screening data, termed “Synergy Maps”, is presented. In a Synergy Map, information about the synergistic interactions of compounds is integrated with information about their properties (chemical structure, physicochemical properties, bioactivity profiles) to produce a single visualization. As a result the relationships between compound and combination properties may be investigated simultaneously, and thus may afford insight into the synergy observed in the screen. An interactive web app implementation, available at http://richlewis42.github.io/synergy-maps, has been developed for public use, which may find use in navigating and filtering larger scale combination datasets. This tool is applied to a recent all-pairs dataset of anti-malarials, tested against Plasmodium falciparum, and a preliminary analysis is given as an example, illustrating the disproportionate synergism of histone deacetylase inhibitors previously described in literature, as well as suggesting new hypotheses for future investigation.

Conclusions

Synergy Maps improve the state of the art in compound combination visualization, by simultaneously representing individual compound properties and their interactions. The web-based tool allows straightforward exploration of combination data, and easier identification of correlations between compound properties and interactions.

Identifying kinase dependency in cancer cells by integrating high-throughput drug screening and kinase inhibition data

MOTIVATION

Targeted kinase inhibitors have dramatically improved cancer treatment, but kinase dependency for an individual patient or cancer cell can be challenging to predict. Kinase dependency does not always correspond with gene expression and mutation status. High-throughput drug screens are powerful tools for determining kinase dependency, but drug polypharmacology can make results difficult to interpret.

RESULTS

We developed Kinase Addiction Ranker (KAR), an algorithm that integrates high-throughput drug screening data, comprehensive kinase inhibition data and gene expression profiles to identify kinase dependency in cancer cells. We applied KAR to predict kinase dependency of 21 lung cancer cell lines and 151 leukemia patient samples using published datasets. We experimentally validated KAR predictions of FGFR and MTOR dependence in lung cancer cell line H1581, showing synergistic reduction in proliferation after combining ponatinib and AZD8055.

AVAILABILITY AND IMPLEMENTATION

KAR can be downloaded as a Python function or a MATLAB script along with example inputs and outputs at: http://tanlab.ucdenver.edu/KAR/.

CONTACT

aikchoon.tan@ucdenver.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network-based approaches for drug response prediction and targeted therapy development in cancer

Signaling pathways implicated in cancer create a complex network with numerous regulatory loops and redundant pathways. This complexity explains frequent failure of one-drug-one-target paradigm of treatment, resulting in drug resistance in patients. To overcome the robustness of cell signaling network, cancer treatment should be extended to a combination therapy approach. Integrating and analyzing patient high-throughput data together with the information about biological signaling machinery may help deciphering molecular patterns specific to each patient and finding the best combinations of candidates for therapeutic targeting. We review state of the art in the field of targeted cancer medicine from the computational systems biology perspective. We summarize major signaling network resources and describe their characteristics with respect to applicability for drug response prediction and intervention targets suggestion. Thus discuss methods for prediction of drug sensitivity and intervention combinations using signaling networks together with high-throughput data. Gradual integration of these approaches into clinical routine will improve prediction of response to standard treatments and adjustment of intervention schemes.