# Protein Surface Characterization Using an Invariant Descriptor

Authors: Zainab Abu Deeb; Donald A. Adjeroh; Bing-Hua Jiang
Journal: International Journal of Biomedical Imaging
 (2011)

---

## Abstract

Aim. To develop a new invariant descriptor for the characterization of protein surfaces, suitable for various analysis tasks, such as protein functional classification, and search and retrieval of protein surfaces over a large database. Methods. We start with a local descriptor of selected circular patches on the protein surface. The descriptor records the distance distribution between the central residue and the residues within the patch, keeping track of the number of particular pairwise residue cooccurrences in the patch. A global descriptor for the entire protein surface is then constructed by combining information from the local descriptors. Our method is novel in its focus on residue-specific distance distributions, and the use of residue-distance co-occurrences as the basis for the proposed protein surface descriptors. Results. Results are presented for protein classification and for retrieval for three protein families. For the three families, we obtained an area under the curve for precision and recall ranging from 0.6494 (without residue co-occurrences) to 0.6683 (with residue co-occurrences). Large-scale screening using two other protein families placed related family members at the top of the rank, with a number of uncharacterized proteins also retrieved. Comparative results with other proposed methods are included.


---

## Body

## 1. Introduction

The Protein Data Bank (http://www.pdb.org/pdb/home/home.do) (PDB) currently has more than 3000 protein structures classified as uncharacterized or as proteins of unknown function. This is about 5% of the total structures in PDB. The Pfam database was recently reported to contain over 2200 gene families with unknown function [1]. It has been argued that there are even more local regions on the protein structures that are not completely characterized, and whose functions are not known [2]. Therefore, with the increasing rate at which protein structures are being generated, the problem of protein function annotation has become a major challenge in the postgenomic era [3–5]. The function of a given protein is largely determined by its three-dimensional structure [6]. The specific shape and orientation of a protein in 3D space are key elements that determine how the protein interacts with its environment, and hence the function of the protein. Although related proteins often have similar functions, it is well known that sequence similarity between proteins does not always lead to functional similarity [7, 8]. Even different functions have been observed for structures with the same fold [9].  Conversely, sequences have been observed with low sequence similarity, but highly structural and functional similarity [10]. The trypsin-like catalytic triad [9] is one example of proteins with different folds, but similar functions. A similar argument can be made between sequence and surface, and between surface and fold. While residues on the protein surface typically make up a small percentage of the total residues in a protein, they often represent the most conserved functional elements of the protein [11]. Therefore, analyzing protein structures using information about their 3D surfaces is essential in the quest for protein function annotation, especially in the study of functional similarities between nonhomologous proteins.At the core of most activities in the analysis of protein structures and protein function is similarity measurement between structures. Such measurements must deal with different levels of structural similarity, arbitrary mutations, deletions, and insertion of residues, local surface similarities, and so forth. When the problem is similarity measurement between protein surfaces, a major issue becomes how the protein surface is represented, and how the representation can be used for the required similarity measurement. Another problem is that of computation. Structure alignment, the basis for most approaches to protein 3D structure analysis is known to be NP-hard [12]. A major difficulty in comparing protein surfaces locally is the problem of matching 3D structures, since structures need to undergo an exhaustive amount of rotation and translation in order to obtain an adequate structural alignment and to perform an accurate matching [8]. Clearly, a method that avoids the step of local structural alignments can have a significant advantage, especially in screening of similar surfaces over a large database.In this paper, we introduce an invariant descriptor for the characterization of protein surfaces. We then use this characterization to study the problem of classifying proteins into their functional families based primarily on their surface characteristics. This is a challenging problem, but one that is important in the quest for functional annotation of proteins, using information from potentially nonhomologous proteins. We also show how we can use such a descriptor in various related analysis activities, such as in effective retrieval of similar protein surfaces from very large databases, such as the Protein Data Bank (PDB).


## 2. Background and Related Work

### 2.1. Protein Sequences, Structure, and Surface

Although proteins could vary significantly in their functions and 3D shapes, they also share a general common structure. Proteins are composed of 20-amino acids that are connected via peptide bonds [13]. Each protein is composed of an ordered sequence of amino acids. The order in which these amino acids are connected is called the protein sequence, or the primary structure of the protein [14]. This primary sequence determines the 3D structure of the protein. All proteins are composed of four common structural types: primary structure, secondary structure, tertiary structure, and quaternary structure. The primary structure is simply the amino acid sequence. The secondary structure is formed by patterns of intermolecular bonding of hydrogen and is determined primarily by the location and the directions of these patterns [14, 15]. This is often described in terms of secondary structural elements (SSEs), such as α-helixes, β-sheets, and turns. The overall 3D shape of the secondary structures determines the tertiary structure of the protein. When two or more chains combine to form a larger molecule, the whole structure is called the quaternary structure. Figure 1 shows an example of some of the common protein structural types (the sequence is not included).Protein structures for a sample protein (PDB id: 2UDI). (a) Secondary structure elements—α-helixes (magenta), β-sheets (gold), and turns (gray); (b) two chains: chain E (blue), chain I (green);  (c) surface and 3D shape for chain E; (d) surface and 3D shape for chain I; (e) quaternary structure for the protein. Figures are produced using PMV [16].
(a)(b)(c)(d)(e)A common method for protein function prediction is by annotation transfer from known homologous proteins [17]. Functions of novel proteins can be determined by sequence comparisons, for instance using sequence alignment. When proteins evolve, the protein structure remains more highly conserved when compared to the sequence. Protein sequences change more easily during evolution due to residue mutations, for instance by substitution, insertion, or deletion. Hence, proteins that belong to the same family (homologous proteins) may not be identified using sequences alone. Orengo et al. [17] reported that proteins related to the same family could share fewer than 15% identical residues. The protein structure retains a significant portion of similarity even between distant homologs. In general, the degree of structural or sequence similarity varies substantially between protein families. Some families can handle more changes than others. This so-called structural plasticity [17] has a considerable impact on the functionality of some proteins, or members of a protein family. A consideration of the protein structure and its variability becomes important in such situations for further analysis of functional similarity between proteins.A classical approach for deriving the protein function is by first determining its 3D structure, which can then provide some ideas about its function [17]. Protein 3D structures provide information about the binding sites, active sites, and how proteins interact with each other, and thus could provide an insight into the function of the protein [17]. How proteins interact with each other and with other molecules (e.g., ligands) is determined primarily by the amino acids on the protein surface [18]. Therefore, knowledge of the protein surface residues could help in a better understanding of what molecules are binding together, and in some cases, why they bind [18]. The protein surface could also provide significant information about protein functions which cannot be easily detected, even in the presence of sequence or fold similarity. Therefore, the analysis of protein surfaces is important in the study of intermolecular interactions. Clearly, advances in our understanding of protein surfaces could have important implications in various biomedical fields, such as personalized medicine, drug discovery, drug design, and so forth.


### 2.2. Protein Surface Characterization Methods

Given the foregoing, it is not surprising that different methods have been proposed to characterize the protein surface.  Popular examples include those based on surface shape distributions [19], Gauss integral [20], Fourier transform [21], spherical harmonics [22, 23], alpha-shapes [2, 24], and Zernike polynomials [7]. Contact maps between protein surfaces were studied in  [25], while similarity networks between surface patches from protein binding sites were studied in [26, 27].  Protein surface similarity using varying resolutions of structural data have also been studied, for instance, using medium-resolution Cryo-EM maps in [26] and low resolution protein structure data in [28].  SHARP [29] provides a mechanism to predict protein-protein interaction by analyzing overlapping protein 3D surface regions. SURFACE [5] is a database of protein surface regions that can be useful for annotation.Much earlier, Jones and Thornton [30] analyzed protein-protein interaction by using surface patches, where patches are defined based on the Cα atoms that have a predetermined accessible surface area, and adhere to defined constraints on the solvent vectors. Each patch is then described using six parameters, namely, solvation potential, residue interface propensity, hydrophobicity, planarity, protrusion, and accessible surface area. Ferrè et al. [5] analyzed locally similar structures by matching surface patches composed of subsets of amino acids. Each residue on the protein surface is represented using a vector joining its Cα atom and the centroid of its side chain atoms. Surface patches are then compared for similarity by comparing the residue vectors for all possible pairs of residues from the query and target surface patch. Matches are determined based on the root mean square distance, and the residue similarity as determined using a standard substitution matrix. The results of using this method on a nonredundant list of protein chains as recorded in the SURFACE database [5], a collection of protein surface regions that can be useful for annotation. Below, we describe three approaches that are more closely related to our work. See [28, 31] for reviews on surface comparison methods.Distance Distributions
Distances, geometry, and topology have for long been used in the analysis of general protein 3D structures [32]. Residue distances have been used in standard texture-based analysis of 2D textures (distance matrices) formed by the distances between residues in a protein structure [33].  The use of topological invariants, as captured using Gaus integrals for the automated analysis and representation of general protein 3D structures was described in [20]. Much earlier, Connolly [19] proposed the analysis of protein surfaces using the notion of surface shape distributions. Essentially, surface shapes correspond to different geometric configurations defined on the protein surface. Binkowski and Joachimiak [11] proposed the use of surface shape signatures (SSSs) as a method to describe protein surfaces by exploiting global shape and geometrical properties of the surfaces. Shape signatures are computed based on the distances measured between each unique atom pairs on the surface. Distances are then sorted based on which their distributions are generated. With the distributions, the problem of matching between two surfaces is now reduced to that of comparing their distributions. Comparison between two distributions is performed using the Kolmogorov-Smirnov (KS) test. The use of the shape distribution is fast and relatively resilient to scale, rotation, and mirroring. However, the discrimination ability is still a problem, as the SSS tends to lose important surface details.Zernike Polynomials
Following earlier work by Canterakis [34] on the use of 3D Zernikes for the analysis of general 3D objects, Sael et al. [7] introduced 3D Zernike to the area of protein structural similarity matching. Here, the protein 3D structure is represented as a series expansion of 3D Zernike functions. The triangulated Connolly surface of the protein is computed, and subsequently the protein is placed into a 3D cubic grid and voxelized. Each voxel has a value of 1 or 0, depending on whether the voxel is on the protein surface or in the interior. The 3D Zernike function is then applied to the voxelized 3D protein shape to obtain the 3D Zernike descriptors. Therefore, the problem of comparison of 3D surfaces is reduced to that of comparing two vectors representing the 3D Zernike descriptors for each protein surface. Several distance measures were tried, such as the Euclidean distance, Manhattan distance, and a correlation-based distance defined as the complement of the correlation coefficient between two Zernike descriptors. Venkatraman et al. [23] studied the use of both spherical harmonics and 3D Zernike descriptors in the retrieval of functionally similar proteins. In a more recent work, Sael and Kihara [28] used the Zernike descriptor to study protein surfaces in low resolution data. Computation of the required Zernike polynomials is, however, known to be a major computational huddle [35]. This problem is even worse for the 3D Zernike polynomials needed for protein surfaces. Thus, the required preprocessing before matching is performed may be a problem for indexing and real-time search of large-scale datasets.Fingerprints
A recent work [8] used the idea of extracting invariant fingerprints from patches on the protein surface. Patches are obtained by generating the dot surface of the protein and constructing a graph to approximate the protein surface.  Afterwards, circular patches are generated as a contiguous surface area from a center point, where the radius of the patch is within a predetermined cutoff. Patches are created for each single point on the surface, after which a fingerprint representation of the patch is computed as a geodesic distance-dependent distribution of directional curvature. Geodesic distances are computed from the central vertex in each patch. Comparisons between fingerprints were performed using the average fingerprint similarity score (AFSS) and the direct fingerprint similarity score (DFSS). Final scores are computed after an alignment procedure based on the AFSS. Clearly, computational complexity will be a major problem here, especially given the computation of the patch representation for each vertex on the surface graph (number of vertices is much more than the number of surface residues). The need for a later stage of alignment for the final computation of matching scores only compounds the computational burden (see [12], for example).The key difference in our method is the use of the local patch descriptors as defined by the distribution of distances between Cα atoms within each surface patch, conditioned on the specific residue at the center of the patch, and the particular residues found within the patch. Our method computes the residue-specific distance distributions, and residue-distance cooccurrences for the protein surface patches using only the Cα atoms on the protein surface. Residues in the interior of the protein are discarded. Unlike the approach in [8], we avoid the time complexity of generating a graph representation of the surface before the surface can be scanned to generate the patches and then compute the distance distribution. Further, ours does not depend on the time-consuming process of initial surface alignment.


## 2.1. Protein Sequences, Structure, and Surface

Although proteins could vary significantly in their functions and 3D shapes, they also share a general common structure. Proteins are composed of 20-amino acids that are connected via peptide bonds [13]. Each protein is composed of an ordered sequence of amino acids. The order in which these amino acids are connected is called the protein sequence, or the primary structure of the protein [14]. This primary sequence determines the 3D structure of the protein. All proteins are composed of four common structural types: primary structure, secondary structure, tertiary structure, and quaternary structure. The primary structure is simply the amino acid sequence. The secondary structure is formed by patterns of intermolecular bonding of hydrogen and is determined primarily by the location and the directions of these patterns [14, 15]. This is often described in terms of secondary structural elements (SSEs), such as α-helixes, β-sheets, and turns. The overall 3D shape of the secondary structures determines the tertiary structure of the protein. When two or more chains combine to form a larger molecule, the whole structure is called the quaternary structure. Figure 1 shows an example of some of the common protein structural types (the sequence is not included).Protein structures for a sample protein (PDB id: 2UDI). (a) Secondary structure elements—α-helixes (magenta), β-sheets (gold), and turns (gray); (b) two chains: chain E (blue), chain I (green);  (c) surface and 3D shape for chain E; (d) surface and 3D shape for chain I; (e) quaternary structure for the protein. Figures are produced using PMV [16].
(a)(b)(c)(d)(e)A common method for protein function prediction is by annotation transfer from known homologous proteins [17]. Functions of novel proteins can be determined by sequence comparisons, for instance using sequence alignment. When proteins evolve, the protein structure remains more highly conserved when compared to the sequence. Protein sequences change more easily during evolution due to residue mutations, for instance by substitution, insertion, or deletion. Hence, proteins that belong to the same family (homologous proteins) may not be identified using sequences alone. Orengo et al. [17] reported that proteins related to the same family could share fewer than 15% identical residues. The protein structure retains a significant portion of similarity even between distant homologs. In general, the degree of structural or sequence similarity varies substantially between protein families. Some families can handle more changes than others. This so-called structural plasticity [17] has a considerable impact on the functionality of some proteins, or members of a protein family. A consideration of the protein structure and its variability becomes important in such situations for further analysis of functional similarity between proteins.A classical approach for deriving the protein function is by first determining its 3D structure, which can then provide some ideas about its function [17]. Protein 3D structures provide information about the binding sites, active sites, and how proteins interact with each other, and thus could provide an insight into the function of the protein [17]. How proteins interact with each other and with other molecules (e.g., ligands) is determined primarily by the amino acids on the protein surface [18]. Therefore, knowledge of the protein surface residues could help in a better understanding of what molecules are binding together, and in some cases, why they bind [18]. The protein surface could also provide significant information about protein functions which cannot be easily detected, even in the presence of sequence or fold similarity. Therefore, the analysis of protein surfaces is important in the study of intermolecular interactions. Clearly, advances in our understanding of protein surfaces could have important implications in various biomedical fields, such as personalized medicine, drug discovery, drug design, and so forth.


## 2.2. Protein Surface Characterization Methods

Given the foregoing, it is not surprising that different methods have been proposed to characterize the protein surface.  Popular examples include those based on surface shape distributions [19], Gauss integral [20], Fourier transform [21], spherical harmonics [22, 23], alpha-shapes [2, 24], and Zernike polynomials [7]. Contact maps between protein surfaces were studied in  [25], while similarity networks between surface patches from protein binding sites were studied in [26, 27].  Protein surface similarity using varying resolutions of structural data have also been studied, for instance, using medium-resolution Cryo-EM maps in [26] and low resolution protein structure data in [28].  SHARP [29] provides a mechanism to predict protein-protein interaction by analyzing overlapping protein 3D surface regions. SURFACE [5] is a database of protein surface regions that can be useful for annotation.Much earlier, Jones and Thornton [30] analyzed protein-protein interaction by using surface patches, where patches are defined based on the Cα atoms that have a predetermined accessible surface area, and adhere to defined constraints on the solvent vectors. Each patch is then described using six parameters, namely, solvation potential, residue interface propensity, hydrophobicity, planarity, protrusion, and accessible surface area. Ferrè et al. [5] analyzed locally similar structures by matching surface patches composed of subsets of amino acids. Each residue on the protein surface is represented using a vector joining its Cα atom and the centroid of its side chain atoms. Surface patches are then compared for similarity by comparing the residue vectors for all possible pairs of residues from the query and target surface patch. Matches are determined based on the root mean square distance, and the residue similarity as determined using a standard substitution matrix. The results of using this method on a nonredundant list of protein chains as recorded in the SURFACE database [5], a collection of protein surface regions that can be useful for annotation. Below, we describe three approaches that are more closely related to our work. See [28, 31] for reviews on surface comparison methods.Distance Distributions
Distances, geometry, and topology have for long been used in the analysis of general protein 3D structures [32]. Residue distances have been used in standard texture-based analysis of 2D textures (distance matrices) formed by the distances between residues in a protein structure [33].  The use of topological invariants, as captured using Gaus integrals for the automated analysis and representation of general protein 3D structures was described in [20]. Much earlier, Connolly [19] proposed the analysis of protein surfaces using the notion of surface shape distributions. Essentially, surface shapes correspond to different geometric configurations defined on the protein surface. Binkowski and Joachimiak [11] proposed the use of surface shape signatures (SSSs) as a method to describe protein surfaces by exploiting global shape and geometrical properties of the surfaces. Shape signatures are computed based on the distances measured between each unique atom pairs on the surface. Distances are then sorted based on which their distributions are generated. With the distributions, the problem of matching between two surfaces is now reduced to that of comparing their distributions. Comparison between two distributions is performed using the Kolmogorov-Smirnov (KS) test. The use of the shape distribution is fast and relatively resilient to scale, rotation, and mirroring. However, the discrimination ability is still a problem, as the SSS tends to lose important surface details.Zernike Polynomials
Following earlier work by Canterakis [34] on the use of 3D Zernikes for the analysis of general 3D objects, Sael et al. [7] introduced 3D Zernike to the area of protein structural similarity matching. Here, the protein 3D structure is represented as a series expansion of 3D Zernike functions. The triangulated Connolly surface of the protein is computed, and subsequently the protein is placed into a 3D cubic grid and voxelized. Each voxel has a value of 1 or 0, depending on whether the voxel is on the protein surface or in the interior. The 3D Zernike function is then applied to the voxelized 3D protein shape to obtain the 3D Zernike descriptors. Therefore, the problem of comparison of 3D surfaces is reduced to that of comparing two vectors representing the 3D Zernike descriptors for each protein surface. Several distance measures were tried, such as the Euclidean distance, Manhattan distance, and a correlation-based distance defined as the complement of the correlation coefficient between two Zernike descriptors. Venkatraman et al. [23] studied the use of both spherical harmonics and 3D Zernike descriptors in the retrieval of functionally similar proteins. In a more recent work, Sael and Kihara [28] used the Zernike descriptor to study protein surfaces in low resolution data. Computation of the required Zernike polynomials is, however, known to be a major computational huddle [35]. This problem is even worse for the 3D Zernike polynomials needed for protein surfaces. Thus, the required preprocessing before matching is performed may be a problem for indexing and real-time search of large-scale datasets.Fingerprints
A recent work [8] used the idea of extracting invariant fingerprints from patches on the protein surface. Patches are obtained by generating the dot surface of the protein and constructing a graph to approximate the protein surface.  Afterwards, circular patches are generated as a contiguous surface area from a center point, where the radius of the patch is within a predetermined cutoff. Patches are created for each single point on the surface, after which a fingerprint representation of the patch is computed as a geodesic distance-dependent distribution of directional curvature. Geodesic distances are computed from the central vertex in each patch. Comparisons between fingerprints were performed using the average fingerprint similarity score (AFSS) and the direct fingerprint similarity score (DFSS). Final scores are computed after an alignment procedure based on the AFSS. Clearly, computational complexity will be a major problem here, especially given the computation of the patch representation for each vertex on the surface graph (number of vertices is much more than the number of surface residues). The need for a later stage of alignment for the final computation of matching scores only compounds the computational burden (see [12], for example).The key difference in our method is the use of the local patch descriptors as defined by the distribution of distances between Cα atoms within each surface patch, conditioned on the specific residue at the center of the patch, and the particular residues found within the patch. Our method computes the residue-specific distance distributions, and residue-distance cooccurrences for the protein surface patches using only the Cα atoms on the protein surface. Residues in the interior of the protein are discarded. Unlike the approach in [8], we avoid the time complexity of generating a graph representation of the surface before the surface can be scanned to generate the patches and then compute the distance distribution. Further, ours does not depend on the time-consuming process of initial surface alignment.


## 3. Methods

We present an invariant descriptor for characterizing protein surfaces. We start with a local descriptor of selected circular patches on the protein surface. For a given surface patch, the local descriptor is computed based on the residue distances from the center of the patch. The descriptor records the distance distribution between the central residue and the residues within the patch, keeping track of the number of particular pairwise residue cooccurrences in the patch. A global descriptor for the entire protein surface is then constructed from the local descriptors by combining information from local descriptors with similar central residues. The proposed descriptor is invariant to rotations of the surface and mirroring.Using a fixed patch size, we obtain a descriptor for the protein surface, independent of the size of the protein structure. Thus, the descriptor can facilitate the rapid matching of protein chains, and will eliminate the need for the exhaustive alignment of the protein 3D structures. For a given protein structure or protein chain from a database, such as the PDB, the proposed method can be summarize in the following steps:(1)
generate the Connolly surface [36] for the protein chain;(2)
generate the surface patches and compute the local invariant descriptor for each patch on the surface;(3)
compute the global invariant surface descriptor for the protein chain, by combining information from the local patch descriptors;(4)
perform surface matching and comparison using the descriptors;(5)
classify the protein into its potential functional family, or perform protein surface retrieval using the invariant descriptors.Figure2 shows a schematic diagram of the general approach. The method has been applied on three protein families: uracil-DNA glycosylase, estrogen receptor, and cell division protein kinase 2. These are the same protein families used in a recently published work [8]. We also tested on epidermal growth factor (EGF) and cyclooxygenase-2  (COX-2), two protein families that are known to play a role in cancer. Below, we provide more details on the steps enumerated above.Figure 2
Schematic diagram for the protein surface characterization using an invariant descriptor. Protein structures in the figure are produced using PMV [3].


### 3.1. Surface Generation

For a given protein, we first generate its Connolly surface [36] at a given atomic radii, using the MSMS program [37], based on which the dot surface is generated. This dot surface is stored in a vertex file. We have used a probe radius of 1.4 Å in all our experiments. Next, MATLAB Bioinformatics Toolbox (Mathworks Inc, Natick, Mass, USA) was used to extract the protein chains and to generate the residue coordinates in each chain. In this step, the chains are extracted while preserving the coordinates of the Cα atoms and their respective residue types by extracting the information from the PDB and the vertex files.


### 3.2. The Invariant Descriptor

#### 3.2.1. Surface Patches

To capture protein structure similarity and to avoid the computational complexity and the time-consuming problem of aligning 3D protein structures, we propose the use of a global rotational-invariant descriptor to represent overlapping patches on the protein surface. A patch is defined as a circular region with a specified radius, centered on the Cα position of a surface residue.  For each residue on the surface (the central residue), we construct a surface patch by recording its residue type, and consider all residues within a certain distance threshold (τp) as part of the patch (see Figure 2). Thus, the proposed surface descriptor is composed of 20 distinct descriptions, one for each protein residue type. For the local descriptor, this is constructed from only information from the patch. For the global descriptor, this is constructed by combining information from patches with the same central residue.The local invariant descriptor for the patch is created by calculating the distribution of distances between the central residue and all other surface residues within the patch. Additionally, the residue cooccurrences within the patch are also recorded as a part of the local descriptor. Each local descriptor is represented in a matrixDA of size (20+1)×(b+1), where the rows correspond to the 20 distinct protein residue types, plus an extra row to describe the summary distance distribution within the patch. The columns represent the individual bins used to capture the distance distributions (total of b bins), plus an extra column to represent the summary of the residue cooccurrences. To reduce the computational time and space requirement, unlike in [8] we define patches only for surface residue positions, rather than for each vertex on the dot surface (the number of vertexes is much more than the number of residues). Therefore, for a given chain, the number of local invariant descriptors will be equal to the number of surface residues. Yet, this number can vary from tens to hundreds and sometimes to thousands of surface residues.  Using a huge number of local invariant descriptors for one chain to perform matching will be very time-consuming. To further reduce the computational requirements, for a given chain, we compute a global rotational-invariant descriptor by combining the 20 distinct residue-specific descriptors. For a given residue type, the global descriptor is constructed by taking the average of all patch descriptors with a given residue type as the central residue (see Figure 2). We consider three ways to represent and use the global surface descriptor, as explained below.


#### 3.2.2. Distance Distribution (DD2)

The basic idea of using the distance distribution is that similar functional proteins should have a similar distribution of distances between the residues on their surfaces. The patch descriptor captures the distribution in two forms. Thefirst form is a detailed distance distribution between the central residue in the surface patch and each of the other residues on the patch. To achieve this, a uniform distribution of the distances is assumed and the total number of bins b is used to estimate the probability distribution of finding a pair of residues at one of the b ranges. The second form is the global distance probability distribution. In this form we estimate the probability of observing any given residue within a patch in a particular distance range from the central residue. In this paper, we study the use of the global distance distribution in identifying similar protein surfaces, and possibly proteins with similar functions. Consequently, the question to be answered is, given a central residue of a specific type, what is the distance distribution for the residues around this central residue? That is, we seek Pr{d∣Rc}, the probability of observing distance d between a central residue of type Rc and any other residue. We expect that the distance distribution should be similar for surface patches from functionally similar proteins.


#### 3.2.3. Residue Cooccurrences (RCs)

Given that surface structures are more conserved than sequence over evolution [15, 17], we expect that functionally similar proteins are likely to have similar surface residues, even though the order of such residues may have changed. This intuition is captured using residue cooccurrences on the protein surface. Using the distance distribution globally provides an idea of how the distances from the central residue are distributed in the protein surface patch. However, there is no constraint on, or indication of, which residues are involved in the formation of these distributions. The co-occurrence of a given residue with the central residue is calculated as the number of times the residue occurs on a patch with the same central residue. Thus, the main problem would be to find the probability of observing residue say, Ri, given a central residue, say Rc. Again, we expect the probability Pr{Ri∣Rc}, to be similar for protein surfaces from functionally similar proteins. We note that the surface co-occurrence does not depend on the specific distance between the residues involved, as far as Ri is within the patch.


#### 3.2.4. Distance-Residue Cooccurrences (DRCs)

The above have considered the distance and the co-occurrence separately. The DRC combines the general distance distribution (represented as a row vector,sumC in matrix DA) and the residue cooccurrences (represented as a column vector, sumR in matrix DA) in describing the protein surface (see Figure 2). The residue-distance co-occurrence vector is defined as follows: DRC=(sumC°sumRT), where ° is the concatenation operator and XT stands for the transpose of X. DRC is used to compute the conditional probability  Pr{d∣Rc,Ri}, that is, the probability of observing the distance d between residues Rc and Ri given that Rc is the central residue in the patch. We expect that the residue co-occurrence (sumR, or RC) should carry more distinctive functionally relevant information than the general distance distribution (sumC, or DD2), since surface residue cooccurrences are likely to be more conserved over evolution. By combining both vectors, we can account for both the geometry of the protein surface and the distribution of specific residues within specific distances on the surface. Using both vectors brings in some biological relevance in the analysis and is likely to lead to improved results in the identification of functionally similar protein surfaces.


### 3.3. Matching and Classification

Given two proteins, say Protein 1 and Protein 2 we characterize them using their global descriptors, sayDg1 and Dg2 respectively. In this work, the global descriptor could be the distance distribution (DD2), residue cooccurrences (RCs), or the distance-residue cooccurrences (DRCs).Distance Distribution
For matching using the distance distribution we create a vectorDg1d that is composed of the 20 global distance distributions represented by all sumC vectors from each descriptor. Dg1d is defined as Dg1d=(Dd1°Dd2°⋯°Dd20), where Dd1,Dd2,…,Dd20 are the distance distributions from each residue type on the surface of Protein 1. Repeat the same process for Protein 2 to create Dg2d. Then we perform matching using the simple Euclidean distance: D12=∑i=1n[Dg1d(i)-Dg2d(i)]2.Residue Cooccurrences
For Protein 1 we create a vectorDg1c that combines the 20 residues co-occurrence vectors (denoted sumR), defined as Dg1c=(Dc1T°Dc2T°⋯°Dc20T), where Dc1T,Dc2T,…,Dc20T represents sumR1T, sumR2T, and sumR20T. Similarly, we compute Dg2c. Matching is performed using the Euclidean distance between Dg1c and Dg2c.Distance-Residue Cooccurrences
Here, we create a vectorDRC that is comprised of all of the distance distributions as well as the residue cooccurrences. For Protein 1, we have DRC1=(Dd1°Dc1T°Dd2°Dc2T°⋯°Dd1°Dc1T). Similarly we obtain DRC2 for Protein 2. Again for simplicity, matching is performed using the Euclidean distance. Clearly, other distance measures could be used.Classification
Having computed the surface descriptors and the distance between protein surfaces using the descriptors, one may be interested in determining whether a given unknown protein belongs to some known protein family. Using some training data, we can compute surface descriptors for the known family, and based on these perform the required classification. Classification is performed using Weka [38, 39], an open-source software for machine leaning that provides a suite of classification algorithms.


## 3.1. Surface Generation

For a given protein, we first generate its Connolly surface [36] at a given atomic radii, using the MSMS program [37], based on which the dot surface is generated. This dot surface is stored in a vertex file. We have used a probe radius of 1.4 Å in all our experiments. Next, MATLAB Bioinformatics Toolbox (Mathworks Inc, Natick, Mass, USA) was used to extract the protein chains and to generate the residue coordinates in each chain. In this step, the chains are extracted while preserving the coordinates of the Cα atoms and their respective residue types by extracting the information from the PDB and the vertex files.


## 3.2. The Invariant Descriptor

### 3.2.1. Surface Patches

To capture protein structure similarity and to avoid the computational complexity and the time-consuming problem of aligning 3D protein structures, we propose the use of a global rotational-invariant descriptor to represent overlapping patches on the protein surface. A patch is defined as a circular region with a specified radius, centered on the Cα position of a surface residue.  For each residue on the surface (the central residue), we construct a surface patch by recording its residue type, and consider all residues within a certain distance threshold (τp) as part of the patch (see Figure 2). Thus, the proposed surface descriptor is composed of 20 distinct descriptions, one for each protein residue type. For the local descriptor, this is constructed from only information from the patch. For the global descriptor, this is constructed by combining information from patches with the same central residue.The local invariant descriptor for the patch is created by calculating the distribution of distances between the central residue and all other surface residues within the patch. Additionally, the residue cooccurrences within the patch are also recorded as a part of the local descriptor. Each local descriptor is represented in a matrixDA of size (20+1)×(b+1), where the rows correspond to the 20 distinct protein residue types, plus an extra row to describe the summary distance distribution within the patch. The columns represent the individual bins used to capture the distance distributions (total of b bins), plus an extra column to represent the summary of the residue cooccurrences. To reduce the computational time and space requirement, unlike in [8] we define patches only for surface residue positions, rather than for each vertex on the dot surface (the number of vertexes is much more than the number of residues). Therefore, for a given chain, the number of local invariant descriptors will be equal to the number of surface residues. Yet, this number can vary from tens to hundreds and sometimes to thousands of surface residues.  Using a huge number of local invariant descriptors for one chain to perform matching will be very time-consuming. To further reduce the computational requirements, for a given chain, we compute a global rotational-invariant descriptor by combining the 20 distinct residue-specific descriptors. For a given residue type, the global descriptor is constructed by taking the average of all patch descriptors with a given residue type as the central residue (see Figure 2). We consider three ways to represent and use the global surface descriptor, as explained below.


### 3.2.2. Distance Distribution (DD2)

The basic idea of using the distance distribution is that similar functional proteins should have a similar distribution of distances between the residues on their surfaces. The patch descriptor captures the distribution in two forms. Thefirst form is a detailed distance distribution between the central residue in the surface patch and each of the other residues on the patch. To achieve this, a uniform distribution of the distances is assumed and the total number of bins b is used to estimate the probability distribution of finding a pair of residues at one of the b ranges. The second form is the global distance probability distribution. In this form we estimate the probability of observing any given residue within a patch in a particular distance range from the central residue. In this paper, we study the use of the global distance distribution in identifying similar protein surfaces, and possibly proteins with similar functions. Consequently, the question to be answered is, given a central residue of a specific type, what is the distance distribution for the residues around this central residue? That is, we seek Pr{d∣Rc}, the probability of observing distance d between a central residue of type Rc and any other residue. We expect that the distance distribution should be similar for surface patches from functionally similar proteins.


### 3.2.3. Residue Cooccurrences (RCs)

Given that surface structures are more conserved than sequence over evolution [15, 17], we expect that functionally similar proteins are likely to have similar surface residues, even though the order of such residues may have changed. This intuition is captured using residue cooccurrences on the protein surface. Using the distance distribution globally provides an idea of how the distances from the central residue are distributed in the protein surface patch. However, there is no constraint on, or indication of, which residues are involved in the formation of these distributions. The co-occurrence of a given residue with the central residue is calculated as the number of times the residue occurs on a patch with the same central residue. Thus, the main problem would be to find the probability of observing residue say, Ri, given a central residue, say Rc. Again, we expect the probability Pr{Ri∣Rc}, to be similar for protein surfaces from functionally similar proteins. We note that the surface co-occurrence does not depend on the specific distance between the residues involved, as far as Ri is within the patch.


### 3.2.4. Distance-Residue Cooccurrences (DRCs)

The above have considered the distance and the co-occurrence separately. The DRC combines the general distance distribution (represented as a row vector,sumC in matrix DA) and the residue cooccurrences (represented as a column vector, sumR in matrix DA) in describing the protein surface (see Figure 2). The residue-distance co-occurrence vector is defined as follows: DRC=(sumC°sumRT), where ° is the concatenation operator and XT stands for the transpose of X. DRC is used to compute the conditional probability  Pr{d∣Rc,Ri}, that is, the probability of observing the distance d between residues Rc and Ri given that Rc is the central residue in the patch. We expect that the residue co-occurrence (sumR, or RC) should carry more distinctive functionally relevant information than the general distance distribution (sumC, or DD2), since surface residue cooccurrences are likely to be more conserved over evolution. By combining both vectors, we can account for both the geometry of the protein surface and the distribution of specific residues within specific distances on the surface. Using both vectors brings in some biological relevance in the analysis and is likely to lead to improved results in the identification of functionally similar protein surfaces.


## 3.2.1. Surface Patches

To capture protein structure similarity and to avoid the computational complexity and the time-consuming problem of aligning 3D protein structures, we propose the use of a global rotational-invariant descriptor to represent overlapping patches on the protein surface. A patch is defined as a circular region with a specified radius, centered on the Cα position of a surface residue.  For each residue on the surface (the central residue), we construct a surface patch by recording its residue type, and consider all residues within a certain distance threshold (τp) as part of the patch (see Figure 2). Thus, the proposed surface descriptor is composed of 20 distinct descriptions, one for each protein residue type. For the local descriptor, this is constructed from only information from the patch. For the global descriptor, this is constructed by combining information from patches with the same central residue.The local invariant descriptor for the patch is created by calculating the distribution of distances between the central residue and all other surface residues within the patch. Additionally, the residue cooccurrences within the patch are also recorded as a part of the local descriptor. Each local descriptor is represented in a matrixDA of size (20+1)×(b+1), where the rows correspond to the 20 distinct protein residue types, plus an extra row to describe the summary distance distribution within the patch. The columns represent the individual bins used to capture the distance distributions (total of b bins), plus an extra column to represent the summary of the residue cooccurrences. To reduce the computational time and space requirement, unlike in [8] we define patches only for surface residue positions, rather than for each vertex on the dot surface (the number of vertexes is much more than the number of residues). Therefore, for a given chain, the number of local invariant descriptors will be equal to the number of surface residues. Yet, this number can vary from tens to hundreds and sometimes to thousands of surface residues.  Using a huge number of local invariant descriptors for one chain to perform matching will be very time-consuming. To further reduce the computational requirements, for a given chain, we compute a global rotational-invariant descriptor by combining the 20 distinct residue-specific descriptors. For a given residue type, the global descriptor is constructed by taking the average of all patch descriptors with a given residue type as the central residue (see Figure 2). We consider three ways to represent and use the global surface descriptor, as explained below.


## 3.2.2. Distance Distribution (DD2)

The basic idea of using the distance distribution is that similar functional proteins should have a similar distribution of distances between the residues on their surfaces. The patch descriptor captures the distribution in two forms. Thefirst form is a detailed distance distribution between the central residue in the surface patch and each of the other residues on the patch. To achieve this, a uniform distribution of the distances is assumed and the total number of bins b is used to estimate the probability distribution of finding a pair of residues at one of the b ranges. The second form is the global distance probability distribution. In this form we estimate the probability of observing any given residue within a patch in a particular distance range from the central residue. In this paper, we study the use of the global distance distribution in identifying similar protein surfaces, and possibly proteins with similar functions. Consequently, the question to be answered is, given a central residue of a specific type, what is the distance distribution for the residues around this central residue? That is, we seek Pr{d∣Rc}, the probability of observing distance d between a central residue of type Rc and any other residue. We expect that the distance distribution should be similar for surface patches from functionally similar proteins.


## 3.2.3. Residue Cooccurrences (RCs)

Given that surface structures are more conserved than sequence over evolution [15, 17], we expect that functionally similar proteins are likely to have similar surface residues, even though the order of such residues may have changed. This intuition is captured using residue cooccurrences on the protein surface. Using the distance distribution globally provides an idea of how the distances from the central residue are distributed in the protein surface patch. However, there is no constraint on, or indication of, which residues are involved in the formation of these distributions. The co-occurrence of a given residue with the central residue is calculated as the number of times the residue occurs on a patch with the same central residue. Thus, the main problem would be to find the probability of observing residue say, Ri, given a central residue, say Rc. Again, we expect the probability Pr{Ri∣Rc}, to be similar for protein surfaces from functionally similar proteins. We note that the surface co-occurrence does not depend on the specific distance between the residues involved, as far as Ri is within the patch.


## 3.2.4. Distance-Residue Cooccurrences (DRCs)

The above have considered the distance and the co-occurrence separately. The DRC combines the general distance distribution (represented as a row vector,sumC in matrix DA) and the residue cooccurrences (represented as a column vector, sumR in matrix DA) in describing the protein surface (see Figure 2). The residue-distance co-occurrence vector is defined as follows: DRC=(sumC°sumRT), where ° is the concatenation operator and XT stands for the transpose of X. DRC is used to compute the conditional probability  Pr{d∣Rc,Ri}, that is, the probability of observing the distance d between residues Rc and Ri given that Rc is the central residue in the patch. We expect that the residue co-occurrence (sumR, or RC) should carry more distinctive functionally relevant information than the general distance distribution (sumC, or DD2), since surface residue cooccurrences are likely to be more conserved over evolution. By combining both vectors, we can account for both the geometry of the protein surface and the distribution of specific residues within specific distances on the surface. Using both vectors brings in some biological relevance in the analysis and is likely to lead to improved results in the identification of functionally similar protein surfaces.


## 3.3. Matching and Classification

Given two proteins, say Protein 1 and Protein 2 we characterize them using their global descriptors, sayDg1 and Dg2 respectively. In this work, the global descriptor could be the distance distribution (DD2), residue cooccurrences (RCs), or the distance-residue cooccurrences (DRCs).Distance Distribution
For matching using the distance distribution we create a vectorDg1d that is composed of the 20 global distance distributions represented by all sumC vectors from each descriptor. Dg1d is defined as Dg1d=(Dd1°Dd2°⋯°Dd20), where Dd1,Dd2,…,Dd20 are the distance distributions from each residue type on the surface of Protein 1. Repeat the same process for Protein 2 to create Dg2d. Then we perform matching using the simple Euclidean distance: D12=∑i=1n[Dg1d(i)-Dg2d(i)]2.Residue Cooccurrences
For Protein 1 we create a vectorDg1c that combines the 20 residues co-occurrence vectors (denoted sumR), defined as Dg1c=(Dc1T°Dc2T°⋯°Dc20T), where Dc1T,Dc2T,…,Dc20T represents sumR1T, sumR2T, and sumR20T. Similarly, we compute Dg2c. Matching is performed using the Euclidean distance between Dg1c and Dg2c.Distance-Residue Cooccurrences
Here, we create a vectorDRC that is comprised of all of the distance distributions as well as the residue cooccurrences. For Protein 1, we have DRC1=(Dd1°Dc1T°Dd2°Dc2T°⋯°Dd1°Dc1T). Similarly we obtain DRC2 for Protein 2. Again for simplicity, matching is performed using the Euclidean distance. Clearly, other distance measures could be used.Classification
Having computed the surface descriptors and the distance between protein surfaces using the descriptors, one may be interested in determining whether a given unknown protein belongs to some known protein family. Using some training data, we can compute surface descriptors for the known family, and based on these perform the required classification. Classification is performed using Weka [38, 39], an open-source software for machine leaning that provides a suite of classification algorithms.


## 4. Results and Discussion

### 4.1. Datasets and Environment

We performed experiments to test the performance of the proposed protein surface descriptor in two protein structure analysis tasks, namely, classifying proteins into their most likely functional groups, and ranking and retrieval of protein surfaces. We used two datasets for the experiments. DATASET-A contained information from three protein families:uracil-DNA glycosylase, cell division protein kinase 2, and estrogen receptor. This was created by scanning the PDB and selecting the protein structures with protein chains belonging to one of the three families. We were able to extract 416 chains that belong to 243 proteins in the PDB. The dataset is distributed as follows: 91 chains from 46 distinct proteins for uracil-DNA glycosylase (Group1), 186 chains from 95 distinct proteins for estrogen receptor family (Group2), and 139 chains from 102 distinct proteins from cell division protein kinase 2 (Group3). We used DATASET-A basically to train the system, and perform initial testing.  DATASET-B contained protein structures from two families, namely cyclooxygenase-2  (COX-2) (51 proteins, 95 chains) and epidermal growth factor  (EGF) (67 proteins, 71 chains). We then extracted protein structures from the PDB that have 10 or less chains and ignored the rest. This resulted in a total of 15,386 protein chains form 6,261 unique proteins. DATASET-B included all structures in DATASET-A. We used DATASET-B for a more comprehensive scan of the PDB, in the quest for potentially novel structures that may be related to the two families. Experiments were performed using a SONY VAIO personal computer, with Intel Core 2 Duo Processor T8100, running at 2.10 GHz, with 2 GB of main memory. Programs were written using Matlab (Mathworks Inc, Natick, Mass, USA) with the Bioinformatics Toolbox. We set probe radius = 1.4 Å and patch distance threshold τp=10 Å. For distance distributions, we used a fixed number of bins, b=5. Classification was performed based on algorithms implemented in Weka [38, 39] version 3-6-4.


### 4.2. Classification Performance

We divide DATASET-A into training and testing sets and apply different classifiers on the different descriptors proposed. In all our experiments, the training sets were kept very separate from the testing sets, with no overlap between the two. Classification performance is measured in terms of classification rate based on the three protein families in the dataset. We tested the method using various classifiers implemented in Weka, such as Naïve Bayes, logistic regression, and simple logistic classifier. We report results mainly for the logistic regression.  First, we explore the impact of the size of the testing set and of the training set on the classification performance using the proposed approach. We varied the size of the training set (from 50 to 300), while keeping the size of the testing set fixed. We then checked the performance using fixed testing sets of size 100, 200, and 300. Figure3 shows the results.Variation of classification rate (CR) with size of training set using the proposed descriptors DD2 (a), RC (b), and DRC (c).  Results are shown for the average over 10 runs, using logistic regression as the classifier.
(a)(b)(c)The figure shows that applying the distance distribution (DD2) alone resulted in the lowest performance accuracy as compared to using the residue cooccurrences (RCs) or distance-residue cooccurrences (DRCs). Yet, our definition of the distance distribution shows encouraging results. A steady improvement in performance with increasing training set size can be observed when using DD2 alone, peaking at about 87%  with a training size of 200 and testing size of 100. The distinctiveness of our approach is the use of residue cooccurrences on the protein surface. This approach assumes that functionally similar surface proteins have similar residue cooccurrences within a small local surface region. Figure3 (middle plot) shows that classification using residue cooccurrences (RCs) provided a significant improvement in the classification rate.  A similar improvement was observed using other classifiers, such as Naïve Bayes. Using the RC descriptor, we can achieve an accuracy rate of 94% using a small training set (50 samples) and six times larger testing set (300 chains). This shows the robustness of the residue cooccurrences, even when using a few training samples. We observe that the performance using DD2 was not as robust (about 81% using small training set, peaking at about 87% using 200 training samples).The use of distance-residue co-occurrence presents a steadier improvement in the classification rate. Using the DRC raised the accuracy rate to 99% using the simple logistic classifier on a training set of 150 and testing set of 100 (data not shown).  We can observe the significant difference between the results of DD2 (which did not use information on residue cooccurrences) and RC and/or DRC (both of which used residue cooccurrences). Figure4 shows a corresponding performance measurement with varying size of the testing set, while keeping the training set size fixed. As expected, there is a general slight decrease in performance with increasing size of the test set. The case of DRC using a training set size of 100 seemed to increase slightly with increasing testing set size. The increase is however within a small range (from 0.91 to 0.93). This shows a steady performance over increasing size of the testing set. Overall trends are similar to Figure 3, with RC and DRC performing much better than DD2.  Similar trends were also observed using other classification algorithms.  The overall classification performance is summarized in Figure 5, which shows the results of the three proposed schemes using n-fold cross validation, for different values of n.Variation of classification rate (CR) with size of testing set using the proposed descriptors DD2 (a), RC (b), and DRC (c).  Results are shown for the average over 10 runs, using logistic regression as the classifier.
(a)(b)(c)Figure 5
Summary classification performance usingn-fold cross validation (the x-axis is for varying n).


### 4.3. Ranking and Retrieval

In this section, we explore the effectiveness of our approach on the problem of search and retrieval of protein surfaces. Given a query protein, we study whether our approach has the robustness to place most of the functionally similar proteins in the top hits of the retrieved surfaces.  Here, a query protein from each of the three groups is used to screen the entire DATASET-A (416 samples) and provide a ranking based on the similarity. Thus, each protein structure is ranked against the query, (from 1 to 416), where a lower rank (smaller distance) implies more similarity to the query. After that, we search over the retrieved proteins to find which ranks the functionally similar proteins (i.e., proteins in the same functional group) have attained. Table1 shows the ranking produced using the proposed descriptor, for three query samples, one for each group. Results are shown only for DRC. RC produced a slightly better ranking (especially for uracil-DNA glycosylase family (Group 1)), while DD2 was worse than both RC and DRC). Overall, for Group 2 and Group 3, the Top 30 ranked proteins belonged to the corresponding family, while Group 1 was more difficult.Table 1
Ranking the screened proteins according to their similarity to the query protein. Results are shown for the Top 30 hits for one query protein from each of the three groups, using DATASET-A (416 protein chains).
(a)
DRC on query protein 1UDI chain I (Group 1)
Protein PDB IDChainRankDistance1UDII102ZHXB22.13061LQMB32.35091LQGC42.45892UUGC52.53531EUIC72.59202ZHXL82.61041UGHI102.63492UGIA152.68091UGIE162.68722ZHXH192.69692ZHXD212.70062ZHXN222.70172ZHXJ232.71411UGIG252.72611EMJA422.77582BOOA452.78081UGID472.78682OWRB502.79522J8XD612.81291LQGD702.82631Q3FA902.85332UUGD992.86751UGIA1012.86892OWRC1102.87492ZHXF1162.8852OWQB1292.89771SSPE1412.91151UGIC1422.91172ZHXA1472.915(b)
DRC on query protein 1QKN chain A (Group 2)
Protein PDB IDChainRankDistance1QKNA102J7XA21.37692J7YA31.57531QKMA41.67931NDEA51.73682GIUA61.73711L2IA71.74601U3RB81.76703ERDA91.76833OS9A101.77152IOGA111.77383LTXC121.77421YIMA131.78543ERTA141.79661U3QD151.80091YY4A161.80902OUZA171.81261YINA181.81521XP6A191.82602AYRA201.82693OS8D211.83112QH6A221.83123OSAA231.83671L2JA241.83852JJ3A251.84381G50A261.84903OS8A271.85092FSZA281.85182QGWA291.86041UOMA301.8683(c)
DRC on query protein 1YKR chain A (Group 3)
Protein PDB IDChainRankDistance1YKRA102UZOA21.03962R3OA31.08103PXYA41.08463PY1A51.09402WMAA61.13892IW6A71.16093NS9A81.20433IGGA91.21412WFYA101.21792C5YA111.22703DDPA121.22802J9MA131.22842R3JA141.23742R3LA151.24023PXRA161.24222DUVA171.25341W8CA181.25863DOGA191.27932V22A201.28222R3PA211.28832V22C221.29633IG7A231.32072JGZA241.32752R64A251.33032WHBA261.33812VTNA271.34583LFNA281.34762WIPA291.35142BKZA301.3580We further measured the performance of our approach using the enrichment plot. The enrichment plot essentially measures how well a given ranking or retrieval system performs, when compared with a random selection of the data samples. At a given percentage of database screening, the enrichment factor is computed as the ratioNobs/Nexp, where Nobs=number of functionally similar proteins observed or retrieved by the system, and Nexp=number of functionally similar proteins expected by random selection. For an effective system, we expect that most of the functionally similar proteins should be observed after a small percentage of screening. That is, the top hits should contain mainly functionally similar proteins, and hence the enrichment factor should be high after a small percentage probe of the database, and gradually decrease towards 1 (which corresponds to random selection). Figure 6(a) shows a plot of the average enrichment factor using 5 queries from Group 3. The enrichment plot shows that our proposed method provides better results as we screen a small percentage of the dataset. In most of the cases, our method retrieved about three times better than the expected random retrieval in the first 10% of screened proteins. As we increase the percent of screening, the retrieval degrades, since we are more likely to have retrieved most, if not all of the similar proteins after a small percentage of the screening. Thus, subsequent retrievals will lead to spurious results.Ranking and retrieval performance for the proposed methods. (a) Enrichment plot for screening protein structures using the proposed descriptors. Results are average for 5 query proteins fromcell division protein kinase 2 family (Group 3), using DATASET-A (416 protein chains).  (b) Average precision and recall for three queries, one for each group in DATASET-A. DD1 corresponds to the distance distribution proposed in [11], as described in Section 2 (see Section 4.5 on comparison with related methods).
(a)(b)


### 4.4. Screening Protein Surfaces in PDB

Encouraged by the results in classification and ranking using the proposed descriptors, we now performed a larger scale experiment, by screening the entire protein structures in PDB, using the protein chains in DATASET-B,with members of the COX-2 and EGF families as the query. The main objective was to see how the proposed descriptors will perform on a large scale, and to see if the methods could predict potentially novel functional linkages between any of the families and other proteins in PDB. For this task, we used only PDB files with 10 or less chains, and ignored the rest. This resulted in a total of 15,386 protein chains from 6,261 unique proteins. Table 2(a) shows the ranking results produced by screening the PDB files based on the proposed descriptors, using a member of the EGF family as a query. Table 2(b) shows corresponding results using a member of COX-2 family. Results are shown only for the DRC descriptor. Generally, similar results were obtained using RC. We can notice that some of the unknown proteins (annotated as “uncharacterized”) were placed in the Top-50 rankings, implying a possible relationship with the respective families.Table 2(a)
Top 50 hits using DRC for a query protein structure from the EGF family on DATASET-B. Annotations in bold correspond to members of the EGF family, predicted proteins, or uncharacterized proteins
ProteinChainDistanceProtein name annotationRank2a2qL0.0000COAGULATION FACTOR VII12firL1.4925COAGULATION FACTOR VII LIGHT CHAIN22zp0L1.5545FACTOR VII LIGHT CHAIN31wtgL1.5628COAGULATION FACTOR VII41wunL1.5844COAGULATION FACTOR VII52b8oL1.6305COAGULATION FACTOR VII LIGHT CHAIN62zwlL1.6379FACTOR VII LIGHT CHAIN72zzuL1.6536FACTOR VII LIGHT CHAIN81wqvL1.6816COAGULATION FACTOR VII92ec9L1.6832COAGULATION FACTOR VII101danL1.7655BLOOD COAGULATION FACTOR VIIA111wssL1.7659COAGULATION FACTOR VII122puqL1.7692COAGULATION FACTOR VII131fakL1.7934PROTEIN (BLOOD COAGULATION FACTOR VIIA)142b7dL1.8024COAGULATION FACTOR VII156acnA1.8061ACONITASE162aerL1.8120COAGULATION FACTOR VII172aeiL1.8164COAGULATION FACTOR VII182flrL1.8196COAGULATION FACTOR VII193elaL1.8668COAGULATION FACTOR VII LIGHT CHAIN201z6jL1.8859COAGULATION FACTOR VII212f9bL1.9027COAGULATION FACTOR VII223phsA1.9169CELL WALL SURFACE ANCHOR FAMILY PROTEIN233n54B1.9263SPORE GERMINATION PROTEIN B3243qbpB1.9264FUMARASE FUM253ma9L1.9367TRANSMEMBRANE GLYCOPROTEIN263mt0A1.9921UNCHARACTERIZED PROTEIN PA1789273lguA2.0004PROTEASE DEGS283m7iA2.0071TRANSKETOLASE291qfkL2.0096PROTEIN (COAGULATION FACTOR VIIA (LIGHT CHAIN))303n9tA2.0169PNPC313pxzA2.0235CELL DIVISION PROTEIN KINASE 2322flbL2.0257COAGULATION FACTOR VII333lh1A2.0289PROTEASE DEGS343nlcA2.0300UNCHARACTERIZED PROTEIN VP0956353no5C2.0302UNCHARACTERIZED PROTEIN363msqC2.0347PUTATIVE UBIQUINONE BIOSYNTHESIS PROTEIN373rykA2.0389DTDP-4-DEHYDRORHSMNOSE 3,5-EPIMERASE383m4aA2.0444PUTATIVE TYPE I TOPOISOMERASE393n3nB2.0444CATALASE-PEROXIDASE402R3GA2.0456CELL DIVISION PROTEIN KINASE 2412R3IA2.0465CELL DIVISION PROTEIN KINASE 2423o0rL2.0473ANTIBODY FAB FRAGMENT LIGHT CHAIN433n3pB2.0490CATALASE-PEROXIDASE443nfhA2.0492DNA-DIRECTED RNA POLYMERASE I SUBUNIT RPA49453qfkA2.0501UNCHARACTERIZED PROTEIN463n5hF2.0517FARNESYL PYROPHOSPHATE SYNTHASE473n3oB2.0529CATALASE-PEROXIDASE483o78B2.0548CHIMERA PROTEIN OF PEPTIDE OF MYOSIN LIGHT CHAIN  SMOOTH MUSCLE, GREEN FLUORESCENT PROTEIN, GREEN FLUORESCENT CALMODULIN493luyA2.0570PROBABLE CHORISMATE MUTASE50(b)
Top 50 hits using DRC for a query protein structure from the COX-2 family onDataset-B. Annotations in bold correspond to members of the COX-2 family, predicted proteins, or uncharacterized proteins
ProteinChainDistanceProtein name annotationRank2zxwB0.0000CYTOCHROME C OXIDASE SUBUNIT 212eilB1.3914CYTOCHROME C OXIDASE SUBUNIT 222eijB1.5853CYTOCHROME C OXIDASE SUBUNIT 233ag4B1.6147CYTOCHROME C OXIDASE SUBUNIT 242dysB1.6991CYTOCHROME C OXIDASE SUBUNIT 253ag1B1.8159CYTOCHROME C OXIDASE SUBUNIT 262eimB1.8824CYTOCHROME C OXIDASE SUBUNIT 273ag2B2.0173CYTOCHROME C OXIDASE SUBUNIT 282occB2.0631CYTOCHROME C OXIDASE93ablB2.1404CYTOCHROME C OXIDASE SUBUNIT 2102eikB2.1413CYTOCHROME C OXIDASE SUBUNIT 2113abmB2.1454CYTOCHROME C OXIDASE SUBUNIT 2121v55B2.1489CYTOCHROME C OXIDASE POLYPEPTIDE II132dyrB2.2044CYTOCHROME C OXIDASE SUBUNIT 2142einB2.2628CYTOCHROME C OXIDASE SUBUNIT 2151v54B2.2976CYTOCHROME C OXIDASE POLYPEPTIDE II163n56B2.4226INSULIN-DEGRADING ENZYME173abkB2.4502CYTOCHROME C OXIDASE SUBUNIT 2183p42A2.4785PREDICTED PROTEIN193r2uB2.4846METALLO-BETA-LACTAMASE FAMILY PROTEIN203msuB2.4920CITRATE SYNTHASE213ag3B2.4950CYTOCHROME C OXIDASE SUBUNIT 2223ntdB2.5052FAD-DEPENDENT PYRIDINE NUCLEOTIDE-DISULPHIDE OXIDOREDUCTASE233ngiA2.5055DNA POLYMERASE247ximB2.5198D-XYLOSE ISOMERASE253mjyA2.5206DIHYDROOROTATE DEHYDROGENASE263nvaB2.5391CTP SYNTHASE273lm3A2.5433UNCHARACTERIZED PROTEIN283ppnB2.5511GLYCINE BETAINE/CARNITINE/CHOLINE-BINDING PROTEIN293o98B2.5557BIFUNCTIONAL GLUTATHIONYLSPERMIDINE SYNTHETASE/AM303pomB2.5558RETINOBLASTOMA-ASSOCIATED PROTEIN313nt6B2.5656FAD-DEPENDENT PYRIDINE NUCLEOTIDE-DISULPHIDE OXIDOREDUCTASE325lymB2.5690LYSOZYME333n1yB2.5728TOLUENE O-XYLENE MONOOXYGENASE COMPONENT341occB2.5772CYTOCHROME C OXIDASE353lxtD2.5860GLUTATHIONE S TRANSFERASE362q70B2.5922ESTROGEN RECEPTOR373l49B2.5930ABC SUGAR (RIBOSE) TRANSPORTER, PERIPLASMIC SUBSTRATE-BINDING SUBUNIT383pvqA2.5944DIPEPTIDYL-PEPTIDASE VI393pufB2.6032RIBONUCLEASE H2 SUBUNIT B403mveB2.6077UPF0255 PROTEIN VV1_0328413ld2B2.6143PUTATIVE ACETYLTRANSFERASE423ne6A2.6160DNA POLYMERASE433qaeA2.61793-HYDROXY-3-METHYLGLUTARYL-COENZYME A REDUCTASE443qh8A2.6197BETA-LACTAMASE-LIKE453m3rA2.6215ALPHA-HEMOLYSIN463nrbB2.6234FORMYLTETRAHYDROFOLATE DEFORMYLASE473n05B2.6240NH(3)-DEPENDENT NAD(+) SYNTHETASE483m2lA2.6324ALPHA-HEMOLYSIN493pnsB2.6357URIDINE PHOSPHORYLASE50


### 4.5. Comparison with Related Methods

The use of distance distributions for protein surface analysis was studied by Binkowski et al. [11]. As earlier discussed, they did not consider the specific residues in constructing the distributions. Their distance distribution (labeled as DD1 in this work) is obtained by removing the reference to the specific residue at the center of the patch (see Figure 2).  Our use of surface residue cooccurrences and combining these with the residue-specific distance distributions are novel methods introduced in this paper. Tables 3(a) and 3(b) compare the overall classification performance using DD1 with those obtained with the proposed descriptors.Table 3(a)
Overall classification rate using different classifiers (300 training samples, 100 testing samples from DATASET-A)
ClassifierDescriptorDD1DD2RCDRCNaïve bayes58%86%94%91%Logistic58%85%99%97%Simple logistic58%89%98%91%(b)
Overall classification rate using different classifiers (100 training samples, 300 testing samples from DATASET-A)
ClassifierDescriptorDD1DD2RCDRCNaïve bayes55%74%94%93%Logistic62%88%89%94%Simple logistic63%85%91%90%Figure6 also shows the comparative performance using both the enrichment plots, and precision and recall. We define precision and recall at a given distance threshold as follows: precision = (number of correct retrievals at the threshold)/(number of total retrievals at the threshold). Recall = (number of correct retrievals at the threshold)/(number of total true matches expected at the threshold).  Here, using the ranked results, for a given query and a given rank, the number of expected true matches will be min{rank, query group size}. This is similar to the definition used in [28]. We performed queries on DATASET-A using query proteins from each of the three groups and computed the average precision and recall for each descriptor. We then computed the area under the curve (AUC) for the average precision-recall plots. The results were as follows:  DD1 (0.501052), DD2 (0.649412), RC (0.668303), and DRC (0.66759). Although the databases used are different, these results compare well with the results reported by Sael and Kihara [28], where they evaluated the retrieval performance of four surface characterization methods, based on the Zernike representation. The maximum AUC reported using standard resolution surfaces was 0.608 (without length filtering) and 0.628 (with length filtering). Yin et al. [8] proposed a fingerprint-based method, using surface alignment on selected surface patches. Their method constructs an initial patch on every vertex on the dot surface, and requires computation of geodesic distances on the surface, two very time-consuming processes. Our method neither requires surface alignment, nor expensive computations on the surface, beyond the surface generation process. Patches are generated only on positions of the surface residues, rather than over all the vertices on the generated protein surface.


### 4.6. Computation Time

The most time consuming part was for preprocessing, as needed to construct the protein surfaces and extract the protein chains. The construction of the protein surface from the original PDB files required about 4.065 seconds, running onCygwin (a version of Linux for Microsoft Windows). Extraction of the protein chains and the Cα atoms was performed using Matlab Bioinformatics Toolbox (Mathworks Inc., Natick, Mass, USA), and required 32 seconds per PDB file. Construction of the descriptors after the above steps took an average of 0.7 seconds per PDB file. Querying DATASET-B (15,386 chains, 6,261 unique PDB files) using the DRC descriptor required an average time of 0.28 seconds for each query PDB file.


## 4.1. Datasets and Environment

We performed experiments to test the performance of the proposed protein surface descriptor in two protein structure analysis tasks, namely, classifying proteins into their most likely functional groups, and ranking and retrieval of protein surfaces. We used two datasets for the experiments. DATASET-A contained information from three protein families:uracil-DNA glycosylase, cell division protein kinase 2, and estrogen receptor. This was created by scanning the PDB and selecting the protein structures with protein chains belonging to one of the three families. We were able to extract 416 chains that belong to 243 proteins in the PDB. The dataset is distributed as follows: 91 chains from 46 distinct proteins for uracil-DNA glycosylase (Group1), 186 chains from 95 distinct proteins for estrogen receptor family (Group2), and 139 chains from 102 distinct proteins from cell division protein kinase 2 (Group3). We used DATASET-A basically to train the system, and perform initial testing.  DATASET-B contained protein structures from two families, namely cyclooxygenase-2  (COX-2) (51 proteins, 95 chains) and epidermal growth factor  (EGF) (67 proteins, 71 chains). We then extracted protein structures from the PDB that have 10 or less chains and ignored the rest. This resulted in a total of 15,386 protein chains form 6,261 unique proteins. DATASET-B included all structures in DATASET-A. We used DATASET-B for a more comprehensive scan of the PDB, in the quest for potentially novel structures that may be related to the two families. Experiments were performed using a SONY VAIO personal computer, with Intel Core 2 Duo Processor T8100, running at 2.10 GHz, with 2 GB of main memory. Programs were written using Matlab (Mathworks Inc, Natick, Mass, USA) with the Bioinformatics Toolbox. We set probe radius = 1.4 Å and patch distance threshold τp=10 Å. For distance distributions, we used a fixed number of bins, b=5. Classification was performed based on algorithms implemented in Weka [38, 39] version 3-6-4.


## 4.2. Classification Performance

We divide DATASET-A into training and testing sets and apply different classifiers on the different descriptors proposed. In all our experiments, the training sets were kept very separate from the testing sets, with no overlap between the two. Classification performance is measured in terms of classification rate based on the three protein families in the dataset. We tested the method using various classifiers implemented in Weka, such as Naïve Bayes, logistic regression, and simple logistic classifier. We report results mainly for the logistic regression.  First, we explore the impact of the size of the testing set and of the training set on the classification performance using the proposed approach. We varied the size of the training set (from 50 to 300), while keeping the size of the testing set fixed. We then checked the performance using fixed testing sets of size 100, 200, and 300. Figure3 shows the results.Variation of classification rate (CR) with size of training set using the proposed descriptors DD2 (a), RC (b), and DRC (c).  Results are shown for the average over 10 runs, using logistic regression as the classifier.
(a)(b)(c)The figure shows that applying the distance distribution (DD2) alone resulted in the lowest performance accuracy as compared to using the residue cooccurrences (RCs) or distance-residue cooccurrences (DRCs). Yet, our definition of the distance distribution shows encouraging results. A steady improvement in performance with increasing training set size can be observed when using DD2 alone, peaking at about 87%  with a training size of 200 and testing size of 100. The distinctiveness of our approach is the use of residue cooccurrences on the protein surface. This approach assumes that functionally similar surface proteins have similar residue cooccurrences within a small local surface region. Figure3 (middle plot) shows that classification using residue cooccurrences (RCs) provided a significant improvement in the classification rate.  A similar improvement was observed using other classifiers, such as Naïve Bayes. Using the RC descriptor, we can achieve an accuracy rate of 94% using a small training set (50 samples) and six times larger testing set (300 chains). This shows the robustness of the residue cooccurrences, even when using a few training samples. We observe that the performance using DD2 was not as robust (about 81% using small training set, peaking at about 87% using 200 training samples).The use of distance-residue co-occurrence presents a steadier improvement in the classification rate. Using the DRC raised the accuracy rate to 99% using the simple logistic classifier on a training set of 150 and testing set of 100 (data not shown).  We can observe the significant difference between the results of DD2 (which did not use information on residue cooccurrences) and RC and/or DRC (both of which used residue cooccurrences). Figure4 shows a corresponding performance measurement with varying size of the testing set, while keeping the training set size fixed. As expected, there is a general slight decrease in performance with increasing size of the test set. The case of DRC using a training set size of 100 seemed to increase slightly with increasing testing set size. The increase is however within a small range (from 0.91 to 0.93). This shows a steady performance over increasing size of the testing set. Overall trends are similar to Figure 3, with RC and DRC performing much better than DD2.  Similar trends were also observed using other classification algorithms.  The overall classification performance is summarized in Figure 5, which shows the results of the three proposed schemes using n-fold cross validation, for different values of n.Variation of classification rate (CR) with size of testing set using the proposed descriptors DD2 (a), RC (b), and DRC (c).  Results are shown for the average over 10 runs, using logistic regression as the classifier.
(a)(b)(c)Figure 5
Summary classification performance usingn-fold cross validation (the x-axis is for varying n).


## 4.3. Ranking and Retrieval

In this section, we explore the effectiveness of our approach on the problem of search and retrieval of protein surfaces. Given a query protein, we study whether our approach has the robustness to place most of the functionally similar proteins in the top hits of the retrieved surfaces.  Here, a query protein from each of the three groups is used to screen the entire DATASET-A (416 samples) and provide a ranking based on the similarity. Thus, each protein structure is ranked against the query, (from 1 to 416), where a lower rank (smaller distance) implies more similarity to the query. After that, we search over the retrieved proteins to find which ranks the functionally similar proteins (i.e., proteins in the same functional group) have attained. Table1 shows the ranking produced using the proposed descriptor, for three query samples, one for each group. Results are shown only for DRC. RC produced a slightly better ranking (especially for uracil-DNA glycosylase family (Group 1)), while DD2 was worse than both RC and DRC). Overall, for Group 2 and Group 3, the Top 30 ranked proteins belonged to the corresponding family, while Group 1 was more difficult.Table 1
Ranking the screened proteins according to their similarity to the query protein. Results are shown for the Top 30 hits for one query protein from each of the three groups, using DATASET-A (416 protein chains).
(a)
DRC on query protein 1UDI chain I (Group 1)
Protein PDB IDChainRankDistance1UDII102ZHXB22.13061LQMB32.35091LQGC42.45892UUGC52.53531EUIC72.59202ZHXL82.61041UGHI102.63492UGIA152.68091UGIE162.68722ZHXH192.69692ZHXD212.70062ZHXN222.70172ZHXJ232.71411UGIG252.72611EMJA422.77582BOOA452.78081UGID472.78682OWRB502.79522J8XD612.81291LQGD702.82631Q3FA902.85332UUGD992.86751UGIA1012.86892OWRC1102.87492ZHXF1162.8852OWQB1292.89771SSPE1412.91151UGIC1422.91172ZHXA1472.915(b)
DRC on query protein 1QKN chain A (Group 2)
Protein PDB IDChainRankDistance1QKNA102J7XA21.37692J7YA31.57531QKMA41.67931NDEA51.73682GIUA61.73711L2IA71.74601U3RB81.76703ERDA91.76833OS9A101.77152IOGA111.77383LTXC121.77421YIMA131.78543ERTA141.79661U3QD151.80091YY4A161.80902OUZA171.81261YINA181.81521XP6A191.82602AYRA201.82693OS8D211.83112QH6A221.83123OSAA231.83671L2JA241.83852JJ3A251.84381G50A261.84903OS8A271.85092FSZA281.85182QGWA291.86041UOMA301.8683(c)
DRC on query protein 1YKR chain A (Group 3)
Protein PDB IDChainRankDistance1YKRA102UZOA21.03962R3OA31.08103PXYA41.08463PY1A51.09402WMAA61.13892IW6A71.16093NS9A81.20433IGGA91.21412WFYA101.21792C5YA111.22703DDPA121.22802J9MA131.22842R3JA141.23742R3LA151.24023PXRA161.24222DUVA171.25341W8CA181.25863DOGA191.27932V22A201.28222R3PA211.28832V22C221.29633IG7A231.32072JGZA241.32752R64A251.33032WHBA261.33812VTNA271.34583LFNA281.34762WIPA291.35142BKZA301.3580We further measured the performance of our approach using the enrichment plot. The enrichment plot essentially measures how well a given ranking or retrieval system performs, when compared with a random selection of the data samples. At a given percentage of database screening, the enrichment factor is computed as the ratioNobs/Nexp, where Nobs=number of functionally similar proteins observed or retrieved by the system, and Nexp=number of functionally similar proteins expected by random selection. For an effective system, we expect that most of the functionally similar proteins should be observed after a small percentage of screening. That is, the top hits should contain mainly functionally similar proteins, and hence the enrichment factor should be high after a small percentage probe of the database, and gradually decrease towards 1 (which corresponds to random selection). Figure 6(a) shows a plot of the average enrichment factor using 5 queries from Group 3. The enrichment plot shows that our proposed method provides better results as we screen a small percentage of the dataset. In most of the cases, our method retrieved about three times better than the expected random retrieval in the first 10% of screened proteins. As we increase the percent of screening, the retrieval degrades, since we are more likely to have retrieved most, if not all of the similar proteins after a small percentage of the screening. Thus, subsequent retrievals will lead to spurious results.Ranking and retrieval performance for the proposed methods. (a) Enrichment plot for screening protein structures using the proposed descriptors. Results are average for 5 query proteins fromcell division protein kinase 2 family (Group 3), using DATASET-A (416 protein chains).  (b) Average precision and recall for three queries, one for each group in DATASET-A. DD1 corresponds to the distance distribution proposed in [11], as described in Section 2 (see Section 4.5 on comparison with related methods).
(a)(b)


## 4.4. Screening Protein Surfaces in PDB

Encouraged by the results in classification and ranking using the proposed descriptors, we now performed a larger scale experiment, by screening the entire protein structures in PDB, using the protein chains in DATASET-B,with members of the COX-2 and EGF families as the query. The main objective was to see how the proposed descriptors will perform on a large scale, and to see if the methods could predict potentially novel functional linkages between any of the families and other proteins in PDB. For this task, we used only PDB files with 10 or less chains, and ignored the rest. This resulted in a total of 15,386 protein chains from 6,261 unique proteins. Table 2(a) shows the ranking results produced by screening the PDB files based on the proposed descriptors, using a member of the EGF family as a query. Table 2(b) shows corresponding results using a member of COX-2 family. Results are shown only for the DRC descriptor. Generally, similar results were obtained using RC. We can notice that some of the unknown proteins (annotated as “uncharacterized”) were placed in the Top-50 rankings, implying a possible relationship with the respective families.Table 2(a)
Top 50 hits using DRC for a query protein structure from the EGF family on DATASET-B. Annotations in bold correspond to members of the EGF family, predicted proteins, or uncharacterized proteins
ProteinChainDistanceProtein name annotationRank2a2qL0.0000COAGULATION FACTOR VII12firL1.4925COAGULATION FACTOR VII LIGHT CHAIN22zp0L1.5545FACTOR VII LIGHT CHAIN31wtgL1.5628COAGULATION FACTOR VII41wunL1.5844COAGULATION FACTOR VII52b8oL1.6305COAGULATION FACTOR VII LIGHT CHAIN62zwlL1.6379FACTOR VII LIGHT CHAIN72zzuL1.6536FACTOR VII LIGHT CHAIN81wqvL1.6816COAGULATION FACTOR VII92ec9L1.6832COAGULATION FACTOR VII101danL1.7655BLOOD COAGULATION FACTOR VIIA111wssL1.7659COAGULATION FACTOR VII122puqL1.7692COAGULATION FACTOR VII131fakL1.7934PROTEIN (BLOOD COAGULATION FACTOR VIIA)142b7dL1.8024COAGULATION FACTOR VII156acnA1.8061ACONITASE162aerL1.8120COAGULATION FACTOR VII172aeiL1.8164COAGULATION FACTOR VII182flrL1.8196COAGULATION FACTOR VII193elaL1.8668COAGULATION FACTOR VII LIGHT CHAIN201z6jL1.8859COAGULATION FACTOR VII212f9bL1.9027COAGULATION FACTOR VII223phsA1.9169CELL WALL SURFACE ANCHOR FAMILY PROTEIN233n54B1.9263SPORE GERMINATION PROTEIN B3243qbpB1.9264FUMARASE FUM253ma9L1.9367TRANSMEMBRANE GLYCOPROTEIN263mt0A1.9921UNCHARACTERIZED PROTEIN PA1789273lguA2.0004PROTEASE DEGS283m7iA2.0071TRANSKETOLASE291qfkL2.0096PROTEIN (COAGULATION FACTOR VIIA (LIGHT CHAIN))303n9tA2.0169PNPC313pxzA2.0235CELL DIVISION PROTEIN KINASE 2322flbL2.0257COAGULATION FACTOR VII333lh1A2.0289PROTEASE DEGS343nlcA2.0300UNCHARACTERIZED PROTEIN VP0956353no5C2.0302UNCHARACTERIZED PROTEIN363msqC2.0347PUTATIVE UBIQUINONE BIOSYNTHESIS PROTEIN373rykA2.0389DTDP-4-DEHYDRORHSMNOSE 3,5-EPIMERASE383m4aA2.0444PUTATIVE TYPE I TOPOISOMERASE393n3nB2.0444CATALASE-PEROXIDASE402R3GA2.0456CELL DIVISION PROTEIN KINASE 2412R3IA2.0465CELL DIVISION PROTEIN KINASE 2423o0rL2.0473ANTIBODY FAB FRAGMENT LIGHT CHAIN433n3pB2.0490CATALASE-PEROXIDASE443nfhA2.0492DNA-DIRECTED RNA POLYMERASE I SUBUNIT RPA49453qfkA2.0501UNCHARACTERIZED PROTEIN463n5hF2.0517FARNESYL PYROPHOSPHATE SYNTHASE473n3oB2.0529CATALASE-PEROXIDASE483o78B2.0548CHIMERA PROTEIN OF PEPTIDE OF MYOSIN LIGHT CHAIN  SMOOTH MUSCLE, GREEN FLUORESCENT PROTEIN, GREEN FLUORESCENT CALMODULIN493luyA2.0570PROBABLE CHORISMATE MUTASE50(b)
Top 50 hits using DRC for a query protein structure from the COX-2 family onDataset-B. Annotations in bold correspond to members of the COX-2 family, predicted proteins, or uncharacterized proteins
ProteinChainDistanceProtein name annotationRank2zxwB0.0000CYTOCHROME C OXIDASE SUBUNIT 212eilB1.3914CYTOCHROME C OXIDASE SUBUNIT 222eijB1.5853CYTOCHROME C OXIDASE SUBUNIT 233ag4B1.6147CYTOCHROME C OXIDASE SUBUNIT 242dysB1.6991CYTOCHROME C OXIDASE SUBUNIT 253ag1B1.8159CYTOCHROME C OXIDASE SUBUNIT 262eimB1.8824CYTOCHROME C OXIDASE SUBUNIT 273ag2B2.0173CYTOCHROME C OXIDASE SUBUNIT 282occB2.0631CYTOCHROME C OXIDASE93ablB2.1404CYTOCHROME C OXIDASE SUBUNIT 2102eikB2.1413CYTOCHROME C OXIDASE SUBUNIT 2113abmB2.1454CYTOCHROME C OXIDASE SUBUNIT 2121v55B2.1489CYTOCHROME C OXIDASE POLYPEPTIDE II132dyrB2.2044CYTOCHROME C OXIDASE SUBUNIT 2142einB2.2628CYTOCHROME C OXIDASE SUBUNIT 2151v54B2.2976CYTOCHROME C OXIDASE POLYPEPTIDE II163n56B2.4226INSULIN-DEGRADING ENZYME173abkB2.4502CYTOCHROME C OXIDASE SUBUNIT 2183p42A2.4785PREDICTED PROTEIN193r2uB2.4846METALLO-BETA-LACTAMASE FAMILY PROTEIN203msuB2.4920CITRATE SYNTHASE213ag3B2.4950CYTOCHROME C OXIDASE SUBUNIT 2223ntdB2.5052FAD-DEPENDENT PYRIDINE NUCLEOTIDE-DISULPHIDE OXIDOREDUCTASE233ngiA2.5055DNA POLYMERASE247ximB2.5198D-XYLOSE ISOMERASE253mjyA2.5206DIHYDROOROTATE DEHYDROGENASE263nvaB2.5391CTP SYNTHASE273lm3A2.5433UNCHARACTERIZED PROTEIN283ppnB2.5511GLYCINE BETAINE/CARNITINE/CHOLINE-BINDING PROTEIN293o98B2.5557BIFUNCTIONAL GLUTATHIONYLSPERMIDINE SYNTHETASE/AM303pomB2.5558RETINOBLASTOMA-ASSOCIATED PROTEIN313nt6B2.5656FAD-DEPENDENT PYRIDINE NUCLEOTIDE-DISULPHIDE OXIDOREDUCTASE325lymB2.5690LYSOZYME333n1yB2.5728TOLUENE O-XYLENE MONOOXYGENASE COMPONENT341occB2.5772CYTOCHROME C OXIDASE353lxtD2.5860GLUTATHIONE S TRANSFERASE362q70B2.5922ESTROGEN RECEPTOR373l49B2.5930ABC SUGAR (RIBOSE) TRANSPORTER, PERIPLASMIC SUBSTRATE-BINDING SUBUNIT383pvqA2.5944DIPEPTIDYL-PEPTIDASE VI393pufB2.6032RIBONUCLEASE H2 SUBUNIT B403mveB2.6077UPF0255 PROTEIN VV1_0328413ld2B2.6143PUTATIVE ACETYLTRANSFERASE423ne6A2.6160DNA POLYMERASE433qaeA2.61793-HYDROXY-3-METHYLGLUTARYL-COENZYME A REDUCTASE443qh8A2.6197BETA-LACTAMASE-LIKE453m3rA2.6215ALPHA-HEMOLYSIN463nrbB2.6234FORMYLTETRAHYDROFOLATE DEFORMYLASE473n05B2.6240NH(3)-DEPENDENT NAD(+) SYNTHETASE483m2lA2.6324ALPHA-HEMOLYSIN493pnsB2.6357URIDINE PHOSPHORYLASE50


## 4.5. Comparison with Related Methods

The use of distance distributions for protein surface analysis was studied by Binkowski et al. [11]. As earlier discussed, they did not consider the specific residues in constructing the distributions. Their distance distribution (labeled as DD1 in this work) is obtained by removing the reference to the specific residue at the center of the patch (see Figure 2).  Our use of surface residue cooccurrences and combining these with the residue-specific distance distributions are novel methods introduced in this paper. Tables 3(a) and 3(b) compare the overall classification performance using DD1 with those obtained with the proposed descriptors.Table 3(a)
Overall classification rate using different classifiers (300 training samples, 100 testing samples from DATASET-A)
ClassifierDescriptorDD1DD2RCDRCNaïve bayes58%86%94%91%Logistic58%85%99%97%Simple logistic58%89%98%91%(b)
Overall classification rate using different classifiers (100 training samples, 300 testing samples from DATASET-A)
ClassifierDescriptorDD1DD2RCDRCNaïve bayes55%74%94%93%Logistic62%88%89%94%Simple logistic63%85%91%90%Figure6 also shows the comparative performance using both the enrichment plots, and precision and recall. We define precision and recall at a given distance threshold as follows: precision = (number of correct retrievals at the threshold)/(number of total retrievals at the threshold). Recall = (number of correct retrievals at the threshold)/(number of total true matches expected at the threshold).  Here, using the ranked results, for a given query and a given rank, the number of expected true matches will be min{rank, query group size}. This is similar to the definition used in [28]. We performed queries on DATASET-A using query proteins from each of the three groups and computed the average precision and recall for each descriptor. We then computed the area under the curve (AUC) for the average precision-recall plots. The results were as follows:  DD1 (0.501052), DD2 (0.649412), RC (0.668303), and DRC (0.66759). Although the databases used are different, these results compare well with the results reported by Sael and Kihara [28], where they evaluated the retrieval performance of four surface characterization methods, based on the Zernike representation. The maximum AUC reported using standard resolution surfaces was 0.608 (without length filtering) and 0.628 (with length filtering). Yin et al. [8] proposed a fingerprint-based method, using surface alignment on selected surface patches. Their method constructs an initial patch on every vertex on the dot surface, and requires computation of geodesic distances on the surface, two very time-consuming processes. Our method neither requires surface alignment, nor expensive computations on the surface, beyond the surface generation process. Patches are generated only on positions of the surface residues, rather than over all the vertices on the generated protein surface.


## 4.6. Computation Time

The most time consuming part was for preprocessing, as needed to construct the protein surfaces and extract the protein chains. The construction of the protein surface from the original PDB files required about 4.065 seconds, running onCygwin (a version of Linux for Microsoft Windows). Extraction of the protein chains and the Cα atoms was performed using Matlab Bioinformatics Toolbox (Mathworks Inc., Natick, Mass, USA), and required 32 seconds per PDB file. Construction of the descriptors after the above steps took an average of 0.7 seconds per PDB file. Querying DATASET-B (15,386 chains, 6,261 unique PDB files) using the DRC descriptor required an average time of 0.28 seconds for each query PDB file.


## 5. Conclusion

We have introduced a novel approach to the description and characterization of protein surfaces. The proposed approach captures the surface structure of the protein by utilizing local patches defined only on the positions of surface residues, rather than over all surface vertices, or over all the surface atoms. We make residue cooccurrences on the surface a central part of the descriptor. The novelty of this approach can be observed by the ease of handling both local and global variation on the surface (using local and global descriptors). Moving from local to global not only reduces the computational problem of matching 3D structures, but also facilitates direct comparison between protein structures of different sizes. By avoiding the construction of the complete 3D surface and retaining only the surface Cα to do the analysis, the need for surface alignment of the 3D structure is eliminated. Further, we do not need to perform any geometrical transformation to insure reliable matching. This is very important for rapid analysis over a large database, such as the PDB.We showed results on the performance of the proposed methods in functional classification of proteins into their putative families, based on the surface information. We further compared the results using enrichment plots, and the standard measures of precision and recall. For the three protein families used, we obtained an area under the curve for precision and recall of 0.6494 (DD2), 0.6683 (RC), and 0.6676 (DRC). A screening of the PDB usingCOX-2 and EGF family members showed that the proposed methods ranked related family members in the Top-20 hits, with a number of uncharacterized proteins also retrieved. It will be interesting to perform further biological lab experiments to verify if any of the retrieved uncharacterized proteins are truly related to the respective families to which they share similar surfaces (as determined by our surface descriptors).



---

*Source: 918978-2011-11-22.xml*