Skip to main content

Quantitatively integrating molecular structure and bioactivity profile evidence into drug-target relationship analysis

Abstract

Background

Public resources of chemical compound are in a rapid growth both in quantity and the types of data-representation. To comprehensively understand the relationship between the intrinsic features of chemical compounds and protein targets is an essential task to evaluate potential protein-binding function for virtual drug screening. In previous studies, correlations were proposed between bioactivity profiles and target networks, especially when chemical structures were similar. With the lack of effective quantitative methods to uncover such correlation, it is demanding and necessary for us to integrate the information from multiple data sources to produce an comprehensive assessment of the similarity between small molecules, as well as quantitatively uncover the relationship between compounds and their targets by such integrated schema.

Results

In this study a multi-view based clustering algorithm was introduced to quantitatively integrate compound similarity from both bioactivity profiles and structural fingerprints. Firstly, a hierarchy clustering was performed with the fused similarity on 37 compounds curated from PubChem. Compared to clustering in a single view, the overall common target number within fused classes has been improved by using the integrated similarity, which indicated that the present multi-view based clustering is more efficient by successfully identifying clusters with its members sharing more number of common targets. Analysis in certain classes reveals that mutual complement of the two views for compound description helps to discover missing similar compound when only single view was applied. Then, a large-scale drug virtual screen was performed on 1267 compounds curated from Connectivity Map (CMap) dataset based on the fused similarity, which obtained a better ranking result compared to that of single-view. These comprehensive tests indicated that by combining different data representations; an improved assessment of target-specific compound similarity can be achieved.

Conclusions

Our study presented an efficient, extendable and quantitative computational model for integration of different compound representations, and expected to provide new clues to improve the virtual drug screening from various pharmacological properties. Scripts, supplementary materials and data used in this study are publicly available at http://lifecenter.sgst.cn/fusion/.

Background

To comprehend relationship between intrinsic characteristics of chemical compound and the compound interaction with protein target is an essential task to evaluate potential protein-binding function for virtual drug screening. Similarity relationship between compounds can be characterized differently, depending on different aspects of features to be measured. The similarity measurement of small molecules has been the focus of essentially every compound-based approach to design or identify novel drug candidates [1]. However, in the process of novel drug screening, the representation of a compound varies dramatically, which results in different similarity measurements. Such similarity difference has given rise to distinct candidate compound similarity ranking lists with only generally about 15% overlap [1]. It is demanding and necessary if information from multiple data sources can be integrated together to produce a comprehensive representation and assessment of similarity relationship between small molecules [2], thus expected to boost the results of virtual drug screening.

Generally, the drug candidates are related to specific targets. The investigation on the nature of target-specific structure–activity relationships of molecules should be based on the available data sources concerning structure, activity and target-binding information from a comprehensive and integrative perspective. Fortunately, public resources are in a rapid growth both in the quantity of data and in the type of data-generating, which provide us a great chance to further mine the relationship between compounds and their targets. Besides the classic representations of small molecules, like various fingerprints characterizing compound chemical structure, public high-throughput experimental data representing bioactivity of compounds are boosting with the development of online database, including PubChem (http://pubchem.ncbi.nlm.nih.gov/) [3], Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) [4] and DrugBank (DrugBank, http://drugbank.ca/) [5] etc., which provides an alternative way for molecule characterization based on bioactivity profiles. Several recent studies on the relationship between different compound features claimed that, correlations were proposed between bioactivity profiles and target networks, especially when chemical structures were similar [2, 68]. By simply combining both public repositories of compound targets and compound bioactivity, these studies indicates that comparison of bioactivity profile can provide insight into the mode of actions (MOA) at the molecular level, which will facilitate the knowledge-based discovery of novel compounds. However although various relationship were found between multiple features, no effective quantitative integrating methods was proposed or evaluated to combine these multi-view features. Inspired by previous works, two important and interesting computational issues are needed to investigate: (1) is there a quantitative relationship between compound features (bioactivity profile and structural feature) and compound target that can be specifically described? (2) Since the former works implicated that an integration of multiple compound features may result in a better measurement of target-specific compound similarity rather than only one specific type was adopted, how such integration can be optimized to quantitatively and automatically combine information from various views of compound representations, i.e., structural features, bioactivity features and other more? Hereby in our study, we refer such multiple features description and integration for compound as a multi-view data representation and learning problem, and we aim at presenting a quantitative relationship between target-specific compound similarity and multi-view representations of compound features in an efficient multi-view learning schema.

It should be noted that the term “multi-view learning” was initially presented from 3D-object recognition by the machine learning and graphic communities [9]. Naturally as implicated by its name, multi-view learning combines models from different aspects of one identical entity to obtain an overall and comprehensive representation for further study. Multi-view learning was classically introduced as co-training, a semi-supervised learning procedure to distinguish webpages using two different types of data [10]. Thereafter the concept of integration of different information sources has been developed for years in the field of information retrieval [1113]. On the other side, as an unsupervised-learning method, multi-view clustering algorithms can be divided into two categories in general [14]: (1) Fusion of similarity data by deriving a convex combination of similarities from different views to minimize a given penalty error [15, 16]. (2) Fusion of clustering decision derived from each view separately [17, 18]. In the clustering process, other techniques like canonical correlation analysis (CCA) [19] and matrix factorization [20] were employed to reduce the feature dimension or reconcile clustering groups. These applications of multi-view learning commonly yield better performance than that of single-view learning. In our study, as both the structure and bioactivity information are two distinguished intrinsic features to describe the small molecule, it is natural to investigate the results with the integration of both the chemical space (molecule structure) and genetic space (bioactivity profile) of molecules for a better evaluation of molecular properties and similarity comparison.

In this study, firstly a data set of 37 compounds (in Additional file 1: Table S1) from previous study based on bioactivity profile similarity [6] were adopted. Two similarity matrix characterizing bioactivity profile and structural similarity were calculated. As we would like to investigate the hierarchical structure of similarity among compounds regarding to multiple data sources, rather than only achieve an integrated ranking decision, a similarity fusion method was employed and modified to automatically optimize the weights of the combination of different similarity data. A hierarchy clustering was produced and discussed based on the fused similarity. Then, in order to evaluate the fusion method on the large scale dataset, Connectivity MAP dataset [21] containing 1267 compounds with their gene expression profile and structure fingerprint representation were used to perform drug virtual screen based on similarity searching. The compound-target interaction in these experiments was also analysed and compared quantitatively to demonstrate the benefits introduced by the integration of multiple data representations.

Materials and methods

Algorithm workflow

The workflow of our analysis is illustrated in Figure 1. The intuition behind this workflow is to automatically identify the weights for two molecule representations in fusion under a mathematical optimization framework. Given two similarity matrix P 1 and P 2 , weights α = α 1 , α 2 were to be optimized for a final similarity matrix p = α 1 p 1 + α 2 p 2 . Initially two similarity matrices of different views were used as input after standardization to the z-value and renormalization. Then a two-step alternative minimization was used to obtain the proper weights for the two similarity matrix in fusion. In the first step, given the initial weights α = α 1 , α 2 Cross-entropy between the input matrices and a combined non-negative factorization was minimized by an EM algorithm. In the second step, given the calculated cross-entropy, the weights were calculated by minimizing the object function, i.e. the cross-entropy and entropy of the weight. The two steps iterate until convergence. The final α was used as an ideal weighing vector that obtains balance between weighted sparseness and informativeness. Details are shown below.

Figure 1
figure 1

Workflow of this study. Initially two similarity matrices of different views were used as input after standardization to the z-value and renormalization. Then a two-step alternative minimization was used to obtain the proper weights for the two similarity matrix in fusion. In the first step, given the initial weights α = α 1 , α 2 . cross-entropy between the input matrices and a combined non-negative factorization was minimized by an EM algorithm. In the second step, given the calculated cross-entropy, the weights were calculated by minimizing the object function, i.e. the cross-entropy and entropy of the weight. The two steps iterate until convergence. The final α was used as an ideal weighing vector that obtains balance between weighted sparseness and informativeness.

Dataset

NCI-60 dataset

In our study, the same data set used in Cheng’s work [6] rather than the up-to-date data is applied for equally comparison purpose, in order to illustrate the superiority of target-relationship analysis with similarity fusion from integration of multi-view information. The NCI-60 data set is available in the PubChem BioAssay Database, derived from the bioassays titled “NCI human tumor cell line growth inhibition assay” with relatively sufficient number of tested compounds (more than 16,000). Finally, filtered through 3 rules as Cheng defined [6], 37 small molecules of eligible quality were curated as the final NCI-60 dataset (in Additional file 1: Table S 1).

CMap dataset

In order to demonstrate the performance of the feature integration on the large-scale dataset, similarity fusion was performed on the well-known Connectivity Map dataset. [21] Justin Lamb, et al. had created the first reference collection of gene-expression profiles from cultured human cells stimulated with bioactive small molecules, together with the pattern-matching algorithm to mine these data. To date, CMap contains approximately 7,100 expression profiles representing 1,309 compounds. Some compounds with only expression profiles of HT_HG-U133A_EA Gene chips were not included in this study due to the lack of chip description information. Compared to the former NCI-60 data, the gene expression profile for a compound can also be viewed as a kind of bioactivity representation.

Methods

Test for NCI-60 dataset

Similarity matrix from two views: Bioactivity profile and molecule structure

The pairwise similarities among the 37 molecules are characterized by two similarity matrices in two views. In the view of bioactivity, similarity between two compounds is measured by the Pearson correlation coefficient of the two bioactivity profiles:

r = n Σ A i B i Σ A i Σ B i n Σ A i 2 Σ A i 2 n Σ B i 2 Σ B i 2
(1)

Where n is 37, A i and B i are the log(GI50) values in the ith NCI-60 cell line for the compound A and B, respectively. In the view of molecule structure, commonly-used path-based 1024-bit fingerprint of each compound is calculated via java CDK library to represent the molecular structure, and the similarity of two compounds is measured by the tanimoto-index of the two structural fingerprints:

t = N AB N A + N B - N AB
(2)

where N A N B is the number of features in compound A(B), and N AB is the number of features common to both A and B. Both of the two similarity measurements are in the interval from 0 to 1. It should be noted that correlation coefficient of bioactivity profile below 0 are assign to 0 for two reasons: (1) only very few compounds pairs have a negative correlation coefficient and the minimum is −0.2, which is not significant as an evidence of negative correlation; (2) regarding to the integration analysis of different similarity information, negative correlation brings in no better information of molecular similarity than noise. Finally, as the input for multi-view fusion [15], the two n ñ n similarity matrices S = S ij were standardized as S = S m e a n s / s d s and renormalized to P = S / Σ ij S ij .

Fusion of similarity matrices by expectation-maximization (EM) algorithm

The fused similarity matrix is a convex combination of the L original matrices weighted by vector αwith Σ i a i = 1 and a i 0 . The advantage of this model lies that it can automatically learn the optimized proper weights for each matrix for fusion rather than arbitrary setting the values. This is achieved by an two-step alternative minimization method introduced by T. Lange and J.M. Buhmann [15]. A brief process of the two alternating steps is summarized as follow:

  1. 1.

    Non-negative matrix factorization:

In this study the Non-Negative Matrix Factorization (NMF) is used as one step of the minimization of cross-entropy. The target fused matrix P 0 , 1 n ñ n can be factorized into a product V H t of the n ñ k matrices of V and H. Here the parameter k was assigned to 6 in accordance with Cheng’s number of clustering [6]. It should be noted that the selection of clustering number in a common cluster algorithm is always a non-trivial problem, however in our study, we just set the same cluster number as in the former study for an equally comparison purpose. The computational model proposed here is well extendable to tune the optimal cluster number if any pre-knowledge are unavailable. Then given the fixed weights α of the similarity matrices (initial α of (0.5, 0.5) is used and the value of α is updated in every iteration), we can obtain estimated and using an EM-process which minimize the cross-entropy between P and by updating V and H iteratively.

  1. 2.

    Optimizations of weights for similarity matrices by minimizing the cross-entropy:

Given the estimated factorized matrices, we minimize the cross-entropy between and V H t

Min αVH Σ l α l C P l | | V H t
(3)

regarding to α subject to Σ i a i = 1 and a i 0 , where C P | | Q denotes the cross-entropy of P and Q, and:

C P | | Q = Σ x p x log q x
(4)

Hence the second step becomes a linear program problem.

Since the solution of the linear program would tend to be too sparse that only one of the data source would be chosen to minimize the object function, which is against our intention to combine multiple data source, it is necessary to modify the object function by introducing the entropy of weight α so that both sparseness and informativeness could be taken into account. Since the information quantity provided by the weights vector could be measured by its entropy [15], the modified object function is: Min αVH Σ l α l C P l | | V H t η H α s.t. Σ α = 1 and α i 0 (5)Where H α denotes the entropy of α:

H α = Σ i p α i log p α i
(5)

Then the second step becomes a NLP problem and can be solved with LINDO API 6.1.

Parameter optimization

The parameter η controls the trade-off between sparseness and informativeness. η 0 indicates that entropy will take few significance in the object function, while η indicates that information are taken as the most important factor in the object function, and the weights of different source are evenly distributed. Tuning η is a non-trivial work. In the previous work the sampling-based assessment of η is not suitable for clustering of small size objects (like the NCI-60 dataset in this study). In our study a leave-one-out stability assessment was used to assign a proper value of η. The ideal η is expected to render better stability when a clustering is performed. For NCI-60 dataset, we performed 37 times leave-one-out sampling for clustering of the whole data, and each time 36 compounds were selected. A series of η ranging from 0.001 to 1000 were used in the fusion model. Two parameters were used as the evaluation of the performance regarding to different η value, as listed in the following. It should be noted that other measurements can also be adopted to tune the η value in clustering, which will be generally consistent to these two measurements and will not be discussed here:

  1. 1.

    Average mean disagreement (AMD)

AMD is defined as the average value of the mean disagreement among the 37 subgroups. Given the clustering results Y 1 , 2 , ... , k n by cutting the clustering tree into k class with by cutting the clustering tree into k class with k 2 , 15 , the AMD is defined as:

AMD = 1 s Σ 1 s m e a n Y , Y ' Y Σ i = 1 n I y i y ' i
(6)

Where S is the number of subgroups, Y , Y ' are two clustering result with equal k value, y i and y ' i are class labels of element i in two results respectively, I A is the indicator function of expression A. The η value with a lower AMD in the clustering result is considered as a good parameter.

  1. 2.

    Average Dunn’s Index (ADI)

ADI is used to describe the partition quality in a clustering result. Dunn’s index is defined as the ratio of the minimal interclass distance and maximal intra-class distance. Higher Dunn’s index indicates better validity of partition. For NCI-60 dataset we calculate the average Dunn’s index regarding to a range of class number k 2 , 15 as defined in AMD. The η that obtains a high Dunn’s index with low variance could be considered as a proper estimation.

Fused similarity matrix

After estimating the sparseness controlling parameter, the alternative minimization steps were repeated on the whole dataset. Final weights vector α was calculated and the fused similarity matrix can be obtained as the convex combination of the two original similarity matrices with the weights α.

Compound-target interaction analysis

A compound-target interaction network can be constructed via target annotation in the PubChem BioAssay database. In general, one compound (identified by CID) was linked to a target protein (identified by NCBI protein ID, or GI) if this compound was tested active in the bioassay which was specified with the protein target. All the target annotation and activity information was retrieved from PubChem BioAssay database via E-Utilities tool. The interaction network was constructed and visualized by using the Cytoscape (version 2.7.0) [22], containing 37 compound nodes and 138 target nodes (in Additional file 1: Table S 2).

We proposed a quantitative method to analysis the relationship between compound similarity and their protein targets. This method is based on the concerning that compounds which have similar features, either structural or biological, tend to share common protein target. Based on such assumption, it is nature to build a connection between the quantitative similarity between compounds and the common target number within a group of compounds. And it is obvious to conclude that the common target number in a cluster derived by a clustering algorithm is an efficient measurement to measure the quality of the similarity adopted for this clustering. The larger common target number obtained in a cluster generally reveals a more reasonable similarity adopted. Followed by this strategy, in our study the compound-target interaction network was modified by taking out all the target nodes by linking two compound nodes together if they have a common protein target. The modifying process was carried out using Pajek [23]. And an average degree within a cluster was presented, which is calculated as an efficient measurement of the common target number in this cluster:

D = 1 n Σ j = 1 n D j
(7)

Where D j is the degree of node j in the graph and n is the number of nodes. The degree analysis was accomplished by using igraph package (version 0.5.1) [24] in R (version 2.12.0) [25].

Similarity fusion on large-scale CMAP dataset

In order to demonstrate the performance of the feature integration on the large-scale dataset, similarity fusion was performed on the CMap dataset. In this study, in order to generate pair-wise relationship among all the compounds, Gene Ontology (GO) fingerprint [26], which is presented in our previous study as a well-defined bioactivity representation, was adopted to combine all the expression profiles of one compound and reduce the high dimensions and noises in the microarray data. This descriptor was used to describe drug in a biological activity view. Similarly, the same structural fingerprint as used for NCI-60 data was used here to describe drug in a compound structure view.

The fusion of structural fingerprint and GO-fingerprint similarity matrices was performed following the same workflow aforementioned for NCI-60 data. And the detailed parameter optimization will not be discussed here. Considering that clustering result of large scale dataset cannot be analysed straightforwardly as the former 37-compound dataset, two typical HDAC and HSP90 inhibitors, which was used as the examples in Lamb’s work, were chosen as the queries to validate our fusion method from the perspective of virtual drug screen. For each query, the ranks of similarity searching derived by the fused similarity were compared to that with only single view, and the targets of top-ranked compounds with similarity above 0.5 to queries were also analysed for further discussion.

Results and discussions

Test results for NCI-60 dataset

Assessment of the sparseness-controlling parameter for NCI-60 data

For the stability assessment, η was chosen in a range as (0.001, 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 100, 200, 500, 1000). It should be noted that if η is smaller than 0.5, extreme large weight would be added on one of the two original similarity matrices (larger than 0.99), while η larger than 10 will generally separate the weight evenly between the two matrices, i.e. 0.5 for each. After the 37 times leave-one-out subgroup clustering, two parameters, AMD and ADI were calculated as the evaluations of the clustering quality (Figure 2).

Figure 2
figure 2

Parameter optimization. Average Mean Disagreement(AMD) and Average Dunn’s Index(ADI) with different η value. When η = 3, AMD and ADI was marked with red color.

As shown in Figure 2A, the Average Mean Disagreement reached the lowest value when η = 3 . Furthermore, the Average Dunn’s Index indicated the validity of the clustering. As shown in Figure 2B, the ADI grew gradually when η increased below 3. The decreasing variance suggested an accretive clustering quality. It should be noted that when η = 3 ADI has a sharp rise, while after that the trend of growing has become attenuated. Later calculation of weights α reveals that η lower than 3 or greater than 100 will tend to give biased weights to the two matrix, i.e. either α = 0 , 1 or α = 0.5 , 0.5 . In summary, given the best value of η in AMD, and a relative high value in ADI, it is reasonable to choose η = 3 as a proper estimation to control the sparseness.

Clustering result

A hierarchy clustering result for the 37 compounds based on fused similarity is shown in Figure 3. It should be noted that there exist several differences on the structure of the hierarchy clustering tree compared to single-view similarity clustering, as shown in Cheng’s work (The previous clustering results of single-view similarity are shown in Additional file 1: Figure S 1). Generally speaking, using correlation of bioactivity profile instead of Euclidean distance helps to find a new member in one cluster (referred to as cluster B in Cheng’s work [6]), and our fused similarity clustering result combines distinct clusters that exist in different single-view clustering separately. It should be noted that clusters with low fused similarity (with a distance above 0.7 in this study) have neither significant structural nor biological resemblance, hence no detailed analysis was presented. In the following part 3 interesting findings for NCI-60 data from the multi-view clustering were discussed compared to the clustering with only bioactivity profiles or structural fingerprints respectively.

Figure 3
figure 3

Clustering result using fused similarity. Numbers represent Pubchem Compound ID (CID). The value of distance d transformed from similarity s:d = 1 s is shown in the left. Three distinct classes was marked with a red box and named as: A [CID: 2723601, 3246652, 6351879]; B [CID: 24360, 60699, 72402, 354677, 5351879]; C [CID: 4212, 5458171].

Overall average common target number

As shown in Figure 4, the x axis is the number of classes the hierarchy tree was cut into; the y axis was the average common compound target number within one class. It is apparent that the common target number would decrease as the member within each class drops when class number increases, since the number of objects in each class decreased. The common target obtained by fused similarity, structural similarity, bioactivity profile similarity and bioactivity profile Euclidean distance were represented by red, green, blue and purple lines respectively. It is interesting to find that in general the common target number obtained by fused similarity is larger than those obtained by the other two single-view similarities, which indicates that the multi-view data representation provides a better similarity measurement and clustering validity in target-specific compound analysis compared to single-view clustering.

Figure 4
figure 4

Average common target number of the clustering result. The hierarchy clustering tree was cut into a range of classes. The common target obtained by fused similarity, structural similarity, bioactivity profile similarity and bioactivity profile Euclidean distance were represented by red, green, blue and purple lines respectively. The value of the default class number 6 was marked with a yellow line.

Figure 5
figure 5

Compound-targets network of clustering. Protein Targets are represented in rectangle shape, and the corresponding compounds are represented in eclipse shape in cluster B. Compounds in previous cluster B were marked in purple. The newly discovered class member using fused similarity was marked in blue.

Highly similar structure as complement of bioactivity profiles

In the hierarchical tree a cluster (Cluster B in Figure 3) with 6 compounds [CID: 72402, 354677, 3246719, 60699, 24360 and 97226] is distinctive in the final clustering result. Among the 6 compounds, 5 of them correspond to the Cluster B in the previous clustering achieved only with bioactivity profiles [6]. It is quite interesting that the one excluded in the single-view clustering was finally introduced into this cluster when the structural information and bioactivity profile information were considered in an integrated way. An insight into the bioactivity profiles reveals that compound [CID: 3246719] was excluded in the former study for a probable reason that its bioactivity profile shifts above the other 5 profiles (Figure 6A), but keeps the similar shape of the profile curve. By further comparing their structures, it is clearly to observe the high similarity among the 6 compounds (Figure 6B). It is possible to reason that the intrinsic similar structures of the 6 compounds results in the similar pattern of bioactivity profiles, i.e. similar chance to function in the compound-target network, and the up-shift dosage of the outlier compound above other bioactivity profiles will influence little on its functions related to specific target. However, with only bioactivity profile distance measurement, such information may be lost by ignoring structural resemblance and corresponding bioactivity correlation. It can be seen that there exist an approximately 1 order of magnitude difference between the GI50 of compound [CID:3246719] and the other 5 compounds. Therefore their bioactivity profiles will varied significantly while the correlation remains relative large.

Figure 6
figure 6

Bioactivity profile(A) and compound structure (B) of the 6 compounds in cluster B. Bioactivity profile and compound structure of the 6 compounds in cluster B were presented in Figure 6A and 6B respectively.

Further compound-target interaction analysis on the protein targets within this cluster shows that the compound CID: 3246719 shares common protein target with other members in the group, while this compound was missing in the previous compound-target network [6] as shown in Additional file 1: Figure S 1A. It can be seen that there’re three protein target linked to compound CID: 3246719 [GI: 119579178, 222080095 and 3287985] (Figure 5). All the three protein targets are shared in the former single-view clustering. It is indicated that by using multi-view similarity analysis a missing group member was discovered by introducing extra structural information. It is evident that compounds sharing similar structural features and bioactivity profiles simultaneously will give bonus to the performance in a similarity-based search. In addition, other two compounds (Cluster C, [CID: 4212 and 5458171]) can be another good examples. These two compounds are significantly similar both in structure and bioactivity profiles (Additional file 1: Figure S 3); hence they get a notably high similarity in the hierarchy tree (Figure 1).

Highly similar bioactivity profiles as complement of moderately similar structure

Another cluster composed of three compounds [CID: 2723601, 3246652 and 5351879] are noteworthy to explain in the fused hierarchy clustering tree. If we only measure the compound similarity with structural information, we can see that there are relatively less similar. However, when combined with the bioactivity information, these three compounds successfully merged into the first cluster during hierarchical clustering [6]. (Additional file 1: Figure S 2). Target interacting analysis reveals that these three compounds share a common target [GI: 4504349], indicating a potential common function in biological process. It is notable that certain fragment of the compounds, thioguanine in this example, instead of the complete structure, is essential in a binding event. Therefore when compounds that bind to a common target exhibit only relatively low overall structural similarity, it could be a good complementary to introduce the bioactivity profiles to suggest a more strong correlation with target binding potent. Such advantage of multi-view similarity assessment could be remarkable when no prior knowledge about either specific functional fragment or target is available.

Drug virtual screen based on fused similarity of CMap dataset

Firstly, trichostatin A (TSA), a typical Histone deacetylases (HDAC) inhibitor [27], was used as the query of similarity searching based on fused similarity, GO fingerprint and structural fingerprint respectively. The top 10 ranking results were listed in Table 1. It is very interesting that among all the ranking compounds, vorinostat and scriptaid, two strong HDAC inhibitors [28, 29] were successfully retrieved in the top 2 (1st and 2nd of the 1267 compounds) using fused similarity. As shown in Table 1, vorinostat and scriptaid ranks as the top two candidates in the view of GO fingerprint, which indicates their similar expression profiles. However in the view of structural fingerprint, vorinostat and scriptaid were ranked 83rd and 295th respectively. Further analysis on the weighing scheme reveals that this is reasonable since the optimized fusion method successfully give a larger weight on the information-rich GO fingerprint.

Table 1 Similarity search based on fused similarity: HDAC inhibitors

More interesting results were discovered on the similarity searching of geldanamycin, an HSP90 inhibitor [30]. (Table 2) It can be seen that tanespimycin and alvespimycin, both the derivants of geldanamycin, which are also typical HSP90 inhibitors [31, 32], were ranked at 1st and 2 rd place using fused similarity. In addition, monorden, which is another common HSP90 inhibitor [33] was ranked in top 10 (rank 8th). However, when using only GO fingerprint and Structure fingerprint, monorden was ranked 10th and 123rd respectively. More interestingly, 15-delta prostaglandin J2 (15d-PGJ2), suggested to exert anti-inflammatory effects in vivo[34], ranked 3rd in the fusion similarity searching. Further literature research on its target indicates that HSP90 is a target for modification by 15d-PGJ2 in renal mesangial cells. [35] This result shows that by assigning balanced weights to the two views, our fusion method successfully picked out this newly discovered HSP90 inhibitor, thus demonstrated the fused similarity provides an effective quantitative assessment of drug-target relationship.

Table 2 Similarity search based on fused similarity: HSP90 inhibitors

Conclusions

A multi-view clustering method was introduced to discover a more robust correlation between fused multi-view similarity and compound-target interacting pattern. By using a similarity-based optimization and fusion model, a hierarchy clustering integrated with both structural and bioactivity profile information was presented on the NCI-60 dataset. It is interesting that comparing to single view analysis, the overall common target number within fused classes has been promoted by integrating information from two views, which indicated a more robust and efficient representation of compound related to specific target. Analysis of compound-target interaction network shows that fusion of data source from different views enhances similar compound discovery, leading to a more comprehensible assessment of target-binding potent. Further analysis in certain classes with high fused similarity shows that the mutual complement of the two views can lead to the discovery of missing similar compound with only one view. A further large-scale similarity searching on the CMap data based on the fused similarity also obtained a better ranking results compared to that of single-view for two inhibitors as queries, thus indicate the potential use of our quantitative similarity fusion in virtual drug screen. In summary, our findings are interesting for the following reasons: Firstly, both the bioactivity profiles and structural fingerprint lack to be a completely direct indicator of interaction, i.e. only partial features instead of overall characterization from either view are essential in a binding event. Hence by integrating potentially correlating features from both views to maximize the utility of available data source, a robust similarity assessment could be achieved without prior knowledge about the detail relationship between target-binding rules and compound features. Secondly, the fusion method in this study provides an extendable framework of integrating multi-view data. Fusion process is applicable to various situations when more than two data sources are available. A comprehensive assessment of the similarity can be achieved in virtual drug screening when various potential pharmacological properties of compounds are integrated.

References

  1. Eckert H, Bojorath J: Molecular similarity analysis in virtual screening: foundations, limitations and novel approaches. Drug Discov Today 2007, 12(5–6):225–233.

    Article  CAS  PubMed  Google Scholar 

  2. Cheng T, Li Q, Wang Y, Bryant SH: Identifying compound-target associations by combining bioactivity profile similarity search and public databases mining. Journal of Chemical Information and Modeling 2011, 51(9):2440–2448. 10.1021/ci200192v

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Bryant SH, Wang YL, Xiao JW, Suzek TO, Zhang J, Wang JY: PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res 2009, 37: W623-W633. 10.1093/nar/gkp456

    Article  PubMed Central  PubMed  Google Scholar 

  4. Barrett T, Troup DB, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, et al.: NCBI GEO: archive for functional genomics data sets-10 years on. Nucleic Acids Res 2011, 39: D1005-D1010. 10.1093/nar/gkq1184

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, Gautam B, Hassanali M: DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res 2008, 36: D901-D906.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Cheng TJ, Wang YL, Bryant SH: Investigating the correlations among the chemical structures, bioactivity profiles and molecular targets of small molecules. Bioinformatics 2010, 26(22):2881–2888. 10.1093/bioinformatics/btq550

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Sheng Z, Kang H, Dai TL, Liu Q, Zhu RX: Complementary study of structure features and gene profile features for chemical compounds. Acta Chim Sinica 2010, 68(23):2395–2400.

    CAS  Google Scholar 

  8. Glover CJ, Rabow AA, Isgor YG, Shoemaker RH, Covell DG: Data mining of NCI's anticancer screening database reveals mitochondrial complex I inhibitors cytotoxic to leukemia cell lines. Biochem Pharmacol 2007, 73(3):331–340. 10.1016/j.bcp.2006.10.005

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Hartley R, Zisserman A: Multiple view geometry in computer vision. Cambridge University Press, The Edinburgh Building, Cambridge, CB2 2RU, UK; 2003.

    Google Scholar 

  10. Blum A, Mitchell T: Combining labeled and unlabeled data with co-training. In Proceedings of the eleventh annual conference on Computational learning theory. Madison, Wisconsin, United States: ACM; 1998: 92–100 1998: 92-100

  11. Ghani R: Combining labeled and unlabeled data for multiclass text categorization. In Proceedings of the Nineteenth International Conference on Machine Learning. Morgan Kaufmann Publishers Inc, San Francisco, CA, USA; 2002: 187–194 2002: 187-194

  12. Brefeld U, Scheffer T: Co-EM support vector learning. In: Proceedings of the twenty-first international conference on Machine learning Banff, Alberta, Canada: ACM; 2004:16. Banff, Alberta, Canada: ACM; 2004:16.

  13. Zhou Z, Li M: Semi-supervised regression with co-training. In: In Proc of Int’l Joint Conf on Artificial Intelligence, Sydney, Australia 2005.

    Google Scholar 

  14. Bruno E, Marchand-Maillet S: Multiview Clustering: a late fusion approach using latent models. Proceedings 32nd Annual International Acm Sigir Conference on Research and Development in Information Retrieval 2009. Boston, MA, USA: 736–737.870. Boston, MA, USA: 736-737.870.

    Google Scholar 

  15. Lange T, Buhman J: Fusion of similarity data in clustering. In Proceeding of Advances in Neural Information Processing Systems 18. MIT Press, Cambridge, MA, USA; 2005.

    Google Scholar 

  16. Long B, Yu PS, Zhang ZM: A general model for multiple view unsupervised learning. In In SDM. Atlanta, Georgia, USA; 2008.

    Google Scholar 

  17. Cowan I, Tesauro G, Virginia R, De Sa: Learning classification with unlabeled data. Proc NIPS93 Neural Information Processing Systems, Morgan Kaufmann Publishers, Denver, CO, USA, 1993, Pages: 112–119.

  18. Bickel S, Scheffer T: Multi-view clustering. In Proceeding of IEEE data mining Conference. Brighton, UK, 2004; 2004.

    Google Scholar 

  19. Chaudhuri K, Kakade SM, Livescu K, Sridharan K: Multi-view clustering via canonical correlation analysis. In Proceedings of the 26th Annual International Conference on Machine Learning. Montreal, Quebec, Canada: ACM; 2009:129–136.

    Google Scholar 

  20. Greene D, Cunningham P: A matrix factorization approach for integrating multiple data views. In Proceedings of the European Conference on Machine Learning and Knowledge Discovery in Databases: Part I. Springer, Bled, Slovenia; 2009:423–438 2009:423-438

  21. Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, et al.: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science 2006, 313(5795):1929–1935. 10.1126/science.1132939

    Article  CAS  PubMed  Google Scholar 

  22. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498–2504. 10.1101/gr.1239303

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Batagelj V, Mrvar A: Pajek - Analysis and visualization of large networks. In: Graph Drawing Software, Mathematics and Visualization.Edited by Jünger, Michael and Mutzel, Petra.Springer Berlin Heidelberg 2004:77–103 2004:77-103

  24. Csardi G, Nepusz T: The igraph software package for complex network research. InterJournal 2006, 1695. Volume: Complex Systems Volume: Complex Systems

    Google Scholar 

  25. Team RDC: R: A language and environment for statistical computing. Austria: R Foundation for Statistical Computing, Vienna; 2009.

    Google Scholar 

  26. Sheng Z, Huang Q, Kang H, Liu Q, Cao ZW, Zhu RX: A new fingerprint of chemical compounds and its application to drugs virtual screening. Acta Chim Sinica 2011, 69(16):1845–1850.

    CAS  Google Scholar 

  27. Vanhaecke T, Papeleu P, Elaut G, Rogiers V: Trichostatin A - like hydroxamate histone deacetylase inhibitors as therapeutic agents: toxicological point of view. Curr Med Chem 2004, 11(12):1629–1643. 10.2174/0929867043365099

    Article  CAS  PubMed  Google Scholar 

  28. Ree AH, Dueland S, Folkvord S, Hole KH, Seierstad T, Johansen M, Abrahamsen TW, Flatmark K: Vorinostat, a histone deacetylase inhibitor, combined with pelvic palliative radiotherapy for gastrointestinal carcinoma: the Pelvic Radiation and Vorinostat (PRAVO) phase 1 study. Lancet Oncol 2010, 11(5):459–464. 10.1016/S1470-2045(10)70058-9

    Article  CAS  PubMed  Google Scholar 

  29. Keen JC, Yan L, Mack KM, Pettit C, Smith D, Sharma D, Davidson NE: A novel histone deacetylase inhibitor, Scriptaid, enhances expression of functional estrogen receptor alpha (ER) in ER negative human breast cancer cells in combination with 5-aza 2 '-deoxycytidine. Breast Cancer Res Tr 2003, 81(3):177–186. 10.1023/A:1026146524737

    Article  CAS  Google Scholar 

  30. Li YH, Tao PZ, Liu YZ, Jiang JD: Geldanamycin, a ligand of heat shock protein 90, inhibits the replication of herpes simplex virus type 1 in vitro. Antimicrob Agents Ch 2004, 48(3):867–872. 10.1128/AAC.48.3.867-872.2004

    Article  CAS  Google Scholar 

  31. Yao JQ, Liu QH, Chen X, Yang Q, Xu ZY, Hu F, Wang L, Li JM: Hsp90 inhibitor 17-allylamino-17-demethoxygeldanamycin inhibits the proliferation of ARPE-19 cells. J Biomed Sci 2010, 17:30.

    Google Scholar 

  32. Ramanathan RK, Egorin MJ, Erlichman C, Remick SC, Ramalingam SS, Naret C, Holleran JL, TenEyck CJ, Ivy SP, Belani CP: Phase I Pharmacokinetic and Pharmacodynamic Study of 17-dimethylaminoethylamino-17-demethoxygeldanamycin, an inhibitor of heat-shock protein 90, in patients with advanced solid tumors. J Clin Oncol 2010, 28(9):1520–1526. 10.1200/JCO.2009.25.0415

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Wang SH, Xu YQ, Maine EA, Wijeratne EMK, Espinosa-Artiles P, Gunatilaka AAL, Molnar I: Functional characterization of the biosynthesis of Radicicol, an Hsp90 inhibitor resorcylic acid lactone from chaetomium chiversii. Chem Biol 2008, 15(12):1328–1338. 10.1016/j.chembiol.2008.10.006

    Article  CAS  PubMed  Google Scholar 

  34. Olmos G, Arenas M, Bienes R, Calzada M, Aragonés J, Garcia-Bermejo M, Landazuri M, Lucio-Cazaña J: 15-Deoxy-Δ<sup>12,14</sup>-prostaglandin-J<sub>2</sub>reveals a new pVHL-independent, lysosomal-dependent mechanism of HIF-1α degradation. Cellular and Molecular Life Sciences 2009, 66(13):2167–2180. 10.1007/s00018-009-0039-x

    Article  CAS  PubMed  Google Scholar 

  35. Stamatakis K, Sanchez-Gomez FJ, Perez-Sala D: Identification of novel protein targets for modification by 15-deoxy-Delta12,14-prostaglandin J2 in mesangial cells reveals multiple interactions with the cytoskeleton. J Am Soc Nephrol 2006, 17(1):89–98.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Funding: This work was supported in part by Project White Magnolia Funding, Shanghai (Grant No. 2010B127), Shanghai Pujiang Talents Funding (Grant No.11PJ1407400), Tongji Excellent Young Scientist Funding (Grant No. 2000219052), National Natural Science Foundation of China (Grant No. 30976611, Grant No.31100956 and Grant No. 61173117), Research Fund for the Doctoral Program of Higher Education of China (Grant No.20100072120050, 20110072120048) and 973 National Key Basic Research Program of China (Grant No. 2010CB833601).

We would like to express our special thanks for Tiejun Cheng and other authors of [6] in NIH, U.S. They kindly helped us for the NCI-60 dataset retrieving.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Qi Liu or Zhiwei Cao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ZC and QL conceived and participated in the design of the study. TX and RZ performed the data collection and processing. TX and QL performed the whole experiments and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2011_5507_MOESM1_ESM.doc

Additional file 1:Additional tables and figures were saved in the word file entitled “supplementary.doc”, containing the 37 compounds CID list, target information and other clustering results for the NCI-60 dataset, etc. Additional tables and figures were saved in the word file entitled “supplementary.doc”, containing the 37 compounds CID list, target information and other clustering results for the NCI-60 dataset, etc.(DOC 978 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 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.

Reprints and permissions

About this article

Cite this article

Xu, T., Zhu, R., Liu, Q. et al. Quantitatively integrating molecular structure and bioactivity profile evidence into drug-target relationship analysis. BMC Bioinformatics 13, 75 (2012). https://doi.org/10.1186/1471-2105-13-75

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-13-75

Keywords