Volume 1 Supplement 1
Linking evolution of protein structures through fragments
© Abeln and Deane; licensee BioMed Central Ltd. 2007
Published: 8 May 2007
At present there is no universal understanding of how proteins can change topology during evolution, and how such pathways can be determined in a systematic way. The ability to create links between fold topologies would have important consequences for structural classification, structure prediction and homology modeling. Several methods based on geometrical measures have been proposed to create links between topologies, e.g. [1, 2]. It has proven difficult, however, to show the evolutionary relevance of such links. Here we use our previously developped age measure for protein superfamilies  to investigate the relationship between structural fragments and protein structure evolution.
Results and discussion
The fragment library generated for this study, contains fragment-pairs of length 10, 15, 20 and 30, with a maximum allowed gap-lengths of 2, 3, 4 and 6 respectively. All fragments are based on pairwise comparisons between structural domain as defined by SCOP. The pairs are scored for similarity purely on structural grounds, using the coordinates of the c-alpha atoms. This is to avoid bias, based on sequence similarity. All possible pairwise fragments between two domains of the given lengths are first screened and aligned using a method similar to the pre-filter used by MAMMOTH . Each fragment pair with an alignment score above a threshold is then superimposed giving the c-alpha RMSD score for the fragment pair.
Age estimates for protein folds or superfamilies are generated using fold recognition of structural domains on a set of completed genomes. The occurrence patterns of such predictions, are analysed with a parsimony algorithm to estimate an age for a superfamily, for more details see . The age of a superfamily is based on a score between [0.0,1.0], with 1.0 indicating the superfamily was estimated to be present at the root of the species tree (oldest), and 0.0 estimating that the superfamily was created at the leaf level (youngest). Here an 'old' fold is defined as a fold with an age of 1.0, and a 'young' fold with an age < 0.5.
Some fragments might be over-represented (e.g. secondary structure is not considered) therefore the number of shared fragments needs to be normalised for the number of times a fragment occurs. Friedberg and Godzik (2005) used a superfamily based normalisation to overcome this problem . We use a similar approach, although the fragment-pairs in this study are based on structural similarity only. (whereas Friedberg and Godzik (2005) used a combination of sequence and structural similarity). A link between two superfamilies (I and J) is established when f(I, J) > 0.1, which is calculated as:
Here Sim(A, B) is the number of shared fragments between two set of domains (e.g. superfamilies), and A is the set of all domains. In this study we do not consider self-similarity of superfamilies.
We show that younger folds have relatively fewer shared fragments with other folds, than old protein folds. This may indicate that evolutionary links above superfamily or fold level could be established, through such shared fragments.
- Taylor WR: A 'periodic table' for protein structures. Nature. 2002, 416 (6881): 657-660. 10.1038/416657aPubMedView Article
- Friedberg I, Godzik A: Fragnostic: walking through protein structure space. Nucleic Acids Res. 2005, W249-W251. 33 Web Server
- Winstanley HF, Abeln S, Deane CM: How old is your fold?. Bioinformatics. 2005, 21 (Suppl 1): i449-i458. 10.1093/bioinformatics/bti1008PubMedView Article
- Ortiz AR, Strauss CE, Olmea O: MAMMOTH (matching molecular models obtained from theory): an automated method for model comparison. Protein Sci. 2002, 11 (11): 2606-2621. 10.1110/ps.0215902PubMedPubMed CentralView Article
This article is published under license to BioMed Central Ltd.