Finding all longest common segments in protein structures efficiently

Research output: Journal Publications and Reviews (RGC: 21, 22, 62)21_Publication in refereed journalpeer-review

2 Scopus Citations
View graph of relations


Related Research Unit(s)


Original languageEnglish
Pages (from-to)644-655
Journal / PublicationIEEE/ACM Transactions on Computational Biology and Bioinformatics
Issue number3
Publication statusPublished - 1 May 2015


The Local/Global Alignment (Zemla, 2003), or LGA, is a popular method for the comparison of protein structures. One of the two components of LGA requires us to compute the longest common contiguous segments between two protein structures. That is, given two structures A = (a1;⋯;an) and B = (b1;⋯;bn) where ak, bk ∈ ℝ3, we are to find, among all the segments f = (ai;⋯; aj) and g = (bi;⋯;bj) that fulfill a certain criterion regarding their similarity, those of the maximum length. We consider the following criteria: (1) the root mean squared deviation (RMSD) between f and g is to be within a given teR; (2) f and g can be superposed such that for each k,i ≤ k ≤ j, kak - bkk ≤t for a given t ∈ ℝ. We give an algorithm of O(n log n + nl) time complexity when the first requirement applies, where l is the maximum length of the segments fulfilling the criterion. We show an FPTAS which, for any t ∈ ℝ, finds a segment of length at least 1, but of RMSD up to (1 + ∈)t, in O(n log n + n=∈) time. We propose an FPTAS which for any given ∈ ∈ ℝ, finds all the segments f and g of the maximum length which can be superposed such that for each k,i ≤k ≤ j, ||ak - bkk <(1 + ∈)t, thus fulfilling the second requirement approximately. The algorithm has a time complexity of O(n log2 n/∈5) when consecutive points in A are separated by the same distance (which is the case with protein structures). These worst-case runtime complexities are verified using C++ implementations of the algorithms, which we have made available at

Research Area(s)

  • Local/Global Alignment (LGA), Longest common structure, PTAS, RMSD

Citation Format(s)