The Cytoplasmic Peptidase DPP9 Is Rate-limiting for Degradation of Proline-containing Peptides

Protein degradation is an essential process that continuously takes place in all living cells. Regulated degradation of most cellular proteins is initiated by proteasomes, which produce peptides of varying length. These peptides are rapidly cleaved to single amino acids by cytoplasmic peptidases. Proline-containing peptides pose a specific problem due to structural constrains imposed by the pyrrolidine ring that prevents most peptidases from cleavage. Here we show that DPP9, a poorly characterized cytoplasmic prolyl-peptidase, is rate-limiting for destruction of proline-containing substrates both in cell extracts and in intact cells. We identified the first natural substrate for DPP9, the RU1_34–42 antigenic peptide (VPYGSFKHV). RU1_34–42 is degraded in vitro by DPP9, and down-regulation of DPP9 in intact cells results in increased presentation of this antigen. Together our findings demonstrate an important role for DPP9 in peptide turnover and antigen presentation.Protein turn-over is an essential process that continuously occurs in all living cells. The ubiquitin-proteasome system is responsible for initiating the regulated degradation of most cellular proteins (1). Proteasome-degradation products are not single amino acids but rather peptides varying in length between 3 and 22 amino acids (2, 3). Cytosolic amino- and endopeptidases rapidly cleave these peptides (4) to allow recycling of amino acids and to prevent accumulation of short peptides, which may be harmful to the cell. In addition, these peptidases also play an important role in the trimming of proteasomal products for antigen presentation on MHC⁴ class I (5–8).Peptides containing proline residues pose a problem for most peptidases due to the pyrrolidine ring of proline that gives it an exceptional conformational rigidity. Only few peptidases are known to cleave after prolines, including the cytoplasmic peptidases prolyl oligopeptidase (POP) and cytoplasmic members of the S9B/DPPIV family (DPP8 and DPP9). POP is a cytosolic endopeptidase of the S9A family, which is broadly distributed with high concentrations in the brain. It has been implicated in the maturation and degradation of peptide hormones and neuropeptides (9, 10).S9B/DPPIV peptidases are a family of exopeptidases that cleave off N-terminal dipeptides from proteins/polypeptides having a proline residue at the second position (Xaa-Pro). The best-characterized member of this family is DPPIV, a membrane protein with a catalytic domain facing the extracellular space. DPPIV knock-out mice show enhanced insulin secretion and improved glucose tolerance (11, 12). This is due to cleavage and, thus, inactivation of the incretin hormones glucagon-like peptide and glucose-dependent insulinotropic polypeptide by DPPIV (13–15). Therefore, DPPIV is used as a drug target for the treatment of diabetes type 2.In contrast, DPP8 and DPP9 are soluble cytoplasmic peptidases of unknown function. They share 60% amino acid identity and are ubiquitously expressed in vertebrate tissues (16–20). Because DPP8 and DPP9 knock-out mice are not available, most studies on these enzymes were done with inhibitors either against the DPPIV family or specifically against DPP8 and -9. Currently two specific DPP8/9 inhibitors are described (21, 22), of which one showed severe effects in animal models (21).Here we show that DPP9 is a rate-limiting enzyme for cytosolic post-proline aminodipeptidase activity. Our work associates an in vivo function with DPP9 in peptide degradation and also suggests that changes in DPP9 expression levels or activity contribute to changes in the repertoire of cytosolic peptides, including those presented by MHC class I.     

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]     

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]     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[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]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[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]     

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]     

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.