Multipolar representation of protein structure
 Apostol Gramada^{1}Email author and
 Philip E Bourne^{1, 2}
DOI: 10.1186/147121057242
© Gramada and Bourne; licensee BioMed Central Ltd. 2006
Received: 22 November 2005
Accepted: 04 May 2006
Published: 04 May 2006
Abstract
Background
That the structure determines the function of proteins is a central paradigm in biology. However, protein functions are more directly related to cooperative effects at the residue and multiresidue scales. As such, current representations based on atomic coordinates can be considered inadequate. Bridging the gap between atomiclevel structure and overall proteinlevel functionality requires parameterizations of the protein structure (and other physicochemical properties) in a quasicontinuous range, from a simple collection of unrelated amino acids coordinates to the highly synergistic organization of the whole protein entity, from a microscopic view in which each atom is completely resolved to a "macroscopic" description such as the one encoded in the threedimensional protein shape.
Results
Here we propose such a parameterization and study its relationship to the standard Euclidian description based on amino acid representative coordinates. The representation uses multipoles associated with residue Cα coordinates as shape descriptors. We demonstrate that the multipoles can be used for the quantitative description of the protein shape and for the comparison of protein structures at various levels of detail. Specifically, we construct a (dis)similarity measure in multipolar configuration space, and show how such a function can be used for the comparison of a pair of proteins. We then test the parameterization on a benchmark set of the protein kinaselike superfamily. We prove that, when the biologically relevant portions of the proteins are retained, it can robustly discriminate between the various families in the set in a way not possible through sequence or conventional structural representations alone. We then compare our representation with the Cartesian coordinate description and show that, as expected, the correlation with that representation increases as the level of detail, measured by the highest rank of multipoles used in the representation, approaches the dimensionality of the fold space.
Conclusion
The results described here demonstrate how a granular description of the protein structure can be achieved using multipolar coefficients. The description has the additional advantage of being immediately generalizable for any residuespecific property therefore providing a unitary framework for the study and comparison of the spatial profile of various protein properties.
Background
The functions of a protein are determined by its three dimensional structure. It is believed that the functional space of all proteins can be spanned by combining a rather small number of structural units, termed folds. The number of different folds is small compared to the total number of proteins, of the order of 1000 for globular, watersoluble proteins [1]. This is in contrast with the exponential complexity of the amino acidlevel configuration space of proteins. Therefore, for some purposes, the description of the configuration based on amino acid coordinates is overdetailed. The number of degrees of freedom needed for the description, comparison and classification of the macroscopic, biologicallyrelevant features of proteins is necessarily much smaller than that associated with the collection of amino acids. Selecting the relevant degrees of freedom and defining methods to compare structures in their reduced space is deemed useful.
The comparison of structures is a central component of many research objectives. For formulations of the problem and a review of various methods of comparison see for example [2]. The problem is particularly important in protein classification. Various classification schemes that have been developed thus far (for example SCOP [3], CATH [4, 5], FSSP [6]) approach the objective of selecting the degrees of freedom by using a combination of various empirical and/or qualitative descriptors of protein conformations and parameters derived directly from coordinates of corresponding amino acids (interdistances for example), recognizing implicitly that, for the purpose of describing its function, the conformation of the proteins is highly unstructured in the microscopic (amino acid) configuration space. Despite this, however, when it comes to quantitatively measure the difference between protein structures, in most cases, the measure of choice remains the root mean square deviation (rmsd) between aligned atomic coordinates.
A disadvantage of the rmsdbased measures is that it assumes a strict onetoone correspondence between at least the C alpha (Cα) atom coordinates of the compared proteins, i.e., it requires an alignment at the amino acid level of detail. Unfortunately, there is no unique formulation for the problem of aligning protein sequences from structure and, as a result, the existing methods produce results which differ in details [7]. Moreover, there exist problems where the shape of a certain region in the molecule is what is analyzed in the context of a biological function. For example, binding sites from different proteins need to be compared for the detection of any similarities [8]. In those cases there is no alignment, or the alignment between amino acids is not relevant and the rmsd cannot be defined. Thus far, from a shape perspective, the comparison has included only the molecular surface in that region [8, 9]. One limitation (induced by its twodimensional nature) is that some methods of surface representation impose restrictions on what kind of surfaces can be studied (for instance only starlike surface, i.e. surfaces described by functions on a unit sphere [8]). Another disadvantage is that it may be the case that other properties of the site, which are threedimensional in nature might be relevant in defining the function. For example, the distribution of charge, hydrophobicity, etc deep beneath the surface of the site. These properties can not be included in a natural way when only the surfaces are compared. In other words, the intrinsic behavior of a protein is a combination of its properties defined through both the internal and external structure of the protein. To more correctly represent a protein warrants the search for alternative parameterizations of protein structure. The approach presented here represents such an alternative.
A fold is characterized by structural features at the multiresidue level. Even though these features are easily recognizable visually in many cases, there is no obvious quantitative way to relate them to the underlying atomic coordinates that exactly describe the structure of the protein. That is, atom coordinates only offer a local description of the structure, while the features defining the fold represent global, shape related properties. Starting with the very rich description given by the coordinates of the atoms that make the protein, one would then need a way to distill from this set of coordinates only the information that is directly associated with these general traits. We do not know of the existence of any systematic approach for the elimination of the redundant information, starting from the initial set of coordinates. Here, we adopt an approach that starts from the other end of the spectrum. Instead of starting from the atom description and discarding nonrelevant features, we start from the global level with a very coarse description and refine it by adding descriptors for more and more detailed features.
The need for methods that use global descriptors for the comparison of protein structures has been recognized before. One such method [10, 11] has been proposed recently that relies on results from knot theory to extract a number of quantitative features of the path of the protein backbone and compare two such paths in the reduced space of this set of features [12]. Here we present another method that is based on a hierarchical set of descriptors for the distribution of atoms in space. The method is general enough so that it can be refined to describe the protein spatial profile with a geometrical level of detail. It can potentially interpolate between the coarsest description of the structure (knowledge of the number of atoms only) and the most detailed description, equivalent in amount of information to the complete set of atom coordinates. This is achieved by using as coordinates the series of multipolar components associated with a given atom property. To be specific, we will refer here to the multipole tensors associated with the mass of Cα atoms. Note, however, that as long as the property is distributed uniformly over the set of atoms (in our case the same mass is assigned to each atom), the nature of the property is not relevant. The multipoles themselves represent, up to a constant multiplicative factor, a set of parameters describing exclusively the spatial configuration of the atom set.
Our approach can be combined with other important parameters, for example, mass of the residue instead of the Cα, charge, hydrophobicity or secondary structure conformation. In those cases, the method will enable the description of the spatial profile of those quantities instead of just the pure geometrical configuration of the molecule. Correspondingly, the method will provide a means to quantitatively compare the proteins according to those properties.
The notion of multipoles originates in physics [13] and is closely related to the representation of functions in terms of spherical harmonics [14]. The use of spherical harmonics in biomolecular research is not new. They have been used for example for the representation and rotation of molecular surface and other properties in an efficient way [15, 16], for the purpose of molecular docking [17], for the comparison of binding sites in molecules [8] or for the display of molecular surfaces in molecular visualization [18]. The lower order (or rank) multipoles (up to the quadrupole order) have been used before as a signature for the electrostatic field in the comparison of small molecules for the purpose of drug design [19]. Here, we take the approach that the whole set of multipoles can be interpreted as an alternative set of coordinates for the description of the structure of the molecule. We then show how their tensorial properties can be used for the definition of a distance function in protein configuration space.
The organization of the paper is as follows. In the Results and Discussion section we first define the multipoles and present a qualitative motivation for their use as an alternative parameterization of the shape of the protein. Then we show how their tensorial properties can be used to define a distance function in the conformational space. Since the multipoles are dependent on the location and orientation of the system of axes, the following subsection is used to define a "canonical" reference frame to be used for the purpose of comparison. The concluding subsection is devoted to testing the method. We show that, given the biological relevant portion of the structure, the method successfully discriminates between the families in a test set of proteins from the protein kinaselike superfamily [20]. We then study the correlation between the multipole and Cartesian coordinates representations for the same test set and show that, as expected, the correlation increases with the level of detail of the multipole representation approaching the dimensionality of the fold space. We conclude the paper with a discussion of the advantages of the method and of the various directions in which it can be generalized.
Results and discussion
Mutlipolar representation
The notion of multipoles comes from physics where they are used to describe the field generated by the spatial distribution of a scalar quantity such as mass (gravitational field) or charge (electrostatic field) density. The potential of the field created by such a quantity, at a given distance outside the region it occupies in space, can be conveniently expressed in the form of a multipole expansion [13]. Each multipole in the series accounts for the contribution of a certain type of deviation of the density field from a spherically symmetric distribution and, in general, the higher the order of the multipole the smaller the spatial scale of the deviation it describes. In this sense, the multipoles can be viewed as descriptors for the shape of the scalar distribution. The use of multipoles as shape descriptors is also closely related to the more general methods of 3D moments used in the field of object recognition in computer science [21–23]. In the above physical example the multipoles of higher orders typically account for relatively small contributions in the force field compared to lower order multipoles, and they may be neglected for many practical purposes. In a similar way, higher order multipoles representing small scale details of the shape, can be ignored in the process of describing protein structures when only a rough comparison is needed. Before giving a formal definition of multipoles in general, let us start by discussing the few lower order multipoles which are more familiar and more widely encountered in the research literature.
Quantitatively, the multipoles associated with the space distribution of a scalar property (density of mass in our example) form a sequence of tensors over the three dimensional position space. The multipole of rank zero, or the monopole, is just the space integral of the scalar property. When the scalar property is the mass density then the monopole is the total mass of the set of atoms.
The multipole of rank one is proportional to the position of the center of mass. We will use it to set the origin of the coordinate system with respect to which all multipoles are calculated. Therefore, in our calculations, the multipole of rank one, (the dipole as it is commonly known), is always going to be a null vector. For completeness, we should mention that for the multipoles of a distribution of charge this can not always be done. If the total charge is zero there may be a nonzero dipole moment that can not be made to vanish by a translation of coordinates. This is however a technical problem and it has been addressed before [19].
The multipole of rank two, or the quadrupole has nine Cartesian components. For our discrete distribution, the components are given by the following expression:
Here, δ_{ i,j }= 1 for i = j and 0 otherwise. The sum runs over all N Cα atoms in the structure, x_{α,i}is the component i of the ${\overrightarrow{x}}_{\alpha}$ position vector (one of the x, y, z Cartesian components) and r_{ α }represents the length of the ${\overrightarrow{x}}_{\alpha}$ position vector. For example, the first diagonal component is $\sum \left(3{x}_{\alpha}^{2}{r}_{\alpha}^{2}\right)={\displaystyle \sum \left(2{x}_{\alpha}^{2}{y}_{\alpha}^{2}{z}_{\alpha}^{2}\right)}$ and one of the nondiagonal terms is ∑x_{ α }y_{ α }.
For higher order multipoles, enumerating the Cartesian components in closed form is not a simple task. Moreover, there is a large number of symmetry properties obeyed by the Cartesian components and therefore keeping all components of a given multipole would be redundant. Instead, more commonly, the irreducible spherical components are defined since they have a compact form when represented in terms of spherical harmonic functions and are independent. This makes them suitable for analytical and numerical calculations. Within the rest of the paper we will use the term multipoles to denote these irreducible components which represent the focus of our attention. We will explicitely name the Cartesian components if needed to distinguish them from the spherical components.
For a discrete set of N atoms the multipoles of rank l are defined as:
where r_{ i }, θ_{ i }, φ_{ i }represent the spherical coordinates of atom i, Y_{ lm }denotes a spherical harmonic function and the * denotes complex conjugation. For the definition and summary properties of these functions, see for example [24]. Here, since we will only consider the Cα atoms, we set the mass of each atom to unity to simplify the notation. For a set of arbitrary atoms, each term in the sum would be weighted by the mass of the atom (or another scalar property in a generalized case).
The rank l can take any integer value from 0 to ∞ and for each given l the number m can take values in the range l...l. Then, the number of irreducible components, specified by the index m, increases linearly with the rank l of the multipole as 2l + 1. When all multipoles with rank from 0 to n are used, the total number of independent components describing the shape of the protein is (n + 1)^{2}. As n increases, this number approaches the number of Cartesian components of the position vectors of the Cα atoms. When the two numbers are equal, the description provided by the set of multipolar components for the structure of the protein is of the same level of detail as the original description offered by the atomic coordinates. We fully recover the amount of information provided by those coordinates. When this happens, from a mathematical standpoint the multipole series is just a coordinate transformation and if not singular, at least in principle, we can transform back and forth from one description to the other.
As a last remark, we will note that the set of multipoles that we use here is only a subset of a larger set which, in its entirety, uniquely describes the potential field surrounding the distribution of charge (for example) [13], for a given set of boundary conditions. The formalism that we present can be extended to include any portion of this complete set of multipoles and this leaves open the possibility for further optimization of a protein comparison process. The reason for retaining this particular subset of multipoles is that they describe the field outside the region occupied by the molecule and therefore they are more likely to correlate with its interaction capabilities and consequently its function.
Constructing a distance function in the protein conformation space
What makes the set of multipoles defined in Eq. (2) a good set of descriptors for comparison purposes is that they form a series of quantities with remarkable symmetry properties. Specifically, for any given rank l, the 2l + 1 components q_{ lm }, m = l, l 1, ... l + 1, l form an irreducible tensorial set of order l [25]. This means that under regular rotations in the threedimensional (3D) physical space, these components are transformed according to a well defined induced rotation matrix (see e.g. [26]), in a way similar to the behavior of a 2l + 1 dimensional vectorial quantity. The immediate benefit is that one can apply the regular operations with vectors to the multipoles of a given rank and one can construct invariant quantities following the known rules from Euclidian vectors. In particular, if we denote by q_{ 1 }the set of all components ${\left\{{q}_{lm}\right\}}_{m=\overline{l,l}}$, we can define the length of the multipole of rank l using the scalar product:
where the last part of the equation follows from the definition of the multipoles and well known symmetry properties of the spherical harmonic functions [24]. This norm can then be used to define a distance between two structures inside the subspace defined by the multipolar components (say q_{ 1 }, ${{q}^{\prime}}_{1}$) of a given rank, provided that either the structures have been previously spatially superimposed, or, a "canonical orientation" of the structures has been set in some consistent manner:
δ_{ l }(q_{ 1 }, q_{ 1 }') = q_{ 1 } q_{ 1 }'. (4)
There is no a priori prescription for combining distances in subspaces of different ranks l to construct a global distance function. Such a prescription needs to be extracted from numerical experimentation with the problem to be modeled and/or from more general principles. The alternative discussed here is the result of our tests of the sensitivity and selectivity in discerning protein structures. Since the dimensionality of multipoles differs with l, in order to combine distances from different subspaces to construct a global metric one has to first define quantities with the same dimensionality. The solution adopted in this paper, is to redefine the distance in all ranks so that it has the same dimensionality, say dimensionality of length. Except eventually for a general factor (with the dimension of mass in our example), the dimensionality of the multipoles is a power of length equal to their rank. The general factor can be rendered dimensionless by a proper rescaling. Then, one can obtain a quantity with dimension of length by taking an appropriate root of the Euclidian distance as follows:
d_{ l }(q_{ 1 }, q_{ 1 }') = q_{ 1 } q_{ 1 }'^{1/l}, l > 0. (5)
Note that, once the multipole components have been calculated, any reference to the original Cartesian coordinates disappears from the representation. As a consequence, unlike the rmsd which requires a onetoone correspondence between the set of atoms in the structures compared, the distance in Eq. (5) is defined for arbitrary structures, without any restriction with respect to their number or sequence of aminoacids. Therefore no alignment is implied. In practice, a normalization with respect to the "size" of the structures involved may still be necessary. For that purpose, each multipole in Eq. (5) can be separately rescaled with a factor inversely dependent on the "size" of the corresponding molecule. Then, instead of Eq. (5), we will use
The notation we use in this formula for the "size" dependent factors (q_{0}, ${{q}^{\prime}}_{0}$) is motivated by the fact that the multipole of rank 0 (monopole) is up to a constant numerical factor the "size" of the molecule (the number of atoms for example). When the two structures have the same size, the rescaling of the multipoles in Eq. (6) reduces to a rescaling of the distance (5) by a factor inversely proportional to the common size of the two molecules. This is qualitatively equivalent to the 1/N factor in the rmsd distance (Eq. (10)). It can be shown that Eq. (6) satisfies the triangle inequality and therefore a global distance function which also satisfies the triangle inequality can be defined by adding the distances (multiplied eventually by a weight factor) for all ranks of the multipoles. Here, we use the following formula:
This function will be used as a dissimilarity measure for proteins in our study. The upper limit in the summation is the maximum rank of the multipoles retained in the representation and determines the dimensionality of the representation and, implicitely, its level of detail.
The interest for reduced representations is manifest in the literature. From a shape perspective, similar to ours, such representations emerge in approaches such as that described in [10, 11]. From a different perspective, starting from individual atomic coordinates and using an averaging approach, an alternative method is presented in [27].
Defining a canonical reference frame
The multipoles behave like vectorial quantities and the numerical values of their components depend on the location and orientation of the reference frame. For the comparison of structures to be meaningful, we need to either minimize the distance in Eq. (7) over all rigid transformations (translations and rotations) of one of the molecules, or to choose a standard for the reference frame with respect to which the multipolar configurations of the molecules are calculated [8]. Since the second approach is much more efficient for large scale computations, we chose to test this second alternative.
The problem of choosing such a standard arises in many research areas where 3D systems are involved [21], and various schemes can be found in the literature, depending of the research field. A common choice is to select a system of axes that is placed at the center of mass and having its three orthogonal directions along the eigenvectors of a suitable, symmetric matrix (principal axes reference frame). While the location at the center of mass is natural, its orientation as described above is ambiguous. This prescription is not appropriate in our case since it does not uniquely define the axes: any combination of permutations and inversions of the versors of a given principal axes frame form also a principal axes frame. Since we are using multipoles of ranks higher than the quadrupole (rank two), which are sensitive to these various orientations, we can not allow an arbitrary choice. We need a prescription that uniquely defines a frame. Our choice is based on the use of lower order vectorial (magnetic) counterparts of the multipoles introduced above. Specifically, we start by defining the following vectors:
The first vector reduces to the relative position of the last amino acid with respect to the first, while the second one is a more complex quantity that is sensitive to the details of the path of the protein backbone. Except for special cases (for example when the two vectors in Eqs. (8, 9) are not welldefined, or they become parallel), these vectors are independent. Then, they can be orthonormalized and the resulting unit vectors will serve as the first two versors of our canonical reference frame. The third one will be their cross product.
The "canonical" reference frame defined by Eqs. (8, 9) is unique by construction. However, other unique definitions can be developed [8]. We are not aware of any rigorous prescription for constructing such a reference frame and therefore our choice remains heuristic.
Testing the multipole representation
To test our representation of protein structure, we performed a number of calculations with the goal of assessing both its discriminatory power and, where meaningful, its correlation with the Cartesian description.
Comparing biologically relevant molecules
As already stated, the use of multipoles opens the possibility of protein shape comparisons without the need for a preexisting amino acid alignment. However, while technically our representation allows for the comparison of arbitrary collections of atoms, in biological applications, such as protein classification, not any comparison will make sense: we need to restrict the comparison to those portions of the proteins which are relevant to the problem, for example, the functional regions. In principle, the multipoles can be used in identifying corresponding domains in structures, however, as of this moment we do not have fully functional tools to do that. Therefore, as a benchmark for testing the method, we use a manual alignment of the catalytic cores from the protein kinaselike superfamily [20]. The set contains 25 typical protein kinases (TPK) and 6 atypical protein kinases (AK) which phosphorylate nonproteins. As has been shown [20] these diverse structures can be traced to a common ancestor, but today their sequence identity is below 15% in some cases and significant structural changes have taken place, particularly in the Cterminal lobe so that a variety of substrates can be phosphorylated using the same ATP gammaphosphate transference mechanism. These structures represent an excellent test case for structure recognition since the accurate hand curated alignment provides a valuable benchmark.
Even at the subfamily level with relatively little shape discrimination, the distance matrix retains some of its discriminatory power. A close examination reveals distinct patterns along the diagonal corresponding to the various groups of kinases in the test set (Figure 1).
The discriminatory power at the family level is limited by unresolved portions of some structures. The lack of coordinates for parts of the polypeptide chain affects the calculated distances both directly (a missing piece of chain is seen as a difference in shape) and indirectly (a missing piece of chain leads to a different canonical reference frame). To reduce these perturbations, we chose to ignore in our calculations any portion of the alignment corresponding to missing parts in at least one of the proteins in the set. Most unresolved portions are relatively short (approximately 20 amino acids) and do not affect the shape dramatically.
Correlation with the Cartesian representation
The multipolar description offers a hierarchical approach to characterizing the shape of a molecule. While at the coarsest level there is no information about the shape, except that defined by the length of the chain, at the most refined level of details (when the number of multipole components is of the same order of magnitude as the original Cartesian coordinates) the description is as rich as the original amino acid coordinate set. At this end of the spatial spectrum we would expect a good correlation with results provided by the Cartesian coordinates. To empirically prove this, we need to devise experiments in which both representations can be applied and then compare the results.
An obvious choice is the comparison of aligned proteins, alignment being necessary for the rmsd to be defined. Since as yet we do not have tools for aligning proteins based on their multipole representation, we used the high quality expert alignments provided by the authors of the benchmark set [20]. We performed an analysis of how well distances calculated based on the multipole representation compare with the ones based on the Cartesian coordinates. Two different cases were considered, each defined by how the Cartesian and multipolar representation were calculated in the two proteins.
Case 1
In the first case, the rmsd distances were calculated based on a prior spatial superposition of the aligned structures (the typical approach for assessing structural similarity). The rmsddistance was calculated using the formula
The two vectors in Eq. (10) denote the coordinates of aligned residues. The multipoles of the aligned portions of the proteins were calculated with the coordinates expressed in the canonical reference frame defined by Eqs. (8, 9).
The multipoles of each protein in a given pair were computed from the Cα coordinates only. The distances were calculated several times, each time retaining a different range of multipoles to analyze the results at different levels of structural detail. The coarsest calculation corresponds to retaining only multipoles up to rank 2 (quadrupole) and the finest one contains all multipoles up to rank 12. For each set of distances obtained in this way we calculated their linear correlation coefficients with the set of rmsd values.
Case 2
It is clear that, while the multipole description differs in some intrafamily details from the typical rmsd results, especially when the latest are calculated with spatial superposition, the results are quite robust in their capacity to discriminate between the families. The level of correlation increases with the level of detail of the representation, i.e. with the number of multipole approaching the dimensionality of the fold space. These latest results, when compared with the previous ones, suggest that the two ways of positioning the structures for the purpose of comparison are not entirely equivalent and the use of a "canonical" positioning produces a more robust clustering of the structures.
Conclusion
In this paper we propose a new parameterization of protein structure which provides a new form of characterization and comparison. The approach uses components of the multipoles of consecutive ranks associated with Cα coordinates.
We have shown:

Once an approximate "superposition" has been calculated using our canonical reference frame, the multipole distance function is capable of discriminating between protein families.

The multipole description allows for the adjustment of the level of detail of the comparison and, implicitly, it provides a systematic method for deriving reduced representations of the protein configuration space.
From a biological perspective, our tests show that the comparison based on multipoles is more robust with respect to intrafamily details and the results are more meaningful biologically. From the comparison tests with the Cartesian description, its robustness appears to be related in part to the use of a "canonical" reference frame for the comparison rather than the spatial superposition of the structures. Also, the visible relationship between the distance matrix in Figure 1 and the family classification of the set, in contrast to the distance based on an exact amino acid alignment, suggest that multiresidue shape features are more important to the biological classification than local variations in the alignment. This supports the idea that evolution is more likely driven by shape optimization as required by molecular recognition. Evolutionary events such as sequence insertions and deletions are merely the means to achieve optimal shape complementarity.
For illustration of the multipole method, we used the mass of the Cα atoms as the relevant physicochemical property. This led to comparison of proteins with respect to their geometrical structure. The method can be refined with minimal changes to use various residuespecific quantities. The only technical difference will be the use of a weight for each term in Eq. (2). The weight is a numerical functions that measure the property of interest, such as hydrophobicity, assigned charge, numerically encoded secondary structure information, etc. It can also represent a composite index reflecting a set of properties assumed to be relevant for investigating a given biological concept.
The use of alternative residuespecific quantities would provide a powerful tool for the comparison of proteins since the residue specific quantities allow an easier discrimination between structures with similar spatial location of the Cα atoms but differing in local properties of the chain, such as secondary structure conformation for example.
 a)
Rigorous definition of the notion of "canonical" reference frame. Our choice, based on features rigidly tied to the set of atoms is inspired by the body reference frames used in physical and engineeringsciences and is intuitive. However, the problem of comparing structures is different and criteria are needed for the identification of "good" reference frames and/or how they affect the protein comparison.
 b)
Algorithms for fast superposition by minimization of the multipolar distance would be needed as an alternative to the use of a "canonical" reference frame.
 c)
The definition of a global metric (Eq. 7) contains coefficients controlling the combination of multipoles of various orders. Further optimization of these coefficients for the purpose of protein comparison can lead to biologically more meaningfull metrics.
As a final remark, our representation allows an estimation of the number of degrees of freedom necessary to describe a given class of properties. The saturation of the correlation with the Cartesian representation marks the maximum number of degrees of freedom necessary to "macroscopically" distinguish structures within that class. Since the structure determines the whole biology of the proteins, one can infer from here that the same number of degrees of freedom describes the whole functional space of that class of proteins. The number obtained from such a correlation curve can be used to adjust the dimensionality of the representations used in protein comparisons.
Methods
The atomic coordinates of the selected members of the protein kinaselike superfamily were obtained from the ASTRAL database [31]. For each pair of proteins, we retain only the Cα atoms of the biologically relevant parts of the proteins as determined in [20]. The calculations presented in the paper were initially prototyped in Mathematica [32]. A Java program was subsequently written to test the performance of the calculations. We tested the program on a notebook computer with a 1.2 GHz Pentium III processor. The comparison of proteins with about 300 residues, with a level of detail corresponding to l_{ max }= 8 takes of the order of 100 ms.
Declarations
Acknowledgements
We are grateful to Marian Anghel and Eric Scheeff for very useful discussions and to Yuting Jia for providing the programs for the spatial superposition of aligned structures. We are grateful for financial support from grant NIGMS GM63208.
Authors’ Affiliations
References
 Chothia C: Proteins. One thousand families for the molecular biologist. Nature 1992, 357: 543–544. 10.1038/357543a0View ArticlePubMedGoogle Scholar
 Eidhammer I, Jonassen I, Taylor WR: Structure Comparison and Structure Patterns. J Comput Biol 2000, 7(5):685–716. 10.1089/106652701446152View ArticlePubMedGoogle Scholar
 Murzin A, Brenner S, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequence and structures. J Mol Biol 1995, 247: 536–540. 10.1006/jmbi.1995.0159PubMedGoogle Scholar
 Orengo C, Michie A, Jones S, Jones D, Swindells M, Thornton J: CATH – A hierarchical classification of protein domain structures. Structure 1997, 5: 1093–1108. 10.1016/S09692126(97)002608View ArticlePubMedGoogle Scholar
 Pearl F, Martin N, Bray J, Buchan D, Harrison A, Lee D, Reeves G, Shepherd A, Sillitoe I, Todd A, Thornton J, Orengo C: A rapid classification protocol for the CATH domain database to support structural genomics. Nucleic Acids Res 2001, 29: 223–227. 10.1093/nar/29.1.223PubMed CentralView ArticlePubMedGoogle Scholar
 Holm L, Sander C: Mapping the protein universe. Science 1996, 273: 595–602.View ArticlePubMedGoogle Scholar
 Godzik A: The structural alignment between two proteins: is there a unique answer? Protein Sci 1996, 5(7):1325–1338.PubMed CentralView ArticlePubMedGoogle Scholar
 Morris RJ, Najmanovich RJ, Abdullah K, Thornton JM: Real spherical harmonic expansion coefficients as 3D shape descriptors for protein binding pocket and ligand comparison. Bioinformatics 2005, 21(10):2347–2355. 10.1093/bioinformatics/bti337View ArticlePubMedGoogle Scholar
 Rosen M, Shuo Liang L, Haim W: Molecular shape comparison in Search for active sites and functional similarity. Protein Eng 1998, 11(4):263–277. 10.1093/protein/11.4.263View ArticlePubMedGoogle Scholar
 Røgen P, Fain B: Automatic Classification of Protein Structures by Gauss Integrals. Proc Natl Acad Sci USA 2003, 100: 119–124. 10.1073/pnas.2636460100PubMed CentralView ArticlePubMedGoogle Scholar
 Røgen P, Bohr H: A New Family of Protein Shape Descriptors. Math Biosci 2003, 182: 167–181. 10.1016/S00255564(02)00216XView ArticlePubMedGoogle Scholar
 BarNatan D: On the Vassiliev Knot Invariants. Topology 1995, 34: 423–472. 10.1016/00409383(95)932372View ArticleGoogle Scholar
 Jackson J: Classical Electrodynamics. third edition. New York: John Wiley & Sons, Inc; 1999.Google Scholar
 Tannoudji CC, Diu B, Laloë F: Quantum Mechanics. New York: John Wiley & Sons, Inc; 1977.Google Scholar
 Ritchie DW, Kemp GJL: Protein docking using spherical polar Fourier correlations. Proteins 2000, 39: 179–194. Publisher Full Text 10.1002/(SICI)10970134(20000501)39:2<178::AIDPROT8>3.0.CO;26View ArticleGoogle Scholar
 Crowther RA: The Fast Rotation Function. In The Molecular Replacement Method: A Collection of Papers on the Use of Noncrystallographic Symmetry. Edited by: Rossmann MG. New York: Gordon and Breach; 1972.Google Scholar
 Ritchie DW, Kemp GJL: Fast computation, rotation, and comparison of low resolution spherical harmonic molecular surfaces. J Comput Chem 1999, 20: 383–395. Publisher Full Text 10.1002/(SICI)1096987X(199903)20:4%3C;383::AIDJCC1%3E;3.0.CO;2MView ArticleGoogle Scholar
 Duncan BS, Olson AJ: Shape analysis of molecular surfaces. Biopolymers 1993, 33: 219–229. 10.1002/bip.360330204View ArticlePubMedGoogle Scholar
 Platt DE, Silverman B: Registration, Orientation, and Similarity of Molecular Electrostatic Potentials through Multipole Matching. J Comput Chem 1996, 17: 358–366. Publisher Full Text 10.1002/(SICI)1096987X(199602)17:3<358::AIDJCC10>3.0.CO;2GView ArticleGoogle Scholar
 Scheeff E, Bourne P: Structural Evolution of the Protein KinaseLike Superfamily. PLoS Comput Biol 2005, 1(5):e49. 10.1371/journal.pcbi.0010049PubMed CentralView ArticlePubMedGoogle Scholar
 Kazhdan M, Funkhouser T, Rusinkiewicz S: Symmetry Descriptors and 3D Shape Matching. Symposium on Geometry Processing 2004.Google Scholar
 Lo C, Don H: 3D moment forms: Their construction and application to object identification and positioning. IEEE Trans Pattern Anal Mach Intell 1989, 11: 1053–1064. 10.1109/34.42836View ArticleGoogle Scholar
 Burel G, Henocq H: Threedimensional invariants and their application to object recognition. Signal Processing 1995, 45: 1–22. 10.1016/01651684(95)00039GView ArticleGoogle Scholar
 Abramovitz M, Stegun I: Handbook of Mathematical Functions. New York: Dover; 1970.Google Scholar
 Fano U: Irreducible Tensorial Sets. New York: Academic Press; 1958.Google Scholar
 Biedenharn L, Louck J: Angular Momentum in Quantum Mechanics, Theory and Applications. AddisonWesley Publ. Co.; 1981.Google Scholar
 Lotan I, Schwarzer F: Approximation of Protein Structure for Fast Similarity Measures. J Comput Biol 2004, 11(2–3):299–317. 10.1089/1066527041410355View ArticlePubMedGoogle Scholar
 Manning G, Plowman G, Hunter T, Sudarsanam S: Evolution of protein kinase signaling from yeast to man. Trends Biochem Sci 2002, 27: 514–520. 10.1016/S09680004(02)021795View ArticlePubMedGoogle Scholar
 Manning G, Whyte D, Martinez R, Hunter T, Sudarsanam S: The protein kinase complement of the human genome. Science 2002, 298: 1912–1928. 10.1126/science.1075762View ArticlePubMedGoogle Scholar
 Ortiz AR, Strauss CE, Olmea O: MAMMOTH (Matching molecular models obtained from theory): An automated method for model comparison. Protein Sci 2002, 11: 2606–2621. 10.1110/ps.0215902PubMed CentralView ArticlePubMedGoogle Scholar
 Chandonia J, Hon G, Walker N, Lo Conte L, Koehl P, Levitt M, Brenner S: The ASTRAL compendium. Nucleic Acids Res 2004, 32: D189D192. 10.1093/nar/gkh034PubMed CentralView ArticlePubMedGoogle Scholar
 Wolfram Research I: Mathematica. Chmpaign, Illinois: Wolfram Research, Inc., version 5.2 edition; 2005.Google Scholar
Copyright
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.