Template modeling score explained
In bioinformatics, the template modeling score or TM-score is a measure of similarity between two protein structures. The TM-score is intended as a more accurate measure of the global similarity of full-length protein structures than the often used RMSD measure. The TM-score indicates the similarity between two structures by a score between
, where 1 indicates a perfect match between two structures (thus the higher the better).
[1] Generally scores below 0.20 corresponds to randomly chosen unrelated proteins whereas structures with a score higher than 0.5 assume roughly the same fold.
[2] A quantitative study
[3] shows that proteins of TM-score = 0.5 have a
posterior probability of 37% in the same
CATH topology family and of 13% in the same
SCOP fold family. The probabilities increase rapidly when TM-score > 0.5. The TM-score is designed to be independent of protein lengths.
The TM-score equation
TM-score between two protein structures (e.g., a template structure and a target structure) is defined by
TM-score=max\left[
| 1 |
1+\left( | di | \right)2 | d0(Ltarget) |
|
\right]
where
is the length of the amino acid sequence of the target protein,and
is the number of residues that appear in both the template and target structures.
is the distance between the
th pair of residues in the template and target structures, and
d0(Ltarget)=1.24\sqrt[3]{Ltarget-15}-1.8
is a distance scale that normalizes distances.
When comparing two protein structures that have the same residue order,
reads from the C-alpha order number of the structure files (i.e., Column 23-26 in
Protein Data Bank (file format)). When comparing two protein structures that have different sequences and/or different residue orders, a
structural alignment is usually performed first, and TM-score is then calculated on the commonly aligned residues from the structural alignment.
Other measures
An often used structural similarity measure is root-mean-square deviation (RMSD). Because RMSD
is calculated as an average of distance error (
) with equal weight over all residue pairs, a large local error on a few residue pairs can result in a quite large RMSD. On the other hand, by putting
in the denominator, TM-score naturally weights smaller distance errors more strongly than larger distance errors. Therefore, TM-score value is more sensitive to the global structural similarity rather than to the local structural errors, compared to RMSD. Another advantage of TM-score is the introduction of the scale
d0(Ltarget)=1.24\sqrt[3]{Ltarget-15}-1.8
which makes the magnitude of TM-score length-independent for random structure pairs, while RMSD and most other measures are length-dependent metrics.
The Global Distance Test (GDT) algorithm, and its GDT TS score to represent "total score", is another measure of similarity between two protein structures with known amino acid correspondences (e.g. identical amino acid sequences) but different tertiary structures.[4] GDT score has the same length-dependence issue as RMSD, because the average GDT score for random structure pairs has a power-law dependence on the protein size.[1]
See also
- RMSD - a different structure comparison measure
- GDT - a different structure comparison measure
- Longest continuous segment (LCS) - A different structure comparison measure
- Global distance calculation (GDC_sc, GDC_all) - Structure comparison measures that use full-model information (not just α-carbon) to assess similarity
- Local global alignment (LGA) - Protein structure alignment program and structure comparison measure
External links
Notes and References
- Zhang Y and Skolnick J . Scoring function for automated assessment of protein structure template quality . Proteins . 57 . 702–710 . 2004 . 15476259 . 10.1002/prot.20264 . 4. 7954787 .
- Zhang Y and Skolnick J . TM-align: a protein structure alignment algorithm based on the TM-score . Nucleic Acids Res . 33 . 7 . 2302–2309 . 2005 . 15849316 . 10.1093/nar/gki524 . 1084323.
- Xu J and Zhang Y . How significant is a protein structure similarity with TM-score = 0.5? . Bioinformatics . 26 . 7 . 889–895 . 2010 . 20164152 . 2913670 . 10.1093/bioinformatics/btq066.
- Zemla A . 2003 . LGA: A method for finding 3D similarities in protein structures . Nucleic Acids Research . 31 . 13 . 3370–3374 . 12824330 . 168977 . 10.1093/nar/gkg571.