 Methodology article
 Open Access
 Published:
A global optimization algorithm for protein surface alignment
BMC Bioinformatics volume 11, Article number: 488 (2010)
Abstract
Background
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 physicochemical 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 proteinligand binding sites and of proteinprotein interfaces but the problem cannot be considered solved.
Results
In this paper we propose a new method for local structural alignment of protein surfaces based on continuous global optimization techniques. Given the threedimensional 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 wellknown Iterative Closest Point (ICP) method for threedimensional (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.
Conclusions
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 nonredundant proteins and applying a clustering technique to the results of all comparisons to classify binding sites.
Background
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 physicochemical 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 proteinligand binding sites and of proteinprotein interfaces. They include hashing techniques [1, 2], graph theoretic methods [3–6], descriptors based on moments [7] and moment invariants [8], 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 wellstudied 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 [17], the Hausdorff distance [18], the bottleneck distance [19].
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 [20], alphashapes [21] and pseudovertices [2]. 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 wellknown Iterative Closest Point (ICP) algorithm [17] 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 [22] 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 [32], and compare its performance with that MolLoc [12].
Notations and Assumptions
In this section we introduce some notations and assumptions that will be used throughout the paper. Given two protein structures $\mathcal{P}$ and $\mathcal{Q}$, let us denote by P and Q the two finite sets of points corresponding to the atoms of the active sites of the two structures $\mathcal{P}$ and $\mathcal{Q}$, respectively. We let
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.
An isometric transformation in threedimensional space can be defined by a unit quaternion a_{ r }= (a_{0}, a_{1}, a_{2}, a_{3})^{⊤} ∈ ℜ^{4} (a_{ r }  = 1) and by a translation vector a_{ t } ∈ℜ^{3}. Let ${a}^{\top}=({a}_{r}^{\top}{a}_{t}^{\top})$ be the transformation defining vector and denote by T_{ a } the corresponding transformation, so that
for every x ∈ ℜ^{3}, where R (a_{ r } ) is the rotation matrix defined by the unit quaternion a_{ r } as follows:
where
Let Θ ⊂ ℜ^{7} be the set of all vectors a ∈ ℜ^{7} defining an isometric transformation in ℜ^{3}. Given a transformation vector a ∈ Θ, let T_{ a } (P) = P_{ a } denote the set of points obtained by applying the transformation T_{ a } to every point of P, that is
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} .
Let ψ ∈ Ψ be a given mapping and a be a vector defining an isometric transformation, then the mean square error function between P and Q is the following
The surface alignment problem consists in finding a mapping ψ* ∈ Ψ of points in P to points in Q and an isometric transformation a* such that
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 p_{1}, p_{2} ∈ P, p_{1} ≠ p_{2} implies φ(p_{1}) ≠ φ (p_{2}).
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).
Results
Algorithm
A wellknown algorithm for shape alignment is the Iterative Closest Point Algorithm [17]. This algorithm stems from the idea that, once a mapping $\overline{\psi}$ ∈ Ψ is fixed, it is possible to compute the isometric transformation a ∈ Θ that minimizes the function f($\overline{\psi}$, a) (a closedform expression for a($\overline{\psi}$) has been given in [33] where we refer the interested reader for the relevant details). Let a($\overline{\psi}$) be the minimizer of f($\overline{\psi}$, a), that is
Hence, the problem implicitly considered by the ICP Algorithm is the following twolevel optimization problem
As it is stated in [17], where ICP has been originally proposed, the method converges to a solution which is a local minimum of the twolevel Problem (2). Futher, in [17] it has been shown that the final transformation $\overline{a}$ and mapping $\overline{\psi}$ 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 innerlevel problem becomes the one defining the mapping function ψ (instead of the transformation a) once the transformation vector a ∈ Θ is fixed in the outer level.
Namely, we consider the following twolevel optimization problem
Problem (3) can be reduced to a onelevel optimization problem by considering that for every vector a ∈ Θ, the innerlevel problem of (3) admits a globally optimal solution, which we denote by
and represents the closest point mapping. Hence, Problem (3) can be equivalently stated as
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 [17] which can be relevant for n and m large. Moreover, due to its definition, g(a) is a nonsmooth (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 nonsmoothness of function g.
On the basis of the above observation we propose the use of a controlled random search method for the solution to Problem (5). In the following we briefly recall the global optimization algorithm that we use and which was originally proposed in [25] and successively improved in [24]. It is a population based algorithm in the sense that, throughout the entire optimization process, a population of points is maintained and iteratively updated in such a way that they cluster around the global minima of the objective function. Roughly speaking, the method is composed of two distinct and consecutive phases: a global phase and a local phase. During the global phase an initial population of points (defining rototranslations in threedimensional space) is generated by randomly sampling a sufficiently large set of points over some feasible domain. Then, at every iteration of the local phase, a new point is generated and the population is updated if this new point improves on the worst point of the population. More in details, the algorithm can be described by the following steps.

1.
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
$${g}_{\mathrm{max}}{g}_{\mathrm{min}}<\u03f5,$$
where
then STOP.

3.
search phase N +1 points are randomly chosen in the set S. Then,

(a)
the weighted centroid a_{ c } of the N + 1 points is computed;

(b)
the new trial point ã is computed by doing a weighted reflection of the centroid onto the worst point among the selected N + 1 points.
Namely, let a^{†} be the worst point, then
where α ∈ (0,1) is a reflection parameter and
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).
The parameter α is iteratively updated during the optimization process [24] in such a way that its value tends to zero as the iteration count increases and the difference g_{max}  g_{min} decreases.

4.
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 24.
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.
In order to avoid this undesirable effect, we consider the function f(φ, a) defined in (1) where φ ∈ Φ is an assignment function and let, for every a, φ (a) be a global solution to problem
Problem (6) can be formulated as a 0, 1optimization problem and is, indeed, the combinatorial optimization problem known as the Asymmetric Assignment Problem (AAP).
In particular, let G (P, Q, E), E ⊂ P × Q, be the bipartite directed graph characterized by the two sets of nodes P and Q and by the edges between all pairs of nodes, one of P and the other of Q. Then, for every pair e = (p, q) ∈ E, define
Let s ∈ {0, 1}^{E}be the edge incidence vector and consider the following minimum cost assignment problem
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.
Clearly, it is
where s* is the optimal solution to Problem (7).
Problem (7), and hence (6), can be solved very efficiently by adhoc codes that have time complexity $O(\sqrt{nm}\mathrm{log}(n\mathcal{C}({T}_{a})))$ where $\mathcal{C}\left({T}_{a}\right)={\text{max}}_{p\in P,q\in Q}\left\{q\u2013{T}_{a}\left(p\right)\right\},$ see for instance [34].
We are now able to define our new dissimilarity measure, 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 physicochemical 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 hydrophobichydrophilic properties play an important role in proteinprotein and proteinligand interactions. Thus, we consider a variant of our approach in which we integrate physicochemical 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 physicochemical properties. According to [35], we say that p and q have the same physicochemical 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.
Testing
We applied our method, referred to as Continuous Optimization (CO) method in the following, to the comparison of binding sites of proteins. We integrated physicochemical 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 [32]. 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 alltoall 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 [12], 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 alltoall 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 [36] 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 nonredundant and the binding sites are not evolutionary related.
The results of alltoall comparisons are illustrated by means of the distance matrix of Figure 1. An entry of the matrix corresponds to a protein pair and contains a value related to the number of aligned atoms of the binding sites of the pair. Namely, in the matrix we report
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 colorcoded 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 [36], 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 [37], 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.
An important aspect of an alignment method is its ability to retrieve, for a given query binding site, those proteins of the dataset binding the same ligand. To evaluate CO in this task we resort to ROC curves. The results of the comparisons of the query with all other proteins are ranked from the best to the worst in terms of the number of aligned atoms. A pairwise comparison in the ranked list is considered correct or true positive if the other protein of the pair binds the same ligand as the query. The results are summarized by the receiver operating characteristic (ROC) curves in Figure 2 that display the fraction of true positives or correct answers vs. the fraction of false positives for all positions of the ranked solutions. The best possible prediction results would yield a curve through the point in the upper left corner or coordinate (0, 1) of the ROC space. A completely random guess would give a point along a diagonal line from the left bottom to the top right corners. We repeated this experiment with each protein of the dataset as query. Each curve in Figure 2 shows the average values obtained on all query proteins of a group. As expected, the curve corresponding to PO4 (in green) deviates the most from the diagonal line being the closest to the topleft corner of the ROC square. Thus CO has a good success in predicting a PO4 binding site. By contrast, the worst performance is achieved for NAD binding proteins with the associated curve in dotted black.
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 [38] where six publicly available programs are evaluated on almost 9 million pairs of proteins. However, a similar largescale 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 pseudocenters. In [36] the points are the spherical sample points derived from the atomic coordinates. MolLoc, on the other hand, uses Connolly's [39] 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 [36], 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 [38], 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 [42] 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 [43].
We performed pairwise comparisons of the active site of the Catalytic Subunit of cAMP dependent ProteinKinase (pdb code 1atp, chain E) with each of the remaining proteins of the input data set. Of the set of proteins only three belong to the same SCOP family as 1atp, namely 1phk, 1csn and 1hck. In Table 1 for each comparison we report the number of aligned atoms along with the RMSD obtained by CO (columns 34) and MolLoc (colums 67). The entries of the table are listed and ranked according to the number of corresponding atoms obtained by CO (column 3). We observe (see Table 1) that both methods correctly rank at the top three positions the proteins in the same family as 1atp, that is 1phk, 1hck and 1csn. Furthermore for the same three proteins the RMSD is typically very low (approx. 1.5 Å). Lower scores are obtained for distantly related proteins, as for instance 1g5t. The table also reports the SAS measure for CO (column 5) and MolLoc (column 8) and their average at the bottom of the same two columns. As a lower SAS value indicates a better match, it follows that CO on average achieves a better quality than MolLoc with respect to this measure.
Figure 3 shows an example superimposition of the binding sites of ligand ATP of proteins 1atp and 1hck after the computed rototraslation is applied.
We conclude this section by reporting that the execution times of CO on average on all 18 pairwise 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 nonredundant 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.
References
 1.
Shatsky M, ShulmanPeleg 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.407
 2.
ShulmanPeleg 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.012
 3.
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.20140
 4.
Chen B, Bryant D, Fofanov V, Kristensen D, Cruess A, Kimmel M, Lichtarge O, Kavraki L: Cavityaware motifs reduce false positives in protein function prediction. Computational System Bioinformatics Conf 2005, 311–323.
 5.
Hofbauer C, Lohninger H, Aszodi A: Surfcomp: A novel graphbased approach to molecular surface comparison. Journal of Chemical Information and Computer Sciences 2004, 44(3):837–847.
 6.
Weskamp N, Kuhn D, Hullermeier E, Klebe G: Efficient similarity search in protein structure databases by kclique hashing. Bioinformatics 2004, 20(10):1522–1526. 10.1093/bioinformatics/bth113
 7.
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.20681
 8.
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/btm503
 9.
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.0145
 10.
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.
 11.
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.018
 12.
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):
 13.
Ausiello G, Gherardini PF, Marcatili P, Tramontano A, Via A, HelmerCitterich M: Funclust: a web server for the identification of structural motifs in a set of nonhomologous protein structures. BMC Bioinformatics 2008, 9(Suppl 2):S2. 10.1186/147121059S2S2
 14.
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/bti645
 15.
Kinoshita N, Furui J, Nakamura H: Identification of protein functions from a molecular surface database, efsite. Journal of Structural and Functional Genomics 2001, 2: 9–22. 10.1023/A:1011318527094
 16.
ShulmanPeleg A, Shatsky M, Nussinov R, Wolfson HJ: MultiBind and MAPPIS: webservers for multiple alignment of protein 3Dbinding sites and their interactions. Nucleic Acids Research 2008., (36 Web server):
 17.
Besl PJ, McKay ND: A method for registration of 3D shapes. IEEE Trans. on Pattern Analysis and Mach. Intelligence 1992, 14: 239–255. 10.1109/34.121791
 18.
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.
 19.
Efrat A, Itai A, Katz MJ: Geometry helps in bottleneck matching and related problems. Algorithmica 2001, 31: 1–28. 10.1007/s0045300100168
 20.
Connolly ML: Analytical molecular surface calculation. Journal of Applied Crystallography 1983, 16: 548–558. 10.1107/S0021889883010985
 21.
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.5560070905
 22.
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_text
 23.
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:1008250020656
 24.
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.
 25.
Price WL: A controlled random search procedure for global optimization. In Towards Global Optimization 2. Edited by: Dixon L, Szego G. Amsterdam: NorthHolland; 1978.
 26.
Daidone A, Parasiliti F, Villani M, Lucidi S: A new method for the design optimization of threephase induction motors. IEEE Trans. on Magnetics 1998, 34: 2932–2935. 10.1109/20.717684
 27.
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/s1010700405285
 28.
Liuzzi G, Lucidi S, Parasiliti F, Villani M: Multiobjective Optimization Techniques for the Design of Induction Motors. IEEE Transactions on Magnetics 2003, 39: 1261–1264. 10.1109/TMAG.2003.810193
 29.
Price WL: Global optimization by controlled random search. Journal of Optimization Theory and Applications 1983, 40: 333–348. 10.1007/BF00933504
 30.
Price WL: Global optimization algorithms for a CAD workstation. Journal of Optimization Theory and Applications 1983, 55: 133–146. 10.1007/BF00939049
 31.
Price WL, Woodhams F: Optimising accelerator for CAD workstations. IEEE proceedings 1988, 135: 214–221.
 32.
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.235
 33.
Horn BKP: Closedform solution of absolute orientation using unit quaternions. Journal of the Optical Society of America 1987, 4: 629–642.
 34.
Goldberg AV, Kennedy R: An Efficient Cost Scaling Algorithm for the Assignment Problem. Mathematical Programming 1995, 71: 153–178.
 35.
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/S00222836(02)008112
 36.
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.086
 37.
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.012
 38.
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.032
 39.
Connolly ML: Analytical molecular surface calculation. Journal of Applied Crystallography 1983, 16: 548–558. 10.1107/S0021889883010985
 40.
Subbiah S, Laurents D, Levitt M: Structural similarity of DNAbinding domains of bacteriophage repressors and the globin core. Current Biology 1993, 3(3):141–148. 10.1016/09609822(93)90255M
 41.
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.5560070226
 42.
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.0045
 43.
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.
Acknowledgements
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"
Author information
Affiliations
Corresponding author
Additional information
Authors' contributions
All three authors participated in the design of the methodology and in set up of the experiments. GL carried out the implementation. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Bertolazzi, P., Guerra, C. & Liuzzi, G. A global optimization algorithm for protein surface alignment. BMC Bioinformatics 11, 488 (2010). https://doi.org/10.1186/1471210511488
Received:
Accepted:
Published:
Keywords
 Iterative Close Point
 Dissimilarity Measure
 Continuous Optimization
 Global Optimization Algorithm
 Iterative Close Point