High quality protein sequence alignment by combining structural profile prediction and profile alignment using SABER-TOOTH.
BACKGROUND
Protein alignments are an essential tool for many bioinformatics analyses. While sequence alignments are accurate for proteins of high sequence similarity, they become unreliable as they approach the so-called 'twilight zone' where sequence similarity gets indistinguishable from random. For such distant pairs, structure alignment is of much better quality. Nevertheless, sequence alignment is the only choice in the majority of cases where structural data is not available. This situation demands development of methods that extend the applicability of accurate sequence alignment to distantly related proteins.


RESULTS
We develop a sequence alignment method that combines the prediction of a structural profile based on the protein's sequence with the alignment of that profile using our recently published alignment tool SABERTOOTH. In particular, we predict the contact vector of protein structures using an artificial neural network based on position-specific scoring matrices generated by PSI-BLAST and align these predicted contact vectors. The resulting sequence alignments are assessed using two different tests: First, we assess the alignment quality by measuring the derived structural similarity for cases in which structures are available. In a second test, we quantify the ability of the significance score of the alignments to recognize structural and evolutionary relationships. As a benchmark we use a representative set of the SCOP (structural classification of proteins) database, with similarities ranging from closely related proteins at SCOP family level, to very distantly related proteins at SCOP fold level. Comparing these results with some prominent sequence alignment tools, we find that SABERTOOTH produces sequence alignments of better quality than those of Clustal W, T-Coffee, MUSCLE, and PSI-BLAST. HHpred, one of the most sophisticated and computationally expensive tools available, outperforms our alignment algorithm at family and superfamily levels, while the use of SABERTOOTH is advantageous for alignments at fold level. Our alignment scheme will profit from future improvements of structural profiles prediction.


CONCLUSIONS
We present the automatic sequence alignment tool SABERTOOTH that computes pairwise sequence alignments of very high quality. SABERTOOTH is especially advantageous when applied to alignments of remotely related proteins. The source code is available at http://www.fkp.tu-darmstadt.de/sabertooth_project/, free for academic users upon request.
Teichert
Florian
F
Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr, Darmstadt, Germany.
Minning
Jonas
J
Bastolla
Ugo
U
Porto
Markus
M
D017510:Protein Folding; D017434:Protein Structure, Tertiary; D011506:Proteins / Q000737:chemistry*; D016415:Sequence Alignment / Q000379:methods*; D020539:Sequence Analysis, Protein / Q000379:methods; D012984:Software
D011506:Proteins
J Mol Biol. 1999 Sep 17;292(2):195-202
10493868
J Mol Biol. 2000 Sep 8;302(1):205-17
10964570
Nucleic Acids Res. 2003 Jul 1;31(13):3370-4
12824330
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D189-92
14681391
Nucleic Acids Res. 2008 Jan;36(Database issue):D234-9
17921502
Nucleic Acids Res. 1997 Sep 1;25(17):3389-402
9254694
Nucleic Acids Res. 1994 Nov 11;22(22):4673-80
7984417
Nucleic Acids Res. 2004 Mar 19;32(5):1792-7
15034147
Bioinformatics. 2005 Apr 1;21(7):951-60
15531603
Trends Biochem Sci. 2002 Mar;27(3):161-4
11893514
BMC Bioinformatics. 2006 Mar 30;7:180
16573808
Proteins. 2007;69 Suppl 8:3-9
17918729
Bioinformatics. 2000 Sep;16(9):776-85
11108700
Proteins. 2005 Jan 1;58(1):22-30
15523667
Protein Sci. 2002 Nov;11(11):2606-21
12381844
Proteins. 2005 Jan 1;58(1):158-65
15523668
Nat Biotechnol. 2004 Aug;22(8):1035-6
15286655
Protein Sci. 1998 Nov;7(11):2469-71
9828015
Bioinformatics. 1999 Jan;15(1):87-8
10068696
BMC Evol Biol. 2006 May 31;6:43
16737532
BMC Bioinformatics. 2006 Sep 05;7:401
16952323
Nucleic Acids Res. 2005 Apr 22;33(7):2302-9
15849316
BMC Bioinformatics. 2007 Oct 31;8:425
17974011
Protein Sci. 1996 Jul;5(7):1325-38
8819165