Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(7): e69513.
Published online 2013 July 23. doi:  10.1371/journal.pone.0069513
PMCID: PMC3720659

Integrated Analysis of Drug-Induced Gene Expression Profiles Predicts Novel hERG Inhibitors

James P Brody, Editor


Growing evidence suggests that drugs interact with diverse molecular targets mediating both therapeutic and toxic effects. Prediction of these complex interactions from chemical structures alone remains challenging, as compounds with different structures may possess similar toxicity profiles. In contrast, predictions based on systems-level measurements of drug effect may reveal pharmacologic similarities not evident from structure or known therapeutic indications. Here we utilized drug-induced transcriptional responses in the Connectivity Map (CMap) to discover such similarities among diverse antagonists of the human ether-à-go-go related (hERG) potassium channel, a common target of promiscuous inhibition by small molecules. Analysis of transcriptional profiles generated in three independent cell lines revealed clusters enriched for hERG inhibitors annotated using a database of experimental measurements (hERGcentral) and clinical indications. As a validation, we experimentally identified novel hERG inhibitors among the unannotated drugs in these enriched clusters, suggesting transcriptional responses may serve as predictive surrogates of cardiotoxicity complementing existing functional assays.


While the single-target approach to drug discovery seeks “silver bullets” that selectively modulate disease-related proteins, recent work has emphasized the often promiscuous interactions of both marketed and candidate therapeutics [13]. The positive impact of such polypharmacology includes the potential to discover novel clinical uses for previously approved medications [46]. However, it also suggests that drugs may share similar and undesirable side effects despite unrelated chemical structures or primary mechanisms-of-action (MOA). While existing quantitative structure activity relationship (QSAR) methods have leveraged structural features of small molecules to predict toxicity, the difficulty of applying such techniques to chemicals that vary substantially from the model inputs has been described, particularly in cases where toxicity is linked to the metabolic by-products of a compound [7,8]. Thus alternative descriptors, such as measurements of drug effects that probe the complex physiology of the cell, may potentially reveal commonalities aiding the prediction of toxicity independent of chemical structure as represented, for example, by conventional chemical fingerprints. Here, we explored similarities in drug-induced transcriptional effects using the Connectivity Map (CMap), a collection of Affymetrix™ microarray profiles generated by treating three independent lineages of cancer cell lines with small molecule drugs [9]. In previous applications, analysis of the CMap has associated transcriptional signatures to known MOAs or disease states, allowing the discovery of novel modulators of autophagy, small cell lung cancer proliferation, and inflammatory bowel disease [5,6,10]. Similarly, computational studies have identified correlations between known drug side effects and transcriptional responses in the CMap [11,12]. Thus, we hypothesized that this data might also be used to predict and verify novel toxicities, which we demonstrate by integrating the CMap with experimentally measured inhibition data for the human ether-à-go-go related (hERG) potassium channel and literature annotations to identify novel antagonists of this important anti-target of many drugs.

Promiscuous inhibition of the hERG channel by therapeutically and structurally diverse drugs prolongs the QT interval quantified by surface electrocardiogram (ECG) [13]. This phenomenon, known as drug-induced Long QT (LQT) syndrome, is a risk factor for sudden cardiac death [13]. To date, the lack of universal chemical patterns and diversity of primary clinical targets among known hERG inhibitors have impeded effective risk assessment of this side effect using computational methods, and experimental evaluation using the “gold standard” of electrophysiology remains an important step in therapeutic development. Such electrophysiological recordings, utilizing recombinantly expressed hERG channels [1416] as well as patient-derived cardiomyocytes [17,18], have afforded valuable experimental opportunities to study the potential LQT side effects of small molecules. More recently, the development of high-throughput electrophysiology platforms has facilitated systematic evaluation of hERG inhibition in large compound collections [19,20]. Concurrently, potential global physiological readouts for channel function are suggested by behavioral assays in model organisms such as C. elegans and D. rerio [21,22], as well as reports linking channel activity to tumor migration and volume [23,24], indicating these phenomena may conceivably be used as ways to probe hERG liability. Computationally, hERG inhibition has also been correlated with the proximity of a drug’s therapeutic target to hERG in a protein–protein interaction network [25].

Our present analysis integrates earlier results in which we have independently profiled over 300,000 compounds (including approximately half of the CMap compounds) in the NIH Molecular Library Small Molecule Repository (MLSMR) for their ability to inhibit hERG current in a high-throughput electrophysiological assay [26]. Combining our database with additional publicly available annotations for LQT side effect allowed us to identify clusters of drugs with similar expression profiles in the CMap enriched for channel inhibitors. Drugs of unknown hERG liability within these clusters, through the principle of ‘guilt by association’, were then experimentally validated using an electrophysiology assay. These results advance the hypothesis that structurally diverse hERG inhibitors mediate similar physiological effects revealed by transcriptional response profiles, even in cell lines not derived from a cardiac lineage and potentially independent of hERG expression. Thus, gene expression signatures may serve as a proxy measurement correlated with reduction of hERG current, and find practical application as a high-throughput platform to complement existing electrophysiological assays. More generally, these analyses suggest that side effect profiles as well as primary MOAs may be predictively correlated with similarities in drug-induced gene expression responses independent of chemical structure.


Microarray normalization and analysis

Our analysis pipeline is outlined in Figure 1. The Connectivity Map (CMap) ( [9] is a collection of Affymetrix™ microarray gene expression profiles representing the responses of three cancer cell lines (breast cancer: MCF7, prostate cancer: PC3, and leukemia: HL60) to small molecule treatments in comparison to dimethyl sulfoxide (DMSO) treated samples used as vehicle-treatment controls for these studies because most drugs are dissolved using DMSO as a solvent. The full dataset currently contains one or more measurements for 1,309 unique substances at varying concentrations. To begin our investigation, we normalized the downloaded CMap data and corrected for batch effects (similarity among arrays correlated with the experimental group in which they were processed rather than biological annotations) by centering the mean expression value per probeset in each batch at zero, as previously proposed [27]. Thus, the resulting profiles represent changes in gene expression dependent on drug effect rather than experimental batch. Following this correction, we also noted a sub-population of arrays representing replicate treatments (the same drug tested against the same cell line at the same concentration) that have approximately zero Pearson correlation (as judged by the correlation of log2 fold changes versus vehicle control). By using thresholds for minimum log2 fold gene expression change among probesets on an array, it is possible to split this population into two groups of correlated and uncorrelated replicates (Figure S1). We propose that these uncorrelated replicates may represent “transcriptionally silent” drugs where the observed response is random signal variation and thus not preserved between repeats, and thus removed them from our sample prior to downstream analysis. Finally, we used the batch-corrected DMSO-treated samples to identify probesets in each drug-induced expression profile whose fold change lies outside the range of variation exhibited by the same probeset in DMSO-treated control samples measured in the same cell background. Table 1 summarizes the number of arrays filtered and unique drug instances filtered during each stage of this analysis, as well as the number of experimental and annotated hERG inhibitors present in each of these samples. The final sample of 673 unique drugs contains 1,033 drug-cell line combinations (‘Merge Replicates’ row in Table 1), as some drugs were profiled in more than one cell background, with 62 experimentally and 57 clinically annotated as hERG inhibitors.

Figure 1
Pipeline for construction and analysis of drug transcriptional response network.
Table 1
Microarray data processing statistics.

Enrichment of structurally diverse hERG inhibitors through transcriptional response similarity

Following pre-processing, we clustered the resulting collection of drug-induced gene expression profiles using affinity propagation [28], an unsupervised learning algorithm that automatically identifies the optimal number of clusters in a dataset using an input of all pairwise similarities (here, the pairwise Pearson correlations between expression profiles). Each cluster generated by this procedure contains an “exemplar”, a single member that best characterizes the pattern shared by the members of the group. To identify higher-level relationships between individual clusters, we further aggregated the data by recursively re-clustering these exemplars to attain a global view of the number of characteristic patterns of drug-induced gene expression changes in this collection. We integrated the gene expression measurements with annotations for hERG inhibition derived from two sources: a previously described dataset of electrophysiology measurements of hERG inhibition ( [26], and lists of drugs that have been clinically linked to LQT side effects ( and Drugs with records in hERGcentral were annotated as inhibitors if they decreased hERG current by 50% or more at 10 µM concentration, representing an IC50 value of approximately 10 µM or less. We selected this threshold as the dataset in hERGCentral contains inhibition measurements at 1 µM and 10 µM, and literature data is frequently annotated with potency (IC50) values, making these two concentrations the most convenient thresholds. We chose a less conservative (10 µM) threshold, as this value correctly identifies 40/53 (76%) of torsades de pointes (TdP)-risk drugs described in a literature survey [29]. The agreement between the hERGcentral measurements (continuous values) and existing LQT drug annotations (binary classifications) is demonstrated by a Wilcoxon rank-sum test comparing the median experimentally determined hERG inhibition values from hERGcentral for drugs with or without previous annotation for LQT side effects (p-value 5.7x10-12, Figure S2). Complete experimental and clinical annotations based on the criteria described above are given in Table S1.

The clusters generated from the drug-induced gene expression profiles derived from the breast cancer cell line MCF7 are displayed in a network diagram in Figure 2A, where nodes represent individual drugs and edge weights represent the magnitude of similarity (Pearson correlation coefficient) between a given drug’s expression profile and the exemplar of its cluster. We found that 2 of the 31 resulting clusters contained an enriched fraction of hERG inhibitors compared to randomized clusters of the same size (using a false discovery rate threshold of 0.2). Even after correction for experimental batch effects described above, clustering of all drug-induced expression profiles demonstrates assortment by cell background, suggesting the existence of cell line-specific drug effects (Figure 2B). However, as the exemplars of these sub-clusters are hierarchically merged, connections emerge between hERG inhibitor-enriched clusters derived from all three different cells lines, indicating the presence of general as well as cell background-specific responses (Figure 2B). This interpretation is supported by a scatterplot of average expression changes in these five hERG inhibitor-enriched clusters in the three cell lines, which indicates that some differentially expressed (DE) genes are shared (Figure S3A), as well as Venn diagrams indicating the overlap of DE genes between these sets (Figure S3, B & C). The clustering results for each of the networks in Figure 2 are given in Table S2.

Figure 2
Network analysis of drug-induced gene expression profiles.

Further, we note that many of the exemplars of these hERG inhibitor-enriched clusters are preserved between cell lines (Figure 2B). The enrichment of hERG blockers with previous experimental or clinical annotation among the five enriched clusters identified in Figure 2B were further quantified with the hypergeometric test, with resulting prediction statistics summarized in Table 2. Quantification of the resulting predictive power in Table 3 suggests an overall accuracy of 82% based on drugs for which experimental or clinical annotation is available (e.g., excluding the grey ‘untested’ drugs in Figure 2B), which is consistent with a test of ‘good’ quality based on previously published metrics [30]. While the pairwise gene expression profile correlations (Pearson coefficients) within the hERG inhibitor-enriched clusters are significantly higher than the correlations between enriched cluster drugs and non-enriched cluster drugs (medians of 0.14 and 0.01, respectively, Wilcoxon rank-sum test p-value < 2.23x10-308) (Figure 3, A & B), the corresponding distributions of pairwise chemical similarities (Tanimoto coefficients) are statistically different (as judged by a Wilcoxon rank-sum test comparing the inter-cluster chemical similarity of the hERG inhibitor-enriched clusters versus their similarity to drugs in other clusters, p-value 1.1510x10-288) yet possess approximately equal medians (0.12 and 0.10, respectively) (Figure 3C). Thus, this analysis highlights correlations in drug-induced gene expression profiles that are not evident from chemical similarity alone. Intriguingly, we also noted that the MCF7 Astemizole-exemplar cluster includes Miconazole and Mefloquine, drugs which have been previously shown to inhibit hERG channels recombinantly expressed in cell lines [31,32], but did not appear in our lists of clinically annotated LQT-causing drugs and were inactive in our high-throughput electrophysiology assay. This suggests that this dataset may contain false negatives which nevertheless cluster with other known inhibitors based on similarity in transcriptional responses. Conversely, compensatory block of other ionic currents in addition to hERG may normalize hERG effects by these drugs, leading to no observable clinical phenotype [33,34]. Examination of the relationship between a drug’s hERG inhibition and maximal Pearson correlation to any of the five hERG blocker enriched cluster exemplars in Figure 2B demonstrates a modest linear correlation which is statistically significant compared to randomized data (Figure S4, A & B). Further, a greater fraction of drugs with high hERG inhibition are present in the enriched clusters of Figure 2B than those with low hERG inhibition (Figure S4C). Evaluation of enriched gene ontology (GO) annotations among genes up and down-regulated in the five hERG inhibitor-enriched clusters indicated positive effects on cholesterol biosynthesis (GO:0006695), isoprenoid biosynthesis (GO:0008299), and the unfolded protein response (GO:0030968), and negative effects on cell cycle checkpoint (GO:0000075), Sphase of mitotic cell cycle (GO:0000084), and DNA replication (GO:0006260). The physiological correlation between hERG block and these processes remains to be investigated though intriguingly, previous reports have linked hERG channel activities to a variety of biological processes in addition to cardiac function [23,24,35,36]. Functional enrichment results for all clusters are given in Table S3.

Table 2
Statistical enrichment of hERG inhibitors in transcriptionally determined clusters.
Table 3
Prediction statistics for transcriptional-signature based hERG inhibitor enrichment.
Figure 3
Expression and structural similarity of hERG inhibitor-enriched clusters.

Experimental validation of predicted hERG inhibitors

To determine whether our analysis could predict novel hERG ligands among the compounds in the inhibitor-enriched clusters, we examined drugs without existing experimental or clinical annotation from the databases used in our analysis in these groups. Figure 4A demonstrates that the Astemizole cluster from the MCF7-tested drug set is the overall center of all hERG inhibitor enriched clusters, representing the most characteristic pattern for these five groups. The structures of the six ‘missing-data’ drugs in the MCF7 Astemizole cluster (Figure 4A) display limited structural similarities, and differences in functional moieties (such as the three chloro-groups of Sulconazole) that could reasonably alter their surface polarities, along with variation in linker group composition (with Fendiline, Sulconazole, and Cloperastine bearing, respectively, a nitrogen, sulfur, and oxygen atom along their carbon backbones). Four of these compounds (Fendiline, Cloperastine, Ethopropazine, and Sulconazole) were untested in our previous electrophysiology data ( [26] and lacked annotation for drug-induced LQT syndrome ( and Though not included in the LQT drug lists we utilized, we found previous literature associating Fendilene with drug-induced LQT [37], but no evidence of its direct inhibition of hERG current. Additionally, Sulconazole has been previously annotated in the bioactivity records in the ChEMBL database as inactive in a hERG binding assay using Astemizole displacement as a functional readout [38]. The remaining compound in the Astemizole cluster, Clomiphene, was not in the databases used to annotate Figure 2B but has been previously shown to inhibit the current of recombinantly expressed hERG channels with an IC50 of 0.18 µM [39]. We could locate no previous data describing Hexetidine, Cloperastine, or Ethopropazine effect on hERG or association with LQT. We therefore evaluated several of these unannotated compounds experimentally for inhibition of hERG current using automated electrophysiology recordings of a Chinese hamster ovary (CHO) stable cell line, a popular expression system due to its low background current from endogenous potassium channels compared to human epithelial kidney 293 (HEK293) cells [40,41]. With the exception of Hexetidine, which appears to disrupt the membrane seal during recording, the remaining four were successfully tested. Figure 4B shows the measured dose–response curves, for which we calculated IC50 values of 18.6 ± 0.8 µM (Ethopropazine, n = 4), 0.36 ± 0.05 µM (Cloperastine, n = 4), 2.6 ± 0.1 µM (Fendiline, n = 4), and 6.1 ± 0.9 µM (Sulconazole, n = 4). Because the Ionworks platform used in our assessment tends to underestimate the potency of compounds [42,43], the effective concentrations could be even lower.

Figure 4
Experimental validation of novel hERG inhibitors.

To assess the statistical significance of the results, we simulated 1,000 random sets of four compounds using the distribution of previously recorded experimental data for 10 µM inhibition of hERG current for the subset of MCF7-tested drugs (those present in the network of Figure 2A) in our high-throughput electrophysiology assay (Figure S5A), finding that the average percent inhibition of the tested compounds (67%, or a hERG activity of 33%) was greater than any random group (an empirical p-value of <0.001) (Figure S5B). Therefore, our clustering analysis significantly enriches for hERG inhibition among previously unannotated compounds.


In this study we identified commonalities in the transcriptional responses of structurally diverse hERG inhibitors, suggesting microarrays as a novel proxy measurement correlated with conduction of potassium currents by hERG and hence liability of channel block. Perhaps as remarkably, the observed change in gene expression is not necessarily a direct result of hERG inhibition. While hERG inhibitor-enriched clusters were observed for profiles generated in all three cell lines, functional evidence for hERG expression has been reported only for MCF7 and HL60 [23,24], indicating that channel expression may not be strictly required for the observed pattern. Indeed, there are also reports that channels including hERG have activities in addition to ion conductance [4446], indicating that independent molecular targets might converge on common signaling pathways or processes also modulated by hERG and leading to the observed correlation. For example, there is a tendency that hERG inhibitors or LQT-causing drugs are also antagonists of the multidrug resistance transporter (MDR) [47]. Alternatively, hERG may be co-expressed with other channels correlated with oncogenesis which possess similar pharmacological profiles, such as hEAG [48], thus confounding causal inference of the relationship between hERG activity and gene expression response.

We also note that the presence of inhibitors in the CMap outside the enriched clusters highlighted in our analysis indicates that this “hERG signature” is not necessarily “dominant” over other expression pattern(s), implying that other such patterns might perturb or mask the signature from being identified for some compounds. In this interpretation, the subset of clustered inhibitors highlighted in our analysis represent drugs for which this signature is dominant over or of equal strength with other expression responses of the compounds. Additionally, we note that a large portion of the compounds in this dataset exhibit silent or weak transcriptional response which prevents profiling for hERG inhibition using the proposed approach. As the signatures in CMap are uniformly generated at a 6 hour time point, it is possible that some compounds may display chronic transcriptional effects at a later time point, and thus be profiled by modifications in the original screening protocol. Indeed, previous studies of time course data from drug-induced gene expression responses have indicated that distinct expression patterns may be detected at different time points [4951], even for frequent measurements such as 3, 6, and 9 hours. We thus hypothesize that some of the ‘silent compounds’ in our study might have detectable signatures at later time points, while the hERG inhibitors outside of enriched clusters may exhibit a dominant ‘hERG signature’ at earlier time points. Taken together, these results suggest that improved sensitivity for this assay might be achieved by using time course instead of single point expression data. Additionally, we note that the sensitivity of our assay may be effected by our choice of a 10 µM IC50 threshold. While this threshold has been used in previous hERG predictive models [52], previous studies have also reported greater accuracy with a lower threshold (e.g., IC50 40 µM) [53]. Thus it may also be possible to improve the sensitivity of our method using inhibition measurements derived from higher drug concentrations.

Research has also suggested that the duration of action potentials at 90% repolarization (APD90), a correlate of clinical LQT which is elongated by hERG inhibition, may be dependent upon multi-channel drugs effects [33,34], and thus the ability of our approach to forecast clinical endpoints may be aided by future integration of high-throughput recording data for other cardiac channels such as Nav1.5. Furthermore, despite the current lack of causal evidence linking the gene-expression profiles of the clustered hERG inhibitors in the CMap with functional modulation of the channel, this analysis does suggest an intriguing possibility that some hERG inhibitors induce a downstream signaling cascade as a consequence of current reduction that is visible as a global change in gene expression. Alternatively, these observations may indicate signaling pathways downstream of potassium channels that are not directly related to their role in conduction [44]. A selection of these hypotheses is diagrammed in Figure 5. Certainly, profiling selective inhibitors of hERG such as E4031 which are not present in the CMap might help clarify these hypotheses, though the large number of transcriptionally silent compounds in the dataset suggests these selective inhibitors may not exhibit a detectable signature at 6 hour time points.

Figure 5
Mechanistic hypotheses for hERG-inhibition correlated gene expression signatures.

From a practical standpoint, the observed similarity of microarray profiles among electrophysiologically confirmed but structurally diverse inhibitors argues for the potential of using such a surrogate as an informative descriptor for hERG liability complementing existing electrophysiological assays. The utility of such a platform is suggested by compounds such as the antidepressant Amoxapine [54] which display slow on-rates beyond the temporal resolution of high throughput electrophysiological systems (which is typically less than 5 minutes), and thus appear as ‘negatives’ in our acute experimental electrophysiology data. In contrast, the microarray data utilized in this study were generated 6 hours following drug treatment, suggesting gene expression measurements may offer complementary temporal resolution not readily accessible by automated electrophysiology data, allowing high-throughput assessment of hERG inhibition in compounds with slow on-rates which have previously required manual patch clamp recordings to resolve. Furthermore, transcriptional signatures may identify false negatives from other assays, such as Sulconazole, which was labeled as inactive in ChemblDB from binding data. Because binding experiments often utilize displacement of a known ligand, they will not identify compounds binding at alternative sites(s) of action.

Gene expression measurements have additional attractive properties compared to other high-throughput technologies. Because microarray profiles represent an integrated output of multiple signaling pathways in the cell, they are potentially more sensitive than biochemical or cellular assays which are commonly designed to test one or a limited number of physiological parameters. Such expression profiles are also certainly more general in terms of measuring diverse signaling pathways and integrated biological events. Thus, assessment of hERG liability may be effectively evaluated in parallel with other endpoints of biological interest, such as inflammatory signaling, oxidative damage response, or metabolic perturbations. Additionally, the fact that our signature utilizes measurements in cancer cells derived from different tissues of origin suggests the attractive possibility of assaying the effects of hERG activity in these oncogenesis models, as previous research has linked hERG expression to tumor migration and cell volume [24,55]. Admittedly, cells with cardiac lineage may be equally or more informative. Indeed, patient-derived induced pluripotent stem cell (iPSC) models of cardiac disease have proven to be attractive disease models in electrophysiology studies [17,18,56], with additional evidence suggesting the potential for cardiac-specific transcriptional activity that may find utility in genomic drug-activity profiles [5759]. Combined with cost savings generated by custom arrays that measure only the subset of differentially expressed genes correlated with hERG risk, these aspects suggest the potential for a novel genetic platform to assess ion channel activity.

More generally, our analysis contributes to growing evidence that systems-level measurements of drug effect reveal connections and similarities often invisible from the perspective of single molecular descriptors or activity measurements [6,9,18,60,61]. These links suggest not only the possibility of mining such connections for predictive purposes, but also that the full pharmacological complexity of even long-standing medications may not yet be appreciated. Integrated analyses are thus poised to illuminate these patterns and suggest possibly novel indications or, as in our study, liabilities of existing drugs.


Gene expression and drug property/activity data

All raw data comprising the Connectivity Map (CMap) build 02 were downloaded as CEL files from the Broad Institute ( We annotated compounds for percentage of hERG inhibition using data from a previously described database of high-throughput electrophysiology measurements [26], and obtained LQT-risk data from online references at and Simplified molecular input line entry system (SMILES) strings representing the chemical structures of all compounds were downloaded from PubChem and ChemBank. Computational filtering of salts, standardization of charge and coordinates, and calculation of functional circular fingerprints (FCFP_6) were performed with Pipeline Pilot Student Edition v 6.1 (Scitegic).

Microarray pre-processing

Our analysis of the Connectivity Map Build 02 data consisted of four steps:

Platform selection and probeset normalization

From the 7,056 cell files in the CMap build 02, we selected 6,029 files generated from HT-HG-U133A arrays, consisting of 5,242 drug treatment instances and 967 DMSO-treated vehicle controls tested on three tumor cell lines of human origin. This selection is performed because unlike the non-parametric processing used in the original analysis [9], the probeset normalization algorithms we employed require a homogenous platform. All selected CEL files were probeset-normalized using GC robust multi-array average (GCRMA) background correction as implemented in the aroma. affymetrix R package [62,63] in R 2.14.2 [64].

Batch correction

Following pre-processing, we sought to remove correlations between arrays due to experimental batch (date of data acquisition) rather than biological similarity by mean-centering probesets across all drugs in each batch, following a previously described pipeline [27]. Since this correction assumes that on average a probeset should not be differentially expressed among an experimental batch of otherwise unrelated drugs, we retained only batches with sufficient numbers (>25) for this assumption to reasonably hold. To provide the most consistent comparison between replicates of the same drug in the same cell background, we retained only arrays representing the single concentration with the most examples for a given drug (the mode). However, we also note that the variation of test concentrations across CMap instances for a given drug is not often large, and thus this source of variability is likely minor. Applying these two criteria left 4,754 drug treatment and 700 control instances.

Filtering Transcriptionally Inactive Drugs

It has previously been reported that some drugs in the CMap are “inactive”, as judged by lack of correlation among replicates (arrays representing transcriptional response to the same drug in the same cell background at the same concentration), but only pairs of drugs tested more than once were considered in previous analysis [65]. We sought a universal filter to apply to all the data to filter these “transcriptionally silent” drugs. In order to filter these “transcriptionally silent” treatment profiles, drug treatment instances without at least 10 probesets exceeding 2 log2 fold units expression change (increase or decrease) compared to batch mean (representing vehicle treated control) and at least 1 probeset exceeding 3 log2 fold units expression change compared to batch mean (increase) were removed, a heuristic criterion that generated the bimodal distribution displayed in Figure S1. This yielded a final set of 1,419 drug treatment arrays representing 673 unique compounds tested on three cancer cells lines. Probeset values for replicate measurements of a given drug-cell line pair were averaged, yielding 1,033 unique combinations of drug and cell background.

Identifying Differentially Expressed Genes

Differentially expressed probesets were determined by calculating the 2.5 and 97.5 percentiles of log2 expression ranges for a given probeset among the DMSO-treated vehicle controls for a given cell line, and setting to 0 all probeset values in drug treatment instances that were not outside this range.

Clustering and enrichment analysis

Statistical analysis of the resulting drug-induced gene expression profiles was performed in MATLAB R2012a (The Mathworks). The similarity between drug treatment profiles in each cell line were quantified using Pearson correlation, and clustered using affinity propagation, a message-passing algorithm that automatically splits datasets into clusters defined by a set of distinct exemplars (centroids) [28], with the probability of each drug becoming a cluster exemplar being set to the median of pairwise Pearson correlations. To aggregate the data further, the exemplars identified in the initial clustering were also grouped using affinity propagation (with the initial probability of being a cluster exemplar held to the original group median used in the first round of clustering), and this process was repeated hierarchically until the resulting clusters could not be further merged, in a manner similar to previous analysis of the CMap [6]. Chemical similarity between compounds was computed using the Tanimoto coefficient (Jaccard coefficient) using FCFP_6 circular fingerprints calculated in Pipeline Pilot (Scitegic).

To calculate cluster enrichment for hERG inhibitors (>50% reduction of activity at 10 µM or LQT side effect) through permutation testing, labels among compounds experimentally or clinically annotated in our databases were randomized 1000 times and the number of times the resulting cluster enrichments (fraction of hERG inhibitors or LQT drugs among all annotated compounds in a cluster) exceeded the observed number of hERG inhibitors and LQT drugs in the clusters in Figure 2B was computed. This count yielded an empirical p-value (number of times out of 1000 permutations that the enrichment of a randomized cluster exceeded the observed enrichment in the clusters of Figure 2B) which was adjusted using the Benjamini–Hochberg procedure [60] to control for multiple hypothesis testing employing a false discovery rate of 0.2. Clusters were visualized using Cytoscape 2.8.2. Gene Ontology (GO) analysis was performed using the topGO package in bioconductor [66], using the Fisher’s exact test and the elim method for the Biological Process ontology. For each cluster, the set of genes that were up or down-regulated in at least half of the cluster members (median greater or less than 0) were tested for GO term enrichment.

Experimental validation of hERG inhibitors

Ethopropazine, Cloperastine, Fendiline, and Sulconazole (Sigma Aldrich) were prepared at 30 µM stock concentration and serially diluted 3-fold for eight-point dose response measurements. Inhibition of hERG current was experimentally assessed using a previously described protocol [19]. Briefly, Chinese hamster ovary (CHO) cells stably expressing the hERG channel were dislodged from tissue culture flasks and dispensed into PPC plates. Background leak currents were estimated by initiating a 100 ms step to -80 mV from an initial holding potential of -70 mV and subtracted from the subsequent current measurement. For each drug concentration, sequential voltage pulses were applied, each using a 100 ms step to -30 mV from a holding potential of -70 mV, a 2 s conditioning step to +45 mV, and a 2 s test step to -30 mV. Small molecule effects on hERG current density were quantified by measuring the peak tail current prior to compound application and dividing by the amplitude following application of each test concentration. Recordings with peak tail current amplitude pre-compound > 0.2 nA, seal resistance > 30 MOhms, and seal resistance drop rate < 25% were retained for subsequent analysis. Data were fit to a sigmoidal dose response curve using Origin 6.0 (Microcal). To assess the statistical enrichment of this result, we simulated 1000 random sets of four compounds drawn from our experimental data for 10 µM hERG inhibition for the MCF7-tested drugs (red or white nodes in Figure 2A) and compared the mean inhibition of these random sets to that observed for the tested compounds.

Supporting Information

Figure S1

Distribution of pairwise correlations among transcriptionally active and silent drugs.

Replicate drug treatments (duplicate microarray instances for the same concentration, cell line, and drug) may be divided into populations of transcriptionally active (red) and silent (black) drugs based on filters for the magnitude of log2 fold expression change compared to batch average (representing vehicle treated control) among probesets of a given drug.


Figure S2

LQT-annotated drugs have statistically enriched hERG inhibition.

Boxplot of distributional difference in experimentally recorded hERG inhibition values for LQT annotated and unannotated drugs.


Figure S3

Cell line dependent and independent expression changes in hERG-inhibitor enriched clusters.

(A) X, Y, and Z (color gradient of scatterplot points) axes denote average log2 fold gene expression changes versus DMSO treated controls for drugs in hERG-inhibitor enriched clusters highlighted in Figure 2 for drugs profiled in PC3 (x axis) HL60 (y axis) and MCF7 cells (z axis, color gradient of scatterplot points). Purple, red, and blue dashed boxes denote regions of cell-line specific transcriptional modulation. Red or blue shaded points within orange dashed box denote genes with cell line-independent transcriptional changes upon drug treatment. (B) Overlap of differentially expressed (DE) genes with average change in expression greater than 0 versus DMSO treated controls for drugs in the hERG-inhibitor-enriched clusters highlighted in Figure 2B profiled in the three cancer cell lines utilized in the CMap. (C) As in (B), for genes with average fold change less than 0 in the highlighted clusters.


Figure S4

Correlation between hERG activity and transcriptional similarity.

(A) Measured hERG activity (%) is plotted for all assayed drugs versus the maximum correlation (Pearson coefficient) judged by gene expression microarray to one of the five drug expression profiles at the centers of the five enriched cluster exemplars in Figure 2B. (B) Comparison of the correlation in (A) to that in 1,000 sets in which drug activities have been randomly permuted. (C) Fraction of drugs in enriched clusters (red) in Figure 2B for all drugs within a given range (window of activity values with width 10) of measured hERG activity (%).


Figure S5

Statistical evaluation of hERG inhibition among validated compounds.

(A) Distribution of experimentally measured hERG inhibition for all compounds tested on the MCF7 cell line in Figure 2A. (B) Mean hERG activity of random sets of 4 drugs selected from the distribution of (A), compared to the set of validated inhibitors in Figure 4.


Table S1

Drug annotations, SMILES strings.


Table S2

Drug network for MCF7, all cells.


Table S3

GO annotations for drug network clusters.



We thank members of the Li laboratory and Joel Bader for valuable discussions and comments on the manuscript.

Funding Statement

This work is supported by grants to M.L. from the National Institutes of Health (GM078579, MH084691). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Keiser MJ, Setola V, Irwin JJ, Laggner C, Abbas AI et al. (2009) Predicting new molecular targets for known drugs. Nature 462: 175-181.10.1038/nature08506 PubMed: 19881490 [PMC free article] [PubMed]
2. Campillos M, Kuhn M, Gavin AC, Jensen LJ, Bork P (2008) Drug target identification using side-effect similarity. Science 321: 263-266.10.1126/science.1158140 PubMed: 18621671 [PubMed]
3. Fliri AF, Loging WT, Thadeio PF, Volkmann RA (2005) Analysis of drug-induced effect patterns to link structure and side effects of medicines. Nat Chem Biol 1: 389-397.10.1038/nchembio747 PubMed: 16370374 [PubMed]
4. Xie L, Kinnings SL, Bourne PE (2012) Novel computational approaches to polypharmacology as a means to define responses to individual drugs. Annu Rev Pharmacol Toxicol 52: 361-379.10.1146/annurev-pharmtox-010611-134630 PubMed: 22017683 [PubMed]
5. Dudley JT, Sirota M, Shenoy M, Pai RK, Roedder S et al. (2011) Computational repositioning of the anticonvulsant topiramate for inflammatory bowel disease. Sci Transl Med 3: 96ra76 [PMC free article] [PubMed]
6. Iorio F, Bosotti R, Scacheri E, Belcastro V, Mithbaokar P et al. (2010) Discovery of drug mode of action and drug repositioning from transcriptional responses. Proc Natl Acad Sci U S A 107: 14621-14626.10.1073/pnas.1000138107 PubMed: 20679242 [PubMed]
7. Tropsha A, Golbraikh A (2007) Predictive QSAR modeling workflow, model applicability domains, and virtual screening. Curr Pharm Des 13: 3494-3504.10.2174/138161207782794257 PubMed: 18220786 [PubMed]
8. Nigsch F, Lounkine E, McCarren P, Cornett B, Glick M et al. (2011) Computational methods for early predictive safety assessment from biological and chemical data. Expert Opin Drug Metab Toxicol 7: 1497-1511.10.1517/17425255.2011.632632 PubMed: 22050465 [PubMed]
9. Lamb J, Crawford ED, Peck D, Modell JW, Blat IC et al. (2006) The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science 313: 1929-1935.10.1126/science.1132939 PubMed: 17008526 [PubMed]
10. Sirota M, Dudley JT, Kim J, Chiang AP, Morgan AA et al. (2011) Discovery and preclinical validation of drug indications using compendia of public gene expression data. Sci Transl Med 3: 96ra77 [PMC free article] [PubMed]
11. Toyoshiba H, Sawada H, Naeshiro I, Horinouchi A (2009) Similar compounds searching system by using the gene expression microarray database. Toxicol Lett 186: 52-57.10.1016/j.toxlet.2008.08.009 PubMed: 18801419 [PubMed]
12. Lee S, Lee K, Song M, Lee D (2011) Building the process-drug-side effect network to discover the relationship between biological Processes and side effects. BMC Bioinformatics 12: S2.10.1186/1471-2105-12-S1-S2 [PMC free article] [PubMed]
13. Sanguinetti MC, Tristani-Firouzi M (2006) hERG potassium channels and cardiac arrhythmia. Nature 440: 463-469.10.1038/nature04710 PubMed: 16554806 [PubMed]
14. Suessbrich H, Waldegger S, Lang F, Busch AE (1996) Blockade of HERG channels expressed in Xenopus oocytes by the histamine receptor antagonists terfenadine and astemizole. FEBS Lett 385: 77-80.10.1016/0014-5793(96)00355-9 PubMed: 8641472 [PubMed]
15. Mohammad S, Zhou Z, Gong Q, January CT (1997) Blockage of the HERG human cardiac K+ channel by the gastrointestinal prokinetic agent cisapride. Am J Physiol 273: H2534-H2538 PubMed: 9374794 [PubMed]
16. Roy M, Dumaine R, Brown AM (1996) HERG, a primary human ventricular target of the nonsedating antihistamine terfenadine. Circulation 94: 817-823.10.1161/01.CIR.94.4.817 PubMed: 8772706 [PubMed]
17. Itzhaki I, Maizels L, Huber I, Zwi-Dantsis L, Caspi O et al. (2011) Modelling the long QT syndrome with induced pluripotent stem cells. Nature 471: 225-229.10.1038/nature09747 PubMed: 21240260 [PubMed]
18. Moretti A, Bellin M, Welling A, Jung CB, Lam JT et al. (2010) Patient-specific induced pluripotent stem-cell models for long-QT syndrome. N Engl J Med 363: 1397-1409.10.1056/NEJMoa0908679 PubMed: 20660394 [PubMed]
19. Zou B, Yu H, Babcock JJ, Chanda P, Bader JS et al. (2010) Profiling diverse compounds by flux- and electrophysiology-based primary screens for inhibition of human Ether-a-go-go related gene potassium channels. Assay Drug Dev Technol 8: 743-754.10.1089/adt.2010.0339 PubMed: 21158688 [PMC free article] [PubMed]
20. Potet F, Lorinc AN, Chaigne S, Hopkins CR, Venkataraman R et al. (2012) Identification and characterization of a compound that protects cardiac tissue from human Ether-a-go-go-related gene (hERG)-related drug-induced arrhythmias. J Biol Chem 287: 39613-39625.10.1074/jbc.M112.380162 PubMed: 23033485 [PMC free article] [PubMed]
21. Petersen CI, McFarland TR, Stepanovic SZ, Yang P, Reiner DJ et al. (2004) In vivo identification of genes that modify ether-a-go-go-related gene activity in Caenorhabditis elegans may also affect human cardiac arrhythmia. Proc Natl Acad Sci U S A 101: 11773-11778.10.1073/pnas.0306005101 PubMed: 15280551 [PubMed]
22. Rubinstein AL (2006) Zebrafish assays for drug toxicity screening. Expert Opin Drug Metab Toxicol 2: 231-240.10.1517/17425255.2.2.231 PubMed: 16866609 [PubMed]
23. Smith GA, Tsui HW, Newell EW, Jiang X, Zhu XP et al. (2002) Functional up-regulation of HERG K+ channels in neoplastic hematopoietic cells. J Biol Chem 277: 18528-18534.10.1074/jbc.M200592200 PubMed: 11893742 [PubMed]
24. Roy J, Vantol B, Cowley EA, Blay J, Linsdell P (2008) Pharmacological separation of hEAG and hERG K+ channel function in the human mammary carcinoma cell line MCF-7. Oncol Rep 19: 1511-1516 PubMed: 18497958 [PubMed]
25. Berger SI, Ma. ayan A, Iyengar R (2010) Systems pharmacology of arrhythmias. Sci Signal 3: 30 [PMC free article] [PubMed]
26. Du F, Yu H, Zou B, Babcock J, Long S et al. (2011) hERGCentral: a large database to store, retrieve, and analyze compound-human Ether-a-go-go related gene channel interactions to facilitate cardiotoxicity assessment in drug development. Assay Drug Dev Technol 9: 580-588.10.1089/adt.2011.0425 PubMed: 22149888 [PMC free article] [PubMed]
27. Iskar M, Campillos M, Kuhn M, Jensen LJ, van Noort V et al. (2010) Drug-Induced Regulation of Target Expression. PLOS Comput Biol 6: e1000925. [PMC free article] [PubMed]
28. Frey BJ, Dueck D (2007) Clustering by passing messages between data points. Science 315: 972-976.10.1126/science.1136800 PubMed: 17218491 [PubMed]
29. Redfern WS, Carlsson L, Davis AS, Lynch WG, MacKenzie I et al. (2003) Relationships between preclinical cardiac electrophysiology, clinical QT interval prolongation and torsade de pointes for a broad range of drugs: evidence for a provisional safety margin in drug development. Cardiovasc Res 58: 32-45.10.1016/S0008-6363(02)00846-5 PubMed: 12667944 [PubMed]
30. Valentin JP, Bialecki R, Ewart L, Hammond T, Leishmann D et al. (2009) A framework to assess the translation of safety pharmacology data to humans. J Pharmacol Toxicol Methods 60: 152-158.10.1016/j.vascn.2009.05.011 PubMed: 19616110 [PubMed]
31. Kang J, Chen XL, Wang L, Rampe D (2001) Interactions of the antimalarial drug mefloquine with the human cardiac potassium channels KvLQT1/minK and HERG. J Pharmacol Exp Ther 299: 290-296 PubMed: 11561091 [PubMed]
32. Kikuchi K, Nagatomo T, Abe H, Kawakami K, Duff HJ et al. (2005) Blockade of HERG cardiac K+ current by antifungal drug miconazole. Br J Pharmacol 144: 840-848.10.1038/sj.bjp.0706095 PubMed: 15778703 [PMC free article] [PubMed]
33. Mirams GR, Cui Y, Sher A, Fink M, Cooper J et al. (2011) Simulation of multiple ion channel block provides improved early prediction of compounds’ clinical torsadogenic risk. Cardiovasc Res 91: 53-61.10.1093/cvr/cvr044 PubMed: 21300721 [PMC free article] [PubMed]
34. Davies MR, Mistry HB, Hussein L, Pollard CE, Valentin JP et al. (2012) An in silico canine cardiac midmyocardial action potential duration model as a tool for early drug safety assessment. Am J Physiol Heart Circ Physiol 302: H1466-H1480.10.1152/ajpheart.00808.2011 PubMed: 22198175 [PubMed]
35. Huffaker SJ, Chen J, Nicodemus KK, Sambataro F, Yang F et al. (2009) A primate-specific, brain isoform of KCNH2 affects cortical physiology, cognition, neuronal repolarization and risk of schizophrenia. Nat Med 15: 509-518.10.1038/nm.1962 PubMed: 19412172 [PMC free article] [PubMed]
36. Hofmann G, Bernabei PA, Crociani O, Cherubini A, Guasti L et al. (2001) HERG K+ channels activation during beta(1) integrin-mediated adhesion to fibronectin induces an up-regulation of alpha(v)beta(3) integrin in the preosteoclastic leukemia cell line FLG. J Biol Chem 29 1: 276: 4923-4931 [PubMed]
37. Shah RR (2002) The significance of QT interval in drug development. Br J Clin Pharmacol 54: 188-202.10.1046/j.1365-2125.2002.01627.x PubMed: 12207642 [PMC free article] [PubMed]
38. Gaulton A, Bellis LJ, Bento AP, Chambers J, Davies M et al. (2012) ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Res 40: D1100-D1107.10.1093/nar/gkr777 PubMed: 21948594 [PMC free article] [PubMed]
39. Yuill KH, Borg JJ, Ridley JM, Milnes JT, Witchel HJ et al. (2004) Potent inhibition of human cardiac potassium (HERG) channels by the anti-estrogen agent clomiphene-without QT interval prolongation. Biochem Biophys Res Commun 318: 556-561.10.1016/j.bbrc.2004.04.063 PubMed: 15120636 [PubMed]
40. Zhou ZF, Gong QM, Ye B, Fan Z, Makielski JC et al. (1998) Properties of HERG channels stably expressed in HEK 293 cells studied at physiological temperature. Biophys J 74: 230-241.10.1016/S0006-3495(98)77782-3 PubMed: 9449325 [PubMed]
41. Yu SP, Kerchner GA (1998) Endogenous voltage-gated potassium channels in human embryonic kidney (HEK293) cells. J Neurosci Res 52: 612-617.10.1002/(SICI)1097-4547(19980601)52:5 PubMed: 9632317 [PubMed]
42. Sorota S, Zhang XS, Margulis M, Tucker K, Priestley T (2005) Characterization of a hERG screen using the IonWorks HT: Comparison to a hERG rubidium efflux screen. Assay Drug Dev Technol 3: 47-57.10.1089/adt.2005.3.47 PubMed: 15798395 [PubMed]
43. Bridal TR, Margulis M, Wang X, Donio M, Sorota S (2010) Comparison of human Ether-a-go-go related gene screening assays based on IonWorks Quattro and thallium flux. Assay Drug Dev Technol 8: 755-765.10.1089/adt.2010.0267 PubMed: 20658944 [PubMed]
44. Kaczmarek LK (2006) Non-conducting functions of voltage-gated ion channels. Nat Rev Neurosci 7: 761-771.10.1038/nrn1988 PubMed: 16988652 [PubMed]
45. Gomez-Ospina N, Tsuruta F, Barreto-Chang O, Hu L, Dolmetsch R (2006) The C terminus of the L-type voltage-gated calcium channel Ca(V)1.2 encodes a transcription factor. Cell 127: 591-606.10.1016/j.cell.2006.10.017 PubMed: 17081980 [PMC free article] [PubMed]
46. Hegle AP, Marble DD, Wilson GF (2006) A voltage-driven switch for ion-independent signaling by ether-a-go-go K+ channels. Proc Natl Acad Sci U S A 103: 2886-2891.10.1073/pnas.0505909103 PubMed: 16477030 [PubMed]
47. Broccatelli F, Mannhold R, Moriconi A, Giuli S, Carosati E (2012) QSAR Modeling and Data Mining Link Torsades de Pointes Risk to the Interplay of Extent of Metabolism, Active Transport, and hERG Liability. Mol Pharm 9: 2290-2301 [PubMed]
48. Huang X, Dubuc AM, Hashizume R, Berg J, He Y et al. (2012) Voltage-gated potassium channel EAG2 controls mitotic entry and tumor growth in medulloblastoma via regulating cell volume dynamics. Genes Dev 26: 1780-1796.10.1101/gad.193789.112 PubMed: 22855790 [PubMed]
49. Nueda MJ, Sebastian P, Tarazona S, Garcia-Garcia F, Dopazo J et al. (2009) Functional assessment of time course microarray data. BMC Bioinformatics 10Suppl 6: S9.10.1186/1471-2105-10-S1-S9 [PMC free article] [PubMed]
50. Tan Y, Shi L, Hussain SM, Xu J, Tong W et al. (2006) Integrating time-course microarray gene expression profiles with cytotoxicity for identification of biomarkers in primary rat hepatocytes exposed to cadmium. Bioinformatics 22: 77-87.10.1093/bioinformatics/bti737 PubMed: 16249259 [PubMed]
51. Takashima K, Mizukawa Y, Morishita K, Okuyama M, Kasahara T et al. (2006) Effect of the difference in vehicles on gene expression in the rat liver--analysis of the control data in the Toxicogenomics Project Database. Life Sci 78: 2787-2796.10.1016/j.lfs.2005.11.010 PubMed: 16360708 [PubMed]
52. Doddareddy MR, Klaasse EC, Shagufta, Ijzerman AP, Bender A (2010) Prospective validation of a comprehensive in silico hERG model and its applications to commercial compound and drug databases. Chemmedchem 5: 716-729.10.1002/cmdc.201000024 PubMed: 20349498 [PubMed]
53. Li Q, Jorgensen FS, Oprea T, Brunak S, Taboureau O (2008) hERG classification model based on a combination of support vector machine method and GRIND descriptors. Mol Pharm 5: 117-127.10.1021/mp700124e PubMed: 18197627 [PubMed]
54. Obers S, Staudacher I, Ficker E, Dennis A, Koschny R et al. (2010) Multiple mechanisms of hERG liability: K+ current inhibition, disruption of protein trafficking, and apoptosis induced by amoxapine. Naunyn Schmiedebergs Arch Pharmacol 381: 385-400.10.1007/s00210-010-0496-7 PubMed: 20229012 [PubMed]
55. Pillozzi S, Brizzi MF, Bernabei PA, Bartolozzi B, Caporale R et al. (2007) VEGFR-1 (FLT-1), beta1 integrin, and hERG K+ channel for a macromolecular signaling complex in acute myeloid leukemia: role in cell migration and clinical outcome. Blood 110: 1238-1250.10.1182/blood-2006-02-003772 PubMed: 17420287 [PubMed]
56. Zhang H, Zou B, Yu H, Moretti A, Wang X et al. (2012) Modulation of hERG potassium channel gating normalizes action potential duration prolonged by dysfunctional KCNQ1 potassium channel. Proc Natl Acad Sci U S A 109: 11866-11871.10.1073/pnas.1205266109 PubMed: 22745159 [PubMed]
57. Yang KC, Ku YC, Lovett M, Nerbonne JM (2012) Combined deep microRNA and mRNA sequencing identifies protective transcriptomal signature of enhanced PI3Kalpha signaling in cardiac hypertrophy. J Mol Cell Cardiol 53: 101-112.10.1016/j.yjmcc.2012.04.012 PubMed: 22580345 [PMC free article] [PubMed]
58. Song HK, Hong SE, Kim T, Kim do H (2012) Deep RNA sequencing reveals novel cardiac transcriptomic signatures for physiological and pathological hypertrophy. PLOS ONE 7: e35552.10.1371/journal.pone.0035552 PubMed: 22523601 [PMC free article] [PubMed]
59. Bagnall RD, Tsoutsman T, Shephard RE, Ritchie W, Semsarian C (2012) Global microRNA profiling of the mouse ventricles during development of severe hypertrophic cardiomyopathy and heart failure. PLOS ONE 7: e44744.10.1371/journal.pone.0044744 PubMed: 23024758 [PMC free article] [PubMed]
60. Hochberg Y, Benjamini Y (1990) More powerful procedures for multiple significance testing. Stat Med 9: 811-818.10.1002/sim.4780090710 PubMed: 2218183 [PubMed]
61. Li H, Zhou H, Wang D, Qiu J, Zhou Y et al. (2012) Versatile pathway-centric approach based on high-throughput sequencing to anticancer drug discovery. Proc Natl Acad Sci U S A 109: 4609-4614.10.1073/pnas.1200305109 PubMed: 22396588 [PubMed]
62. Bengtsson H, Simpson K, Bullard J, Hansen K (2008) aroma.affymetrix: A generic framework in R for analyzing small to very large Affymetrix data sets in bounded memory.
63. Wu ZJ, Irizarry RA, Gentleman R, Martinez-Murillo F, Spencer F (2004) A model-based background adjustment for oligonucleotide expression arrays. J Am Stat Assoc 99: 909-917.10.1198/016214504000000683
64. R Development Core Team (2012) R: A Language and Environment for Statistical Computing.
65. Ma C, Chen H-IH, Huang Y, Chen Y (2010) Constructing a compound mode-of-action network for personalized drug effectiveness prediction. Proceedings of the First ACM International Conference on Bioinformatics and Computational Biology; Niagara Falls, New York: ACM. 520-528
66. Alexa A, Rahnenfuhrer J (2010) topGO: topGO: Enrichment analysis for Gene Ontology

Articles from PLoS ONE are provided here courtesy of Public Library of Science