Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Genomics**|**v.14(Suppl 2); 2013**|**PMC3582452

Formats

Article sections

- Abstract
- Background
- Methods
- Testing results
- Discussion
- Competing interests
- Authors' contributions
- References

Authors

Related links

BMC Genomics. 2013; 14(Suppl 2): S1.

Published online 2013 February 15. doi: 10.1186/1471-2164-14-S2-S1

PMCID: PMC3582452

Cody Ashby,^{1} Daniel Johnson,^{1} Karl Walker,^{2} Iyad A Kanj,^{3} Ge Xia,^{}^{4} and Xiuzhen Huang^{}^{5}

Cody Ashby: ude.etatsa.liams@ybhsa.ydoc; Daniel Johnson: ude.etatsa.liams@nosnhojleinad; Karl Walker: ude.rlau@1reklawak; Iyad A Kanj: ude.luaped.mdc@jnaki; Ge Xia: ude.etteyafal@gaix; Xiuzhen Huang: ude.etatsa@gnauhx

Selected articles from ISCB-Asia 2012

Victor Jin and Paul Horton

Publication of this supplement was funded by the authors.

ISCB-Asia 2012

17-19 December 2012

Shenzhen, China

Copyright ©2013 Ashby et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Protein structure comparison and classification is an effective method for exploring protein structure-function relations. This problem is computationally challenging. Many different computational approaches for protein structure comparison apply the secondary structure elements (SSEs) representation of protein structures.

We study the complexity of the protein structure comparison problem based on a mixed-graph model with respect to different computational frameworks. We develop an effective approach for protein structure comparison based on a novel independent set enumeration algorithm. Our approach (named: ePC, **e**fficient **e**numeration-based **P**rotein structure **C**omparison) is tested for general purpose protein structure comparison as well as for specific protein examples. Compared with other graph-based approaches for protein structure comparison, the theoretical running-time *O*(1.47^{rn}n^{2}) of our approach ePC is significantly better, where *n *is the smaller number of SSEs of the two proteins, *r *is a parameter of small value.

Through the enumeration algorithm, our approach can identify different substructures from a list of high-scoring solutions of biological interest. Our approach is flexible to conduct protein structure comparison with the SSEs in sequential and non-sequential order as well. Supplementary data of additional testing and the source of ePC will be available at http://bioinformatics.astate.edu/.

Protein structure comparison is an effective method for exploring protein structure-function relations and for studying evolutionary relations of different species. It can also be applied to identify the active sites of carrier proteins, the binding sites of antibodies, the inhibition sites of enzymes, and the common structural motifs of proteins, which has significant applications in biological and biomedical research.

The computational methods for protein structure comparison usually represent a protein structure by atomic coordinates in the Euclidean space, as a distance matrix [1] whose entries represent the distances between two residues of the protein, or as a contact map [2], where a binary matrix is used to represent the distances between the residue pairs. A structure graph representation of a protein tertiary structure was first defined in [3] for protein structure prediction. In this current work, we adopt the structure graph representation in [3]. We develop a very efficient graph-based approach for protein structure comparison. Our approach transforms the comparison problem to an independent set problem in an auxiliary graph, and then applies a novel enumeration algorithm to identify the best out of a set of good comparison candidates.

We first show the problem of comparing a query structure to another structure is intractable with respect to several computational frameworks. For example, we show that the problem is $\mathcal{N}\mathcal{P}-\mathsf{\text{hard}}$ (even for very restricted instances), cannot be approximated to a ratio ${n}^{\frac{1}{2}-\epsilon}$, for any *ε *> 0, unless $\mathcal{P}=\mathcal{N}\mathcal{P}$, and is *W *[1]- complete with respect to the framework of parameterized complexity. We also show that a useful case of the problem is solvable in polynomial time by reducing it to the 2-CNF-SATISFIABILITY problem.

Whereas the above results are negative hinting at the challenging nature of the problem, the graph-based approach we use allows us to model the problem as a maximum independent set problem, for which a repertoire of effective exact algorithms exist in the literature. We use an algorithm developed by (some of) the authors [4] to enumerate the top-*K *maximum independent sets in a graph in time *O*(1.47^{n}n^{2}), where *n *is the number of vertices in the graph (Note that the algorithm in [4] enumerates the top-*K *minimum vertex covers in a graph, but obviously can be used to enumerate the top-*K *maximum independent sets in a graph using the standard reduction between vertex cover and independent set); this enumeration algorithm allows us to sift through the top SSE alignments for the protein structure comparison problem, looking for the best amongst them in terms of accuracy. Compared with other graph-base approaches, the theoretical running-time *O*(1.47^{rn}n^{2}) of our approach ePC is the current best, where n is the smaller number of SSEs of the two proteins, *r *is an introduced parameter of small values.

Many different approaches for protein structure comparison apply the secondary structure elements (SSEs) representation and database searching, such as deconSTRUCT [5], SSM [6], GANGSTA [7], MASS [8,9], VAST [10], TOPS [11] and approaches in [12-19]. Our approach ePC utilizes the SSE-based representation of the protein structure, and takes into consideration the global 3D structural arrangements of the SSEs of the proteins. We compare our approach with two other SSE-based approaches: deconSTRUCT, an approach for general purpose protein structure comparison and database search, and SSM, a high-resolution structure comparison approach. Our approach has comparable performance as deconSTRUCT. With a more general and simplified representation and a unified graph enumeration algorithm, our approach could detect a substructure or motif structure in a set of large structures, more than one common substructure shared by a set of proteins. It is very flexible. Our approach could use a wide range of evaluation functions for protein structure comparison. It could be applied to handle sequential and non-sequential order of SSE alignment and be extended to handle challenging protein multiple structure alignment and protein subset alignment.

A mixed graph for a protein structure is constructed from the PDB file as follows: each vertex represents a core/secondary structure element (i.e., an alpha helix element, or, a beta strand element), each undirected edge represents the interaction between two cores, and each directed edge (arc) represents the loop between two consecutive cores (from the N-terminal to the C-terminal). A mixed graph representation is used for protein structure prediction in [3]. The DSSP program [20,21] was used for the assignments of secondary structure elements for the protein entries from the Protein Data Bank (PDB). Refer to the protein structure and the corresponding mixed graph representation in Figures Figures11 for protein with ID: 6ldh. Alpha helix elements are represented by circles and beta strand elements are represented by squares. Therefore, a mixed graph can be represented as a triple *G *= (*V *(*G*), *A *(*G*), *E *(*G*)), where *V *(*G*) is the vertex-set of *G*, *E *(*G*) is the set of undirected edges of *G*, and *A *(*G*) is the set of directed edges of *G*, which induces a directed path spanning all vertices of *G*, thus defining a linear order among the vertices of *V *(*G*). The aforementioned mixed graph representation incorporates the SSE type, the sequential order of the SSEs, and the interactions of the SSEs. When comparing two protein structures, the problem could now be reduced to finding the common subgraph of the two mixed graph.

Goldman et al. [2] studied the protein comparison problem using the notion of *contact maps*. Contact maps are undirected graphs whose vertices are linearly ordered. Goldman et al. [2] formulated the protein comparison problem as a CONTACT MAP OVERLAP problem, in which we are given two contact maps and we need to identify a subset of vertices *S *in the first contact map, a subset of vertices ${S}^{\text{'}}$ in the second with $\left|S\right|=\phantom{\rule{2.77695pt}{0ex}}|{S}^{\text{'}}$, and an order preserving (w.r.t. linear ordering) bijection $f:S\to {S}^{\text{'}}$, such that the number of edges in *S *(i.e., between the vertices in *S*) that correspond (under *f*) to edges in ${S}^{\text{'}}$ is maximized. In [2], the authors proved that the CONTACT MAP OVERLAP problem is MAXSNP-complete even when both contact maps have maximum degree one.

Song et al. [3] studied the problem of mixed-graph comparison, when each vertex *v *in the first mixed-graph is associated with a subset of vertices *S _{v }*in the second mixed-graph, and the bijection

In this section, we study the complexity of the mixed graph embedding problem, which corresponds to the problem of identifying the query protein structure (e.g., a motif structure) as a substructure in a larger protein structure.

We define the GRAPH EMBEDDING problem as follows:

GRAPH EMBEDDING

Given two mixed graphs *G *= (*V *(*G*), *A*(*G*), *E*(*G*)) and *H *= (*V *(*H*), *A*(*H*), *E*(*H*)), where *H *is referred to as the *host graph*, such that each vertex *v * *V *(*G*) has a list *L*(*v*) *V *(*H*) of vertices in *H *that it can be mapped to, decide if there exists an injection *f*: *V *(*G*) → *V *(*H*) such that:

(i) *f *(*v*) *L*(*v*) for every *v * *V *(*G*);

(ii) for any two vertices *v*, ${\nu}^{\text{'}}\in V\left(G\right)$, there is a directed path from *v *to ${\nu}^{\text{'}}$ in G if and only if there is a directed path from f(*v*) to $f\left({\nu}^{\text{'}}\right)$ in *H*; and

(iii) for any two vertices *v*, ${\nu}^{\text{'}}\in V\left(G\right)$, if $\nu {\nu}^{\text{'}}\in E\left(G\right)$ then $f\left(\nu \right)f\left({\nu}^{\text{'}}\right)\in E\left(H\right)$.

We shall call an injective embedding *f *satisfying properties (i)-(iii) above a *valid embedding*.

Informally speaking, the GRAPH EMBEDDING problem asks if we can embed *G *into *H *in such a way that the precedence order determined by the arcs of *G *is respected by this embedding, and the undirected edges of *G *are respected by this embedding.

We define the restriction of the GRAPH EMBEDDING problem, denoted *r*-GRAPH EMBEDDING, where *r *is positive integer, by restricting the cardinality of the set *L*(*v*) to be at most *r*, for every *v * *V *(*G*); that is, in the restrictions of the problems, a vertex in *V *(*G*) can be mapped to at most *r *vertices in *H*.

If one cannot embed the whole graph *G *into *H*, it is natural to seek an embedding that embeds the maximum number of vertices in *G *into *H*, while respecting conditions (i)-(iii) above. Therefore, we define a version of GRAPH EMBEDDING, denoted GRAPH EMBEDDING_{≥}, by introducing a nonnegative parameter *k*, and asking whether there exists a subset *S * *V *(*G*) with |*S*|≥ *k*, and an injection *f*: *S *→ *V *(*H*) such that:

(i) *f*(*v*) *L*(*v*) for every *v * *S*;

(ii) for any two vertices *v*, ${\nu}^{\text{'}}\in S$, if there is a directed path from *v *to ${\nu}^{\text{'}}$ in *G *then there is a directed path from *f*(*v*) to $f\left({\nu}^{\text{'}}\right)$ in *H*; and

(iii) for any two vertices *v*, ${\nu}^{\text{'}}\in S$, if $\nu {\nu}^{\text{'}}\in E\left(G\right)$ then $f\left(\nu \right)f\left({\nu}^{\text{'}}\right)\in E\left(H\right)$.

The optimization/maximization version of the GRAPH EMBEDDING_{≥ }problem, denoted MAXIMUM GRAPH EMBEDDING, asks for a set *S *of maximum cardinality that satisfies conditions (i)-(iii) above. Similarly, we can define the problems *r*-GRAPH EMBEDDING_{≥ }and MAXIMUM r-GRAPH EMBEDDING.

It was shown in [3] that a more general problem than *r*-GRAPH EMBEDDING, in which the set of edges *A*(*G*) do not necessarily induce a path, is $\mathcal{N}\mathcal{P}-\phantom{\rule{2.77695pt}{0ex}}\mathsf{\text{complete}}$ for any *r *≥ 3. The same proof actually shows that the *r*-GRAPH EMBEDDING problem is $\mathcal{N}\mathcal{P}-\phantom{\rule{2.77695pt}{0ex}}\mathsf{\text{complete}}$ for any *r *≥ 3. We show next that the 2-GRAPH EMBEDDING is solvable in polynomial time.

**Theorem 0.1 ***The *2-GRAPH EMBEDDING *problem is solvable in polynomial time*.

PROOF. We reduce the problem to 2-CNF-SATISFIABILITY, which is solvable in polynomial time (for example, see [22]. Recall that in the 2-CNF-SATISFIABILITY problem we are given a Boolean formula in the *conjunctive normal form *(CNF) (i.e., the formula is the conjunction of clauses, and each clause is the disjunction of a literals, which are variables or negations of variables), in which each clause contains at most two literals, and we are asked to decide whether or not the formula is satisfiable. Let (*G*, *H*) be an instance of 2-GRAPH EMBEDDING satisfying |*L*(*v*)| ≤ 2, for every *v * V (G). We show how to construct in polynomial time an instance *F *of 2-CNF-SATISFIABILITY such that *G *has a valid embedding into *H *if and only if *F *is satisfiable.

For every vertex *v * *G*: if $L\left(\nu \right)=\left\{{\nu}^{\text{'}}\right\}$ we add a variable ${x}_{\nu {\nu}^{\text{'}}}$ and add the clause $\left\{{x}_{\nu {\nu}^{\text{'}}}\right\}$ to *F*; and if $L\left(\nu \right)=\left\{{\nu}^{\text{'}},{\nu}^{\u2033}\right\}$ we add the two variables ${x}_{\nu {\nu}^{\text{'}}}$, ${x}_{\nu {\nu}^{\u2033}}$and the two clauses $\left\{{x}_{\nu {\nu}^{\text{\'}}},{x}_{\nu {\nu}^{"}}\right\}$, $\left\{{\overline{x}}_{\nu {\nu}^{\text{'}}},{\overline{x}}_{\nu {\nu}^{\u2033}}\right\}$ to *F*. This ensures that every vertex *v *in *G *is mapped to one and only one vertex in *H *(i.e., the map is a well-defined function). (We assume that |L(*v*)| ≠ 0; otherwise, the instance can be rejected.)

For every two vertices *v *and *u *in *G *such that there is a directed path from *v *to *u *in *G *(i.e., *v *appears before *u *in the directed path in *G*), and for very ${\nu}^{\text{'}}\in L\left(\nu \right)$ and ${u}^{\text{'}}\in L\left(u\right)$ such that ${\nu}^{\text{'}}={u}^{\text{'}}$ or ${u}^{\text{'}}$appears before *u *in the directed path in *H*, we add the clause $\left\{{\overline{x}}_{\nu {\nu}^{\text{'}}},{\overline{x}}_{u{u}^{\text{'}}}\right\}$ to *F*. This ensures that the desired mapping is injective, and ensures that the mapping respects the precedence order among the vertices in *G *that is defined by the directed path in *G *(property (ii)).

For every two vertices *v *and *u *in G such that *vu * *E*(*G*), and for very ${\nu}^{\text{'}}\in L\left(\nu \right)$ and ${u}^{\text{'}}\in L\left(u\right)$ such that ${\nu}^{\text{'}}{u}^{\text{'}}\notin E\left(H\right)$, we add the clause $\left\{{\overline{x}}_{\nu {\nu}^{\text{'}}},{\overline{x}}_{u{u}^{\text{'}}}\right\}$ to *F*. This ensures that the desired mapping respects the undirected edges of *G *(property (iii)).

This completes the construction of *F*. Clearly, this construction can be carried out in polynomial time.

It is not difficult to verify that (*G*, *H*) is a yes-instance of 2-GRAPH EMBEDDING if and only if *F *is a yes-instance of 2-CNF-SATISFIABILITY. This implies that 2-GRAPH EMBEDDING is polynomial-time solvable. □

The above theorem, together with the result in [3], provides a complete characterization of the complexity (NP-hardness) of *r*-GRAPH EMBEDDING with respect to *r*.

If we consider the *r*-GRAPH EMBEDDING parameterized by *r*, the fact that the problem is NP-complete for *r *≥ 3 [3] implies that the problem is not solvable in time *O*(*n ^{r}*) unless $\mathcal{P}=\mathcal{N}\mathcal{P}$, and hence, with respect to the parameterized complexity framework, the problem is not in the class $\mathcal{X}\mathcal{P}$. Therefore, there is not much hope behind seeking parameterized algorithms (with respect to

We investigate next the complexity of the *r*- GRAPH EMBEDDING_{≥ }problem.

**Theorem 0.2 ***The r*-GRAPH EMBEDDING_{≥ }*problem is *$\mathcal{N}\mathcal{P}-complete$, *for any r *≥ 1.

PROOF. It suffices to prove the $\mathcal{N}\mathcal{P}-\mathsf{\text{completeness}}$ of the 1-GRAPH EMBEDDING_{≥ }problem. We only prove the NP-hardness, as it is very easy to show the membership of the problem in $\mathcal{N}\mathcal{P}$. The proof is via a reduction from the CLIQUE problem: Given a graph and a nonnegative integer k, determine if the graph has a clique (complete subgraph) of size *k*.

Let $\left({G}^{\text{'}},k\right)$ be an instance of CLIQUE, where $V\left({G}^{\text{'}}\right)=\left\{{\nu}_{1}^{\text{'}},...,{\nu}_{n}^{\text{'}}\right\}$. We construct the instance (*G*, *H*, *k*) of 1- GRAPH EMBEDDING_{≥ }as follows. The set of vertices *V *(*G*) = {*v*_{1}, ... ,*v _{n}*} and

It is not difficult to verify that $\left({G}^{\text{'}},k\right)$ is a yes-instance of CLIQUE if and only if (*G*, *H*, *k*) is a yes-instance of 1- GRAPH EMBEDDING_{≥}. This completes the proof. □

The reduction in the above theorem is an fpt-reduction, from the CLIQUE problem to 1- GRAPH EMBEDDING_{≥}, where the parameter is the size of the subgraph sought *k*. Since CLIQUE is known to be *W *[1]-hard in the parameterized complexity hierarchy, we obtain:

**Theorem 0.3 **The *r*- GRAPH EMBEDDING_{≥ }*problem is W *[1]-*complete*, *for any r *≥ 1. *(Note that membership in W *[1]*follows from the results in the next section.*)

Finally, we observe that the same reduction in Theorem 0.2 provides an *L*-reduction [24] (i.e., approximation-preserving reduction) from MAXIMUM CLIQUE (the problem of computing a clique of maximum cardinality in a graph) to MAXIMUM 1-GRAPH EMBEDDING. It is well known that, unless $\mathcal{P}=\mathcal{N}\mathcal{P}$, MAXIMUM CLIQUE cannot be approximated to a ratio ${n}^{\frac{1}{2}-\epsilon}$ for any *ε *> 0 [25]. It follows that:

**Theorem 0.4 ***Unless *$\mathcal{P}=\mathcal{N}\mathcal{P}$, *the *MAXIMUM *r*-GRAPH EMBEDDING *problem cannot be approximated to a **ratio *${n}^{\frac{1}{2}-\epsilon}$*for any **ε *> 0.

In this section we show that the MAXIMUM *r*-GRAPH EMBEDDING problem can be modeled as an MAXIMUM INDEPENDENT SET problem. Recall that an *independent **set *in a graph is set of vertices such that no two of them are adjacent, and the MAXIMUM INDEPENDENT SET problem asks for an independent set of maximum cardinality in a graph.

Let (*G*, *H*) be an instance of MAXIMUM *r*-GRAPH EMBEDDING. Suppose that *V *(*G*) = {*g*_{1}, *g*_{2}, ... ,*g _{n}*} with directed edges from

**Theorem 0.5 ***If *MAXIMUM INDEPENDENT SET *is solvable in time *2* ^{cn}*,

PROOF. Create an auxiliary graph *X *as follows. For each possible choice mapping *g _{i }*to

1. *i *= *k *or *j *= *l*.

2. *i *<*k *and *j *>*l*, or *i *>*k *and *j *<*l*.

3. There is an undirected edge between *g _{i }*and

Note that Condition 2 could be removed when the order of the mapped vertices are not required to be the same for the two graphs.

It is clear that any independent set of *X *corresponds to a common subgraph of *G *and *H *of the same size. So the problem of finding a maximum common subgraph of *G *and *H *is reduced to the problem of finding a maximum independent set of *X*, which has *rn *vertices. In particular, to find if *G *is a subgraph of *H *it suffices to find an independent set of size *n*. Therefore if MAXIMUM INDEPENDENT SET is solvable in time 2* ^{cn}*, then MAXIMUM

If we use the current-best exact algorithm for MAXIMUM INDEPENDENT SET by Robson [26] that runs in time *O*(2^{n}^{/4}), we conclude that:

**Theorem 0.6 ***The *MAXIMUM *r*-GRAPH EMBEDDING *problem is solvable in time O*(2^{rn}^{/4}), *where n is the number of vertices in graph **G*.

The problem of protein structure comparison could be modeled as finding an independent set problem of an auxiliary graph. When aligning two protein structures, the auxiliary graph *X *is created as is in the proof of Theorem 2.5. Note that when aligning three and more protein structures, the auxiliary graph *X *could be created similarly.

Refer to the following for the outline of the algorithm for protein structure comparison.

1 *(Preprocessing)*. Generate the two structure graphs for the two proteins, based on both their secondary structure information (local structure) and tertiary structure (global structure) information.

2 *(Auxiliary graph)*. Build the auxiliary graph from the two structure graphs;

3 *(Top K independent sets)*. Generate the top *K *maximum independent sets of the auxiliary graph by applying the enumeration algorithm developed in [4].

4 *(Matched SSEs)*. Evaluate the generated top *K *maximum independent sets and generate the SSE pairs with the best score of the two proteins.

We analyze the time complexity of the algorithm:

Step 1: The algorithm processes the two proteins to generate the corresponding two structure graphs, where each vertex of a graph represents an SSE of the corresponding protein. Suppose the number of the vertices of each structure graph is bounded by *n*.

Step 2: We introduce a parameter *r *as the maximum number of pairs associated with each vertex of the structure graphs. The number of vertices of the auxiliary graph is bounded *rn*.

Step 3: Through calling the enumeration algorithm develop in [4], it takes time *O*(1.47* ^{rn}*) to generate the top

Step 4: It takes time *O*(1.47^{rn}n^{2}) to evaluate the generated independent sets and identify the independent set, which corresponding to the SSE pairs with the best score of the two proteins.

Refer to [27] for a discussion of the theoretical running times of several other graph-based approaches for protein structure comparison, which are of *O*((*mn*)* ^{n}*) or

Our approach ePC is designed for general-purpose protein structure comparison. In this section we test our approach for this purpose using SABmark-sup and SABmark-twi [28], and specific novel folds studied in the literature. Our approach is implemented using C++. The testing is mainly performed on a regular Macbook (8GB Mem). The running-time testing is conducted on a Dell server (PowerEdge 2950III, 32GB Mem). Due to the space limit, some testing results are not presented.

Given two proteins, *A *and *B*, the score of the a SSE pair is the sum of the *L _{ij }*of the residues for the SSE pair.

Let *S *be the sum of the scores of all the aligned SSEs. The normalized score ${S}_{n}=S/\sqrt{\left({l}_{A}*{l}_{B}\right)}$, where *l _{A }*and

There are two important parameters of our algorithm *r *and *K*, where *r *is the maximum number of SSE pairs associated with each SSE of the structure graphs, and *K *is the number of enumerated independent sets. Note that the score *L _{ij }*of the SSE pairs is the criteria for identifying the associated

We present our testing results for accuracy (using the score *S *as a criteria) and running-time of our approach with different parameter *r *values. We have conducted the testing of 200 protein pairs from SABmark-sup database with different parameter *r *values, where each SSE from one protein is matched with the top *r *SSEs from the other protein. Refer to Figure Figure22 for the average scores of 200 protein pairs from Sup database, when testing our approach with different parameter *r *values, *r *= 2, 3, 4, 5, 6, 7, 8, 9. Our testing results indicate that when the parameter *r *value increases, the score has increased. Refer to Figure Figure33 for the average running times of 200 protein pairs from Sup database, when testing our approach with different parameter *r *values. When *r *increases, the running time of our approach increases in general. However note that the running times when *r *= 5, 6, 7, 8, 9 are very similar; this is because trimming has been applied to reduce the sizes of the auxiliary graphes before the enumeration of the independent sets, and also because the impact of the parameter *K *on the running time. Especially the running time when *r *= 2 is significantly lower than the other cases, which matches our theoretical result that for *r *= 2 the *r*-GRAPH EMBEDDING problem is in P.

For the enumeration of independent sets, we have introduced a parameter *K*, which is the bound of the number of enumerated independent sets. Here we present our testing results for accuracy and running-time of our approach with different parameter K values (See Table Table1).1). Similar as the testing for the parameter *r*, we have conducted the testing of 200 protein pairs from SABmark-sup database with different parameter *K *values, *K *= 125, 250, 500, 1000. Our testing results indicate that when the parameter *K *value increases, the score has increased and the running time also increases.

*Self-querying in a large database of structures*. As pointed in [5], a necessary condition for a approach to be of practical value for structure comparison and classification, it should be able to find the query itself in a database of protein structures. To test this property of our approach, 1000 protein structures from the SABmark-sup database. Our approach with the normalized score function can identify the query structure with ranking No. 1 with 100% accuracy.

*Detecting a substructure in a set of larger structures*. Our approach can detect a smaller query structure (or, a motif structure) within a larger target structure. We use the test set from the previous test and required for each domain to be matched to the target domain embedded in the original full-protein structure. Our approach with the normalized score can identify the substructure with ranking No. 1 with 100% accuracy.

*Protein family classification*. We compare the performance of our approach for protein family classification with deconSTRUCT, which is also an SSE-based method and designed for protein structure database filtering. We have tested 1000 proteins pairs of the SABmark [28]. Due to the space limit, we only discuss some of the representative testing result. We align protein d1a6m (core size:7; AAs: 151; from SABmark [28]) to proteins from 10 different families of the twi database with each family 10 proteins. Of the proteins in the top 10 ranking, 7 proteins identified through our approach are the proteins from the same family as protein d1a6m. For deconSTRUCT, 7 proteins of the identified 10 proteins (without ranking) are the proteins from the same family as protein d1a6m. Form the testing results, our approach has comparable performance with DeconSTRUCT for the general purpose protein structure comparison and structure classification. The mixed graph representation of our approach ePC is much simpler compared with deconSTRUCT. Our approach ePC is more flexible than deconSTRUCT in that ePC can handle SSE alignments with and without respect to the order of SSEs, which will be discussed in the next section for specific examples.

We test our approach on specific examples for common substructures and novel folds which share common substructures with non-sequential SSEs.

*Detection of several different common substructures*. We test our approach ePC using the four protein structures (PDB codes: 1a02N, 1iknA, 1nfiA, and 1a3qA) studied in [8,9]. The proteins share two common domains: "p53-like transcription factors" and "E set domains". In [8,9] two different common substructures were detected, one for each domain. The first common substructure is part of the "p53-like transcription factors" domain. It consists of 114 residues, and it forms a sandwich of nine *beta*-strands. The second common substructure is part of the "E set domains" domain. It consists of 87 residues, and it forms a sandwich of seven *beta*-strands.

Please refer to the following testing results of our approaches, when 1a02N is compared with: 1iknA, 1nfiA, and 1a3qA. Our testing results match the results in [8,9]. Especially for the second common substructure that is part of the "E set domains" domain with conserved matched SSEs: 12, 13, 14, 15, 16, 17 of 1a02N. Please refer to Figure Figure44 for its 3D structure and the two domains.

*Three novel folds*. The three novel folds were discussed in [7] to study the unique feature of GANGSTA+ to conduct non-sequential SSE alignment. Note that the protein structures that are structurally similar to the listed three new folds were detected through scanning the ASTRAL40 database by GANGSTA+. The detected similar protein structures have non-sequential SSE alignments with the three novel folds respectively. Please refer to our testing result in Table Table22 Figure Figure55 and and66.

We use an SSE-based graph model for general purpose protein structure comparison. We presented the computational complexity results related to the protein structure comparison problem. An effective algorithm is developed integrating a novel enumeration of independent sets and parameterized computation for the problem. Our approach is tested for protein structure comparison using benchmark testing sets. Compared with other SSE-based approaches, our approach has comparable performance for the general purpose protein structure comparison. We also demonstrate that our approach could be applied to identify common substructure with non-sequential SSEs and proteins sharing more than one common substructure.

The authors declare that they have no competing interests.

XH, IK and GX carried out the study on the complexity and the design of the approach for the protein structure comparison problem, and drafted the manuscript. CA, DJ and KW participated in the implementation and the testing of the algorithm. All authors have approved the final manuscript.

This research is supported by the National Institute of Health grants from the National Center for Research Resources (5P20RR016460-11) and the National Institute of General Medical Sciences (8P20GM103429-11).

The publication costs for this article were funded by the corresponding author's institution.

This article has been published as part of *BMC Genomics *Volume 14 Supplement 2, 2013: Selected articles from ISCB-Asia 2012. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/14/S2.

- Holm L, Sander C. Protein structure comparison by alignment of distance matrices. J of Molecular Biology. 1993;233:123–138. doi: 10.1006/jmbi.1993.1489. [PubMed] [Cross Ref]
- Goldman D, Istrail S, Papadimitriou CH. Algorithmic Aspects of Protein Structure Similarity. FOCS. 1999. pp. 512–522.
- Song Y, Liu C, Huang X, Malmberg RL, Xu Y, Cai L. Efficient parameterized algorithms for biopolymer structuresequence alignment. IEEE/ACM Trans Comput Biology Bioinform. 2006;3(4):423–432. [PubMed]
- Chen J, Kanj I, Meng J, Xia G, Zhang F. On the effective enumerability of NP problems. Proceedings of the 2nd InternationalWorkshop on Parameterized and Exact Computation, volume 4169 of Lecture Notes in Computer Science. 2006. pp. 215–226.
- Zhang ZH, Bharatham K, Sherman WA, Mihalek I. deconSTRUCT: general purpose protein database search on the substructure level. Nucleic Acids Research. 2010;38(Web Server):W590–W594. doi: 10.1093/nar/gkq489. [PMC free article] [PubMed] [Cross Ref]
- Krissinel E, Henrick K. Secondary-structure matching (PDBeFold), a new tool for fast protein structure alignment in three dimensions. Acta Cryst D60. 2004. pp. 2256–2268. [PubMed]
- Guerler, Knapp. Novel Folds and their Nonsequential Structural Analogs. Protein Science. 2008;17:8:1374–1382. [PubMed]
- Dror O, Benyamini H, Nussinov R, Wolfson H. MASS: Multiple structural alignment by secondary structures. Bioinformatics. 2003;19(Suppl 1):i95–i104. doi: 10.1093/bioinformatics/btg1012. [PubMed] [Cross Ref]
- Dror O, Benyamini H, Nussinov R, Wolfson H. Multiple structural alignment by secondary structures: algorithm and applications. Protein Science. 2003;12:2492–2507. [PubMed]
- Gibrat JF, Madej T, Bryant SH. Surprising similarities in structure comparison. Curr Opin Struct Biol. 1996;6(3):377–385. doi: 10.1016/S0959-440X(96)80058-3. [PubMed] [Cross Ref]
- Michalopoulos I, Torrance GM, Gilbert DR, Westhead DR. TOPS: an enhanced database of protein structural topology. Nucleic Acids Research. 2004;32:251–254. doi: 10.1093/nar/gkh060. [PMC free article] [PubMed] [Cross Ref]
- Alesker V, Nussinov R, Wolfson H. Detection of non-topological motifs in protein structures. Protein Eng. 1996;9:1103–1119. doi: 10.1093/protein/9.12.1103. [PubMed] [Cross Ref]
- Alexandrov N, Fischer D. Analysis of topological and nontopological structural similarities in the PDB: New examples with old structures. Proteins. 1996;25:354–365. doi: 10.1002/(SICI)1097-0134(199607)25:3<354::AID-PROT7>3.3.CO;2-W. [PubMed] [Cross Ref]
- Grindley H, Artymiuk P, Rice D, Willett P. Identification of tertiary structure resemblance in proteins using a maximal common subgraph isomorphism algorithm. J Mol Biol. 1993;229:707–721. doi: 10.1006/jmbi.1993.1074. [PubMed] [Cross Ref]
- Holm L, Sander C. 3-D lookup: Fast protein structure database searches at 90% reliability. The Third International Conference on Intelligent Systems for Molecular Biology. 1995. pp. 179–187. [PubMed]
- Koch I, Lengauer T, Wanke E. An algorithm for finding maximal common subtopologies in a set of proteins. J Comp Biol. 1996;3:289–306. doi: 10.1089/cmb.1996.3.289. [PubMed] [Cross Ref]
- Lu G. TOP: A new method for protein structure comparisons and similarity searches. J Appl Crystallogr. 2000;33:176–183. doi: 10.1107/S0021889899012339. [Cross Ref]
- Mitchel E, Artymiuk P, Rice D, Willet P. Use of techniques derived from graph theory to compare secondary structure motifs in proteins. J Mol Biol. 1990;212:151–166. doi: 10.1016/0022-2836(90)90312-A. [PubMed] [Cross Ref]
- Yang AS, Honig B. An integrated approach to the analysis and modeling of protein sequences and structures. I. Protein structural alignment and a quantitative measure for protein structural distance. J Mol Biol. 2000;301:65–678. [PubMed]
- Joosten RP, Te Beek TAH, Krieger E, Hekkelman ML, Hooft RWW, Schneider R, Sander C, Vriend G. A series of PDB related databases for everyday needs. NAR. 2010. doi: 10.1093/nar/gkq1105. [PMC free article] [PubMed]
- Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983;22:2577–2637. doi: 10.1002/bip.360221211. [PubMed] [Cross Ref]
- Papadimitriou CH. Computational Complexity. Addison-Wesley; 1994.
- Impagliazzo R, Paturi R, Zane F. Which problems have strongly exponential complexity? Journal of Computer and System Sciences. 2001;63(4):512–530. doi: 10.1006/jcss.2001.1774. [Cross Ref]
- Papadimitriou CH, Yannakakis M. Optimization, approximation, and complexity classes. J Comput Syst Sci. 1991;43(3):425–440. doi: 10.1016/0022-0000(91)90023-X. [Cross Ref]
- Håstad Johan. Clique is Hard to Approximate Within n
^{1-epsilon}. Proceedings of the 37th Annual Symposium on Foundations of Computer Science. 1996. pp. 627–636. - Robson JM. Technical Report. LaBRI, Universite Bordeaux I; Finding a maximum independent set in time
*O*(2/^{n}), 2001; pp. 1251–01.^{4} - Krissinel E, Henrick K. Protein structure comparison service Fold at European Bioinformatics Institute. http://www.ebi.ac.uk/msd-srv/ssm
- Van Walle I. et al. SABmark: a benchmark for sequence alignment that covers the entire known fold space. Bioinformatics. 2005;21:1267–1268. doi: 10.1093/bioinformatics/bth493. [PubMed] [Cross Ref]
- Zhu J, Weng Z. FAST: a novel protein structure alignment algorithm. Proteins. 2005;58(3):618–627. [PubMed]

Articles from BMC Genomics are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |