The construction and use of log-odds substitution scores for multiple sequence alignment

PLoS Computational Biology

Volumn 6, Issue 7, 2010, Pages 11-

  Altschul, Stephen F ^a   Wootton, John C ^a   Zaslavsky, Elena ^b   Yu, Yi Kuo ^a  

^a NATIONAL INSTITUTES OF HEALTH   (United States)

^b MOUNT SINAI SCHOOL OF MEDICINE   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

BAYESIAN NETWORKS; BIOINFORMATICS;

AMINO-ACIDS; BAYESIAN; BAYESIAN METHODS; GIBBS SAMPLING; INDIVIDUAL SEQUENCES; MULTIPLE ALIGNMENT; MULTIPLE SEQUENCE ALIGNMENTS; PAIRWISE ALIGNMENT; PAIRWISE SEQUENCE ALIGNMENT; PRIOR DISTRIBUTION;

HIDDEN MARKOV MODELS;

ACCURACY; ARTICLE; BAYES THEOREM; COMPUTER PREDICTION; COMPUTER PROGRAM; CONTROLLED STUDY; DNA BINDING; HIDDEN MARKOV MODEL; INFORMATION PROCESSING; INFORMATION RETRIEVAL; NONHUMAN; PLASMODIUM FALCIPARUM; PROBABILITY; PROTEIN STRUCTURE; SCORING SYSTEM; SEQUENCE ALIGNMENT; TOXOPLASMA GONDII; ALGORITHM; AMINO ACID SEQUENCE; AUTOMATED PATTERN RECOGNITION; BIOLOGY; CHEMISTRY; CONSENSUS SEQUENCE; GENETIC DATABASE; GENETICS; METHODOLOGY; NUCLEOTIDE SEQUENCE; PLASMODIUM; PROTEIN TERTIARY STRUCTURE; STATISTICAL MODEL; TOXOPLASMA;

PLASMODIUM FALCIPARUM; TOXOPLASMA GONDII;

DNA BINDING PROTEIN;

ALGORITHMS; AMINO ACID SEQUENCE; BASE SEQUENCE; BAYES THEOREM; COMPUTATIONAL BIOLOGY; CONSENSUS SEQUENCE; DATABASES, GENETIC; DNA-BINDING PROTEINS; MODELS, STATISTICAL; PATTERN RECOGNITION, AUTOMATED; PLASMODIUM; PROTEIN STRUCTURE, TERTIARY; SEQUENCE ALIGNMENT; TOXOPLASMA;

EID: 77957880647     PISSN: 1553734X     EISSN: 15537358     Source Type: Journal    
DOI: 10.1371/journal.pcbi.1000852     Document Type: Article

References (122)

1. Smith TF, Waterman MS (1981) Identification of common molecular subsequences. J Mol Biol 147: 195-197.
2. Sellers PH (1984) Pattern recognition in genetic sequences by mismatch density. Bull Math Biol 46: 501-514.
3. Pearson WR, Lipman DJ (1988) Improved tools for biological sequence comparison. Proc Natl Acad Sci USA 85: 2444-2448.
4. Altschul SF, Madden TL, Schä ffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389-3402.
5. Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. (2009) CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res 37: D205-D210.
6. Kann MG, Sheetlin SL, Park Y, Bryant SH, Spouge JL (2007) The identification of complete domains within protein sequences using accurate e-values for semi-global alignment. Nucleic Acids Res 35: 4678-4685.
7. Dayhoff MO, Schwartz RM, Orcutt BC (1978) A model of evolutionary change in proteins. In: Dayhoff MO, ed. Atlas of Protein Sequence and Structure. Washington, DC: Natl. Biomed. Res. Found., volume 5, suppl. 3. pp 345-352.
8. Schwartz RM, Dayhoff MO (1978) Matrices for detecting distant relationships. In: Dayhoff MO, ed. Atlas of Protein Sequence and Structure. Washington, DC: Natl. Biomed. Res. Found., volume 5, suppl. 3. pp 353-358.
9. Feng DF, Johnson MS, Doolittle RF (1985) Aligning amino acid sequences: comparison of commonly used methods. J Mol Evol 21: 112-125.
10. Taylor WR (1986) The classification of amino acid conservation. J Theor Biol 119: 205-218.
11. Rao JKM (1987) New scoring matrix for amino acid residue exchanges based on residue characteristic physical parameters. Int J Peptide Protein Res 29: 276-281.
12. Risler JL, Delorme MO, Delacroix H, Henaut A (1988) Amino acid substitutions in structurally related proteins. J Mol Biol 204: 1019-1029.
13. Gonnet GH, Cohen MA, Benner SA (1992) Exhaustive matching of the entire protein sequence database. Science 256: 1443-1445.
14. Henikoff S, Henikoff JG (1992) Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci USA 89: 10915-10919.
15. Overington J, Donnelly D, Johnson MS, Sali A, Blundell TL (1992) Environment-specific amino acid substitution tables: Tertiary templates and prediction of protein folds. Prot Sci 1: 216-226.
16. Jones DT, Taylor WR, Thornton JM (1992) The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 8: 275-282.
17. Kann M, Qian B, Goldstein RA (2000) Optimization of a new score function for the detection of remote homologs. Proteins 41: 498-503.
18. Ng PC, Henikoff JG, Henikoff S (2000) PHAT: a transmembrane-specific substitution matrix. Bioinformatics 16: 760-766.
19. Müller T, Rahmann S, Rehmsmeier M (2001) Non-symmetric score matrices and the detection of homologous transmembrane proteins. Bioinformatics 17, Suppl. 1: S182-S189.
20. Goonesekere NC, Lee B (2008) Context-specific amino acid substitution matrices and their use in the detection of protein homologs. Proteins 71: 910-919.
21. States DJ, Gish W, Altschul SF (1991) Improved sensitivity of nucleic acid database searches using application-specific scoring matrices. Methods 3: 66-70.
22. Chiaromonte F, Yap VB, Miller W (2002) Scoring pairwise genomic sequence alignments. In: Altman R, Dunker AK, Hunter L, Lauderdale K, Klein TE, eds. Proc. Pacific Symp. Biocomput. Mountain View, CA: World Scientific. pp 115-126.
23. Karlin S, Altschul SF (1990) Methods for assessing the statistical significance of molecular sequence features by using general scoring schemes. Proc Natl Acad Sci USA 87: 2264-2268.
24. Dembo A, Karlin S, Zeitouni O (1994) Limit distribution of maximal nonaligned two-sequence segmental score. Ann Prob 22: 2022-2039.
25. Sankoff D (1975) Minimal mutation trees of sequences. SIAM J Appl Math 28: 35-42.
26. Sankoff D, Cedergren RJ (1983) Simultaneous comparison of three or more sequences related by a tree. In: Sankoff D, Kruskal JB, eds. Time Warps, String Edits and Macromolecules: The Theory and Practice of Sequence Comparison. Reading, MA: Addison-Wesley. pp 253-263.
27. Murata M, Richardson JS, Sussman JL (1985) Simultaneous comparison of three protein sequences. Proc Natl Acad Sci USA 82: 3073-3077.
28. Bacon DJ, Anderson WF (1986) Multiple sequence alignment. J Mol Biol 191: 153-161.
29. Schneider TD, Stormo GD, Gold L, Ehrenfeucht A (1986) Information content of binding sites on nucleotide sequences. J Mol Biol 188: 415-431.
30. Altschul SF (1991) Amino acid substitution matrices from an information theoretic perspective. J Mol Biol 219: 555-565.
31. Brown M, Hughey R, Krogh A, Mian IS, Sjölander K, et al. (1993) Using Dirichlet mixture priors to derive hidden Markov models for protein families. In: Hunter L, Searls D, Shavlik J, eds. Proc. First Int. Conf. on Intelligent System for Mol. Biol. Menlo Park, CA: AAAI Press. pp 47-55.
32. Sjölander K, Karplus K, Brown M, Hughey R, Krogh A, et al. (1996) Dirichlet mixtures: a method for improved detection of weak but significant protein sequence homology. Comput Appl Biosci 12: 327-345.
33. Grünwald PD (2007) The Minimum Description Length Principle. Cambridge, MA: MIT Press.
34. Xing EP, Karp RM (2004) MotifPrototyper: a Bayesian profile model for motif families. Proc Natl Acad Sci USA 101: 10523-10528.
35. Lunter G, Miklós I, Drummond A, Jensen JL, Hein J (2005) Bayesian coestimation of phylogeny and sequence alignment. BMC Bioinformatics 6: 83.
36. Bradley RK, Roberts A, Smoot M, Juvekar S, Do J, et al. (2009) Fast statistical alignment. PLoS Comput Biol 5: e1000392.
37. Satija R, Novák A, Miklós I, Lyngsø R, Hein J (2009) BigFoot: Bayesian alignment and phylogenetic footprinting with MCMC. BMC Evol Biol 9: 217.
38. Feng DF, Doolittle RF (1987) Progressive sequence alignment as a prerequisite to correct phylogenetic trees. J Mol Evol 25: 351-360.
39. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673-4680.
40. Pietrokovski S (1996) Searching databases of conserved sequence regions by aligning protein multiple-alignments. Nucleic Acids Res 24: 3836-3845.
41. Rychlewski L, Jaroszewski L, Li W, Godzik A (2000) Comparison of sequence profiles. strategies for structural predictions using sequence information. Protein Sci 9: 232-241.
42. Yona G, Levitt M (2002) Within the twilight zone: a sensitive profile-profile comparison tool based on information theory. J Mol Biol 315: 1257-1275.
43. Edgar RC, Sjölander K (2003) SATCHMO: sequence alignment and tree construction using hidden markov models. Bioinformatics 19: 1404-1411.
44. Panchenko AR (2003) Finding weak similarities between proteins by sequence profile comparison. Nucleic Acids Res 31: 683-689.
45. Sadreyev R, Grishin N (2003) COMPASS: a tool for comparison of multiple protein alignments with assessment of statistical significance. J Mol Biol 326: 317-336.
46. Edgar RC, Sjölander K (2004) A comparison of scoring functions for protein sequence profile alignment. Bioinformatics 20: 1301-1308.
47. Wang G, Dunbrack RL, Jr. (2004) Scoring profile-to-profile sequence alignments. Protein Sci 13: 1612-1626.
48. Söding J (2005) Protein homology detection by HMM-HMM comparison. Bioinformatics 21: 951-960.
49. MacKay DJC (2003) Information Theory, Inference, and Learning Algorithms. New York, NY: Cambridge University Press.
50. Altschul SF, Carroll RJ, Lipman DJ (1989) Weights for data related by a tree. J Mol Biol 207: 647-653.
51. Sibbald PR, Argos P (1990) Weighting aligned protein or nucleic acid sequences to correct for unequal representation. J Mol Biol 216: 813-818.
52. Sander C, Schneider R (1991) Database of homology-derived protein structures and the structural meaning of sequence alignment. Proteins 9: 56-68.
53. Vingron M, Sibbald PR (1993) Weighting in sequence space: a comparison of methods in terms of generalized sequences. Proc Natl Acad Sci USA 90: 8777-8781.
54. Gerstein M, Sonnhammer ELL, Chothia C (1994) Volume changes in protein evolution. Appendix: A method to weight protein sequences to correct for unequal representation. J Mol Biol 236: 1067-1078.
55. Henikoff S, Henikoff JG (1994) Position-based sequence weights. J Mol Biol 243: 574-578.
56. Thompson JD, Higgins DG, Gibson TJ (1994) Improved sensitivity of profile searches through the use of sequence weights and gap excision. Comput Appl Biosci 10: 19-29.
57. Eddy SR, Mitchison G, Durbin R (1995) Maximum discrimination hidden Markov models of sequence consensus. J Comput Biol 2: 9-23.
58. Gotoh O (1995) A weighting system and algorithm for aligning many phylogenetically related sequences. Comput Appl Biosci 11: 543-551.
59. Krogh A, Mitchison G (1995) Maximum entropy weighting of aligned sequences of protein or DNA. In: Rawlings C, Clark D, Altman R, Hunter L, Lengauer T, et al. (1995) Proc. Third Int. Conf. on Intelligent System for Mol. Biol. Menlo Park, CA: AAAI Press. pp 215-221.
60. Bailey TL, Gribskov M (1996) The megaprior heuristic for discovering protein sequence patterns. In: States D, Agarwal P, Gaasterland T, Hunter L, Smith R, eds. Proc. Fourth Int. Conf. on Intelligent System for Mol. Biol. pp 15-24.
61. Sunyaev SR, Eisenhaber F, Rodchenkov IV, Eisenhaber B, Tumanyan VG, et al. (1999) PSIC: profile extraction from sequence alignments with positionspecific counts of independent observations. Protein Eng 12: 387-394.
62. Brown DP, Krishnamurthy N, Sjö lander K (2007) Automated protein subfamily identification and classification. PLoS Comput Biol 3: e160.
63. Altschul SF, Gertz EM, Agarwala R, Schäffer AA, Yu YK (2009) PSI-BLAST pseudocounts and the minimum description length principle. Nucleic Acids Res 37: 815-824.
64. Yu YK, Wootton JC, Altschul SF (2003) The compositional adjustment of amino acid substitution matrices. Proc Natl Acad Sci USA 100: 15688-15693.
65. Yu YK, Altschul SF (2005) The construction of amino acid substitution matrices for the comparison of proteins with non-standard compositions. Bioinformatics 21: 902-911.
66. Jeffreys H (1946) An invariant form of the prior probability in estimation problems. Proc Royal Soc London Series A 186: 453-461.
67. Nishida K, Frith MC, Nakai K (2009) Pseudocounts for transcription factor binding sites. Nucleic Acids Res 37: 939-944.
68. Vingron M, Waterman MS (1994) Sequence alignment and penalty choice. Review of concepts, case studies and implications. J Mol Biol 235: 1-12.
69. Lawrence CE, Altschul SF, Boguski MS, Liu JS, Neuwald AF, et al. (1993) Detecting subtle sequence signals: A Gibbs sampling strategy for multiple alignment. Science 262: 208-214.
70. Cover TM, Thomas JA (1991) Elements of Information Theory. New York, NY: Wiley.
71. Sjölander K (1998) Phylogenetic inference in protein superfamilies: analysis of SH2 domains. In: Glasgow J, Littlejohn T, Major F, Lathrop R, Sankoff D, et al. (1998) Proc. Sixth Int. Conf. on Intelligent System for Mol. Biol. Menlo Park, CA: AAAI Press. pp 165-174.
72. Brown DP (2008) Efficient functional clustering of protein sequences using the Dirichlet process. Bioinformatics 24: 1765-1771.
73. Altschul SF (1989) Gap costs for multiple sequence alignment. J Theor Biol 138: 297-309.
74. Thorne JL, Kishino H, Felsenstein J (1991) An evolutionary model for maximum likelihood alignment of DNA sequences. J Mol Evol 33: 114-124.
75. Thorne JL, Kishino H, Felsenstein J (1992) Inching toward reality: an improved likelihood model of sequence evolution. J Mol Evol 34: 3-16.
76. Tanaka H, Ishikawa M, Asai K, Konagaya A (1993) Hidden Markov models and iterative aligners: study of their equivalence and possibilities. In: Hunter L, Searls D, Shavlik J, eds. Proc. First Int. Conf. on Intelligent System for Mol. Biol. Menlo Park, CA: AAAI Press. pp 395-401.
77. Baldi P, Chauvin Y, Hunkapiller T, McClure MA (1994) Hidden Markov models of biological primary sequence information. Proc Natl Acad Sci USA 91: 1059-1063.
78. Krogh A, Brown M, Mian IS, Sjölander K, Haussler D (1994) Hidden Markov models in computational biology. Applications to protein modeling. J Mol Biol 235: 1501-1531.
79. Durbin R, Eddy S, Krogh A, Mitchison G (1998) Biological sequence analysis. Probabilistic models of proteins and nucleic acids. Cambridge, England: Cambridge University Press.
80. Eddy SR (1998) Profile hidden Markov models. Bioinformatics 14: 755-763.
81. Karplus K, Barrett C, Hughey R (1998) Hidden Markov models for detecting remote protein homologies. Bioinformatics 14: 846-856.
82. Neuwald AF, Liu JS (2004) Gapped alignment of protein sequence motifs through Monte Carlo optimization of a hidden Markov model. BMC Bioinformatics 5: 157.
83. Gotoh O (1982) An improved algorithm for matching biological sequences. J Mol Biol 162: 705-708.
84. Fitch WM, Smith TF (1983) Optimal sequence alignments. Proc Natl Acad Sci USA 80: 1382-1386.
85. Altschul SF, Erickson BW (1986) Optimal sequence alignment using affine gap costs. Bull Math Biol 48: 603-616.
86. Waterman MS, Smith TF, Beyer WA (1976) Some biological sequence metrics. Adv Math 20: 367-387.
87. Miller W, Myers EW (1988) Sequence comparison with concave weighting functions. Bull Math Biol 50: 97-120.
88. Benner SA, Cohen MA, Gonnet GH (1993) Empirical and structural models for insertions and deletions in the divergent evolution of proteins. J Mol Biol 229: 1065-1082.
89. Goonesekere NC, Lee B (2004) Frequency of gaps observed in a structurally aligned protein pair database suggests a simple gap penalty function. Nucleic Acids Res 32: 2838-2843.
90. Ladurner AG, Fersht AR (1997) Glutamine, alanine or glycine repeats inserted into the loop of a protein have minimal effects on stability and folding rates. J Mol Biol 273: 330-337.
91. Scalley-Kim M, Minard P, Baker D (2003) Low free energy cost of very long loop insertions in proteins. Protein Sci 12: 197-206.
92. Erickson BW, Sellers PH (1983) Recognition of patterns in genetic sequences. In: Sankoff D, Kruskal JB, eds. Time Warps, String Edits and Macromolecules: The Theory and Practice of Sequence Comparison. Reading, MA: Addison- Wesley. pp 55-91.
93. Rocke E, Tompa M (1998) An algorithm for finding novel gapped motifs in dna sequences. In: RECOMB '98: Proceedings of the second annual international conference on Computational molecular biology. pp 228-233.
94. Wareham HT, Jiang T, Zhang X, Trendall CG (2000) Stochastic heuristic algorithms for target motif identification (extended abstract). Pac Symp Biocomput. pp 392-403.
95. Thompson JD, Plewniak F, Poch O (1999) BAliBASE: a benchmark alignment database for the evaluation of multiple alignment programs. Bioinformatics 15: 87-88.
96. Subramanian AR, Weyer-Menkhoff J, Kaufmann M, Morgenstern B (2005) DIALIGN-T: an improved algorithm for segment-based multiple sequence alignment. BMC Bioinformatics 6: 66.
97. Stoye J, Evers D, Meyer F (1998) Rose: generating sequence families. Bioinformatics 14: 157-163.
98. Subramanian AR, Kaufmann M, Morgenstern B (2008) DIALIGN-TX: greedy and progressive approaches for segment-based multiple sequence alignment. Algorithms Mol Biol 3: 6.
99. Papadopoulos JS, Agarwala R (2007) COBALT: constraint-based alignment tool for multiple protein sequences. Bioinformatics 23: 1073-1079.
100. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673-4680.
101. Pei J, Sadreyev R, Grishin NV (2003) PCMA: fast and accurate multiple sequence alignment based on profile consistency. Bioinformatics 19: 427-428.
102. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792-1797.
103. Edgar RC (2004) MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5: 113.
104. Do CB, Mahabhashyam MS, Brudno M, Batzoglou S (2005) ProbCons: Probabilistic consistency-based multiple sequence alignment. Genome Res 15: 330-340.
105. Sonnhammer EL, Eddy SR, Durbin R (1997) Pfam: a comprehensive database of protein domain families based on seed alignments. Proteins 28: 405-420.
106. Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R (1998) Pfam: multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res 26: 320-322.
107. Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, et al. (2008) The Pfam protein families database. Nucleic Acids Res 36: D281-288.
108. Yada T, Ishikawa M, Tanaka H, Asai K (1996) Extraction of hidden Markov model representations of signal patterns in DNA sequences. Pac Symp Biocomput. pp 686-696.
109. Won KJ, Prugel-Bennett A, Krogh A (2004) Training HMM structure with genetic algorithm for biological sequence analysis. Bioinformatics 20: 3613-3619.
110. Won KJ, Sandelin A, Marstrand TT, Krogh A (2008) Modeling promoter grammars with evolving hidden Markov models. Bioinformatics 24: 1669-1675.
111. Mott R (1999) Local sequence alignments with monotonic gap penalties. Bioinformatics 15: 455-462.
112. Balaji S, Babu MM, Iyer LM, Aravind L (2005) Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res 33: 3994-4006.
113. Magnani E, Sjölander K, Hake S (2004) From endonucleases to transcription factors: evolution of the AP2 DNA binding domain in plants. Plant Cell 16: 2265-2277.
114. Wuitschick JD, Lindstrom PR, Meyer AE, Karrer KM (2004) Homing endonucleases encoded by germ line-limited genes in Tetrahymena thermophila have APETELA2 DNA binding domains. Eukaryotic Cell 3: 685-694.
115. De Silva EK, Gehrke AR, Olszewski K, León I, Chahal JS, et al. (2008) Specific DNA-binding by apicomplexan AP2 transcription factors. Proc Natl Acad Sci USA 105: 8393-8398.
116. Yuda M, Iwanaga S, Shigenobu S, Mair GR, Janse CJ, et al. (2009) Identification of a transcription factor in the mosquito-invasive stage of malaria parasites. Mol Microbiol 71: 1402-1414.
117. Phuong TM, Do CB, Edgar RC, Batzoglou S (2006) Multiple alignment of protein sequences with repeats and rearrangements. Nucleic Acids Res 34: 5932-5942.
118. Raphael B, Zhi D, Tang H, Pevzner P (2004) A novel method for multiple alignment of sequences with repeated and shuffled elements. Genome Res 14: 2336-2346.
119. Wootton JC (1994) Non-globular domains in protein sequences: automated segmentation using complexity measures. Comput Chem 18: 269-285.
120. Schneider TD, Stephens RM (1990) Sequence logos: a new way to display consensus sequences. Nucleic Acids Res 18: 6097-6100.
121. Allen MD, Yamasaki K, Ohme-Takagi M, Tateno M, Suzuki M (1998) A novel mode of DNA recognition by a beta-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA. EMBO J 17: 5484-5496.
122. Lindner SE, De Silva EK, Keck JL, Llinás M (2010) Structural determinants of DNA binding by a P. falciparum ApiAP2 transcriptional regulator. J Mol Biol 395: 558-567.