Towards a reliable objective function for multiple sequence alignments.
Multiple sequence alignment is a fundamental tool in a number of different domains in modern molecular biology, including functional and evolutionary studies of a protein family. Multiple alignments also play an essential role in the new integrated systems for genome annotation and analysis. Thus, the development of new multiple alignment scores and statistics is essential, in the spirit of the work dedicated to the evaluation of pairwise sequence alignments for database searching techniques. We present here norMD, a new objective scoring function for multiple sequence alignments. NorMD combines the advantages of the column-scoring techniques with the sensitivity of methods incorporating residue similarity scores. In addition, norMD incorporates ab initio sequence information, such as the number, length and similarity of the sequences to be aligned. The sensitivity and reliability of the norMD objective function is demonstrated using structural alignments in the SCOP and BAliBASE databases. The norMD scores are then applied to the multiple alignments of the complete sequences (MACS) detected by BlastP with E-value<10, for a set of 734 hypothetical proteins encoded by the Vibrio cholerae genome. Unrelated or badly aligned sequences were automatically removed from the MACS, leaving a high-quality multiple alignment which could be reliably exploited in a subsequent functional and/or structural annotation process. After removal of unreliable sequences, 176 (24 %) of the alignments contained at least one sequence with a functional annotation. 103 of these new matches were supported by significant hits to the Interpro domain and motif database.
Thompson
J D
JD
Laboratoire de Biologie et Génomique Structurales, Institut de Génétique et de Biologie Moléculaire et Cellulaire, (CNRS/INSERM/ULP), Illkirch Cedex, 67404, France.
Plewniak
F
F
Ripp
R
R
Thierry
J C
JC
Poch
O
O
D020816:Amino Acid Motifs; D000595:Amino Acid Sequence; D019843:Archaeal Proteins / Q000737:chemistry / Q000235:genetics; D001426:Bacterial Proteins / Q000737:chemistry* / Q000235:genetics / Q000378:metabolism*; D019295:Computational Biology / Q000379:methods*; D030541:Databases, Genetic; D005057:Eukaryotic Cells / Q000378:metabolism; D016680:Genome, Bacterial; D023281:Genomics / Q000379:methods; D008958:Models, Molecular; D008969:Molecular Sequence Data; D017434:Protein Structure, Tertiary; D015203:Reproducibility of Results; D012107:Research Design; D012680:Sensitivity and Specificity; D016415:Sequence Alignment / Q000379:methods*; D012984:Software; D014734:Vibrio cholerae / Q000737:chemistry* / Q000235:genetics / Q000472:pathogenicity
D019843:Archaeal Proteins; D001426:Bacterial Proteins