Quantitative Proteomics with siRNA Screening Identifies Novel Mechanisms of Trastuzumab Resistance in HER2 Amplified Breast Cancers

HER2 is a receptor tyrosine kinase that is overexpressed in 20% to 30% of human breast cancers and which affects patient prognosis and survival. Treatment of HER2-positive breast cancer with the monoclonal antibody trastuzumab (Herceptin) has improved patient survival, but the development of trastuzumab resistance is a major medical problem. Many of the known mechanisms of trastuzumab resistance cause changes in protein phosphorylation patterns, and therefore quantitative proteomics was used to examine phosphotyrosine signaling networks in trastuzumab-resistant cells. The model system used in this study was two pairs of trastuzumab-sensitive and -resistant breast cancer cell lines. Using stable isotope labeling, phosphotyrosine immunoprecipitations, and online TiO₂ chromatography utilizing a dual trap configuration, ∼1700 proteins were quantified. Comparing quantified proteins between the two cell line pairs showed only a small number of common protein ratio changes, demonstrating heterogeneity in phosphotyrosine signaling networks across different trastuzumab-resistant cancers. Proteins showing significant increases in resistant versus sensitive cells were subjected to a focused siRNA screen to evaluate their functional relevance to trastuzumab resistance. The screen revealed proteins related to the Src kinase pathway, such as CDCP1/Trask, embryonal Fyn substrate, and Paxillin. We also identify several novel proteins that increased trastuzumab sensitivity in resistant cells when targeted by siRNAs, including FAM83A and MAPK1. These proteins may present targets for the development of clinical diagnostics or therapeutic strategies to guide the treatment of HER2+ breast cancer patients who develop trastuzumab resistance.HER2 is a member of the epidermal growth factor receptor (EGFR)/ErbB family of receptor tyrosine kinases. Under normal physiologic conditions, HER2 tyrosine kinase signaling is tightly regulated spatially and temporally by the requirement for it to heterodimerize with a ligand bound family member, such as EGFR, HER3/ErbB3, or HER4/ErbB4 (1). However, in 20% to 30% of human breast cancer cases, HER2 gene amplification is present, resulting in a high level of HER2 protein overexpression and unregulated, constitutive HER2 tyrosine kinase signaling (2, 3). HER2 gene amplified breast cancer, also termed HER2-positive breast cancer, carries a poor prognosis, but the development of the HER2 targeted monoclonal antibody trastuzumab (Herceptin) has significantly improved patient survival (2). Despite the clinical effectiveness of trastuzumab, the development of drug resistance significantly increases the risk of patient death. This poses a major medical problem, as most metastatic HER2-positive breast cancer patients develop trastuzumab resistance over the course of their cancer treatment (4). The treatment approach for HER2+ breast cancer patients after trastuzumab resistance develops is mostly a trial-and-error process that subjects the patient to increased toxicity. Therefore, there is a substantial medical need for strategies to overcome trastuzumab resistance.Multiple trastuzumab-resistance mechanisms have been identified, and they alter signaling networks and protein phosphorylation patterns in either a direct or an indirect manner. These mechanisms can be grouped into three categories. The first category is the activation of a parallel signaling network by other tyrosine kinases. These kinases include the receptor tyrosine kinases, EGFR, IGF1R, Her3, Met, EphA2, and Axl, as well as the erythropoietin-receptor-mediated activation of the cytoplasmic tyrosine kinases Jak2 and Src (5–11). The second category is the activation of downstream signaling proteins. Multiple studies have demonstrated activation of the phosphatidylinositol-3-kinase (PI3K)/AKT pathway in trastuzumab resistance, which occurs either via deletion of the PTEN lipid phosphatase or mutation of the PI3K genes (12, 13). Activation of Src family kinases or overexpression of cyclin E, which increases the cyclin E–cyclin-dependent kinase 2 signaling pathway, has also been reported (14). The third category includes mechanisms that maintain HER2 signaling even in the presence of trastuzumab. The production of a truncated isoform of HER2, p95HER2, which lacks the trastuzumab binding site, causes constitutive HER2 signaling (15, 16). Overexpression of the MUC4 sialomucin complex inhibits trastuzumab binding to HER2 and thereby maintains HER2 signaling (17, 18).Given that multiple trastuzumab-resistance mechanisms alter signaling networks and protein phosphorylation patterns, we reasoned that mapping phosphotyrosine signaling networks using quantitative proteomics would be a highly useful strategy for analyzing known mechanisms and identifying novel mechanisms of trastuzumab resistance. Quantitative proteomics and phosphotyrosine enrichment approaches have been extensively used to study the EGFR signal networks (19–23). We and others have used these approaches to map the HER2 signaling network (22, 24, 25). Multiple other tyrosine kinase signaling networks were analyzed using quantitative proteomics, including Ephrin receptor, EphB2 (26–28), platelet-derived growth factor receptor (PDGFR) (21), insulin receptor (29, 30), and the receptor for hepatocyte growth factor, c-MET (31).The goal of this study is to identify, quantify, and functionally screen proteins that might be involved in trastuzumab resistance. We used two pairs of HER2 gene amplified trastuzumab-sensitive (parental, SkBr3 and BT474) and -resistant (SkBr3^R and BT474^R) human breast cancer cell lines as models for trastuzumab resistance. These cell lines and their trastuzumab-resistant derivatives have been extensively characterized and highly cited in the breast cancer literature (32, 33). Using stable isotope labeling of amino acids in cell culture (SILAC),¹ phosphotyrosine immunoprecipitations, and online TiO₂ chromatography with dual trap configuration, we quantified the changes in phosphotyrosine containing proteins and interactors between trastuzumab-sensitive and -resistant cells. Several of the known trastuzumab-resistance mechanisms were identified, which serves as a positive control and validation of our approach, and large protein ratio changes were measured in proteins that had not been previously connected with trastuzumab resistance. We then performed a focused siRNA screen targeting the proteins with significantly increased protein ratios. This screen functionally tested the role of the identified proteins and identifies which proteins might have the largest effect on reversing trastuzumab resistance.     

Delineation of a Carcinogenic Helicobacter pylori Proteome

Helicobacter pylori is the strongest known risk factor for gastric adenocarcinoma, yet only a fraction of infected persons ever develop cancer. The extensive genetic diversity inherent to this pathogen has precluded comprehensive analyses of constituents that mediate carcinogenesis. We previously reported that in vivo adaptation of a non-carcinogenic H. pylori strain endowed the output derivative with the ability to induce adenocarcinoma, providing a unique opportunity to identify proteins selectively expressed by an oncogenic H. pylori strain. Using a global proteomics DIGE/MS approach, a novel missense mutation of the flagellar protein FlaA was identified that affects structure and function of this virulence-related organelle. Among 25 additional differentially abundant proteins, this approach also identified new proteins previously unassociated with gastric cancer, generating a profile of H. pylori proteins to use in vaccine development and for screening persons infected with strains most likely to induce severe disease.Helicobacter pylori is a Gram-negative bacterial species that selectively colonizes gastric epithelium and induces an inflammatory response within the stomach that persists for decades (1, 2). Biological costs incurred by the long term relationship between H. pylori and humans include an increased risk for distal gastric adenocarcinoma (3–8), and eradication of this pathogen significantly decreases cancer risk among infected individuals without premalignant lesions (9). However, only a fraction of colonized persons ever develop neoplasia, and enhanced cancer risk is related to H. pylori strain differences, inflammatory responses governed by host genetic diversity, and/or specific interactions between host and microbial determinants (10).H. pylori strains are remarkably diverse (11–15), and the genetic composition of strains can change over time within an individual colonized stomach (16, 17). Despite this diversity, several genetic loci have been identified that augment disease risk. The cag pathogenicity island encodes a type IV bacterial secretion system, and the product of the terminal gene in this island, CagA, is translocated into host epithelial cells by the cag secretion system following adherence (18–20). Within the host cell, CagA undergoes Src- and Abl-dependent tyrosine phosphorylation (21) and activates the eukaryotic phosphatase SHP-2, leading to dephosphorylation of host cell proteins and cellular morphological changes (19–21). CagA also dysregulates β-catenin signaling (22, 23) and apical-junctional complexes (24), events linked to increased cell motility and oncogenic transformation in several models (25, 26). Another H. pylori constituent linked to gastric cancer is the cytotoxin VacA, encoded by the gene vacA, which is present in virtually all H. pylori strains (27). In vitro, VacA induces the formation of intracellular vacuoles (27) and can induce apoptosis (28), and vacuolating activity is significantly associated with the presence of the cag pathogenicity island (3).Approximately 20% of H. pylori bind to gastric epithelial cells in vivo (29), and sequence analysis has revealed that the H. pylori genome contains an unusually high number of ORFs relative to its genome size that are predicted to encode outer membrane proteins (15). BabA, a member of a family of highly conserved outer membrane proteins and encoded by the strain-specific gene babA2, binds the Lewis^b histo-blood group antigen on gastric epithelial cells (30, 31), and H. pylori babA2⁺ strains are associated with an increased risk for gastric cancer (30). However, not all persons infected with cag⁺ babA2⁺ toxigenic strains develop gastric cancer, indicating that additional H. pylori constituents are important in carcinogenesis.We recently identified a strain of H. pylori, 7.13, that reproducibly induces gastric cancer in two rodent models of gastritis, Mongolian gerbils and hypergastrinemic INS-GAS mice (22). This strain was derived via in vivo adaptation of a clinical H. pylori strain, B128, which induces inflammation, but not cancer, in rodent gastric mucosa. The oncogenic 7.13 phenotype is not due to an enhanced ability of strain 7.13 to colonize as there were no significant differences in gastric colonization density or efficiency between strains B128 and 7.13 as assessed by either quantitative culture or histology. However, carcinogenic strain 7.13 binds more avidly to gastric epithelial cells in vitro than does strain B128, suggesting that the two strains may variably express different outer membrane proteins.To define proteins that may mediate the development of H. pylori-induced gastric cancer, we performed two-dimensional (2D)¹ DIGE coupled with MS to identify differentially abundant membrane-associated and cytosolic proteins from non-carcinogenic H. pylori strain B128 and its carcinogenic derivative, strain 7.13 (22). DIGE/MS is a well established proteomics technology based on conventional 2D gel protein separations whereby prelabeling samples with spectrally resolvable fluorescent dyes and multiplexing samples onto a series of gels that contain a mixture of all experimental samples (internal standard) provide quantitative data on abundance changes for thousands of intact proteins from multiple experimental conditions, each measured in replicate for statistical confidence (32–36). Techniques including DIGE/MS have recently been utilized to robustly define differences in protein abundance profiles between bacterial strains and to compare expression patterns of proteins harvested from bacteria maintained under different growth conditions (37, 38).Utilizing DIGE/MS, we detected and identified 26 proteins with statistically significant differences between strains B128 and 7.13, including a novel cysteine-to-arginine mutation in the H. pylori flagellar protein FlaA. We demonstrate that this FlaA mutation results in structural and functional aberrations. Application of this technique to two genetically related bacterial strains that induce distinct phenotypes also identified several novel candidate H. pylori virulence factors, providing a framework for studies targeting the pathogenesis of microbially induced cancer.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Identification of Tumor-associated Autoantigens for the Diagnosis of Colorectal Cancer in Serum Using High Density Protein Microarrays

There is a mounting evidence of the existence of autoantibodies associated to cancer progression. Antibodies are the target of choice for serum screening because of their stability and suitability for sensitive immunoassays. By using commercial protein microarrays containing 8000 human proteins, we examined 20 sera from colorectal cancer (CRC) patients and healthy subjects to identify autoantibody patterns and associated antigens. Forty-three proteins were differentially recognized by tumoral and reference sera (p value <0.04) in the protein microarrays. Five immunoreactive antigens, PIM1, MAPKAPK3, STK4, SRC, and FGFR4, showed the highest prevalence in cancer samples, whereas ACVR2B was more abundant in normal sera. Three of them, PIM1, MAPKAPK3, and ACVR2B, were used for further validation. A significant increase in the expression level of these antigens on CRC cell lines and colonic mucosa was confirmed by immunoblotting and immunohistochemistry on tissue microarrays. A diagnostic ELISA based on the combination of MAPKAPK3 and ACVR2B proteins yielded specificity and sensitivity values of 73.9 and 83.3% (area under the curve, 0.85), respectively, for CRC discrimination after using an independent sample set containing 94 sera representative of different stages of progression and control subjects. In summary, these studies confirmed the presence of specific autoantibodies for CRC and revealed new individual markers of disease (PIM1, MAPKAPK3, and ACVR2B) with the potential to diagnose CRC with higher specificity and sensitivity than previously reported serum biomarkers.Colorectal cancer (CRC)¹ is the second most prevalent cancer in the western world. The development of this disease takes decades and involves multiple genetic events. CRC remains a major cause of mortality in developed countries because most of the patients are diagnosed at advanced stages because of the reluctance to use highly invasive diagnostic tools like colonoscopy. Actually only a few proteins have been described as biomarkers in CRC (carcinoembryonic antigen (CEA), CA19.9, and CA125 (1–3)), although none of them is recommended for clinical screening (4). Proteomics analysis is actively used for the identification of new biomarkers. In previous studies, the use of two-dimensional DIGE and antibody microarrays allowed the identification of differentially expressed proteins in CRC tissue, including isoforms and post-translational modifications responsible for modifications in signaling pathways (5–8). Both approaches resulted in the identification of a collection of potential tumoral tissue biomarkers that is currently being investigated.However, the implementation of simpler, non-invasive methods for the early detection of CRC should be based on the identification of proteins or antibodies in serum or plasma (9–13). There is ample evidence of the existence of an immune response to cancer in humans as demonstrated by the presence of autoantibodies in cancer sera. Self-proteins (autoantigens) altered before or during tumor formation can elicit an immune response (13–19). These tumor-specific autoantibodies can be detected at early cancer stages and prior to cancer diagnosis revealing a great potential as biomarkers (14, 15, 20). Tumor proteins can be affected by specific point mutations, misfolding, overexpression, aberrant glycosylation, truncation, or aberrant degradation (e.g. p53, HER2, NY-ESO1, or MUC1 (16, 21–25)). In fact, a number of tumor-associated autoantigens (TAAs) were identified previously in different studies involving autoantibody screening in CRC (26–28).Several approaches have been used to identify TAAs in cancer, including natural protein arrays prepared with fractions obtained from two-dimensional LC separations of tumoral samples (29, 30) or protein extracts from cancer cells or tissue (9, 31) probed by Western blot with patient sera, cancer tissue peptide libraries expressed as cDNA expression libraries for serological screening (serological analysis of recombinant cDNA expression libraries (SEREX)) (22, 32), or peptides expressed on the surface of phages in combination with microarrays (17, 18, 33, 34). However, these approaches suffer from several drawbacks. In some cases TAAs have to be isolated and identified from the reactive protein lysate by LC-MS techniques, or in the phage display approach, the reactive TAA could be a mimotope without a corresponding linear amino acid sequence. Moreover, cDNA libraries might not be representative of the protein expression levels in tumors as there is a poor correspondence between mRNA and protein levels.Protein arrays provide a novel platform for the identification of both autoantibodies and their respective TAAs for diagnostic purposes in cancer serum patients. They present some advantages. Proteins printed on the microarray are known “a priori,” avoiding the need for later identifications and the discovery of mimotopes. There is no bias in protein selection as the proteins are printed at a similar concentration. This should result in a higher sensitivity for biomarker identification (13, 35, 36).In this study, we used commercially available high density protein microarrays for the identification of autoantibody signatures and tumor-associated antigens in colorectal cancer. We screened 20 CRC patient and control sera with protein microarrays containing 8000 human proteins to identify the CRC-associated autoantibody repertoire and the corresponding TAAs. Autoantibody profiles that discriminated the different types of CRC metastasis were identified. Moreover a set of TAAs showing increased or decreased expression in cancer cell lines and paired tumoral tissues was found. Finally an ELISA was set up to test the ability of the most immunoreactive proteins to detect colorectal adenocarcinoma. On the basis of the antibody response, combinations of three antigens, PIM1, MAPKAPK3, and ACVR2B, showed a great potential for diagnosis.     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]