Systems analysis of protein signatures predicting cetuximab responses in KRAS, NRAS, BRAF and PIK3CA wild-type patient-derived xenograft models of metastatic colorectal cancer.

Antibodies targeting the human epidermal growth factor receptor (EGFR) are used for the treatment of RAS wild-type metastatic colorectal cancer. A significant proportion of patients remains unresponsive to this therapy. Here, we performed a reverse-phase protein array-based (phospho)protein analysis of 63 KRAS, NRAS, BRAF and PIK3CA wild-type metastatic CRC tumours. Responses of tumours to anti-EGFR therapy with cetuximab were recorded in patient-derived xenograft (PDX) models. Unsupervised hierarchical clustering of pretreatment tumour tissue identified three clusters, of which Cluster C3 was exclusively composed of responders. Clusters C1 and C2 exhibited mixed responses. None of the three protein clusters exhibited a significant correlation with transcriptome-based subtypes. Analysis of protein signatures across all PDXs identified 14 markers that discriminated cetuximab-sensitive and cetuximab-resistant tumours: PDK1 (S241), caspase-8, Shc (Y317), Stat3 (Y705), p27, GSK-3ß (S9), HER3, PKC-α (S657), EGFR (Y1068), Akt (S473), S6 ribosomal protein (S240/244), HER3 (Y1289), NF-κB-p65 (S536) and Gab-1 (Y627). Least absolute shrinkage and selection operator and binominal logistic regression analysis delivered refined protein signatures for predicting response to cetuximab. (Phospo-)protein analysis of matched pretreated and posttreated models furthermore showed significant reduction of Gab-1 (Y627) and GSK-3ß (S9) exclusively in responding models, suggesting novel targets for treatment.

BCL-2 system analysis identifies high-risk colorectal cancer patients.

OBJECTIVE: The mitochondrial apoptosis pathway is controlled by an interaction of multiple BCL-2 family proteins, and plays a key role in tumour progression and therapy responses. We assessed the prognostic potential of an experimentally validated, mathematical model of BCL-2 protein interactions (DR_MOMP) in patients with stage III colorectal cancer (CRC). DESIGN: Absolute protein levels of BCL-2 family proteins were determined in primary CRC tumours collected from n=128 resected and chemotherapy-treated patients with stage III CRC. We applied DR_MOMP to categorise patients as high or low risk based on model outputs, and compared model outputs with known prognostic factors (T-stage, N-stage, lymphovascular invasion). DR_MOMP signatures were validated on protein of n=156 patients with CRC from the Cancer Genome Atlas (TCGA) project. RESULTS: High-risk stage III patients identified by DR_MOMP had an approximately fivefold increased risk of death compared with patients identified as low risk (HR 5.2, 95% CI 1.4 to 17.9, p=0.02). The DR_MOMP signature ranked highest among all molecular and pathological features analysed. The prognostic signature was validated in the TCGA colon adenocarcinoma (COAD) cohort (HR 4.2, 95% CI 1.1 to 15.6, p=0.04). DR_MOMP also further stratified patients identified by supervised gene expression risk scores into low-risk and high-risk categories. BCL-2-dependent signalling critically contributed to treatment responses in consensus molecular subtypes 1 and 3, linking for the first time specific molecular subtypes to apoptosis signalling. CONCLUSIONS: DR_MOMP delivers a system-based biomarker with significant potential as a prognostic tool for stage III CRC that significantly improves established histopathological risk factors.

HGF induces epithelial-to-mesenchymal transition by modulating the mammalian hippo/MST2 and ISG15 pathways.

Epithelial to mesenchymal transition (EMT) is a fundamental cell differentiation/dedifferentiation process which is associated with dramatic morphological changes. Formerly polarized and immobile epithelial cells which form cell junctions and cobblestone-like cell sheets undergo a transition into highly motile, elongated, mesenchymal cells lacking cell-to-cell adhesions. To explore how the proteome is affected during EMT we profiled protein expression and tracked cell biological markers in Madin-Darby kidney epithelial cells undergoing hepatocyte growth factor (HGF) induced EMT. We were able to identify and quantify over 4000 proteins by mass spectrometry. Enrichment analysis of this revealed that expression of proteins associated with the ubiquitination machinery was induced, whereas expression of proteins regulating apoptotic pathways was suppressed. We show that both the mammalian Hippo/MST2 and the ISG15 pathways are regulated at the protein level by ubiquitin ligases. Inhibition of the Hippo pathway by overexpression of either ITCH or A-Raf promotes HGF-induced EMT. Conversely, ISG15 overexpression is sufficient to induce cell scattering and an elongated morphology without external stimuli. Thus, we demonstrate for the first time that the Hippo/MST2 and ISG15 pathways are regulated during growth-factor induced EMT.

5'ValCAC tRNA fragment generated as part of a protective angiogenin response provides prognostic value in amyotrophic lateral sclerosis.

Loss-of-function mutations in the ribonuclease angiogenin are associated with amyotrophic lateral sclerosis. Angiogenin has been shown to cleave transfer RNAs during stress to produce 'transfer-derived stress-induced RNAs'. Stress-induced tRNA cleavage is preserved from single-celled organisms to humans indicating it represents part of a highly conserved stress response. However, to date, the role of tRNA cleavage in amyotrophic lateral sclerosis remains to be fully elucidated. To this end, we performed small RNA sequencing on a human astrocytoma cell line to identify the complete repertoire of tRNA fragments generated by angiogenin. We found that only a specific subset of tRNAs is cleaved by angiogenin and identified 5'ValCAC transfer-derived stress-induced RNA to be secreted from neural cells. 5'ValCAC was quantified in spinal cord and serum from SOD1G93A amyotrophic lateral sclerosis mouse models where we found it to be significantly elevated at symptom onset correlating with increased angiogenin expression, imbalanced protein translation initiation factors and slower disease progression. In amyotrophic lateral sclerosis patient serum samples, we found 5'ValCAC to be significantly higher in patients with slow disease progression, and interestingly, we find 5'ValCAC to hold prognostic value for amyotrophic lateral sclerosis patients. Here, we report that angiogenin cleaves a specific subset of tRNAs and provide evidence for 5'ValCAC as a prognostic biomarker in amyotrophic lateral sclerosis. We propose that increased serum 5'ValCAC levels indicate an enhanced angiogenin-mediated stress response within motor neurons that correlates with increased survival. These data suggest that the previously reported beneficial effects of angiogenin in SOD1G93A mice may result from elevated levels of 5'ValCAC transfer RNA fragment.

Elevation in plasma tRNA fragments precede seizures in human epilepsy.

Transfer RNAs (tRNAs) are a major class of noncoding RNA. Stress-induced cleavage of tRNA is highly conserved and results in tRNA fragments. Here we find specific tRNA fragments in plasma are associated with epilepsy. Small RNA sequencing of plasma samples collected during video-EEG monitoring of focal epilepsy patients identified significant differences in three tRNA fragments (5', 5'AlaTGC, and 5'GluCTC) from controls. Levels of these tRNA fragments were higher in pre-seizure than post-seizure samples, suggesting they may serve as biomarkers of seizure risk in epilepsy patients. In vitro studies confirmed that production and extracellular release of tRNA fragments was lower after epileptiform-like activity in hippocampal neurons. We designed PCR-based assays to quantify tRNA fragments in a cohort of pre- and post-seizure plasma samples from focal epilepsy patients and healthy controls (n = 32/group). Receiver operating characteristic analysis indicated that tRNA fragments potently distinguished pre- from post-seizure patients (area under the curve of 0.8-0.95). Elevated tRNA fragments levels were not detected in patients with psychogenic non-epileptic seizures, and did not result from medication tapering. This study identifies a novel class of epilepsy biomarker and reveals the potential existence of prodromal molecular patterns in blood that could be used to predict seizure risk.

Implementing Reverse Phase Protein Array Profiling as a Sensitive Method for the Early Pre-Clinical Detection of Off-Target Toxicities Associated with Sunitinib Malate.

PURPOSE: The tyrosine kinase inhibitor (TKI) sunitinib is a multi-targeted agent approved across multiple cancer indications. Nevertheless, since approval, data has emerged to describe a worrisome side effect profile including hypertension, hand-foot syndrome, fatigue, diarrhea, mucositis, proteinuria, and (rarely) congestive heart failure. It has been hypothesized that the observed multi-parameter toxicity profile is related to "on-target" kinase inhibition in "off-target" tissues. EXPERIMENTAL DESIGN: To interrogate off-target effects in pre-clinical studies, a reverse phase protein array (RPPA) approach is employed. Mice are treated with sunitinib (40 mg kg-1 ) for 4 weeks, following which critical organs are removed. The Zeptosens RPPA platform is employed for protein expression analysis. RESULTS: Differentially expressed proteins associated with damage and/or stress are found in the majority of organs from treated animals. Proteins differentially expressed in the heart are associated with myocardial hypertrophy, ischaemia/reperfusion, and hypoxia. However, hypertrophy is not evidenced on histology. Mild proteinuria is observed; however, no changes in renal glomerular structure are visible via electron microscopy. In skin, proteins associated with cutaneous inflammation, keratinocyte hyper-proliferation, and increased inflammatory response are differentially expressed. CONCLUSIONS AND CLINICAL RELEVANCE: It is posited that pre-clinical implementation of a combined histopathological/RPPA approach provides a sensitive method to mechanistically elucidate the early manifestation of TKI on-target/organ off-target toxicities.

The BAX/BAK-like protein BOK is a prognostic marker in colorectal cancer.

The intrinsic or mitochondrial apoptosis pathway is controlled by the interaction of antiapoptotic and pro-apoptotic members of the BCL-2 protein family. Activation of this death pathway plays a crucial role in cancer progression and chemotherapy responses. The BCL-2-related ovarian killer (BOK) possesses three BCL-2 homology domains and has been proposed to act in a similar pro-apoptotic pathway as the pro-apoptotic proteins BAX and BAK. In this study, we showed that stage II and III colorectal cancer patients possessed decreased levels of BOK protein in their tumours compared to matched normal tissue. BOK protein levels in tumours were also prognostic of clinical outcome but increased BOK protein levels surprisingly associated with earlier disease recurrence and reduced overall survival. We found no significant association of BOK protein tumour levels with ER stress markers GRP78 or GRP94 or with cleaved caspase-3. In contrast, BOK protein levels correlated with Calreticulin. These data indicate BOK as a prognostic marker in colorectal cancer and suggest that different activities of BOK may contribute to cancer progression and prognosis.

Computational Analysis of AMPK-Mediated Neuroprotection Suggests Acute Excitotoxic Bioenergetics and Glucose Dynamics Are Regulated by a Minimal Set of Critical Reactions.

Loss of ionic homeostasis during excitotoxic stress depletes ATP levels and activates the AMP-activated protein kinase (AMPK), re-establishing energy production by increased expression of glucose transporters on the plasma membrane. Here, we develop a computational model to test whether this AMPK-mediated glucose import can rapidly restore ATP levels following a transient excitotoxic insult. We demonstrate that a highly compact model, comprising a minimal set of critical reactions, can closely resemble the rapid dynamics and cell-to-cell heterogeneity of ATP levels and AMPK activity, as confirmed by single-cell fluorescence microscopy in rat primary cerebellar neurons exposed to glutamate excitotoxicity. The model further correctly predicted an excitotoxicity-induced elevation of intracellular glucose, and well resembled the delayed recovery and cell-to-cell heterogeneity of experimentally measured glucose dynamics. The model also predicted necrotic bioenergetic collapse and altered calcium dynamics following more severe excitotoxic insults. In conclusion, our data suggest that a minimal set of critical reactions may determine the acute bioenergetic response to transient excitotoxicity and that an AMPK-mediated increase in intracellular glucose may be sufficient to rapidly recover ATP levels following an excitotoxic insult.

Substrate-Trapped Interactors of PHD3 and FIH Cluster in Distinct Signaling Pathways.

Amino acid hydroxylation is a post-translational modification that regulates intra- and inter-molecular protein-protein interactions. The modifications are regulated by a family of 2-oxoglutarate- (2OG) dependent enzymes and, although the biochemistry is well understood, until now only a few substrates have been described for these enzymes. Using quantitative interaction proteomics, we screened for substrates of the proline hydroxylase PHD3 and the asparagine hydroxylase FIH, which regulate the HIF-mediated hypoxic response. We were able to identify hundreds of potential substrates. Enrichment analysis revealed that the potential substrates of both hydroxylases cluster in the same pathways but frequently modify different nodes of signaling networks. We confirm that two proteins identified in our screen, MAPK6 (Erk3) and RIPK4, are indeed hydroxylated in a FIH- or PHD3-dependent mechanism. We further determined that FIH-dependent hydroxylation regulates RIPK4-dependent Wnt signaling, and that PHD3-dependent hydroxylation of MAPK6 protects the protein from proteasomal degradation.

Bistability in the Rac1, PAK, and RhoA Signaling Network Drives Actin Cytoskeleton Dynamics and Cell Motility Switches.

Dynamic interactions between RhoA and Rac1, members of the Rho small GTPase family, play a vital role in the control of cell migration. Using predictive mathematical modeling, mass spectrometry-based quantitation of network components, and experimental validation in MDA-MB-231 mesenchymal breast cancer cells, we show that a network containing Rac1, RhoA, and PAK family kinases can produce bistable, switch-like responses to a graded PAK inhibition. Using a small chemical inhibitor of PAK, we demonstrate that cellular RhoA and Rac1 activation levels respond in a history-dependent, bistable manner to PAK inhibition. Consequently, we show that downstream signaling, actin dynamics, and cell migration also behave in a bistable fashion, displaying switches and hysteresis in response to PAK inhibition. Our results demonstrate that PAK is a critical component in the Rac1-RhoA inhibitory crosstalk that governs bistable GTPase activity, cell morphology, and cell migration switches.

Network-based identification of feedback modules that control RhoA activity and cell migration.

Cancer cell migration enables metastatic spread causing most cancer deaths. Rho-family GTPases control cell migration, but being embedded in a highly interconnected feedback network, the control of their dynamical behavior during cell migration remains elusive. To address this question, we reconstructed the Rho-family GTPases signaling network involved in cell migration, and developed a Boolean network model to analyze the different states and emergent rewiring of the Rho-family GTPases signaling network at protrusions and during extracellular matrix-dependent cell migration. Extensive simulations and experimental validations revealed that the bursts of RhoA activity induced at protrusions by EGF are regulated by a negative-feedback module composed of Src, FAK, and CSK. Interestingly, perturbing this module interfered with cyclic Rho activation and extracellular matrix-dependent migration, suggesting that CSK inhibition can be a novel and effective intervention strategy for blocking extracellular matrix-dependent cancer cell migration, while Src inhibition might fail, depending on the genetic background of cells. Thus, this study provides new insights into the mechanisms that regulate the intricate activation states of Rho-family GTPases during extracellular matrix-dependent migration, revealing potential new targets for interfering with extracellular matrix-dependent cancer cell migration.

Searching for synergies: matrix algebraic approaches for efficient pair screening.

Functionally interacting perturbations, such as synergistic drugs pairs or synthetic lethal gene pairs, are of key interest in both pharmacology and functional genomics. However, to find such pairs by traditional screening methods is both time consuming and costly. We present a novel computational-experimental framework for efficient identification of synergistic target pairs, applicable for screening of systems with sizes on the order of current drug, small RNA or SGA (Synthetic Genetic Array) libraries (>1000 targets). This framework exploits the fact that the response of a drug pair in a given system, or a pair of genes' propensity to interact functionally, can be partly predicted by computational means from (i) a small set of experimentally determined target pairs, and (ii) pre-existing data (e.g. gene ontology, PPI) on the similarities between targets. Predictions are obtained by a novel matrix algebraic technique, based on cyclical projections onto convex sets. We demonstrate the efficiency of the proposed method using drug-drug interaction data from seven cancer cell lines and gene-gene interaction data from yeast SGA screens. Our protocol increases the rate of synergism discovery significantly over traditional screening, by up to 7-fold. Our method is easy to implement and could be applied to accelerate pair screening for both animal and microbial systems.

Extracellular signal-regulated kinase regulates RhoA activation and tumor cell plasticity by inhibiting guanine exchange factor H1 activity.

In certain Ras mutant cell lines, the inhibition of extracellular signal-regulated kinase (ERK) signaling increases RhoA activity and inhibits cell motility, which was attributed to a decrease in Fra-1 levels. Here we report a Fra-1-independent augmentation of RhoA signaling during short-term inhibition of ERK signaling. Using mass spectrometry-based proteomics, we identified guanine exchange factor H1 (GEF-H1) as mediating this effect. ERK binds to the Rho exchange factor GEF-H1 and phosphorylates it on S959, causing inhibition of GEF-H1 activity and a consequent decrease in RhoA activity. Knockdown experiments and expression of a nonphosphorylatable S959A GEF-H1 mutant showed that this site is crucial in regulating cell motility and invasiveness. Thus, we identified GEF-H1 as a critical ERK effector that regulates motility, cell morphology, and invasiveness.

Comparative drug pair screening across multiple glioblastoma cell lines reveals novel drug-drug interactions.

BACKGROUND: Glioblastoma multiforme (GBM) is the most aggressive brain tumor in adults, and despite state-of-the-art treatment, survival remains poor and novel therapeutics are sorely needed. The aim of the present study was to identify new synergistic drug pairs for GBM. In addition, we aimed to explore differences in drug-drug interactions across multiple GBM-derived cell cultures and predict such differences by use of transcriptional biomarkers. METHODS: We performed a screen in which we quantified drug-drug interactions for 465 drug pairs in each of the 5 GBM cell lines U87MG, U343MG, U373MG, A172, and T98G. Selected interactions were further tested using isobole-based analysis and validated in 5 glioma-initiating cell cultures. Furthermore, drug interactions were predicted using microarray-based transcriptional profiling in combination with statistical modeling. RESULTS: Of the 5 × 465 drug pairs, we could define a subset of drug pairs with strong interaction in both standard cell lines and glioma-initiating cell cultures. In particular, a subset of pairs involving the pharmaceutical compounds rimcazole, sertraline, pterostilbene, and gefitinib showed a strong interaction in a majority of the cell cultures tested. Statistical modeling of microarray and interaction data using sparse canonical correlation analysis revealed several predictive biomarkers, which we propose could be of importance in regulating drug pair responses. CONCLUSION: We identify novel candidate drug pairs for GBM and suggest possibilities to prospectively use transcriptional biomarkers to predict drug interactions in individual cases.