An optimized TOPS+ comparison method for enhanced TOPS models.  
Jump to Full Text  
MedLine Citation:

PMID: 20236520 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: Although methods based on highly abstract descriptions of protein structures, such as VAST and TOPS, can perform very fast protein structure comparison, the results can lack a high degree of biological significance. Previously we have discussed the basic mechanisms of our novel method for structure comparison based on our TOPS+ model (Topological descriptions of Protein Structures Enhanced with Ligand Information). In this paper we show how these results can be significantly improved using parameter optimization, and we call the resulting optimised TOPS+ method as advanced TOPS+ comparison method i.e. advTOPS+. RESULTS: We have developed a TOPS+ string model as an improvement to the TOPS 123 graph model by considering loops as secondary structure elements (SSEs) in addition to helices and strands, representing ligands as first class objects, and describing interactions between SSEs, and SSEs and ligands, by incoming and outgoing arcs, annotating SSEs with the interaction direction and type. Benchmarking results of an allagainstall pairwise comparison using a large dataset of 2,620 nonredundant structures from the PDB40 dataset 4 demonstrate the biological significance, in terms of SCOP classification at the superfamily level, of our TOPS+ comparison method. CONCLUSIONS: Our advanced TOPS+ comparison shows better performance on the PDB40 dataset 4 compared to our basic TOPS+ method, giving 90% accuracy for SCOP alpha+beta; a 6% increase in accuracy compared to the TOPS and basic TOPS+ methods. It also outperforms the TOPS, basic TOPS+ and SSAP comparison methods on the ChewKedem dataset 5, achieving 98% accuracy. SOFTWARE AVAILABILITY: The TOPS+ comparison server is available at http://balabio.dcs.gla.ac.uk/mallika/WebTOPS/. 
Authors:

Mallika Veeramalai; David Gilbert; Gabriel Valiente 
Publication Detail:

Type: Comparative Study; Journal Article; Research Support, NonU.S. Gov't Date: 20100317 
Journal Detail:

Title: BMC bioinformatics Volume: 11 ISSN: 14712105 ISO Abbreviation: BMC Bioinformatics Publication Date: 2010 
Date Detail:

Created Date: 20100422 Completed Date: 20101027 Revised Date: 20130715 
Medline Journal Info:

Nlm Unique ID: 100965194 Medline TA: BMC Bioinformatics Country: England 
Other Details:

Languages: eng Pagination: 138 Citation Subset: IM 
Affiliation:

Joint Center for Molecular Modeling, SanfordBurnham Medical Research Institute, La Jolla, CA 92037, USA. mallikav@burnham.org 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Algorithms Computational Biology / methods* Databases, Protein Ligands Models, Molecular Protein Conformation Protein Folding Proteins / chemistry* Software* 
Chemical  
Reg. No./Substance:

0/Ligands; 0/Proteins 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): BMC Bioinformatics ISSN: 14712105 Publisher: BioMed Central 
Article Information Download PDF Copyright ?2010 Veeramalai et al; licensee BioMed Central Ltd. openaccess: Received Day: 11 Month: 9 Year: 2009 Accepted Day: 17 Month: 3 Year: 2010 collection publication date: Year: 2010 Electronic publication date: Day: 17 Month: 3 Year: 2010 Volume: 11First Page: 138 Last Page: 138 ID: 2858036 Publisher Id: 1471210511138 PubMed Id: 20236520 DOI: 10.1186/1471210511138 
An optimized TOPS+ comparison method for enhanced TOPS models  
Mallika Veeramalai1  Email: mallikav@burnham.org 
David Gilbert2  Email: david.gilbert@brunel.ac.uk 
Gabriel Valiente3  Email: valiente@lsi.upc.edu 
1Joint Center for Molecular Modeling, SanfordBurnham Medical Research Institute, La Jolla, CA 92037, USA 

2School of Information Systems, Computing and Mathematics, Brunel University, Uxbridge, Middlesex UB8 3PH, UK 

3Algorithms, Bioinformatics, Complexity and Formal Methods Research Group, Technical University of Catalonia, E08034 Barcelona, Spain 
The structural genomics consortium [^{6}] aims to populate protein fold space using highthroughput experimental technologies. As a result the number of known structures in the Protein Data Bank (PDB) [^{7}] is increasing rapidly every year and currently holds 59,790 structures (August 25, 2009). This highlights the importance of the need for fast and reliable protein structure comparison methods. There are various methods which use detailed 3D structures for comparison; SSAP [^{8},^{9}] uses a double dynamic programming method that takes into account several different features of protein structure including phi/psi angles, accessibility and interresidue vectors to align two protein structures. Other approaches include STAMP [^{10}], DALI [^{11}] and the Combinatorial Extension method [^{12}]. On the other hand abstract level structural comparison methods are based on topological/vector models of secondary structure elements (SSEs) and their relationships. VAST is a vector based protein structure comparison method [^{13},^{14}]. GRATH [^{15}] is a graphbased algorithm that compares the axial vectors of alpha helices and beta strands of two proteins, together with the distances, angles and chirality between these vectors. It is based on a method by Grindley et al. [^{16}]. Earlier work by Koch et al. [^{17}] uses a graph method to find maximal common SSEs in a pair of proteins. TOPS is a graphbased method applied to the topological representation of the protein structures [^{3}]. Although these methods perform very fast protein structure comparison in most cases the results have significantly less biological interpretation due to the abstract nature of the protein model. Moreover, the functional annotation problem is made much more complex by the fact that the number of protein folds is limited while their range of functions is very diverse. For example, the current version of the SCOP database classified the (single) TIM barrels protein fold into 33 distinct functional superfamilies.
This motivated our research to develop a novel topological model for protein structures, enhanced with structural and biochemical features, such as ligand interaction information and amino acid sequence length of the secondary structures, in order to permit better, more biologically significant comparison methods. Previously, we have discussed the basic mechanisms of our novel TOPS+ comparison method for novel topological models. We compute the edit distance between two proteins based on TOPS+ strings elements using a dynamic programming approach. We have benchmarked our method with an allagainstall pairwise comparison using a large dataset of 2,620 nonredundant structures from the PDB40 and the results were validated using the standard SCOP superfamily classification numbers. We have also compared our method against other methods and showed that it is faster than SSAP, FATCAT, DALI and TOPS and that it has a comparable performance against TOPS [^{18}]. Recently we developed the TOPS++FATCAT system that exploits the TOPS+ strings comparison method to speedup the FATCAT protein structural alignment program for fast flexible structural alignment, while preserving the accuracy of the original FATCAT method [^{19}]. These promising results have facilitated the introduction of further constraints on ligandarc matching.
In this paper, we show how the above results can be significantly improved using parameter optimization at two stages of the TOPS+ method: (i) in the generation of the dynamic programming table and (ii) in the computation of the comparison score using a compression measure. The dynamic programming algorithm includes weight tables for matching TOPS+ strings elements, the match scores take into account not only the SSEType, orientation but in addition they include scores for total in/out/ligand arcs together with their arc types such as right and left chiralities, and parallel and antiparallel hydrogenbonds. This research work involved (a) generating the TOPS descriptions enhanced with in/out/ligand arc information for a large set of proteins; (b) designing the weight tables; (c) optimization of weights in the table; (d) designing a pairwise comparison metric based on a compression measure and optimizing different parameters to take into account the variability on both components of the topological and ligand interaction features. The optimization of our advanced TOPS+ comparison method was carried out on the PDB40 representative dataset. Furthermore, we assess the biological significance of our method against existing protein structure comparison methods based on cluster analysis and validation using an Fmeasure calculation [^{20},^{21}] on the ChewKedem dataset [^{4},^{5}].
Figure 1 shows the ROC curves and Table 1 gives AUC values for SCOP classes allalpha, allbeta, alpha/beta and alpha+beta on the PDB40 subset dataset obtained from the TOPS, TOPS+ and advanced TOPS+ (advTOPS+) methods. The results show that the advTOPS+ method is superior to TOPS+ on classes allbeta, alpha/beta and alpha+beta with 82%, 77% and 90% accuracy (see Table 1); while it gives similar results on the allalpha class with an accuracy of 82%. When we compared our advTOPS+ method with TOPS we have better performance on alpha+beta and allalpha classes with accuracy level increased by 6%; the result is comparable in the case of classes allbeta and alpha/beta protein domains. The TOPS method relies on arc information and in alpha rich proteins there are no hydrogenbond arcs and few chirality arcs, hence it performs poorly on the allalpha class. Our TOPS+ methods has a better performance compared to TOPS in allalpha class of proteins, because in our TOPS+ model we have included additional biochemical features such as loops, SSEligand interactions, and SSE segment length. Moreover in the allalpha class, most of the proteins have structuredependent ligand interactions such as DNAbinding proteins with (HelixLoopHelix = HLH, HelixTurnHelix = HTH) and metalbinding proteins (like HEMbinding proteins, etc.). Thus our method can recognize those proteins more efficiently compared to TOPS.
In alpha+beta class our advTOPS+ method has a 90% accuracy, which is superior when compared with both TOPS and our basic TOPS+ method, which have only 84% accuracy (see Table 1). Because these proteins are composed of segregated alpha and beta regions the structuredependent ligand interactions and additional chiral, hydrogen bonds are also present. Thus our parameter optimization can handle all arcs more efficiently.
On the other hand, the alpha/beta class of proteins contains mixed alpha and beta secondary structures; more importantly although the protein domains from these classes have ligand interactions, they may not be structuredependent ligand interactions. In these classes for most of the protein superfamilies the ligands have a tendency to bind the clefts or binding pocket which have appropriate physiochemical properties and the correct conformational geometry of the amino acids. Furthermore it is important to note that in our TOPS+ and advTOPS+ comparison methods we have considered only the total number of ligandarcs rather than the actual ligand property match, thus we have false positives in some SCOP classes. In the case of allbeta class proteins our advTOPS+ method has comparable performance against TOPS with 87% accuracy (see Table 1); in this class proteins contain a significant number of hydrogen bond and chiral arcs, and thus parameter optimization is performed more efficiently. From the Fmeasure statistical evaluation analysis (we used the same cutoff value of 0.35 for all three methods) we found that the advTOPS+ method appears to always do better than TOPS and TOPS+ except for the alpha/beta class of proteins (see Table 1).
The overall results show that our advTOPS+ method exhibits substantial improvement compared to basic TOPS+. It has better performance for allalpha and alpha+beta proteins compared to TOPS. On the other SCOP classes the performance is comparable with TOPS. Since our method considers only the total number of ligand arcs rather than the actual ligand property this leads to false positives to some extent. Our advTOPS+ method can efficiently recognize structuredependent ligand interactions appropriately in the case of DNAbinding proteins and metal binding proteins.
Our advanced TOPS+ comparison method outperforms all the other methods (TOPS, basic TOPS+ and SSAP) and groups 36 representative proteins from five fold families into biologically significant clusters (see Figure 2). In Figure 2 each protein domain is represented in the following format "domainName_foldFamilies". The tags of the "foldFamilies" represents protein structures from five different protein folds; where tb = TIM barrel, g = Globins, ab = alpha beta, b = all beta, and a = all alpha protein families. The clusters obtained from our advanced TOPS+ method show that all of the protein domains are grouped according to their structural fold and biological significance except d1ct9a1 (see explanations below). The basic TOPS+ method also groups most of the protein domains into correct fold families except for an all beta protein d1cdb__ and two protein domains d2hbg__, d1hlm__ from globins which are grouped together with all alpha proteins (see the supplementary material page for clusters from other methods and additional information at http://balabio.dcs.gla.ac.uk/mallika/WebTOPS/optTOPSplusresults.html). In comparison both the TOPS and SSAP methods produce clusters in which there are many wrongly grouped domains (see supplementary material page for clusters from other methods and additional information at http://balabio.dcs.gla.ac.uk/mallika/WebTOPS/optTOPSplusresults.html). The quantitative analysis using the Fmeasure calculation results in more than 98% accuracy for our advTOPS+ method, which is higher than SSAP (96%), TOPS (95%) and basic TOPS+ (93%) (see Table 2). When we compare our advTOPS+ with our basic TOPS+ method we have achieved a 5% improvement based on parameter optimization.
One interesting example from our cluster analysis (see Figure 2) is that of the Asparagine synthetase B, Cterminal domain d1ct9a1 (alphabeta fold) which is grouped by the advanced TOPS+ method with the TIM barrel domains (see Figure 3). Table 3 provides advanced TOPS+ comparison scores for d1ct9a1 (SSE_Ln = 53) against those protein domains belonging to alphabeta and TIMbarrel fold families, together with their SSEs length (SSE_Ln) and LCS pattern length (LCS PAT_Ln). Interestingly, d1ct9a1 has a smaller distance score against all of the TIMbarrel proteins compared to alphabeta protein domains; specifically, for the xylose isomerase protein d6xia__. When we closely checked the 3D structure of d1ct9a1, we found that it has alpha helices in its Nterminal domain and it contains half TIMbarrel like structures in its Cterminal domain (see Figure 3). This suggests that a subdomain of the d1ct9a1 exhibits a structural drift between the alphabeta fold and the TIMbarrel fold. The structural drift [^{22}] is a special case of gregariousness as described by [^{23}]. Since our TOPS+ method is based on an abstraction of protein structure in the form of SSEs and topological features without any geometrical properties it is able to perform matching at the fold level and to include most of the SSEs which are common to two protein structures.
In this paper we have reported the generation of TOPS+ and TOPS+ strings models for large datasets and have presented an improved TOPS+ comparison method using parameter optimization both for the computation of the dynamic programming table and the computation of the comparison score using a compression metric. Through our evaluation analysis we have showed that our advanced TOPS+ comparison method has a substantial improvement on all the SCOP classes compared to our basic TOPS+ method. Our advanced TOPS+ method has better performance compared to TOPS on alpha+beta and allalpha and is comparable on allbeta and alpha/beta. On the ChewKedem dataset our advanced TOPS+ comparison outperforms all the other methods.
This demonstrates that our TOPS+ and TOPS+ strings models can find more biologically significant results and has led to interesting new directions to incorporate ligandpattern discovery in TOPS+ comparison [^{24}]. Our method is faster than TOPS and SSAP because it has time complexity O(n^{2}), where n represents the number of SSEs in the protein domains. This research opens new doors to an exciting improvement to our TOPS+ models and advanced TOPS+ comparison method by the addition of features such as aminoacid sequences, biochemical properties of the proteinligand interaction at atomic level, and arc scores (both topological level and ligand level) for each SSE. Moreover we can improve the comparison process with additional statistical scoring values for each TOPS+ strings element match both at the micro (atomicdetails of proteinligand interaction information) and the macro level (abstract level).
Furthermore our novel TOPS+ models, TOPS+ strings and comparison approaches could be applicable to different problem areas such as RNA secondary structure comparison and prediction. Most of the drugdiscovery process starts with insilico chemical compound screening which is computationally expensive. Our TOPS+ comparison approach could be applied as an initial step to prune the search space and filter the proteins into same folds interacting with similar or different ligands and different folds interacting with similar or different ligands.
The TOPS model [^{1},^{2},^{25},^{26}] represented protein structures at the fold level by a graph where the nodes stand for SSEs(up or down) alphahelices and betastrandsand (nondirected) edges represent right or lefthanded chirality and parallel or antiparallel hydrogenbond relationships. In addition, there is a total ordering over the nodes, corresponding to the backbone of the protein. Our TOPS+ model enhances the original TOPS graph model with structural and biochemical features such as ligand interaction information and amino acid sequence length of the secondary structures. We have added extra nodes for loops (represented as a first class objectSSE) and ligands as well as maintaining the existing nodes for alphahelices and betastrands.
Further, we have designed a string model based on our TOPS+ graph model where the longrange and shortrange interactions between the SSEs are converted into incoming and outgoing arcs for each SSE, which maintain the directions and arc type properties. All relevant SSE nodes are enhanced with SSEligand interaction information which includes loopligand interaction information. We abstract away from the ligands themselves, to give a linear model called TOPS+ string which preserves the essential biochemical information whilst permitting more efficient and nonheuristic algorithms for comparison.
In detail, each SSE node is enriched with SSEType, SSE segment details are indicated by PSSL where PS is the PDB start number and SL is the SSE segment length, total incoming arcs (InArc) and total outgoing arcs (OutArc), total number of ArcTypes, and total number of ligand arcs (LigArc). The SSEType is given by {E, e, H, h, U, u} where, 'E' and 'e' represent the 'up' and 'down' oriented betastrands; 'H' and 'h' indicate the 'up' and 'down' oriented alphahelices; 'U' and 'u' represent ligandbound and ligandfree loops. The InArcType is represented as a {R, L, P, A}, where 'R' and 'L' represent right and left chiralities; 'P' and 'A' represent parallel and antiparallel hydrogenbonds respectively. The OutArcType is represented in a similar manner by {R', L', P', A'}. Ligand arcs are indicated by LT = AA where LT is the ligand type and AA is the PDB number. For example, Figure 4(a) and 4(b) shows the representation of the TOPS+ model and the abstract TOPS+ string representation for the protein domain 1fnb01 respectively. Here the triangles represent the beta strands; red curves represent the alpha helix; circles indicate loop regions and green arcs indicate hydrogen bonds between two beta strands, called the antiparallel beta sheet. In the 3D cartoon of 1fnb01 shown in Figure 4(c), the ligand molecules 'FAD' and 'SO4' are indicated by spheres. We have computed the topological information for a given PDB based on either SCOP or CATH domain definitions by using the TOPS cartoon generation method [^{27}]. The InterCal program (Wallace and Michalopoulos, personal communication) provides the proteinligand interaction information. Thus, by combining topological details and the proteinligand interaction information, we have constructed our TOPS+ and TOPS+ strings models for a given protein domain; for more details refer to Veeramalai and Gilbert [^{18}]. We have generated the TOPS+ model and the TOPS+ strings representation of 28 976 and 28 298 protein domains including 16 163 and 14 887 ligandbound protein domains corresponding to the CATH 2.4 and the SCOP 1.61, respectively.
Our TOPS+ comparison method computes a comparison score between two proteins based on edit distance using a dynamic programming approach. The Levenshtein distance or edit distance [^{28},^{29}] gives a measure (the cost) of the minimum number of elementary edit operations (insertions, deletions and substitutions of characters) necessary to transform one sequence into the other. In this research we have improved our existing method using parameter optimization in the dynamic programming table computation and also in the computation of the comparison score.
We have added some additional functions to the standard edit distance algorithm in order to compare the TOPS+ strings models of the two protein structures. Our optimized comparison method comprises five major steps as follows, where steps (1) and (5) take parameter tuning tables and penalty weights for arc information (both topological arcs and ligand arcs) into account:
1. Recursive definition of the optimal dissimilarity score for match and mismatch between TOPS+ strings elements (this process is based on the advanced_SSEArc+Match function, which incorporates the parameter optimization process using parameter tuning table).
2. Construction of the Edit Distance (ED) matrix (dynamic programming table).
3. Traceback on the ED matrix (dynamic programming table).
4. Obtain the LCS (Longest Common Substring), which is equivalent to the largest common structural core.
5. Computation of the comparison score based on the compression measure which is optimized with penalty weights for arc information (at both topological arcs and ligand arcs).
In our optimized TOPS+ comparison method, the computation of the edit distance matrix M is an important process, in which the advanced_SSEArc+Match function plays a key role in assigning dissimilarity scores for each TOPS+ strings element match or mismatch between the target t_{i }? T and the source s_{j }? S. This function handles the parameter optimization process while computing the construction of the editdistance matrix using a dynamic programming approach. It takes the basic parameter list P_{b }supplied together with the input and constructs the parameter tuning table PT with 12 weights (w_{1 }to w_{12}) and integrates these weights with the absolute arc differences (D_{1 }to D_{12}) between the TOPS+ strings elements t_{i }? T and s_{j }? S, computing the final normalized dissimilarity score for match or mismatch between the t_{i }? T and s_{j }? S. In each step the advanced_SSEArc+Match function performs the following processes in order to obtain the dissimilarity scores between each pair of TOPS+ strings elements of T and S and to construct the dynamic programming table:
? Construct the parameter tuning table PT based on the basic parameter list P_{b }and this process performed once.
? Compute the absolute difference for the arc features such in/out/ligand arc between t_{i }and s_{j }of T and S respectively.
? Compute the optimized dissimilarity score for t_{i }and s_{j }match using equations (1) and (2) below.
? Construct the dynamic programming table.
Algorithm 1 (Edit distance between TOPS+ strings) A function call ComputeEditDistance(T, S) computes the edit distance matrix M, the backtrace pointer matrix P, the edit distance value ed, and the longest common subsequence lcs of two TOPS+ strings T and S.
???function ComputeEditDistance(T = t_{1}, ..., t_{n}, S = s_{1}, ..., s_{m})
??????M [0, 0] ? 0
??????for i ? 1, ..., n do
?????????M [i, 0] ? i
??????for j ? 1, ..., m do
?????????M [0, j] ? j
??????for i ? 1, ... n do
?????????for j ? 1, ..., m do
????????????A ? SSEArc+Match(t_{i}, s_{j})
????????????M [i, j] ? min{M [i, j  1] + 1, M [i, j  1] + 1, M [i  1, j  1] + A}
????????????if M [i, j] = M [i  1, j  1] + A then
???????????????P [i, j] ? 'm'???? match or mismatch of s_{i }to t_{j}
????????????else if M [i, j] = M [i, j  1] + 1 then
???????????????P [i, j] ? 'i'???? insertion of s_{j }into t
????????????else
???????????????P [i, j] ? 'd'???? deletion of t_{i }from t
??????ed ? M [n, m]
??????lcs ? BuildLCS(M, P, T, S)
??????return ?M, P, ed, lcs?
???function BuildLCS(M, P, T, S)
??????lcs ? ????? empty sequence
??????k ? 0
??????i ? n ? length of T
??????j ? m ? length of S
??????while i > 0 or j > 0 do
?????????if P [i, j] = 'm' then
????????????lcs ? lcs ? t_{j1 }? match or mismatch of s_{i1 }to t_{j1}
????????????k ? k + 1
????????????i ? i  1
????????????j ? j  1
?????????else if P [i, j] = 'd' then
????????????i ? i  1???? deletion of t_{i1 }from t
?????????else
????????????j ? j  1???? insertion of s_{j1 }into t
??????return lcs
???function SSEArc+Match(t_{i}, s_{j})
??????mS ? 0
??????Parse(t_{i}, t_{sk}, t_{I}, t_{O}, t_{L}, t_{IR}, t_{IL}, t_{IP}, t_{IA}, t_{OR}, t_{OL}, t_{OP}, t_{OA})
??????Parse(s_{j}, s_{sk}, s_{I}, s_{O}, s_{L}, s_{IR}, s_{IL}, s_{IP}, s_{IA}, s_{OR}, s_{OL}, s_{OP}, s_{OA})
??????if MatchSSEArc+features(t_{i}, s_{j}) then
?????????mS ? mS + 1
??????return mS
???procedure Parse(t_{i}, t_{sk}, t_{I}, t_{O}, t_{L}, t_{IR}, t_{IL}, t_{IP}, t_{IA}, t_{OR}, t_{OL}, t_{OP}, t_{OA})
??????t_{sk }? secondary structure length of t_{i}
??????t_{I}, t_{O}, t_{L }? total number of incoming, outgoing, ligand arcs of t_{i}
??????t_{IR}, t_{IL}, t_{IP}, t_{IA }? total number of incoming arcs of type R, L, P, A of t_{i}
??????t_{OR}, t_{OL}, t_{OP}, t_{OA }? total number of outgoing arcs of type R, L, P, A of t_{i}
The time complexity is O(n^{2}) where n is the length of the string of SSEs. The current version of our TOPS+ method performs global alignment [^{30}]. However, local alignment [^{31}] can be applied to find the local structural similarity or patterns such as similar SSEligand interactions at local level across different folds.
We performed parameter tuning/optimization in order to obtain the optimal approximate match between two protein structures. In general, at the superfamily level, only core structures are conserved throughout evolution across the members of protein families. Studies have shown that the number of SSE insertions and deletions is variable for different sequence families or organisms [^{32}]. This implies that variable numbers of 'indels' are applicable to the ArcsTypes and SSETypes across protein families from various organisms within a superfamily. Thus, it is important to develop a cost matrix with an additional penalty scoring function for such an approximate matching process. In the following sections we discuss the development of the parameter tuning table and the computation of the absolute difference between ArcTypes types and SSETypes. Subsequently, we explain the main parameter optimization process involved in the computation of the dynamic programming table, which exploits the computation of a normalized dissimilarity score for TOPS+ strings element match.
The parametertuning table PT contains a list of weights for all the SSEArc+ features given by w_{1}, ..., w_{12}, where, w_{1 }is applied to SSEType; w_{2 }to total incoming arcs; w_{3}, w_{4}, w_{5 }and w_{6 }to total incoming arcs of type R, L, P and A respectively; w_{7 }to total outgoing arcs; w_{8}, w_{9}, w_{10 }and w_{11 }to total outgoing arcs of type R, L, P and A respectively; and w_{12 }to total ligands arcs. Every element 'y' in PT is in the range 1..9. The SSEType weight w_{1 }is obtained from the ISM or DSM scoring matrices (see Table 4). All the other weights w_{2 }to w_{12 }are computed using a basic parameter set P_{b }which is supplied with the input. The basic parameter P_{b }consists of a list of five positive integers in the range 1..9 given by P_{b }= (p, q, r, s, t). The parameter values p and q correspond to the parallel and the antiparallel hydrogenbond arcs respectively; r and s correspond to righthanded and lefthanded chiralities respectively; t is for the ligand arcs. Based on these conditions, we have constructed a total of 1,134 unique lists of basic parameters P_{b }used for tuning the parameter table PT with weights as follows:
? w_{3 }= w_{8 }= r (for incoming and outgoing arc type_R)
? w_{4 }= w_{9 }= s (for incoming and outgoing arc type_L)
? w_{5 }= w_{10 }= p (for incoming and outgoing arc type_P)
? w_{6 }= w_{11 }= q (for incoming and outgoing arc type_A)
? w_{2 }= w_{7 }= r + s + p + q (for total incoming and outgoing arcs)
? w_{12 }= t (for total ligand arcs)
The normalized dissimilarity score for SSE is obtained from the absolute difference of the SSE length (amino acid sequence length of the SSE) between the TOPS+ strings elements t_{i }? T and s_{j }? S given by,
where, t_{sk }and s_{sk }are the secondary structure length for t_{i }? T and s_{j }? S respectively. The absolute differences of the total number of incoming, outgoing and ligand arcs between the TOPS+ strings elements t_{i }? T and s_{j }? S are computed based on the equations given in Table 5. This process includes the computation of the absolute arc differences for the incoming and outgoing arcs according to their ArcTypes {R, L, P, A} between t_{i }? T and s_{j }? S of the TOPS+ strings elements.
The absolute differences for total incoming arcs are given by,
where t_{I }and s_{I }are the total number of incoming arcs of the TOPS+ strings elements t_{i }? T and s_{j }? S respectively. The absolute arc differences between incoming arcs, outgoing arcs and their arc types are calculated based on Table 5.
We have combined the parameter tuning table PT and the absolute arc differences to compute the final normalized dissimilarity score mSn for the matching or mismatching of the TOPS+ strings elements t_{i }? T and s_{j }? S. The general formula for calculating the normalized dissimilarity score is given in equation (3), where all the weights w_{i }are obtained from the parameter tuning table PT; all the values for absolute arc differences D_{i }are computed based on equation (2) and Table 5. Since the SSETypes {H, h} do not have any Hbond arcs {P, A}, we have substituted zero weights for corresponding arc weights as follows: w_{IP }= w_{IA }= w_{OP }= w_{OA }= 0. Note that this condition is applicable only when both t_{i }? T and s_{j }? S have SSEType {H, h}. The final normalized dissimilarity score for the helix SSETypes helix_mSn is calculated using equation (5). A similar condition is also applicable to the loops; they do not have any other arcs except their ligand arcs and their normalized similarity score loop_mSn is given by equation (6). The normalized dissimilarity score strand_mSn for the beta strands is calculated using equation (4). For all the other nonmatch SSEType we have considered the weights obtained from the ISM or DSM scoring matrices (see Table 4). When we applied these dissimilarity scores to loops we consider both 'u' (loop without ligand interaction) and 'U' (loop with ligand interaction) as identical. However when we perform ligand pattern discovery in a different context (for instance, similar ligand interaction on different fold types) for each SSE we have considered them as different SSETypes [^{24}].
We have computed a pairwise comparison score based on a compression measure (itself a metric) to evaluate the goodness of patterns for a set of TOPS+ strings models of the proteins. This procedure was adapted from Brazma et al. [^{33}] and has been employed in the TOPS comparison and pattern discovery methods [^{1},^{2},^{26}]. We have calculated the compression value based on the total number of proteins in the given set (in this case there are two proteins), the total number of SSEs in the proteins in the given input set and total number of SSEs in the common pattern. Similarly, we also take into account arc information such as total number of in/out/ligand arcs in the given set of proteins as well as the same information for the LCS pattern. The raw compression value over SSEs is given as:
where, S_{i} is the total number of SSEs in a protein i, n is the total number of proteins in the target set and S_{p} is the total number of SSEs in the LCS pattern. The normalized compression score over SSEs, which varies from 1 (best) to 0 (worst), is computed by:
Similarly we compute the normalized compression (nC) score for InArcs, OutArcs and LigandArcs. The overall normalized TOPS+ comparison score for TOPS+ strings models is computed by (9), combining the compression for the SSEs and arc type, optimized by different weights k_{1}, ..., k_{4 }(see below).
We have computed 17 different combinations of compression values based on ED and LCS together with or without different levels of SSEArc+ features information. Supplementary Table 1 (see supplementary material page at http://balabio.dcs.gla.ac.uk/mallika/WebTOPS/optTOPSplusresults.html) gives all the 17 normalized compression scores we have calculated based on ED and LCS from our advanced TOPS+ comparison with (output order of the results) and description.
We have performed training and analysis of our advanced TOPS+ comparison method with the parameter tuning table. Our method incorporates parameter optimization at two levels, both in the computation of the dynamic programming table and in the computation of the normalized compression measure. We have tested our method with 1,134 unique basic parameter lists on the training dataset of 7,000 random protein domain pairs from the PDB40 dataset, which contain both ligandbound and ligandfree proteins. We validated our results with the SCOP superfamily classification numbers and obtained the ROC and AUC values corresponding to each basic parameter list. The experimental testing and evaluation analysis involved the following steps:
? Perform advanced TOPS+ comparison based on the advanced_SSEArc+Match function for all basic parameters in list P_{b}.
? Compute ROC (Receiver Operating Characteristic) curve analysis for all 7,000 results, and for each parameter list.
? Calculate the AUC values corresponding to the 17 different nC scores.
We have obtained the 17 ? 1, 134 AUC values and we plot the results according to the SCOP classes 14, which correspond to allalpha, allbeta, alpha/beta and alpha+beta, and all classes together. In Figure 5, the xaxis represents the 1,134 parameters and the yaxis denotes the AUC values corresponding to 17 nC scores and for each of the SCOP classes allalpha, allbeta, alpha/beta, alpha+beta, and all classes (17 ? 5 = 85 points). The color codes adjacent to the actual graph indicate the range of AUC values between 0 and 1. This corresponds to the percentage of accuracy where 0% is represented by blue and 100% by maroon. The first 250 basic parameter lists give the best performance for all the SCOP classes, and the AUC values based on nCE again give consistently better performance for all the parameter values compared to nC scores based on LCS. Specifically, nCEnA(6), nCEnL(8) and nCESAnL(10) give higher AUC values for all SCOP classes and for all the parameter values; while nCLnL(14) and nCLSAnL(16) give higher performance for the allalpha class. From the parameter tuning table evaluation analysis we have found that the basic parameter P_{b }[P92:3,1,1,1,1] always provides the best result for our training dataset. So we have selected this basic parameter as a default parameter for our advTOPS+ method and performed further analysis and produced the results for PDB40 and ChewKedem datasets.
We have considered the PDB40 representative subset of 4,261 protein domain pairs (which excludes the domain pairs used in our training dataset) corresponding to SCOP 1.61, see Table 6. These proteins contain both ligandbound and ligandfree protein domains. We have compared our advanced TOPS+ method (based on basic parameter value [P92:3,1,1,1,1]) with our basic TOPS+ and TOPS [^{3},^{26}] methods. The basic TOPS+ comparison method computes the edit distance between two proteins based on TOPS+ strings elements using a dynamic programming approach.
We have considered the Chew and Kedem dataset [^{4},^{5}] to assess the biological significance of our advanced TOPS+ comparison method. This dataset contains 36 medium size representative proteins of five different families: globins (17 entries), alphabeta (6 entries), timbarrels (4 entries), allalpha (2 entries), and allbeta (7 entries) proteins. We compared our method against the SSAP structure alignment program [^{9},^{34},^{35}] and TOPS [^{3},^{26}] and validated our results based on computation of the Fmeasure [^{20},^{21}].
For the PDB40 dataset we have performed evaluation analysis as given below:
? Obtain the pairwise comparison score from the protein comparison method for a given dataset.
? Assignment of Homolog (TP, true positive) and nonHomolog (FP, false positive) based on the SCOP superfamily numbers for each protein domain (see below) and rank them according to the comparison score.
? Perform the Receiver Operating Characteristics (ROC) curve analysis for equation (2). For all the ROC curves we have computed the AUC (Area Under the ROC Curve) values.
We have considered the assignment of homologous or nonhomologous information of a protein domain pair, based on the standard SCOP classification numbers at superfamily level as an indication of structural homology. If both protein domains belong to same superfamily then they are homologous, otherwise they are nonhomologous.
We have obtained allagainstall comparison scores from all the comparison methods and based on these scores, for each method, we performed pairwise hierarchical clustering using the OC program [^{36}]. We have evaluated the biological significance of the clusters obtained from different protein structure comparison methods based on Fmeasure calculations [^{20},^{21}].
We performed all the analyses using a RedHat 7.2 linux environment with an Intel Pentium IV 1.6 GHz processor. The methods SSAP, TOPS, TOPS+ and advTOPS+ took 9139 s, 75 s, 21 s and 1805 s (s = seconds) respectively to complete 630 pairwise comparisons.
MV developed the advanced TOPS+ comparison method, performed the calculations and prepared the figures. All authors prepared the manuscript, contributed to the discussion, and have approved the final manuscript.
We would like to thank the TOPS Project for TOPS data resources and Juris Viksna for advice on our method. MV has been supported by a PhD studentship from the University of Glasgow.
References
Gilbert D,Westhead DR,Viksna J,Thornton J,A Computer System to Perform Structure Comparison using TOPS Representations of Protein StructureJ Comput ChemYear: 200126233010.1016/S00978485(01)000961  
Gilbert D,Westhead DR,Viksna J,Frasconi P, Shamir RTechniques for Comparison, Pattern Matching and Pattern Discovery: From Sequences to Protein TopologyArtificial Intelligence and Heuristic Methods in Bioinformatics, Volume 183 of NATO Science Series: Computer & Systems SciencesYear: 2003IOS Press128147  
Viksna J,Gilbert D,Pattern Matching and Pattern Discovery Algorithms for Protein TopologiesAlgorithms in BioInformatics, Volume 2149 of Lecture Notes in Comput. SciYear: 2001SpringerVerlag98111  
Krasnogor N,Pelta DA,Measuring the Similarity of Protein Structures by Means of the Universal Similarity MetricBioinformaticsYear: 20042071015102110.1093/bioinformatics/bth03114751983  
Chew LP,Kedem K,Finding the Consensus Shape for a Protein FamilyAlgorithmicaYear: 20033811512910.1007/s0045300310452  
GoldsmithFischman S,Honig B,Structural Genomics: Computational Methods for Structure AnalysisProtein SciYear: 20031291813182110.1110/ps.024290312930981  
Berman HM,Westbrook J,Feng Z,Gilliland G,Bhat TN,Weissig H,Shindyalov IN,Bourne PE,The Protein Data BankNucleic Acids ResYear: 20002823524210.1093/nar/28.1.23510592235  
Orengo CA,Taylor WR,A Rapid Method for Protein Structure AlignmentJ Theor BiolYear: 1990147451755110.1016/S00225193(05)8026322074728  
Taylor WR,Orengo CA,Protein Structure AlignmentJ Mol BiolYear: 198920812210.1016/00222836(89)9008432769748  
Russell RB,Barton GJ,Multiple Protein Sequence Alignment from Tertiary Structure Comparison: Assignment of Global and Residue Confidence LevelsProteinsYear: 199214230932310.1002/prot.3401402161409577  
Holm L,Sander C,Protein Structure Comparison by Alignment of Distance MatricesJ Mol BiolYear: 199323312313810.1006/jmbi.1993.14898377180  
Shindyalov IN,Bourne PE,Protein Structure Alignment by Incremental Combinatorial Extension (CE) of the Optimal PathProtein EngineeringYear: 199811973974710.1093/protein/11.9.7399796821  
Madej T,Mossing MC,Hamiltonians for Protein Tertiary Structure Prediction Based on Threedimensional Environment PrinciplesJ Mol BiolYear: 1993233348048710.1006/jmbi.1993.15257692069  
Madej T,Gibrat JF,Bryant SH,Threading a Database of Protein CoresProteinsYear: 199523335636910.1002/prot.3402303098710828  
Harrison A,Pearl F,Sillitoe I,Slidel T,Mott R,Thornton J,Orengo C,Recognizing the Fold of a Protein StructureBioinformaticsYear: 200319141748175910.1093/bioinformatics/btg24014512345  
Grindley HM,Artymiuk PJ,Rice DW,Willett P,Identification of Tertiary Structure Resemblance in Proteins Using a Maximal Common Subgraph Isomorphism AlgorithmJ Mol BiolYear: 1993229370772110.1006/jmbi.1993.10748381875  
Koch I,Lengauer T,Wanke E,An Algorithm for Finding Maximal Common Subtopologies in a Set of Protein StructuresJ Comput BiolYear: 19963228930610.1089/cmb.1996.3.2898811488  
Veeramalai M,Gilbert D,A Novel Method for Comparing Topological Models of Protein Structures Enhanced with Ligand InformationBioinformaticsYear: 200824232698270510.1093/bioinformatics/btn51818842602  
Veeramalai M,Ye Y,Godzik A,TOPS++FATCAT: fast flexible structural alignment using constraints derived from TOPS+ Strings ModelBMC BioinformaticsYear: 2008935810.1186/14712105935818759993  
Handl J,Knowles J,Kell DB,Computational Cluster Validation in PostGenomic Data AnalysisBioinformaticsYear: 200521153201321210.1093/bioinformatics/bti51715914541  
van Rijsbergen CJ,Information RetrievalYear: 19792London: Butterworths  
Krishna SS,Grishin NV,Structural Drift: A Possible Path to Protein Fold ChangeBioinformaticsYear: 20052181308131010.1093/bioinformatics/bti22715604105  
Harrison A,Pearl F,Mott R,Thornton J,Orengo C,Quantifying the Similarities within Fold SpaceJ Mol BiolYear: 2002323590992610.1016/S00222836(02)00992012417203  
Veeramalai M,A Novel Method for Comparing Topological Models of Protein Structures Enhanced with Ligand InformationPhD thesisYear: 2005University of Glasgow  
Michalopoulos I,Torrance GM,Gilbert D,Westhead DR,TOPS: An Enhanced Database of Protein Structural TopologyNucleic Acids ResYear: 200332D251254  
Torrance GM,Gilbert D,Michalopoulos I,Westhead DR,Protein Structure Topological Comparison, Discovery and Matching ServiceBioinformaticsYear: 200521102537253810.1093/bioinformatics/bti33115741246  
Westhead D,Slidel T,Flores T,Thornton J,Protein structural topology: automated analysis and diagrammatic representationProtein ScienceYear: 1999889790410211836  
Levenshtein VI,Binary Codes Capable of Correcting Deletions, Insertions, and ReversalsSoviet Physics DokladyYear: 1966108707710  
Valiente G,Combinatorial Pattern Matching Algorithms in Computational Biology using Perl and RYear: 2009Taylor & Francis/CRC Press  
Needleman SB,Wunsch CD,A General Method applicable to the Search for Similarities in the Amino Acid Sequence of two ProteinsJ Mol BiolYear: 197048344345310.1016/00222836(70)9005745420325  
Smith TF,Waterman MS,Identification of Common Molecular SubsequencesJ Mol BiolYear: 198114719519710.1016/00222836(81)9008757265238  
Mizuguchi K,Blundell TL,Analysis of conservation and substitutions of secondary structure elements within protein superfamiliesBioinformaticsYear: 200016121111111910.1093/bioinformatics/16.12.111111159330  
Brazma A,Jonassen I,Vilo J,Ukkonen E,Pattern Discovery in BiosequencesProc. 4th Int. Coll. Grammatical Inference, Volume 1433 of Lecture Notes in Comput. SciYear: 1998SpringerVerlag257270  
Orengo CA,Brown NP,Taylor WR,Fast Structure Alignment for Protein Databank SearchingProteinsYear: 199214213916710.1002/prot.3401402031409565  
Orengo CA,Taylor WR,SSAP: Sequential Structure Alignment Program for Protein Structure ComparisonMethods EnzymolYear: 1996266617635 full_text. 8743709  
Barton GJ,OCA Cluster Analysis ProgramYear: 2002http://www.compbio.dundee.ac.uk/downloads/oc/ 
Figures
Tables
ROC curve and Fmeasure analysis of structural homology for the PDB40 dataset.
SCOP Class  TOPS  TOPS+  advTOPS+  

1  All alpha  0.76/0.79  0.83/0.85  0.82/0.88 
2  All beta  0.89/0.85  0.85/0.83  0.87/0.86 
3  Alpha/beta  0.82/0.75  0.75/0.70  0.77/0.70 
4  Alpha+beta  0.84/0.75  0.84/0.74  0.90/0.81 
AUC/Fmeasure values for all the comparison methods for the PDB40 dataset.
Biological significance of protein domain clusters for the ChewKedem dataset.
Method  Fmeasure 

SSAP  0.966 
TOPS  0.955 
TOPS+  0.931 
advTOPS+  0.985 
Fmeasure of the clusters obtained from different protein structure comparison methods.
advTOPS+ comparison scores for the ChewKedem dataset.
Protein Fold  Domain  SSE Ln  LCS PAT Ln  Adv TOPS+ Score  LCS SSE PATTERN 

Alphabeta  d1aa9__  23  20  0.49  uEUhuEuhuEuhuEUhuEhu 
d1gnp__  25  20  0.51  uEUhuEhuEuhuuEUhuEhu  
d6q21a_  21  16  0.59  uEUhuEuhuEUhuEhu  
d1qraa_  25  20  0.51  ueUHueHueuHuueUHueHu  
d5p21__  25  20  0.51  ueUHueHueuHuueUHueHu  


TIMbarrel  d6xia__  60  40  0.30  uhuhuHuHueuHueuHuhuehhuHuHueuuuHuhuhuuhu 
d2mnr_1  37  29  0.37  uuuuHueuHueuHueuHueuHuuhuuuhu  
d1chra1  41  31  0.35  uuuHuHueuHueuHueuHuHueuHuhuuuhu  
d4enl_1  55  35  0.37  uhuhuHuueuuuHuhueuHuHuuhuuHuuHuHuhu 
advTOPS+ comparison scores for alphabeta and TIMbarrel folds in the ChewKedem dataset.
Identity and dissimilarity scoring matrices for TOPS+ diagrams.
ISM  DSM  

SSE  E  e  H  h  U  u  E  e  H  h  U  u 


E  0  1  1  1  1  1  0  1  2  2  2  2 
e  1  0  1  1  1  1  1  0  2  2  2  2 
H  1  1  0  1  1  1  2  2  0  1  2  2 
h  1  1  1  0  1  1  2  2  1  0  2  2 
U  1  1  1  1  0  1  2  2  2  2  0  1 
u  1  1  1  1  1  0  2  2  2  2  1  0 
Identity Scoring Matrix (ISM) and Dissimilarity Scoring Matrix (DSM) of secondary structure elements used for matching TOPS+ diagrams.
Normalized similarity score between secondary structure elements.
Absolute Differences  Equations  Description 

total incoming arcs  D_{2 }= t_{I }s_{I}  t_{I }and s_{I }are the total number of incoming arcs of the TOPS+ strings elements t_{i }? T and s_{j }? S respectively 


total incoming arcs type_R  D_{3 }= t_{IR }s_{IR}  t_{IR }and s_{IR }indicate the total number of incoming arcs type_R for t_{i }? T and s_{j }? S respectively 


total incoming arcs type_L  D_{4 }= t_{IL }s_{IL}  t_{IL }and s_{IL }indicate the total number of incoming arcs type_L for t_{i }? T and s_{j }? S respectively 


total incoming arcs type_P  D_{5 }= t_{IP }s_{IP}  t_{IP }and s_{IP }indicate the total number of incoming arcs type_P for t_{i }? T and s_{j }? S respectively 


total incoming arcs type_A  D_{6 }= t_{IA }s_{IA}  t_{IA }and s_{IA }indicate the total number of incoming arcs type_A for t_{i }? T and s_{j }? S respectively 


total outgoing arcs  D_{7 }= t_{O }s_{O}  t_{O }and s_{O }are the total number of outgoing arcs for the TOPS+ strings elements t_{i }? T and s_{j }? S respectively 


total outgoing arcs type_R  D_{8 }= t_{OR }s_{OR}  t_{OR }and s_{OR }indicate the total number of outgoing arcs type_R for t_{i }? T and s_{j }? S respectively 


total outgoing arcs type_L  D_{9 }= t_{OL }s_{OL}  t_{OL }and s_{OL }indicate the total number of outgoing arcs type_L for t_{i }? T and s_{j }? S respectively 


total outgoing arcs type_P  D_{10 }= t_{OP }s_{OP}  t_{OP }and s_{OP }indicate the total number of outgoing arcs type_P for t_{i }? T and s_{j }? S respectively 


total outgoing arcs type_A  D_{11 }= t_{OA }s_{OA}  t_{OA }and s_{OA }indicate the total number of outgoing arcs type_A for t_{i }? T and s_{j }? S respectively 


total ligand arcs  D_{12 }= t_{L }s_{L}  t_{L }and s_{L }are the total number of ligand arcs for the TOPS+ strings elements t_{i }? T and s_{j }? S respectively 
Equations used for computing the normalized similarity score between secondary structure elements of TOPS+ strings.
Structural homology of protein domains for the PDB40 dataset.
SCOP Class  Hom  %  NonHom  %  Total  % 

All alpha  129  69  58  31  187  10 
All beta  219  68  102  32  321  18 
Alpha/beta  452  48  487  52  939  52 
Alpha+beta  167  46  193  54  360  20 


Total  967  54  840  46  1807  100 
Superfamily homologous (Hom) and nonhomologous (NonHom) statistics of the PDB40 dataset.
Article Categories:

Previous Document: Effects of prolonging administration gonadotropin on unexpectedly poor ovarian responders undergoing...
Next Document: Comparison of heterologous xylose transporters in recombinant Saccharomyces cerevisiae.