The cancer cell proteome and transcriptome predicts sensitivity to targeted and cytotoxic drugs

Mapping the proteogenomic landscape enables prediction of drug response in acute myeloid leukemia

Acute myeloid leukemia is a poor-prognosis cancer commonly stratified by genetic aberrations, but these mutations are often heterogeneous and fail to consistently predict therapeutic response. Here, we combine transcriptomic, proteomic, and phosphoproteomic datasets with ex vivo drug sensitivity data to help understand the underlying pathophysiology of AML beyond mutations. We measure the proteome and phosphoproteome of 210 patients and combine them with genomic and transcriptomic measurements to identify four proteogenomic subtypes that complement existing genetic subtypes. We build a predictor to classify samples into subtypes and map them to a "landscape" that identifies specific drug response patterns. We then build a drug response prediction model to identify drugs that target distinct subtypes and validate our findings on cell lines representing various stages of quizartinib resistance. Our results show how multiomics data together with drug sensitivity data can inform therapy stratification and drug combinations in AML.

Proteomic and phosphoproteomic measurements enhance ability to predict ex vivo drug response in AML

Acute Myeloid Leukemia (AML) affects 20,000 patients in the US annually with a five-year survival rate of approximately 25%. One reason for the low survival rate is the high prevalence of clonal evolution that gives rise to heterogeneous sub-populations of leukemic cells with diverse mutation spectra, which eventually leads to disease relapse. This genetic heterogeneity drives the activation of complex signaling pathways that is reflected at the protein level. This diversity makes it difficult to treat AML with targeted therapy, requiring custom patient treatment protocols tailored to each individual's leukemia. Toward this end, the Beat AML research program prospectively collected genomic and transcriptomic data from over 1000 AML patients and carried out ex vivo drug sensitivity assays to identify genomic signatures that could predict patient-specific drug responses. However, there are inherent weaknesses in using only genetic and transcriptomic measurements as surrogates of drug response, particularly the absence of direct information about phosphorylation-mediated signal transduction. As a member of the Clinical Proteomic Tumor Analysis Consortium, we have extended the molecular characterization of this cohort by collecting proteomic and phosphoproteomic measurements from a subset of these patient samples (38 in total) to evaluate the hypothesis that proteomic signatures can improve the ability to predict response to 26 drugs in AML ex vivo samples. In this work we describe our systematic, multi-omic approach to evaluate proteomic signatures of drug response and compare protein levels to other markers of drug response such as mutational patterns. We explore the nuances of this approach using two drugs that target key pathways activated in AML: quizartinib (FLT3) and trametinib (Ras/MEK), and show how patient-derived signatures can be interpreted biologically and validated in cell lines. In conclusion, this pilot study demonstrates strong promise for proteomics-based patient stratification to assess drug sensitivity in AML.

Gene expression-based prediction of pazopanib efficacy in sarcoma

BACKGROUND

The multi-receptor tyrosine kinase inhibitor pazopanib is approved for the treatment of advanced soft-tissue sarcoma and has also shown activity in other sarcoma subtypes. However, its clinical efficacy is highly variable, and no reliable predictors exist to select patients who are likely to benefit from this drug.

PATIENTS AND METHODS

We analysed the molecular profiles and clinical outcomes of patients with pazopanib-treated sarcoma enrolled in a prospective observational study by the German Cancer Consortium, DKTK MASTER, that employs whole-genome/exome sequencing and transcriptome sequencing to inform the care of young adults with advanced cancer across histology and patients with rare cancers.

RESULTS

Among 109 patients with available whole-genome/exome sequencing data, there was no correlation between clinical parameters, specific genetic alterations or mutational signatures and clinical outcome. In contrast, the analysis of a subcohort of 62 patients who underwent molecular analysis before pazopanib treatment and had transcriptome sequencing data available showed that mRNA levels of NTRK3 (hazard ratio [HR] = 0.53, p = 0.021), IGF1R (HR = 1.82, p = 0.027) and KDR (HR = 0.50, p = 0.011) were independently associated with progression-free survival (PFS). Based on the expression of these multi-receptor tyrosine kinase genes, i.e. the features NTRK3-high, IGF1R-low and KDR-high, we developed a pazopanib efficacy predictor that stratified patients into three groups with significantly different PFS (p < 0.0001). Application of the pazopanib efficacy predictor to an independent cohort of patients with pazopanib-treated sarcoma from DKTK MASTER (n = 43) confirmed its potential to separate patient groups with significantly different PFS (p = 0.02), whereas no such association was observed in patients with sarcoma from DKTK MASTER (n = 97) or The Cancer Genome Atlas sarcoma cohort (n = 256) who were not treated with pazopanib.

CONCLUSION

A score based on the combined expression of NTRK3, IGF1R and KDR allows the identification of patients with sarcoma and with good, intermediate and poor outcome following pazopanib therapy and warrants prospective investigation as a predictive tool to optimise the use of this drug in the clinic.

MECHANISTIC AND DATA-DRIVEN MODELS OF CELL SIGNALING: TOOLS FOR FUNDAMENTAL DISCOVERY AND RATIONAL DESIGN OF THERAPY

A full understanding of cell signaling processes requires knowledge of protein structure/function relationships, protein-protein interactions, and the abilities of pathways to control phenotypes. Computational models offer a valuable framework for integrating that knowledge to predict the effects of system perturbations and interventions in health and disease. Whereas mechanistic models are well suited for understanding the biophysical basis for signal transduction and principles of therapeutic design, data-driven models are particularly suited to distill complex signaling relationships among samples and between multivariate signaling changes and phenotypes. Both approaches have limitations and provide incomplete representations of signaling biology, but their careful implementation and integration can provide new understanding for how manipulating system variables impacts cellular decisions.

Integrative molecular characterization of sarcomatoid and rhabdoid renal cell carcinoma

Sarcomatoid and rhabdoid (S/R) renal cell carcinoma (RCC) are highly aggressive tumors with limited molecular and clinical characterization. Emerging evidence suggests immune checkpoint inhibitors (ICI) are particularly effective for these tumors, although the biological basis for this property is largely unknown. Here, we evaluate multiple clinical trial and real-world cohorts of S/R RCC to characterize their molecular features, clinical outcomes, and immunologic characteristics. We find that S/R RCC tumors harbor distinctive molecular features that may account for their aggressive behavior, including BAP1 mutations, CDKN2A deletions, and increased expression of MYC transcriptional programs. We show that these tumors are highly responsive to ICI and that they exhibit an immune-inflamed phenotype characterized by immune activation, increased cytotoxic immune infiltration, upregulation of antigen presentation machinery genes, and PD-L1 expression. Our findings build on prior work and shed light on the molecular drivers of aggressivity and responsiveness to ICI of S/R RCC.

Sarcomatoid and rhabdoid tumours are highly aggressive forms of renal cell carcinoma that are also responsive to immunotherapy. In this study, the authors perform a comprehensive molecular characterization of these tumours discovering an enrichment of specific alterations and an inflamed phenotype.

A modular master regulator landscape controls cancer transcriptional identity

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

Discovering the anti-cancer potential of non-oncology drugs by systematic viability profiling

Anti-cancer uses of non-oncology drugs have occasionally been found, but such discoveries have been serendipitous. We sought to create a public resource containing the growth inhibitory activity of 4,518 drugs tested across 578 human cancer cell lines. We used PRISM, a molecular barcoding method, to screen drugs against cell lines in pools. An unexpectedly large number of non-oncology drugs selectively inhibited subsets of cancer cell lines in a manner predictable from the cell lines' molecular features. Our findings include compounds that killed by inducing PDE3A-SLFN12 complex formation; vanadium-containing compounds whose killing depended on the sulfate transporter SLC26A2; the alcohol dependence drug disulfiram, which killed cells with low expression of metallothioneins; and the anti-inflammatory drug tepoxalin, which killed via the multi-drug resistance protein ABCB1. The PRISM drug repurposing resource (https://depmap.org/repurposing) is a starting point to develop new oncology therapeutics, and more rarely, for potential direct clinical translation.

Discovering the anti-cancer potential of non-oncology drugs by systematic viability profiling

Anti-cancer uses of non-oncology drugs have occasionally been found, but such discoveries have been serendipitous. We sought to create a public resource containing the growth inhibitory activity of 4,518 drugs tested across 578 human cancer cell lines. We used PRISM, a molecular barcoding method, to screen drugs against cell lines in pools. An unexpectedly large number of non-oncology drugs selectively inhibited subsets of cancer cell lines in a manner predictable from the cell lines' molecular features. Our findings include compounds that killed by inducing PDE3A-SLFN12 complex formation; vanadium-containing compounds whose killing depended on the sulfate transporter SLC26A2; the alcohol dependence drug disulfiram, which killed cells with low expression of metallothioneins; and the anti-inflammatory drug tepoxalin, which killed via the multi-drug resistance protein ABCB1. The PRISM drug repurposing resource (https://depmap.org/repurposing) is a starting point to develop new oncology therapeutics, and more rarely, for potential direct clinical translation.