title: Detection of conserved segments in proteins: iterative scanning of sequence databases with alignment blocks.
abstract: We describe an approach to analyzing protein sequence databases that, starting from a single uncharacterized sequence or group of related sequences, generates blocks of conserved segments. The procedure involves iterative database scans with an evolving position-dependent weight matrix constructed from a coevolving set of aligned conserved segments. For each iteration, the expected distribution of matrix scores under a random model is used to set a cutoff score for the inclusion of a segment in the next iteration. This cutoff may be calculated to allow the chance inclusion of either a fixed number or a fixed proportion of false positive segments. With sufficiently high cutoff scores, the procedure converged for all alignment blocks studied, with varying numbers of iterations required. Different methods for calculating weight matrices from alignment blocks were compared. The most effective of those tested was a logarithm-of-odds, Bayesian-based approach that used prior residue probabilities calculated from a mixture of Dirichlet distributions. The procedure described was used to detect novel conserved motifs of potential biological importance.
authors.affiliation: National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894.
authors.forename: R L
authors.initials: RL
authors.lastname: Tatusov
authors.forename: S F
authors.initials: SF
authors.lastname: Altschul
authors.forename: E V
authors.initials: EV
authors.lastname: Koonin
mesh_terms: D000595:Amino Acid Sequence; D001419:Bacteria / Q000201:enzymology / Q000235:genetics; D005075:Biological Evolution; D016384:Consensus Sequence; D017124:Conserved Sequence; D004264:DNA Topoisomerases, Type I / Q000737:chemistry / Q000235:genetics; D016208:Databases, Factual; D008962:Models, Theoretical; D008969:Molecular Sequence Data; D011506:Proteins / Q000737:chemistry* / Q000235:genetics*; D012441:Saccharomyces cerevisiae / Q000201:enzymology / Q000235:genetics; D013223:Statistics as Topic
chemical_list: D011506:Proteins; D004264:DNA Topoisomerases, Type I