 Research article
 Open access
 Published:
Votingbased consensus clustering for combining multiple clusterings of chemical structures
Journal of Cheminformatics volume 4, Article number: 37 (2012)
Abstract
Background
Although many consensus clustering methods have been successfully used for combining multiple classifiers in many areas such as machine learning, applied statistics, pattern recognition and bioinformatics, few consensus clustering methods have been applied for combining multiple clusterings of chemical structures. It is known that any individual clustering method will not always give the best results for all types of applications. So, in this paper, three voting and graphbased consensus clusterings were used for combining multiple clusterings of chemical structures to enhance the ability of separating biologically active molecules from inactive ones in each cluster.
Results
The cumulative votingbased aggregation algorithm (CVAA), clusterbased similarity partitioning algorithm (CSPA) and hypergraph partitioning algorithm (HGPA) were examined. The Fmeasure and Quality Partition Index method (QPI) were used to evaluate the clusterings and the results were compared to the Ward’s clustering method. The MDL Drug Data Report (MDDR) dataset was used for experiments and was represented by two 2D fingerprints, ALOGP and ECFP_4. The performance of votingbased consensus clustering method outperformed the Ward’s method using Fmeasure and QPI method for both ALOGP and ECFP_4 fingerprints, while the graphbased consensus clustering methods outperformed the Ward’s method only for ALOGP using QPI. The Jaccard and Euclidean distance measures were the methods of choice to generate the ensembles, which give the highest values for both criteria.
Conclusions
The results of the experiments show that consensus clustering methods can improve the effectiveness of chemical structures clusterings. The cumulative votingbased aggregation algorithm (CVAA) was the method of choice among consensus clustering methods.
Background
Chemoinformatics, as defined by Brown [1], is the collection, representation and organisation of chemical data in order to create chemical information, which is applied to create chemical knowledge. It has been used for the process of drug discovery and design, especially in the lead identification and optimisation process, which is known as HighThroughput Screening (HTS).
According to Brown and Martin [2], the advent of highthroughput biological screening methods has given pharmaceutical companies the ability to screen many thousands of compounds in a short time. However, there are many hundreds of thousands of compounds available both inhouse and from commercial vendors. Whilst it may be feasible to screen many, or all, of the compounds available, this is undesirable for reasons of cost and time and may be unnecessary if it results in the production of some redundant information. Therefore, there has been a great deal of interest in the use of compound clustering techniques to aid in the selection of a representative subset of all the compounds available [3]. Given a clustering method, which can group structurally similar compounds together, and application of the similar property principle[4], which states that structurally similar molecules will exhibit similar physicochemical and biological properties, implies that the selection, or synthesis, and testing of representatives from each cluster produced from a set of compounds should be sufficient to understand the structureactivity relationships of the whole set, without the need to test them all. An appropriate clustering method will, ideally, cluster all similar compounds together whilst separating active and inactive compounds into different sets of clusters [2].
The main objective of clustering is to subdivide data objects into smaller groups known as clusters so that each group exhibits a high degree of intracluster similarity and intercluster dissimilarity [5]. Many different types of clustering techniques for chemical structures have been used in the literature [6–13].
Brown and Martin [2] considered the Ward’s clustering to be the most efficient method in clusterbased compound selection. However, as it is known, there is no clustering method capable of correctly finding the underlying structure for all data sets. So, the idea of combining different clustering results (consensus clustering) is considered as an alternative approach for improving the quality of the individual clustering algorithms [14].
Consensus clustering involves two main steps: (i) partitions generation and (ii) consensus function. In the first step, many partitions will be generated (the collection of partitions is called ensemble). There are no constraints about how the partitions must be generated. In the partitions generation step, many mechanisms can be applied including the using of: (i) different data representations, (ii) different individual clustering methods, (iii) different parameters initialisation for clustering methods and (iv) data resampling. In the second step, there are two main approaches, i.e. the objects cooccurrencebased and the median partitionbased approaches. Voting and graph based consensus clusterings are widely used for the first approach.
Topchy et al.[15] and Fred and Jain [16] summarised the main advantages of using consensus clustering in the following terms:

Robustness: The combination process must have better average performance than the single clustering algorithms.

Consistency: The result of the combination should be somehow, very similar to all combined single clustering algorithm results.

Novelty: Cluster ensembles must allow finding solutions unattainable by single clustering algorithms.

Stability: Results with lower sensitivity to noise and outliers.
In chemoinformatics, it is most unlikely that any single method will yield the best classification under all circumstances, even if attention is restricted to a single type of application [17]. Furthermore, the consensus scoring (data fusion) has been successfully used in chemoinformatics and, in particular, for virtual screening [18–25].
Over the last few years, data fusion has become accepted as a simple way of enhancing the performance of existing systems for ligandbased virtual screening by combining the results of two or more screening methods. In some cases, the fused search may even be better than the best individual screening method when averaged over large numbers of searches [20].
Chu et al.[17] used consensus clustering methods on sets of chemical compounds represented by 2D fingerprints (ECFP_4), and concluded that consensus methods can outperform the Ward's method, the current standard clustering method for chemoinformatics applications. However, based on the implemented methods, it was not the case if the clustering is restricted to a single consensus method. In this paper, we examined the use of voting and graphbased consensus clustering methods for combining multiple clusterings of chemical structures with different distance measures in order to improve the effectiveness of chemical structures clustering.
Experimental
Molecular fingerprints
For the clustering experiments, two molecular fingerprints were developed by Scitegic’s Pipeline Pilot software [26]. The first one was 120bit ALOGP, which includes octanolwater partitioning coefficient based on Ghose and Crippen’s method [27, 28]. ALOGP atom type code is generated based on the molecular hydrophobicity (lipophilicity), usually quantified as log P (the logarithm of 1octanol/water partition coefficient), which is an important molecular characteristic in drug discovery [28]. The second descriptor was the Scitegic extendedconnectivity fingerprints (1024 ECFP_4). The first character E in the fingerprint name denotes the atom abstraction method used to assign initial atom code which is derived from the number of connections to an atom, the element type, the charge and the atomic mass [29].
Dataset
Experiments were conducted over the most popular chemoinformatics databases: the MDL Drug Data Report database [30]. This database consists of 102516 molecules which was described by Hert et al.[21] and used for consensus scoring that combined the results of different similarity searches of a chemical database. Also, it has been used for many virtual screening experiments [31–33]. According to Hert et al.[26], the subset dataset was chosen from MDDR database which is quite disparate in nature, some of the molecules being structurally homogeneous (e.g., rennin and HIV1 protease inhibitors) while others were structurally diverse (e.g., cyclooxygenase and protein kinase C inhibitors); the diversity was estimated by the mean pairwise Tanimoto similarity across each set of active molecules (activity class). The calculations of pairwise Tanimoto similarity were conducted using Pipeline Pilot software. The MDDR dataset contains eleven activity classes (8294 molecules) and the details of this dataset are listed in Table 1. Each row in the table contains an activity class, the number of molecules belonging to the class and the diversity of the class.
Ensemble generation
Every consensus clustering method is made up of two steps: (i) partitions generation and (ii) consensus functions. For the purposes of this paper, the partitions were generated by a single run of multiple individual clustering algorithms (singlelinkage, complete linkage, average linkage, weighted average distance, Ward and Kmeans methods). Every individual clustering used six distance measures in order to generate different ensembles. The thresholds of 500, 600, 700, 800, 900 and 1000 were used to generate partitions with different sizes (number of clusters). The same process was done for each 2D fingerprint in order to study the effectiveness of consensus clusterings on different molecular representations. The distance measures that were used with each clustering technique were Correlation, Cosine, Euclidean, Hamming, Jaccard and Manhattan.
Methods
Graphbased consensus clustering
Two graphbased consensus clustering algorithms, proposed by Strehl and Gosh [34], were used to obtain the consensus partition from ensembles generated in the previous step. The two algorithms were developed based on transforming the set of clusterings into a hypergraph representation. The first algorithm is Clusterbased Similarity Partitioning Algorithm (CSPA) in which a clustering signifies a relationship between objects in the same cluster and can thus be used to establish a measure of pairwise similarity. Because of this similarity measure, CSPA is also categorized under consensus similarity matrix methods. The second algorithm is the HyperGraph Partitioning Algorithm (HGPA) in which the cluster ensemble problem is formulated as partitioning the hypergraph by cutting a minimal number of hyperedges. Both algorithms were coded by the published cluster ensemble package that is available on (http://www.strehl.com).
For CSPA, the similarity matrix is generated so that each two objects have a similarity of 1 if they are in the same cluster and 0 otherwise. The process is repeated for each clustering method. A n× n binary similarity matrix S is created where n is the total number of objects in the dataset. The entries of S are divided by r, which is the number of clustering methods. Then, the similarity matrix is used to recluster the objects using any reasonable similaritybased clustering algorithm. Here, we view the similarity matrix as graph (vertex = object, edge weight = similarity) and cluster it using graph partitioning algorithm METIS [35] because of its robust and scalable properties in order to obtain the consensus partition.
The HGPA portions the hypergraph directly. This is done by removing the lower number of hyperedges. All hyperedges have the same weight and are searched by cutting the minimum possible number of hyperedges that partition the hypergraph in k connected components of approximately the same dimension. For the implementation of this method, the hypergraphs partitioning package HMETIS [36] was used.
Votingbased consensus clustering
The cumulative votingbased aggregation algorithm consists of two steps; the first one is to obtain the optimal relabeling for all partitions, which is known as the voting problem. Then, the votingbased aggregation algorithm is used to obtain the aggregated (consensus) partition. The votingbased aggregation algorithm described by Ayed and Kamel [37, 38] is modified to be used in this paper.
Let χ denote a set of n data objects, and let a partition of χ into k clusters be represented by an n×k matrix U such that ∑ _{q = 1}^{k}ujq = 1, for ∀ j. Let u = {U_{ i }}_{i = 1}^{b} denote an ensemble of partitions. The votingbased aggregation problem is concerned with searching for an optimal relabeling for each partition V^{i} with respect to representative partition U^{0} (with k^{0} clusters) and for a central aggregated partition denoted as Ū that summarises the ensemble partitions. The matrix of coefficients W^{i}, which is a k^{i} × k^{0} matrix of w_{ lq }^{i} coefficients, is used to obtain the optimal relabeling for ensemble partitions.
In this paper, the fixedreference approach is used, whereby an initial reference partition is used as a common representative partition for all the ensemble partitions and remains unchanged throughout the aggregation process. Instead of selecting random partition, the partition that is generated by the method, which showed high ability to separate active from inactive molecules in our experiments, is suggested to be the reference partition U^{0}; and this method is the Ward’s clustering (the current standard clustering method for Chemoinformatics applications). The cumulative votingbased aggregation algorithm is described as follows:
Cumulative Votingbased Aggregation Algorithm

1.
select a partition U^{i} ∈u which is generated by the Ward’s method and assign to U^{0}

2.
for i=1 to b do

3.
W^{i} = (U^{iT}U^{i})^{− 1}U^{i}^{T}U^{0}

4.
V^{i} = U^{i}W^{i}

5.
{\mathrm{U}}^{0}=\frac{i1}{i}{\mathrm{U}}^{0}+\frac{1}{i}{\mathrm{V}}^{i}

6.
end for

7.
Ū = U ^{0}.
Performance evaluation
The results were evaluated based on the effectiveness of the methods to separate active from inactive molecules using two measures: the Fmeasure [39] and Quality Partition Index (QPI) measure [40]. As defined by [17], if the cluster contains n compounds, that a of these are active and that there is a total of A compounds with the chosen Activity. The precision, P, and the recall, R, for that cluster are:
This calculation is carried out on each cluster and the Fmeasure is the maximum value across all clusters.
In addition, according to [13], an active cluster can be defined as a nonsingleton cluster for which the percentage of active molecules in the cluster is greater than the percentage of active molecules in the dataset as a whole. Let p be the number of actives in active clusters, q be the number of inactives in active clusters, r be the number of actives in inactive clusters (i.e., clusters that are not active clusters) and s be the number of singleton actives. The high value occurs when the actives are clustered tightly together and separated from the inactive molecules. The QPI is defined to be:
Results and discussion
The ensembles were generated by running the six individual clusterings, each with the six distance measures. Then, the ensembles were combined using voting and graphbased consensus clustering methods: CVAA, CSPA and HGPA. This process was repeated for each fingerprint (ALOGP and ECFP_4).
The mean of Fmeasure and QPI values were averaged over the eleven activity classes of the dataset. Tables 2, 3, 4 and 5 show the effectiveness of MDDR dataset clustering for ALOGP and ECFP_4 fingerprints. In all tables, the best values for Fmeasure and QPI of consensus clustering methods for each column were boldfaced for ease of reference.
Visual inspection of Fmeasure and QPI values in Tables 2, 3, 4 and 5 enables comparisons to be made between the effectiveness of three consensus clustering methods and the Ward’s method. In addition, two fingerprints were used for the experiments in order to study the effectiveness of consensus clustering on different representations of molecular dataset.
For clustering of MDDR dataset which is represented by ALOGP fingerprint (Tables 2 and 4), the performance of votingbased consensus method (CVAA) outperformed the Ward’s method and the graphbased consensus methods (CSPA and HGPA) using the two criteria: Fmeasure and QPI. The highest Fmeasure values were obtained by using Euclidean distance measure with individual clustering methods in the ensemble generation step. While, using the QPI measure, the highest values were obtained by using the Jaccard distance measure. Moreover, the ensembles generated by the other distance measures showed a better performance of CVAA than Ward and graphbased consensus clustering methods using both criteria.
Similarly, the results in Tables 3 and 5 show that, when ECFP_4 fingerprint is used, the CVAA consensus clustering performed very well and outperformed Ward and graphbased consensus clustering methods using F and QPI measures. The Jaccard distance measure was the method of choice to generate the ensembles, which gives the highest values for both criteria.
Some statistical significance tests (Ttests) were performed to show the improvements achieved by the consensus clustering methods, as shown in Tables 6 and 7. It was found that the performance of votingbased consensus method is statistically significant when using both criteria.
Tables 6 and 7 display a number of parameters: mean value, standard deviation, standard error and significance values for the pairs of the best Fmeasure and QPI values of clustering methods which are ((CVAA, Ward’s method), (CVAA, CSPA) and (CVAA, HGPA)) compared in the paired samples ttest procedure. The pairedsamples ttest procedure compares the means of two variables that represent the same group at different cluster size. Since the paired samples ttest compares the means for the two variables, it is quite useful to know what the mean values are. A low significance value for the ttest (typically less than 0.05) indicates that there is a significant difference that was satisfied between two variables. In Tables 6 and 7, it was noted that the significance field (Sig. (2tailed)) in terms of Fmeasure for AlogP is: CVAAWard’s method (0.000004), CVAACSPA (0.000003) and CVAAHGPA (0.000002), and for ECFP_4: CVAAWard’s method (0.000589), CVAACSPA (0.000391) and CVAAHGPA (0.000357). Similarly, the significance field (Sig. (2tailed)) in terms of QPI values for AlogP is: CVAAWard’s method (0.000199), CVAACSPA (0.004290) and CVAAHGPA (0.013842), and for ECFP_4 is: CVAA Ward’s method (0.000301), CVAACSPA (0.000005) and CVAAHGPA (0.000010). In addition, the significance value is low in Fmeasure and QPI values and the confidence interval for the mean difference does not contain zero. It is then concluded that the consensus clustering method, CVAA, obtained significant results by using Fmeasure and QPI values compared to Ward and graphbased consensus clustering methods. Moreover, the CVAA method is more efficient because the computational complexity of CVAA is O (n) which is better than other consensus clustering methods such as CSPA that has complexity of O (n^{2}), where n is the number of data objects [14].
Conclusions
The experiments results show that consensus clustering methods can improve the effectiveness of chemical structures clusterings. The cumulative votingbased aggregation algorithm CVAA was the method of choice among consensus clustering methods. The performance of CVAA consensus clustering significantly outperforms Ward and graphbased consensus clustering methods (CSPA and HGPA) using F and QPI measures for both ALOGP and ECFP_4 fingerprints, while the graphbased consensus methods outperform the Ward’s method only for ALOGP using QPI measure. The experiments reported here suggest that votingbased consensus clustering can perform very well when the partitions are generated by a single run of multiple individual clusterings that use Jaccard and Euclidean distance measures in the ensemble generation process.
References
Brown FK: Chemoinformatics: what is it and how does it impact drug discovery. Annu Rep Med Chem. 1998, 33: 375384.
Brown RD, Martin YC: Use of structureactivity data to compare structurebased clustering methods and descriptors for use in compound selection. J Chem Inf Comput Sci. 1996, 36: 572584. 10.1021/ci9501047.
Willett P, Winterman V, Bawden D: Implementation of nonhierarchic cluster analysis methods in chemical information systems: selection of compounds for biological testing and substructure search output. J Chem Inf Comput Sci. 1986, 26: 109118. 10.1021/ci00051a005.
Johnson M, Maggiora GM: Concepts and Applications of Molecular Similarity. 1990, New York: Wiley
Everitt BS, Landau S, Leese M: Cluster Analysis. 2001, London: Edward Arnold, 4
Adamson GW, Bush JA: A method for the automatic classification of chemical structures. Information Storage and Retrieval. 1973, 9: 561568. 10.1016/00200271(73)900594.
Downs GM, Barnard JM: Clustering of chemical structures on the basis of twodimensional similarity measures. J Chem Inf Comput Sci. 1992, 32: 644649. 10.1021/ci00010a010.
Willett P: Similarity and Clustering in Chemical Information Systems. 1987, Letchworth: Research Studies Press
Downs GM, Willett P, Fisanick W: Similarity searching and clustering of chemicalstructure databases using molecular property data. J Chem Inf Comput Sci. 1994, 34: 10941102. 10.1021/ci00021a011.
Brown RD, Martin YC: The information content of 2D and 3D structural descriptors relevant to ligand–receptor binding. J Chem Inf Comput Sci. 1997, 37: 19. 10.1021/ci960373c.
Schuffenhauer A, Brown N, Ertl P, Jenkins JL, Selzer P, Hamon J: Clustering and rulebased classifications of chemical structures evaluated in the biological activity space. J Chem Inf Model. 2007, 47 (2): 325336. 10.1021/ci6004004.
Holliday JD, Rodgers SL, Willet P: Clustering files of chemical structures using the fuzzy kmeans clustering method. J Chem Inf Comput Sci. 2004, 44: 894902. 10.1021/ci0342674.
Varin T, Bureau R, Mueller C, Willett P: Clustering files of chemical structures using the Székely–Rizzo generalization of Ward's method. J Mol Graph Model. 2009, 28 (2): 187195. 10.1016/j.jmgm.2009.06.006.
VegaPons S, RuizSchulcloper J: A survey of clustering ensemble algorithms. Int J Pattern Recognit Artificial Intelligence. 2011, 25 (Issue 3): 337372.
Topchy A, Jain AK, Punch W: A mixture model of clustering ensembles. Proc. SIAM Intl. Conf. on Data Mining. 2004
Fred ALN, Jain AK: Combining multiple clustering using evidence accumulation. IEEE Trans Patt Anal Mach Intell. 2005, 27: 835850. 10.1109/TPAMI.2005.113.
Chu CW, Holliday J, Willett P: Combining multiple classifications of chemical structures using consensus clustering. Bioorgan Med Chem. 2012, 20 (18): 53665371. 10.1016/j.bmc.2012.03.010.
Feher M: Consensus scoring for proteinligand interactions. Drug Discov Today. 2006, 11: 421428. 10.1016/j.drudis.2006.03.009.
Salim N, Holliday JD, Willett P: Combination of fingerprintbased similarity coefficients using data FusionJ. Chem Inf Comput Sci. 2003, 43: 435442. 10.1021/ci025596j.
Willet P: Enhancing the effectiveness of ligandbased virtual screening using data fusion. QSAR Comb Sci. 2006, 25: 11431152. 10.1002/qsar.200610084.
Hert J, Willett P, Wilton DJ, Acklin P, Azzaoui K, Jacoby E, Schuffenhauer A: New methods for ligandbased virtualscreening: use of data fusion and machinelearning to enhancethe effectiveness of similarity searching. J Chem Inf Model. 2006, 46: 462470. 10.1021/ci050348j.
Whittle M, Gillet VJ, Willett P: Analysis of data fusion methods in virtual screening: Similarity and group fusion. J Chem Inf Model. 2006, 6: 22062219.
Chen B, Mueller C, Willett P: Combination rules for GroupFusion in similaritybased virtual screening. MolInf. 2010, 29: 533541.
RiveraBorroto OM, MarreroPonce Y, García de la Vega JM, GrauÁbalo RC: Comparison of combinatorial clustering methods on pharmacological data sets represented by machine learningselected real molecular descriptors. J Chem Inf Model. 2011, 51 (12): 30363049. 10.1021/ci2000083.
Svensson F, Karlen A, Skold C: Virtual Screening DataFusion Using Both Structure and LigandBased Methods. 2011, Model: J. Chem. Inf
Pipeline Pilot software: SciTegic Accelrys Inc. 2008, San Diego: Accelrys Inc website, http://www.accelrys.com/,
Ghose AK, Crippen GM: Atomic physicochemical parameters for threedimensional structuredirected quantitative structureactivity relationships 1. Partition coefficients as a Measure of hydrophobicity. J Comput Chem. 1986, 7: 565577. 10.1002/jcc.540070419.
Ghose AK, Viswanadhan VN, Wendoloski JJ: Prediction of hydrophobic (lipophilic) properties of small organic molecules using fragmental methods: An analysis of ALOGP and CLOGP methods. J Phys Chem A. 1998, 102: 37623772. 10.1021/jp980230o.
Chen L, Li Y, Zaho Q, Peng H, Hou T: ADME evaluation in drug discovery. 10. Predictions of Pglycoprotein inhibitors using recursive partitioning and naive Bayesian classification techniques. Mol Pharm. 2011, 8: 889900. 10.1021/mp100465q.
Sci Tegic Accelrys Inc., the MDL Drug Data Report (MDDR). database is available from at http://www.accelrys.com/ (accessed 1st of November 2012)
Abdo A, Chen B, Mueller C, Salim N, Willett P: Ligandbased virtual screening using bayesian networks. J Chem Inf Model. 2010, 50: 10121020. 10.1021/ci100090p.
Abdo A, Salim N: New fragment weighting scheme for the bayesian inference network in ligandbased virtual screening. J Chem Inf Model. 2011, 51: 2532. 10.1021/ci100232h.
Abdo A, Saeed F, Hentabli H, Ali A, Salim N, Ahmed A: Ligand expansion in ligandbased virtual screening using relevance feedback. J ComputAided Mol Des. 2012, 26: 279287. 10.1007/s1082201295434.
Strehl A, Ghosh J: Cluster ensembles—a knowledge reuse framework for combining multiple partitions. J Machine Learning Res. 2002, 3: 583617.
Karypis G, Kumar V: A fast and high quality multilevel scheme for partitioning irregular graphs. SIAM J Scient Comput. 1998, 20: 359392. 10.1137/S1064827595287997.
Karypis G, Aggarwal R, Kumar V, Shekhar S: Multilevel hypergraph partitioning: Application in VLSI domain. Proceedings of the 34th annual Design Automation Conference. 1997, 526529. ACM
Ayad HG, Kamel MS: Cumulative voting consensus method for partitions with a variable number of clusters. IEEE Trans Pattern Anal Mach Intell. 2008, 30 (1): 160173. January
Ayad HG, Kamel MS: On votingbased consensus of cluster ensembles. Patt Recogn. 2010, 43: 19431953. 10.1016/j.patcog.2009.11.012.
Van Rijsbergen CJ: Information Retrieval. 1979, London: Butterworths, 2
Varin T, Saettel N, Villain J, Lesnard A, Dauphin F, Bureau R, Rault SJ: 3D Pharmacophore, hierarchical methods, and 5HT4 receptor binding data. Enzyme Inhib Med Chem. 2008, 23: 593603. 10.1080/14756360802204748.
Acknowledgements
This work is supported by the Ministry of Higher Education (MOHE) and Research Management Centre (RMC) at the Universiti Teknologi Malaysia (UTM) under Research University Grant Category (VOT Q.J130000.7128.00H72). We also would like to thank MISMOHE for sponsoring the first author.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
FS is a PhD candidate and performed the experiments under the supervision of NS and AA. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Saeed, F., Salim, N. & Abdo, A. Votingbased consensus clustering for combining multiple clusterings of chemical structures. J Cheminform 4, 37 (2012). https://doi.org/10.1186/17582946437
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/17582946437