A synthetic lethal dependency on casein kinase 2 in response to replication-perturbing therapeutics in RB1-deficient cancer cells.

Resistance to therapy commonly develops in patients with high-grade serous ovarian carcinoma (HGSC) and triple-negative breast cancer (TNBC), urging the search for improved therapeutic combinations and their predictive biomarkers. Starting from a CRISPR knockout screen, we identified that loss of RB1 in TNBC or HGSC cells generates a synthetic lethal dependency on casein kinase 2 (CK2) for surviving the treatment with replication-perturbing therapeutics such as carboplatin, gemcitabine, or PARP inhibitors. CK2 inhibition in RB1-deficient cells resulted in the degradation of another RB family cell cycle regulator, p130, which led to S phase accumulation, micronuclei formation, and accelerated PARP inhibition-induced aneuploidy and mitotic cell death. CK2 inhibition was also effective in primary patient-derived cells. It selectively prevented the regrowth of RB1-deficient patient HGSC organoids after treatment with carboplatin or niraparib. As about 25% of HGSCs and 40% of TNBCs have lost RB1 expression, CK2 inhibition is a promising approach to overcome resistance to standard therapeutics in large strata of patients.

Development of HDAC Inhibitors Exhibiting Therapeutic Potential in T-Cell Prolymphocytic Leukemia.

Epigenetic targeting has emerged as an efficacious therapy for hematological cancers. The rare and incurable T-cell prolymphocytic leukemia (T-PLL) is known for its aggressive clinical course. Current epigenetic agents such as histone deacetylase (HDAC) inhibitors are increasingly used for targeted therapy. Through a structure-activity relationship (SAR) study, we developed an HDAC6 inhibitor KT-531, which exhibited higher potency in T-PLL compared to other hematological cancers. KT-531 displayed strong HDAC6 inhibitory potency and selectivity, on-target biological activity, and a safe therapeutic window in nontransformed cell lines. In primary T-PLL patient cells, where HDAC6 was found to be overexpressed, KT-531 exhibited strong biological responses, and safety in healthy donor samples. Notably, combination studies in T-PLL patient samples demonstrated KT-531 synergizes with approved cancer drugs, bendamustine, idasanutlin, and venetoclax. Our work suggests HDAC inhibition in T-PLL could afford sufficient therapeutic windows to achieve durable remission either as stand-alone or in combination with targeted drugs.

Regression plane concept for analysing continuous cellular processes with machine learning.

Biological processes are inherently continuous, and the chance of phenotypic discovery is significantly restricted by discretising them. Using multi-parametric active regression we introduce the Regression Plane (RP), a user-friendly discovery tool enabling class-free phenotypic supervised machine learning, to describe and explore biological data in a continuous manner. First, we compare traditional classification with regression in a simulated experimental setup. Second, we use our framework to identify genes involved in regulating triglyceride levels in human cells. Subsequently, we analyse a time-lapse dataset on mitosis to demonstrate that the proposed methodology is capable of modelling complex processes at infinite resolution. Finally, we show that hemocyte differentiation in Drosophila melanogaster has continuous characteristics.

Multi-modal meta-analysis of cancer cell line omics profiles identifies ECHDC1 as a novel breast tumor suppressor.

Molecular and functional profiling of cancer cell lines is subject to laboratory-specific experimental practices and data analysis protocols. The current challenge therefore is how to make an integrated use of the omics profiles of cancer cell lines for reliable biological discoveries. Here, we carried out a systematic analysis of nine types of data modalities using meta-analysis of 53 omics studies across 12 research laboratories for 2,018 cell lines. To account for a relatively low consistency observed for certain data modalities, we developed a robust data integration approach that identifies reproducible signals shared among multiple data modalities and studies. We demonstrated the power of the integrative analyses by identifying a novel driver gene, ECHDC1, with tumor suppressive role validated both in breast cancer cells and patient tumors. The multi-modal meta-analysis approach also identified synthetic lethal partners of cancer drivers, including a co-dependency of PTEN deficient endometrial cancer cells on RNA helicases.

Improved detection of differentially represented DNA barcodes for high-throughput clonal phenomics.

Cellular DNA barcoding has become a popular approach to study heterogeneity of cell populations and to identify clones with differential response to cellular stimuli. However, there is a lack of reliable methods for statistical inference of differentially responding clones. Here, we used mixtures of DNA-barcoded cell pools to generate a realistic benchmark read count dataset for modelling a range of outcomes of clone-tracing experiments. By accounting for the statistical properties intrinsic to the DNA barcode read count data, we implemented an improved algorithm that results in a significantly lower false-positive rate, compared to current RNA-seq data analysis algorithms, especially when detecting differentially responding clones in experiments with strong selection pressure. Building on the reliable statistical methodology, we illustrate how multidimensional phenotypic profiling enables one to deconvolute phenotypically distinct clonal subpopulations within a cancer cell line. The mixture control dataset and our analysis results provide a foundation for benchmarking and improving algorithms for clone-tracing experiments.

SynToxProfiler: An interactive analysis of drug combination synergy, toxicity and efficacy.

Drug combinations are becoming a standard treatment of many complex diseases due to their capability to overcome resistance to monotherapy. In the current preclinical drug combination screening, the top combinations for further study are often selected based on synergy alone, without considering the combination efficacy and toxicity effects, even though these are critical determinants for the clinical success of a therapy. To promote the prioritization of drug combinations based on integrated analysis of synergy, efficacy and toxicity profiles, we implemented a web-based open-source tool, SynToxProfiler (Synergy-Toxicity-Profiler). When applied to 20 anti-cancer drug combinations tested both in healthy control and T-cell prolymphocytic leukemia (T-PLL) patient cells, as well as to 77 anti-viral drug pairs tested in Huh7 liver cell line with and without Ebola virus infection, SynToxProfiler prioritized as top hits those synergistic drug pairs that showed higher selective efficacy (difference between efficacy and toxicity), which offers an improved likelihood for clinical success.

JAK/STAT-Activating Genomic Alterations Are a Hallmark of T-PLL.

T-cell prolymphocytic leukemia (T-PLL) is a rare and poor-prognostic mature T-cell leukemia. Recent studies detected genomic aberrations affecting JAK and STAT genes in T-PLL. Due to the limited number of primary patient samples available, genomic analyses of the JAK/STAT pathway have been performed in rather small cohorts. Therefore, we conducted-via a primary-data based pipeline-a meta-analysis that re-evaluated the genomic landscape of T-PLL. It included all available data sets with sequence information on JAK or STAT gene loci in 275 T-PLL. We eliminated overlapping cases and determined a cumulative rate of 62.1% of cases with mutated JAK or STAT genes. Most frequently, JAK1 (6.3%), JAK3 (36.4%), and STAT5B (18.8%) carried somatic single-nucleotide variants (SNVs), with missense mutations in the SH2 or pseudokinase domains as most prevalent. Importantly, these lesions were predominantly subclonal. We did not detect any strong association between mutations of a JAK or STAT gene with clinical characteristics. Irrespective of the presence of gain-of-function (GOF) SNVs, basal phosphorylation of STAT5B was elevated in all analyzed T-PLL. Fittingly, a significant proportion of genes encoding for potential negative regulators of STAT5B showed genomic losses (in 71.4% of T-PLL in total, in 68.4% of T-PLL without any JAK or STAT mutations). They included DUSP4, CD45, TCPTP, SHP1, SOCS1, SOCS3, and HDAC9. Overall, considering such losses of negative regulators and the GOF mutations in JAK and STAT genes, a total of 89.8% of T-PLL revealed a genomic aberration potentially explaining enhanced STAT5B activity. In essence, we present a comprehensive meta-analysis on the highly prevalent genomic lesions that affect genes encoding JAK/STAT signaling components. This provides an overview of possible modes of activation of this pathway in a large cohort of T-PLL. In light of new advances in JAK/STAT inhibitor development, we also outline translational contexts for harnessing active JAK/STAT signaling, which has emerged as a 'secondary' hallmark of T-PLL.

Chemogenomic Analysis of the Druggable Kinome and Its Application to Repositioning and Lead Identification Studies.

Owing to the intrinsic polypharmacological nature of most small-molecule kinase inhibitors, there is a need for computational models that enable systematic exploration of the chemogenomic landscape underlying druggable kinome toward more efficient kinome-profiling strategies. We implemented VirtualKinomeProfiler, an efficient computational platform that captures distinct representations of chemical similarity space of the druggable kinome for various drug discovery endeavors. By using the computational platform, we profiled approximately 37 million compound-kinase pairs and made predictions for 151,708 compounds in terms of their repositioning and lead molecule potential, against 248 kinases simultaneously. Experimental testing with biochemical assays validated 51 of the predicted interactions, identifying 19 small-molecule inhibitors of EGFR, HCK, FLT1, and MSK1 protein kinases. The prediction model led to a 1.5-fold increase in precision and 2.8-fold decrease in false-discovery rate, when compared with traditional single-dose biochemical screening, which demonstrates its potential to drastically expedite the kinome-specific drug discovery process.

Network pharmacology modeling identifies synergistic Aurora B and ZAK interaction in triple-negative breast cancer.

Cancer cells with heterogeneous mutation landscapes and extensive functional redundancy easily develop resistance to monotherapies by emerging activation of compensating or bypassing pathways. To achieve more effective and sustained clinical responses, synergistic interactions of multiple druggable targets that inhibit redundant cancer survival pathways are often required. Here, we report a systematic polypharmacology strategy to predict, test, and understand the selective drug combinations for MDA-MB-231 triple-negative breast cancer cells. We started by applying our network pharmacology model to predict synergistic drug combinations. Next, by utilizing kinome-wide drug-target profiles and gene expression data, we pinpointed a synergistic target interaction between Aurora B and ZAK kinase inhibition that led to enhanced growth inhibition and cytotoxicity, as validated by combinatorial siRNA, CRISPR/Cas9, and drug combination experiments. The mechanism of such a context-specific target interaction was elucidated using a dynamic simulation of MDA-MB-231 signaling network, suggesting a cross-talk between p53 and p38 pathways. Our results demonstrate the potential of polypharmacological modeling to systematically interrogate target interactions that may lead to clinically actionable and personalized treatment options.

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.

Computational-experimental approach to drug-target interaction mapping: A case study on kinase inhibitors.

Due to relatively high costs and labor required for experimental profiling of the full target space of chemical compounds, various machine learning models have been proposed as cost-effective means to advance this process in terms of predicting the most potent compound-target interactions for subsequent verification. However, most of the model predictions lack direct experimental validation in the laboratory, making their practical benefits for drug discovery or repurposing applications largely unknown. Here, we therefore introduce and carefully test a systematic computational-experimental framework for the prediction and pre-clinical verification of drug-target interactions using a well-established kernel-based regression algorithm as the prediction model. To evaluate its performance, we first predicted unmeasured binding affinities in a large-scale kinase inhibitor profiling study, and then experimentally tested 100 compound-kinase pairs. The relatively high correlation of 0.77 (p < 0.0001) between the predicted and measured bioactivities supports the potential of the model for filling the experimental gaps in existing compound-target interaction maps. Further, we subjected the model to a more challenging task of predicting target interactions for such a new candidate drug compound that lacks prior binding profile information. As a specific case study, we used tivozanib, an investigational VEGF receptor inhibitor with currently unknown off-target profile. Among 7 kinases with high predicted affinity, we experimentally validated 4 new off-targets of tivozanib, namely the Src-family kinases FRK and FYN A, the non-receptor tyrosine kinase ABL1, and the serine/threonine kinase SLK. Our sub-sequent experimental validation protocol effectively avoids any possible information leakage between the training and validation data, and therefore enables rigorous model validation for practical applications. These results demonstrate that the kernel-based modeling approach offers practical benefits for probing novel insights into the mode of action of investigational compounds, and for the identification of new target selectivities for drug repurposing applications.

The impact of low-frequency and rare variants on lipid levels.

Using a genome-wide screen of 9.6 million genetic variants achieved through 1000 Genomes Project imputation in 62,166 samples, we identify association to lipid traits in 93 loci, including 79 previously identified loci with new lead SNPs and 10 new loci, 15 loci with a low-frequency lead SNP and 10 loci with a missense lead SNP, and 2 loci with an accumulation of rare variants. In six loci, SNPs with established function in lipid genetics (CELSR2, GCKR, LIPC and APOE) or candidate missense mutations with predicted damaging function (CD300LG and TM6SF2) explained the locus associations. The low-frequency variants increased the proportion of variance explained, particularly for low-density lipoprotein cholesterol and total cholesterol. Altogether, our results highlight the impact of low-frequency variants in complex traits and show that imputation offers a cost-effective alternative to resequencing.