- Methodology article
- Open Access
A global optimization algorithm for protein surface alignment
© Bertolazzi et al; licensee BioMed Central Ltd. 2010
Received: 7 October 2009
Accepted: 29 September 2010
Published: 29 September 2010
A relevant problem in drug design is the comparison and recognition of protein binding sites. Binding sites recognition is generally based on geometry often combined with physico-chemical properties of the site since the conformation, size and chemical composition of the protein surface are all relevant for the interaction with a specific ligand. Several matching strategies have been designed for the recognition of protein-ligand binding sites and of protein-protein interfaces but the problem cannot be considered solved.
In this paper we propose a new method for local structural alignment of protein surfaces based on continuous global optimization techniques. Given the three-dimensional structures of two proteins, the method finds the isometric transformation (rotation plus translation) that best superimposes active regions of two structures. We draw our inspiration from the well-known Iterative Closest Point (ICP) method for three-dimensional (3D) shapes registration. Our main contribution is in the adoption of a controlled random search as a more efficient global optimization approach along with a new dissimilarity measure. The reported computational experience and comparison show viability of the proposed approach.
Our method performs well to detect similarity in binding sites when this in fact exists. In the future we plan to do a more comprehensive evaluation of the method by considering large datasets of non-redundant proteins and applying a clustering technique to the results of all comparisons to classify binding sites.
The function of a protein typically depends on the structure of specific binding sites located at the surface of the protein where the interaction with a ligand takes place. The identification of protein binding sites, their classification and analysis is of much interest for drug design and treatment of diseases. Binding sites recognition is generally based on geometry often combined with physico-chemical properties of the site since the conformation, size and chemical composition of the protein surface are all relevant for the interaction with a specific ligand.
In this paper we address the problem of optimally aligning protein surfaces, i.e. of finding atom pairs on two protein surfaces that occupy spatially equivalent positions. Our computational method integrates geometry with chemical properties of the matched atoms. It can be applied to the comparison of binding sites as well as of any other surface patches, such as cavities, that may be of interest.
Although the literature in protein surface alignment is not as vast as the one on complete structure or fold alignment, nevertheless several matching strategies have been designed for the recognition of protein-ligand binding sites and of protein-protein interfaces. They include hashing techniques [1, 2], graph theoretic methods [3–6], descriptors based on moments  and moment invariants , shape descriptors such as spin images [9–11]. A few web servers have recently become available [12–16].
Most of the proposed methods require the solution to a 3D matching problem which is a well-studied problem also in computer vision and robotics. It can be formulated as follows: given two sets A and B of points, find two possibly large subsets A' of A and B' of B with high degree of similarity. There are various ways of defining the similarity between two point sets in 3D space leading to the proposal of different distance functions and associated algorithms; they include the root mean square distance, the closest point distance , the Hausdorff distance , the bottleneck distance .
An important aspect of the matching is the choice of a suitable surface representation; in the literature common ways of representing a surface are Connolly's representation , alpha-shapes  and pseudo-vertices . In our approach, we represent the surface as a cloud of points, each corresponding to a surface atom. Thus, the protein surface alignment problem is the same as the aforementioned 3D matching problem.
One possible way to solve the surface alignment problem is by using the well-known Iterative Closest Point (ICP) algorithm  from which we draw our inspiration. The ICP algorithm, originally introduced in the area of computer vision for image registration, has been used in bioinformatics  for the alignment of complete protein structures. Indeed, we take a similar approach for surface alignment, namely we search for the isometric transformation which best superimposes two given protein structures.
Our main contribution is in the adoption of a different, more efficient global optimization approach along with a new dissimilarity measure. The global optimization algorithm we design belongs to the class of controlled random search methods [23–25]. These methods, although heuristic in nature, are very efficient and reliable for the global minimization of nonlinear multivariate functions of several variables. In the past years, controlled random search algorithms have been successfully used to solve many real world problems, see for instance [26–31]. The dissimilarity measure we propose is based on the solution to an "Asymmetric Assignment Problem" on a bipartite graph associated to the matching problem. Our method is capable of generating very accurate local alignments. We benchmark it on various sets of protein structures from the PDB , and compare its performance with that MolLoc .
Notations and Assumptions
and assume, without loss of generality, that n ≤ m. The set P is conventionally representative of a query shape while Q defines a reference model shape.
Let: ψ: P → Q denote a point to point mapping that associates to every point of P a point of Q. Since, as assumed above, P and Q are finite sets, the class Ψ of all mappings ψ has finite cardinality given by |Ψ| = m n .
for all ψ ∈ Ψ and a ∈ Θ.
The problem can also be formulated in terms of the following definition of an assignment. A function φ ∈ Ψ is an assignment from P to Q if, by definition, it is injective, that is for every p1, p2 ∈ P, p1 ≠ p2 implies φ(p1) ≠ φ (p2).
Let us denote by Φ ⊆ Ψ the class of all possible assignments from P to Q. Obviously, since P and Q are finite sets, Φ is finite as well and its cardinality is |Φ| = m (m - 1) ... (m - n + 1).
As it is stated in , where ICP has been originally proposed, the method converges to a solution which is a local minimum of the two-level Problem (2). Futher, in  it has been shown that the final transformation and mapping obtained by Algorithm ICP heavily depend on the initial relative positioning of sets P and Q.
In this section we discuss the use of a continuous global optimization algorithm for the solution of the shape alignment problem. To this aim, it is necessary to reformulate the shape alignment problem in a complementary way with respect to Problem (2). More in particular, the inner-level problem becomes the one defining the mapping function ψ (instead of the transformation a) once the transformation vector a ∈ Θ is fixed in the outer level.
where g(a) = f(ψ(a), a). Every global solution a* of (5) is, by definition, a solution such that f (ψ (a*), a*) ≤ f (ψ (a), a), for all a ∈ Θ.
Observe that the computation of function g requires the computation of the optimal mapping ψ (a), that is, the solution to Problem (4). This latter problem can be solved with a time complexity O(nm) in the worst case  which can be relevant for n and m large. Moreover, due to its definition, g(a) is a non-smooth (Lipschitz) continuous function and its derivatives are not available. Indeed, for the minimization of function g(a) we can neither directly use its derivatives nor approximate them through finite differences since this would require too much time and produce numerical derivatives which are unreliable because of the non-smoothness of function g.
initialization Let N = 6 and choose an integer M ≫ N. The objective function is sampled on a set S of M points a ∈ ℜ7 randomly chosen within the feasible domain Θ strictly containing the global minimizer.
- 2.stopping criterion If the maximum and minimum values of the objective function over S are sufficiently close to each other, namely
search phase N +1 points are randomly chosen in the set S. Then,
the weighted centroid a c of the N + 1 points is computed;
the new trial point ã is computed by doing a weighted reflection of the centroid onto the worst point among the selected N + 1 points.
The normalization matrix W is introduced to ensure that the first four components of the resulting point represent a unit quaternion (i.e. a rotation).
updating phase If the objective function value on the new point ã improves on the maximum function value over S, then the set S is updated by adding the new point and discarding the worst one. Otherwise the set S is left unchanged and the new point is discarded. The algorithm continues iterating through steps 2-4.
The algorithm starts by randomly choosing M ≫ N points over the feasible set Θ. In the literature, a typically accepted value of M is 25N[25, 29]. This value is able to convey to the algorithm sufficient ability to find the global minimum point without excessively slowing down the convergence.
A new dissimilarity measure
In this section we propose a new dissimilarity measure between two given sets of points of two proteins. This measure is based on a distance other than the closest point distance.
In particular, it can be noted that, using the closest point distance, it can happen that different points of set T a (P) are mapped to the same point of set Q. This, in turn can yield a distance value which is small just because many points are all mapped to the same closest point.
Problem (6) can be formulated as a 0, 1-optimization problem and is, indeed, the combinatorial optimization problem known as the Asymmetric Assignment Problem (AAP).
where, δ+ (p) and δ- (q) are the sets of edges leaving node p and, respectively, entering in node q.
Note that the constraints of Problem (7) require each node p ∈ P to be assigned to exactly one node q ∈ Q and each node q ∈ Q to be assigned to at most one node p ∈ P, which is why Problem (7) is known as Asymmetric Assignment Problem.
where s* is the optimal solution to Problem (7).
Problem (7), and hence (6), can be solved very efficiently by ad-hoc codes that have time complexity where see for instance .
We are now able to define our new dissimilarity mea-sure, that we call matching distance.
Definition 1 Given an isometric transformation a ∈ Θ and two distinct sets of points P and Q, the matching distance between T a (P ) and Q is given by f(φ(a), a).
Reasoning as in the preceding section, we can now search for a global solution to problem (5) where now g(a) = f(φ(a), a).
Integration of physico-chemical properties
Up to this point, the discussed approach is based on geometry only. However, as is well known in biology, there are other properties that affect the binding of molecules. For instance, electrostatic as well as hydrophobic-hydrophilic properties play an important role in protein-protein and protein-ligand interactions. Thus, we consider a variant of our approach in which we integrate physico-chemical properties. Specifically in the graph G(P, Q, E) we assume that the edge e = (p, q) is present only if the two atoms p ∈ P and q ∈ Q have the same physico-chemical properties. According to , we say that p and q have the same physico-chemical properties if they are both Acceptor (ACC), Donor (DO), Acceptor/Donor (AD), Aliphatic (ALI) or Aromatic (PI). Furthermore, we assume that, for every p ∈ P at least a node q ∈ Q exists such that (p, q) ∈ E.
We applied our method, referred to as Continuous Optimization (CO) method in the following, to the comparison of binding sites of proteins. We integrated physico-chemical properties in our method, as discussed in the previous section. The structures of the proteins in complex with specific ligands are taken from the PDB . The binding sites are extracted by a simple algorithm that finds all protein atoms within a certain distance (4.0 Å) from an atom of the ligand. We run our algorithm on pairs of binding sites producing in output the list of matched atoms on the two binding sites, the rigid transformation that best superimposes them, and the RMSD after superposition.
We benchmarked CO on a dataset of 100 proteins in complex with 9 ligands that differ in chemical composition as well as in size and shape. The results of all-to-all pairwise comparisons are visualized by means of a distance matrix and by the ROC curves. The goal is to evaluate the ability of CO in assigning a binding site to the correct group of proteins, i.e. those binding the same ligand.
We then present more detailed results on a set of 19 binding sites of proteins in complex with the ligand ATP with the goal of judging the quality of the alignments. For each comparison we report the number of aligned atoms as well as their RMSD after superposition. The results on this dataset are compared with those of another method, MolLOC , which derives the same two measures, i.e. number of aligned atoms and RMSD.
Classification of proteins according to their bound ligand
In the first experiment we perform all-to-all comparisons on a dataset of 100 proteins in complex with one of 9 ligands: AMP, ATP, FAD, FMN, GLC, HEME, NAD, PO4, and Steroid. This dataset was used in  for an analysis of shape variation in protein binding sites. The proteins were carefully selected, with a number of criteria, so that the dataset is non-redundant and the binding sites are not evolutionary related.
where n and m are the numbers of atoms of the two binding sites. The proteins are listed along the rows and columns of the matrix so that proteins binding the same ligand are grouped together. Horizontal and vertical black lines on the matrix separate different groups of proteins. The matrix is color-coded from 0 to 1, with red corresponding to high number of aligned atoms and therefore high similarity in the shape of the binding sites and blue to the lowest degree of similarity. A good classification of sites based on bound ligands implies the presence of mostly red areas around the main diagonal, corresponding to pairwise comparisons within the same group of proteins, i.e. in complex with one specific ligand. This can be in fact observed in the image matrix although with different degrees for the different groups of proteins. As it is known , ligand PO4 tends to be rigid, exhibiting little conformational variability in the binding. Not surprisingly, the corresponding area is the one showing the highest degree of similarity. The method CO appears to perform well also in distinguishing the PO4 group from any other group, as PO4 binding sites are more similar to themselves than to binding sites of other groups. Similar considerations apply to steroid and GLC. A good performance is also obtained for the HEME group, although the discriminating power with the NAD group is not clear. As noted in , ligand ATP has great variation in its conformation when binding different proteins: it can be in an extended conformation or in a compact one, resulting in different sizes and shapes of the binding regions. This is reflected in our experiments, as can be seen from the distance matrix where blue or green areas are present.
From the above sets of experiments we can conclude that CO has a good accuracy in the retrieval of similarity information: for a given query binding site the highest scoring solutions are generally the binding sites of the dataset in complex with the same ligand as the query. Furthermore, when a good similarity in the binding is expected because of the relative rigidity of the ligands, CO is able to capture such a similarity, as shown in the distance matrix.
Comparing CO with other alignment methods on ATP binding sites
Several studies have been conducted to evaluate and compare different methods for determining the structural similarity of proteins. For instance, a comprehensive assessment of structural alignment methods is presented in  where six publicly available programs are evaluated on almost 9 million pairs of proteins. However, a similar large-scale experiment is not available for the related problems of aligning protein surfaces and binding sites, despite the growing number of methods and web servers available. There are several factors that contribute to the difficulty of the comparison. First, different methods solve different instances of the matching problem: some methods compare binding sites, while others recognize binding sites in cavities or even entire surfaces. Second, the methods differ in the input representations and scoring functions. For instance, in CO the input points are the atom centers, in Multibind a reduced set of points, the pseudo-centers. In  the points are the spherical sample points derived from the atomic coordinates. MolLoc, on the other hand, uses Connolly's  points and a richer surface representation based on local shape descriptors of surface points. As for the scoring function, although most methods produce the RMSD of the superimposed structures, some methods have a different native scoring function that cannot be easily derived by other methods.
As a comprehensive evaluation of all the techniques is beyond the scope of this paper, only MOlLoc will be considered in comparison with CO. The reason for choosing MolLoc is that both methods judge the quality of the alignments by the number of aligned atoms and their RMSD after superposition. Such measures are available at MolLOC website. As Multibind does not report the RMSD of two aligned structures at its website it will not be considered here. Moreover, the method in , based on spherical harmonics and benchmarked on the same dataset of 100 proteins, is not used in our evaluation because it computes a measure of similarity of two shapes without an alignment.
As observed in , although the ROC curves are a valid tool for assessing the quality of a classification approach they are often of limited value in comparing different methods; in fact such curves take into account only the ranking of the alignments not their quality. For this reason, since we want to assess the quality of the alignments we choose the geometric measure SAS [38, 40, 41]. Clearly, a better match has a higher number of aligned atoms and smaller RMSD. Since the two measures are not independent SAS combines them into a single expression: SAS = (RMSD × 100)/(num. aligned atoms).
We run both programs on the set of 19 proteins used in  for a related although different problem, that is binding site recognition within a cavity. The proteins all bind ligand ATP and are from different families according to the structural classification SCOP .
Comparison of CO with MolLoc.
N. corresp atoms
N. corresp atoms
We conclude this section by reporting that the execution times of CO on average on all 18 pair-wise comparisons considered in this experiment was 14.06 s for a total of 253.1 s on an Intel Pentium IV processor running at 2.66Ghz with 1Gb main memory. As we mentioned before, the low computational complexity of our proposed approach is one of the key points of our design. We do not report the execution times of MolLoc since they are not available from the web server interface.
Discussion and Conclusions
The main challenge for a method that compares and classifies binding sites is to be able to cluster the binding sites in groups according to the type of ligands they bind while at the same time allowing some conformational variability within the same group, as is often observed for binding sites of different proteins complexed with the same ligand. The difficulty arises because of the variety of ways in which a ligand can bind proteins. Although we expect a computational method to be able to distinguish among different types of ligands relatively well, there are obviously cases when only experimental methods can determine the binding affinity of two molecules.
Our proposed method, CO, performs well to detect similarity in binding sites when this in fact exists. In the future we plan to do a more comprehensive evaluation of the method by considering large datasets of non-redundant proteins and applying a clustering technique to the results of all comparisons to classify binding sites. A systematic evaluation of CO with other existing methods will be done through the introduction of a common scoring function that will overcome the problem that the available methods use native scoring functions difficult to export to other methods.
The authors are indebted to two anonymous Referees whose many helpful comments and suggestions greatly helped improving the paper. The work of C. Guerra was supported by the "Progetto di Ateneo, University of Padova" and by "Fondazione Cariparo"
- Shatsky M, Shulman-Peleg A, Nussinov R, Wolfson HJ: The multiple common point set problem and its application to molecule binding pattern detection. Journal of Computational Biology 2006, 13: 407–428. 10.1089/cmb.2006.13.407View ArticlePubMedGoogle Scholar
- Shulman-Peleg A, Nussinov R, Wolfson HJ: Recognition of functional sites in protein structures. Journal of Molecular Biology 2004, 339: 607–633. 10.1016/j.jmb.2004.04.012View ArticlePubMedGoogle Scholar
- Artymiuk P, Spriggs R, Willett P: Graph theoretic methods for the analysis of structural relationships in biological macromolecules. Journal of the American Society for Information Science and Technology 2005, 56(5):518–528. 10.1002/asi.20140View ArticleGoogle Scholar
- Chen B, Bryant D, Fofanov V, Kristensen D, Cruess A, Kimmel M, Lichtarge O, Kavraki L: Cavity-aware motifs reduce false positives in protein function pre-diction. Computational System Bioinformatics Conf 2005, 311–323.Google Scholar
- Hofbauer C, Lohninger H, Aszodi A: Surfcomp: A novel graph-based approach to molecular surface comparison. Journal of Chemical Information and Computer Sciences 2004, 44(3):837–847.PubMedGoogle Scholar
- Weskamp N, Kuhn D, Hullermeier E, Klebe G: Efficient similarity search in protein structure databases by k-clique hashing. Bioinformatics 2004, 20(10):1522–1526. 10.1093/bioinformatics/bth113View ArticlePubMedGoogle Scholar
- Ballester P, Richards W: Ultrafast shape recognition to search compound databases for similar molecular shapes. Journal of Computational Chemistry 2007, 28(10):1711–1723. 10.1002/jcc.20681View ArticlePubMedGoogle Scholar
- Sommer I, Müller O, Domingues F, Sander O, Weickert J, Lengauer T: Moment invariants as shape recognition technique for comparing protein binding sites. Bioinformatics 2007, 23: 3139–3146. 10.1093/bioinformatics/btm503View ArticlePubMedGoogle Scholar
- Bock M, Garutti C, Guerra C: Discovery of similar regions on protein surfaces. Journal of Computational Biology 2007, 14(3):285–299. 10.1089/cmb.2006.0145View ArticlePubMedGoogle Scholar
- Bock M, Garutti C, Guerra C: Effective labeling of molecular surface points for cavity detection and location of putative binding sites. Proc of the VI Int Conf on Computational Systems Bioinformatics; San Diego 2007, 263–274.View ArticleGoogle Scholar
- Bock M, Garutti C, Guerra C: Cavity detection and matching for binding site recognition. Theoretical Computer Science 2008, 408(2–3):151–162. 10.1016/j.tcs.2008.08.018View ArticleGoogle Scholar
- Angaran S, Bock M, Garutti C, Guerra C: MolLoc: a web tool for the local structural alignment of molecular surfaces. Nucleic Acids Research 2009., (Web server):Google Scholar
- Ausiello G, Gherardini PF, Marcatili P, Tramontano A, Via A, Helmer-Citterich M: Funclust: a web server for the identification of structural motifs in a set of non-homologous protein structures. BMC Bioinformatics 2008, 9(Suppl 2):S2. 10.1186/1471-2105-9-S2-S2View ArticlePubMedPubMed CentralGoogle Scholar
- Jambon M, Olivier A, Combet C, Deleage G, Delfaud F, Geourjon C: The SuMo server: 3D search for protein functional sites. Bioinformatics 2005, 21(20):3929–3930. 10.1093/bioinformatics/bti645View ArticlePubMedGoogle Scholar
- Kinoshita N, Furui J, Nakamura H: Identification of protein functions from a molecular surface database, ef-site. Journal of Structural and Functional Genomics 2001, 2: 9–22. 10.1023/A:1011318527094View ArticleGoogle Scholar
- Shulman-Peleg A, Shatsky M, Nussinov R, Wolfson HJ: MultiBind and MAPPIS: webservers for multiple alignment of protein 3D-binding sites and their interactions. Nucleic Acids Research 2008., (36 Web server):Google Scholar
- Besl PJ, McKay ND: A method for registration of 3-D shapes. IEEE Trans. on Pattern Analysis and Mach. Intelligence 1992, 14: 239–255. 10.1109/34.121791View ArticleGoogle Scholar
- Goodrich MT, Mitchell JSB, Orletsky MW: Practical methods for approximate geometric pettern matching under rigid motions. In Proc of the 10th Ann Symp on Computational Geometry Edited by: Press A. 1994, 103–112.Google Scholar
- Efrat A, Itai A, Katz MJ: Geometry helps in bottle-neck matching and related problems. Algorithmica 2001, 31: 1–28. 10.1007/s00453-001-0016-8View ArticleGoogle Scholar
- Connolly ML: Analytical molecular surface calculation. Journal of Applied Crystallography 1983, 16: 548–558. 10.1107/S0021889883010985View ArticleGoogle Scholar
- Liang J, Edelsbrunner H, Woodward C: Anatomy of protein pockets and caviteis: measurement of binding site geometry and implications for ligand design. Protein Science 1998, 7: 1884–1897. 10.1002/pro.5560070905View ArticlePubMedPubMed CentralGoogle Scholar
- Xu D, Li H, Gu T: Protein Structure Superposition by Curve Moment Invariants and Iterative Closest Point. The 1st Int Conf on Bioinformatics and Biomedical Engineering 2007, 6: 25–28. full_textGoogle Scholar
- Brachetti P, De Felice Ciccoli M, Di Pillo G, Lucidi S: A new version of the Price's algorithm for global optimization. Journal of Global Optimization 1997, 10: 165–184. 10.1023/A:1008250020656View ArticleGoogle Scholar
- Cirio L, Lucidi S, Parasiliti F, Villani M: A global optimization approach for the synchronous motors design by finite element analysis. Journal of Applied Electromagnetics and Mechanics 2002, 16: 13–27.Google Scholar
- Price WL: A controlled random search procedure for global optimization. In Towards Global Optimization 2. Edited by: Dixon L, Szego G. Amsterdam: North-Holland; 1978.Google Scholar
- Daidone A, Parasiliti F, Villani M, Lucidi S: A new method for the design optimization of three-phase induction motors. IEEE Trans. on Magnetics 1998, 34: 2932–2935. 10.1109/20.717684View ArticleGoogle Scholar
- Liuzzi G, Lucidi S, Piccialli V, Sotgiu A: A magnetic resonance device designed via global optimization techniques. Mathematical Programming 2004, 101: 339–364. 10.1007/s10107-004-0528-5View ArticleGoogle Scholar
- Liuzzi G, Lucidi S, Parasiliti F, Villani M: Multi-objective Optimization Techniques for the Design of Induction Motors. IEEE Transactions on Magnetics 2003, 39: 1261–1264. 10.1109/TMAG.2003.810193View ArticleGoogle Scholar
- Price WL: Global optimization by controlled random search. Journal of Optimization Theory and Applications 1983, 40: 333–348. 10.1007/BF00933504View ArticleGoogle Scholar
- Price WL: Global optimization algorithms for a CAD workstation. Journal of Optimization Theory and Applications 1983, 55: 133–146. 10.1007/BF00939049View ArticleGoogle Scholar
- Price WL, Woodhams F: Optimising accelerator for CAD workstations. IEEE proceedings 1988, 135: 214–221.Google Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Research 2000, 28: 235–242. 10.1093/nar/28.1.235View ArticlePubMedPubMed CentralGoogle Scholar
- Horn BKP: Closed-form solution of absolute orientation using unit quaternions. Journal of the Optical Society of America 1987, 4: 629–642.View ArticleGoogle Scholar
- Goldberg AV, Kennedy R: An Efficient Cost Scaling Algorithm for the Assignment Problem. Mathematical Programming 1995, 71: 153–178.Google Scholar
- Schmitt S, Kuhn D, Klebe G: A New Method to Detect Related Function Among Proteins Independent of Sequence and Fold Homology. Journal of Molecular Biology 2002, 323: 387–406. 10.1016/S0022-2836(02)00811-2View ArticlePubMedGoogle Scholar
- Kahraman A, Morris RJ, Laskowski RA, Thornton JM: Shape Variation in Protein Binding Pockets and their Ligands. Journal of Molecular Biology 2007, 368: 283–301. 10.1016/j.jmb.2007.01.086View ArticlePubMedGoogle Scholar
- Stockwell GR, Thornton JM: Conformational diversity of ligands bound to proteins. Journal of Molecular Biology 2006, 356: 928–944. 10.1016/j.jmb.2005.12.012View ArticlePubMedGoogle Scholar
- Kolodny R, Koehl P, Levitt M: Comprehensive Evaluation of Protein Structure Alignment Methods: Scoring by Geometric Measures. Journal of Molecular Biology 2005, 346: 1173–1188. 10.1016/j.jmb.2004.12.032View ArticlePubMedPubMed CentralGoogle Scholar
- Connolly ML: Analytical molecular surface calculation. Journal of Applied Crystallography 1983, 16: 548–558. 10.1107/S0021889883010985View ArticleGoogle Scholar
- Subbiah S, Laurents D, Levitt M: Structural similarity of DNA-binding domains of bacteriophage repressors and the globin core. Current Biology 1993, 3(3):141–148. 10.1016/0960-9822(93)90255-MView ArticlePubMedGoogle Scholar
- Gerstein M, Levitt M: Comprehensive assessment of automatic structural alignment against a manual standard, the scop classification of proteins. Protein Science 1998, 7(2):445–456. 10.1002/pro.5560070226View ArticlePubMedPubMed CentralGoogle Scholar
- Comin M, Guerra C, Dellaert F: Binding Balls: Fast detection of Binding Sites using a property of Spherical Fourier Transform. Journal of Computational Biology 2009, 16(11):1577–1591. 10.1089/cmb.2009.0045View ArticlePubMedGoogle Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. Journal of Molecular Biology 1995, 247: 536–540.PubMedGoogle Scholar
This article is published under license to 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.