Comprehensive assessment of automatic structural alignment against a manual standard, the scop classification of proteins
Tóm tắt
We apply a simple method for aligning protein sequences on the basis of a 3D structure, on a large scale, to the proteins in the scop classification of fold families. This allows us to assess, understand, and improve our automatic method against an objective, manually derived standard, a type of comprehensive evaluation that has not yet been possible for other structural alignment algorithms. Our basic approach directly matches the backbones of two structures, using repeated cycles of dynamic programming and least‐squares fitting to determine an alignment minimizing coordinate difference. Because of simplicity, our method can be readily modified to take into account additional features of protein structure such as the orientation of side chains or the location‐dependent cost of opening a gap. Our basic method, augmented by such modifications, can find reasonable alignments for all but 1.5% of the known structural similarities in scop, i.e., all but 32 of the 2,107 superfamily pairs. We discuss the specific protein structural features that make these 32 pairs so difficult to align and show how our procedure effectively partitions the relationships in scop into different categories, depending on what aspects of protein structure are involved (e.g., depending on whether or not consideration of side‐chain orientation is necessary for proper alignment). We also show how our pairwise alignment procedure can be extended to generate a multiple alignment for a group of related structures. We have compared these alignments in detail with corresponding manual ones culled from the literature. We find good agreement (to within 95% for the core regions), and detailed comparison highlights how particular protein structural features (such as certain strands) are problematical to align, giving somewhat ambiguous results. With these improvements and systematic tests, our procedure should be useful for the development of scop and the future classification of protein folds. Supplementary material is available at http://bioinfo.mbb.yale.edu/align.
Từ khóa
Tài liệu tham khảo
BrennerS ChothiaC HubbardT.1997. Assessing sequence comparison methods.Proc Natl Acad Sci USA. Forthcoming.
CohenGH.1997. ALIGN: A program to superimpose protein coordinates accounting for insertions and deletions.J Appl Crystallogr. Forthcoming.
Doolittle RF, 1987, Of Urfs and Orfs
Gerstein M, 1995, A structurally invariant core for the globins, CABIOS, 11, 633
Gerstein M, 1996, Proc Fourth Int Conf on Intell Sys Mol Biol, 59
GersteinM LevittM.1997. A structural census of the current population of protein sequences.Proc Natl Acad Sci USA. Forthcoming.
Godzik A, 1994, Flexible algorithm for direct multiple alignment of protein structures and sequences, CABIOS, 10, 587
Gribskov M, 1992, Sequence analysis primer
Holm L, 1994, The FSSP database of structurally aligned protein fold families, Nucleic Acid Res, 22, 3600
LevittM GersteinM.1998. A unified statistical framework for sequence comparison and structure comparison.Proc Natl Acad Sci USA. In press.
Overington JP, 1993, Molecular recognition in protein families: A database of aligned three‐dimensional structures of related proteins, Biochem Soc Transact, 3, 597, 10.1042/bst0210597
Rossmann MG, 1975, A comparison of the heme binding pocket in globins and cytochrome b5, J Biol Chem, 250, 7525, 10.1016/S0021-9258(19)40974-5