chemical shifts are now routinely used as a source of local conformational restraints in the structure determination of proteins by nmr , due mostly to the widespread use of programs such as talos ( cornilescu et al .
1999 ) and shiftor / preditor ( neal et al . 2006 ; berjanskii et al . 2006 ) .
these programs share a common approach and output similar data ; both search a database that correlates local patterns of chemical shifts with local conformation , and both provide backbone dihedral angle restraints for individual residues .
this approach has been very successful , but has some limitations in the stringent criteria needed for selecting proteins or protein fragments to populate the database , i.e. only those with both highly reliable chemical shift and structural data can be included .
although this may seem adequate  for example , the talos database contains 186 proteins subdivided into over 24,000 tripeptide fragments ( http://spin.niddk.nih.gov/nmrpipe/talos/)the sequence / conformation search space is very large , and the database coverage is unevenly distributed . as a result , rare combinations of amino - acid and conformation may be under - represented in the database , leading to significant under - prediction and even to errors outside of the heavily populated regions of the ramachandran map .
we have adopted an alternative approach to extracting structural data from chemical shifts based on our simshiftdb algorithm .
the original simshift was designed to test for structural similarities between proteins in a pair wise manner using chemical shifts to supplement sequence data .
simshift showed improved ability to detect distant structural relationships when compared to state - of - the - art methods based on the sequence alone .
a natural further development of pairwise comparison was to adapt the simshift algorithm for database searching , resulting in simshiftdb ( ginzinger et al .
, simshiftdb provides a list of matching proteins in the database , scored by a measure of statistical significance . in effect
it searches a synthetic chemical shift database of 13,000 proteins based on the astral library ( chandonia et al .
the matching sequence can be of any length , and structurally similar regions can be found ranging from small , locally similar fragments up to full domains . in principle , any structural alignment method can also be used to make predictions of local conformation by extracting torsion angles from matching regions of the target proteins , and it is this implementation of the simshiftdb algorithm we present here .
we benchmark the program against talos as a standard for current methods and hhpred , a sequence search method based on hidden markov models ( sding et al .
we show that simshiftdb can significantly increase the amount of information that can be derived from chemical shifts .
2007a ) routine for standardizing chemical shift referencing to produce a pipeline for analysis of chemical shift data .
full details of the simshiftdb algorithm have been published previously ( see ginzinger et al .
2007b ) , but it is worthwhile giving a brief overview here ( see fig .  1 ) .
given a target protein , simshiftdb analyzes each possible pairing of the target protein with one of the template protein structures in the database .
here we use the proteins from the astral database ( version 1.71 , filtered for 95% sequence identity ) to build the simshiftdb template database . as simshiftdb
is based on the comparison of chemical shift data , the chemical shifts for all entries in the template database are back - calculated from the three - dimensional structure using shiftx ( neal et al . 2003 ) . for each combination of the target protein with a template from the database ,
top : an example for a set of local similarities for a target protein ( s ) and a template protein ( t ) .
bottom : combination of a subset of local similarities to yield a self - consistent final alignmentstep 1 : local similarities are found by looking for high scoring combinations of parts of the target protein sequence ( s ) with parts of the template protein sequence ( t ) .
for example , block b in the figure shows that the chemical shifts of the target protein sequence from index xmin to xmax are similar to the chemical shifts of the template protein sequence from index ymin to ymax .
the similarity is calculated by summing the pairwise scores of the residues in the similar region in analogy to a pairwise sequence alignment .
the pairwise similarity scores are given by the so - called chemical shift substitution matrices , which give a score for each combination of two residues with associated chemical shifts ( for more details see ginzinger 2008).step 2 : the set of local similarities from step 1 is taken as an input for step 2 , where the most significant combination of blocks is identified , according to a statistical model of alignment scores ( karlin and altschul 1993 ) .
additionally , two blocks have to fulfill two constraints for their combination to be considered : an explanation of the simshift algorithm .
top : an example for a set of local similarities for a target protein ( s ) and a template protein ( t ) .
bottom : combination of a subset of local similarities to yield a self - consistent final alignment step 1 : local similarities are found by looking for high scoring combinations of parts of the target protein sequence ( s ) with parts of the template protein sequence ( t ) . fig . 
for example , block b in the figure shows that the chemical shifts of the target protein sequence from index xmin to xmax are similar to the chemical shifts of the template protein sequence from index ymin to ymax .
the similarity is calculated by summing the pairwise scores of the residues in the similar region in analogy to a pairwise sequence alignment .
the pairwise similarity scores are given by the so - called chemical shift substitution matrices , which give a score for each combination of two residues with associated chemical shifts ( for more details see ginzinger 2008 ) .
step 2 : the set of local similarities from step 1 is taken as an input for step 2 , where the most significant combination of blocks is identified , according to a statistical model of alignment scores ( karlin and altschul 1993 ) .
additionally , two blocks have to fulfill two constraints for their combination to be considered : blocks may not overlap in the target or in the template protein ; this would otherwise result in an ambiguous alignment.as the three - dimensional structure of the template protein is known , we further require that the euclidean distance between the end of the first block and the beginning of the second block may be bridged ( according to chemical restraints ) by the relevant sequence of amino acids in the target protein .
blocks may not overlap in the target or in the template protein ; this would otherwise result in an ambiguous alignment . as the three - dimensional structure of the template protein
is known , we further require that the euclidean distance between the end of the first block and the beginning of the second block may be bridged ( according to chemical restraints ) by the relevant sequence of amino acids in the target protein . finally , we calculate an e - value for the optimal combination of blocks .
this e - value represents the number of alignments of equal or better quality , which are expected to occur by chance , given the distribution of the amino acids with associated chemical shifts in the target protein and the template database .
additionally , the e - value takes the size of the template database into account . according to the following evaluation , an e - value of < 10 guarantees a high quality alignment .
to test the performance of simshiftdb , a benchmark set has to be defined for which both chemical shifts and the three - dimensional structure of the protein are available .
the bmrb ( seavey et al . 1991 ) is the main public repository for chemical shift data .
however , there is no consistent mapping to the structural databases , making it difficult to relate structural with chemical shift information reliably .
therefore a mapping between bmrb and astral is calculated based on amino acid sequence similarity .
every entry in the benchmark set has to fulfill the following constraints : a 100% sequence match to an astral entry.at least 100 residues with associated chemical shifts ( to exclude very short protein fragments ; e.g. single helices ) .
at least 100 residues with associated chemical shifts ( to exclude very short protein fragments ; e.g. single helices ) . to identify protein structures corresponding to the respective bmrb entries , a blast - search ( altschul et al . 1990 ) against the sequences from the astral database is conducted for each bmrb entry . if the full bmrb sequence can be matched without gaps against an astral sequence , the corresponding astral structure is assigned to the bmrb entry .
as some entries in bmrb match more than one sequence in astral , one representative structure has to be chosen .
2004 ) provided for each astral entry , thereby selecting the structure with the best resolution . through this procedure
a benchmark set containing 144 entries was derived . when calculating the similarity score for two residues , simshiftdb is restricted to at most three chemical shifts from the following list : h , h ,
n , c , c , and c. thus it is important to select a combination of shifts to extract maximum information , and a priority for replacing missing shifts . to identify the most successful strategy , we tested all possible priorities for the six atom types , resulting in 6 ! 
the most successful priority was : c  >  c  >  h  >  c  >  h  >  n. this is the default priority , and is used in the following analysis . to evaluate the prediction accuracy
, we applied the program to all entries in the benchmark set , using all proteins from the simshiftdb template database as potential templates .
subsequently , we used all alignments that achieved an e - value better than 10 to infer torsion angles for the target residue from the associated template residues . if residues of the target were mapped multiple times , we based the prediction on the highest scoring alignment .
it is important to evaluate the performance of simshiftdb as a function of sequence similarity ; therefore 9 evaluations were performed with the maximum allowed sequence similarity in the evaluated alignments set between 20 and 100% .
about 70% of all torsion angles predicted using simshiftdb have a high accuracy ( 15 error ) .
another 10 to 20% of the predicted angles have an error of less than 30. therefore simshiftdb yields accurate results in 85 to 90% of the evaluated predictions . through the use of the chemical shifts ,
this performance is nearly independent of the percentage of sequence identity in the respective alignments.fig . 
2evaluation of torsion angle predictions calculated using simshiftdb , separated according to the maximum allowed sequence identity between query and target .
the boxes show the percentage of dihedral angle predictions with an error of 30. the red part of each box shows the percentage of predictions with an error of 15. the brown line represents the percentage of residues from the test set for which a prediction exists evaluation of torsion angle predictions calculated using simshiftdb , separated according to the maximum allowed sequence identity between query and target .
the boxes show the percentage of dihedral angle predictions with an error of 30. the red part of each box shows the percentage of predictions with an error of 15. the brown line represents the percentage of residues from the test set for which a prediction exists also of interest is the performance of simshiftdb on different types of secondary structure .
figures  3 and 4 show the difference between the secondary structure content in all predictions versus that in high quality and erroneous predictions , respectively .
it can be seen that predictions for sheet ( for both high quality and erroneous predictions ) match well with the percentage observed in all predictions , and this match is largely independent of sequence similarity .
in contrast , the percentage of high quality predictions in helix increases with decreasing sequence similarity , whereas the corresponding percentage in coil regions decreases . for erroneous predictions
this seems logical , as the structures for coil regions are less reliable , and predictions are clearly harder to make than for secondary structure .
this test shows empirically that simshiftdb has no significant bias when comparing performance in predicting helix versus sheet .
the independence of sheet predictions from sequence similarity indicates that the chemical shift data is more diagnostic for sheet than for helix.fig . 
3difference between secondary structure content in high quality predictions versus the content for all predictions .
the bars show the difference in helix , sheet and coil content ( red , yellow and grey , respectively ) .
4difference between secondary structure content in erroneous predictions versus the content for all predictions using the same scheme as in fig . 
3 difference between secondary structure content in high quality predictions versus the content for all predictions .
the bars show the difference in helix , sheet and coil content ( red , yellow and grey , respectively ) .
the lines show the overall secondary structure content in all predictions difference between secondary structure content in erroneous predictions versus the content for all predictions using the same scheme as in fig .  3 to show empirically that simshiftdb uses the information in the chemical shift data to yield more sensitive alignments , especially in the case of low sequence similarity , we compare simshiftdb to hhsearch ( sding 2005 ) .
hhsearch , a sensitive search tool based on hidden markov models , calculates alignments between proteins using the primary sequence complemented by sequence - based predictions of secondary structure .
2005 ) , a protein structure prediction method based on hhsearch alignments , ranked second best in the casp7 ( battey et al .
additionally , it is freely available for download and gives the user the possibility to define arbitrary template databases . therefore it is perfectly suited to serve as a reference for purely sequence - based methods . to compare the performance of simshiftdb and hhsearch we used the benchmark set defined in the previous section . for both methods we ran each target protein against the simshiftdb template database and used alignments achieving an e - value better than 10 to predict torsion angles for the residues in the target protein .
if residues were mapped multiple times , the prediction was based on the highest scoring alignment .
the following notation is used for the presentation of the results :- a simshiftdb prediction is called better if it has an error of 30 and the corresponding hhsearch prediction has an error which is worse by more than 5.two predictions are called equal if both have an error of 30 and the difference between the errors is less than 5 , or both predictions have an error > 30.missing predictions are treated as predictions with an error > 30. a simshiftdb prediction is called better if it has an error of 30 and the corresponding hhsearch prediction has an error which is worse by more than 5. two predictions are called equal if both have an error of 30 and the difference between the errors is less than 5 , or both predictions have an error > 30. missing predictions are treated as predictions with an error > 30. figure  5 shows the results of this evaluation for alignments with maximal sequence identities ranging from 10 to 100% .
there is a clear trend for simshiftdb to outperform hhsearch as the sequence identity decreases , demonstrating that simshiftdb uses the structural information in chemical shifts to improve alignments.fig . 
5simshiftdb predictions compared to the respective hhsearch predictions , separated according to the maximum allowed sequence identity between query and target .
the red line shows the percentage of predictions where simshiftdb performs better , the yellow line show the percentage where simshiftdb is either better or gives a result of equal quality .
the brown line corresponds to the number of torsion angles predicted ( right axis ) simshiftdb predictions compared to the respective hhsearch predictions , separated according to the maximum allowed sequence identity between query and target .
the red line shows the percentage of predictions where simshiftdb performs better , the yellow line show the percentage where simshiftdb is either better or gives a result of equal quality .
the brown line corresponds to the number of torsion angles predicted ( right axis ) it is important to compare simshiftdb to the most prominent method for predicting torsion angles from chemical shifts , namely talos ( cornilescu et al .
again we used the benchmark set defined earlier , and made torsion angle predictions based on simshiftdb alignments that achieved an e - value better than 10 .
we then compared the quality of each torsion angle prediction to the corresponding talos prediction .
figure  6 shows that simshiftdb outperforms talos in at least 30% of all cases .
it should be noted that there is likely to be a significant bias towards talos in these results , as many members of the benchmark set will have been calculated using talos restraints.fig . 
6simshiftdb predictions compared to the respective talos predictions using the same notation as in fig . 
5 simshiftdb predictions compared to the respective talos predictions using the same notation as in fig .  5
to test the performance of simshiftdb , a benchmark set has to be defined for which both chemical shifts and the three - dimensional structure of the protein are available .
the bmrb ( seavey et al . 1991 ) is the main public repository for chemical shift data .
however , there is no consistent mapping to the structural databases , making it difficult to relate structural with chemical shift information reliably .
therefore a mapping between bmrb and astral is calculated based on amino acid sequence similarity .
every entry in the benchmark set has to fulfill the following constraints : a 100% sequence match to an astral entry.at least 100 residues with associated chemical shifts ( to exclude very short protein fragments ; e.g. single helices ) .
at least 100 residues with associated chemical shifts ( to exclude very short protein fragments ; e.g. single helices ) . to identify protein structures corresponding to the respective bmrb entries , a blast - search ( altschul et al . 1990 ) against the sequences from the astral database is conducted for each bmrb entry . if the full bmrb sequence can be matched without gaps against an astral sequence , the corresponding astral structure is assigned to the bmrb entry .
as some entries in bmrb match more than one sequence in astral , one representative structure has to be chosen .
2004 ) provided for each astral entry , thereby selecting the structure with the best resolution . through this procedure
when calculating the similarity score for two residues , simshiftdb is restricted to at most three chemical shifts from the following list : h , h , n , c , c , and c. thus it is important to select a combination of shifts to extract maximum information , and a priority for replacing missing shifts . to identify the most successful strategy , we tested all possible priorities for the six atom types , resulting in 6 !  =  720 evaluations .
the most successful priority was : c  >  c  >  h  >  c  >  h  >  n. this is the default priority , and is used in the following analysis . to evaluate the prediction accuracy
, we applied the program to all entries in the benchmark set , using all proteins from the simshiftdb template database as potential templates .
subsequently , we used all alignments that achieved an e - value better than 10 to infer torsion angles for the target residue from the associated template residues . if residues of the target were mapped multiple times , we based the prediction on the highest scoring alignment .
it is important to evaluate the performance of simshiftdb as a function of sequence similarity ; therefore 9 evaluations were performed with the maximum allowed sequence similarity in the evaluated alignments set between 20 and 100% .
about 70% of all torsion angles predicted using simshiftdb have a high accuracy ( 15 error ) .
another 10 to 20% of the predicted angles have an error of less than 30. therefore simshiftdb yields accurate results in 85 to 90% of the evaluated predictions . through the use of the chemical shifts ,
this performance is nearly independent of the percentage of sequence identity in the respective alignments.fig . 
2evaluation of torsion angle predictions calculated using simshiftdb , separated according to the maximum allowed sequence identity between query and target .
the boxes show the percentage of dihedral angle predictions with an error of 30. the red part of each box shows the percentage of predictions with an error of 15. the brown line represents the percentage of residues from the test set for which a prediction exists evaluation of torsion angle predictions calculated using simshiftdb , separated according to the maximum allowed sequence identity between query and target .
the boxes show the percentage of dihedral angle predictions with an error of 30. the red part of each box shows the percentage of predictions with an error of 15. the brown line represents the percentage of residues from the test set for which a prediction exists also of interest is the performance of simshiftdb on different types of secondary structure .
figures  3 and 4 show the difference between the secondary structure content in all predictions versus that in high quality and erroneous predictions , respectively .
it can be seen that predictions for sheet ( for both high quality and erroneous predictions ) match well with the percentage observed in all predictions , and this match is largely independent of sequence similarity .
in contrast , the percentage of high quality predictions in helix increases with decreasing sequence similarity , whereas the corresponding percentage in coil regions decreases . for erroneous predictions
this seems logical , as the structures for coil regions are less reliable , and predictions are clearly harder to make than for secondary structure .
this test shows empirically that simshiftdb has no significant bias when comparing performance in predicting helix versus sheet .
the independence of sheet predictions from sequence similarity indicates that the chemical shift data is more diagnostic for sheet than for helix.fig . 
3difference between secondary structure content in high quality predictions versus the content for all predictions .
the bars show the difference in helix , sheet and coil content ( red , yellow and grey , respectively ) .
4difference between secondary structure content in erroneous predictions versus the content for all predictions using the same scheme as in fig . 
3 difference between secondary structure content in high quality predictions versus the content for all predictions .
the bars show the difference in helix , sheet and coil content ( red , yellow and grey , respectively ) .
the lines show the overall secondary structure content in all predictions difference between secondary structure content in erroneous predictions versus the content for all predictions using the same scheme as in fig .  3
to show empirically that simshiftdb uses the information in the chemical shift data to yield more sensitive alignments , especially in the case of low sequence similarity , we compare simshiftdb to hhsearch ( sding 2005 ) .
hhsearch , a sensitive search tool based on hidden markov models , calculates alignments between proteins using the primary sequence complemented by sequence - based predictions of secondary structure .
2005 ) , a protein structure prediction method based on hhsearch alignments , ranked second best in the casp7 ( battey et al . 2007 ) experiment .
additionally , it is freely available for download and gives the user the possibility to define arbitrary template databases .
therefore it is perfectly suited to serve as a reference for purely sequence - based methods . to compare the performance of simshiftdb and hhsearch we used the benchmark set defined in the previous section .
for both methods we ran each target protein against the simshiftdb template database and used alignments achieving an e - value better than 10 to predict torsion angles for the residues in the target protein .
if residues were mapped multiple times , the prediction was based on the highest scoring alignment .
the following notation is used for the presentation of the results :- a simshiftdb prediction is called better if it has an error of 30 and the corresponding hhsearch prediction has an error which is worse by more than 5.two predictions are called equal if both have an error of 30 and the difference between the errors is less than 5 , or both predictions have an error > 30.missing predictions are treated as predictions with an error > 30. a simshiftdb prediction is called better if it has an error of 30 and the corresponding hhsearch prediction has an error which is worse by more than 5. two predictions are called equal if both have an error of 30 and the difference between the errors is less than 5 , or both predictions have an error > 30. missing predictions are treated as predictions with an error > 30. figure  5 shows the results of this evaluation for alignments with maximal sequence identities ranging from 10 to 100% .
there is a clear trend for simshiftdb to outperform hhsearch as the sequence identity decreases , demonstrating that simshiftdb uses the structural information in chemical shifts to improve alignments.fig . 
5simshiftdb predictions compared to the respective hhsearch predictions , separated according to the maximum allowed sequence identity between query and target .
the red line shows the percentage of predictions where simshiftdb performs better , the yellow line show the percentage where simshiftdb is either better or gives a result of equal quality .
the brown line corresponds to the number of torsion angles predicted ( right axis ) simshiftdb predictions compared to the respective hhsearch predictions , separated according to the maximum allowed sequence identity between query and target .
the red line shows the percentage of predictions where simshiftdb performs better , the yellow line show the percentage where simshiftdb is either better or gives a result of equal quality .
it is important to compare simshiftdb to the most prominent method for predicting torsion angles from chemical shifts , namely talos ( cornilescu et al . 1999 ) .
again we used the benchmark set defined earlier , and made torsion angle predictions based on simshiftdb alignments that achieved an e - value better than 10 .
we then compared the quality of each torsion angle prediction to the corresponding talos prediction .
figure  6 shows that simshiftdb outperforms talos in at least 30% of all cases .
it should be noted that there is likely to be a significant bias towards talos in these results , as many members of the benchmark set will have been calculated using talos restraints.fig . 
6simshiftdb predictions compared to the respective talos predictions using the same notation as in fig . 
5 simshiftdb predictions compared to the respective talos predictions using the same notation as in fig .  5
we have presented simshiftdb and shown that the program is able to sensitively extract structural information from chemical shift data .
this information is to a certain extent complementary to that from currently available tools . on one hand
simshiftdb shows its strength especially in cases of low sequence similarity , which underlines the advantage of including chemical shift information in the alignment algorithm . on the other hand
, we were able to show that one - third of the predictions by simshiftdb clearly have a higher quality than the corresponding talos predictions , and this is largely independent of sequence similarity .
the main advantage of simshiftdb is derived from its superior coverage of the search space , due to the large and quickly adaptable template database .
simshiftdb outperforms talos especially in those cases where talos finds no predictions classified as  good  according to its selection criteria .
simshiftdb and talos are therefore complementary , and can be used in parallel to increase the number of available predictions .
ph1500c is a 78-residue homo - hexameric domain currently under investigation in our laboratories , and was chosen because it shows no significant sequence similarity to proteins of known structure .
we used simshiftdb to search for templates matching ph1500c , using several different chemical shift priorities .
the search identifies several templates at e - values around 10 , which correspond to the region g17-f40 of ph1500c and consist of three - stranded -meander linked by two tight turns ( fig . 
the consensus of simshiftdb predictions shows the first turn to be a standard type i -turn , while the second is a less commonly observed 5-residue -turn .
an example template agrees very well with all structural data available in this region ( fig . 
7 ) , and suggests local conformational details such as sidechain positions and local hydrogen bonding networks.fig . 
the left panel shows a representative template structure found in a simshiftdb search for ph1500c ( 1p22 , residues a392-w415 ) .
residues shown in green are predicted by talos to be in the -region of ramachandran space and the residue in yellow in the + -region , while for residues in white there is no prediction .
only one talos prediction is made for the six residues outside of canonical secondary structure , i.e. the two turns and the -bulge in the first strand .
the right panel shows the structure calculated with all available nmr data , showing very good agreement to the template .
the coloring is as above , with the addition of blue residues for residues in the -region of ramachandran space .
a conserved sidechain hydrogen - bond acceptor ( t101 in ph1500c ) is shown in grey , highlighting the ability of simshiftdb searches to reveal fine structural detailsfig . 
8the simshiftdb web - interface at the center for applied molecular engineering an example simshiftdb search result .
the left panel shows a representative template structure found in a simshiftdb search for ph1500c ( 1p22 , residues a392-w415 ) .
residues shown in green are predicted by talos to be in the -region of ramachandran space and the residue in yellow in the + -region , while for residues in white there is no prediction .
only one talos prediction is made for the six residues outside of canonical secondary structure , i.e. the two turns and the -bulge in the first strand .
the right panel shows the structure calculated with all available nmr data , showing very good agreement to the template .
the coloring is as above , with the addition of blue residues for residues in the -region of ramachandran space .
a conserved sidechain hydrogen - bond acceptor ( t101 in ph1500c ) is shown in grey , highlighting the ability of simshiftdb searches to reveal fine structural details the simshiftdb web - interface at the center for applied molecular engineering this example highlights the major difference in the simshiftdb and talos approaches , i.e. the length of the template structures found by simshiftdb when compared to the tripeptides used to make talos predictions .
the second difference is the use of the e - value as a continuous measure of quality , rather than a discrete selection criterion based on a consensus of the ten best hits . in some cases there
may be only one or two templates found for any region of the protein , but low e - value scores can nevertheless allow predictions with high confidence .
we have established an accuracy of above 85% for simshiftdb predictions , based on our benchmark set of proteins .
this may at first glance compare poorly to talos , where an accuracy of 9798% is reported .
however , it must be considered that this value is based on single simshiftdb predictions , rather than the consensus of 10 predictions .
also , it is worth noting that talos is very accurate within secondary structure , and therefore the 23% of errors must be concentrated in the smaller fraction of other predictions . in our experience
, these errors often result from predictions made out of structural context ; e.g. for a residue in a -turn based on tripeptides from a helix .
the wider context provided by simshiftdb results should therefore add both to the confidence of its predictions and those from talos .
the optimum chemical shift priority found for simshiftdb searches is somewhat surprising in that it contains h , which is not generally regarded as containing much structural information .
perhaps this is due to some complementarities of the information from h and that from other shifts .
it is worth noting , though , that the difference between the best priorities is small , and it may be worth testing a range of priorities .
this is easily possible ; although the program searches a database of 13,000 protein structures , an average simshiftdb run takes only 30 seconds on a standard laptop ( intel t2500 , 2.0  ghz , 1  gb ram ) .
the different results are comparable using the calculated e - values , thereby enabling the user to select the most promising result .
 we present simshiftdb , a new program to extract conformational data from protein chemical shifts using structural alignments . 
 the alignments are obtained in searches of a large database containing 13,000 structures and corresponding back - calculated chemical shifts . 
 simshiftdb makes use of chemical shift data to provide accurate results even in the case of low sequence similarity , and with even coverage of the conformational search space . 
 we compare simshiftdb to hhsearch , a state - of - the - art sequence - based search tool , and to talos , the current standard tool for the task . 
 we show that for a significant fraction of the predicted similarities , simshiftdb outperforms the other two methods . 
 particularly , the high coverage afforded by the larger database often allows predictions to be made for residues not involved in canonical secondary structure , where talos predictions are both less frequent and more error prone . 
 thus simshiftdb can be seen as a complement to currently available methods . 
Introduction
Implementation
Results
The benchmark set
Evaluation of prediction accuracy
Comparison to HHsearch
Comparison to TALOS
Discussion