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.

Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct.

The extent to which genetic risk factors are shared between childhood-onset (COA) and adult-onset (AOA) asthma has not been estimated. On the basis of data from the UK Biobank study (n = 447,628), we found that the variance in disease liability explained by common variants is higher for COA (onset at ages between 0 and 19 years; h2g = 25.6%) than for AOA (onset at ages between 20 and 60 years; h2g = 10.6%). The genetic correlation (rg) between COA and AOA was 0.67. Variation in age of onset among COA-affected individuals had a low heritability (h2g = 5%), which we confirmed in independent studies and also among AOA-affected individuals. To identify subtype-specific genetic associations, we performed a genome-wide association study (GWAS) in the UK Biobank for COA (13,962 affected individuals) and a separate GWAS for AOA (26,582 affected individuals) by using a common set of 300,671 controls for both studies. We identified 123 independent associations for COA and 56 for AOA (37 overlapped); of these, 98 and 34, respectively, were reproducible in an independent study (n = 262,767). Collectively, 28 associations were not previously reported. For 96 COA-associated variants, including five variants that represent COA-specific risk factors, the risk allele was more common in COA- than in AOA-affected individuals. Conversely, we identified three variants that are stronger risk factors for AOA. Variants associated with obesity and smoking had a stronger contribution to the risk of AOA than to the risk of COA. Lastly, we identified 109 likely target genes of the associated variants, primarily on the basis of correlated expression quantitative trait loci (up to n = 31,684). GWAS informed by age of onset can identify subtype-specific risk variants, which can help us understand differences in pathophysiology between COA and AOA and so can be informative for drug development.

Identification of selective cytotoxic and synthetic lethal drug responses in triple negative breast cancer cells.

BACKGROUND: Triple negative breast cancer (TNBC) is a highly heterogeneous and aggressive type of cancer that lacks effective targeted therapy. Despite detailed molecular profiling, no targeted therapy has been established. Hence, with the aim of gaining deeper understanding of the functional differences of TNBC subtypes and how that may relate to potential novel therapeutic strategies, we studied comprehensive anticancer-agent responses among a panel of TNBC cell lines. METHOD: The responses of 301 approved and investigational oncology compounds were measured in 16 TNBC cell lines applying a functional profiling approach. To go beyond the standard drug viability effect profiling, which has been used in most chemosensitivity studies, we utilized a multiplexed readout for both cell viability and cytotoxicity, allowing us to differentiate between cytostatic and cytotoxic responses. RESULTS: Our approach revealed that most single-agent anti-cancer compounds that showed activity for the viability readout had no or little cytotoxic effects. Major compound classes that exhibited this type of response included anti-mitotics, mTOR, CDK, and metabolic inhibitors, as well as many agents selectively inhibiting oncogene-activated pathways. However, within the broad viability-acting classes of compounds, there were often subsets of cell lines that responded by cell death, suggesting that these cells are particularly vulnerable to the tested substance. In those cases we could identify differential levels of protein markers associated with cytotoxic responses. For example, PAI-1, MAPK phosphatase and Notch-3 levels associated with cytotoxic responses to mitotic and proteasome inhibitors, suggesting that these might serve as markers of response also in clinical settings. Furthermore, the cytotoxicity readout highlighted selective synergistic and synthetic lethal drug combinations that were missed by the cell viability readouts. For instance, the MEK inhibitor trametinib synergized with PARP inhibitors. Similarly, combination of two non-cytotoxic compounds, the rapamycin analog everolimus and an ATP-competitive mTOR inhibitor dactolisib, showed synthetic lethality in several mTOR-addicted cell lines. CONCLUSIONS: Taken together, by studying the combination of cytotoxic and cytostatic drug responses, we identified a deeper spectrum of cellular responses both to single agents and combinations that may be highly relevant for identifying precision medicine approaches in TNBC as well as in other types of cancers.

Relevance Rank Platform (RRP) for Functional Filtering of High Content Protein-Protein Interaction Data.

High content protein interaction screens have revolutionized our understanding of protein complex assembly. However, one of the major challenges in translation of high content protein interaction data is identification of those interactions that are functionally relevant for a particular biological question. To address this challenge, we developed a relevance ranking platform (RRP), which consist of modular functional and bioinformatic filters to provide relevance rank among the interactome proteins. We demonstrate the versatility of RRP to enable a systematic prioritization of the most relevant interaction partners from high content data, highlighted by the analysis of cancer relevant protein interactions for oncoproteins Pin1 and PME-1. We validated the importance of selected interactions by demonstration of PTOV1 and CSKN2B as novel regulators of Pin1 target c-Jun phosphorylation and reveal previously unknown interacting proteins that may mediate PME-1 effects via PP2A-inhibition. The RRP framework is modular and can be modified to answer versatile research problems depending on the nature of the biological question under study. Based on comparison of RRP to other existing filtering tools, the presented data indicate that RRP offers added value especially for the analysis of interacting proteins for which there is no sufficient prior knowledge available. Finally, we encourage the use of RRP in combination with either SAINT or CRAPome computational tools for selecting the candidate interactors that fulfill the both important requirements, functional relevance, and high confidence interaction detection.

From drug response profiling to target addiction scoring in cancer cell models.

Deconvoluting the molecular target signals behind observed drug response phenotypes is an important part of phenotype-based drug discovery and repurposing efforts. We demonstrate here how our network-based deconvolution approach, named target addiction score (TAS), provides insights into the functional importance of druggable protein targets in cell-based drug sensitivity testing experiments. Using cancer cell line profiling data sets, we constructed a functional classification across 107 cancer cell models, based on their common and unique target addiction signatures. The pan-cancer addiction correlations could not be explained by the tissue of origin, and only correlated in part with molecular and genomic signatures of the heterogeneous cancer cells. The TAS-based cancer cell classification was also shown to be robust to drug response data resampling, as well as predictive of the transcriptomic patterns in an independent set of cancer cells that shared similar addiction signatures with the 107 cancers. The critical protein targets identified by the integrated approach were also shown to have clinically relevant mutation frequencies in patients with various cancer subtypes, including not only well-established pan-cancer genes, such as PTEN tumor suppressor, but also a number of targets that are less frequently mutated in specific cancer types, including ABL1 oncoprotein in acute myeloid leukemia. An application to leukemia patient primary cell models demonstrated how the target deconvolution approach offers functional insights into patient-specific addiction patterns, such as those indicative of their receptor-type tyrosine-protein kinase FLT3 internal tandem duplication (FLT3-ITD) status and co-addiction partners, which may lead to clinically actionable, personalized drug treatment developments. To promote its application to the future drug testing studies, we have made available an open-source implementation of the TAS calculation in the form of a stand-alone R package.

Quantitative analysis of PPT1 interactome in human neuroblastoma cells.

Mutations in the CLN1 gene that encodes Palmitoyl protein thioesterase 1 (PPT1) or CLN1, cause Infantile NCL (INCL, MIM#256730). PPT1 removes long fatty acid chains such as palmitate from modified cysteine residues of proteins. The data shown here result from isolated protein complexes from PPT1-expressing SH-SY5Y stable cells that were subjected to single step affinity purification coupled to mass spectrometry (AP-MS). Prior to the MS analysis, we utilised a modified filter-aided sample preparation (FASP) protocol. Based on label free quantitative analysis of the data by SAINT, 23 PPT1 interacting partners (IP) were identified. A dense connectivity in PPT1 network was further revealed by functional coupling and extended network analyses, linking it to mitochondrial ATP synthesis coupled protein transport and thioester biosynthetic process. Moreover, the terms: inhibition of organismal death, movement disorders and concentration of lipid were predicted to be altered in the PPT1 network. Data presented here are related to Scifo et al. (J. Proteomics, 123 (2015) 42-53).

Systematic Mapping of Kinase Addiction Combinations in Breast Cancer Cells by Integrating Drug Sensitivity and Selectivity Profiles.

Chemical perturbation screens offer the possibility to identify actionable sets of cancer-specific vulnerabilities. However, most inhibitors of kinases or other cancer targets result in polypharmacological effects, which complicate the identification of target dependencies directly from the drug-response phenotypes. In this study, we developed a chemical systems biology approach that integrates comprehensive drug sensitivity and selectivity profiling to provide functional insights into both single and multi-target oncogenic signal addictions. When applied to 21 breast cancer cell lines, perturbed with 40 kinase inhibitors, the subtype-specific addiction patterns clustered in agreement with patient-derived subtypes, while showing considerable variability between the heterogeneous breast cancers. Experimental validation of the top predictions revealed a number of co-dependencies between kinase targets that led to unexpected synergistic combinations between their inhibitors, such as dasatinib and axitinib in the triple-negative basal-like HCC1937 cell line.

Proteomic analysis of the palmitoyl protein thioesterase 1 interactome in SH-SY5Y human neuroblastoma cells.

Neuronal ceroid lipofuscinoses (NCL) are a group of inherited progressive childhood disorders, characterized by early accumulation of autofluorescent storage material in lysosomes of neurons or other cells. Clinical symptoms of NCL include: progressive loss of vision, mental and motor deterioration, epileptic seizures and premature death. CLN1 disease (MIM#256730) is caused by mutations in the CLN1 gene, which encodes palmitoyl protein thioesterase 1 (PPT1). In this study, we utilised single step affinity purification coupled to mass spectrometry (AP-MS) to unravel the in vivo substrates of human PPT1 in the brain neuronal cells. Protein complexes were isolated from human PPT1 expressing SH-SY5Y stable cells, subjected to filter-aided sample preparation (FASP) and analysed on a Q Exactive Hybrid Quadrupole-Orbitrap mass spectrometer. A total of 23 PPT1 interacting partners (IP) were identified from label free quantitation of the MS data by SAINT platform. Three of the identified PPT1 IP, namely CRMP1, DBH, and MAP1B are predicted to be palmitoylated. Our proteomic analysis confirmed previously suggested roles of PPT1 in axon guidance and lipid metabolism, yet implicates the enzyme in novel roles including: involvement in neuronal migration and dopamine receptor mediated signalling pathway. BIOLOGICAL SIGNIFICANCE: The significance of this work lies in the unravelling of putative in vivo substrates of human CLN1 or PPT1 in brain neuronal cells. Moreover, the PPT1 IP implicate the enzyme in novel roles including: involvement in neuronal migration and dopamine receptor mediated signalling pathway.

Quantitative scoring of differential drug sensitivity for individually optimized anticancer therapies.

We developed a systematic algorithmic solution for quantitative drug sensitivity scoring (DSS), based on continuous modeling and integration of multiple dose-response relationships in high-throughput compound testing studies. Mathematical model estimation and continuous interpolation makes the scoring approach robust against sources of technical variability and widely applicable to various experimental settings, both in cancer cell line models and primary patient-derived cells. Here, we demonstrate its improved performance over other response parameters especially in a leukemia patient case study, where differential DSS between patient and control cells enabled identification of both cancer-selective drugs and drug-sensitive patient sub-groups, as well as dynamic monitoring of the response patterns and oncogenic driver signals during cancer progression and relapse in individual patient cells ex vivo. An open-source and easily extendable implementation of the DSS calculation is made freely available to support its tailored application to translating drug sensitivity testing results into clinically actionable treatment options.

Naturally occurring human phosphorylcholine antibodies are predominantly products of affinity-matured B cells in the adult.

Phosphorylcholine (PC) is a classic T-independent Ag that is exposed on apoptotic cells, oxidized phospholipids, and bacterial polysaccharides. Experimental as well as epidemiological studies have over the past decade implicated Abs against PC (anti-PC) as anti-inflammatory and a strong protective factor in cardiovascular disease. Although clinically important, little is known about the development of anti-PC in humans. This study was conceived to dissect the human anti-PC repertoire and generate human mAbs. We designed a PC-specific probe to identify, isolate, and characterize PC-reactive B cells from 10 healthy individuals. The donors had all mounted somatically mutated Abs toward PC using a broad variety of Ig genes. PC-reactive B cells were primarily found in the IgM(+) memory subset, although significant numbers also were detected among naive, IgG(+), and CD27(+)CD43(+) B cells. Abs from these subsets were clonally related, suggesting a common origin. mAbs derived from the same donors exhibited equivalent or higher affinity for PC than the well-characterized murine T-15 clone. These results provide novel insights into the cellular and molecular ontogeny of atheroprotective PC Abs, thereby offering new opportunities for Ab-based therapeutic interventions.

Making sense of large-scale kinase inhibitor bioactivity data sets: a comparative and integrative analysis.

We carried out a systematic evaluation of target selectivity profiles across three recent large-scale biochemical assays of kinase inhibitors and further compared these standardized bioactivity assays with data reported in the widely used databases ChEMBL and STITCH. Our comparative evaluation revealed relative benefits and potential limitations among the bioactivity types, as well as pinpointed biases in the database curation processes. Ignoring such issues in data heterogeneity and representation may lead to biased modeling of drugs' polypharmacological effects as well as to unrealistic evaluation of computational strategies for the prediction of drug-target interaction networks. Toward making use of the complementary information captured by the various bioactivity types, including IC50, K(i), and K(d), we also introduce a model-based integration approach, termed KIBA, and demonstrate here how it can be used to classify kinase inhibitor targets and to pinpoint potential errors in database-reported drug-target interactions. An integrated drug-target bioactivity matrix across 52,498 chemical compounds and 467 kinase targets, including a total of 246,088 KIBA scores, has been made freely available.

Chk1 targeting reactivates PP2A tumor suppressor activity in cancer cells.

Checkpoint kinase Chk1 is constitutively active in many cancer cell types and new generation Chk1 inhibitors show marked antitumor activity as single agents. Here we present a hitherto unrecognized mechanism that contributes to the response of cancer cells to Chk1-targeted therapy. Inhibiting chronic Chk1 activity in cancer cells induced the tumor suppressor activity of protein phosphatase protein phosphatase 2A (PP2A), which by dephosphorylating MYC serine 62, inhibited MYC activity and impaired cancer cell survival. Mechanistic investigations revealed that Chk1 inhibition activated PP2A by decreasing the transcription of cancerous inhibitor of PP2A (CIP2A), a chief inhibitor of PP2A activity. Inhibition of cancer cell clonogenicity by Chk1 inhibition could be rescued in vitro either by exogenous expression of CIP2A or by blocking the CIP2A-regulated PP2A complex. Chk1-mediated CIP2A regulation was extended in tumor models dependent on either Chk1 or CIP2A. The clinical relevance of CIP2A as a Chk1 effector protein was validated in several human cancer types, including neuroblastoma, where CIP2A was identified as an NMYC-independent prognostic factor. Because the Chk1-CIP2A-PP2A pathway is driven by DNA-PK activity, functioning regardless of p53 or ATM/ATR status, our results offer explanative power for understanding how Chk1 inhibitors mediate single-agent anticancer efficacy. Furthermore, they define CIP2A-PP2A status in cancer cells as a pharmacodynamic marker for their response to Chk1-targeted therapy.

Individualized systems medicine strategy to tailor treatments for patients with chemorefractory acute myeloid leukemia.

UNLABELLED: We present an individualized systems medicine (ISM) approach to optimize cancer drug therapies one patient at a time. ISM is based on (i) molecular profiling and ex vivo drug sensitivity and resistance testing (DSRT) of patients' cancer cells to 187 oncology drugs, (ii) clinical implementation of therapies predicted to be effective, and (iii) studying consecutive samples from the treated patients to understand the basis of resistance. Here, application of ISM to 28 samples from patients with acute myeloid leukemia (AML) uncovered five major taxonomic drug-response subtypes based on DSRT profiles, some with distinct genomic features (e.g., MLL gene fusions in subgroup IV and FLT3-ITD mutations in subgroup V). Therapy based on DSRT resulted in several clinical responses. After progression under DSRT-guided therapies, AML cells displayed significant clonal evolution and novel genomic changes potentially explaining resistance, whereas ex vivo DSRT data showed resistance to the clinically applied drugs and new vulnerabilities to previously ineffective drugs. SIGNIFICANCE: Here, we demonstrate an ISM strategy to optimize safe and effective personalized cancer therapies for individual patients as well as to understand and predict disease evolution and the next line of therapy. This approach could facilitate systematic drug repositioning of approved targeted drugs as well as help to prioritize and de-risk emerging drugs for clinical testing.

Target inhibition networks: predicting selective combinations of druggable targets to block cancer survival pathways.

A recent trend in drug development is to identify drug combinations or multi-target agents that effectively modify multiple nodes of disease-associated networks. Such polypharmacological effects may reduce the risk of emerging drug resistance by means of attacking the disease networks through synergistic and synthetic lethal interactions. However, due to the exponentially increasing number of potential drug and target combinations, systematic approaches are needed for prioritizing the most potent multi-target alternatives on a global network level. We took a functional systems pharmacology approach toward the identification of selective target combinations for specific cancer cells by combining large-scale screening data on drug treatment efficacies and drug-target binding affinities. Our model-based prediction approach, named TIMMA, takes advantage of the polypharmacological effects of drugs and infers combinatorial drug efficacies through system-level target inhibition networks. Case studies in MCF-7 and MDA-MB-231 breast cancer and BxPC-3 pancreatic cancer cells demonstrated how the target inhibition modeling allows systematic exploration of functional interactions between drugs and their targets to maximally inhibit multiple survival pathways in a given cancer type. The TIMMA prediction results were experimentally validated by means of systematic siRNA-mediated silencing of the selected targets and their pairwise combinations, showing increased ability to identify not only such druggable kinase targets that are essential for cancer survival either individually or in combination, but also synergistic interactions indicative of non-additive drug efficacies. These system-level analyses were enabled by a novel model construction method utilizing maximization and minimization rules, as well as a model selection algorithm based on sequential forward floating search. Compared with an existing computational solution, TIMMA showed both enhanced prediction accuracies in cross validation as well as significant reduction in computation times. Such cost-effective computational-experimental design strategies have the potential to greatly speed-up the drug testing efforts by prioritizing those interventions and interactions warranting further study in individual cancer cases.

Drafting the CLN3 protein interactome in SH-SY5Y human neuroblastoma cells: a label-free quantitative proteomics approach.

Neuronal ceroid lipofuscinoses (NCL) are the most common inherited progressive encephalopathies of childhood. One of the most prevalent forms of NCL, Juvenile neuronal ceroid lipofuscinosis (JNCL) or CLN3 disease (OMIM: 204200), is caused by mutations in the CLN3 gene on chromosome 16p12.1. Despite progress in the NCL field, the primary function of ceroid-lipofuscinosis neuronal protein 3 (CLN3) remains elusive. In this study, we aimed to clarify the role of human CLN3 in the brain by identifying CLN3-associated proteins using a Tandem Affinity Purification coupled to Mass Spectrometry (TAP-MS) strategy combined with Significance Analysis of Interactome (SAINT). Human SH-SY5Y-NTAP-CLN3 stable cells were used to isolate native protein complexes for subsequent TAP-MS. Bioinformatic analyses of isolated complexes yielded 58 CLN3 interacting partners (IP) including 42 novel CLN3 IP, as well as 16 CLN3 high confidence interacting partners (HCIP) previously identified in another high-throughput study by Behrends et al., 2010. Moreover, 31 IP of ceroid-lipofuscinosis neuronal protein 5 (CLN5) were identified (18 of which were in common with the CLN3 bait). Our findings support previously suggested involvement of CLN3 in transmembrane transport, lipid homeostasis and neuronal excitability, as well as link it to G-protein signaling and protein folding/sorting in the ER.