Laura M. Heiser

Basal Subtype and MAPK/ERK Kinase (MEK)-Phosphoinositide 3-Kinase Feedback Signaling Determine Susceptibility of Breast Cancer Cells to MEK Inhibition

Specific inhibitors of mitogen-activated protein kinase/extracellular signal-regulated kinase (ERK) kinase (MEK) have been developed that efficiently inhibit the oncogenic RAF-MEK-ERK pathway. We used a systems-based approach to identify breast cancer subtypes particularly susceptible to MEK inhibitors and to understand molecular mechanisms conferring resistance to such compounds. Basal-type breast cancer cells were found to be particularly susceptible to growth inhibition by small-molecule MEK inhibitors. Activation of the phosphatidylinositol 3-kinase (PI3K) pathway in response to MEK inhibition through a negative MEK-epidermal growth factor receptor-PI3K feedback loop was found to limit efficacy. Interruption of this feedback mechanism by targeting MEK and PI3K produced synergistic effects, including induction of apoptosis and, in some cell lines, cell cycle arrest and protection from apoptosis induced by proapoptotic agents. These findings enhance our understanding of the interconnectivity of oncogenic signal transduction circuits and have implications for the design of future clinical trials of MEK inhibitors in breast cancer by guiding patient selection and suggesting rational combination therapies.

Subtype and pathway specific responses to anticancer compounds in breast cancer

Breast cancers are comprised of molecularly distinct subtypes that may respond differently to pathway-targeted therapies now under development. Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways, and drug response. In a test of 77 therapeutic compounds, nearly all drugs showed differential responses across these cell lines, and approximately one third showed subtype-, pathway-, and/or genomic aberration-specific responses. These observations suggest mechanisms of response and resistance and may inform efforts to develop molecular assays that predict clinical response.

A community effort to assess and improve drug sensitivity prediction algorithms

Predicting the best treatment strategy from genomic information is a core goal of precision medicine. Here we focus on predicting drug response based on a cohort of genomic, epigenomic and proteomic profiling data sets measured in human breast cancer cell lines. Through a collaborative effort between the National Cancer Institute (NCI) and the Dialogue on Reverse Engineering Assessment and Methods (DREAM) project, we analyzed a total of 44 drug sensitivity prediction algorithms. The top-performing approaches modeled nonlinear relationships and incorporated biological pathway information. We found that gene expression microarrays consistently provided the best predictive power of the individual profiling data sets; however, performance was increased by including multiple, independent data sets. We discuss the innovations underlying the top-performing methodology, Bayesian multitask MKL, and we provide detailed descriptions of all methods. This study establishes benchmarks for drug sensitivity prediction and identifies approaches that can be leveraged for the development of new methods.

Metrics other than potency reveal systematic variation in responses to cancer drugs

Large-scale analysis of cellular response to anticancer drugs typically focuses on variation in potency (half-maximum inhibitory concentration, (IC50)), assuming that it is the most important difference between effective and ineffective drugs or sensitive and resistant cells. We took a multiparametric approach involving analysis of the slope of the dose-response curve, the area under the curve and the maximum effect (Emax). We found that some of these parameters vary systematically with cell line and others with drug class. For cell-cycle inhibitors, Emax often but not always correlated with cell proliferation rate. For drugs targeting the Akt/PI3K/mTOR pathway, dose-response curves were unusually shallow. Classical pharmacology has no ready explanation for this phenomenon, but single-cell analysis showed that it correlated with significant and heritable cell-to-cell variability in the extent of target inhibition. We conclude that parameters other than potency should be considered in the comparative analysis of drug response, particularly at clinically relevant concentrations near and above the IC50.

Modeling precision treatment of breast cancer

BackgroundFirst-generation molecular profiles for human breast cancers have enabled the identification of features that can predict therapeutic response; however, little is known about how the various data types can best be combined to yield optimal predictors. Collections of breast cancer cell lines mirror many aspects of breast cancer molecular pathobiology, and measurements of their omic and biological therapeutic responses are well-suited for development of strategies to identify the most predictive molecular feature sets.ResultsWe used least squares-support vector machines and random forest algorithms to identify molecular features associated with responses of a collection of 70 breast cancer cell lines to 90 experimental or approved therapeutic agents. The datasets analyzed included measurements of copy number aberrations, mutations, gene and isoform expression, promoter methylation and protein expression. Transcriptional subtype contributed strongly to response predictors for 25% of compounds, and adding other molecular data types improved prediction for 65%. No single molecular dataset consistently out-performed the others, suggesting that therapeutic response is mediated at multiple levels in the genome. Response predictors were developed and applied to TCGA data, and were found to be present in subsets of those patient samples.ConclusionsThese results suggest that matching patients to treatments based on transcriptional subtype will improve response rates, and inclusion of additional features from other profiling data types may provide additional benefit. Further, we suggest a systems biology strategy for guiding clinical trials so that patient cohorts most likely to respond to new therapies may be more efficiently identified.

A community computational challenge to predict the activity of pairs of compounds

Recent therapeutic successes have renewed interest in drug combinations, but experimental screening approaches are costly and often identify only small numbers of synergistic combinations. The DREAM consortium launched an open challenge to foster the development of in silico methods to computationally rank 91 compound pairs, from the most synergistic to the most antagonistic, based on gene-expression profiles of human B cells treated with individual compounds at multiple time points and concentrations. Using scoring metrics based on experimental dose-response curves, we assessed 32 methods (31 community-generated approaches and SynGen), four of which performed significantly better than random guessing. We highlight similarities between the methods. Although the accuracy of predictions was not optimal, we find that computational prediction of compound-pair activity is possible, and that community challenges can be useful to advance the field of in silico compound-synergy prediction.

Bayesian Inference of Signaling Network Topology in a Cancer Cell Line

MOTIVATION Protein signaling networks play a key role in cellular function, and their dysregulation is central to many diseases, including cancer. To shed light on signaling network topology in specific contexts, such as cancer, requires interrogation of multiple proteins through time and statistical approaches to make inferences regarding network structure. RESULTS In this study, we use dynamic Bayesian networks to make inferences regarding network structure and thereby generate testable hypotheses. We incorporate existing biology using informative network priors, weighted objectively by an empirical Bayes approach, and exploit a connection between variable selection and network inference to enable exact calculation of posterior probabilities of interest. The approach is computationally efficient and essentially free of user-set tuning parameters. Results on data where the true, underlying network is known place the approach favorably relative to existing approaches. We apply these methods to reverse-phase protein array time-course data from a breast cancer cell line (MDA-MB-468) to predict signaling links that we independently validate using targeted inhibition. The methods proposed offer a general approach by which to elucidate molecular networks specific to biological context, including, but not limited to, human cancers. AVAILABILITY http://mukherjeelab.nki.nl/DBN (code and data).

Inferring causal molecular networks: empirical assessment through a community-based effort

It remains unclear whether causal, rather than merely correlational, relationships in molecular networks can be inferred in complex biological settings. Here we describe the HPN-DREAM network inference challenge, which focused on learning causal influences in signaling networks. We used phosphoprotein data from cancer cell lines as well as in silico data from a nonlinear dynamical model. Using the phosphoprotein data, we scored more than 2,000 networks submitted by challenge participants. The networks spanned 32 biological contexts and were scored in terms of causal validity with respect to unseen interventional data. A number of approaches were effective, and incorporating known biology was generally advantageous. Additional sub-challenges considered time-course prediction and visualization. Our results suggest that learning causal relationships may be feasible in complex settings such as disease states. Furthermore, our scoring approach provides a practical way to empirically assess inferred molecular networks in a causal sense.

Tumor-derived Cell Lines as Molecular Models of Cancer Pharmacogenomics

Compared with normal cells, tumor cells have undergone an array of genetic and epigenetic alterations. Often, these changes underlie cancer development, progression, and drug resistance, so the utility of model systems rests on their ability to recapitulate the genomic aberrations observed in primary tumors. Tumor-derived cell lines have long been used to study the underlying biologic processes in cancer, as well as screening platforms for discovering and evaluating the efficacy of anticancer therapeutics. Multiple -omic measurements across more than a thousand cancer cell lines have been produced following advances in high-throughput technologies and multigroup collaborative projects. These data complement the large, international cancer genomic sequencing efforts to characterize patient tumors, such as The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). Given the scope and scale of data that have been generated, researchers are now in a position to evaluate the similarities and differences that exist in genomic features between cell lines and patient samples. As pharmacogenomics models, cell lines offer the advantages of being easily grown, relatively inexpensive, and amenable to high-throughput testing of therapeutic agents. Data generated from cell lines can then be used to link cellular drug response to genomic features, where the ultimate goal is to build predictive signatures of patient outcome. This review highlights the recent work that has compared -omic profiles of cell lines with primary tumors, and discusses the advantages and disadvantages of cancer cell lines as pharmacogenomic models of anticancer therapies. Mol Cancer Res; 14(1); 3–13. ©2015 AACR.

Integrated analysis of breast cancer cell lines reveals unique signaling pathways

BackgroundCancer is a heterogeneous disease resulting from the accumulation of genetic defects that negatively impact control of cell division, motility, adhesion and apoptosis. Deregulation in signaling along the EgfR-MAPK pathway is common in breast cancer, though the manner in which deregulation occurs varies between both individuals and cancer subtypes.ResultsWe were interested in identifying subnetworks within the EgfR-MAPK pathway that are similarly deregulated across subsets of breast cancers. To that end, we mapped genomic, transcriptional and proteomic profiles for 30 breast cancer cell lines onto a curated Pathway Logic symbolic systems model of EgfR-MAPK signaling. This model was composed of 539 molecular states and 396 rules governing signaling between active states. We analyzed these models and identified several subtype-specific subnetworks, including one that suggested Pak1 is particularly important in regulating the MAPK cascade when it is over-expressed. We hypothesized that Pak1 over-expressing cell lines would have increased sensitivity to Mek inhibitors. We tested this experimentally by measuring quantitative responses of 20 breast cancer cell lines to three Mek inhibitors. We found that Pak1 over-expressing luminal breast cancer cell lines are significantly more sensitive to Mek inhibition compared to those that express Pak1 at low levels. This indicates that Pak1 over-expression may be a useful clinical marker to identify patient populations that may be sensitive to Mek inhibitors.ConclusionsAll together, our results support the utility of symbolic system biology models for identification of therapeutic approaches that will be effective against breast cancer subsets.