Skip to main content

Pairwise statistical significance of local sequence alignment using multiple parameter sets and empirical justification of parameter set change penalty

Abstract

Background

Accurate estimation of statistical significance of a pairwise alignment is an important problem in sequence comparison. Recently, a comparative study of pairwise statistical significance with database statistical significance was conducted. In this paper, we extend the earlier work on pairwise statistical significance by incorporating with it the use of multiple parameter sets.

Results

Results for a knowledge discovery application of homology detection reveal that using multiple parameter sets for pairwise statistical significance estimates gives better coverage than using a single parameter set, at least at some error levels. Further, the results of pairwise statistical significance using multiple parameter sets are shown to be significantly better than database statistical significance estimates reported by BLAST and PSI-BLAST, and comparable and at times significantly better than SSEARCH. Using non-zero parameter set change penalty values give better performance than zero penalty.

Conclusion

The fact that the homology detection performance does not degrade when using multiple parameter sets is a strong evidence for the validity of the assumption that the alignment score distribution follows an extreme value distribution even when using multiple parameter sets. Parameter set change penalty is a useful parameter for alignment using multiple parameter sets. Pairwise statistical significance using multiple parameter sets can be effectively used to determine the relatedness of a (or a few) pair(s) of sequences without performing a time-consuming database search.

Background

Local sequence alignment plays a major role in the analysis of DNA and protein sequences [1–3]. It is the basic step of many other applications like detecting homology, finding protein structure and function, deciphering evolutionary relationships, etc. There exist several local sequence alignment programs that use well-known algorithms [4, 5] or their heuristic versions [3, 6, 7]. Database search is a special case of pairwise local sequence alignment where the second sequence is a database in itself consisting of many sequences. Recently, there have been many enhancements in alignment program features [8, 9] using difference blocks and multiple scoring matrices, in an attempt to incorporate more biological features in the alignment algorithm.

Why statistical significance?

The local sequence alignment programs report alignment scores for the alignments constructed, and related (homologous) sequences will have higher alignment scores. But the definition of high depends strongly on the alignment score distribution, which gives importance to the concept of statistical significance. An alignment score is considered statistically significant if it has a low probability of occurring by chance. Since the alignment score distribution depends on various factors like alignment program, scoring scheme, sequence lengths, sequence compositions [10], it implies that it is possible to have two alignments of different sequence pairs with scores x and y with x <y, but x more significant than y. Therefore, instead of using the alignment score for detecting homology, the statistical significance of an alignment score is more widely accepted as a metric to comment on the relatedness of the two sequences being aligned. Of course, it is important to emphasize here that although statistical significance is a good preliminary indicator of biological significance, it does not necessarily imply biological significance [10, 11].

Accurate statistical theory for the ungapped alignment score distribution is available [12]. However, no precise statistical theory currently exists for the gapped alignment score distribution and for score distributions from alignment programs using additional features. Accurate estimation of statistical significance of gapped sequence alignment scores has attracted a lot of attention in the recent years [13–21].

Although there exists some understanding of the statistics of gapped alignment score distributions for simple scoring schemes [22, 23], but a complete mathematical description of the optimal score distribution remains far from reach [23]. There exist some excellent reviews on statistical significance in sequence comparison in the literature [10, 24–26].

Database statistical significance versus pairwise statistical significance

Recently, a study of pairwise statistical significance and its comparison with database statistical significance [27] was conducted. In summary, the database statistical significance reported by most database search programs like SSEARCH, FASTA, PSI-BLAST is database-dependent, and hence, the same alignment of two sequences with the same alignment score can be evaluated as having different significance values in database searches with different databases, and even with the same database at different times, as the database size can be variable. On the other hand, pairwise statistical significance is specific to the sequence pair being aligned, and is database-independent. In [27], various approaches to estimate pairwise statistical significance were compared to find that maximum likelihood fitting with censoring left of peak is the most accurate method for estimating pairwise statistical significance. Further, this method was compared with database statistical significance in a homology detection experiment to find that pairwise statistical significance performs comparably to and sometimes significantly better than database statistical significance.

Accurate statistical significance estimates for pairwise alignments can be very useful to comment on the relatedness of a pair of sequences aligned by an alignment program independent of any database. And thus, it can also be used to compare different combination of alignment parameters – like the alignment program itself, substitution matrices, gap costs. A comparison of different gap opening penalties for four commonly used BLOSUM matrices using pairwise statistical significance was presented in [28]. In addition to the standard local alignment algorithms [4, 5], some recent algorithms have been developed [8, 9] that take into account other desirable biological features in addition to gaps – like difference blocks or the use of multiple parameter sets (substitution matrices, gap penalties). These features of the alignment programs enhance the sequence alignment of real sequences by better suiting to different conservation rates at different spatial locations of the sequences. As pointed out earlier, accurate statistical theory for alignment score distribution is available only for ungapped alignment, and not even for its simplest extension, i.e., alignment with gaps. Accurate statistics of the alignment score distribution from newer and more sophisticated alignment programs therefore is not expected to be straightforward. For comparing the performance of newer alignment programs, accurate estimates of pairwise statistical significance can be very useful. Further, quick and accurate estimates of pairwise statistical significance can also be helpful for applications like multiple sequence alignment and phylogenetic tree construction which are based on pairwise sequence alignment to select most related pairs of sequences, for example, in a progressive multiple sequence alignment.

The extreme value distribution for ungapped and gapped alignments

Just as the distribution of the sum of a large number of independent identically distributed (i.i.d.) random variables tends to a normal distribution (central limit theorem), the distribution of the maximum of a large number of i.i.d. random variables tends to an extreme value distribution (EVD) [29]. This is an important and useful fact, because in principle it allows us to fit an EVD to the score distribution from any local alignment program and use it for estimating statistical significance of scores from that program. The distribution of Smith-Waterman local alignment score between random, unrelated sequences is approximately a Gumbel-type EVD [12]. In the limit of sufficiently large sequence lengths m and n, the statistics of HSP (High-scoring Segment Pairs which correspond to ungapped local alignments) scores are characterized by two parameters, K and λ. The probability (P-value) that the optimal local alignment score S exceeds x is estimated by:

Pr(S > x) ~1 - e-E,

where E is the E-value and is given by

E = Kmne-λx.

For E-values less than 0.01, both E-value and P-values are very close to each other. The above formulae are valid for ungapped alignments [12], and the parameters K and λ can be computed analytically from the substitution scores and sequence compositions. For gapped alignments, no perfect statistical theory has yet been developed, although there exist some good starting points for the problem as mentioned before [22, 23]. Recently, researchers have also looked closely at the low probability tail distribution, and the work in [30] applied a rare-event sampling technique earlier used in [31] and suggested a Gaussian correction to the Gumbel distribution to better describe the rare event tail, resulting in a considerable change in the reported significance values. However, for most practical purposes, the original Gumbel distribution has been widely used to describe gapped alignment score distribution [9, 13–15, 17, 27, 32, 33]. From an empirically generated score distribution, we can directly observe the E-value E for a particular score x, by counting the number of times a score x or higher was attained. Since this number would be different for different number of random shuffles N (or number of sequences in the database in case of database search), a normalized E-value is defined as

E n o r m a l i z e d = E N MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyrau0aaSbaaSqaaiabd6gaUjabd+gaVjabdkhaYjabd2gaTjabdggaHjabdYgaSjabdMgaPjabdQha6jabdwgaLjabdsgaKbqabaGccqGH9aqpjuaGdaWcaaqaaiabdweafbqaaiabd6eaobaaaaa@3F0C@

In theory, this normalized E-value is same as the P-value (for large N).

Contributions

In this paper, we extend the existing work on pairwise statistical significance [27] to incorporate in it the use of multiple parameter sets, and evaluate it on an important knowledge discovery application-homology detection. We conducted similar experiments as reported in [34], and later in [27] on a subset of the CATH 2.3 database to compare pairwise statistical significance with single and multiple parameter sets. This benchmark database was earlier created in [34] to evaluate seven protein structure comparison methods and two sequence comparison programs: SSEARCH and PSI-BLAST. SSEARCH uses the original Smith-Waterman algorithm [4], and is considered the most sensitive algorithm in terms of retrieval accuracy, better than the heuristic versions like BLAST and FASTA [35, 36]. PSI-BLAST is a modification to the BLAST program, where position specific scoring matrices are constructed over multiple iterations of BLAST algorithm. Comparison of pairwise statistical significance results using multiple parameter sets with pairwise statistical significance using a single parameter set shows that at least for some error levels, using multiple parameter sets is significantly better than using a single parameter set. This is because sequences can have different conservation rates at different spatial locations, which can be better aligned using multiple parameter sets (substitution matrices, gap penalties, etc.). Comparison with database statistical significance results show that pairwise statistical significance with multiple parameter sets gives significantly better performance than the statistical significance estimates reported by BLAST and PSI-BLAST, and comparable and at times significantly better performance than the SSEARCH program. Further, the results also give concrete evidence that for the practical application of homology detection, the score distribution from alignment program using multiple parameter sets can also be assumed to follow an extreme value distribution. This is an important and useful finding since it is in general difficult to accurately determine statistics of alignment scores from enhanced alignment programs. Finally, experiments with different values of parameter set change penalties indicate that it is indeed important to use a non-zero parameter set change penalty while performing alignment using multiple parameter sets.

Methods

Pairwise statistical significance estimation

Consider the pairwise statistical significance described in [27] to be obtainable by the following function: PairwiseStatSig(Seq 1, Seq 2, SC, N) where Seq 1 is the first sequence, Seq 2 is the second sequence, SC is the scoring scheme (substitution matrix, gap penalties), and N is the number of shuffles. The function PairwiseStatSig, therefore generates a score distribution by aligning Seq 1 with N shuffled versions of Seq 2, fits the distribution to an extreme value distribution using censored maximum likelihood fitting to obtain the statistical parameters K and λ, and returns the pairwise statistical significance estimate of the pairwise alignment score between Seq 1 and Seq 2 using the parameters K and λ in the P-value formula. More details on pairwise statistical significance can be found in [27]. In this paper, we dynamically use multiple parameter sets instead of a single scoring scheme SC for estimation of pairwise statistical significance.

Dynamic use of multiple parameter sets in sequence alignment

Usually, pairwise sequence alignment is done with a single parameter set (substitution matrix, gap penalties). But to suit the different levels of conservation between sequences, there exists an algorithm [9] which can dynamically use multiple parameter sets and generate a single optimal alignment with possibly different parameter sets used in different regions of the alignment. The algorithm is implemented in a program named GAP4. The algorithm uses a dynamic programming approach as explained in [9]. Consider alignment of two sequences A = a1, a2,..., a m and B = b1, b2,... b n using p parameter sets P1, P2,..., P p . Let A i and B j be the subsequences a1, a2,... a i and b1, b2,..., b j respectively. For each alignment position (i, j) and each parameter set P k , the algorithm keeps track of the optimal alignment score of the subsequences A i and B j where the last component (substitution, gap, or difference block) is scored using P k . Dynamic programming is used to get optimal alignment for progressive alignment positions, until i becomes m and j becomes n. Appropriate modification of the algorithm also allows it to calculate the optimal local alignment. More details about using multiple parameter sets for pairwise sequence alignment can be found in [9].

Evaluation methodology

To evaluate the performance of pairwise statistical significance using multiple parameter sets, we used a non-redundant subset of the CATH 2.3 database (Class, Architecture, Topology, and Hierarchy, [37]) provided by [34] and available at ftp://ftp.ebi.ac.uk/pub/software/unix/fasta/prot_sci_04/. This database was selected in [34] to evaluate seven structure comparison programs and two sequence comparison programs. As described in [34], this dataset consists of 2771 domain sequences and includes 86 selected test query sequences. This domain set is considered as a valid benchmark for testing protein comparison algorithms [38].

We used this database and query set for experimenting with pairwise statistical significance using multiple parameter sets. For each of the 86 × 2771 comparisons, we used the maximum likelihood method with censoring left of peak with 1000 shuffles to fit the score distribution from the GAP4 program with substitution matrices BLOSUM45, BLOSUM50, BLOSUM62, BLOSUM80, and their all possible combinations (24 - 1 = 15 in number). All matrices were used in 1/3 bit scale. The gap opening penalties for each of these matrices was set to the values empirically determined to be the best for this database in [28]. These are listed in Table 1. The gap extension penalties were set to 2 for all the four matrices. Following [27, 34], Error per Query (EPQ) versus Coverage plots were used to present the results. To create these plots, the list of pairwise comparisons was sorted based on decreasing statistical significance (increasing P-values). Going down the list, the coverage count is increased by one if the two sequences of the pair are homologs, and the error count is increased by one if they are not. At a given point in the list, Errors Per Query (EPQ) is the total number of errors incurred so far, divided by the number of queries; and coverage is the fraction of total homolog pairs so far detected. In the ideal case, the curve would go from 0% to 100% coverage, without incurring any errors, which would correspond to a straight line on the x-axis. Therefore, the more the curve is towards the right, the better the curve is.

Table 1 Effective gap opening penalties for commonly used BLOSUM matrices determined for the benchmark database used

Just as gap opening and gap extension penalties are dynamically charged during the alignment process whenever a gap is inserted and extended respectively, the GAP4 [9] program allows the use of a parameter set change penalty, which is dynamically charged whenever the parameter set mapping is changed during the alignment process. To see the effect of parameter set change penalty on the coverage performance, we conducted a series of homology detection experiments with one of the substitution matrix combinations (BLOSUM45 and BLOSUM62) with different parameter set change penalties. Coverage vs. parameter set change penalty curves were plotted at different error levels to find the usefulness of the parameter set change penalty, as reported in the next section.

Results

Comparison with pairwise statistical significance using single parameter set

Out of the 15 substitution matrix combinations, 4 are using single parameter sets, 6 are using two parameter sets, 4 are using three parameter sets, and 1 is using all four parameter sets. The EPQ vs. Coverage curves using pairwise statistical significance with two, three, and four parameter sets are presented in Figures 1, 2, and 3 respectively. For comparison purposes, the EPQ vs. Coverage curves using corresponding single parameter sets are also presented in the same figures. The y-axis (error-axis) in all these graphs is in log-scale, and hence there is more information in the upper part of the graphs. These figures suggest that pairwise statistical significance using multiple parameter sets performs comparably to and sometimes significantly better than pairwise statistical significance using a single parameter set for most instances of using a single parameter set, and at most error levels.

Figure 1
figure 1

Pairwise statistical significance using two parameter sets. Errors per Query vs. Coverage plot for pairwise statistical significance using two parameter sets, along with the curves using corresponding single parameter sets. (a) BLOSUM45, BLOSUM50; (b) BLOSUM45, BLOSUM62; (c) BLOSUM45, BLOSUM80; (d) BLOSUM50, BLOSUM62; (e) BLOSUM50, BLOSUM80; (f) BLOSUM62, BLOSUM80. In 5 panels (b) through (f) out of 6, using two parameter sets leads to better coverage than using a single parameter set at most error levels.

Figure 2
figure 2

Pairwise statistical significance using three parameter sets. Errors per Query vs. Coverage plot for pairwise statistical significance using three parameter sets, along with the curves using corresponding single parameter sets. (a) BLOSUM45, BLOSUM50, BLOSUM62; (b) BLOSUM45, BLOSUM50, BLOSUM80; (c) BLOSUM45, BLOSUM62, BLOSUM80; (d) BLOSUM50, BLOSUM62, BLOSUM80. In all 4 panels, using three parameter sets leads to better coverage than using a single parameter set at most error levels for at least two instances of using a single parameter set.

Figure 3
figure 3

Pairwise statistical significance using four parameter sets. Errors per Query vs. Coverage plot for pairwise statistical significance using four parameter sets BLOSUM45, BLOSUM50, BLOSUM62, BLOSUM80 along with the curves using corresponding single parameter sets. Using four parameter sets results in better coverage than using a single parameter set at most error levels for at least three instances of using a single parameter set.

Comparison with database statistical significance

Since the EPQ vs. Coverage curves on the complete dataset can be distorted due to poor performance by one or two queries (if those queries produce many errors at low coverage levels) [34], to compare the performance of pairwise statistical significance using multiple parameter sets with database statistical significance, we examined the performance of the methods with individual queries, following the work in [34]. The coverage of each of the 86 queries at the 1st, 3rd, 10th, 30th, and 100th error was recorded, and the median coverage for each error level across the 86 queries was plotted to obtain EPQ vs. Coverage curves for the sequence comparison method to be evaluated. Figure 4 shows the median coverage level at the 1st, 3rd, 10th, 30th, and 100th false positive for homologs (i.e. 43 of the queries have worse coverage, and 43 have better coverage). The curve for SSEARCH in Figure 4 is derived from the figure 2A in [34]. The curves for BLAST and PSI-BLAST were obtained by experimentation. It is clear that the proposed pairwise statistical significance using multiple parameter sets performs significantly better than BLAST and PSI-BLAST at all error levels, comparable to SSEARCH at low error levels, and significantly better than SSEARCH at higher error levels.

Figure 4
figure 4

Comparison with database statistical significance. Comparison of pairwise statistical significance (using multiple parameter sets) and database statistical significance. All parameter combinations are significantly better than database statistical significance estimates reported by BLAST and PSI-BLAST at all error levels, and better than SSEARCH especially at higher error levels.

Empirical justification of parameter set change penalty

The coverage vs. parameter set change penalty plot for the substitution matrix combination of BLOSUM45 and BLOSUM62 is illustrated in Figure 5. The curve shows a poor coverage performance for the case when the parameter set change penalty is not charged, i.e., when the alignment algorithm is freely allowed to change the parameter set during alignment without charging any penalty. This can be explained by the fact that the algorithm would try to mathematically maximize the alignment score by changing the parameter set as frequently as possible, which may produce more biologically irrelevant alignments. A similar phenomenon is also observed when very low gap penalty is used [28]. The coverage performance clearly improves for non-zero values of parameter set change penalty, which provides its empirical justification.

Figure 5
figure 5

Empirical justification of parameter set change penalty. Coverage vs. Parameter Set Change Penalty plots at different errors per query for the substitution matrix combination of BLOSUM45 and BLOSUM62. Poor coverage is obtained if the parameter set change penalty is zero. The coverage is better and steady for non-zero values of parameter set change penalty.

Discussion

As pointed out earlier, SSEARCH employs the original Smith-Waterman algorithm for alignment, and is considered more sensitive than its heuristic implementations like BLAST and FASTA. PSI-BLAST uses an iterative approach with query-specific substitution matrices, and its performance mainly depends on the quality of position-specific scoring matrices (PSSMs) constructed iteratively. The results show that PSI-BLAST gave poorer performance than pairwise statistical significance using multiple parameter sets, even with PSSMs constructed against the benchmark CATH database used in our experiments. However, using PSSMs derived against BLAST non-redundant protein database has been shown to give better results [34] as it uses much more information than just a pair of sequences. Comparable and at times significantly better results than SSEARCH using pairwise statistical significance with multiple parameter sets implies that statistical significance estimates at least as good as database statistical significance can be obtained by pairwise statistical significance using multiple parameter sets without having to do a time-consuming database search. This can be very useful to estimate accurate pairwise statistical significance of two (or a few) sequences, which is a common scenario in many pairwise alignment based applications like phylogenetic tree construction, progressive multiple sequence alignment.

It is important to note that the proposed method is not a database search method but statistical significance estimation method for pairwise local alignment, and the comparison with database search programs like BLAST, SSEARCH, and PSI-BLAST is of their statistical significance estimation strategies. The proposed method, as of now is not scalable to a database search, but can be used to refine the results from a fast database search program like BLAST.

Since pairwise alignment using multiple parameter sets takes more computational time than using a single parameter set, pairwise statistical significance estimation using multiple parameter sets also takes more time than pairwise statistical significance estimation using a single parameter set. In general, using k parameter sets increases the computation time by a factor little more than k. Therefore, faster methods for significance estimation can be very helpful.

Conclusion

This paper extends the work on pairwise statistical significance by incorporating in it the use of multiple parameter sets (substitution matrices, gap penalties, etc.), and compares it with database statistical significance for the knowledge discovery application of homology detection. The results show that pairwise statistical significance using multiple parameter sets performs better than pairwise statistical significance using a single parameter set. It also performs significantly better than database statistical significance using BLAST and PSI-BLAST, and comparable and at times significantly better than database statistical significance using SSEARCH. Further, an empirical justification of the use of parameter set change penalty is provided.

Since PSI-BLAST results can be improved by using better quality PSSMs derived from larger protein databases, we believe that the performance of pairwise statistical significance can also be improved using sequence-specific/position-specific substitution matrices, which is a significant part of our future work. Another important contribution can be to estimate the pairwise statistical significance accurately in less time, since using multiple parameter sets increases the significance estimation time. Faster methods for determining pairwise statistical significance would be very useful.

References

  1. Pearson WR, Lipman DJ: Improved Tools for Biological Sequence Comparison. Proc Natl Acad Sci U S A. 1988, 85 (8): 2444-2448.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool. Journal of Molecular Biology. 1990, 215 (3): 403-410.

    Article  CAS  PubMed  Google Scholar 

  3. 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 Research. 1997, 25 (17): 3389-3402.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Smith TF, Waterman MS: Identification of Common Molecular Subsequences. Journal of Molecular Biology. 1981, 147: 195-197.

    Article  CAS  PubMed  Google Scholar 

  5. Sellers PH: Pattern Recognition in Genetic Sequences by Mismatch Density. Bulletin of Mathematical Biology. 1984, 46 (4): 501-514. [http://www.springerlink.com/content/2v4477481102w030]

    Article  CAS  Google Scholar 

  6. Pearson WR: Effective Protein Sequence Comparison. Methods in Enzymology. 1996, 266: 227-259.

    Article  CAS  PubMed  Google Scholar 

  7. Pearson WR: Flexible Sequence Similarity Searching with the FASTA3 Program Package. Methods in Molecular Biology. 2000, 132: 185-219.

    CAS  PubMed  Google Scholar 

  8. Huang X, Chao KM: A Generalized Global Alignment Algorithm. Bioinformatics. 2003, 19 (2): 228-233.

    Article  CAS  PubMed  Google Scholar 

  9. Huang X, Brutlag DL: Dynamic Use of Multiple Parameter Sets in Sequence Alignment. Nucleic Acids Research. 2007, 35 (2): 678-686.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Mott R: Alignment: Statistical Significance. Encyclopedia of Life Sciences. 2005, [Doi.wiley.com/10.1038/npg.els.0005264].

    Google Scholar 

  11. Altschul SF, Boguski MS, Gish W, Wootton JC: Issues in searching molecular sequence databases. Nature Genetics. 1994, 6 (2): 119-129.

    Article  CAS  PubMed  Google Scholar 

  12. Karlin S, Altschul SF: Methods for Assessing the Statistical Significance of Molecular Sequence Features by Using General Scoring Schemes. Proc Natl Acad Sci USA. 1990, 87 (6): 2264-2268.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Waterman MS, Vingron M: Rapid and Accurate Estimates of Statistical Significance for Sequence Database Searches. Proc Natl Acad Sci U S A. 1994, 91 (11): 4625-4628.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Altschul SF, Gish W: Local Alignment Statistics. Methods in Enzymology. 1996, 266: 460-80.

    Article  CAS  PubMed  Google Scholar 

  15. Pearson WR: Empirical Statistical Estimates for Sequence Similarity Searches. Journal of Molecular Biology. 1998, 276: 71-84.

    Article  CAS  PubMed  Google Scholar 

  16. Mott R, Tribe R: Approximate Statistics of Gapped Alignments. Journal of Computational Biology. 1999, 6: 91-112.

    Article  CAS  PubMed  Google Scholar 

  17. Mott R: Accurate Formula for P-values of Gapped Local Sequence and Profile Alignments. Journal of Molecular Biology. 2000, 300: 649-659.

    Article  CAS  PubMed  Google Scholar 

  18. Altschul SF, Bundschuh R, Olsen R, Hwa T: The estimation of statistical parameters for local alignment score distributions. Nucleic Acids Research. 2001, 29 (2): 351-361.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Schäffer AA, Aravind L, Madden TL, Shavirin S, Spouge JL, Wolf YI, Koonin EV, Altschul SF: Improving the Accuracy of PSI-BLAST Protein Database Searches with Composition-based Statistics and Other Refinements. Nucleic Acids Research. 2001, 29 (14): 2994-3005.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Sheetlin S, Park Y, Spouge JL: The Gumbel Pre-factor k for Gapped Local Alignment can be Estimated From Simulations of Global Alignment. Nucleic Acids Research. 2005, 33 (15): 4987-4994.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Yu YK, Gertz EM, Agarwala R, Schäffer AA, Altschul SF: Retrieval Accuracy, Statistical Significance and Compositional Similarity in Protein Sequence Database Searches. Nucleic Acids Research. 2006, 34 (20): 5966-5973.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Kschischo M, Lässig M, Yuc YK: Toward an Accurate Statistics of Gapped Alignments. Bulletin of Mathematical Biology. 2004, 67: 169-191.

    Article  Google Scholar 

  23. Grossmann S, Yakir B: Large Deviations for Global Maxima of Independent Superadditive Processes with Negative Drift and an Application to Optimal Sequence Alignments. Bernoulli. 2004, 10 (5): 829-845.

    Article  Google Scholar 

  24. Pagni M, Jongeneel CV: Making Sense of Score Statistics for Sequence Alignments. Briefings in Bioinformatics. 2001, 2: 51-67.

    Article  CAS  PubMed  Google Scholar 

  25. Pearson WR, Wood TC: Statistical Significance in Biological Sequence Comparison. Handbook of Statistical Genetics. Edited by: Balding DJ, Bishop M, Cannings C. 2001, Chichester, UK: Wiley, 39-66.

    Google Scholar 

  26. Mitrophanov AY, Borodovsky M: Statistical Significance in Biological Sequence Analysis. Briefings in Bioinformatics. 2006, 7: 2-24.

    Article  CAS  PubMed  Google Scholar 

  27. Agrawal A, Brendel V, Huang X: Pairwise Statistical Significance Versus Database Statistical Significance for Local Alignment of Protein Sequences. Bioinformatics Research and Applications. 2008, LNCS(LNBI), Springer Berlin/Heidelberg, 4983: 50-61.

    Chapter  Google Scholar 

  28. Agrawal A, Brendel V, Huang X: Pairwise Statistical Significance and Empirical Determination of Effective Gap Opening Penalties for Protein Local Sequence Alignment. International Journal of Computational Biology and Drug Design. 2008, 1 (4): 347-367.

    Article  CAS  PubMed  Google Scholar 

  29. Kotz S, Nadarajah S: Extreme Value Distributions: Theory and Applications. 2000, Imperial College Press, London, UK, chap 1: 3-4.

    Book  Google Scholar 

  30. Wolfsheimer S, Burghardt B, Hartmann AK: Local Sequence Alignments Statistics: Deviations from Gumbel Statistics in the Rare-event Tail. Algorithms Mol Biol. 2007, 2: 9-

    Article  PubMed Central  PubMed  Google Scholar 

  31. Hartmann AK: Sampling Rare Events: Statistics of Local Sequence Alignments. Physical Review E. 2002, 65 (5): 056102-

    Article  Google Scholar 

  32. Olsen R, Bundschuh R, Hwa T: Rapid Assessment of Extremal Statistics for Gapped Local Alignment. Proc of the Seventh International Conference on Intelligent Systems for Molecular Biology. 1999, AAAI Press, 211-222.

    Google Scholar 

  33. Mott RF: Maximum-likelihood Estimation of the Statistical Distribution of SmithWaterman Local Sequence Similarity Scores. Bulletin of Mathematical Biology. 1992, 54: 59-75.

    Article  CAS  PubMed  Google Scholar 

  34. Sierk ML, Pearson WR: Sensitivity and Selectivity in Protein Structure Comparison. Protein Science. 2004, 13 (3): 773-785.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Brenner SE: Practical database searching. Trends in Biotechnology. 1998, 16: 9-12.

    Article  Google Scholar 

  36. Bucher P, Hofmann K: A Sequence Similarity Search Algorithm Based on a Probabilistic Interpretation of an Alignment Scoring System. Proceedings of the Fourth International Conference on Intelligent Systems for Molecular Biology. 1996, AAAI Press, 44-51.

    Google Scholar 

  37. Orengo CA, Michie AD, Jones S, Jones DT, Swindells MB, Thornton JM: CATH – A Hierarchic Classification of Protein Domain Structures. Structure. 1997, 28: 1093-1108.

    Article  Google Scholar 

  38. Rocha J, Rosselló F, Segura J: Compression Ratios Based on the Universal Similarity Metric Still Yield Protein Distances far from CATH Distances. CoRR. 2006, abs/q-bio/0603007.

    Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Volker Brendel for helpful discussions and providing links to the data.

This article has been published as part of BMC Bioinformatics Volume 10 Supplement 3, 2009: Second International Workshop on Data and Text Mining in Bioinformatics (DTMBio) 2008. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/10?issue=S3.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ankit Agrawal.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AA conceived the study based on the initial idea given by XH. Both AA and XH did the programming, AA carried out the experiments and analysis, and drafted the initial manuscript. Both AA and XH read and approved the final version of the manuscript.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://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

Agrawal, A., Huang, X. Pairwise statistical significance of local sequence alignment using multiple parameter sets and empirical justification of parameter set change penalty. BMC Bioinformatics 10 (Suppl 3), S1 (2009). https://doi.org/10.1186/1471-2105-10-S3-S1

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-10-S3-S1

Keywords