 Methodology
 Open Access
 Published:
Automated shapebased clustering of 3D immunoglobulin protein structures in chronic lymphocytic leukemia
BMC Bioinformatics volume 19, Article number: 414 (2018)
Abstract
Background
Although the etiology of chronic lymphocytic leukemia (CLL), the most common type of adult leukemia, is still unclear, strong evidence implicates antigen involvement in disease ontogeny and evolution. Primary and 3D structure analysis has been utilised in order to discover indications of antigenic pressure. The latter has been mostly based on the 3D models of the clonotypic B cell receptor immunoglobulin (BcR IG) amino acid sequences. Therefore, their accuracy is directly dependent on the quality of the model construction algorithms and the specific methods used to compare the ensuing models. Thus far, reliable and robust methods that can group the IG 3D models based on their structural characteristics are missing.
Results
Here we propose a novel method for clustering a set of proteins based on their 3D structure focusing on 3D structures of BcR IG from a large series of patients with CLL. The method combines techniques from the areas of bioinformatics, 3D object recognition and machine learning. The clustering procedure is based on the extraction of 3D descriptors, encoding various properties of the local and global geometrical structure of the proteins. The descriptors are extracted from aligned pairs of proteins. A combination of individual 3D descriptors is also used as an additional method. The comparison of the automatically generated clusters to manual annotation by experts shows an increased accuracy when using the 3D descriptors compared to plain bioinformaticsbased comparison. The accuracy is increased even more when using the combination of 3D descriptors.
Conclusions
The experimental results verify that the use of 3D descriptors commonly used for 3D object recognition can be effectively applied to distinguishing structural differences of proteins. The proposed approach can be applied to provide hints for the existence of structural groups in a large set of unannotated BcR IG protein files in both CLL and, by logical extension, other contexts where it is relevant to characterize BcR IG structural similarity. The method does not present any limitations in application and can be extended to other types of proteins.
Background
The concept of molecular similarity underlies a methodology where molecules are grouped together based on their biological effects, physicochemical properties and threedimensional structures [1]. Considering that the threedimensional (3D) protein structure plays a pivotal role in protein functional characterization [2], the comparison of the threedimensional (3D) molecular structures is a key technique in a variety of applications such as protein function prediction, computer aided molecular design, rational drug design and protein docking [3].
In the absence of known structure, alternative approaches such as comparative modeling can provide a 3D model of a protein, related to at least one experimentally determined protein structure. The most comprehensive examples of these approaches are SCOP [4] and CATH [5], protein structure classification databases that were established to address the evolutionary relationships between protein structures. They are widely used as a benchmark for novel protein structure comparison methods and as a training dataset for machine learning algorithms focused on protein structure classification and prediction [6]. Their rationale is that protein structures are conserved during evolution and the existence of a protein family would facilitate the identification of related proteins through similarities in their structures [7].
Techniques that define similarity between 3D structures can be classified into three categories, i.e. (1) superposition of protein structures where alignment between equivalent residues in not given a priori [8], (2) feature representation of protein spatial profile in multidimensional vectors [9] and (3) time series formed by the alteration of the protein tertiary structure [10].
In the first category, the structural similarity is determined by scaling, rotation, transformation and then superpositioning [11]. Numerous scoring functions have been proposed towards the definition of the positional deviations of equivalent atoms upon rigidbody superimposition. Aligners were implemented with the ability to identify similarities between proteins with large conformational changes. Various metrics for comparing and scoring identity between two protein structures are employed but the most commonly used are pvalues [12] and root mean square deviation (RMSD) [2]. Highlighted aligners in this category are represented in Table 1. Although this type of approach is very effective, it is a computationally expensive and time consuming method.
The second approach includes all the shapebased methods. In shapebased approaches, the protein is treated as a 3D object and represented by a multidimensional vector that uniquely characterizes the object. Consequently the comparison between feature vectors is characterized by lower complexity and higher accuracy. Similarity search is committed through global or local features. The global features are computed by the transformation of euclidean space into a metric space that measures the pairwise distances between the points of the 3D objects. The global features are invariant to the deformations of the 3D object. The local features are computed on each keypoint of the surface by accumulating pairwise relations among oriented surface points into a local histogram [3].
The last method is related to the comparison between time series. According to this type of methodology, protein structures are translated into polygonal chains [13]. The aforementioned transformation of the 3D object to a feature vector reduces object complexity and it can be handled as a time series [10]. Protein tertiary structures, such as the alphacarbon atoms along the backbone of a protein, essentially form a 3D polygonal chain and a natural measure for comparing the geometric similarity estimates the similarity between the structures.
The current work aims to categorize chronic lymphocytic leukemia (CLL) patients based on their 3dimensional protein structures of the clonotypic B cell receptor immunoglobulin (BcR IG) amino acid (IG) sequences following on the shapebased approach. In the paragraphs that follow, the corresponding stateoftheart analysis is presented.
Geometrical descriptor vectors able to achieve very fast comparisons especially for applications of virtual screening are described in [14, 15]. Spin Images [16] and Shape Histograms [17] are methodologies in molecular surface representation [18]. The former is related to local 2D descriptors while the later exploits the global geometric properties of the molecules. Efforts for the implementation of multiview methods [19, 20] of molecular surface representation were proved insufficient in proteins in cases without symmetries. Computational approaches that address the significance of small variations between 3D protein structure of high similarity are conducted in the methodologies of [21–23].
Pattern recognition establishes approaches that extract moments from the 3D object. Zernike moments were applied in [24] for a feature representation on a PositionSpecific Scoring Matrix (PSSM). Zernike moments descriptors were utilized to extract features in each protein PSSM forming a 42dimensional feature vector. Finally, machine learning methods called PCVM where applied to accomplish classification. A similar method that implements Legendre moments to predict proteinprotein interactions is described in [25]. Zernike descriptors provide a rotation invariant ability to the protein shape comparison as they do not necessitate structural alignment. Additionally, they allow other characteristics of a protein surface, such as electrostatic potentials, to be incorporated into the descriptor vector [26]. Zernike moments can be applied in several problems related to protein structures with satisfactory results. Regarding moment extraction, 2D PolarFourier coefficients [27] and 2D Krawtchouk moments [28] create a rotation invariant feature vector by taking as input the volume of the 3D object. Spherical Harmonics [29] are widely used for a large scale of structural similarity comparison. Although the formation of an orthonormal 1D vector allows fast comparisons, they present an inaccuracy in the results that is connected to the alignment parameters.
Rigid object methodologies are inherently limited by ignoring the flexibility of the molecule. To overcome this limitation, approaches that respect the shape deformation of molecules were utilized. The two main categories on feature extraction of nonrigid approaches are: a) globalshapebased [30, 31] and b) localshapebased methods [32–34]. The former usually is used to create a metric where Euclidean space or Euclidean metrics are transformed into pairwise distances between points of the 3D object surface. The aforementioned points are invariant to deformations of the 3D object. The final descriptor vector is formed by the feature histograms of the distances. Localshapebased methods sample the surface and extract descriptors for each of the sampled local regions. When the local descriptors are extracted then a featurebased methodology is implemented in order to translate them into global. Besides the discriminative ability between the proper local shape descriptors, they also satisfy significant criteria such as fast descriptor extraction, compactness and rotation invariance.
The motivation of this study was to examine the abilities of global and local 3D descriptors for pairwise distance calculation, instead of applying bioinformaticsspecific similarity scores. We hypothesized that their high capability in describing general 3D structures could be applied to the comparison of 3D protein structures. The structures for analysis emerged from the primary sequences of the clonotypic BcR IG of patients with CLL.
CLL is the most common adult leukemia, with still unclear etiology. That said, primary and 3D structurebased reasoning strongly implicates antigen selection in disease ontogeny and evolution [35, 36]. Molecular categorization of CLL patients based on BcR IG sequence similarity has so far been addressed using bioinformatics methods of structural similarity calculation [35, 37]. The novelty of the approach proposed herein lies in the combination of current stateoftheart bioinformatics methods with the extraction of features arising from 3D object recognition methods. The most uptodate 3D prediction structure algorithms were implemented to construct patients’ models. The proposed combined methodology achieves an efficient grouping of CLL patients in accordance to their biological and clinical features, especially in light of the recently identified stereotyped subsetsb [37, 38]. The results confirm the original hypothesis that the combination of bioinformaticsspecific techniques, such as TMalign, and generalpurpose 3D descriptors achieves a high discriminative power compared to using only bioinformaticsspecific methods.
Methods
In Fig. 1, a pipeline diagram presents the layers that compose the proposed methodology. The method is separated in three main levels, i.e (1) primary data collection, (2) creation of the 3D protein structures and (3) 3D protein structure comparison method. A baseline method raised from pure bioinformatics approaches was utilized as a benchmark for validation.
Extensively, the baseline benchmark consists of the following steps. First, the original BcR IG heavy and light amino acid sequences of 925 CLL patients were transformed to 3D protein models, using stateoftheart protein structure prediction tools. This resulted in a set of protein models in PDB format. For each pair of models, structural alignment was carried out and a structural similarity score was computed. The TMalign and TMscore algorithms were used for the alignment and the similarity score computation, respectively. A similarity matrix was formed to consider the similarity scores between every pair of models. The matrix was was ultimately used to organize the proteins into clusters, using various existing clustering methods.
The proposed method modified this baseline approach by replacing the bioinformatics related TMscore to a similarity metric between a pair of models with 3D descriptors originated from the field of 3D object recognition. Specifically, the Fast Point Feature Histograms (FPFH) [39], 3D Shape Context (3DSC) [40], Radiusbased Surface Descriptors (RSD) [41], Viewpoint Feature Histograms (VFH) [42] and a combination of the above local descriptors where applied to the 3D structures in order to extract the appropriate features for the comparison. The descriptors were compared using the root mean square deviation (RMSD) distance metric. The distance for every pair of models was computed and all results were incorporated into a distance matrix, which was finally used to model clustering.
Structural alignment
The TMalign algorithm [43] identifies the best structural alignment between protein pairs by combining the TMscore rotation matrix and Dynamic Programming (DP). When comparing two protein structures, the second model is rotated and translated appropriately, until the maximum alignment between the two structures is achieved [44]. The Kabsch algorithm, a method for calculating the optimal rotation matrix that minimizes the RMSD between two paired sets of points performed the structure comparison.
The alignment process includes alignment over the secondary structures of the BcR IG sequences, based on the gapless matching of the IG sequences, and alignment using an equally weighted combination of the previously extracted results. Heuristic iterations were applied, meaning that the steps were repeated until the alignment became stable and the highest TMscore was achieved. The TMscore is calculated as:
where, L_{Target} is the length of the target protein that the source protein is aligned to, L_{align} is the number of aligned parts, d_{i} is the distance between the i_{th} pair of residues, and d_{0} (L_{Target}) is given by the equation:
3D descriptors used
The proposed method is based on the extraction of 3D descriptors from the raw PDB files. Considering a pair of PDB files, they are first aligned in space and then 3D descriptors are extracted from them. The descriptors capture specific geometric properties of the distribution of the atom positions. The following four types of descriptors have been considered:

Fast Point Feature Histogram (FPFH  local descriptor)

3D Shape Context (3DSC  local descriptor)

Radiusbased Surface Descriptor (RSD  local descriptor)

Viewpoint Feature Histogram (VFH  global descriptor)
The above descriptors are briefly described in the following sections. In order to apply the descriptors, the proteins are considered as 3D point clouds, where each point of the cloud corresponds to a protein atom. For the conducted experiments presented in this paper, the implementations contained in the PCL C++ library [45] were used for these descriptors.
Fast point feature histogram (FPFH)
The Fast Point Feature Histogram (FPFH) [39] represents the angular variations in the neighborhood of each point in the input point set. It is a modification of the Point Feature Histogram (PFH) descriptor [46], towards faster computation and smaller size. The PFH descriptor considers the kneighborhood around each query point i.e. atom. For each pair of points in the neighborhood, a coordinate system is formed by the pair of points and the surface normals at these points, and 3 measures of angular variation between the points and the normals are computed. The histogram of these measures computed for all pairs of points in the neighborhood constitutes the PFH descriptor for the query point.
In the FPFH descriptor, instead of considering all pairs of points within the neighborhood of the query point, only the relationships between the query point and its neighbors are initially considered. The resulting histograms constitute an intermediate Simplified Point Feature Histogram (SPFH) descriptor for each point. The final FPFH descriptor is computed by adding the SPFH descriptors of the neighbors of the query point, weighted by their inverse distance to the query point:
where p is the query point, k is the number of its nearest neighbors, p_{i} is the ith nearest neighbor point, and d_{i} is the distance between p and p_{i}.
The resulting FPFH descriptor for a query point is a 33dimensional vector, being more compact than the original PFH descriptor, as well as faster to extract. The FPFH descriptor is invariant to position, scale and orientation, as the original PFH descriptor. However, due to numerical limitations, FPFH may not be completely invariant to rotation in certain cases. For this reason, the protein files are aligned prior to computing the FPFH descriptors, in order to ensure that there is no rotational variance.
3D shape context (3DSC)
The 3D Shape Context (3DSC) [40] descriptor describes the distribution of points around each point of a point cloud. Considering each point p of the point cloud, the support region of p is defined as a sphere centered at p. The orientation of the sphere is determined by the surface normal at p. The support region is divided into bins, determined by uniform divisions in the two angular dimensions (azimuth and elevation) and logarithmic divisions along the radial dimension.
The 3D division of the support region defines the bins of a histogram of point counts within this region. The value of the (i,j,k) bin, where i, j and k are indexes of the angular and radial position of the bin, is computed by counting the points of the point cloud that fall within the bin. However, the contribution of each point p_{i} within the support region is weighted by the following factor:
where V(i,j,k) is the volume of bin (i,j,k) and ρ_{i} is the density of the points around p_{i}. The density is computed by counting the points within a sphere of radius δ around p_{i}.
Radiusbased surface descriptor (RSD)
The Radiusbased Surface Descriptor (RSD) [41] describes the region around a point by approximating the local surfaces with spheres and estimating the minimum and maximum radii of the fitted spheres. Considering a point p and a point of its neighborhood p_{i}, the two points can be thought as lying on a sphere of radius r. Since infinite spheres pass through the two points, the one that also respects the surface normals at the two points is selected. If the distance between the two points is d and the angle between their surface normals is α, then the radius of the fitted sphere can be computed using the formula for determining the length of a sphere chord:
For computational efficiency reasons, the radius is computed from d and α using just the first terms of the Taylor expansion of the above equation. After computing the radii of the spheres for every neighbor of point p, the minimum and maximum ones are kept and used as the descriptor for p, thus obtaining a very compact yet discriminating descriptor.
Viewpoint feature histogram (VFH)
The three above descriptors are local descriptors, i.e. they compute a vector representation for each point in the point cloud of the 3D object considered. The Viewpoint Feature Histogram (VFH) [42] is a global descriptor, describing the whole point cloud. The VFH descriptor is conceptually based on the FPFH descriptor. Instead of computing FPFH descriptors for each point of the point cloud and its neighborhood, a single FPFH descriptor is computed for the object’s centroid, considering all the points of the point cloud as its neighbors.
This central descriptor constitutes one part of the VFH descriptor. The other part considers the histogram of the angles between the normals at each point of the point cloud and a fixed direction, determined by a fixed viewpoint, outside the point cloud. First, the vector from the viewpoint to the object’s centroid is computed and then the angles between this vector and each of the normals of the point cloud are used to construct an angle histogram. This histogram constitutes the second part of the VFH descriptor. Considering a fixed direction from the viewpoint to the centroid, instead of considering the direction from the viewpoint to each point in the point cloud ensures scale invariance. However, the VFH descriptor is not rotation invariant, since it depends on the object’s pose relative to the viewpoint. However, the advantage of the VFH descriptor is the compactness offered by a global descriptor, as it represents the whole object with a single vector.
Distance matrix calculation
After extracting the 3D descriptors from a pair of aligned proteins, a distance measure can be calculated between them, quantifying their structural differences in the corresponding descriptor space. In the general case, a descriptor extracted from protein i, is a set of vectors F_{i}={f_{i,1},f_{i,2},…,f_{i,L}}, where L is the number of points in the protein’s 3D model. The feature vector for point k of protein i is a vector \(\mathbf {f}_{i,k} \in \mathbb {R}^{D}\), where D is the descriptor dimensionality, which is generally different for different descriptor types. The above formulation fits well with the local 3D descriptors, such as FPFH, 3DSC and RSD, since they consist of a feature vector for each point in the 3D point cloud. However, the same formulation can be used for global descriptors, such as VFH, as well, if the global descriptor is considered as a local descriptor extracted from only a single point. Thus, in the following, the same formulation is used for all types of descriptors.
In order to compare between the descriptors of two proteins, a distance measure that can assess the difference in space between two sets of points needs to be used. The Root Mean Square Deviation (RMSD) distance metric has been used hereby for this purpose. The RMSD metric is commonly used for the comparison between protein structures, by computing an average of the pointtopoint differences among the protein atoms. However, hereby it is not used to compare the actual 3D coordinates of the atoms, but instead the highdimensional coordinates of the extracted feature vectors for each pair of points. The RMSD distance metric between proteins i and j is defined as follows:
where · denotes the Euclidean distance. In case that the descriptor type is a global one, such as VFH, i.e. L=1, the RMSD metric is reduced to the Euclidean distance between the descriptors vectors. The smaller the RMSD metric, the closer the corresponding proteins are, in terms of their similarity with respect to the corresponding descriptor type.
The RMSD distance is computed between every pair of proteins in the considered dataset, so that a square symmetric distance matrix is computed. This distance matrix can then be provided as input to clustering algorithms.
Clustering methods
In order to cluster the proteins in groups of similar characteristics with respect to the various types of descriptors considered, the following clustering methods have been used:

kmedoids

hierarchical agglomerative

DBScan
The kmedoids method is similar to the kmeans clustering method, with the difference that it uses a distance matrix as input instead of vectorial representations of the objects to cluster. The kmeans method proceeds by guessing at the cluster center positions within the feature space, and gradually updating them, in order to fit better with the available data. In kmedoids, no feature space is defined, so there is no notion of cluster center positions. Instead, only a distance matrix is provided as input, encoding the similarities and differences among objects. In such a case, the role of cluster centers is played by representative objects from the set of objects themselves, called medoids. At the beginning, the set of medoids is selected arbitrarily, e.g. randomly, from the set of objects in the dataset. The objects are grouped by assigning each object to the closest one, in terms of the distance matrix used. Then, through an iterative procedure, the medoids are updated, selecting more representative objects, so that they are better fitted with the other objects. Considering each group of objects, the object with the smallest sum of distances from the other objects in the group is selected as the group medoid.
The hierarchical agglomerative clustering method is a bottomup clustering approach, building gradually larger clusters of data, in a hierarchical manner. At the beginning, each data point is considered as a separate cluster. At each iteration, the two clusters that are nearest to each other are merged into a single cluster. The procedure continues until all clusters are merged into a single cluster, containing all data points. The distance between two clusters is hereby defined as the mean value of the pairwise distances between each pair of points in the two clusters. Hierarchical clustering results in a treelike representation of the data. Cutting the tree at a specified height results in clustering the data at different granularities. In this paper, the height of the cutting is determined by specifying the number of desired clusters.
DBScan is a densitybased clustering method, which does not require the number of clusters to be known from the beginning. It starts from an arbitrarily chosen point and considers its ε neighborhood. If it is in a dense part, it forms a cluster, also cotaining the ε neighborhoods of its neighbor points. The procedure is repeated for each point in the cluster, until no other point can be considered as being near the points of the cluster. Then, another point is selected, to begin a new cluster. DBScan uses two parameters: the neighborhood size ε and the minimum number of points in the neighborhood, in order to characterize the neighborhood as dense.
All the above clustering method do not require the objects to be represented by vectors; they only need distances to be defined between objects. This fits well with the local 3D descriptors used hereby, since a protein is not represented by a single feature vector, but rather by a set of feature vectors, one for each point in the protein model. However, the representation as sets of vectors does allow the definition of distances among proteins, e.g. using the RMSD measure, as described above, which makes the above methods suitable. Another approach would be to use the distance matrix as the input to methods such as multidimensional scaling, in order to map the proteins to points in a lowdimensional space, before performing traditional clustering methods such as kmeans. However, this could potentially lead to information loss, if the selected space dimensionality does not correspond to the underlying intrinsic dimensionality of the points. Using the distance matrix directly as input overcomes this issue.
The kmedoids and the hierarchical agglomerative methods require the number of clusters to be a priori provided as an input parameter. However, in the exploratory task of examining a set of proteins for clusters, the number of clusters to be discovered is not known. A method that can also determine the number of clusters in the data is needed. In this paper, this issue is overcome by performing the clustering several times, considering a range of number of clusters and selecting the one that maximizes a certain clustering quality criterion. The average silhouette width has been used hereby as this clustering quality criterion. Considering a protein i, let α_{i} be the average distance of protein i to all other proteins of the same cluster. Let also b_{i} be the minimum of the average distance of protein i to the proteins of all other clusters. The silhouette width for protein i is defined as:
The silhouette width takes values in the range from 1 to 1. Values close to 1 mean that b_{i} is large and α_{i} is small, which means that object i is very close to the other objects of its cluster, while at the same time it is far away from the objects of the other clusters. This in turn means that object i has been correctly clustered. On the other hand, values close to 1 mean that the object would be more properly assigned to another cluster. Considering the average silhouette width for all objects provides a measure of the clustering quality. Large values of the average silhouette width mean that the clustering produced compact and clearly divisible clusters. Thus, seeking for the number of clusters that achieves the largest average silhouette width is equivalent to seeking for the number of clusters that is most appropriate to describe the underlying data.
Combination of descriptors based on their clustering performance
In addition to the individual descriptor types (FPFH, 3DSC, RSD and VFH), a combination of them has also been considered. The various descriptor types are diverse in form, since they may be local or global and also contain vectors of different dimensionalities. This makes the process of combining them not straightforward. However, since the input to the clustering algorithm is not the descriptors themselves but rather the distance matrices produced from them, a natural way to combine the multiple descriptors is by merging their corresponding distance matrices.
Let D_{m} be the distance matrix associated with descriptor type m. The combined distance matrix D is defined as the weighted sum:
where the sum is over all descriptor types considered. Two different approaches have been examined for the definition of the weights used in the sum. The first approach is to consider them all equal to 1. This creates a distance matrix where the distance for a pair of proteins is the sum of the distances computed for this pair using the various descriptor types. Taking this sum implicitly considers that all descriptor types are equal in terms of discriminating power. In reality, some descriptor types may be more suitable for clustering the protein datasets than others. In order to handle such differences, the second approach is to consider unequal weights for the sum, ones that reflect the discriminating power of the descriptors. Hereby, the average silhouette widths of the clusterings produced by the different types of descriptors have been used as the weights. The higher the average silhouette width of a clustering, the more fitting the clustering is for the underlying data, i.e. the more descriptive the corresponding descriptor may be for the protein set. Thus, using a high average silhouette width as a weight for the sum, would mean that more importance is given to the corresponding descriptor while computing the combined distance matrix.
Note that the average silhouette width of a clustering also depends on the number of clusters. The maximum average silhouette width computed for a descriptor type, using a range of cluster numbers has been used as the weight for this descriptor.
Results
For the experimental evaluation of the proposed methodology, two different protein model datasets were formed. The clusters were evaluated externally and internally in order to obtain cluster accuracy and quality respectively.
Datasets
The datasets derive from BcR IG sequences obtained from 925 CLL cases diagnosed according to the iwCLL criteria [47]. Following established bioinformatics methods, 137/925 cases were found to belong to subsets with stereotyped i.e. (quasi)identical BcR IG, hereafter designated as stereotyped subsets [48, 49]. As a first step, we examined a dataset including BcR IG sequences from cases belonging to six wellcharacterized stereotyped subsets. As a second step, we analyzed BcR IG sequences from all cases, stereotyped and nonstereotyped. The first dataset was deployed as ground truth for the evaluation of the proposed method. Hence the second dataset formed for unsupervised clustering applications. The number of CLL sequences in each dataset is described in Table 2. The subset size distribution of the dataset consisting of only the stereotyped BcR IG is summarized in Table 3.
Regarding the BcR IG structure prediction, the CPIGS method was used, which is based on the Prediction of Immunoglobulin Structure (PIGS) [50] approach. In the PIGS web server, antibody VL and VH framework regions were used as input parameters. The sequence identity within both chains was examined with the threshold of 80%. If the aforementioned criterion was not satisfied then the two templates with the highest sequence similarity, measured by the Blosum24 score [51] of both the light and heavy chain were utilized. The H3 loop was always modeled using the template with the best sequence similarity; the other loops were modeled using a different template only if the corresponding loop in the framework template did not display the same length and canonical structure of the target loop. Finally, the CPIGS models were built by remodeling the H3 loop of the PIGS models using the template identified by the approach developed in [52]. An outstanding study of the customized PIGS methodology is described by [53].
Clustering of annotated proteins
In order to evaluate the accuracy of clustering using the 3D descriptors, a first round of experiments was conducted, using the annotated dataset. Each type of 3D descriptor was used to cluster the data into 6 clusters, as many as the ground truth subsets in the annotated dataset. The resulting clusterings were compared to the ground truth clustering, i.e. the one where each cluster corresponds to the established protein subset. The same procedure was also performed using the TMscore for clustering.
The kmedoids, Agglomerative Hierarchical clustering and Densitybased spatial clustering of applications with noise (DBSCAN) methods were used to cluster the data using the extracted descriptors. After extracting descriptors from each pair of aligned proteins, the RMSD distances between each pair of descriptors is computed, forming a square distance matrix. This matrix is used as the input to the clustering methods.
In order to compare the resulting clusterings to the ground truth clustering, the Rand index was used. The Rand index measures the number of agreements between the two compared clustering, over all pairs of points. Considering a set of N objects S={o_{1},o_{2},…,o_{N}}, e.g. a set of proteins, and two partitionings of this set, X={X_{1},X_{2},…,X_{m}}, partitioning the objects into m groups, and Y={Y_{1},Y_{2},…,Y_{n}}, partitioning the objects into n groups, the Rand index is computed as follows:
where a is the number of pairs in S that are grouped in the same cluster in X and in the same cluster in Y, while b is the number of pairs in S that are grouped in different clusters in X and in different clusters in Y. The denominator is the number of pairs of objects in S, and is equal to \(\frac {n(n1)}{2}\). In other words, the Rand index measures the percentage of pairs that have been clustered in the same way in both clusterings, over all possible pairs of objects. The Rand index takes values from 0 to 1, with 1 meaning that the two clusterings are the same, while 0 means that the two clusterings are completely different.
The results of the multiple clusterings are summarized in Table 4. With the exception of the VFH descriptor, all individual 3D descriptors manage to produce clusterings that are very close to the ground truth one, achieving accuracies from 77 to 89.5%. The accuracies of the 3D descriptors are larger than the accuracy achieved with the method using the TMscore for clustering. This shows that exploiting the 3D structural information encoded in 3D descriptors commonly used in the area of object recognition achieves an improvement compared to using traditional structural information. Regarding the individual 3D descriptors, the RSD descriptor achieves the highest accuracy. The VFH descriptor achieves the least accuracy compared to the other descriptors. This can be attributed to the fact that it is a global type of descriptor, hence a lot of information regarding local variations of the points in the protein is discarded, thus losing discriminating capacity.
The table also includes the results produced using the combined distance matrix, considering either equal weights for the combination or weights based on the silhouette width. The combination of descriptors achieves the maximum accuracy for all types of clustering considered, reaching 89.7% compared to the ground truth for the kmedoids clustering. This demonstrates that combining 3D descriptors manages to produce clusterings that are more accurate than using the individual descriptors. This is an expected result, since inaccuracies of one clustering can be filtered out when considering multiple clusterings at once.
The comparative analysis of the clustering methods also demonstrated that the agglomerative and DBScan methods achieved lower accuracy than the kmedoids method. For this reason, the kmedoids method is selected for the rest of the present study.
In the results of Table 4, a fixed number of 6 clusters was considered for clustering using any of the descriptor types; it was chosen as this is the number of stereotyped subsets existing in the annotated dataset. However, some subsets may be further subdivided into smaller categories, due to finer differences between their BcR IGs, which were not reflected in the annotation by experts. In order to compensate for this, the same experiments were repeated, but using an optimal number of clusters for each type of descriptor. This optimal number of clusters was computed based on an internal measure of cluster compactness, namely the average silhouette width, as described in the “Methods” section. Considering a single type of descriptor, several clusterings were computed, using a range of number of clusters from 3 up to 25, in order to determine the number of clusters that achieves the largest average silhouette width. An example of such a determination of the optimal number of clusters is depicted in Fig. 2, for the FPFH descriptor, where the number of clusters achieving the largest average silhouette width is 9.
Table 5 summarizes the clustering accuracy results for the annotated dataset, when using the optimal number of clusters, computed separately for each descriptor type. Allowing the number of clusters to vary provides more freedom to the clustering algorithm to cluster the data based on their intrinsic clusters. This allows smaller clusters to emerge, resulting in larger numbers of clusters than the ground truth ones, ranging from 7 to 9. However, this freedom allows the clustering algorithms to group the proteins more accurately, thus resulting in higher values for the Rand index. Using this type of analysis with the silhouettebased combination of the descriptors achieves a 92.2% accuracy with respect to the ground truth protein separation.
As an illustration of the clustering performance, Fig. 3 presents the computed clusters graphically. Each vertical bar corresponds to one of the computed clusters. Each bar is constructed from small rectangles, each one representing a single protein. The higher a vertical bar, the more proteins are contained in the corresponding cluster. The colors assigned to each protein correspond to the different ground truth subsets. It can be observed that, most of the proteins have been correctly clustered, with few exceptions. Moreover, the clustering method discovered 7 clusters, instead of 6, splitting stereotyped subset #2 (greenblue color) into two clusters (indexes 4 and 7). The reason behind this separation is probably the pattern of somatic mutations in the immunoglobulin heavychain variable region gene (IGHV).
Clustering of all BcR IGs
The procedure followed for clustering the annotated dataset was repeated, this time using the whole BcR IG protein dataset, including both stereotyped (annotated) and nonstereotyped (non annotated) cases. For each type of descriptor, the optimal number of clusters was first determined, using the maximum average silhouette width method. Then, the proteins were clustered using the kmedoids method with the optimal number of clusters.
The performance of the various clusterings was evaluated using two types of measures. The first is the average silhouette width itself, which is a measure of the cluster compactness and separation. In general, clustering is based on the assumption that the underlying data form compact clusters of similar characteristics. Larger average silhouette width means that the result of a clustering algorithm consists of compact clusters which are well separated from each other, i.e. probably close to the actual data distribution. A small average silhouette width means e.g. that one of the clusters discovered by the clustering algorithm could be separated in two clusters, or that some of the discovered clusters could be merged together. The average silhouette width is an internal evaluation measure, in the sense that it uses only information contained in the dataset, without assuming any knowledge of ground truth class labels or clusterings.
The second type of evaluation measure is the Rand index, which is an external measure, in the sense that it makes use of ground truth knowledge. The evaluation using the Rand index is similar to the evaluation of the annotated dataset in the previous section, by comparing the produced clusterings to the ground truth clustering. However, only the annotated BcR IG were used for the computation of the Rand index. In other words, after computing a clustering of all proteins, both annotated and unannotated, we wanted to evaluate how well they have been clustered by examining the clustering distribution of the annotated ones, within the whole clustering. The results of the clustering are summarized in Table 6.
To evaluate the quality of the applied clustering methods of BcR IG 3D models, we included in our cohort 137 cases originating from 6 CLL stereotyped subsets namely subsets #1, #2, #4, #6, #7 and #8. The reason for this approach was that (i) stereotyped, highly homologous BcR IG primary sequences are anticipated to produce overall similar 3D structures, hence providing a reference for evaluating the developed workflow; and, (ii) these subsets are well characterized in terms of both biological and clinical properties [38]. Subset size distribution was as described in Table 3.
Next, we assessed the reasons behind the separation of subsets #1 and #2 into 2 different clusters each. As mentioned above, stereotyped subset #1 sequences express different IGHV genes that belong to the same phylogenetic clan. Indeed, the utilization of a different IGHV gene was most likely the reason behind the separation of subset #1 models. In more detail, 7/8 (88%) of subset #1 models assigned to cluster 3 expressed the IGHV5101 gene, whereas 27/29 (93%) of subset #1 models in cluster 6 utilized genes belonging to the IGHV1 gene subgroup. In regard to subset #2 models, the reason behind the separation was probably the pattern of somatic mutations within the IGHV gene. More specifically, 9/15 (66%) of subset #2 models from cluster 4 exhibited the presence of somatic mutations within the FR1 region of the Ig heavy chain sequence. In contrast, none of the subset #2 models in cluster 9 carried such a mutation.
We assessed the efficacy of each individual clustering algorithm as well as the “combined” method regarding their potential of biological significance through evaluating the distribution of these 137 stereotyped BcR IG 3D models across different clusters. At the level of individual descriptors, the best results were observed in the case of the 3DSC clustering method, yet, the most robust results were obtained through the combined approach. In detail, the combined method led to the assignment of stereotyped BcR IG 3D models in 9/15 clusters with the following distribution, as also illustrated in Fig. 4:

cluster 1 contained 18/22 (82%) subset #4 models,

cluster 2 contained 11/11 (100%) subset #8 models along with a single (4%) subset #4 model,

cluster 3 contained 8/38 (21%) subset #1 models,

cluster 4 contained 15/42 (36%) subset #2 models,

cluster 5 contained 9/12 (75%) subset #7 models,

cluster 6 contained 29/38 (76%) subset #1 models,

cluster 7 contained 12/12 (100%) subset #6 models and a single (8%) subset #7 model,

cluster 8 contained 2/12 (17%) subset #7 models, and finally,

cluster 9 contained 27/42 (64%) subset #2 models, a single (3%) subset #1 model as well 3/22 (14%) subset #4 models.
Therefore, we focused our analysis on the results obtained with the combined method. Relevant to mention, our cohort of nonstereotyped IG models was representative of CLL in terms of BcR IG heavy and light gene repertoire properties and, thus, largely informative. In specific, IGHV3 gene subgroup cases predominated (395/788, 50.1%) followed by IGHV4 (186/788, 23.6%) and IGHV1 cases (151/788, 19.2%). In regard to the IG light chain expression, around 2/3 of cases expressed kappa light chains (507/788, 64.3%), as reported for CLL. As mentioned above, the "combined" clustering method ended up with the formation of 15 clusters. The distribution of BcR IG 3D models was not equal among different clusters whose size ranged from 26 to 98 models (average: 52.5, median: 53 models). Interestingly, cluster 9 consisted exclusively of 31 stereotyped models, whereas clusters 1015 did not contain any stereotyped Ig 3D models.
In detail, cluster 1 besides containing the majority (82%) of subset #4 models included mostly (38/40, 95%) nonstereotyped BcR IG 3D models of the IGHV4 gene subgroup. From these, most models (25/38, 66%) expressed the IGHV434 gene, as in the case of subset #4 models. Cluster 2 comprised all subset #8 models and a single (4%) #4 model. This was not unexpected, since both subsets express IGHV4 subgroup genes (IGHV439 versus IGHV434, respectively) and carry the IgG heavy chain isotype, in itself a rarity for CLL [54]. Besides stereotyped models, cluster 2 contained IGHV4 models (51/59, 86%) from which the majority expressed the IGHV439 gene, which is characteristic for subset #8. Clusters 3 and 6 contained subset #1 models as well as nonstereotyped models utilizing IGHV genes from the same phylogenetic clan (Clan I: IGHV gene subgroups 1, 5 and 7) with the frequencies being 100% and 84%, respectively. Clusters 4 and 9 consisted of subset #2 models. Nonstereotyped models assigned to cluster 4 utilized IGHV genes from the IGHV3 subgroup (85%) with the most frequent gene being IGHV321 (7/51, 14%) as in subset #2, whereas cluster 9 did not contain any nonstereotyped models. Cluster 5 contained the majority of subset #7 models (75%). Nonstereotyped models in this cluster expressed BcR IG utilizing IGHV1 subgroup genes, predominated by IGHV169, the hallmark of subset #7. Cluster 7 contained all subset #6 models and a single model of subset #7: a possible explanation for this is that both subsets express the IGHV169 gene. Again, nonstereotyped models mostly expressed the IGHV169 gene. Cluster 8 contained a small fraction (17%) of #7 models. In this case, the same trend was not followed and most nonstereotyped models expressed IGHV3 genes with the most prominent being the IGHV323 gene. Clusters 10, 11 and 1315 comprised IGHV3expressing models yet different genes predominated in each cluster: IGHV323, IGHV37, IGHV330, IGHV348 and IGHV330, respectively. Finally, cluster 12 contained models of the IGHV4 subgroup with the most frequent gene being IGHV434. In terms of light chain usage, we observed a dominance of either the kappa (clusters 13, 58, 1012, 14) or the lambda (clusters 4, 13, 15) light chain.
According to our results, the clustering of BcR IG 3D models reflected to very great extent the classification of IG molecules based on the primary sequences of both the heavy and the light chains of the Ig molecule. Indeed, each individual cluster was characterized by the predominance of a single IGHV gene subgroup and a specific light chain isotype.
Discussion
In this work, a novel method of clustering BcR IG protein 3D structures is introduced. The method underlines the significance of the combination between classic bioinformatics methods with 3D descriptors that goes beyond the realm of bioinformatics. The proposed methodology relies on a combination approach based on local descriptors. Two approaches of combination have been used, one using equal weights and another using silhouette widths as weights. Both performances perform better than using the individual descriptors. However, a clear comparison between the two approaches cannot be made yet. When clustering the 137 annotated protein structures, the silhouettebased approach achieved higher accuracy, while the equal weights approach achieved slightly higher accuracy when clustering the 925 structures. Further investigation on the suitability of each approach for different datasets is considered as future work.
Methodology evaluation succeeded through the dataset separation to ground truth and test set. The ground truth was formed by the wellestablished CLL stereotyped subsets and the test set from the unlabeled structures that result from CPIGS methodology. As final analysis level, a clustering of significant external accuracy and internal quality resulted. Overall, our findings support that the innovative workflow described here enables robust clustering of 3D models produced from BcR IG primary sequences from patients with CLL. Furthermore, they indicate that CLL classification based on stereotypy of BcR IG primary sequences is likely also verified at the IG 3D structural level. More generally, this approach can be implemented to the analysis of BcR IG amino acid sequences in any domain of immunology ranging from normal, autoreactive and malignant B cell populations.
Conclusions
A novel BcR IG protein 3D structure comparison technique is proposed for determining the local structural similarity between the 3D models. The method’s generalizability was demonstrated by applying it to two different datasets: one labeled, formed by 137 protein structures that belong to six wellestablished CLL stereotyped subsets, and one mainly unlabeled, formed by 925 (including the cases that belong to stereotyped subsets) protein structures.
Local and global 3D descriptors where tested and the optimal combination of the local descriptors was selected, based on their performance regarding the average silhouette width. The combination of the localbased descriptors derived from the structurally aligned parts is used to compute an overall distance matrix, which is then used as input for the clustering procedure. The combined descriptor presented Rand Index 89.7% and 92.2% in clustering the labeled data to six and the optimal number of clusters, respectively. The higher accuracy in the optimal cluster number is justified by the biological meaning in data. Additionally, the clustering results of the unlabeled data revealed 90.8% accuracy in the optimal cluster numbers. These results support that the innovative workflow described here enables robust clustering of 3D models produced from BcR IG sequences from patients with CLL. The established methodology can be expanded in different types of 3D protein structures.
The selection of the appropriate 3D descriptor is an issue worth studying further. In future work, methods of time series analysis in combination to 3D descriptors will be examined. More specifically the implementation of Fréchet distance and Dynamic Time Warping will be examined on estimating the distance between the 3D models in combination to 3D descriptors.
Abbreviations
 3DSC:

3D shape context
 CLL:

Chronic lymphocytic leukemia
 FPFH:

Fast point feature histogram
 IGHV:

Immunoglobulin heavy chain variable region
 PIGS:

Prediction of immunoglobulin structure
 RSD:

Radiusbased surface descriptor
 VFH:

Viewpoint Feature Histogram
References
 1
Bender A, Glen RC. Molecular similarity: a key technique in molecular informatics. Org Biomol Chem. 2004; 2(22):3204–18.
 2
Webb B, Sali A. Comparative protein structure modeling using MODELLER. Curr Protoc Bioinforma. 2014; 47(1):5–6.
 3
Axenopoulos A, Rafailidis D, Papadopoulos G, Houstis EN, Daras P. Similarity search of flexible 3d molecules combining local and global shape descriptors. IEEE/ACM Trans Comput Biol Bioinforma. 2016; 13(5):954–70.
 4
Murzin AG, Brenner SE, Hubbard T, Chothia C. Scop: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol. 1995; 247(4):536–40.
 5
Knudsen M, Wiuf C. The cath database. Hum Genomics. 2010; 4(3):207.
 6
Csaba G, Birzele F, Zimmer R. Systematic comparison of scop and cath: a new gold standard for protein structure analysis. BMC Struct Biol. 2009; 9(1):23.
 7
Sillitoe I, Dawson N, Thornton J, Orengo C. The history of the cath structural classification of protein domains. Biochimie. 2015; 119:209–17.
 8
Li Z, Natarajan P, Ye Y, Hrabe T, Godzik A. Posa: a userdriven, interactive multiple protein structure alignment server. Nucleic Acids Res. 2014; 42(W1):240–5.
 9
Liu YS, Li Q, Zheng GQ, Ramani K, Benjamin W. Using diffusion distances for flexible molecular shape comparison. BMC Bioinformatics. 2010; 11(1):480.
 10
Arenas AF, Salcedo GE, Montoya AM, GomezMarin JE. Msca: a spectral comparison algorithm between time series to identify proteinprotein interactions. BMC Bioinformatics. 2015; 16(1):152.
 11
Srivastava S, Lal SB, Mishra D, Angadi U, Chaturvedi K, Rai SN, Rai A. An efficient algorithm for protein structure comparison using elastic shape analysis. Algoritm Mol Biol. 2016; 11(1):27.
 12
Kamburov A, Lawrence MS, Polak P, Leshchiner I, Lage K, Golub TR, Lander ES, Getz G. Comprehensive assessment of cancer missense mutation clustering in protein structures. Proc Natl Acad Sci. 2015; 112(40):5486–95.
 13
Jiang M, Xu Y, Zhu B. Protein structure–structure alignment with discrete fréchet distance. J Bioinforma Comput Biol. 2008; 6(01):51–64.
 14
Ballester PJ, Richards WG. Ultrafast shape recognition to search compound databases for similar molecular shapes. J Comput Chem. 2007; 28(10):1711–23.
 15
Ballester PJ, Richards WG. Ultrafast shape recognition for similarity search in molecular databases. In: Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. 463, 2081. The Royal Society;2007.1307–1321.
 16
Bock M, Cortelazzo G, Ferrari C, Guerra C. Identifying similar surface patches on proteins using a spinimage surface representation. In: Combinatorial Pattern Matching. Springer: 2005. p. 29–99.
 17
Ankerst M, Kastenmüller G, Kriegel HP, Seidl T. 3d shape histograms for similarity search and classification in spatial databases. In: International Symposium on Spatial Databases. Springer: 1999. p. 207–26.
 18
Kinoshita K, Nakamura H. Identification of protein biochemical functions by similarity search using the molecular surface database efsite. Protein Sci. 2003; 12(8):1589–95.
 19
Furuya T, Ohbuchi R. Dense sampling and fast encoding for 3d model retrieval using bagofvisual features. In: Proceedings of the ACM International Conference on Image and Video Retrieval. ACM: 2009. p. 26.
 20
Chen DY, Tian XP, Shen YT, Ouhyoung M. On visual similarity based 3d model retrieval, vol. 22. In: Computer Graphics Forum. Wiley Online Library: 2003. p. 223–32.
 21
Chen BY, Honig B. Vasp: a volumetric analysis of surface properties yields insights into proteinligand binding specificity. PLoS Comput Biol. 2010; 6(8):1000881.
 22
Chen BY. Vaspe: Specificity annotation with a volumetric analysis of electrostatic isopotentials. PLoS Comput Biol. 2014; 10(8):1003792.
 23
Amin SR, Erdin S, Ward RM, Lua RC, Lichtarge O. Prediction and experimental validation of enzyme substrate specificity in protein structures. Proc Natl Acad Sci. 2013; 110(45):4195–202.
 24
Wang Y, You Z, Li X, Chen X, Jiang T, Zhang J. Pcvmzm: Using the probabilistic classification vector machines model combined with a zernike moments descriptor to predict protein–protein interactions from protein sequences. Int J Mol Sci. 2017; 18(5):1029.
 25
Wang YB, You ZH, Li LP, Huang YA, Yi HC. Detection of interactions between proteins by using legendre moments descriptor to extract discriminatory information embedded in pssm. Molecules. 2017; 22(8):1366.
 26
Sael L, Li B, La D, Fang Y, Ramani K, Rustamov R, Kihara D. Fast protein tertiary structure retrieval based on global surface shape similarity. Proteins Struct Funct Bioinforma. 2008; 72(4):1259–73.
 27
Ritchie DW, Venkatraman V. Ultrafast fft protein docking on graphics processors. Bioinformatics. 2010; 26(19):2398–405.
 28
Sit A, Kihara D. Comparison of image patches using local moment invariants. IEEE Trans Image Process. 2014; 23(5):2369–79.
 29
Eck S, Wörz S, MüllerOtt K, Hahn M, Biesdorf A, Schotta G, Rippe K, Rohr K. A spherical harmonics intensity model for 3d segmentation and 3d shape analysis of heterochromatin foci. Med Image Anal. 2016; 32:18–31.
 30
Li Z, Geng C, He P, Yao Y. A novel method of 3d graphical representation and similarity analysis for proteins. MATCH Commun Math Comput Chem. 2014; 71:213–26.
 31
Fang Y, Liu YS, Ramani K. Three dimensional shape comparison of flexible proteins using the localdiameter descriptor. BMC Struct Biol. 2009; 9(1):29.
 32
Li B, Lu Y, Li C, Godil A, Schreck T, Aono M, Burtscher M, Chen Q, Chowdhury NK, Fang B, et al. A comparison of 3d shape retrieval methods based on a largescale benchmark supporting multimodal queries. Comp Vision Image Underst. 2015; 131:1–27.
 33
Can T, Wang YF. Ctss: a robust and efficient method for protein structure alignment based on local geometrical and biological features. In: Bioinformatics Conference, 2003. CSB 2003. Proceedings of the 2003 IEEE. IEEE: 2003. p. 169–79.
 34
Mrozek D, BroŻek M, MałysiakMrozek B. Parallel implementation of 3d protein structure similarity searches using a gpu and the cuda. J Mol Model. 2014; 20(2):2067.
 35
Marcatili P, Ghiotto F, Tenca C, Chailyan A, Mazzarello AN, Yan XJ, Colombo M, Albesiano E, Bagnara D, Cutrona G, et al. Igs expressed by chronic lymphocytic leukemia b cells show limited bindingsite structure variability. J Immunol. 2013; 190(11):5771–8.
 36
Sutton LA, Agathangelidis A, Belessi C, Darzentas N, Davi F, Ghia P, Rosenquist R, Stamatopoulos K. Antigen selection in bcell lymphomas—tracing the evidence. vol. 23. In: Seminars in Cancer Biology. Elsevier: 2013. p. 399–409.
 37
Agathangelidis A, Darzentas N, Hadzidimitriou A, Brochet X, Murray F, Yan XJ, Davis Z, van GastelMol EJ, Tresoldi C, Chu CC, et al. Stereotyped bcell receptors in onethird of chronic lymphocytic leukemia: a molecular classification with implications for targeted therapies. Blood. 2012; 119(19):4467–75.
 38
Stamatopoulos K, Agathangelidis A, Rosenquist R, Ghia P. Antigen receptor stereotypy in chronic lymphocytic leukemia. Leukemia. 2017; 31(2):282. Nature Publishing Group.
 39
Rusu RB, Blodow N, Beetz M. Fast point feature histograms (fpfh) for 3d registration. In: Robotics and Automation, 2009. ICRA’09. IEEE International Conference On. IEEE: 2009. p. 3212–7.
 40
Frome A, Huber D, Kolluri R, Bülow T, Malik J. Recognizing objects in range data using regional point descriptors. Comp Vision ECCV 2004. 2004:224–37.
 41
Marton ZC, Pangercic D, Blodow N, Kleinehellefort J, Beetz M. General 3d modelling of novel objects from a single view. In: Intelligent Robots and Systems (IROS), 2010 IEEE/RSJ International Conference On. IEEE: 2010. p. 3700–5.
 42
Rusu RB, Bradski G, Thibaux R, Hsu J. Fast 3d recognition and pose using the viewpoint feature histogram. In: Intelligent Robots and Systems (IROS), 2010 IEEE/RSJ International Conference2 On. IEEE: 2010. p. 2155–62.
 43
Zhang Y, Skolnick J. Tmalign: a protein structure alignment algorithm based on the tmscore. Nucleic Acids Res. 2005; 33(7):2302–9.
 44
Zhang Y, Skolnick J. Scoring function for automated assessment of protein structure template quality. Proteins Struct Funct Bioinforma. 2004; 57(4):702–10.
 45
Rusu RB, Cousins S. 3D is here: Point Cloud Library (PCL). In: IEEE International Conference on Robotics and Automation (ICRA). Shanghai: 2011.
 46
Rusu RB, Marton ZC, Blodow N, Beetz M. Learning informative point classes for the acquisition of object model maps. In: Control, Automation, Robotics and Vision, 2008. ICARCV 2008. 10th International Conference On. IEEE: 2008. p. 643–650.
 47
Hallek M, Cheson BD, Catovsky D, CaligarisCappio F, Dighiero G, Döhner H, Hillmen P, Keating MJ, Montserrat E, Rai KR, et al. Guidelines for the diagnosis and treatment of chronic lymphocytic leukemia: a report from the international workshop on chronic lymphocytic leukemia updating the national cancer institute–working group 1996 guidelines. Blood. 2008; 111(12):5446–56.
 48
Darzentas N, Stamatopoulos K. The significance of stereotyped bcell receptors in chronic lymphocytic leukemia. Hematol Oncol Clin N Am. 2013; 27(2):237–50.
 49
Bystry V, Agathangelidis A, Bikos V, Sutton LA, Baliakas P, Hadzidimitriou A, Stamatopoulos K, Darzentas N. Arrest/assignsubsets: a novel application for robust subclassification of chronic lymphocytic leukemia based on b cell receptor ig stereotypy. Bioinformatics. 2015; 31(23):3844–6.
 50
Marcatili P, Olimpieri PP, Chailyan A, Tramontano A. Antibody modeling using the prediction of immunoglobulin structure (pigs) web server. Nat Protoc. 2014; 9(12):2771–83.
 51
Henikoff S, Henikoff JG. Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci. 1992; 89(22):10915–9.
 52
Messih MA, Lepore R, Marcatili P, Tramontano A. Improving the accuracy of the structure prediction of the third hypervariable loop of the heavy chains of antibodies. Bioinformatics. 2014; 30(19):2733–40.
 53
Marcatili P, Mochament K, Agathangelidis A, Moschonas P, Sutton LA, Yan XJ, Bikos V, Vardi A, Chailyan A, Stavroyianni N, et al.Automated clustering analysis of immunoglobulin sequences in chronic lymphocytic leukemia based on 3D structural descriptors. Blood. 2016; 128(22).
 54
Vardi A, Agathangelidis A, Sutton LA, Chatzouli M, Scarfò L, Mansouri L, Douka V, Anagnostopoulos A, Darzentas N, Rosenquist R, et al.Iggswitched cll has a distinct immunogenetic signature from the common md variant: ontogenetic implications. Clin Cancer Res. 2014; 20(2):323–30.
 55
Ortiz AR, Strauss CE, Olmea O. Mammoth (matching molecular models obtained from theory): an automated method for model comparison. Protein Sci. 2002; 11(11):2606–21.
 56
Zemla A. Lga: a method for finding 3d similarities in protein structures. Nucleic Acids Res. 2003; 31(13):3370–4.
 57
Wrabl JO, Grishin NV. Statistics of random protein superpositions: pvalues for pairwise structure alignment. J Comput Biol. 2008; 15(3):317–55.
 58
Kolodny R, Koehl P, Levitt M. Comprehensive evaluation of protein structure alignment methods: scoring by geometric measures. J Mol Biol. 2005; 346(4):1173–88.
 59
Pandit SB, Skolnick J. Frtmalign: a new protein structural alignment method based on fragment alignments and the tmscore. BMC Bioinformatics. 2008; 9(1):531.
 60
Aung Z, Tan KL. Matalign: precise protein structure comparison by matrix alignment. J Bioinforma Comput Biol. 2006; 4(06):1197–216.
 61
Martínez L, Andreani R, Martínez JM. Convergent algorithms for protein structural alignment. BMC Bioinformatics. 2007; 8(1):306.
 62
Krissinel E, Henrick K. Secondarystructure matching (ssm), a new tool for fast protein structure alignment in three dimensions. Acta Crystallogr D Biol Crystallogr. 2004; 60(12):2256–68.
Funding
This work was supported in part by the H2020 MEDGENET project (H2020TWINN2015 No:692298).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 19 Supplement 14, 2018: Selected articles from the 5th International WorkConference on Bioinformatics and Biomedical Engineering: bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume19supplement14.
Author information
Affiliations
Contributions
KS, KM, IK, EP and DT conceived of the presented idea. KM, IK and EP developed the theory and performed the computations. AA and KS investigated the biological aspects and supervised the findings of this work. LS, XY, VB, AV, NC, CB, RR and PG contributed sequences and/or were involved in the preparation of the final dataset. PM, AC and NO carried out the BcR IG structure prediction. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The author(s) declare(s) that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Polychronidou, E., Kalamaras, I., Agathangelidis, A. et al. Automated shapebased clustering of 3D immunoglobulin protein structures in chronic lymphocytic leukemia. BMC Bioinformatics 19, 414 (2018). https://doi.org/10.1186/s1285901823811
Published:
Keywords
 CLL protein clustering
 3D protein descriptors
 descriptor fusion