Activation of the Anaphase Promoting Complex Reverses Multiple Drug Resistant Cancer in a Canine Model of Multiple Drug Resistant Lymphoma.

Like humans, canine lymphomas are treated by chemotherapy cocktails and frequently develop multiple drug resistance (MDR). Their shortened clinical timelines and tumor accessibility make canines excellent models to study MDR mechanisms. Insulin-sensitizers have been shown to reduce the incidence of cancer in humans prescribed them, and we previously demonstrated that they also reverse and delay MDR development in vitro. Here, we treated canines with MDR lymphoma with metformin to assess clinical and tumoral responses, including changes in MDR biomarkers, and used mRNA microarrays to determine differential gene expression. Metformin reduced MDR protein markers in all canines in the study. Microarrays performed on mRNAs gathered through longitudinal tumor sampling identified a 290 gene set that was enriched in Anaphase Promoting Complex (APC) substrates and additional mRNAs associated with slowed mitotic progression in MDR samples compared to skin controls. mRNAs from a canine that went into remission showed that APC substrate mRNAs were decreased, indicating that the APC was activated during remission. In vitro validation using canine lymphoma cells selected for resistance to chemotherapeutic drugs confirmed that APC activation restored MDR chemosensitivity, and that APC activity was reduced in MDR cells. This supports the idea that rapidly pushing MDR cells that harbor high loads of chromosome instability through mitosis, by activating the APC, contributes to improved survival and disease-free duration.

Synthetic lethal interactions of DEAD/H-box helicases as targets for cancer therapy.

DEAD/H-box helicases are implicated in virtually every aspect of RNA metabolism, including transcription, pre-mRNA splicing, ribosomes biogenesis, nuclear export, translation initiation, RNA degradation, and mRNA editing. Most of these helicases are upregulated in various cancers and mutations in some of them are associated with several malignancies. Lately, synthetic lethality (SL) and synthetic dosage lethality (SDL) approaches, where genetic interactions of cancer-related genes are exploited as therapeutic targets, are emerging as a leading area of cancer research. Several DEAD/H-box helicases, including DDX3, DDX9 (Dbp9), DDX10 (Dbp4), DDX11 (ChlR1), and DDX41 (Sacy-1), have been subjected to SL analyses in humans and different model organisms. It remains to be explored whether SDL can be utilized to identity druggable targets in DEAD/H-box helicase overexpressing cancers. In this review, we analyze gene expression data of a subset of DEAD/H-box helicases in multiple cancer types and discuss how their SL/SDL interactions can be used for therapeutic purposes. We also summarize the latest developments in clinical applications, apart from discussing some of the challenges in drug discovery in the context of targeting DEAD/H-box helicases.

PIIKA 2.5: Enhanced quality control of peptide microarrays for kinome analysis.

Peptide microarrays consisting of defined phosphorylation target sites are an effective approach for high throughput analysis of cellular kinase (kinome) activity. Kinome peptide arrays are highly customizable and do not require species-specific reagents to measure kinase activity, making them amenable for kinome analysis in any species. Our group developed software, Platform for Integrated, Intelligent Kinome Analysis (PIIKA), to enable more effective extraction of meaningful biological information from kinome peptide array data. A subsequent version, PIIKA2, unveiled new statistical tools and data visualization options. Here we introduce PIIKA 2.5 to provide two essential quality control metrics and a new background correction technique to increase the accuracy and consistency of kinome results. The first metric alerts users to improper spot size and informs them of the need to perform manual resizing to enhance the quality of the raw intensity data. The second metric uses inter-array comparisons to identify outlier arrays that sometimes emerge as a consequence of technical issues. In addition, a new background correction method, background scaling, can sharply reduce spatial biases within a single array in comparison to background subtraction alone. Collectively, the modifications of PIIKA 2.5 enable identification and correction of technical issues inherent to the technology and better facilitate the extraction of meaningful biological information. We show that these metrics demonstrably enhance kinome analysis by identifying low quality data and reducing batch effects, and ultimately improve clustering of treatment groups and enhance reproducibility. The web-based and stand-alone versions of PIIKA 2.5 are freely accessible at via http://saphire.usask.ca.

The KH domain facilitates the substrate specificity and unwinding processivity of DDX43 helicase.

The K-homology (KH) domain is a nucleic acid-binding domain present in many proteins. Recently, we found that the DEAD-box helicase DDX43 contains a KH domain in its N-terminus; however, its function remains unknown. Here, we purified recombinant DDX43 KH domain protein and found that it prefers binding ssDNA and ssRNA. Electrophoretic mobility shift assay and NMR revealed that the KH domain favors pyrimidines over purines. Mutational analysis showed that the GXXG loop in the KH domain is involved in pyrimidine binding. Moreover, we found that an alanine residue adjacent to the GXXG loop is critical for binding. Systematic evolution of ligands by exponential enrichment, chromatin immunoprecipitation-seq, and cross-linking immunoprecipitation-seq showed that the KH domain binds C-/T-rich DNA and U-rich RNA. Bioinformatics analysis suggested that the KH domain prefers to bind promoters. Using 15N-heteronuclear single quantum coherence NMR, the optimal binding sequence was identified as TTGT. Finally, we found that the full-length DDX43 helicase prefers DNA or RNA substrates with TTGT or UUGU single-stranded tails and that the KH domain is critically important for sequence specificity and unwinding processivity. Collectively, our results demonstrated that the KH domain facilitates the substrate specificity and processivity of the DDX43 helicase.

MDIPA: a microRNA-drug interaction prediction approach based on non-negative matrix factorization.

MOTIVATION: Evidence has shown that microRNAs, one type of small biomolecule, regulate the expression level of genes and play an important role in the development or treatment of diseases. Drugs, as important chemical compounds, can interact with microRNAs and change their functions. The experimental identification of microRNA-drug interactions is time-consuming and expensive. Therefore, it is appealing to develop effective computational approaches for predicting microRNA-drug interactions. RESULTS: In this study, a matrix factorization-based method, called the microRNA-drug interaction prediction approach (MDIPA), is proposed for predicting unknown interactions among microRNAs and drugs. Specifically, MDIPA utilizes experimentally validated interactions between drugs and microRNAs, drug similarity and microRNA similarity to predict undiscovered interactions. A path-based microRNA similarity matrix is constructed, while the structural information of drugs is used to establish a drug similarity matrix. To evaluate its performance, our MDIPA is compared with four state-of-the-art prediction methods with an independent dataset and cross-validation. The results of both evaluation methods confirm the superior performance of MDIPA over other methods. Finally, the results of molecular docking in a case study with breast cancer confirm the efficacy of our approach. In conclusion, MDIPA can be effective in predicting potential microRNA-drug interactions. AVAILABILITY AND IMPLEMENTATION: All code and data are freely available from https://github.com/AliJam82/MDIPA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

From Beef to Bees: High-Throughput Kinome Analysis to Understand Host Responses of Livestock Species to Infectious Diseases and Industry-Associated Stress.

Within human health research, the remarkable utility of kinase inhibitors as therapeutics has motivated efforts to understand biology at the level of global cellular kinase activity (the kinome). In contrast, the diminished potential for using kinase inhibitors in food animals has dampened efforts to translate this research approach to livestock species. This, in our opinion, was a lost opportunity for livestock researchers given the unique potential of kinome analysis to offer insight into complex biology. To remedy this situation, our lab developed user-friendly, cost-effective approaches for kinome analysis that can be readily incorporated into most research programs but with a specific priority to enable the technology to livestock researchers. These contributions include the development of custom software programs for the creation of species-specific kinome arrays as well as comprehensive deconvolution and analysis of kinome array data. Presented in this review are examples of the application of kinome analysis to highlight the utility of the technology to further our understanding of two key complex biological events of priority to the livestock industry: host immune responses to infectious diseases and animal stress responses. These advances and examples of application aim to provide both mechanisms and motivation for researchers, particularly livestock researchers, to incorporate kinome analysis into their research programs.

Humanized yeast genetic interaction mapping predicts synthetic lethal interactions of FBXW7 in breast cancer.

BACKGROUND: Synthetic lethal interactions (SLIs) that occur between gene pairs are exploited for cancer therapeutics. Studies in the model eukaryote yeast have identified ~ 550,000 negative genetic interactions that have been extensively studied, leading to characterization of novel pathways and gene functions. This resource can be used to predict SLIs that can be relevant to cancer therapeutics. METHODS: We used patient data to identify genes that are down-regulated in breast cancer. InParanoid orthology mapping was performed to identify yeast orthologs of the down-regulated genes and predict their corresponding SLIs in humans. The predicted network graphs were drawn with Cytoscape. CancerRXgene database was used to predict drug response. RESULTS: Harnessing the vast available knowledge of yeast genetics, we generated a Humanized Yeast Genetic Interaction Network (HYGIN) for 1009 human genes with 10,419 interactions. Through the addition of patient-data from The Cancer Genome Atlas (TCGA), we generated a breast cancer specific subnetwork. Specifically, by comparing 1009 genes in HYGIN to genes that were down-regulated in breast cancer, we identified 15 breast cancer genes with 130 potential SLIs. Interestingly, 32 of the 130 predicted SLIs occurred with FBXW7, a well-known tumor suppressor that functions as a substrate-recognition protein within a SKP/CUL1/F-Box ubiquitin ligase complex for proteasome degradation. Efforts to validate these SLIs using chemical genetic data predicted that patients with loss of FBXW7 may respond to treatment with drugs like Selumitinib or Cabozantinib. CONCLUSIONS: This study provides a patient-data driven interpretation of yeast SLI data. HYGIN represents a novel strategy to uncover therapeutically relevant cancer drug targets and the yeast SLI data offers a major opportunity to mine these interactions.

EPHB6 augments both development and drug sensitivity of triple-negative breast cancer tumours.

Triple-negative breast cancer (TNBC) tumours that lack expression of oestrogen, and progesterone receptors, and do not overexpress the HER2 receptor represent the most aggressive breast cancer subtype, which is characterised by the resistance to therapy in frequently relapsing tumours and a high rate of patient mortality. This is likely due to the resistance of slowly proliferating tumour-initiating cells (TICs), and understanding molecular mechanisms that control TICs behaviour is crucial for the development of effective therapeutic approaches. Here, we present our novel findings, indicating that an intrinsically catalytically inactive member of the Eph group of receptor tyrosine kinases, EPHB6, partially suppresses the epithelial-mesenchymal transition in TNBC cells, while also promoting expansion of TICs. Our work reveals that EPHB6 interacts with the GRB2 adapter protein and that its effect on enhancing cell proliferation is mediated by the activation of the RAS-ERK pathway, which allows it to elevate the expression of the TIC-related transcription factor, OCT4. Consistent with this, suppression of either ERK or OCT4 activities blocks EPHB6-induced pro-proliferative responses. In line with its ability to trigger propagation of TICs, EPHB6 accelerates tumour growth, potentiates tumour initiation and increases TIC populations in xenograft models of TNBC. Remarkably, EPHB6 also suppresses tumour drug resistance to DNA-damaging therapy, probably by forcing TICs into a more proliferative, drug-sensitive state. In agreement, patients with higher EPHB6 expression in their tumours have a better chance for recurrence-free survival. These observations describe an entirely new mechanism that governs TNBC and suggest that it may be beneficial to enhance EPHB6 action concurrent with applying a conventional DNA-damaging treatment, as it would decrease drug resistance and improve tumour elimination.

NovoExD: De novo Peptide Sequencing for ETD/ECD Spectra.

De novo peptide sequencing using tandem mass spectrometry (MS/MS) data has become a major computational method for sequence identification in recent years. With the development of new instruments and technology, novel computational methods have emerged with enhanced performance. However, there are only a few methods focusing on ECD/ETD spectra, which mainly contain variants of c -ions and z-ions. Here, a de novo sequencing method for ECD/ETD spectra, NovoExD, is presented. NovoExD applies a new form of spectrum graph with multiple edge types (called a GMET), considers multiple peptide tags, and integrates amino acid combination (AAC) and fragment ion charge information. Its performance is compared with another successful de novo sequencing method, pNovo+, which has an option for ECD/ETD spectra. Experiments conducted on three different datasets show that the average full length peptide identification accuracy of NovoExD is as high as 88.70 percent, and that NovoExD's average accuracy is more than 20 percent greater on all datasets than that of pNovo+.

De novo peptide sequencing using CID and HCD spectra pairs.

In tandem mass spectrometry (MS/MS), there are several different fragmentation techniques possible, including, collision-induced dissociation (CID) higher energy collisional dissociation (HCD), electron-capture dissociation (ECD), and electron transfer dissociation (ETD). When using pairs of spectra for de novo peptide sequencing, the most popular methods are designed for CID (or HCD) and ECD (or ETD) spectra because of the complementarity between them. Less attention has been paid to the use of CID and HCD spectra pairs. In this study, a new de novo peptide sequencing method is proposed for these spectra pairs. This method includes a CID and HCD spectra merging criterion and a parent mass correction step, along with improvements to our previously proposed algorithm for sequencing merged spectra. Three pairs of spectral datasets were used to investigate and compare the performance of the proposed method with other existing methods designed for single spectrum (HCD or CID) sequencing. Experimental results showed that full-length peptide sequencing accuracy was increased significantly by using spectra pairs in the proposed method, with the highest accuracy reaching 81.31%.

Rapamycin reduces fibroblast proliferation without causing quiescence and induces STAT5A/B-mediated cytokine production.

Rapamycin is a well-known inhibitor of the Target of Rapamycin (TOR) signaling cascade; however, the impact of this drug on global genome function and organization in normal primary cells is poorly understood. To explore this impact, we treated primary human foreskin fibroblasts with rapamycin and observed a decrease in cell proliferation without causing cell death. Upon rapamycin treatment chromosomes 18 and 10 were repositioned to a location similar to that of fibroblasts induced into quiescence by serum reduction. Although similar changes in positioning occurred, comparative transcriptome analyses demonstrated significant divergence in gene expression patterns between rapamycin-treated and quiescence-induced fibroblasts. Rapamycin treatment induced the upregulation of cytokine genes, including those from the Interleukin (IL)-6 signaling network, such as IL-8 and the Leukemia Inhibitory Factor (LIF), while quiescent fibroblasts demonstrated up-regulation of genes involved in the complement and coagulation cascade. In addition, genes significantly up-regulated by rapamycin treatment demonstrated increased promoter occupancy of the transcription factor Signal Transducer and Activator of Transcription 5A/B (STAT5A/B). In summary, we demonstrated that the treatment of fibroblasts with rapamycin decreased proliferation, caused chromosome territory repositioning and induced STAT5A/B-mediated changes in gene expression enriched for cytokines.

Recent Developments in Computational Methods for De Novo Peptide Sequencing from Tandem Mass Spectrometry (MS/MS).

Tandem mass spectrometry (MS/MS) has emerged as a major technology for peptide sequencing. Typically, there are three kinds of methods for the peptide sequencing: database searching, peptide tagging, and de novo sequencing. De novo sequencing has drawn increasing attention because of its independence from existing protein databases and potential for identifying new proteins, proteins resulting from mutations, proteins with unexpected modifications and so on. Recently, with the improvements in the accuracy of MS/MS and development of alternative fragmentation modes of MS/MS, many new de novo sequencing methods have been formulated. This paper reviews these recently developed sequencing methods including those for alternative MS/MS spectra. The paper first introduces background knowledge on peptide sequencing and mass spectrometry, and then reviews de novo peptide sequencing methods for traditional CID spectra. After that, it focuses on the recent development of de novo methods for alternative MS/MS spectra. In addition, methods using multiple spectra from the same peptide are surveyed. Finally, conclusions and some directions of future work are discussed.

Characterization of the host response to pichinde virus infection in the Syrian golden hamster by species-specific kinome analysis.

The Syrian golden hamster has been increasingly used to study viral hemorrhagic fever (VHF) pathogenesis and countermeasure efficacy. As VHFs are a global health concern, well-characterized animal models are essential for both the development of therapeutics and vaccines as well as for increasing our understanding of the molecular events that underlie viral pathogenesis. However, the paucity of reagents or platforms that are available for studying hamsters at a molecular level limits the ability to extract biological information from this important animal model. As such, there is a need to develop platforms/technologies for characterizing host responses of hamsters at a molecular level. To this end, we developed hamster-specific kinome peptide arrays to characterize the molecular host response of the Syrian golden hamster. After validating the functionality of the arrays using immune agonists of defined signaling mechanisms (lipopolysaccharide (LPS) and tumor necrosis factor (TNF)-α), we characterized the host response in a hamster model of VHF based on Pichinde virus (PICV(1)) infection by performing temporal kinome analysis of lung tissue. Our analysis revealed key roles for vascular endothelial growth factor (VEGF), interleukin (IL) responses, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling, and Toll-like receptor (TLR) signaling in the response to PICV infection. These findings were validated through phosphorylation-specific Western blot analysis. Overall, we have demonstrated that hamster-specific kinome arrays are a robust tool for characterizing the species-specific molecular host response in a VHF model. Further, our results provide key insights into the hamster host response to PICV infection and will inform future studies with high-consequence VHF pathogens.

NovoHCD: de novo peptide sequencing from HCD spectra.

In recent years, de novo peptide sequencing from mass spectrometry data has developed as one of the major peptide identification methods with the emergence of new instruments and advanced computational methods. However, there are still limitations to this method; for example, the typically used spectrum graph model cannot represent all the information and relationships inherent in tandem mass spectra (MS/MS spectra). Here, we present a new method named NovoHCD which applies a spectrum graph model with multiple types of edges (called a multi-edge graph), and integrates into it amino acid combination (AAC) information and peptide tags. In addition, information on immonium ions observed particularly in higher-energy collisional dissociation (HCD) spectra is incorporated. Comparisons between NovoHCD and another successful de novo peptide sequencing method for HCD spectra, pNovo, were performed. Experiments were conducted on five HCD spectral datasets. Results show that NovoHCD outperforms pNovo in terms of full length peptide identification accuracy; specifically, the accuracy increases 13%-21% over the five datasets.

TFPI1 mediates resistance to doxorubicin in breast cancer cells by inducing a hypoxic-like response.

Thrombin and hypoxia are important players in breast cancer progression. Breast cancers often develop drug resistance, but mechanisms linking thrombin and hypoxia to drug resistance remain unresolved. Our studies using Doxorubicin (DOX) resistant MCF7 breast cancer cells reveals a mechanism linking DOX exposure with hypoxic induction of DOX resistance. Global expression changes between parental and DOX resistant MCF7 cells were examined. Westerns, Northerns and immunocytochemistry were used to validate drug resistance and differentially expressed genes. A cluster of genes involved in the anticoagulation pathway, with Tissue Factor Pathway Inhibitor 1 (TFPI1) the top hit, was identified. Plasmids overexpressing TFPI1 were utilized, and 1% O2 was used to test the effects of hypoxia on drug resistance. Lastly, microarray datasets from patients with drug resistant breast tumors were interrogated for TFPI1 expression levels. TFPI1 protein levels were found elevated in 3 additional DOX resistant cells lines, from humans and rats, indicating evolutionarily conservation of the effect. Elevated TFPI1 in DOX resistant cells was active, as thrombin protein levels were coincidentally low. We observed elevated HIF1α protein in DOX resistant cells, and in cells with forced expression of TFPI1, suggesting TFPI1 induces HIF1α. TFPI1 also induced c-MYC, c-SRC, and HDAC2 protein, as well as DOX resistance in parental cells. Growth of cells in 1% O2 induced elevated HIF1α, BCRP and MDR-1 protein, and these cells were resistant to DOX. Our in vitro results were consistent with in vivo patient datasets, as tumors harboring increased BCRP and MDR-1 expression also had increased TFPI1 expression. Our observations are clinically relevant indicating that DOX treatment induces an anticoagulation cascade, leading to inhibition of thrombin and the expression of HIF1α. This in turn activates a pathway leading to drug resistance.

Kinotypes: stable species- and individual-specific profiles of cellular kinase activity.

BACKGROUND: Recently, questions have been raised regarding the ability of animal models to recapitulate human disease at the molecular level. It has also been demonstrated that cellular kinases, individually or as a collective unit (the kinome), play critical roles in regulating complex biology. Despite the intimate relationship between kinases and health, little is known about the variability, consistency and stability of kinome profiles across species and individuals. RESULTS: As a preliminary investigation of the existence of species- and individual-specific kinotypes (kinome signatures), peptide arrays were employed for the analysis of peripheral blood mononuclear cells collected weekly from human and porcine subjects (n = 6) over a one month period. The data revealed strong evidence for species-specific signalling profiles. Both humans and pigs also exhibited evidence for individual-specific kinome profiles that were independent of natural changes in blood cell populations. CONCLUSIONS: Species-specific kinotypes could have applications in disease research by facilitating the selection of appropriate animal models or by revealing a baseline kinomic signature to which treatment-induced profiles could be compared. Similarly, individual-specific kinotypes could have implications in personalized medicine, where the identification of molecular patterns or signatures within the kinome may depend on both the levels of kinome diversity and temporal stability across individuals.

Systems kinomics demonstrates Congo Basin monkeypox virus infection selectively modulates host cell signaling responses as compared to West African monkeypox virus.

Monkeypox virus (MPXV) is comprised of two clades: Congo Basin MPXV, with an associated case fatality rate of 10%, and Western African MPXV, which is associated with less severe infection and minimal lethality. We thus postulated that Congo Basin and West African MPXV would differentially modulate host cell responses and, as many host responses are regulated through phosphorylation independent of transcription or translation, we employed systems kinomics with peptide arrays to investigate these functional host responses. Using this approach we have demonstrated that Congo Basin MPXV infection selectively down-regulates host responses as compared with West African MPXV, including growth factor- and apoptosis-related responses. These results were confirmed using fluorescence-activated cell sorting analysis demonstrating that West African MPXV infection resulted in a significant increase in apoptosis in human monocytes as compared with Congo Basin MPXV. Further, differentially phosphorylated kinases were identified through comparison of our MPXV data sets and validated as potential targets for pharmacological inhibition of Congo Basin MPXV infection, including increased Akt S473 phosphorylation and decreased p53 S15 phosphorylation. Inhibition of Akt S473 phosphorylation resulted in a significant decrease in Congo Basin MPXV virus yield (261-fold) but did not affect West African MPXV. In addition, treatment with staurosporine, an apoptosis activator resulted in a 49-fold greater decrease in Congo Basin MPXV yields as compared with West African MPXV. Thus, using a systems kinomics approach, our investigation demonstrates that West African and Congo Basin MPXV differentially modulate host cell responses and has identified potential host targets of therapeutic interest.

Computational prediction of eukaryotic phosphorylation sites.

MOTIVATION: Kinase-mediated phosphorylation is the central mechanism of post-translational modification to regulate cellular responses and phenotypes. Signaling defects associated with protein phosphorylation are linked to many diseases, particularly cancer. Characterizing protein kinases and their substrates enhances our ability to understand and treat such diseases and broadens our knowledge of signaling networks in general. While most or all protein kinases have been identified in well-studied eukaryotes, the sites that they phosphorylate have been only partially elucidated. Experimental methods for identifying phosphorylation sites are resource intensive, so the ability to computationally predict potential sites has considerable value. RESULTS: Many computational techniques for phosphorylation site prediction have been proposed, most of which are available on the web. These techniques differ in several ways, including the machine learning technique used; the amount of sequence information used; whether or not structural information is used in addition to sequence information; whether predictions are made for specific kinases or for kinases in general; and sources of training and testing data. This review summarizes, categorizes and compares the available methods for phosphorylation site prediction, and provides an overview of the challenges that are faced when designing predictors and how they have been addressed. It should therefore be useful both for those wishing to choose a phosphorylation site predictor for their particular biological application, and for those attempting to improve upon established techniques in the future. CONTACT: brett.trost@usask.ca.

Codon number shapes peptide redundancy in the universal proteome composition.

The proteomes catalogued in the UniRef100 database were collected into a single proteome set and examined for actual versus theoretical pentapeptide occurrences. We found a highly diversified degree of pentapeptide redundancy. Numerically, 953 pentamers are expressed only once in the protein world, whereas 103 pentamers occur more than 50,000 times. Moreover, it seems that 417 potentially possible pentapeptides are not present in the protein world. On the whole, tracing the redundancy profile of the protein world as a function of pentapeptide occurrences reveals a quasi-Gaussian curve, with tails representing scarcely and repeatedly occurring 5-mers. Analysis of physico-chemical-biological parameters shows that codon number is the main factor influencing and favoring specific pentapeptide frequencies in the universal proteome composition. That is, when compared to the set of never-expressed 5-mers, the pentapeptides frequently represented in the universal proteome are endowed with a higher number of multi-codonic amino acids. In contrast, the bulkiness degree and the hydrophobicity level play a smaller role. Unexpectedly, the heat of formation of pentapeptide appears to have the least influence.

Rare peptide segments are found significantly more often in proto-oncoproteins than control proteins: implications for immunology and oncology.

There is some evidence to suggest that peptide segments that are found rarely or never in the host proteome play a role in the immune response to disease-related proteins, both those derived from microbes and those derived from the host itself. We conjecture that this pattern may extend to human proto-oncoproteins. Our hypothesis in this study is that the frequency of rare peptide segments in sets of human proto-oncoproteins is significantly higher than in sets of control proteins, and we show that this is the case. Possible immunological implications of this observation are discussed.