Breeze 2.0: an interactive web-tool for visual analysis and comparison of drug response data.

Functional precision medicine (fPM) offers an exciting, simplified approach to finding the right applications for existing molecules and enhancing therapeutic potential. Integrative and robust tools ensuring high accuracy and reliability of the results are critical. In response to this need, we previously developed Breeze, a drug screening data analysis pipeline, designed to facilitate quality control, dose-response curve fitting, and data visualization in a user-friendly manner. Here, we describe the latest version of Breeze (release 2.0), which implements an array of advanced data exploration capabilities, providing users with comprehensive post-analysis and interactive visualization options that are essential for minimizing false positive/negative outcomes and ensuring accurate interpretation of drug sensitivity and resistance data. The Breeze 2.0 web-tool also enables integrative analysis and cross-comparison of user-uploaded data with publicly available drug response datasets. The updated version incorporates new drug quantification metrics, supports analysis of both multi-dose and single-dose drug screening data and introduces a redesigned, intuitive user interface. With these enhancements, Breeze 2.0 is anticipated to substantially broaden its potential applications in diverse domains of fPM.

Evaluation of statistical approaches for association testing in noisy drug screening data.

BACKGROUND: Identifying associations among biological variables is a major challenge in modern quantitative biological research, particularly given the systemic and statistical noise endemic to biological systems. Drug sensitivity data has proven to be a particularly challenging field for identifying associations to inform patient treatment. RESULTS: To address this, we introduce two semi-parametric variations on the commonly used concordance index: the robust concordance index and the kernelized concordance index (rCI, kCI), which incorporate measurements about the noise distribution from the data. We demonstrate that common statistical tests applied to the concordance index and its variations fail to control for false positives, and introduce efficient implementations to compute p-values using adaptive permutation testing. We then evaluate the statistical power of these coefficients under simulation and compare with Pearson and Spearman correlation coefficients. Finally, we evaluate the various statistics in matching drugs across pharmacogenomic datasets. CONCLUSIONS: We observe that the rCI and kCI are better powered than the concordance index in simulation and show some improvement on real data. Surprisingly, we observe that the Pearson correlation was the most robust to measurement noise among the different metrics.

Computational Pipeline for Rational Drug Combination Screening in Patient-Derived Cells.

In many complex diseases, such as cancers, resistance to monotherapies easily occurs, and longer-term treatment responses often require combinatorial therapies as next-line regimens. However, due to a massive number of possible drug combinations to test, there is a need for systematic and rational approaches to finding safe and effective drug combinations for each individual patient. This protocol describes an ecosystem of computational methods to guide high-throughput combinatorial screening that help experimental researchers to identify optimal drug combinations in terms of synergy, efficacy, and/or selectivity for further preclinical and clinical investigation. The methods are demonstrated in the context of combinatorial screening in primary cells of leukemia patients, where the translational aim is to identify drug combinations that show not only high synergy but also maximal cancer-selectivity. The mechanism-agnostic and cost-effective computational methods are widely applicable to various cancer types, which are amenable to drug testing, as the computational methods take as input only the phenotypic measurements of a subset of drug combinations, without requiring target information or genomic profiles of the patient samples.

Crowdsourced mapping of unexplored target space of kinase inhibitors.

Despite decades of intensive search for compounds that modulate the activity of particular protein targets, a large proportion of the human kinome remains as yet undrugged. Effective approaches are therefore required to map the massive space of unexplored compound-kinase interactions for novel and potent activities. Here, we carry out a crowdsourced benchmarking of predictive algorithms for kinase inhibitor potencies across multiple kinase families tested on unpublished bioactivity data. We find the top-performing predictions are based on various models, including kernel learning, gradient boosting and deep learning, and their ensemble leads to a predictive accuracy exceeding that of single-dose kinase activity assays. We design experiments based on the model predictions and identify unexpected activities even for under-studied kinases, thereby accelerating experimental mapping efforts. The open-source prediction algorithms together with the bioactivities between 95 compounds and 295 kinases provide a resource for benchmarking prediction algorithms and for extending the druggable kinome.

Breeze: an integrated quality control and data analysis application for high-throughput drug screening.

SUMMARY: High-throughput screening (HTS) enables systematic testing of thousands of chemical compounds for potential use as investigational and therapeutic agents. HTS experiments are often conducted in multi-well plates that inherently bear technical and experimental sources of error. Thus, HTS data processing requires the use of robust quality control procedures before analysis and interpretation. Here, we have implemented an open-source analysis application, Breeze, an integrated quality control and data analysis application for HTS data. Furthermore, Breeze enables a reliable way to identify individual drug sensitivity and resistance patterns in cell lines or patient-derived samples for functional precision medicine applications. The Breeze application provides a complete solution for data quality assessment, dose-response curve fitting and quantification of the drug responses along with interactive visualization of the results. AVAILABILITY AND IMPLEMENTATION: The Breeze application with video tutorial and technical documentation is accessible at https://breeze.fimm.fi; the R source code is publicly available at https://github.com/potdarswapnil/Breeze under GNU General Public License v3.0. CONTACT: swapnil.potdar@helsinki.fi. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Patient-Customized Drug Combination Prediction and Testing for T-cell Prolymphocytic Leukemia Patients.

The molecular pathways that drive cancer progression and treatment resistance are highly redundant and variable between individual patients with the same cancer type. To tackle this complex rewiring of pathway cross-talk, personalized combination treatments targeting multiple cancer growth and survival pathways are required. Here we implemented a computational-experimental drug combination prediction and testing (DCPT) platform for efficient in silico prioritization and ex vivo testing in patient-derived samples to identify customized synergistic combinations for individual cancer patients. DCPT used drug-target interaction networks to traverse the massive combinatorial search spaces among 218 compounds (a total of 23,653 pairwise combinations) and identified cancer-selective synergies by using differential single-compound sensitivity profiles between patient cells and healthy controls, hence reducing the likelihood of toxic combination effects. A polypharmacology-based machine learning modeling and network visualization made use of baseline genomic and molecular profiles to guide patient-specific combination testing and clinical translation phases. Using T-cell prolymphocytic leukemia (T-PLL) as a first case study, we show how the DCPT platform successfully predicted distinct synergistic combinations for each of the three T-PLL patients, each presenting with different resistance patterns and synergy mechanisms. In total, 10 of 24 (42%) of selective combination predictions were experimentally confirmed to show synergy in patient-derived samples ex vivo The identified selective synergies among approved drugs, including tacrolimus and temsirolimus combined with BCL-2 inhibitor venetoclax, may offer novel drug repurposing opportunities for treating T-PLL.Significance: An integrated use of functional drug screening combined with genomic and molecular profiling enables patient-customized prediction and testing of drug combination synergies for T-PLL patients. Cancer Res; 78(9); 2407-18. ©2018 AACR.

C-SPADE: a web-tool for interactive analysis and visualization of drug screening experiments through compound-specific bioactivity dendrograms.

The advent of polypharmacology paradigm in drug discovery calls for novel chemoinformatic tools for analyzing compounds' multi-targeting activities. Such tools should provide an intuitive representation of the chemical space through capturing and visualizing underlying patterns of compound similarities linked to their polypharmacological effects. Most of the existing compound-centric chemoinformatics tools lack interactive options and user interfaces that are critical for the real-time needs of chemical biologists carrying out compound screening experiments. Toward that end, we introduce C-SPADE, an open-source exploratory web-tool for interactive analysis and visualization of drug profiling assays (biochemical, cell-based or cell-free) using compound-centric similarity clustering. C-SPADE allows the users to visually map the chemical diversity of a screening panel, explore investigational compounds in terms of their similarity to the screening panel, perform polypharmacological analyses and guide drug-target interaction predictions. C-SPADE requires only the raw drug profiling data as input, and it automatically retrieves the structural information and constructs the compound clusters in real-time, thereby reducing the time required for manual analysis in drug development or repurposing applications. The web-tool provides a customizable visual workspace that can either be downloaded as figure or Newick tree file or shared as a hyperlink with other users. C-SPADE is freely available at http://cspade.fimm.fi/.

Comprehensive Drug Testing of Patient-derived Conditionally Reprogrammed Cells from Castration-resistant Prostate Cancer.

BACKGROUND: Technology development to enable the culture of human prostate cancer (PCa) progenitor cells is required for the identification of new, potentially curative therapies for PCa. OBJECTIVE: We established and characterized patient-derived conditionally reprogrammed cells (CRCs) to assess their biological properties and to apply these to test the efficacies of drugs. DESIGN, SETTING, AND PARTICIPANTS: CRCs were established from seven patient samples with disease ranging from primary PCa to advanced castration-resistant PCa (CRPC). The CRCs were characterized by genomic, transcriptomic, protein expression, and drug profiling. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The phenotypic quantification of the CRCs was done based on immunostaining followed by image analysis with Advanced Cell Classifier using Random Forest supervised machine learning. Copy number aberrations (CNAs) were called from whole-exome sequencing and transcriptomics using in-house pipelines. Dose-response measurements were used to generate multiparameter drug sensitivity scores using R-statistical language. RESULTS AND LIMITATIONS: We generated six benign CRC cultures which all had an androgen receptor-negative, basal/transit-amplifying phenotype with few CNAs. In three-dimensional cell culture, these cells could re-express the androgen receptor. The CRCs from a CRPC patient (HUB.5) displayed multiple CNAs, many of which were shared with the parental tumor. We carried out high-throughput drug-response studies with 306 emerging and clinical cancer drugs. Using the benign CRCs as controls, we identified the Bcl-2 family inhibitor navitoclax as the most potent cancer-specific drug for the CRCs from a CRPC patient. Other drug efficacies included taxanes, mepacrine, and retinoids. CONCLUSIONS: Comprehensive cancer pharmacopeia-wide drug testing of CRCs from a CRPC patient highlighted both known and novel drug sensitivities in PCa, including navitoclax, which is currently being tested in clinical trials of CRPC. PATIENT SUMMARY: We describe an approach to generate patient-derived cancer cells from advanced prostate cancer and apply such cells to discover drugs that could be applied in clinical trials for castration-resistant prostate cancer.

Identification of drug candidates and repurposing opportunities through compound-target interaction networks.

INTRODUCTION: System-wide identification of both on- and off-targets of chemical probes provides improved understanding of their therapeutic potential and possible adverse effects, thereby accelerating and de-risking drug discovery process. Given the high costs of experimental profiling of the complete target space of drug-like compounds, computational models offer systematic means for guiding these mapping efforts. These models suggest the most potent interactions for further experimental or pre-clinical evaluation both in cell line models and in patient-derived material. AREAS COVERED: The authors focus here on network-based machine learning models and their use in the prediction of novel compound-target interactions both in target-based and phenotype-based drug discovery applications. While currently being used mainly in complementing the experimentally mapped compound-target networks for drug repurposing applications, such as extending the target space of already approved drugs, these network pharmacology approaches may also suggest completely unexpected and novel investigational probes for drug development. EXPERT OPINION: Although the studies reviewed here have already demonstrated that network-centric modeling approaches have the potential to identify candidate compounds and selective targets in disease networks, many challenges still remain. In particular, these challenges include how to incorporate the cellular context and genetic background into the disease networks to enable more stratified and selective target predictions, as well as how to make the prediction models more realistic for the practical drug discovery and therapeutic applications.

Impact of normalization methods on high-throughput screening data with high hit rates and drug testing with dose-response data.

MOTIVATION: Most data analysis tools for high-throughput screening (HTS) seek to uncover interesting hits for further analysis. They typically assume a low hit rate per plate. Hit rates can be dramatically higher in secondary screening, RNAi screening and in drug sensitivity testing using biologically active drugs. In particular, drug sensitivity testing on primary cells is often based on dose-response experiments, which pose a more stringent requirement for data quality and for intra- and inter-plate variation. Here, we compared common plate normalization and noise-reduction methods, including the B-score and the Loess a local polynomial fit method under high hit-rate scenarios of drug sensitivity testing. We generated simulated 384-well plate HTS datasets, each with 71 plates having a range of 20 (5%) to 160 (42%) hits per plate, with controls placed either at the edge of the plates or in a scattered configuration. RESULTS: We identified 20% (77/384) as the critical hit-rate after which the normalizations started to perform poorly. Results from real drug testing experiments supported this estimation. In particular, the B-score resulted in incorrect normalization of high hit-rate plates, leading to poor data quality, which could be attributed to its dependency on the median polish algorithm. We conclude that a combination of a scattered layout of controls per plate and normalization using a polynomial least squares fit method, such as Loess helps to reduce column, row and edge effects in HTS experiments with high hit-rates and is optimal for generating accurate dose-response curves. CONTACT: john.mpindi@helsinki.fi. AVAILABILITY AND IMPLEMENTATION: Supplementary information: R code and Supplementary data are available at Bioinformatics online.

Novel Lignan and stilbenoid mixture shows anticarcinogenic efficacy in preclinical PC-3M-luc2 prostate cancer model.

Prostate cancer is the most common cancer of men in the Western world, and novel approaches for prostate cancer risk reduction are needed. Plant-derived phenolic compounds attenuate prostate cancer growth in preclinical models by several mechanisms, which is in line with epidemiological findings suggesting that consumption of plant-based diets is associated with low risk of prostate cancer. The objective of this study was to assess the effects of a novel lignan-stilbenoid mixture in PC-3M-luc2 human prostate cancer cells in vitro and in orthotopic xenografts. Lignan and stilbenoid -rich extract was obtained from Scots pine (Pinus sylvestris) knots. Pine knot extract as well as stilbenoids (methyl pinosylvin and pinosylvin), and lignans (matairesinol and nortrachelogenin) present in pine knot extract showed antiproliferative and proapoptotic efficacy at ≥ 40 µM concentration in vitro. Furthermore, pine knot extract derived stilbenoids enhanced tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) induced apoptosis already at ≥ 10 µM concentrations. In orthotopic PC-3M-luc2 xenograft bearing immunocompromized mice, three-week peroral exposure to pine knot extract (52 mg of lignans and stilbenoids per kg of body weight) was well tolerated and showed anti-tumorigenic efficacy, demonstrated by multivariate analysis combining essential markers of tumor growth (i.e. tumor volume, vascularization, and cell proliferation). Methyl pinosylvin, pinosylvin, matairesinol, nortrachelogenin, as well as resveratrol, a metabolite of pinosylvin, were detected in serum at total concentration of 7-73 µM, confirming the bioavailability of pine knot extract derived lignans and stilbenoids. In summary, our data indicates that pine knot extract is a novel and cost-effective source of resveratrol, methyl pinosylvin and other bioactive lignans and stilbenoids. Pine knot extract shows anticarcinogenic efficacy in preclinical prostate cancer model, and our in vitro data suggests that compounds derived from the extract may have potential as novel chemosensitizers to TRAIL. These findings promote further research on health-related applications of wood biochemicals.

A comparative evaluation of software for the analysis of liquid chromatography-tandem mass spectrometry data from isotope coded affinity tag experiments.

The options available for processing quantitative data from isotope coded affinity tag (ICAT) experiments have mostly been confined to software specific to the instrument of acquisition. However, recent developments with data format conversion have subsequently increased such processing opportunities. In the present study, data sets from ICAT experiments, analysed with liquid chromatography/tandem mass spectrometry (MS/MS), using an Applied Biosystems QSTAR Pulsar quadrupole-TOF mass spectrometer, were processed in triplicate using separate mass spectrometry software packages. The programs Pro ICAT, Spectrum Mill and SEQUEST with XPRESS were employed. Attention was paid towards the extent of common identification and agreement of quantitative results, with additional interest in the flexibility and productivity of these programs. The comparisons were made with data from the analysis of a specifically prepared test mixture, nine proteins at a range of relative concentration ratios from 0.1 to 10 (light to heavy labelled forms), as a known control, and data selected from an ICAT study involving the measurement of cytokine induced protein expression in human lymphoblasts, as an applied example. Dissimilarities were detected in peptide identification that reflected how the associated scoring parameters favoured information from the MS/MS data sets. Accordingly, there were differences in the numbers of peptides and protein identifications, although from these it was apparent that both confirmatory and complementary information was present. In the quantitative results from the three programs, no statistically significant differences were observed.