Drug combination discovery depends on reliable synergy metrics but no consensus exists on the correct synergy criterion to characterize combined interactions. The fragmented state of the field confounds analysis, impedes reproducibility, and delays clinical translation of potential combination treatments. Here we present a mass-action based formalism to quantify synergy. With this formalism, we clarify the relationship between the dominant drug synergy principles, and present a mapping of commonly used frameworks onto a unified synergy landscape. From this, we show how biases emerge due to intrinsic assumptions which hinder their broad applicability and impact the interpretation of synergy in discovery efforts. Specifically, we describe how traditional metrics mask consequential synergistic interactions, and contain biases dependent on the Hill-slope and maximal effect of single-drugs. We show how these biases systematically impact synergy classification in large combination screens, potentially misleading discovery efforts. Thus the proposed formalism can provide a consistent, unbiased interpretation of drug synergy, and accelerate the translatability of synergy studies.

An understanding of how cells respond to perturbation is essential for biological applications; however, most approaches for profiling cellular response are limited in scope to pre-established targets. Global analysis of molecular mechanism will advance our understanding of the complex networks constituting cellular perturbation and lead to advancements in areas, such as infectious disease pathogenesis, developmental biology, pathophysiology, pharmacology, and toxicology. We have developed a high-throughput multiomics platform for comprehensive, de novo characterization of cellular mechanisms of action. Platform validation using cisplatin as a test compound demonstrates quantification of over 10 000 unique, significant molecular changes in less than 30 days. These data provide excellent coverage of known cisplatin-induced molecular changes and previously unrecognized insights into cisplatin resistance. This proof-of-principle study demonstrates the value of this platform as a resource to understand complex cellular responses in a high-throughput manner.

BackgroundThere is a lack of biomarkers to predict outcome with targeted therapy in metastatic clear cell renal cancer (mccRCC). This may be because dynamic molecular changes occur with therapy.ObjectiveTo explore if dynamic, targeted-therapy-driven molecular changes correlate with mccRCC outcome.Design, setting, and participantsMultiple frozen samples from primary tumours were taken from sunitinib-naïve (n = 22) and sunitinib-treated mccRCC patients (n = 23) for protein analysis. A cohort (n = 86) of paired, untreated and sunitinib/pazopanib-treated mccRCC samples was used for validation. Array comparative genomic hybridisation (CGH) analysis and RNA interference (RNAi) was used to support the findings.InterventionThree cycles of sunitinib 50 mg (4 wk on, 2 wk off).Outcome measurements and statistical analysisReverse phase protein arrays (training set) and immunofluorescence automated quantitative analysis (validation set) assessed protein expression.Results and limitationsDifferential expression between sunitinib-naïve and treated samples was seen in 30 of 55 proteins (p < 0.05 for each). The proteins B-cell CLL/lymphoma 2 (BCL2), mutL homolog 1 (MLH1), carbonic anhydrase 9 (CA9), and mechanistic target of rapamycin (mTOR) (serine/threonine kinase) had both increased intratumoural variance and significant differential expression with therapy. The validation cohort confirmed increased CA9 expression with therapy. Multivariate analysis showed high CA9 expression after treatment was associated with longer survival (hazard ratio: 0.48; 95% confidence interval, 0.26–0.87; p = 0.02). Array CGH profiles revealed sunitinib was associated with significant CA9 region loss. RNAi CA9 silencing in two cell lines inhibited the antiproliferative effects of sunitinib. Shortcomings of the study include selection of a specific protein for analysis, and the specific time points at which the treated tissue was analysed.ConclusionsCA9 levels increase with targeted therapy in mccRCC. Lower CA9 levels are associated with a poor prognosis and possible resistance, as indicated by the validation cohort.Patient summaryDrug treatment of advanced kidney cancer alters molecular markers of treatment resistance. Measuring carbonic anhydrase 9 levels may be helpful in determining which patients benefit from therapy.

Purpose: The aim of this study was to investigate the effect of VEGF-targeted therapy (sunitinib) on molecular intratumoral heterogeneity (ITH) in metastatic clear cell renal cancer (mccRCC).Experimental Design: Multiple tumor samples (n ¼ 187 samples) were taken from the primary renal tumors of patients with mccRCC who were sunitinib treated (n ¼ 23, SuMR clinical trial) or untreated (n ¼ 23, SCOTRRCC study). ITH of pathologic grade, DNA (aCGH), mRNA (Illumina Beadarray) and candidate proteins (reverse phase protein array) were evaluated using unsupervised and supervised analyses (driver mutations, hypoxia, and stromal-related genes). ITH was analyzed using intratumoral protein variance distributions and distribution of individual patient aCGH and gene-expression clustering.Results: Tumor grade heterogeneity was greater in treated compared with untreated tumors (P ¼ 0.002). In unsupervised analysis, sunitinib therapy was not associated with increased ITH in DNA or mRNA. However, there was an increase in ITH for the driver mutation gene signature (DNA and mRNA) as well as increasing variability of protein expression with treatment (P < 0.05). Despite this variability, significant chromosomal and transcript changes to key targets of sunitinib, such as VHL, PBRM1, and CAIX, occurred in the treated samples.Conclusions: These findings suggest that sunitinib treatment has significant effects on the expression and ITH of key tumor and treatment specific genes/proteins in mccRCC. The results, based on primary tumor analysis, do not support the hypothesis that resistant clones are selected and predominate following targeted therapy.

Tissue microarrays (TMAs) allow multiplexed analysis of tissue samples and are frequently used to estimate biomarker protein expression in tumour biopsies. TMA Navigator (www.tmanavigator.org) is an open access web application for analysis of TMA data and related information, accommodating categorical, semi-continuous and continuous expression scores. Non-biological variation, or batch effects, can hinder data analysis and may be mitigated using the ComBat algorithm, which is incorporated with enhancements for automated application to TMA data. Unsupervised grouping of samples (patients) is provided according to Gaussian mixture modelling of marker scores, with cardinality selected by Bayesian information criterion regularization. Kaplan–Meier survival analysis is available, including comparison of groups identified by mixture modelling using the Mantel-Cox log-rank test. TMA Navigator also supports network inference approaches useful for TMA datasets, which often constitute comparatively few markers. Tissue and cell-type specific networks derived from TMA expression data offer insights into the molecular logic underlying pathophenotypes, towards more effective and personalized medicine. Output is interactive, and results may be exported for use with external programs. Private anonymous access is available, and user accounts may be generated for easier data management.

SummaryA major barrier to the practical utilization of large, complex models of biochemical systems is the lack of open-source computational tools to evaluate model behaviors over high-dimensional parameter spaces. This is due to the high computational expense of performing thousands to millions of model simulations required for statistical analysis. To address this need, we have implemented a user-friendly interface between cupSODA, a GPU-powered kinetic simulator, and PySB, a Python-based modeling and simulation framework. For three example models of varying size, we show that for large numbers of simulations PySB/cupSODA achieves order-of-magnitude speedups relative to a CPU-based ordinary differential equation integrator.Availability and implementationThe PySB/cupSODA interface has been integrated into the PySB modeling framework (version 1.4.0), which can be installed from the Python Package Index (PyPI) using a Python package manager such as pip. cupSODA source code and precompiled binaries (Linux, Mac OS/X, Windows) are available at github.com/aresio/cupSODA (requires an Nvidia GPU; developer.nvidia.com/cuda-gpus). Additional information about PySB is available at pysb.org.Supplementary information Supplementary data are available at Bioinformatics online.

Advances in microscopy imaging technologies have enabled the visualization of live-cell dynamic processes using time-lapse microscopy imaging. However, modern methods exhibit several limitations related to the training phases and to time constraints, hindering their application in the laboratory practice. In this work, we present a novel method, named Automated Cell Detection and Counting (ACDC), designed for activity detection of fluorescent labeled cell nuclei in time-lapse microscopy. ACDC overcomes the limitations of the literature methods, by first applying bilateral filtering on the original image to smooth the input cell images while preserving edge sharpness, and then by exploiting the watershed transform and morphological filtering. Moreover, ACDC represents a feasible solution for the laboratory practice, as it can leverage multi-core architectures in computer clusters to efficiently handle large-scale imaging datasets. Indeed, our Parent-Workers implementation of ACDC allows to obtain up to a 3.7× speed-up compared to the sequential counterpart. ACDC was tested on two distinct cell imaging datasets to assess its accuracy and effectiveness on images with different characteristics. We achieved an accurate cell-count and nuclei segmentation without relying on large-scale annotated datasets, a result confirmed by the average Dice Similarity Coefficients of 76.84 and 88.64 and the Pearson coefficients of 0.99 and 0.96, calculated against the manual cell counting, on the two tested datasets.

Cell identity is governed by gene expression, regulated by transcription factor (TF) binding at cis-regulatory modules. Decoding the relationship between TF binding patterns and gene regulation is nontrivial, remaining a fundamental limitation in understanding cell decision-making. We developed the NetNC software to predict functionally active regulation of TF targets; demonstrated on nine datasets for the TFs Snail, Twist, and modENCODE Highly Occupied Target (HOT) regions. Snail and Twist are canonical drivers of epithelial to mesenchymal transition (EMT), a cell programme important in development, tumour progression and fibrosis. Predicted “neutral” (non-functional) TF binding always accounted for the majority (50% to 95%) of candidate target genes from statistically significant peaks and HOT regions had higher functional binding than most of the Snail and Twist datasets examined. Our results illuminated conserved gene networks that control epithelial plasticity in development and disease. We identified new gene functions and network modules including crosstalk with notch signalling and regulation of chromatin organisation, evidencing networks that reshape Waddington’s epigenetic landscape during epithelial remodelling. Expression of orthologous functional TF targets discriminated breast cancer molecular subtypes and predicted novel tumour biology, with implications for precision medicine. Predicted invasion roles were validated using a tractable cell model, supporting our approach.