High-order combination effects and biological robustness

Integrated omics profiling and network pharmacology uncovers the prognostic genes and multi-targeted therapeutic bioactives to combat lung cancer

Non-small cell lung cancer (NSCLC) is the frequent subtype of lung cancer and the currently used treatment methods, diagnosis, and chemoresistance are relatively ineffective. Determining the pharmacological targets from active biomolecules of medicinal plants has become a frontiers era for biomedical research to develop novel therapies. In view of these scenarios, this pilot study, network pharmacology, cheminformatics, integrative omics, molecular docking and in vitro anti-cancer analysis were performed to unveil the multi-targeted treatment mechanisms of novel plant bioactives to treat lung cancer. Bioactive molecules from medicinal plants were compiled from PubChem. Network pharmacology approach revealed that 29 compounds efficiently target the 390 human and lung cancer associated genes. In addition, comparative analysis was performed and identified the 7 bioactive molecules significantly targeting 8 lung cancer genes. The integrative omics analysis discovered unique genes between the lung cancer and normal lung tissues. These genes were further validated through protein-protein interaction, gene ontology, gene functional and pathway enrichment, boxplot and overall survival analyses to understand the function of unique genes and their involvement in cancer signaling pathways. Survival heatmap analyses identified the significant prognostic genes. Docking results revealed that, lupeol and p-coumaric acid displayed high binding affinities with MIF, CCNB1, FABP4. Hence, we selected these two bioactives for in vitro analysis. Furthermore, these selected bioactives were showed concentration dependent cytotoxicity against the lung adenocarcinoma cells (A549). This holistic study has opened up novel avenues and unravels the cancer prognostic genes which could serve as druggable target and bioactives with anti-cancerous efficacy. Further functional validations are prerequisites to deciphering these bioactives as commercial drug candidates.

Intelligent host engineering for metabolic flux optimisation in biotechnology

Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.

ASN007 is a selective ERK1/2 inhibitor with preferential activity against RAS-and RAF-mutant tumors

Inhibition of the extracellular signal-regulated kinases ERK1 and ERK2 (ERK1/2) offers a promising therapeutic strategy in cancers harboring activated RAS/RAF/MEK/ERK signaling pathways. Here, we describe an orally bioavailable and selective ERK1/2 inhibitor, ASN007, currently in clinical development for the treatment of cancer. In preclinical studies, ASN007 shows strong antiproliferative activity in tumors harboring mutations in BRAF and RAS (KRAS, NRAS, and HRAS). ASN007 demonstrates activity in a BRAFV600E mutant melanoma tumor model that is resistant to BRAF and MEK inhibitors. The PI3K inhibitor copanlisib enhances the antiproliferative activity of ASN007 both in vitro and in vivo due to dual inhibition of RAS/MAPK and PI3K survival pathways. Our data provide a rationale for evaluating ASN007 in RAS/RAF-driven tumors as well as a mechanistic basis for combining ASN007 with PI3K inhibitors.

Genetic interaction analysis in microbial pathogens: unravelling networks of pathogenesis, antimicrobial susceptibility and host interactions

Genetic interaction (GI) analysis is a powerful genetic strategy that analyzes the fitness and phenotypes of single- and double-gene mutant cells in order to dissect the epistatic interactions between genes, categorize genes into biological pathways, and characterize genes of unknown function. GI analysis has been extensively employed in model organisms for foundational, systems-level assessment of the epistatic interactions between genes. More recently, GI analysis has been applied to microbial pathogens and has been instrumental for the study of clinically important infectious organisms. Here, we review recent advances in systems-level GI analysis of diverse microbial pathogens, including bacterial and fungal species. We focus on important applications of GI analysis across pathogens, including GI analysis as a means to decipher complex genetic networks regulating microbial virulence, antimicrobial drug resistance and host-pathogen dynamics, and GI analysis as an approach to uncover novel targets for combination antimicrobial therapeutics. Together, this review bridges our understanding of GI analysis and complex genetic networks, with applications to diverse microbial pathogens, to further our understanding of virulence, the use of antimicrobial therapeutics and host-pathogen interactions. .

SPADIS: An Algorithm for Selecting Predictive and Diverse SNPs in GWAS

Phenotypic heritability of complex traits and diseases is seldom explained by individual genetic variants identified in genome-wide association studies (GWAS). Many methods have been developed to select a subset of variant loci, which are associated with or predictive of the phenotype. Selecting connected SNPs on SNP-SNP networks have been proven successful in finding biologically interpretable and predictive SNPs. However, we argue that the connectedness constraint favors selecting redundant features that affect similar biological processes and therefore does not necessarily yield better predictive performance. In this paper, we propose a novel method called SPADIS that favors the selection of remotely located SNPs in order to account for their complementary effects in explaining a phenotype. SPADIS selects a diverse set of loci on a SNP-SNP network. This is achieved by maximizing a submodular set function with a greedy algorithm that ensures a constant factor approximation to the optimal solution. We compare SPADIS to the state-of-the-art method SConES, on a dataset of Arabidopsis Thaliana with continuous flowering time phenotypes. SPADIS has better average phenotype prediction performance in 15 out of 17 phenotypes when the same number of SNPs are selected and provides consistent improvements across multiple networks and settings on average. Moreover, it identifies more candidate genes and runs faster.

Multi-parametric characterization of drug effects on cells

We present here a novel multi-parametric approach for the characterization of multiple cellular features, using images acquired by high-throughput and high-definition light microscopy. We specifically used this approach for deep and unbiased analysis of the effects of a drug library on five cultured cell lines. The presented method enables the acquisition and analysis of millions of images, of treated and control cells, followed by an automated identification of drugs inducing strong responses, evaluating the median effect concentrations and those cellular properties that are most highly affected by the drug. The tools described here provide standardized quantification of multiple attributes for systems level dissection of complex functions in normal and diseased cells, using multiple perturbations. Such analysis of cells, derived from pathological samples, may help in the diagnosis and follow-up of treatment in patients.

In silico authentication of amygdalin as a potent anticancer compound in the bitter kernels of family Rosaceae

Amygdalin a naturally occurring compound, predominantly in the bitter kernels of apricot, almond, apple and other members of Rosaceae family. Though, amygdalin is used as an alternative therapy to treat various types of cancer but its role in cancer pathways has rarely been explored yet. Therefore, present study was intended with the aim to investigate the alleged anti-cancerous effects of amygdalin specifically on PI3K-AKT-mTOR and Ras pathways of cancer in human body. Computational modelling and simulation techniques were used to assess the effect of amygdalin on PI3K-AKT-mTOR and Ras pathways using different level of dosage. It was observed that amygdalin had direct and substantial contribution to regulate PI3K-mTOR activities on threshold levels while the other caner pathways were effected indirectly. Consequently, amygdalin is a down-regulator of a cancer within a specified amount and contribute considerably to reduce various types of cancer in human. Furthermore, in-vitro and in-vivo analyses of amygdalin could be of helpful to authenticate its pharmacological effects.

Antimicrobial and Antioxidant Properties of Four Lycopus Taxa and an Interaction Study of Their Major Compounds

The compositions of leaf infusions of three genotypes of Lycopus europaeus L. with origins in central Europe, namely L. europaeus A (LeuA), L. europaeus B (LeuB), and L. europaeus C (LeuC), and one genotype of L. exaltatus (Lex), were examined by LC-MS-DAD (Liquid Chromatography Mass Spectrometry and Diode Array Detection) analysis. This revealed the presence of thirteen compounds belonging to the groups of phenolic acids and flavonoids, with a predominance of rosmarinic acid (RA) and luteolin-7-O-glucuronide (LGlr). The antimicrobial activity of leaf infusions was tested on the collection strains of Gram-positive and Gram-negative bacteria, and on the clinical Staphylococcus aureus strains. We detected higher activity against Gram-positive bacteria, of which the most susceptible strains were those of Staphylococcus aureus, including methicillin-resistant and poly-resistant strains. Furthermore, we examined the antioxidant activity of leaf infusions using 2,2-diphenyl-1-picrylhydrazyl (DPPH) and 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) methods, and on NIH/3T3 cell lines using dichlorodihydrofluorescein diacetate (DCFH-DA). We also studied the mutual interactions between selected infusions, namely RA and/or LGlr. In the mixtures of leaf infusion and RA or LGlr, we observed slight synergism and a high dose reduction index in most cases. This leads to the beneficial dose reduction at a given antioxidant effect level in mixtures compared to the doses of the parts used alone. Therefore, our study draws attention to further applications of the Lycopus leaves as a valuable alternative source of natural antioxidants and as a promising topical antibacterial agent for medicinal use.

Identification of a Synergistic Multi-Drug Combination Active in Cancer Cells via the Prevention of Spindle Pole Clustering

A major limitation of clinically used cancer drugs is the lack of specificity resulting in toxicity. To address this, we performed a phenotypically-driven screen to identify optimal multidrug combinations acting with high efficacy and selectivity in clear cell renal cell carcinoma (ccRCC). The search was performed using the Therapeutically Guided Multidrug Optimization (TGMO) method in ccRCC cells (786-O) and nonmalignant renal cells and identified a synergistic low-dose four-drug combination (C2) with high efficacy and negligible toxicity. We discovered that C2 inhibits multipolar spindle pole clustering, a survival mechanism employed by cancer cells with spindle abnormalities. This phenotype was also observed in 786-O cells resistant to sunitinib, the first line ccRCC treatment, as well as in melanoma cells with distinct percentages of supernumerary centrosomes. We conclude that C2-treatment shows a high efficacy in cells prone to form multipolar spindles. Our data suggest a highly effective and selective C2 treatment strategy for malignant and drug-resistant cancers.

Response envelope analysis for quantitative evaluation of drug combinations

MOTIVATION

The concept of synergy between two agents, over a century old, is important to the fields of biology, chemistry, pharmacology and medicine. A key step in drug combination analysis is the selection of an additivity model to identify combination effects including synergy, additivity and antagonism. Existing methods for identifying and interpreting those combination effects have limitations.

RESULTS

We present here a computational framework, termed response envelope analysis (REA), that makes use of 3D response surfaces formed by generalized Loewe Additivity and Bliss Independence models of interaction to evaluate drug combination effects. Because the two models imply two extreme limits of drug interaction (mutually exclusive and mutually non-exclusive), a response envelope defined by them provides a quantitatively stringent additivity model for identifying combination effects without knowing the inhibition mechanism. As a demonstration, we apply REA to representative published data from large screens of anticancer and antibiotic combinations. We show that REA is more accurate than existing methods and provides more consistent results in the context of cross-experiment evaluation.

AVAILABILITY AND IMPLEMENTATION

The open-source software package associated with REA is available at: https://github.com/4dsoftware/rea.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In silico perturbation of drug targets in pan-cancer analysis combining multiple networks and pathways

The knowledge of cancer cell response to conventional therapies is crucial in order to choose the correct therapy of patients affected by cancer. The major problem is generally attributed to the lack of specific biological processes able to predict the therapy efficacy. Here, we optimized a computational method for the analysis of gene networks able to detect and quantify the effects of a drug in a pan-cancer study. Overall, our method, using several network topological measures has identified a cancer gene network with a key role in biological processes. The gene network, able to classify with a good performance cancer vs normal samples, was modulated in silico to evaluate the effects of new or approved drugs. This computational model could offer an interesting hint to decipher molecular mechanisms contributing to resistance or inefficacy of drugs.

Diagonal Method to Measure Synergy Among Any Number of Drugs

A synergistic drug combination has a higher efficacy compared to the effects of individual drugs. Checkerboard assays, where drugs are combined in many doses, allow sensitive measurement of drug interactions. However, these assays are costly and do not scale well for measuring interaction among many drugs. Several recent studies have reported drug interaction measurements using a diagonal sampling of the traditional checkerboard assay. This alternative methodology greatly decreases the cost of drug interaction experiments and allows interaction measurement for combinations with many drugs. Here, we describe a protocol to measure the three pairwise interactions and one three-way interaction among three antibiotics in duplicate, in five days, using only three 96-well microplates and standard laboratory equipment. We present representative results showing that the three-antibiotic combination of Levofloxacin + Nalidixic Acid + Penicillin G is synergistic. Our protocol scales up to measure interactions among many drugs and in other biological contexts, allowing for efficient screens for multi-drug synergies against pathogens and tumors.

Evaluation of drug combination effect using a Bliss independence dose-response surface model

To test the anticancer effect of combining two drugs targeting different biological pathways, the popular way to show synergistic effect of drug combination is a heat map or surface plot based on the percent excess the Bliss prediction using the average response measures at each combination dose. Such graphs, however, are inefficient in the drug screening process and it doesn't give a statistical inference on synergistic effect. To make a statistically rigorous and robust conclusion for drug combination effect, we present a two-stage Bliss independence response surface model to estimate an overall interaction index (τ) with 95% confidence interval (CI). By taking into all data points account, the overall τ with 95% CI can be applied to determine if the drug combination effect is synergistic overall. Using some example data, the two-stage model was compared to a couple of classic models following Bliss rule. The data analysis results obtained from our model reflect the pattern shown from other models. The application of overall τ helps investigators to make decision easier and accelerate the preclinical drug screening.

Modeling of the jasmonate signaling pathway in Arabidopsis thaliana with respect to pathophysiology of Alternaria blight in Brassica

The productivity of Oilseed Brassica, one of the economically important crops of India, is seriously affected by the disease, Alternaria blight. The disease is mainly caused by two major necrotrophic fungi, Alternaria brassicae and Alternaria brassicicola which are responsible for significant yield losses. Till date, no resistant source is available against Alternaria blight, hence plant breeding methods can not be used to develop disease resistant varieties. Jasmonate mediated signalling pathway, which is known to play crucial role during defense response against necrotrophs, could be strengthened in Brassica plants to combat the disease. Since scanty information is available in Brassica-Alternaria pathosystems at molecular level therefore, in the present study efforts have been made to model jasmonic acid pathway in Arabidopsis thaliana to simulate the dynamic behaviour of molecular species in the model. Besides, the developed model was also analyzed topologically for investigation of the hubs node. COI1 is identified as one of the promising candidate genes in response to Alternaria and other linked components of plant defense mechanisms against the pathogens. The findings from present study are therefore informative for understanding the molecular basis of pathophysiology and rational management of Alternaria blight for securing food and nutritional security.

Niche harmony search algorithm for detecting complex disease associated high-order SNP combinations

Genome-wide association study is especially challenging in detecting high-order disease-causing models due to model diversity, possible low or even no marginal effect of the model, and extraordinary search and computations. In this paper, we propose a niche harmony search algorithm where joint entropy is utilized as a heuristic factor to guide the search for low or no marginal effect model, and two computationally lightweight scores are selected to evaluate and adapt to diverse of disease models. In order to obtain all possible suspected pathogenic models, niche technique merges with HS, which serves as a taboo region to avoid HS trapping into local search. From the resultant set of candidate SNP-combinations, we use G-test statistic for testing true positives. Experiments were performed on twenty typical simulation datasets in which 12 models are with marginal effect and eight ones are with no marginal effect. Our results indicate that the proposed algorithm has very high detection power for searching suspected disease models in the first stage and it is superior to some typical existing approaches in both detection power and CPU runtime for all these datasets. Application to age-related macular degeneration (AMD) demonstrates our method is promising in detecting high-order disease-causing models.

Enhancing Drug Efficacy and Therapeutic Index through Cheminformatics-Based Selection of Small Molecule Binary Weapons That Improve Transporter-Mediated Targeting: A Cytotoxicity System Based on Gemcitabine

The transport of drug molecules is mainly determined by the distribution of influx and efflux transporters for which they are substrates. To enable tissue targeting, we sought to develop the idea that we might affect the transporter-mediated disposition of small-molecule drugs via the addition of a second small molecule that of itself had no inhibitory pharmacological effect but that influenced the expression of transporters for the primary drug. We refer to this as a “binary weapon” strategy. The experimental system tested the ability of a molecule that on its own had no cytotoxic effect to increase the toxicity of the nucleoside analog gemcitabine to Panc1 pancreatic cancer cells. An initial phenotypic screen of a 500-member polar drug (fragment) library yielded three “hits.” The structures of 20 of the other 2,000 members of this library suite had a Tanimoto similarity greater than 0.7 to those of the initial hits, and each was itself a hit (the cheminformatics thus providing for a massive enrichment). We chose the top six representatives for further study. They fell into three clusters whose members bore reasonable structural similarities to each other (two were in fact isomers), lending strength to the self-consistency of both our conceptual and experimental strategies. Existing literature had suggested that indole-3-carbinol might play a similar role to that of our fragments, but in our hands it was without effect; nor was it structurally similar to any of our hits. As there was no evidence that the fragments could affect toxicity directly, we looked for effects on transporter transcript levels. In our hands, only the ENT1-3 uptake and ABCC2,3,4,5, and 10 efflux transporters displayed measurable transcripts in Panc1 cultures, along with a ribonucleoside reductase RRM1 known to affect gemcitabine toxicity. Very strikingly, the addition of gemcitabine alone increased the expression of the transcript for ABCC2 (MRP2) by more than 12-fold, and that of RRM1 by more than fourfold, and each of the fragment “hits” served to reverse this. However, an inhibitor of ABCC2 was without significant effect, implying that RRM1 was possibly the more significant player. These effects were somewhat selective for Panc cells. It seems, therefore, that while the effects we measured were here mediated more by efflux than influx transporters, and potentially by other means, the binary weapon idea is hereby fully confirmed: it is indeed possible to find molecules that manipulate the expression of transporters that are involved in the bioactivity of a pharmaceutical drug. This opens up an entirely new area, that of chemical genomics-based drug targeting.

Detecting High-Order Epistasis in Nonlinear Genotype-Phenotype Maps

High-order epistasis has been observed in many genotype-phenotype maps. These multi-way interactions between mutations may be useful for dissecting complex traits and could have profound implications for evolution. Alternatively, they could be a statistical artifact. High-order epistasis models assume the effects of mutations should add, when they could in fact multiply or combine in some other nonlinear way. A mismatch in the “scale” of the epistasis model and the scale of the underlying map would lead to spurious epistasis. In this article, we develop an approach to estimate the nonlinear scales of arbitrary genotype-phenotype maps. We can then linearize these maps and extract high-order epistasis. We investigated seven experimental genotype-phenotype maps for which high-order epistasis had been reported previously. We find that five of the seven maps exhibited nonlinear scales. Interestingly, even after accounting for nonlinearity, we found statistically significant high-order epistasis in all seven maps. The contributions of high-order epistasis to the total variation ranged from 2.2 to 31.0%, with an average across maps of 12.7%. Our results provide strong evidence for extensive high-order epistasis, even after nonlinear scale is taken into account. Further, we describe a simple method to estimate and account for nonlinearity in genotype-phenotype maps.

Dissecting Cell-Fate Determination Through Integrated Mathematical Modeling of the ERK/MAPK Signaling Pathway

The past three decades have witnessed an enormous progress in the elucidation of the ERK/MAPK signaling pathway and its involvement in various cellular processes. Because of its importance and complex wiring, the ERK pathway has been an intensive subject for mathematical modeling, which facilitates the unraveling of key dynamic properties and behaviors of the pathway. Recently, however, it became evident that the pathway does not act in isolation but closely interacts with many other pathways to coordinate various cellular outcomes under different pathophysiological contexts. This has led to an increasing number of integrated, large-scale models that link the ERK pathway to other functionally important pathways. In this chapter, we first discuss the essential steps in model development and notable models of the ERK pathway. We then use three examples of integrated, multipathway models to investigate how crosstalk of ERK signaling with other pathways regulates cell-fate decision-making in various physiological and disease contexts. Specifically, we focus on ERK interactions with the phosphoinositide-3 kinase (PI3K), c-Jun N-terminal kinase (JNK), and β-adrenergic receptor (β-AR) signaling pathways. We conclude that integrated modeling in combination with wet-lab experimentation have been and will be instrumental in gaining an in-depth understanding of ERK signaling in multiple biological contexts.

Prediction of multidimensional drug dose responses based on measurements of drug pairs

Finding potent multidrug combinations against cancer and infections is a pressing therapeutic challenge; however, screening all combinations is difficult because the number of experiments grows exponentially with the number of drugs and doses. To address this, we present a mathematical model that predicts the effects of three or more antibiotics or anticancer drugs at all doses based only on measurements of drug pairs at a few doses, without need for mechanistic information. The model provides accurate predictions on available data for antibiotic combinations, and on experiments presented here on the response matrix of three cancer drugs at eight doses per drug. This approach offers a way to search for effective multidrug combinations using a small number of experiments.

Cancer network activity associated with therapeutic response and synergism

Background

Cancer patients often show no or only modest benefit from a given therapy. This major problem in oncology is generally attributed to the lack of specific predictive biomarkers, yet a global measure of cancer cell activity may support a comprehensive mechanistic understanding of therapy efficacy. We reasoned that network analysis of omic data could help to achieve this goal.

Methods

A measure of “cancer network activity” (CNA) was implemented based on a previously defined network feature of communicability. The network nodes and edges corresponded to human proteins and experimentally identified interactions, respectively. The edges were weighted proportionally to the expression of the genes encoding for the corresponding proteins and relative to the number of direct interactors. The gene expression data corresponded to the basal conditions of 595 human cancer cell lines. Therapeutic responses corresponded to the impairment of cell viability measured by the half maximal inhibitory concentration (IC₅₀) of 130 drugs approved or under clinical development. Gene ontology, signaling pathway, and transcription factor-binding annotations were taken from public repositories. Predicted synergies were assessed by determining the viability of four breast cancer cell lines and by applying two different analytical methods.

Results

The effects of drug classes were associated with CNAs formed by different cell lines. CNAs also differentiate target families and effector pathways. Proteins that occupy a central position in the network largely contribute to CNA. Known key cancer-associated biological processes, signaling pathways, and master regulators also contribute to CNA. Moreover, the major cancer drivers frequently mediate CNA and therapeutic differences. Cell-based assays centered on these differences and using uncorrelated drug effects reveals novel synergistic combinations for the treatment of breast cancer dependent on PI3K-mTOR signaling.

Conclusions

Cancer therapeutic responses can be predicted on the basis of a systems-level analysis of molecular interactions and gene expression. Fundamental cancer processes, pathways, and drivers contribute to this feature, which can also be exploited to predict precise synergistic drug combinations.

Electronic supplementary material

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

NLLSS: Predicting Synergistic Drug Combinations Based on Semi-supervised Learning

Fungal infection has become one of the leading causes of hospital-acquired infections with high mortality rates. Furthermore, drug resistance is common for fungus-causing diseases. Synergistic drug combinations could provide an effective strategy to overcome drug resistance. Meanwhile, synergistic drug combinations can increase treatment efficacy and decrease drug dosage to avoid toxicity. Therefore, computational prediction of synergistic drug combinations for fungus-causing diseases becomes attractive. In this study, we proposed similar nature of drug combinations: principal drugs which obtain synergistic effect with similar adjuvant drugs are often similar and vice versa. Furthermore, we developed a novel algorithm termed Network-based Laplacian regularized Least Square Synergistic drug combination prediction (NLLSS) to predict potential synergistic drug combinations by integrating different kinds of information such as known synergistic drug combinations, drug-target interactions, and drug chemical structures. We applied NLLSS to predict antifungal synergistic drug combinations and showed that it achieved excellent performance both in terms of cross validation and independent prediction. Finally, we performed biological experiments for fungal pathogen Candida albicans to confirm 7 out of 13 predicted antifungal synergistic drug combinations. NLLSS provides an efficient strategy to identify potential synergistic antifungal combinations.

Drug combinations represent a promising strategy for overcoming fungal drug resistance and treating complex diseases. There is an urgent need to establish powerful computational methods for systematic prediction of synergistic drug combination on a large scale. Based on the assumption that principal drugs which obtain synergistic effect with similar adjuvant drugs are often similar and vice versa, NLLSS was developed to predict potential synergistic drug combinations by integrating known synergistic drug combinations, unlabeled drug combinations, drug-target interactions, and drug chemical structures. NLLSS has obtained the reliable performance in the cross validation and experimental validations, which indicated that NLLSS has an excellent performance of identifying potential synergistic drug combinations. Out of 13 predicted antifungal synergistic drug combinations, 7 candidates were experimentally confirmed. It is anticipated that NLLSS would be an important and useful resource by providing a new strategy to identify potential synergistic antifungal combinations, explore new indications of existing drugs, and provide useful insights into the underlying molecular mechanisms of synergistic drug combinations.

Systems Pharmacology and Rational Polypharmacy: Nitric Oxide-Cyclic GMP Signaling Pathway as an Illustrative Example and Derivation of the General Case

Impaired nitric oxide (NO˙)-cyclic guanosine 3', 5'-monophosphate (cGMP) signaling has been observed in many cardiovascular disorders, including heart failure and pulmonary arterial hypertension. There are several enzymatic determinants of cGMP levels in this pathway, including soluble guanylyl cyclase (sGC) itself, the NO˙-activated form of sGC, and phosphodiesterase(s) (PDE). Therapies for some of these disorders with PDE inhibitors have been successful at increasing cGMP levels in both cardiac and vascular tissues. However, at the systems level, it is not clear whether perturbation of PDE alone, under oxidative stress, is the best approach for increasing cGMP levels as compared with perturbation of other potential pathway targets, either alone or in combination. Here, we develop a model-based approach to perturbing this pathway, focusing on single reactions, pairs of reactions, or trios of reactions as targets, then monitoring the theoretical effects of these interventions on cGMP levels. Single perturbations of all reaction steps within this pathway demonstrated that three reaction steps, including the oxidation of sGC, NO˙ dissociation from sGC, and cGMP degradation by PDE, exerted a dominant influence on cGMP accumulation relative to other reaction steps. Furthermore, among all possible single, paired, and triple perturbations of this pathway, the combined perturbations of these three reaction steps had the greatest impact on cGMP accumulation. These computational findings were confirmed in cell-based experiments. We conclude that a combined perturbation of the oxidatively-impaired NO˙-cGMP signaling pathway is a better approach to the restoration of cGMP levels as compared with corresponding individual perturbations. This approach may also yield improved therapeutic responses in other complex pharmacologically amenable pathways.

Developing drugs for a well-defined biochemical or molecular pathway has conventionally been approached by optimizing the inhibition (or activation) of a single target by a single pharmacologic agent. On occasion, drug combinations have been used that generally target multiple pathways affecting a common phenotype, again by optimizing the extent of inhibition of individual targets, semi-empirically adjusting their doses to minimize toxicities as they are manifest. Here, we present a computational approach for identifying optimal combinations of agents that can affect (inhibit) a well-defined biochemical pathway, doing so at minimal combined concentrations, thereby potentially minimizing dose-dependent toxicities. This approach is illustrated computationally and experimentally with a well-known pathway, the nitric oxide-cyclic GMP pathway, but is readily generalizable to rational polypharmacy.

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.

Unlocking the Potential of High-Throughput Drug Combination Assays Using Acoustic Dispensing

Assessment of synergistic effects of drug combinations in vitro is a critical part of anticancer drug research. However, the complexities of dosing and analyzing two drugs over the appropriate range of doses have generally led to compromises in experimental design that restrict the quality and robustness of the data. In particular, the use of a single dose response of combined drugs, rather than a full two-way matrix of varying doses, has predominated in higher-throughput studies. Acoustic dispensing unlocks the potential of high-throughput dose matrix analysis. We have developed acoustic dispensing protocols that enable compound synergy assays in a 384-well format. This experimental design is considerably more efficient and flexible with respect to time, reagent usage, and labware than is achievable using traditional serial-dilution approaches. Data analysis tools integrated in Genedata Screener were used to efficiently deconvolute the combination compound mapping scheme and calculate compound potency and synergy metrics. We have applied this workflow to evaluate interactions among drugs targeting different nodes of the mitogen-activated protein kinase pathway in a panel of cancer cell lines.

Searching High-Order SNP Combinations for Complex Diseases Based on Energy Distribution Difference

Single nucleotide polymorphisms, a dominant type of genetic variants, have been used successfully to identify defective genes causing human single gene diseases. However, most common human diseases are complex diseases and caused by gene-gene and gene-environment interactions. Many SNP-SNP interaction analysis methods have been introduced but they are not powerful enough to discover interactions more than three SNPs. The paper proposes a novel method that analyzes all SNPs simultaneously. Different from existing methods, the method regards an individual's genotype data on a list of SNPs as a point with a unit of energy in a multi-dimensional space, and tries to find a new coordinate system where the energy distribution difference between cases and controls reaches the maximum. The method will find different multiple SNPs combinatorial patterns between cases and controls based on the new coordinate system. The experiment on simulated data shows that the method is efficient. The tests on the real data of age-related macular degeneration (AMD) disease show that it can find out more significant multi-SNP combinatorial patterns than existing methods.

Network pharmacology strategies toward multi-target anticancer therapies: from computational models to experimental design principles

Polypharmacology has emerged as novel means in drug discovery for improving treatment response in clinical use. However, to really capitalize on the polypharmacological effects of drugs, there is a critical need to better model and understand how the complex interactions between drugs and their cellular targets contribute to drug efficacy and possible side effects. Network graphs provide a convenient modeling framework for dealing with the fact that most drugs act on cellular systems through targeting multiple proteins both through on-target and off-target binding. Network pharmacology models aim at addressing questions such as how and where in the disease network should one target to inhibit disease phenotypes, such as cancer growth, ideally leading to therapies that are less vulnerable to drug resistance and side effects by means of attacking the disease network at the systems level through synergistic and synthetic lethal interactions. Since the exponentially increasing number of potential drug target combinations makes pure experimental approach quickly unfeasible, this review depicts a number of computational models and algorithms that can effectively reduce the search space for determining the most promising combinations for experimental evaluation. Such computational-experimental strategies are geared toward realizing the full potential of multi-target treatments in different disease phenotypes. Our specific focus is on system-level network approaches to polypharmacology designs in anticancer drug discovery, where we give representative examples of how network-centric modeling may offer systematic strategies toward better understanding and even predicting the phenotypic responses to multi-target therapies.

Searching for synergies: matrix algebraic approaches for efficient pair screening

Functionally interacting perturbations, such as synergistic drugs pairs or synthetic lethal gene pairs, are of key interest in both pharmacology and functional genomics. However, to find such pairs by traditional screening methods is both time consuming and costly. We present a novel computational-experimental framework for efficient identification of synergistic target pairs, applicable for screening of systems with sizes on the order of current drug, small RNA or SGA (Synthetic Genetic Array) libraries (>1000 targets). This framework exploits the fact that the response of a drug pair in a given system, or a pair of genes' propensity to interact functionally, can be partly predicted by computational means from (i) a small set of experimentally determined target pairs, and (ii) pre-existing data (e.g. gene ontology, PPI) on the similarities between targets. Predictions are obtained by a novel matrix algebraic technique, based on cyclical projections onto convex sets. We demonstrate the efficiency of the proposed method using drug-drug interaction data from seven cancer cell lines and gene-gene interaction data from yeast SGA screens. Our protocol increases the rate of synergism discovery significantly over traditional screening, by up to 7-fold. Our method is easy to implement and could be applied to accelerate pair screening for both animal and microbial systems.

Optimization and model reduction in the high dimensional parameter space of a budding yeast cell cycle model

Background

Parameter estimation from experimental data is critical for mathematical modeling of protein regulatory networks. For realistic networks with dozens of species and reactions, parameter estimation is an especially challenging task. In this study, we present an approach for parameter estimation that is effective in fitting a model of the budding yeast cell cycle (comprising 26 nonlinear ordinary differential equations containing 126 rate constants) to the experimentally observed phenotypes (viable or inviable) of 119 genetic strains carrying mutations of cell cycle genes.

Results

Starting from an initial guess of the parameter values, which correctly captures the phenotypes of only 72 genetic strains, our parameter estimation algorithm quickly improves the success rate of the model to 105–111 of the 119 strains. This success rate is comparable to the best values achieved by a skilled modeler manually choosing parameters over many weeks. The algorithm combines two search and optimization strategies. First, we use Latin hypercube sampling to explore a region surrounding the initial guess. From these samples, we choose ∼20 different sets of parameter values that correctly capture wild type viability. These sets form the starting generation of differential evolution that selects new parameter values that perform better in terms of their success rate in capturing phenotypes. In addition to producing highly successful combinations of parameter values, we analyze the results to determine the parameters that are most critical for matching experimental outcomes and the most competitive strains whose correct outcome with a given parameter vector forces numerous other strains to have incorrect outcomes. These “most critical parameters” and “most competitive strains” provide biological insights into the model. Conversely, the “least critical parameters” and “least competitive strains” suggest ways to reduce the computational complexity of the optimization.

Conclusions

Our approach proves to be a useful tool to help systems biologists fit complex dynamical models to large experimental datasets. In the process of fitting the model to the data, the tool identifies suggestive correlations among aspects of the model and the data.

Navigating personalized medicine dependent on modular flexibility

Deconstructing networks and rewiring alterable modules in a rational way is critical to optimize drug discovery and develop personalized medicine.

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

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

An Integrative Platform of TCM Network Pharmacology and Its Application on a Herbal Formula, Qing-Luo-Yin

The scientific understanding of traditional Chinese medicine (TCM) has been hindered by the lack of methods that can explore the complex nature and combinatorial rules of herbal formulae. On the assumption that herbal ingredients mainly target a molecular network to adjust the imbalance of human body, here we present a-self-developed TCM network pharmacology platform for discovering herbal formulae in a systematic manner. This platform integrates a set of network-based methods that we established previously to catch the network regulation mechanism and to identify active ingredients as well as synergistic combinations for a given herbal formula. We then provided a case study on an antirheumatoid arthritis (RA) formula, Qing-Luo-Yin (QLY), to demonstrate the usability of the platform. We revealed the target network of QLY against RA-related key processes including angiogenesis, inflammatory response, and immune response, based on which we not only predicted active and synergistic ingredients from QLY but also interpreted the combinatorial rule of this formula. These findings are either verified by the literature evidence or have the potential to guide further experiments. Therefore, such a network pharmacology strategy and platform is expected to make the systematical study of herbal formulae achievable and to make the TCM drug discovery predictable.

Insights into the regenerative property of plant cells and their receptivity to transgenesis: wheat as a research case study

From a holistic perspective, the discovery of cellular plasticity, a very interesting property of totipotency, underlies many topical issues in biology with important medical applications, while transgenesis is a core research tool in biology. Partially known, some basic mechanisms involved in the regenerative property of cells and in their receptivity to transgenesis are common to plant and animal cells and highlight the principle of the unity of life. Transgenesis provides an important investigative instrument in plant physiology and is regarded as a valuable tool for crop improvement. The economic, social, cultural and scientific importance of cereals has led to a rich stream of research into their genetics, biology and evolution. Sustained efforts to achieve the results obtained in the fields of genetic engineering and applied biotechnology reflect this deep interest. Difficulties encountered in creating genetically modified cereals, especially wheat, highlighted the central notions of tissue culture regeneration and transformation competencies. From the perspective of combining or encountering these competencies in the same cell lineage, this reputedly recalcitrant species provides a stimulating biological system in which to explore the physiological and genetic complexity of both competencies. The former involves two phases, dedifferentiation and redifferentiation. Cells undergo development switches regulated by extrinsic and intrinsic factors. The re-entry into the cell division cycle progressively culminates in the development of organized structures. This is achieved by global chromatin reorganization associated with the reprogramming of the gene expression pattern. The latter is linked with surveillance mechanisms and DNA repair, aimed at maintaining genome integrity before cells move into mitosis, and with those mechanisms aimed at genome expression control and regulation. In order to clarify the biological basis of these two physiological properties and their interconnectedness, we look at both competencies at the core of defense/adaptive mechanisms and survival, between undifferentiated cell proliferation and organization, constituting a transition phase between two different dynamic regimes, a typical feature of critical dynamic systems. Opting for a candidate-gene strategy, several gene families could be proposed as relevant targets for investigating this hypothesis at the molecular level.

Exploiting genomic knowledge in optimising molecular breeding programmes: algorithms from evolutionary computing

Comparatively few studies have addressed directly the question of quantifying the benefits to be had from using molecular genetic markers in experimental breeding programmes (e.g. for improved crops and livestock), nor the question of which organisms should be mated with each other to best effect. We argue that this requires in silico modelling, an approach for which there is a large literature in the field of evolutionary computation (EC), but which has not really been applied in this way to experimental breeding programmes. EC seeks to optimise measurable outcomes (phenotypic fitnesses) by optimising in silico the mutation, recombination and selection regimes that are used. We review some of the approaches from EC, and compare experimentally, using a biologically relevant in silico landscape, some algorithms that have knowledge of where they are in the (genotypic) search space (G-algorithms) with some (albeit well-tuned ones) that do not (F-algorithms). For the present kinds of landscapes, F- and G-algorithms were broadly comparable in quality and effectiveness, although we recognise that the G-algorithms were not equipped with any ‘prior knowledge’ of epistatic pathway interactions. This use of algorithms based on machine learning has important implications for the optimisation of experimental breeding programmes in the post-genomic era when we shall potentially have access to the full genome sequence of every organism in a breeding population. The non-proprietary code that we have used is made freely available (via Supplementary information).

High-throughput analysis of peptide-binding modules

Modular protein interaction domains (PIDs) that recognize linear peptide motifs are found in hundreds of proteins within the human genome. Some PIDs such as SH2, 14-3-3, Chromo, and Bromo domains serve to recognize posttranslational modification (PTM) of amino acids (such as phosphorylation, acetylation, methylation, etc.) and translate these into discrete cellular responses. Other modules such as SH3 and PSD-95/Discs-large/ZO-1 (PDZ) domains recognize linear peptide epitopes and serve to organize protein complexes based on localization and regions of elevated concentration. In both cases, the ability to nucleate-specific signaling complexes is in large part dependent on the selectivity of a given protein module for its cognate peptide ligand. High-throughput (HTP) analysis of peptide-binding domains by peptide or protein arrays, phage display, mass spectrometry, or other HTP techniques provides new insight into the potential protein-protein interactions prescribed by individual or even whole families of modules. Systems level analyses have also promoted a deeper understanding of the underlying principles that govern selective protein-protein interactions and how selectivity evolves. Lastly, there is a growing appreciation for the limitations and potential pitfalls associated with HTP analysis of protein-peptide interactomes. This review will examine some of the common approaches utilized for large-scale studies of PIDs and suggest a set of standards for the analysis and validation of datasets from large-scale studies of peptide-binding modules. We will also highlight how data from large-scale studies of modular interaction domain families can provide insight into systems level properties such as the linguistics of selective interactions.

Predicting and characterizing selective multiple drug treatments for metabolic diseases and cancer

BACKGROUND

In the field of drug discovery, assessing the potential of multidrug therapies is a difficult task because of the combinatorial complexity (both theoretical and experimental) and because of the requirements on the selectivity of the therapy. To cope with this problem, we have developed a novel method for the systematic in silico investigation of synergistic effects of currently available drugs on genome-scale metabolic networks.

RESULTS

The algorithm finds the optimal combination of drugs which guarantees the inhibition of an objective function, while minimizing the side effect on the other cellular processes. Two different applications are considered: finding drug synergisms for human metabolic diseases (like diabetes, obesity and hypertension) and finding antitumoral drug combinations with minimal side effect on the normal human cell. The results we obtain are consistent with some of the available therapeutic indications and predict new multiple drug treatments. A cluster analysis on all possible interactions among the currently available drugs indicates a limited variety on the metabolic targets for the approved drugs.

CONCLUSION

The in silico prediction of drug synergisms can represent an important tool for the repurposing of drugs in a realistic perspective which considers also the selectivity of the therapy. Moreover, for a more profitable exploitation of drug-drug interactions, we have shown that also experimental drugs which have a different mechanism of action can be reconsider as potential ingredients of new multicompound therapeutic indications. Needless to say the clues provided by a computational study like ours need in any case to be thoroughly evaluated experimentally.

The crossregulation between ERK and PI3K signaling pathways determines the tumoricidal efficacy of MEK inhibitor

MEK inhibitor has been highlighted as a promising anti-tumor drug but its effect has been reported as varying over a wide range depending on patho-physiological conditions. In this study, we employed a systems approach by combining biochemical experimentation with in silico simulations to investigate the resistance mechanism and functional consequences of MEK inhibitor. To this end, we have developed an extended integrative model of ERK and PI3K signaling pathways by considering the crosstalk between Ras and PI3K, and analyzed the resistance mechanism to the MEK inhibitor under various mutational conditions. We found that the phospho-Akt level under the Raf mutation was remarkably augmented by MEK inhibitor, while the phospho-ERK level was almost completely repressed. These results suggest that bypassing of the ERK signal to the PI3K signal causes the resistance to the MEK inhibitor in a complex oncogenic signaling network. We further investigated the underlying mechanism of the drug resistance and revealed that the MEK inhibitor disrupts the negative feedback loops from ERK to SOS and GAB1, but activates the positive feedback loop composed of GAB1, Ras, and PI3K, which induces the bypass of the ERK signal to the PI3K signal. Based on these core feedback circuits, we suggested promising candidates for combination therapy and examined the improved inhibitory effects.

High-order SNP combinations associated with complex diseases: efficient discovery, statistical power and functional interactions

There has been increased interest in discovering combinations of single-nucleotide polymorphisms (SNPs) that are strongly associated with a phenotype even if each SNP has little individual effect. Efficient approaches have been proposed for searching two-locus combinations from genome-wide datasets. However, for high-order combinations, existing methods either adopt a brute-force search which only handles a small number of SNPs (up to few hundreds), or use heuristic search that may miss informative combinations. In addition, existing approaches lack statistical power because of the use of statistics with high degrees-of-freedom and the huge number of hypotheses tested during combinatorial search. Due to these challenges, functional interactions in high-order combinations have not been systematically explored. We leverage discriminative-pattern-mining algorithms from the data-mining community to search for high-order combinations in case-control datasets. The substantially improved efficiency and scalability demonstrated on synthetic and real datasets with several thousands of SNPs allows the study of several important mathematical and statistical properties of SNP combinations with order as high as eleven. We further explore functional interactions in high-order combinations and reveal a general connection between the increase in discriminative power of a combination over its subsets and the functional coherence among the genes comprising the combination, supported by multiple datasets. Finally, we study several significant high-order combinations discovered from a lung-cancer dataset and a kidney-transplant-rejection dataset in detail to provide novel insights on the complex diseases. Interestingly, many of these associations involve combinations of common variations that occur in small fractions of population. Thus, our approach is an alternative methodology for exploring the genetics of rare diseases for which the current focus is on individually rare variations.

Multi-biomarker pattern for tumor identification and prognosis

In last decades, the basic, clinical, and translational research efforts have been directed to the identification of standard biomarkers associated with the degree of malignancy. There is an increasingly public health concern for earlier detection of cancer development at stages in which successful treatments can be achieved. To meet this urgent clinical demand, early stage cancer biomarkers supported by reliable and robust experimental data that can be readily applicable in the clinical practice, are required. In the current standard protocols, when one or two of the canonical proliferating index biomarkers are analyzed, contradictory results are frequently reached leading to incorrect cancer diagnostic and unsuccessful therapies. Therefore, the identification of other cellular characteristics or signatures present in the tumor cells either alone or in combination with the well-established proliferation markers emerge as an alternative strategy in the improvement of cancer diagnosis and treatment. Because it is well known that several pathways and processes are altered in tumor cells, the concept of "single marker" in cancer results incorrect. Therefore, this review aims to analyze and discuss the proposal that the molecular profile of different genes or proteins in different altered tumor pathways must be established to provide a better global clinical pattern for cancer detection and prognosis.

An update on the strategies in multicomponent activity monitoring within the phytopharmaceutical field

Background

To-date modern drug research has focused on the discovery and synthesis of single active substances. However, multicomponent preparations are gaining increasing importance in the phytopharmaceutical field by demonstrating beneficial properties with respect to efficacy and toxicity.

Discussion

In contrast to single drug combinations, a botanical multicomponent therapeutic possesses a complex repertoire of chemicals that belong to a variety of substance classes. This may explain the frequently observed pleiotropic bioactivity spectra of these compounds, which may also suggest that they possess novel therapeutic opportunities. Interestingly, considerable bioactivity properties are exhibited not only by remedies that contain high doses of phytochemicals with prominent pharmaceutical efficacy, but also preparations that lack a sole active principle component. Despite that each individual substance within these multicomponents has a low molar fraction, the therapeutic activity of these substances is established via a potentialization of their effects through combined and simultaneous attacks on multiple molecular targets. Although beneficial properties may emerge from such a broad range of perturbations on cellular machinery, validation and/or prediction of their activity profiles is accompanied with a variety of difficulties in generic risk-benefit assessments. Thus, it is recommended that a comprehensive strategy is implemented to cover the entirety of multicomponent-multitarget effects, so as to address the limitations of conventional approaches.

Summary

An integration of standard toxicological methods with selected pathway-focused bioassays and unbiased data acquisition strategies (such as gene expression analysis) would be advantageous in building an interaction network model to consider all of the effects, whether they were intended or adverse reactions.

Chemical-genomic profiling: systematic analysis of the cellular targets of bioactive molecules

Chemical-genomic (CG) profiling of bioactive compounds is a powerful approach for drug target identification and mode of action studies. Within the last decade, research focused largely on the development and application of CG approaches in the model yeast Saccharomyces cerevisiae. The success of these methods has sparked interest in transitioning CG profiling to other biological systems to extend clinical and evolutionary relevance. Additionally, CG profiling has proven to enhance drug-synergy screens for developing combinatorial therapies. Herein, we briefly review CG profiling, focusing on emerging cross-species technologies and novel drug-synergy applications, as well as outlining needs within the field.

Parsimonious discovery of synergistic drug combinations

Combination therapies that enhance efficacy or permit reduced dosages to be administered have seen great success in a variety of therapeutic applications. More fundamentally, the discovery of epistatic pathway interactions not only informs pharmacologic intervention but can be used to better understand the underlying biological system. There is, however, no systematic and efficient method to identify interacting activities as candidates for combination therapy and, in particular, to identify those with synergistic activities. We devised a pooled, self-deconvoluting screening paradigm for the efficient comprehensive interrogation of all pairs of compounds in 1000-compound libraries. We demonstrate the power of the method to recover established synergistic interactions between compounds. We then applied this approach to a cell-based screen for anti-inflammatory activities using an assay for lipopolysaccharide/interferon-induced acute phase response of a monocytic cell line. The described method, which is >20 times as efficient as a naïve approach, was used to test all pairs of 1027 bioactive compounds for interleukin-6 suppression, yielding 11 pairs of compounds that show synergy. These 11 pairs all represent the same two activities: β-adrenergic receptor agonists and phosphodiesterase-4 inhibitors. These activities both act through cyclic AMP elevation and are known to be anti-inflammatory alone and to synergize in combination. Thus we show proof of concept for a robust, efficient technique for the identification of synergistic combinations. Such a tool can enable qualitatively new scales of pharmacological research and chemical genetics.