|SimShiftDB; local conformational restraints derived from chemical shift similarity searches on a large synthetic database.|
|Jump to Full Text|
|PMID: 19224375 Owner: NLM Status: MEDLINE|
|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.|
|Simon W Ginzinger; Murray Coles|
Related Documents :
|20472545 - Identifying duplicate content using statistically improbable phrases.
12870245 - An investigation of the suitability of the efqm excellence model for a pharmacy departm...
21068445 - Comparative performance of supertree algorithms in large data sets using the soapberry ...
20064025 - Memetic algorithms for continuous optimisation based on local search chains.
9278725 - Coevolutionary chase on exploiter-victim systems with polygenic characters.
25209335 - Cranial morphological diversity of early, middle, and late holocene brazilian groups: i...
|Type: Journal Article Date: 2009-02-18|
|Title: Journal of biomolecular NMR Volume: 43 ISSN: 1573-5001 ISO Abbreviation: J. Biomol. NMR Publication Date: 2009 Mar|
|Created Date: 2009-02-24 Completed Date: 2009-05-08 Revised Date: 2013-06-02|
Medline Journal Info:
|Nlm Unique ID: 9110829 Medline TA: J Biomol NMR Country: Netherlands|
|Languages: eng Pagination: 179-85 Citation Subset: IM|
|Department of Molecular Biology, Division of Bioinformatics, Center of Applied Molecular Engineering, University of Salzburg, Hellbrunnerstr. 34/3.OG, 5020, Salzburg, Austria. email@example.com|
|APA/MLA Format Download EndNote Download BibTex|
Nuclear Magnetic Resonance, Biomolecular*
Proteins / chemistry*
Sequence Homology, Amino Acid
Structural Homology, Protein
Journal ID (nlm-ta): J Biomol NMR
Publisher: Springer Netherlands, Dordrecht
? Springer Science+Business Media B.V. 2009
Received Day: 8 Month: 9 Year: 2008
Accepted Day: 7 Month: 1 Year: 2009
Electronic publication date: Day: 18 Month: 2 Year: 2009
pmc-release publication date: Day: 18 Month: 2 Year: 2009
Print publication date: Month: 3 Year: 2009
Volume: 43 Issue: 3
First Page: 179 Last Page: 185
PubMed Id: 19224375
Publisher Id: 9301
|SimShiftDB; local conformational restraints derived from chemical shift similarity searches on a large synthetic database|
|Simon W. Ginzinger1||
1Department of Molecular Biology, Division of Bioinformatics, Center of Applied Molecular Engineering, University of Salzburg, Hellbrunnerstr. 34/3.OG, 5020 Salzburg, Austria
2Department of Protein Evolution, Max-Planck-Institute for Developmental Biology, Spemannstrasse. 35, 72076 T?bingen, Germany
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. This restricts current databases to less than a few hundred proteins. 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. Experimental query shifts were compared to those back-calculated from the target. 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. 2007b). Given a target sequence and shifts, 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. 2004). 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 (S?ding et al. 2005), as a standard for purely sequenced based methods. We show that SimShiftDB can significantly increase the amount of information that can be derived from chemical shifts. We combine SimShiftDB with our CheckShift (Ginzinger et al. 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, a pairwise alignment is calculated. The alignment algorithm works in two steps:
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.?1 shows a depiction of a set of local similarities. 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.
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?3 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).
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. This is accomplished by using the AEROSPACI score (Chandonia et al. 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: 1H?, 1HN, 15N, 13C?, 13C?, and 13C?. 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: 13C??>?13C??>?1HN?>?13C??>?1H??>?15N. 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?3 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%. Figure?2 presents the results of this analysis. 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.
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 the inverse effect is observed. 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.
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 (S?ding 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. HHpred (S?ding et al. 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?3 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. We then compared all SimShiftDB predictions to the respective HHsearch predictions. 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?.
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.
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?3. We then compared the quality of each torsion angle prediction to the corresponding TALOS prediction. The results are presented using the same notation as in the previous section. 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.
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 we have compared SimShiftDB to a sequence-based method. 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.
The functional differences between SimShiftDB and TALOS are best illustrated by an example. 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?3, which correspond to the region G17-F40 of Ph1500C and consist of three-stranded ?-meander linked by two tight turns (Fig.?7). 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.
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 97?98% 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 2?3% 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 1HN, which is not generally regarded as containing much structural information. Perhaps this is due to some complementarities of the information from 1HN 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.
Most of the work by Simon W. Ginzinger was done while being a research fellow at the Department of Bioinformatics, Ludwig-Maximilians-Universit?t Munich (LMU), Amalienstrasse 17, 80333 Munich, Germany.
We thank Prof. Horst Kessler and his group at the Technical University of Munich for access to data on Ph1500C and other test cases and Johannes S?ding for helpful discussions and the integration of the SimShiftDB server in the MPI Bioinformatics Toolkit. Additionally the first author would like to thank Robert Konrat and Manfred J. Sippl for first introducing him to the problem of searching similarities in chemical shift sequences.
Availability SimShiftDB is available via a web server (http://simshiftdb.services.came.sbg.ac.at), see Fig.?8 for a sample screenshot. This server also provides a variety of functions for analyzing the results of a SimShiftDB Search interactively. Additionally, SimShiftDB will be included in the MPI Bioinformatics Toolkit (http://toolkit.tuebingen.mpg.de).
|Altschul FS,Gish W,Miller W,Myers EW,Lipman DJ. Basic local alignment search toolJ Mol BiolYear: 19902154034102231712|
|Battey JND,Kopp J,Bordoli L,Read RJ,Clarke ND,Schwede T. Automated server predictions in CASP7ProteinsYear: 2007698688210.1002/prot.2176117894354|
|Berjanskii M, Neal S, Wishart D (2006) PREDITOR: a web server for predicting protein torsion angle restraints. Nucleic Acids Res. 34 (Web Server issue): W63|
|Chandonia J-M,Hon G,Walker NS,Conte LL,Koehl P,Levitt M,Brenner SE. The ASTRAL compendium in 2004Nucleic Acids ResYear: 200432Database issueD189D19210.1093/nar/gkh03414681391|
|Cornilescu G,Delaglio F,Bax A. Protein backbone angle restraints from searching a database for chemical shift and sequence homologyJ Biomol NMRYear: 19991328930210.1023/A:100839240574010212987|
|Ginzinger SW (2008) Bioinformatics methods for NMR chemical shift data. PhD thesis, Ludwig-Maximilians Universit?t M?nchen URL: http://edoc.ub.uni-muenchen.de/8077/1/Ginzinger_Simon_Wolfgang.pdf|
|Ginzinger SW,Fischer J. SimShift: identifying structural similarities from NMR chemical shiftsBioinformaticsYear: 20062246046510.1093/bioinformatics/bti80516317071|
|Ginzinger SW,Gerick F,Coles M,Heun V. Checkshift: automatic correction of inconsistent chemical shift referencingJ Biomol NMRYear: 20073922322710.1007/s10858-007-9191-517899394|
|Ginzinger SW, Gr?upl T, Heun V (2007b) SimShiftDB: Chemical-shift-based homology modeling. In: Proceedings of the First Conference on Bioinformatics Research and Development, Lecture Notes in Bioinformatics, vol 4414, pp 357?370|
|Karlin S,Altschul SF. Applications and statistics for multiple high-scoring segments in molecular sequencesProc Natl Acad SciYear: 1993905873587710.1073/pnas.90.12.58738390686|
|Neal S,Nip AM,Zhang H,Wishart DS. Rapid and accurate calculation of protein 1H, 13C and 15N chemical shiftsJ Biomol NMRYear: 20032621524010.1023/A:102381293028812766419|
|Neal S,Berjanskii M,Zhang H,Wishart D. Accurate prediction of protein torsion angles using chemical shifts and sequence homologyMagn Reson ChemYear: 200644S158S16710.1002/mrc.183216823900|
|Seavey B,Farr E,Westler W,Markley J. A relational database for sequence specific protein NMR dataJ Biomol NMRYear: 1991121723610.1007/BF018755161841696|
|S?ding J. Protein homology detection by HMM-HMM comparisonBioinformaticsYear: 20052195196010.1093/bioinformatics/bti12515531603|
|S?ding J,Biegert A,Lupas AN. The HHpred interactive server for protein homology detection and structure predictionNucleic Acids ResYear: 200533Web Server issueW244W24810.1093/nar/gki40815980461|
Keywords: Keywords Chemical shift, Homology modeling, Alignment, Database search.
Previous Document: Relaxation-optimised Hartmann-Hahn transfer using a specifically Tailored MOCCA-XY16 mixing sequence...
Next Document: Species as explanatory hypotheses: refinements and implications.