BioSysBio:abstracts/2007/Sanne Abeln

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
(Linking Folds)
(Linking Folds)
Line 58: Line 58:
A link between two superfamilies (I and J) is established when <math> f(I,J) > 0.1</math>, which is calculated as:
A link between two superfamilies (I and J) is established when <math> f(I,J) > 0.1</math>, which is calculated as:
-
<math> f(I,J) = \frac{Sim(I,J)}{ min\left \{ Sim(A-I,I),Sim(A-J,J))\right \} \mbox{ if }\ I \neq  J </math>  
+
<math> f(I,J) = \frac{Sim(I,J)}{ min( Sim(A-I,I),Sim(A-J,J)) \mbox{ if }\ I \neq  J </math>  
Here <math> Sim(A,B) </math> is the number of shared fragments between two set of domains (e.g. superfamilies), and A is the set of all domains. In this studies we do not consider selfsimilarity of superfamilies.  
Here <math> Sim(A,B) </math> is the number of shared fragments between two set of domains (e.g. superfamilies), and A is the set of all domains. In this studies we do not consider selfsimilarity of superfamilies.  

Revision as of 06:17, 29 September 2006

Linking evolution of protein structures through fragments

Author(s): Sanne Abeln, Charlotte M. Deane
Affiliations: University of Oxford
Contact:email: abeln@stats.ox.ac.uk
Keywords: 'protein structure' 'evolution' 'fragments' 'completed genomes'


Summary

Here we use a strucutural fragment library to investigate evolutionary links between protein folds. We show that 'older' folds have relatively more such links than 'younger' folds.


Motivation

At present there is no universal understanding 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 modelling. It has been proven difficult however to show the evolutionary relevance of such links between topologies based on geometrical measures. Here we use our a previously determined age measure for protein folds or superfamilies [1] to investigate the effect of structural fragments on protein structure evolution .

Results

  • When comparing the number of links (=degree), young and old superfamilies make


  • Figure 2 shows that the difference in links becomes stronger when we consider larger fragment lengths.
  • show figure with links between set of old folds and new folds
  • need brief discussion of evolutionary model
Figure 1: Density chart of the distribution of links between superfamilies for "old" folds and "yound" folds based on fragment length of 30. This histogram also shows that the distribution of links per superfamily is not a normally distributed
Figure 1: Density chart of the distribution of links between superfamilies for "old" folds and "yound" folds based on fragment length of 30. This histogram also shows that the distribution of links per superfamily is not a normally distributed
Figure 2: Fragment length versus W score of Wilcoxon's signed-rank test. Wilcoxon singed-rank tests were performed on paired data: the number of links each superfamily has to the group of young superfamilies, and old superfamilies, normalised for the size of the age groups.  Each test shows that superfamilies have significantly fewer links to young superfamilies (p-values < 2.2e-16). Since the number of compared folds in each test set are identical, the W scores can be compared directly.
Figure 2: Fragment length versus W score of Wilcoxon's signed-rank test. Wilcoxon singed-rank tests were performed on paired data: the number of links each superfamily has to the group of young superfamilies, and old superfamilies, normalised for the size of the age groups. Each test shows that superfamilies have significantly fewer links to young superfamilies (p-values < 2.2e-16). Since the number of compared folds in each test set are identical, the W scores can be compared directly.


Methods

Fragments

The fragment library generated for this study, contains fragment-pairs of length 10,15,20 and 30, with a maximum allowed gaplength of 2,3,4,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 in order to avoid dependencies between fragments and age estimates, which are generated through fold recognition techniques using 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 prefilter used by MAMMOTH [2]. Each fragment pair with an alignment score above a threshold is then superimposed to create an RMSD score for the fragment pair.

Age estimates

Age estimates for protein folds or superfamilies are generated using fold recongnition 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 or fold, for more details see [1].

The age of a fold or superfamily is based on a score between [0.0,1.0] with 0.0 indicating a last common recent ancestor at the leafs (youngest), and 1.0 indicating present at the root of the species tree (oldest). Here an 'old' fold is defined as a fold with an age of 1.0, and a 'young' fold with an age < 0.5

Linking Folds

Since no consideration of secondary structure is taken into account, the amount of shared fragments needs to be normalised for the amount a fragment occurs in general. Friedberg and Godzik (2005) used a used a superfamily based normalisation to overcome this problem [3]. We use a similar approach, although the fragment-pairs in this study are based on structural similarity only, whereas riedberg 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:

Failed to parse (syntax error): f(I,J) = \frac{Sim(I,J)}{ min( Sim(A-I,I),Sim(A-J,J)) \mbox{ if }\ I \neq J


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 studies we do not consider selfsimilarity of superfamilies.



Conclusion

We show that younger folds have relatiely fewer shared fragments with other fold, than old protein fold. This might indicate that evolutionary links above superfamily or fold level could be established, through such shared fragments.

References

  1. Winstanley HF, Abeln S, and Deane CM. . pmid:15961490. PubMed HubMed [Winstanley-2005]
  2. Ortiz AR, Strauss CE, and Olmea O. . pmid:12381844. PubMed HubMed [Ortiz-2002]
  3. Friedberg I and Godzik A. . pmid:15980462. PubMed HubMed [Friedberg-2005]
All Medline abstracts: PubMed HubMed
Personal tools