Skip to main content

The oligodeoxynucleotide sequences corresponding to never-expressed peptide motifs are mainly located in the non-coding strand

Abstract

Background

We study the usage of specific peptide platforms in protein composition. Using the pentapeptide as a unit of length, we find that in the universal proteome many pentapeptides are heavily repeated (even thousands of times), whereas some are quite rare, and a small number do not appear at all. To understand the physico-chemical-biological basis underlying peptide usage at the proteomic level, in this study we analyse the energetic costs for the synthesis of rare and never-expressed versus frequent pentapeptides. In addition, we explore residue bulkiness, hydrophobicity, and codon number as factors able to modulate specific peptide frequencies. Then, the possible influence of amino acid composition is investigated in zero- and high-frequency pentapeptide sets by analysing the frequencies of the corresponding inverse-sequence pentapeptides. As a final step, we analyse the pentadecamer oligodeoxynucleotide sequences corresponding to the never-expressed pentapeptides.

Results

We find that only DNA context-dependent constraints (such as oligodeoxynucleotide sequence location in the minus strand, introns, pseudogenes, frameshifts, etc.) provide a coherent mechanistic platform to explain the occurrence of never-expressed versus frequent pentapeptides in the protein world.

Conclusions

This study is of importance in cell biology. Indeed, the rarity (or lack of expression) of specific 5-mer peptide modules implies the rarity (or lack of expression) of the corresponding n-mer peptide sequences (with n > 5), so possibly modulating protein compositional trends. Moreover the data might further our understanding of the role exerted by rare pentapeptide modules as critical biological effectors in protein-protein interactions.

Background

Proteins comprise subsets of all plausible amino acid sequences, i.e. peptide motifs that occur in different quantitative percentages and with different qualitative significance at the proteomic level. To understand the correspondence between structure and function, we must understand the rules dictating the modular arrangement of proteins. We chose the pentapeptide as a basic structural/functional unit to analyse the compositional distribution of peptide sequences. Indeed, pentapeptides appear to be minimal biological units exerting a central role in fundamental cellular processes such as inhibition/stimulation of cell growth, hormone activity, regulation of transcript expression, enzyme activity, and immune recognition [1]. Following a robust set of experimental protein analyses [2–9], we determined that, as a rule, amino acid stretches with low/no proteomic redundancy alternate with portions of high proteomic redundancy along protein primary structures [2], independently of the protein length [3, 4], whether the protein is derived from microbial or mammalian organisms [3–9], and the proteome under analysis [5–9]. Preliminarily to any evolutionary/functional/physio-pathological considerations, the data prompt a fundamental question: what makes one pentapeptide occur more frequently than another in the protein world? In this paper, we undertake a large-scale analysis of the physico-(bio)chemical factors that theoretically might account for the modular peptide composition of proteins, and examine a total of 20991 pentapeptides, divided into eleven sets characterized by frequencies ranging from zero to 2500.

Methods

The complete UniRef100, UniRef90 and UniRef50 databases (http://www.uniprot.org/downloads) were downloaded as single proteomes and analysed for internal peptide redundancy using 5-mers sequentially overlapping by four residues. The scans were performed using standard UNIX/LINUX commands and custom programs written in Perl [10].

The proteins were manipulated and analysed as follows. All the protein sequences were decomposed in silico to a set of 5-mers (including all duplicates). Any 5-mers containing ambiguous amino acids (i.e., denoted by the letters B, X, or Z, which respectively represent ambiguity between N and D, ambiguity between Q and E, and an unknown amino acid) or non-standard amino acid codes (i.e., -, U, *, O, denoting gaps, selenocysteine residues, stop codons, etc.) were eliminated. Since there are only 3200000 possible 5-mers, a simple linear scan was used to determine the counts of occurrences and 5-mers that do not occur. That is, for each pentamer, the UniRef100 (or UniRef90 or UniRef50) proteome was searched for instances of that pentamer. Any such occurrence was termed a match. The number of matches defines the proteomic frequency of each pentapeptide.

Eleven peptide sets with zero, low, medium and high frequencies (i.e., from zero to 2500 matches) were selected from UniRef100 (hereafter called the "universal proteome") for physico-(bio)chemical analyses. Specifically, the frequencies defining the eleven sets were: 0, 1, 4, 5, 50, 100, 341, 500, 1000, 1368 and 2500. The pentapeptide sets were screened by starting with the UniRef100 database and then using the Perfect Peptide Match program at the Protein Information Resource (PIR) website (http://pir.georgetown.edu/pirwww)[11] to eliminate repeated sequences and fragments. The protein entries containing the 5-mer under analysis were further filtered using the UniProtKB resources (http://www.uniprot.org) to eliminate obsolete entries.

Analysis of the energetics was carried out for each pentaptide using Spartan'06 software (from Wavefunction Inc, Irvine, CA) and applying the semi-empirical method. The peptide bulkiness degree was measured using the ProtScale program available at http://www.expasy.ch/tools[12]. The hydrophobicity level was determined using the scale described by Takano and Yutani [13]. The codon number per pentapeptide was calculated by summing the number of codons of each amino acid forming the 5-mer. One-way analysis of variance (ANOVA, F-test) was used to derive a p-value indicating whether the means of the measurements for the different sets were all equal.

To analyse DNA constraints, we analysed the oligodeoxynucleotide coding sequences corresponding to the pentameric amino acid sequences. The Sequence Manipulation Suite Reverse Translate program (http://www.bioinformatics.org/sms2/)[14] was used to generate a DNA sequence representing the most likely, optimized coding sequence. Additionally, Reverse Translate a Protein (http://www.vivo.colostate.edu/molkit/rtranslate/index.html), a program that uses the standard genetic code and does not consider differences in codon usage, was used in order to obtain all the possible degenerate oligodeoxynucleotide coding frames for each pentapeptide under analysis.

The pentadecameric oligodeoxynucleotide sequences so obtained were the subject of nucleotide-nucleotide BLAST (blastn) analysis at NCBI (http://blast.ncbi.nlm.nih.gov) to find and localize regions of 100% similarity (i.e. with no gaps allowed) in the entire nucleotide collection (nr/nt) comprehending genomic and transcript sequences [15].

Results

Pentapeptide redundancy and ΔG°

The biosynthesis of the peptide bond from amino acids involves an increase in free energy and must therefore depend on energy yielding reactions. We reasoned that, if a substantial fraction of energy is needed to convert starting amino acids into peptides, then the pentapeptide composition of proteins expressed in the proteomes should be biased toward less energetically costly pentapeptides. Theoretically, the extent to which pentapeptide composition is biased to reduce metabolic costs should positively correlate with the pentapeptide redundancy at the proteomic level.

Consequently, we analysed rare versus frequent pentapeptides for the standard enthalpy (or standard heat of formation) associated with the synthesis of the peptide bond. This quantity is highly variable, with the heat generated or absorbed during the formation of a peptide bond depending on the amino acids involved. As an example, Figure 1A reports the frequency distribution of the 400 dipeptides present in the protein world, and, in parallel, the heat of formation in kJ/mol as determined using Spartan '06 software (Figure 1B). It can be seen that the standard heat of formation of the semi-empirically optimised dipeptide structures varies widely from the highly exothermic value of DE dipeptide formation (-944.34 kJ/mol) to the endothermic CP dipeptide formation (1062.47 kJ/mol) (Figure 1B). Moreover, Figure 1 shows that a negative correlation exists between dipeptide redundancy (Figure 1A) and ΔG° level (Figure 1B). As a synthetic datum, mean ΔG° values equal to -219.06 ± 237.83 kJ/mol and 56.34 ± 249.51 kJ/mol characterize the 50 most frequent dipeptides and the 50 less frequent ones, respectively.

Figure 1
figure 1

Correlation between frequency distribution and standard heat of peptide formation of the 400 dipeptides present in the protein world. Panel A: The frequency distribution of the 400 dipeptides present in the protein world widely varies from a maximum of 21,927,296 times (LL dipeptide) to a minimum of 407,573 (WC dipeptide). Panel B: The heat of formation in kJ/mol as determined by the Spartan'06 software for the 400 dipeptides varies widely from the highly exothermic value of DE dipeptide formation (-944.34 kJ/mol) to the endothermic CP dipeptide formation (1062.47 kJ/mol). Mean ΔG° values equal to -219.06 ± 237.83 kJ/mol and 56.34 ± 249.51 kJ/mol characterize the 50 most frequent dipeptides and the 50 less frequent ones, respectively.

Therefore, we reasoned that the same ΔG° variability would apply even more strongly to longer peptide units. Based on this rationale, we calculated the heats of formation for pentapeptide sets with different frequencies in the universal proteome (i.e., from zero to 2500 occurrences). As a universal proteome database, we used UniRef100, which representsone of the most comprehensive non-redundant protein sequence datasets available ([16–18], see also http://www.ebi.ac.uk/uniref/). To control for existing bias and redundancies in the UniRef100 database, the protein entries containing the 5-mers under analysis were filtered for repeated sequences, fragments, and obsolete entries.

Figure 2 reports the distribution of pentapeptide frequencies in the universal proteome as the log of the occurrence count versus the number of 5-mers with that count. The same trend in the quantitative pentapeptide composition of the protein world was observed using UniRef90 and UniRef50 protein datasets (not shown). Then, we selected pentapeptide sets for physico-chemical analyses along the distribution curve of pentapeptide frequencies shown in Figure 2. The frequencies of the different 5-mer sets elected for analysis correspond to 1, 4, 5, 50, 100, 341, 500, 1000, 1368 and 2500 occurrences (as indicated by the lettered arrows b to k, plus a 5-mer set having zero occurrences, namely a). That is, we selected: peptides occurring just once for the obvious reason that such peptides are expected to be "interesting"; the occurrence count of 50 was chosen because the maximum is reached at this point; 341 was chosen because it is a median value (i.e. half the peptides have occurrence counts less than this value, and half the peptides have occurrence counts more than this value); high occurrence counts (e.g. 2500) were chosen to represent the "tail" of the distribution. However, extremely high counts (e.g. 5000 or more) were not chosen because the number of pentapeptides with these frequencies tended to be too small to give results in which we were confident. Finally, other occurrence counts were chosen so as to broaden this sampling.

Figure 2
figure 2

Location of the 5-mer sets selected for physico-chemical analyses along the distribution curve of pentapeptide frequencies in the universal proteome. UniRef100, the most comprehensive protein dataset available [[16–18], see also http://www.ebi.ac.uk/uniref/], was used. The arrows, lettered from b to k, indicate the frequencies of the different 5-mer sets corresponding, in the order, to 1, 4, 5, 50, 100, 341, 500, 1000, 1368 and 2500 occurrences and selected for physico-chemical analyses. A further set a, corresponding to the set of never-occurring pentapeptides, was also chosen.

Afterwards, we calculated the relationship between metabolic costs of pentapeptide biosynthesis (as estimated from heat of peptide bond formation data) and pentapeptide redundancy (as estimated by the number of occurrences). The histograms reported in Figure 3, Panels A to E, refer to the energetic profiles of pentapeptide sets with the following different frequencies in the universal proteome: A) never expressed, B) expressed only once, C) occurring 100 times, D) occurring 341 times, and E) occurring 2500 times. It can be seen that the range in heat of formation values varies considerably across the five sets of pentapeptides. For instance, many among the high frequency pentapeptides have extremely high or low heat-of-formation values (panels C, D and E), while the absent or rare pentapeptides fall into an energetically narrower window (panels A and B). Because of this large variance, considering the central tendency in each panel does not seem to allow one to distinguish among the sets.

Figure 3
figure 3

Energetic cost of pentapeptides with different frequencies in the universal proteome. Panels A to E indicate pentapeptide sets that, in the universal proteome: A) are absent, B) are expressed only once, C) occur 100 times, D) occur 341 times, and E) occur 2500 times.

This result is even more clear in the boxplot diagram reported in Figure 4, where the analysis of the distribution of ΔG° scores is extended to eleven pentapeptide sets occurring with different frequencies in the universal proteome (see Figure 2). It is evident that the never-expressed pentapeptides are confined to restricted energy levels, i.e. have smaller variance, while, on the contrary, many of the pentapeptides occurring repeatedly in the universal proteome have higher energetic costs. Moreover, specifically and importantly, the boxplot diagram shows that outliers are usually associated with high frequency pentapeptides rather than rare ones. Figure 4 clearly shows that the heat of formation has no stringent influence on pentapeptide frequency.

Figure 4
figure 4

Statistical characterization of the energetic cost of pentapeptide sets with different frequencies in the universal proteome. The boxplots show the distribution of ΔG° values for each set of pentapeptides. The line within each box represents the median value. The top and bottom of each box represent the 75th and 25th percentile, respectively. The whiskers show the range of values that are not considered to be outliers. Outliers are plotted individually as plus signs. The p-value was 0.008, indicating that the means of the different sets are different, though clearly the magnitude of the differences is small.

The relationship between pentapeptide redundancy and hydrophobicity, bulkiness, and codon number

We analysed the relationship between pentapeptide frequencies and the following physico-(bio)chemical parameters: side-chain bulkiness, hydrophobicity and amino acid codon number. The results are reported in Figure 5: it can be seen that the pentapeptide redundancy appears to be shaped by, in order of importance, the amino acid codon number (panel C), residue hydrophobicity (panel A), and residue bulkiness (panel B). However, in all instances many values are outliers, indicating a non-stringent relationship between the physico-chemical factors analysed and the distribution of pentapeptide redundancy in the universal proteome.

Figure 5
figure 5

Statistical characterization of hydrophobicity (A), bulkiness (B), and amino acid codon number (C) for pentapeptide sets with different frequencies in the universal proteome. The boxplots show the distribution of the values of each physico-biochemical factor for each set of pentapeptides. The line within each box represents the median value. The top and bottom of each box represent the 75th and 25th percentile, respectively. The whiskers show the range of values that are not considered to be outliers. Outliers are plotted individually as plus signs. The p-values among the different classes of 5-mers for hydrophobicity, bulkiness, and amino acid codon number were all less than 0.001, indicating in each case that the means of the different sets are different.

Pentapeptide redundancy and amino acid composition

Figures 4 and 5 indicate almost no relationship between pentapeptide frequencies and physico-chemical factors such as hydrophobicity and bulkiness. On the other hand, the analyses reported in Figure 3 suggest that rare pentapeptides are formed primarily by Trp, Tyr, and Met, i.e. by essential low-concentration amino acids endowed with high values of hydrophobicity and residue bulkiness. This raises the question: might amino acid frequencies affect pentapeptide frequency?

To analyse the relationship between pentapeptide frequency and amino acid composition, we used the pentapeptide set with zero occurrences and investigated the frequency of the corresponding inverse amino acid sequences. We reasoned that if the factor dictating the rarity/frequency of a certain pentapeptide was specific amino acid composition, then inverting the order of those amino acids but keeping constant the amino acid composition would have little or no effect on pentapeptide occurrence. Panels A and B of Figure 6 show that the inverse sequences of the never-expressed pentapeptides occur in the universal proteome as many as 50 times. Hence, amino acid composition does not represent the factor precluding the expression of the zero-frequency pentapeptide set. Similar results were obtained using the set of pentapeptides with 2500 occurrences in the universal proteome: Figure 6 panel D shows that the inverse amino acid sequences occur in the universal proteome with a wide variety of frequencies. As a further control the frequency of pentapeptides uniquely formed by the rare W, Y and M amino acids was determined. We found that the highly structured WWWWW, YYYYY and MMMMM pentapeptides occur 112, 972 and 1568 times, respectively, in the universal proteome. I.e., pentapeptides formed by rare, mono-codonic, highly hydrophobic, and bulky amino acid residues can even fall in the category of the "highly repeated" pentamers.

Figure 6
figure 6

Effect of amino acid composition on pentapeptide frequency in the universal proteome. Frequency of: A) pentapeptides never-expressed in the universal proteome, and B) their inverse sequences. Frequency of: C) pentapeptides with 2500 occurrences in the proteome, and D) their inverse sequences.

Taken together, these data indicate that amino acid composition appears to modulate at some extent, but does not dictate, the pentapeptide composition of the universal proteome.

Analysing the never-expressed pentapeptides at the DNA level

After obtaining the results above, we postulated that the lack of occurrence of the pentapeptides never found in the universal proteome could be ascribed to a lack of the corresponding pentadecameric oligodeoxynucleotides in the DNA coding sequence. Therefore, a search was conducted for occurrences of the oligodeoxynucleotide sequences coding for the pentapeptides never expressed in the universal proteome using the standard nucleotide-nucleotide BLAST (blastn) program as described under Methods.

The two examples reported in Table 1 synthetically illustrate that while all of the pentadecameric oligodeoxynucleotide sequences corresponding to the zero-frequency pentapeptides are present in a number of different organisms, they are mainly located in the DNA minus strand, introns, frameshifts, or pseudogenes, i.e. in untranslatable DNA positions/structures. The data from Table 1 are further confirmed by the data given in Additional file 1, Table S1, where analysis of the most likely and degenerate oligodeoxynucleotide coding frames for each pentapeptide sequence is reported.

Table 1 The oligodeoxynucleotide sequences corresponding to never-expressed peptide motifs are mainly located in the non-coding strand

From this we conclude that DNA context-dependent constraints (e.g., oligodeoxynucleotide sequence location in the minus strand, introns, splicing-dependent frameshifts, etc.) are the main factors limiting/preventing the expression of the corresponding amino acid sequences in the universal proteome.

Discussion

The factors acting on the amino acid composition of proteins have been thoroughly investigated with particular attention to the habitat of the organisms (e.g., growth temperature and salinity) [19–22], sub-cellular localization (e.g., cytosolic, membrane or nuclear) [23], physical properties (e.g., mass and charge) [24], translational constraints [25], and the metabolic costs of amino acid biosynthesis [26]. In contrast, less attention has been dedicated to the structural and functional constraints acting on the peptide composition of proteins. Clearly, the empirical distribution of pentapeptide frequencies has, one way or another, an impact upon protein expression as well as on function/structure, and it is important to understand and define the physico-chemical-biological factors that correlate with pentapeptide frequencies in the protein world.

We already reported preliminary data showing that certain short sequences of amino acids (i.e. pentapeptides) are very common, whereas some are quite rare, and a small number do not appear at all in the collection of all known proteins [27]. Here we report the results of a comprehensive study of the influence of physico-(bio)chemical parameters (energetic cost, bulkiness, hydrophobicity and amino acid codon number), amino acid composition, and DNA constraints on pentapeptide expression in the protein world.

First, we observe a definite (although not determining) role of, in descending order of importance, amino acid codon number, hydrophobicity and bulkiness in modulating pentapeptide frequency in the universal proteome. On the other hand, we find that ΔG° has little influence in defining the pentapeptide composition of the universal proteome. This result is relevant and deserves to be emphasized. We explored in detail whether variations in the peptide bond energetical cost might explain the extent of the pentapeptide compositional bias in the universal proteome based on the following rationale. The data reported for protein amino acid composition indicate increases in the abundance of less energetically costly amino acids in highly expressed proteins [26]. Accordingly and further supported by the correlation existing between dipeptide redundancy (Figure 1A) and ΔG° level (Figure 1B), we expected that energetically costly pentapeptides would be rare, whereas more frequent pentapeptides would have a low energetic cost. In conflict with this theoretical expectation, the experimental data obtained in this study and reported in Figures 3 and 4 clearly demonstrate that there is no such correlation at the pentapeptide level. Surprisingly we found that high energies of formation are associated with moderately or highly frequent pentapeptides.

A second unexpected finding is that amino acid composition is a marginal factor in determining pentapeptide rarity: although enriched in hydropathic, rare amino acids such as Trp, Tyr, and Met, the inverse sequences of never-expressed pentapeptides are indeed expressed in the universal proteome.

Third, and as a logical consequence of the previous two points, we show that the constraints acting on pentapeptide expression mainly lie at the nucleotide sequence level. Once we excluded possible limitations due to Trp, Met, and Tyr rarity [28] (see Figure 6), we had to suppose that other constraints are active in defining the proteomic pentapeptide frequencies. Effectively, as demonstrated in Table 1 (see also Additional file 1, Table S1), we found that never-expressed pentapeptides correspond to untranslatable, frameshifted or mistranslated oligodeoxynucleotide sequences. In other words, allocation of the coding oligodeoxynucleotide in pseudogenes/minus strand/untranslated regions/introns as well the shift of the reading frame are the main factors determining the distribution of pentapeptide frequencies throughout the protein world.

Conclusions

The results above are of importance both in the biochemical and functional cellular context. Indeed, as already described [1, 3–10, 29], it seems that rare pentapeptides are basic to control functions [30], whereas possibly frequent modules are preferentially involved in structure definition. In this regard, it is worth noting that multicodonic Leu, Ser, Pro, Ala, and Gly residues are the most common ones in high-frequency, low-complexity peptides whose function, in many cases, is the spacing of structural/functional domains [31]. Conversely, the mono/di-codonic amino acids Asn, Cys, Tyr, Met, Phe, and Trp are relatively rare in highly-frequent, low-complexity peptides and characterize functionally critical proteins such as proto-oncogenes [29]. In this case, specific usage of mono/di-codonic amino acids would allow the control of the proto-oncogene product at the transcriptional level. Moreover, during the last decade one of us proposed and demonstrated the association between rare pentapeptides and immunogenic potential [32–39]. Hence, understanding the mechanisms by which peptide platforms are used in the protein world not only is of biochemical interest but also proves of practical importance for biotechnology, e.g. vaccines, expression vectors and peptide therapy approaches [40] with the relevant advantage of effectiveness [41] without adverse side-effects [42–45].

References

  1. Lucchese G, Stufano A, Trost B, Kusalik A, Kanduc D: Peptidology: short amino acid modules in cell biology and immunology. Amino Acids 2007, 33: 703–707. 10.1007/s00726-006-0458-z

    Article  CAS  PubMed  Google Scholar 

  2. Kanduc D, Capone G, Delfino VP, Losa G: The fractal dimension of protein information. Adv Stud Biol 2010, 2: 53–62.

    Google Scholar 

  3. Kanduc D, Lucchese A, Mittelman A: Individuation of monoclonal anti-HPV16 E7 antibody linear peptide epitope by computational biology. Peptides 2001, 22: 1981–1985. 10.1016/S0196-9781(01)00539-3

    Article  CAS  PubMed  Google Scholar 

  4. Mittelman A, Tiwari R, Lucchese G, Willers J, Dummer R, Kanduc D: Identification of monoclonal anti-HMW-MAA antibody linear peptide epitope by proteomic database mining. J Invest Dermatol 2004, 123: 670–675. 10.1111/j.0022-202X.2004.23417.x

    Article  CAS  PubMed  Google Scholar 

  5. Mittelman A, Lucchese A, Sinha AA, Kanduc D: Monoclonal and polyclonal humoral immune response to EC HER-2/NEU peptides with low similarity to the host's proteome. Int J Cancer 2002, 98: 741–747. 10.1002/ijc.10259

    Article  CAS  PubMed  Google Scholar 

  6. Lucchese A, Mittelman A, Lin MS, Kanduc D, Sinha AA: Epitope definition by proteomic similarity analysis: identification of the linear determinant of the anti-Dsg3 MAb 5H10. J Transl Med 2004, 2: 43. 10.1186/1479-5876-2-43

    Article  PubMed  PubMed Central  Google Scholar 

  7. Lucchese A, Willers J, Mittelman A, Kanduc D, Dummer R: Proteomic scan for tyrosinase peptide antigenic pattern in vitiligo and melanoma: role of sequence similarity and HLA-DR1 affinity. J Immunol 2005, 175: 7009–7020.

    Article  CAS  PubMed  Google Scholar 

  8. Willers J, Lucchese A, Mittelman A, Dummer R, Kanduc D: Definition of anti-tyrosinase MAb T311 linear determinant by proteome-based similarity analysis. Exp Dermatol 2005, 14: 543–550. 10.1111/j.0906-6705.2005.00327.x

    Article  CAS  PubMed  Google Scholar 

  9. Stufano A, Kanduc D: Proteome-based epitopic peptide scanning along PSA. Exp Mol Pathol 2009, 86: 36–40. 10.1016/j.yexmp.2008.11.009

    Article  CAS  PubMed  Google Scholar 

  10. Gusfield D: Algorithms on strings, trees, and sequences: computer science and computational biology. Cambridge: Cambridge University Press; 1997.

    Book  Google Scholar 

  11. Wu CH, Huang H, Arminski L, Castro-Alvear J, Chen Y, Hu ZZ, Ledley RS, Lewis KC, Mewes HW, Orcutt BC, Suzek BE, Tsugita A, Vinayaka CR, Yeh LS, Zhang J, Barker WC: The Protein Information Resource: an integrated public resource of functional annotation of proteins. Nucleic Acids Res 2002, 30: 35–37. 10.1093/nar/30.1.35

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zimmerman JM, Eliezer N, Simha R: The characterization of amino acid sequences in proteins by statistical methods. J Theor Biol 1968, 21: 170–201. 10.1016/0022-5193(68)90069-6

    Article  CAS  PubMed  Google Scholar 

  13. Takano K, Yutani K: A new scale for side-chain contribution to protein stability based on the empirical stability analysis of mutant proteins. Protein Eng 2001, 14: 525–528. 10.1093/protein/14.8.525

    Article  CAS  PubMed  Google Scholar 

  14. Stothard P: The Sequence Manipulation Suite: JavaScript programs for analyzing and formatting protein and DNA sequences. Biotechniques 2000, 28: 1102–1104.

    CAS  PubMed  Google Scholar 

  15. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Huang H, Shukla HD, Wu C, Saxena S: Challenges and Solutions in Proteomics. Curr Genomics 2007, 8: 21–28. 10.2174/138920207780076910

    Article  CAS  Google Scholar 

  17. The UniProt Consortium: The Universal Protein Resource (UniProt) 2009. Nucleic Acids Res 2009, 37: D169–174. 10.1093/nar/gkn664

    Article  PubMed Central  Google Scholar 

  18. The UniProt Consortium: The Universal Protein Resource (UniProt) in 2010. Nucleic Acids Res 2010, 38: D142–148. 10.1093/nar/gkp846

    Article  PubMed Central  Google Scholar 

  19. Kreil DP, Ouzounis CA: Identification of thermophilic species by the amino acid compositions deduced from their genomes. Nucleic Acids Res 2001, 29: 1608–1615. 10.1093/nar/29.7.1608

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Slesarev AI, Mezhevaya KV, Makarova KS, Polushin NN, Shcherbinina OV, Shakhova VV, Belova GI, Aravind L, Natale DA, Rogozin IB, Tatusov RL, Wolf YI, Stetter KO, Malykh AG, Koonin EV, Kozyavkin SA: The complete genome of hyperthermophile Methanopyrus kandleri AV19 and monophyly of archaeal methanogens. Proc Natl Acad Sci USA 2002, 99: 4644–4649. 10.1073/pnas.032671499

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Lobry JR, Chessel D: Internal correspondence analysis of codon and amino-acid usage in thermophilic bacteria. J Appl Genet 2003, 44: 235–261.

    PubMed  Google Scholar 

  22. Peer I, Felder CE, Man O, Silman I, Sussman JL, Beckmann JS: Proteomic signatures: Amino acid and oligopeptide compositions differentiate among phyla. Proteins Struct, Funct Bioinform 2004, 54: 20–40. 10.1002/prot.10559

    Article  CAS  Google Scholar 

  23. Schwartz R, Ting CS, King J: Whole proteome pI values correlate with subcellular localizations of proteins for organisms within the three domains of life. Genome Res 2001, 11: 703–709. 10.1101/gr.GR-1587R

    Article  CAS  PubMed  Google Scholar 

  24. Knight CG, Kassen R, Hebestreit H, Rainey PB: Global analysis of predicted proteomes: Functional adaptation of physical properties. Proc Natl Acad Sci USA 2004, 101: 8390–8395. 10.1073/pnas.0307270101

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Akashi H: Gene expression and molecular evolution. Curr Opin Genet Dev 2001, 11: 660–666. 10.1016/S0959-437X(00)00250-1

    Article  CAS  PubMed  Google Scholar 

  26. Akashi H, Gojobori T: Metabolic efficiency and amino acid composition in the proteomes of Escherichia Coli and Bacillus subtilis. Proc Natl Acad Sci USA 2002, 99: 3695–3700. 10.1073/pnas.062526999

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Kusalik A, Trost B, Bickis M, Fasano C, Capone G, Kanduc D: Codon number shapes peptide redundancy in the universal proteome composition. Peptides 2009, 10: 1940–1944. 10.1016/j.peptides.2009.06.035

    Article  Google Scholar 

  28. Brooks DJ, Fresco JR, Lesk AM, Singh M: Evolution of amino acid frequencies in proteins over deep time: inferred order of introduction of amino acids into the genetic code. Mol Biol Evol 2002, 19: 1645–1655.

    Article  CAS  PubMed  Google Scholar 

  29. Trost B, Kanduc D, Kusalik A: Rare peptide segments are found significantly more often in proto-oncoproteins than control proteins: implications for immunology and oncology. J R Soc Interface 2009, 6: 123–127. 10.1098/rsif.2008.0320

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kanduc D: Protein information content resides in rare peptide segments. Peptides 2010, 31: 983–988. 10.1016/j.peptides.2010.02.003

    Article  CAS  PubMed  Google Scholar 

  31. Wootton JC: Sequences with 'unusual' amino acid composition. Curr Opin Struct Biol 1994, 4: 413–421. 10.1016/S0959-440X(94)90111-2

    Article  CAS  Google Scholar 

  32. Willers J, Lucchese A, Kanduc D, Ferrone S: Molecular mimicry of phage displayed peptides mimicking GD3 ganglioside. Peptides 1999, 20: 1021–1026. 10.1016/S0196-9781(99)00095-9

    Article  CAS  PubMed  Google Scholar 

  33. Natale C, Giannini T, Lucchese A, Kanduc D: Computer-assisted analysis of molecular mimicry between HPV16 E7 oncoprotein and human protein sequences. Immunol Cell Biol 2000, 78: 580–585. 10.1046/j.1440-1711.2000.00949.x

    Article  CAS  PubMed  Google Scholar 

  34. Kanduc D: Peptimmunology: immunogenic peptides and sequence redundancy. Curr Drug Discov Technol 2005, 2: 239–244. 10.2174/157016305775202946

    Article  CAS  PubMed  Google Scholar 

  35. Kanduc D: Defining peptide sequences: from antigenicity to immunogenicity through redundancy. Curr Pharmacogenomics 2006, 4: 33–37. 10.2174/157016006776055374

    Article  CAS  Google Scholar 

  36. Kanduc D: Correlating low-similarity peptide sequences and allergenic epitopes. Curr Pharm Des 2008, 14: 289–295. 10.2174/138161208783413257

    Article  CAS  PubMed  Google Scholar 

  37. Kanduc D: Immunogenicity in peptide-immunotherapy: from self/nonself to similar/dissimilar sequences. In Multichain Immune Recognition Receptor Signaling: From Spatiotemporal Organization to Human Disease. Landes Biosci. Edited by: Sigalov A. Austin, TX, USA; 2008:198–207. full_text

    Chapter  Google Scholar 

  38. Kanduc D: Self-nonself peptides in the design of vaccines. Curr Pharm Des 2009, 15: 3283–3289. 10.2174/138161209789105135

    Article  CAS  PubMed  Google Scholar 

  39. Lucchese G, Stufano A, Kanduc D: Proteome-guided search for influenza A B-cell epitopes. FEMS Immunol Med Microbiol 2009, 57: 88–92. 10.1111/j.1574-695X.2009.00582.x

    Article  CAS  PubMed  Google Scholar 

  40. Lucchese A, Serpico R, Crincoli V, Shoenfeld Y, Kanduc D: Sequence uniqueness as a molecular signature of HIV-1-derived B-cell epitopes. Int J Immunopathol Pharmacol 2009, 22: 639–646.

    CAS  PubMed  Google Scholar 

  41. Kanduc D: Epitopic peptides with low similarity to the host proteome: towards biological therapies without side effects. Expert Opin Biol Ther 2009, 9: 45–53. 10.1517/14712590802614041

    Article  CAS  PubMed  Google Scholar 

  42. Mandavilli A: When the vaccine causes disease. Nat Med 2007, 13: 274. 10.1038/nm0307-274b

    Article  CAS  PubMed  Google Scholar 

  43. Kanduc D: Penta- and hexapeptide sharing between HPV16 and Homo sapiens proteomes. Int J Med Sci 2009, 1: 387.

    Google Scholar 

  44. Kanduc D: Quantifying the possible cross-reactivity risk of an HPV16 vaccine. J Exp Ther Oncol 2009, 8: 65–76.

    CAS  PubMed  Google Scholar 

  45. Ricco R, Kanduc D: Hepatitis B virus and Homo sapiens proteome-wide analysis: A profusion of viral peptide overlaps in neuron-specific human proteins. Biologics 2010, 4: 75–81.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Partially funded by: the National Sciences and Engineering Research Council of Canada (BT and AK) and Ministry of University, Italy (DK). GN and CF are PhD students of the course "Analytical Morphometry and Models of Molecular Medicine".

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Darja Kanduc.

Additional information

Authors' contributions

GC and CF performed the physico-(bio)chemical analyses; GN performed the standard nucleotide-nucleotide BLAST searches as reported in Table 1 and Additional file 1, Table S1; BT and MB performed the mathematical analyses and the statistical treatment of the data. AK performed the pentapeptide decomposition of the universal proteomes and determined the occurrence counts of the pentapeptides. DK proposed the original idea, interpreted the data, developed the research project, and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2010_3840_MOESM1_ESM.DOC

Additional file 1: Table S1. The oligodeoxynucleotide sequences corresponding to never-expressed peptide motifs are mainly located in the non-coding strand. The additional Table shows that the pentadecameric oligodeoxynucleotide sequences coding for the never-expressed pentapeptides correspond to untranslatable, frameshifted or mistranslated oligodeoxynucleotide sequences. (DOC 870 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Capone, G., Novello, G., Fasano, C. et al. The oligodeoxynucleotide sequences corresponding to never-expressed peptide motifs are mainly located in the non-coding strand. BMC Bioinformatics 11, 383 (2010). https://doi.org/10.1186/1471-2105-11-383

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-11-383

Keywords