title: Evolutionary profiles from the QR factorization of multiple sequence alignments.
abstract: We present an algorithm to generate complete evolutionary profiles that represent the topology of the molecular phylogenetic tree of the homologous group. The method, based on the multidimensional QR factorization of numerically encoded multiple sequence alignments, removes redundancy from the alignments and orders the protein sequences by increasing linear dependence, resulting in the identification of a minimal basis set of sequences that spans the evolutionary space of the homologous group of proteins. We observe a general trend that these smaller, more evolutionarily balanced profiles have comparable and, in many cases, better performance in database searches than conventional profiles containing hundreds of sequences, constructed in an iterative and computationally intensive procedure. For more diverse families or superfamilies, with sequence identity <30%, structural alignments, based purely on the geometry of the protein structures, provide better alignments than pure sequence-based methods. Merging the structure and sequence information allows the construction of accurate profiles for distantly related groups. These structure-based profiles outperformed other sequence-based methods for finding distant homologs and were used to identify a putative class II cysteinyl-tRNA synthetase (CysRS) in several archaea that eluded previous annotation studies. Phylogenetic analysis showed the putative class II CysRSs to be a monophyletic group and homology modeling revealed a constellation of active site residues similar to that in the known class I CysRS.
authors.affiliation: Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
authors.forename: Anurag
authors.initials: A
authors.lastname: Sethi
authors.forename: Patrick
authors.initials: P
authors.lastname: O'Donoghue
authors.forename: Zaida
authors.initials: Z
authors.lastname: Luthey-Schulten
mesh_terms: D000465:Algorithms; D000595:Amino Acid Sequence; D000619:Aminohydrolases / Q000235:genetics; D001105:Archaea / Q000235:genetics; D019143:Evolution, Molecular; D008969:Molecular Sequence Data; D005810:Multigene Family; D016415:Sequence Alignment
chemical_list: C087227:imidazole glycerol phosphate synthase; D000619:Aminohydrolases
references.citation: Protein Sci. 1995 Feb;4(2):268-74
references.pmid: 7757015
references.citation: FEBS Lett. 2002 Mar 6;514(1):34-6
references.pmid: 11904177
references.citation: Proteins. 1997 Jul;28(3):405-20
references.pmid: 9223186
references.citation: Structure. 1997 Aug 15;5(8):1093-108
references.pmid: 9309224
references.citation: Science. 1997 Nov 7;278(5340):1119-22
references.pmid: 9353192
references.citation: J Mol Biol. 2005 Feb 25;346(3):875-94
references.pmid: 15713469
references.citation: Nat Struct Biol. 2002 Apr;9(4):257-62
references.pmid: 11887185
references.citation: Nucleic Acids Res. 2003 Jan 1;31(1):365-70
references.pmid: 12520024
references.citation: J Mol Biol. 2003 Apr 11;327(5):911-7
references.pmid: 12662918
references.citation: Microbiol Mol Biol Rev. 2003 Dec;67(4):550-73
references.pmid: 14665676
references.citation: J Bacteriol. 2004 Jan;186(1):8-14
references.pmid: 14679218
references.citation: Nucleic Acids Res. 2004 Jan 1;32(Database issue):D189-92
references.pmid: 14681391
references.citation: J Mol Biol. 1995 Apr 7;247(4):536-40
references.pmid: 7723011
references.citation: Biochimie. 1999 Nov;81(11):1037-9
references.pmid: 10575359
references.citation: Nucleic Acids Res. 2000 Mar 15;28(6):1481-8
references.pmid: 10684945
references.citation: Microbiol Mol Biol Rev. 2000 Mar;64(1):202-36
references.pmid: 10704480
references.citation: Comput Chem. 2001 Jul;25(4):329-39
references.pmid: 11459349
references.citation: J Struct Biol. 2001 May-Jun;134(2-3):257-68
references.pmid: 11551184
references.citation: Proc Natl Acad Sci U S A. 2001 Dec 4;98(25):14292-7
references.pmid: 11717392
references.citation: Science. 1996 Aug 23;273(5278):1058-73
references.pmid: 8688087