- Methodology article
- Open Access
Protein structural similarity search by Ramachandran codes
BMC Bioinformatics volume 8, Article number: 307 (2007)
Protein structural data has increased exponentially, such that fast and accurate tools are necessary to access structure similarity search. To improve the search speed, several methods have been designed to reduce three-dimensional protein structures to one-dimensional text strings that are then analyzed by traditional sequence alignment methods; however, the accuracy is usually sacrificed and the speed is still unable to match sequence similarity search tools. Here, we aimed to improve the linear encoding methodology and develop efficient search tools that can rapidly retrieve structural homologs from large protein databases.
We propose a new linear encoding method, SARST (S tructural similarity search A ided by R amachandran S equential T ransformation). SARST transforms protein structures into text strings through a Ramachandran map organized by nearest-neighbor clustering and uses a regenerative approach to produce substitution matrices. Then, classical sequence similarity search methods can be applied to the structural similarity search. Its accuracy is similar to Combinatorial Extension (CE) and works over 243,000 times faster, searching 34,000 proteins in 0.34 sec with a 3.2-GHz CPU. SARST provides statistically meaningful expectation values to assess the retrieved information. It has been implemented into a web service and a stand-alone Java program that is able to run on many different platforms.
As a database search method, SARST can rapidly distinguish high from low similarities and efficiently retrieve homologous structures. It demonstrates that the easily accessible linear encoding methodology has the potential to serve as a foundation for efficient protein structural similarity search tools. These search tools are supposed applicable to automated and high-throughput functional annotations or predictions for the ever increasing number of published protein structures in this post-genomic era.
The number of proteins found in structural databases has increased at such an unprecedented rate in recent years that achieving speed and accuracy simultaneously in protein structure similarity searches has become a formidable task. During evolution, three-dimensional (3D) structures are more conserved than amino acid sequences , and protein homologs that share highly conserved 3D structures may have unrecognizable sequence homology . Amino acid sequence search tools are fast; however, they have proven to be insufficient for detection of remote homology in structural databases . Structure alignment using delicate geometric algorithms is much more accurate than amino acid sequence comparisons, especially when the sequence homology is low . Many brilliant pairwise comparison tools have been created, such as Distance Alignment Tools, DALI , Combinatorial Extension, CE , and FAST Alignment and Search Tool, FAST , but still there is a demand for rapid similarity search tools because protein structure databases have outgrown the utility of pairwise-based searches.
Protein structures are not fully flexible; there are physical constraints on polypeptide conformation [7–11]. It is believed that the 3D structure can be reduced to a simpler form while maintaining the intrinsic structural information [12–24]. With the reduced data, a similarity search can become much easier and faster. A number of methods have been designed based on this idea to one-dimensionalize the 3D protein structure. For instance, Levine et al. (1984) compared 3D protein structures using the sequence of dihedral angles (ϕ, ψ) in a pairwise manner . Lesk (1998) modified Efimov's dissections on the Ramachandran diagram  and combined them with a reduced amino acid alphabet to linearly encode protein structures . Martin (2000) developed TOPSCAN, which uses topology strings to represent protein structures . However, most of these methods could not reach the accuracy comparable to conventional 3D structural comparison methods, and furthermore the implementation of some of them were limited because their methodology could not conveniently analyze fragments with missing residues . Consequently, linear encoding methods have long been considered to compromise accuracy for speed in protein structure comparisons . Nevertheless, there are advantages of the one-dimensional (1D) representation of protein structure, such as its easy applicability to multiple structural alignments , fold-recognition and genome annotations [19, 20]; besides, local backbone structure prediction has long been using linear encoding methodologies [24–27].
In recent years, linear encoding has been applied to large scale structural database searches. Methods like YAKUSA  and 3D-BLAST  can scan thousands of proteins thousands of times as rapid as CE with good performance in searching accuracy. In this post-genomic era when protein structural data increase exponentially, we believe that linear encoding methodology is capable of serving as the foundation for efficient protein structural similarity search tools and that there is still much valuable room left for its improvement.
Here we propose a linear encoding algorithm, Ramachandran sequential transformation, and introduce an efficient protein structural similarity search method, SARST (S tructural similarity search A ided by R amachandran S equential T ransformation) . SARST improves the linear encoding methodology and achieves higher search speed with less sacrifice of accuracy than previous methods. SARST converts 3D protein structures into two-dimensional Ramachandran maps  and further to 1D sequences by predefined assignments of regions to text letters (Ramachandran codes). Finally, conventional sequence similarity search methods can be applied to retrieve homologous proteins from structure databases. These approaches are illustrated in Figure 1.
SARST, using structurally meaningful Ramachandran strings, converts structural similarity search problems into sequence similarity searches. Besides inheriting the speed advantages of sequence-based methods, it provides a ranked hit list with similarity scores and statistically meaningful expectation values (E-values) to assess the reliability of the retrieved information.
Because SARST is aimed to be a database search method, information retrieval (IR) techniques, which have been widely used in the document, image, spatial database, and 3D protein structure database fields [30, 31], were used to evaluate it. SARST can detect remote homology and overcome structural incompleteness; we also report its performances on different structural classes (all alpha, all beta, alpha/beta and alpha+beta).
Algorithm – R amachandran s equential t ransformation (RST)
One thousand domains, each composed of a single polypeptide chain without missing residues, were randomly selected as the training set from the ASTRAL SCOP 1.67 40% identities (ID) subset [32–34] [see Additional file 1]. For every residue r n in the training set, the torsion angle phi (ϕ) formed by atoms of rn-1and r n , and psi (ψ) formed by r n and rn+1were calculated to convey the two-residue-long backbone conformation involving three consecutive residues. All the torsion angle pairs (ϕ, ψ) were mapped onto the Ramachandran (RM) plot, and their distribution was analyzed by counting the pairs contained in each 10° × 10° unit cell. There were 36 × 36 = 1,296 cells on the RM map, each with a known number of (ϕ, ψ) spots. These cells were clustered into 22 groups based on a parameter, RSAD (Root Square Angular Distance), defined to represent the "distances" among cells:
where -180° < Δϕ < 180° and -180° < Δψ < 180°. They represent the differences in ϕ and ψ angles between a pair of cells. S c is a scaling constant assisting in restricting the number of clustered groups.
The 1,296 cells were first ranked in descending order by their spot numbers and assigned as x1–x1296; then, each cell was assigned a representative angle pair (ϕ i , ψ i ), where ϕ i and ψ i stood for the central ϕ and ψ angles of x i , respectively. We defined n i as the spot number of x i and set up a distance matrix D by assigning each element, D ij , the RSAD between x i and x j . With this matrix, a nearest-neighbor clustering algorithm  was performed following the steps below:
Set x i = x1 and g = 1.
Assign x i to cluster C g and let it be the center of C g . Now, N g , the total number of spots in C g , is n i .
Find the nearest neighbor of x i . Let x j denote it.
For those x j with n j > 3:
If D ij is smaller than T D , the threshold of distance, N g + n j is smaller than T N , the threshold of spot number, and x j has not been clustered, then assign x j to C g .
If D ij <T D , N g + n j <T N , and x j has been clustered in C m , compare D ij and the distance between x j and the center of C m . If D ij is the smaller, reassign x j to C g .
For those x j with n j ≤ 3: if D ij < 0.5 × T D , assign x j to C g ; otherwise, simply assign it to C22.
Repeat steps (4)–(5) with the next nearest neighbor while D ij <T D .
If every cell has been clustered, then stop. Otherwise, find the cell possessing the most spots from those not yet clustered, let it be the new x i and set g = 2, then go to step (2).
In this procedure, we were able to adjust T D , T N and S c in formula (1) so as to cluster all the cells into 22 groups. Finally, each group was assigned an English letter. As shown in Figure 2, these assigned letters represented specific regions of the Ramachandran map, and were called "Ramachandran codes". According to these codes, the coordinates of a protein could be transformed into a text sequence in the order of residue serial numbers. If a chain contained missing or internal (ϕ, ψ) incalculable residues, those positions would be labelled as "X"s. The "sequence" generated by RST algorithm is structurally meaningful and very different from the amino acid sequence in nature; therefore, we call it Ramachandran sequence or Ramachandran string.
Building scoring matrices – a regenerative approach
Because RM codes differ from amino acids, suitable scoring matrices were created to perform RM sequence alignment searches. We developed a "regenerative approach", which started with a primitive (and trial) matrix and enabled us to produce scoring matrices generation after generation until the quality was acceptable:
(1) The densest cell of each RM code region has been assigned as the representative center during RST. Code regions with smaller RSAD would be spatially close on the map, which we believed could be given higher scores. Based on this concept, we first calculated the average RSAD () of all the representative centers and then built the "primitive scoring matrix I" (PSMI) using the formula:
for i = j, the scores were uniformly appointed as 10.
(2) Using PSMI, all-against-all RM sequence alignments were performed between the training set and the ASTRAL SCOP 1.67 50% ID subset by blast [36, 37]. FAST  was used as a filter to pick the pairs with alignment lengths larger than 50% to form a "primitive pair database".
(3) The algorithm of BLOSUM matrices  was implemented to this pair database to build primitive scoring matrix II (PSMII).
(4) After performing recursive all-against-all RM sequence alignment on the ASTRAL SCOP 1.67 50% ID subset using PSMII, the pairs with FAST alignment lengths larger than 80% were picked to form a "50% ID primitive pair database", which then generated PSMIII.
(5) With PSMIII, we repeated step (4) for ASTRAL SCOP 1.67 10, 20, 30, 40, 50, 70, 90, and 100% ID subsets and produced the SARST Scoring Matrix (SARSTSM) 10, 20, 30, 40, 50, 70, 90 and 100, respectively.
Optimization of the scoring matrix
A. Selection of the scoring matrix
In 2004, Aung collected 34,055 proteins covering about 90% of the ASTRAL SCOP 1.59 dataset to form a large target database, from which 108 query proteins were selected . To assess the applicability of the SARSTSM, we adopted this database as well as the following parameters commonly used in the information retrieval experiments:
A protein is regarded as "relevant" if it belongs to the same SCOP family classification as the query. These two parameters always had opposite tendencies; when attempts were made to ask for higher recalls with the same query, the precision would decrease. Because of this property, to judge the quality of IR experiments, the F-measure  was also used:F = 2 × Recall × Precision/(Recall + Precision)(5)
For every query protein, RM sequence searches were performed asking for 50 to 5,000 retrievals. We observed that SATSTSM20 outperformed other matrices in most of the cases [see Additional file 2].
B. Determination of the score scaling factor
BLOSUM matrices were generated using Henikoffs' formula :Score ij = f s × log2(q ij /e ij )(6)
where q ij is the observed and e ij is the expected probability of the occurrence for each i, j pair, and f s is a scaling factor. In their study, f s was appointed as 2. To optimize our scoring matrix, we selected the "20% ID pair database", which produced SARSTSM20, and then adjusted the scaling factor. The highest average F-measure (70.0%) after the retrieval of 500 proteins was determined with f s = 1.78. Accordingly, the matrix produced from the 20% ID pair database with f s = 1.78 was chosen as the standard scoring matrix for SARST (Table 1).
C. Determination of the X scores
In RM strings, code "X"s stood for missing residues, residues with incomplete backbone coordinates, or those providing insufficient information for the calculation of torsion angles. We supposed the X scores should be zero to exert a minimum effect on the accuracy of SARST. After the retrieval of 500 proteins with integer X score ranging from -3 to 3, the highest average F-measures (70.0%) was found at zero X score, in agreement with our supposition.
Evaluation of speed
Aung and Tan have used their large database to assess the performances of ProtDex2 and several other methods . We adopted their system and added our assessments to CE, FAST, YAKUSA , 3D-BLAST , BLAST  and SARST.
As shown in Table 2, when using a single 3.2-GHz CPU to search this large database (34,055 proteins), SARST registered an average running time of 0.34 second, almost as rapid as BLAST (0.30 second). The SARST running time is approximately 243,500, 18,400, 250, 105, 27 and over 2 times faster than CE, FAST, TOPSCAN, YAKUSA, 3D-BLAST and ProtDex2, respectively. In a multi-processor system, SARST is capable of distributing the calculation work. If dual 3.2-GHz hyperthreading processors were used, its average running time would be 0.16 second, about 5 times faster than ProtDex2 and 517,400 times faster than CE, which itself could not recruit multiple processors.
Evaluation of accuracy
The goal of SARST is to create an efficient database search method, information retrieval techniques that have been widely used in many database search and management fields were used to evaluate its accuracy. As shown in Figure 3, FAST was the most accurate method. SARST was the third most accurate, and had a higher accuracy when compared with YAKUSA, 3D-BLAST, TOPSCAN, BLAST, and ProtDex2, the former three of which are linear encoding methods.
Implementation: Performance using different structural classes
The 108 query proteins from Aung  were composed of SCOP entries belonging to the four major classes with an average family size of ~80. To examine the performances of SARST using different structural classes, a measure known as "fallout" was calculated after the retrieval of 80 proteins.
Fallout is a measure of the false positive rate; it is the probability of retrieving an irrelevant protein . Accordingly, an effective retrieval system will yield lower fallout. SARST generated lower fallout values when compared with recent linear encoding database search methods, YAKUSA and 3D-BLAST (Figure 4). The fallout rates of SARST are close to those of CE. Unlike BLAST, SARST and other structure-based algorithms, inclusive of linear encoding ones, had limited bias among the four structural classes.
Performance on incomplete structures
Linear encoding methods may have a weakness in transforming structural information for proteins with incomplete backbone coordinates or missing residues , which constituted about one-fifth of Aung's query proteins  and the entire ASTRAL SCOP dataset. The fallout values of SARST for query proteins with incomplete structures were compared with those of several other methods. As illustrated in the right part of Figure 4, SARST generated fewer false positives than other linear encoding methods and, more interestingly, CE (see Additional file 5 for further information). These data indicate that SARST has achieved improved performance on incomplete structures.
Effects of low sequence identities
To more precisely assess the efficiency of searching methods challenged with low sequence identities, a well-organized, non-redundant target database in which all query proteins have remote homologs is necessary. Since Aung's databases do not satisfy this purpose , we have generated a new database from the ASTRAL SCOP 1.69 dataset. Query proteins were selected from the ASTRAL 100% ID subset following these criteria: (1) belonging to the four major classes (all-α, all-β, α/β and α+β), (2) having a family size between 30 and 140 proteins, (3) sharing 10% or less sequence identities, (4) having at least two family members in the 10% ID subset, (5) having no missing residues or incomplete backbone coordinates, and (6) being able to reach 100% recall with all of the assessed tools. The proteins meeting criteria (1)–(5) were grouped according to the family classifications and ranked by their lengths; then, the median length from each family was chosen and tested for criterion (6). The 83 query proteins that met these criteria were subtracted from the original subset, yielding a target database of 24,337 proteins [see Additional file 3].
Using this new target database, IR experiments were performed to examine the effects of low sequence identities. Various identity subsets of the target database were searched. As shown in Figure 5, the precision of SARST decreased as it encountered proteins with low sequences identities but was not as negatively affected as the precision of BLAST, which decreased substantially when the sequence identities fell below 30%. In comparison with recent linear encoding methods like YAKUSA and 3D-BLAST, the precision of SARST was generally improved. It could be observed that, when tested with these non-redundant datasets, the accuracy of linear encoding methods was substantially lower than geometric algorithms like FAST and CE. We propose that this is because of the unavoidable loss of structural information in the process of 3D to 1D transformation, a phenomena discussed in the latter part of this article.
Reliability of searching results
For every retrieved structure, SARST provides not only a similarity score but an expectation value (E-value) to assess the significance of the score (see Discussion). A lower E-value correlates to a higher significance of the score. IR experiments were done to test the reliability of the E-value. As shown in Table 3, low E-values gave high precisions and low fallouts at both the superfamily and family levels. When E-values were below 10-10, for instance, the average precision was greater than 92% and the average fallout was lower than 0.04%. Thus, the rate of negative answers' being retrieved (as positives) was at most 0.04% by SARST in this particular database search. The reliability of the E-value lends greater significance to the structural, functional and even evolutionary relatedness information retrieved by SARST.
Normalization of SARST scores
According to our observations, larger proteins generated higher SARST scores, which did not always translate into smaller root mean square distances (RMSD) in the actual structural superimpositions. For this reason, in some situations, it would be better to normalize SARST scores. The same formula used to normalize FAST scores  was used to normalize SARST scores.
where S is the raw score and is the normalized one. M and N are the RM string lengths for two proteins. The precision increased when the hit list was rearranged in a descending order according to the value. For example, when SARST was run with Aung's database  under the recommended parameter settings (see Methods), the average precision increased from 84.1 to 86.3% after normalization.
The normalized scores were more sensitive to global structural similarities and thus more likely to retrieve SCOP family members, which were mainly clustered according to their overall structural similarities. However, local similarities, which measure the structural relatedness of substructures such as domains, are also important in many situations. Hence, the score normalization is adjustable by the user in the SARST web service.
Distantly related homologs retrieved by SARST – two examples
We selected two pairs of proteins to demonstrate how SARST could detect remote homology from a large structure database. These protein pairs were retrieved from the ASTRAL SCOP 1.69 dataset. The coordinates of the positive positions aligned by RM strings were extracted to perform superimposition before calculation of their minimum RMSD.
In the first example (Figure 6a), [SCOP:d1b3aa_] was the query protein and [SCOP:d1tvxa_] was one of its relevant retrievals. Both of these proteins are interleukin 8-like human chemokines. Their amino acid sequence identity was only 17.2% over a small alignment length (29 residues), whereas they were structurally very similar (minimum RMSD: 1.68 Å) with a much larger RM string alignment length (51 positions). This example indicates that SARST could successfully identify protein homologs sharing highly conserved 3D structures but low overall sequence homology (also seen in Figure 5). In the second example (Figure 6b), [SCOP:d1p3ca_], a Bacillus intermedius glutamyl endopeptidase, was a high score irrelevant retrieval of the query protein [SCOP:d1tpo__], trypsin from cow (Bos taurus). These two proteases exhibited only a 22% amino acid sequence identity. They had similar structures, and the catalytic triads were well aligned by SARST even though they belong to different families in the SCOP classification. There were several missing residues in the query protein, and there were major differences in length for some of the secondary structure elements (SSE), which would normally cause some failure to previous linear encoding methods . SARST successfully identified the structural and functional similarities using suitable "X" scores and gap penalties. (Note that SARST is a database search tool that aims to rapidly distinguish high from low similarities but not to give optimum pairwise structural alignments. The RM sequence alignments shown in Figure 6 demonstrate how SARST works on protein homologs sharing low amino acid sequence identity but does not guarantee the best way to superimpose protein structures.)
SARST  transformed structural information into text strings through the Ramachandran plot, and converted complex geometric superimposition problems into relatively simple sequence similarity search problems. Therefore, SARST compared favorably with conventional structure alignment search methods in terms of speed.
Because SARST uses a relatively simple scoring scheme and an optimized scoring matrix, it ran remarkably faster than previous linear encoding methods, like TOPSCAN, YAKUSA and 3D-BLAST. There are several structure similarity search tools that can run at impressively high speeds by searching databases stored with pre-analyzed structural information, such as ProtDex2. SARST was not only faster than ProtDex2 but also much more accurate.
Given reasonable thresholds and a single CPU, SARST could run more than two hundred thousand times faster than CE. In a multi-processor system, SARST can automatically distribute the calculation work and run even faster. For example, when we used 2 hyperthreading CPUs to run SARST, it executed over half a million times faster than CE, which could not support multiple processors unless it were run by multi-thread scripts programmed by the user. If SARST is run in a clustered environment, mpiBLAST [41–43] can be used as the search engine, and the increase in its running speed would be even more impressive.
Although the current version of SARST could not match FAST in terms of precision, its accuracy is only slightly lower than that of CE. Additionally, SARST could achieve much higher precision than common IR-based tools, such as ProtDex2. In comparison with other linear encoding methods like TOPSCAN, YAKUSA and 3D-BLAST, the accuracy of SARST is improved. FAST is reportedly more accurate than CE . SARST alone can serve as an efficient protein structural database search method; furthermore, a good web service can be developed if we combine SARST and FAST through a filter-and-refine strategy .
After examining the high score hits, we found that the irrelevant retrievals obtained by SARST were largely due to common substructures shared by proteins with different overall structures. In fact, because SARST uses BLAST (basic local alignment search tool) as the core search method, it is suitable for local structural alignment searches. We have suggested a method to normalize SARST scores, which could promote accuracy in the situation that users are more concerned about the overall structural similarities. Parameters such as gap penalties could be adjusted to achieve higher accuracy according to the needs of the user.
Alpha-helices are the most abundant form of regular secondary structure, and therefore the alpha helix-related codes inevitably have the highest occurrence in linear encoding methods . Because there is high probability that two alpha helix-related codes could be aligned by chance, one may expect that SARST, as well as many other linear encoding methods, would produce more false positives in searching structural homologs for all-alpha proteins. However, our results indicated that SARST and other recent linear encoding algorithms had fairly even performance for different structural classes (Figure 4) in comparison with traditional sequence alignment method. In the case of SARST, which outperforms other linear encoding methods, this improvement may result from two factors. (1) The substitution matrices were generated with Henikoffs' algorithm , which calculates similarity scores as the logarithm of the odds (lod) ratio of the observed versus expected probabilities of every code pair. The most abundantly occurring helix-related RM code pairs thus do not have high lod scores, preventing the overweighting of helical SSEs. (2) The introduction of T N , the threshold of group size, into Ramachandran sequential transformation, resulted in fine dissections of the helix-like region of the RM plot. There were nine helix-like RM codes, e.g. ABCDETKVP, enabling SARST to detect minor structural differences between two helical SSEs, reducing the false positive rate.
Missing residues can reduce the performance and accuracy in protein structural similarity searches, as reported in the SA-Search linear encoding system . SARST, however, uses "X" codes to represent missing residues and, given suitable X scores and optimum gap penalties, it suppresses the effects of structural incompleteness (Figure 4).
The precision of SARST was higher than TOPSCAN, YAKUSA and 3D-BLAST probably due to torsion angle properties. Torsion angles are too local to describe long-range residue-residue interactions and may be insufficient for the development of structure "alignment" methods; however, in developing "similarity search" methods through linear encoding, this regional property may have advantages. We hypothesize that linear encoding methods lose structural information in the transformation process, and thus the more information to be encoded the more likely it would be lost. As shown in Table 2 and Figure 3, SARST, encoding two-residue-long conformations by torsion angles, was faster and more accurate than YAKUSA, 3D-BAST and TOPSCAN. The YAKUSA algorithm uses alpha angle to convey four-residue-long interactions , 3D-BLAST uses alpha and kappa angles to describe five-residue-long backbone conformations , and, TOPSCAN considers even longer topological changes of secondary structural elements . These results may imply that the "encoding ratio" has an inverse relationship to the range of interactions and may play a major role in protein structural linearization methodologies.
We plan to modify the SARST algorithm to preserve more structural information in the transformation process and to achieve higher accuracy. Future versions of SARST may consult the hidden Markov model and the methods of the Camproux-derived structural alphabet 27 [44, 45]. About the alphabet size of Ramachandran codes, we had made many preliminary tests ranging from 13 to 23 prior to the choice of 23. It was found that, at least in this range, a larger alphabet size gave a higher precision; however, the performance of the current version of SARST is limited by its search engine such that a maximum of only 23 symbols can be used to compose RM sequences. We hypothesize that if more symbols could be used, the dissection of the Ramachandran plot would be finer, thereby increasing the accuracy.
Significance of SARST score
After the database search, SARST produces a list of hits ordered by a score measuring the structural similarity. Additionally, SARST provides the statistically meaningful E-value to assess the significance of the score (S). The E-value is the number of different alignments with scores equivalent to or better than S that are expected to occur by chance in a database search [36, 37, 46]. Thus, lower E-values yield more significant scores. This statistical significance is transferable to structural relatedness and functional classifications. For instance, as shown in Table 3, the retrievals with E-values lower than 10-25 almost all belong to the same family as the query protein (average precision > 99%; average fallout < 3.1 × 10-5). A score with an E-value lower than 10-13 can be regarded to have a superfamily-level significance as the average precision at the superfamily level is higher than 99% and the average fallout is lower than 7.8 × 10-5 under this E-value threshold. Thus, the chance that one retrieved protein with such a low E-value belongs to a different superfamily from the query is at most 7.8 × 10-5 in this particular database search.
Proteins in the same SCOP family have a clear evolutionary relationship and those sharing the same superfamily most likely have a common evolutionary origin . Automated procedures like Classification by Optimization (CO)  have been developed to link the Z-score, a measure of the statistical significance of the result relative to an alignment of random structures [4, 5, 49], to SCOP classifications and thereby predict protein evolutionary relationships, to which we hypothesize that the E-value provided by SARST is also transferable.
Expected applications of SARST
The primary advantage of SARST is its speed. SARST provides a high search speed without substantially compromising the accuracy. Identification of distantly related protein homologs from a large structural database may prove difficult for sequence search methods or be time-consuming when using conventional structural alignment methods; however, SARST can accomplish this task within one second. In addition, this methodology is easy to implement, and multiple parameters can be adjusted by users to meet their research preferences. Moreover, the stand-alone version of SARST is written in Java and can run on many different platforms, turning a personal computer into an efficient instrument for a protein structural similarity searches.
Because of its high efficiency and portability, we hypothesize that SARST will be useful in automated and high-throughput functional annotations or predictions of the rapidly increasing protein structures produced by structural genomics researches. Because SARST describes protein structures as 1D strings, it can work together with multiple sequence alignment tools such as CLUSTAL W [50, 51] to perform rapid structural, functional or evolutionary clustering of proteins. In addition, fold-recognition and backbone structure predication have used the one-dimensionalization of protein structures for years [19, 20, 24–27] and may be applicable fields for SARST.
We have introduced a new protein structure similarity search method, SARST (S tructural similarity search A ided by R amachandran S equential T ransformation), which transforms 3D protein structures into 1D strings through a clustered Ramachandran map . This technique uses a regenerative approach to produce improved substitution matrices and recruits classical sequence alignment search methods to perform structural similarity searches. As a hybrid, SARST combines the speed advantages of sequence-based methods and accuracy advantages of structural comparisons. Its precision is only slightly lower than CE, and SARST executes hundreds of thousand times faster, almost as rapid as BLAST. In addition, SARST provides E-values to assess the reliability of the retrieved information.
SARST can detect remote homology that escapes a typical amino acid sequence alignment search. Its performance among different structural classes is similar to that of CE, without the normal bias shown by BLAST. Compared with previous linear encoding methods, SARST suppresses the problems caused by structural incompleteness by utilizing "X" codes and major differences in SSEs between homologous structures by using suitable gap penalties: it also achieves higher search speed and precision.
The fact that most linear encoding methods could not match conventional structure alignment methods in accuracy indicates that linear encoding might not be the best solution to protein structural comparisons; however, SARST demonstrates that it still has the potential to develop efficient structural similarity search tools. Protein structural data is increasing exponentially; thus, we hypothesize that efficient, easily accessible and highly portable similarity search methods like SARST will be the basic tool for post-genomic era researches.
The operating system was linux (Fedora Core 4) and, PHP (v.5.0.4) and Java 2 (v.1.4.2) were used to develop programs. The blast method described by Altschul et al.[36, 37] was used as the SARST search engine. All structures presented in the figures were drawn using PyMol .
Optimization of the search engine parameters
Because blastall (v.2.2.13) was recruited as the search engine, its parameter settings would affect the performance of SARST. Based on our early experience, the query sequence filter must be disabled (parameter setting: -F F) to achieve better search results. In addition, three other parameters were optimized: word size (W), gap-opening penalty (G), and gap-extension penalty (E).
There were two W values (2 and 3) allowed by blastall. We had used the small database developed by Aung  to determine their effects. It was found that the word size had limited effects on the precision of SARST, but the speed of SARST running under W = 3 was 3.4 times as rapid as that under W = 2. To meet speed requirements, size 3 was adopted.
At setting W = 3, the effects of all allowed combinations of G and E values were analyzed after the retrieval of 500 proteins. As shown in the additional material [see Additional file 4], SARST yielded the highest IR quality when G = 9 and E = 2, and it ran fastest when G = 25 and E = 2. Therefore, these are the recommended settings for SARST.
Practical parameter settings for SARST
Generally speaking, a well-developed search tool, such as NCBI's BLAST , would offer many parameters freely adjustable by the user to satisfy individual research preferences; however, a set of default values should also be provided to meet the common needs of users and to ensure high performance. To determine the practical parameter settings of SARST, Aung's large database  was used to compare the precision and speed of SARST under various v (number of database sequences to show one-line descriptions) and e (expectation value, or E-value) thresholds. When the v threshold was 250, the average recall was over 80%; thus, higher values seemed unnecessary. When the E-value threshold was above 10-7, the average precision fell below 80%; thus, higher thresholds appeared impractical. As such, we suggest that the combination of v = 250 and e = 10-7 would satisfy common needs. Running under these settings, the average recall of SARST was 76.0% and the average precision was 84.1%.
Assessment of speed and precision
Among the tools assessed in this report, the stand-alone version of CE and FAST could only perform pairwise comparisons; hence, the database searches were achieved using numerous pairwise comparisons with script programs. However, to make fair assessments only the actual running times were considered. The calculation times for parsing and sorting the results were omitted. The time consumed in parsing the outcomes of BLAST, ProtDex2, YAKUSA, 3D-BLAST and SARST was also omitted.
CE was good at local alignment, and therefore its output might contain many pairwise alignments of polypeptide fragments . In such cases, the alignment with the greatest length was selected as the final result.
FAST was designed to align two single polypeptide chains . Because many SCOP domains were composed of multiple fragments from different chains, they would cause FAST to function improperly. Thus, before any PDB file of SCOP domains was entered, it had to be "unified" first – all the chain IDs were changed to "A" regardless of the original labels, and all the residues were re-numbered consecutively.
Chothia C, Lesk AM: The relation between the divergence of sequence and structure in proteins. Embo J. 1986, 5 (4): 823-826.
Murzin AG: How far divergent evolution goes in proteins. Curr Opin Struct Biol. 1998, 8 (3): 380-387. 10.1016/S0959-440X(98)80073-0.
Sauder JM, Arthur JW, Dunbrack RL: Large-scale comparison of protein sequence alignment algorithms with structure alignments. Proteins. 2000, 40 (1): 6-22. 10.1002/(SICI)1097-0134(20000701)40:1<6::AID-PROT30>3.0.CO;2-7.
Holm L, Sander C: Protein structure comparison by alignment of distance matrices. J Mol Biol. 1993, 233 (1): 123-138. 10.1006/jmbi.1993.1489.
Shindyalov IN, Bourne PE: Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng. 1998, 11 (9): 739-747. 10.1093/protein/11.9.739.
Zhu J, Weng Z: FAST: a novel protein structure alignment algorithm. Proteins. 2005, 58 (3): 618-627. 10.1002/prot.20331.
Chou KC, Carlacci L: Energetic approach to the folding of alpha/beta barrels. Proteins. 1991, 9 (4): 280-295. 10.1002/prot.340090406.
Lasters I, Wodak SJ, Alard P, van Cutsem E: Structural principles of parallel beta-barrels in proteins. Proc Natl Acad Sci U S A. 1988, 85 (10): 3338-3342. 10.1073/pnas.85.10.3338.
Lasters I, Wodak SJ, Pio F: The design of idealized alpha/beta-barrels: analysis of beta-sheet closure requirements. Proteins. 1990, 7 (3): 249-256. 10.1002/prot.340070306.
Murzin AG, Lesk AM, Chothia C: Principles determining the structure of beta-sheet barrels in proteins. I. A theoretical analysis. J Mol Biol. 1994, 236 (5): 1369-1381. 10.1016/0022-2836(94)90064-7.
Richardson JS, Richardson DC: Principles and patterns of protein conformation. Prediction of protein structure and the principles of protein conformations. 1989, New York , Plenum, 1-98.
Sheridan RP, Dixon JS, Venkataraghavan R, Kuntz ID, Scott KP: Amino acid composition and hydrophobicity patterns of protein domains correlate with their structures. Biopolymers. 1985, 24 (10): 1995-2023. 10.1002/bip.360241011.
Mitchell EM, Artymiuk PJ, Rice DW, Willett P: Use of techniques derived from graph theory to compare secondary structure motifs in proteins. J Mol Biol. 1990, 212 (1): 151-166. 10.1016/0022-2836(90)90312-A.
Levine M, Stuart D, Williams J: A method for the systematic comparison of the three-dimensional structures of proteins and some results. Acta Cryst. 1984, A40: 600-610.
Efimov AV: Standard structures in proteins. Prog Biophys Mol Biol. 1993, 60 (3): 201-239. 10.1016/0079-6107(93)90015-C.
Lesk AM: Application of sequence alignment methods to multiple structural alignment and superposition. Proceedings of Prague Stringology Club Workshop '98. 1998, Prague , 95-100.
Martin AC: The ups and downs of protein topology; rapid comparison of protein structure. Protein Eng. 2000, 13 (12): 829-837. 10.1093/protein/13.12.829.
Guyon F, Camproux AC, Hochez J, Tuffery P: SA-Search: a web tool for protein structure mining based on a Structural Alphabet. Nucleic Acids Res. 2004, 32 (Web Server issue): W545-8. 10.1093/nar/gkh467.
Shi J, Blundell TL, Mizuguchi K: FUGUE: sequence-structure homology recognition using environment-specific substitution tables and structure-dependent gap penalties. J Mol Biol. 2001, 310 (1): 243-257. 10.1006/jmbi.2001.4762.
Kelley LA, MacCallum RM, Sternberg MJ: Enhanced genome annotation using structural profiles in the program 3D-PSSM. J Mol Biol. 2000, 299 (2): 499-520. 10.1006/jmbi.2000.3741.
Carpentier M, Brouillet S, Pothier J: YAKUSA: a fast structural database scanning method. Proteins. 2005, 61 (1): 137-151. 10.1002/prot.20517.
Yang JM, Tung CH: Protein structure database search and evolutionary classification. Nucleic Acids Res. 2006, 34 (13): 3646-3659. 10.1093/nar/gkl395.
Tyagi M, Gowri VS, Srinivasan N, de Brevern AG, Offmann B: A substitution matrix for structural alphabet based on structural alignment of homologous proteins and its applications. Proteins. 2006, 65 (1): 32-39. 10.1002/prot.21087.
de Brevern AG, Benros C, Gautier R, Valadie H, Hazout S, Etchebest C: Local backbone structure prediction of proteins. In Silico Biol. 2004, 4 (3): 381-386.
de Brevern AG, Etchebest C, Hazout S: Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks. Proteins. 2000, 41 (3): 271-287. 10.1002/1097-0134(20001115)41:3<271::AID-PROT10>3.0.CO;2-Z.
Etchebest C, Benros C, Hazout S, de Brevern AG: A structural alphabet for local protein structures: improved prediction methods. Proteins. 2005, 59 (4): 810-827. 10.1002/prot.20458.
De Brevern AG, Etchebest C, Benros C, Hazout S: "Pinning strategy": a novel approach for predicting the backbone structure in terms of protein blocks from sequence. J Biosci. 2007, 32 (1): 51-70. 10.1007/s12038-007-0006-3.
Structural similarity search Aided by Ramachandran Sequential Transformation. [http://sarst.life.nthu.edu.tw/sarst]
Ramachandran GN, Sasisekharan V: Conformation of polypeptides and proteins. Adv Protein Chem. 1968, 23: 283-438.
Aung Z, Tan KL: Rapid 3D protein structure database searching using information retrieval techniques. Bioinformatics. 2004, 20 (7): 1045-1052. 10.1093/bioinformatics/bth036.
Bertino E, Ooi BC, Sacks-Davis R, Tan KL, Zobel J, Shidlovsky B, Catania B: Indexing techniques for advanced database systems. 1997, Kluwer Academic Publishers
Brenner SE, Koehl P, Levitt M: The ASTRAL compendium for protein structure and sequence analysis. Nucleic Acids Res. 2000, 28 (1): 254-256. 10.1093/nar/28.1.254.
Chandonia JM, Hon G, Walker NS, Lo Conte L, Koehl P, Levitt M, Brenner SE: The ASTRAL Compendium in 2004. Nucleic Acids Res. 2004, 32 (Database issue): D189-92. 10.1093/nar/gkh034.
Chandonia JM, Walker NS, Lo Conte L, Koehl P, Levitt M, Brenner SE: ASTRAL compendium enhancements. Nucleic Acids Res. 2002, 30 (1): 260-263. 10.1093/nar/30.1.260.
Jain AK, Dubes RC: Algorithms for clustering data. 1988, New Jersey , Prentice Hall
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci U S A. 1992, 89 (22): 10915-10919. 10.1073/pnas.89.22.10915.
Hripcsak G, Rothschild AS: Agreement, the f-measure, and reliability in information retrieval. J Am Med Inform Assoc. 2005, 12 (3): 296-298. 10.1197/jamia.M1733.
Salton G, McGill MJ: Retrieval evaluation. Introduction to modern information retrieval. 1983, New York , McGraw-Hill, 174-177.
Darling A, Carey L, Feng W: The design, implementation, and evaluation of mpiBLAST. 4th International Conference on Linux Clusters: The HPC Revolution 2003 in conjunction with the ClusterWorld Conference & Expo. 2003, San Jose, CA
Feng W: Green Destiny + mpiBLAST = Bioinfomagic. 10th International Conference on Parallel Computing (ParCo). 2003, Dresden, Germany
Lin H, Ma X, Chandramohan P, Geist A, Samatova N: Efficient data access for parallel BLAST. IEEE International Parallel & Distributed Processing Symposium. 2005, Denver, CO
Camproux AC, Gautier R, Tuffery P: A hidden markov model derived structural alphabet for proteins. J Mol Biol. 2004, 339 (3): 591-605. 10.1016/j.jmb.2004.04.005.
Camproux AC, Tuffery P, Chevrolat JP, Boisvieux JF, Hazout S: Hidden Markov model approach for identifying the modular framework of the protein backbone. Protein Eng. 1999, 12 (12): 1063-1073. 10.1093/protein/12.12.1063.
Altschul SF, Gish W: Local alignment statistics. Methods Enzymol. 1996, 266: 460-480.
Park J, Karplus K, Barrett C, Hughey R, Haussler D, Hubbard T, Chothia C: Sequence comparisons using multiple sequences detect three times as many remote homologues as pairwise methods. J Mol Biol. 1998, 284 (4): 1201-1210. 10.1006/jmbi.1998.2221.
Getz G, Vendruscolo M, Sachs D, Domany E: Automated assignment of SCOP and CATH protein structure classifications from FSSP scores. Proteins. 2002, 46 (4): 405-415. 10.1002/prot.1176.
Holm L, Sander C: The FSSP database of structurally aligned protein fold families. Nucleic Acids Res. 1994, 22 (17): 3600-3609.
Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996, 266: 383-402.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
DeLano WL: The PyMOL molecular graphics system. 2002, San Carlos, CA, USA , DeLano Scientific
Protein-protein BLAST. [http://www.ncbi.nlm.nih.gov/blast/Blast.cgi]
Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM: AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996, 8 (4): 477-486. 10.1007/BF00228148.
This work was supported by the National Science Council, Taiwan, R.O.C. (NSC grant numbers: 94-3112-B-007-004-Y, 95-2752-B-007-003-PAE). We thank Dr. Yi-Shin Chen, Department of Computer Science, for insightful discussions.
WCL performed the clustering of the Ramachandran map, carried out the information retrieval experiments and participated in the design of the study. PJH participated in the design of the study and the production of scoring matrices. CHC participated in the development of the web service. PCL has conceived and coordinated the study. All authors contributed to writing the manuscript.
Electronic supplementary material
Additional file 1: The training set. A list of the SCOP domains in the training set. (XLS 94 KB)
Additional file 2: Performances of SARST scoring matrices. The F-measure and running time for SARST using different scoring matrices. (XLS 16 KB)
Additional file 3: The query and target proteins. A list of the 83 query proteins and the large target database containing 24,337 proteins collected from the ASTRAL SCOP 1.69 100% identity subset. (XLS 3 MB)
Additional file 4: Effects of gap penalties. The average recall, precision, F-measure and running time for SARST under various gap penalties. (XLS 14 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Authors’ original file for figure 1
Authors’ original file for figure 2
Authors’ original file for figure 3
Authors’ original file for figure 4
Authors’ original file for figure 5
Authors’ original file for figure 6
Authors’ original file for figure 7
Authors’ original file for figure 8
About this article
Cite this article
Lo, WC., Huang, PJ., Chang, CH. et al. Protein structural similarity search by Ramachandran codes. BMC Bioinformatics 8, 307 (2007). https://doi.org/10.1186/1471-2105-8-307
- Query Protein
- Linear Encode
- Glutamyl Endopeptidase
- Combinatorial Extension
- Scop Family