Characterization of drug-induced transcriptional modules:
towards drug repositioning and functional understanding

Murat Iskar, Georg Zeller, Peter Blattmann, Monica Campillos, Michael Kuhn, Katarzyna Kaminska,
Heiko Runz, Anne-Claude Gavin, Rainer Pepperkok, Vera van Noort, Peer Bork

Drug‐induced transcriptional modules (biclusters) were identified and annotated in three human cell lines and rat liver. Analysis of drugs clustered into modules can give new insights into their mechanisms of action and provide leads for drug repositioning.

Correspondence to: Peer Bork,

Supplementary materials for the article :
Characterization of drug-induced transcriptional modules:
towards drug repositioning and functional understanding
Link to paper: PDF   Molecular Systems Biology (2013)


In pharmacology, it is crucial to understand the complex biological responses that drugs elicit in the human organism and how well they can be inferred from model organisms. We therefore identified a large set of drug-induced transcriptional modules from genome-wide microarray data of drug-treated human cell lines and rat liver and first characterized their conservation. Over 70% of these modules were common to multiple cell lines and 15% were conserved between the human in vitro and the rat in vivo system. We then illustrate the utility of conserved and cell-type specific drug-induced modules by predicting and experimentally validating i) gene functions, e.g. ten novel regulators of cellular cholesterol homeostasis and ii) new mechanisms of action for existing drugs thereby providing a starting point for drug repositioning, e.g. novel cell-cycle inhibitors and new modulators of α-adrenergic receptor, peroxisome proliferator-activated receptor and estrogen receptor. Taken together, the identified modules reveal the conservation of transcriptional responses towards drugs across cell types and organisms, and improve our understanding of both the molecular basis of drug action and human biology.

•   Biclustering of drug‐induced gene expression profiles resulted in modules of drugs and genes, which were enriched in both drug and gene annotations.
•   Identifying drug‐induced transcriptional modules separately in three human cell lines and rat liver allows assessment of their conservation across model systems. About 70% of modules are conserved across cell lines, a lower bound of 15% was estimated for their conservation across organisms, and between the in vitro and in vivo systems.
•   Drug‐induced transcriptional modules can predict novel gene functions. A conserved module associated with (chole)sterol metabolism revealed novel regulators of cellular cholesterol homeostasis; 10 of them were validated in functional imaging assays.
•   Analysis of drugs clustered into modules can give new insights into their mechanisms of action and provide leads for drug repositioning. We predicted and experimentally validated novel cell cycle inhibitors and modulators of PPARγ, estrogen and adrenergic receptors, with potential for developing new therapies against diabetes and cancer.

Conserved Drug-Induced Transcriptional Modules

Search Gene or Drug name in CODIMs: (examples: LDLR, Methotrexate)

Complete list of conserved drug-induced transcriptional modules (CODIMs):
CODIM Functional enrichment of genes (GO biological process) Most prominent drug mechanism of action
CODIM1 Cell cycle process, M phase Cell cycle blockers
CODIM2 Sterol biosynthetic process, cholesterol metabolic process Psycholeptics
CODIM3 Citrullination, nucleosome assembly, chromatin assembly Protein synthesis inhibitors
CODIM4 Response to wounding, defense/inflammatory response Corticosteroids
CODIM5 Response to unfolded protein and protein stimulus Heat shock protein inhibitors
CODIM6 Enrichment of LIM domain Histone deactylase inhibitors
CODIM7 Pyrimidine metabolism, (m)RNA processing Flavonoids
CODIM8 Inner ear development, sensory organ development, cell adhesion
CODIM9 Protein folding; keywords: transcriptional regulation
CODIM10 Enrichment of WD40 repeat
CODIM11 Intracellular transport, cellular protein localization
CODIM13 Enrichment of metallothionein, chelation Serpine1
CODIM14 Intracellular protein transport, protein localization
CODIM15 L-serine metabolic process, Glycine, serine and threonine metabolism, Aminoacyl-tRNA biosynthesis SLC6A2
CODIM16 ncRNA metabolic process HRH1
CODIM17 Enrichment of WD40 repeat NA+/K+ pump inhibitors
CODIM18 Enrichment for transmembrane region and signal peptide
CODIM20 Vesicle-mediated transport
CODIM21 Defense/inflammatory response, cell adhesion
CODIM22 Enrichment for glycoprotein, signal peptide
CODIM23 Enrichment Krueppel-associated box domain

MCF7 - Drug-induced modules Functional enrichment of genes Most prominent drug MOA
MCF7-1 (CODIM15) L-serine metabolic process G01; Gynecological antiinfectives and antiseptics
MCF7-2 (CODIM2) Sterol biosynthetic process N05; Psycholeptics
MCF7-3 (CODIM13)
MCF7-4 Retinoid metabolic process
MCF7-5 (CODIM8) Chemotaxis
MCF7-6 (CODIM9) Regulation of RNA metabolic process
MCF7-7 (CODIM3) Positive regulation of transcription from RNA polymerase II promoter CYP2D6
MCF7-8 (CODIM5) Response to unfolded protein  
MCF7-9 G-protein coupled receptor protein signaling pathway ESR1;estrogen receptor 1
MCF7-10 (CODIM10)    
MCF7-11 (CODIM6)    
MCF7-12 (CODIM11)    
MCF7-13 Pathogenic Escherichia coli infection  
MCF7-14 (CODIM1) Cell cycle phase  
MCF7-15 (CODIM4) Leukocyte migration  
MCF7-16 (CODIM12) Protein biosynthesis  
MCF7-18 (CODIM16)    
MCF7-19 (CODIM17) Tube morphogenesis C01; Cardiac therapy
MCF7-20 (CODIM18) Coagulation  
MCF7-21 (CODIM19) Ubiquitin-dependent protein catabolic process TTR;transthyretin
MCF7-22 (CODIM7) RNA processing  
MCF7-23 (CODIM20) Regulation of nucleocytoplasmic transport  
MCF7-24 (CODIM21) Response to wounding  
MCF7-25 (CODIM22) Membrane organization  

PC3 - Drug-induced modules

Functional enrichment of genes Most prominent drug MOA
PC3-1 (CODIM3) Apoptosis  
PC3-2 (CODIM15) Aminoacyl-tRNA biosynthesis HRH1;histamine receptor H1
PC3-3 (CODIM14) Response to unfolded protein R06; Antihistamines for systemic use
PC3-4 (CODIM8) Chemotaxis  
PC3-5 Response to hypoxia  
PC3-6 (CODIM16) ncRNA processing SERPINE1;serpin peptidase inhibitor, clade E, member 1
PC3-7 (CODIM4) Response to wounding D07; Corticosteroids, dermatological preparations
PC3-8 (CODIM1) Mitotic cell cycle  
PC3-9   PPAR; Peroxisome proliferator-activated receptor
PC3-10 (CODIM10) Actin cytoskeleton organization  
PC3-11 Complement and coagulation cascades ADRB2; Adrenergic, beta-2-, receptor, surface
PC3-12 (CODIM23) Regulation of RNA metabolic process  
PC3-13 (CODIM7) mRNA processing  
PC3-14 (CODIM2) Steroid metabolic process  
PC3-15 (CODIM5) Response to unfolded protein  
PC3-16 (CODIM6)    
PC3-17 (CODIM9) Nuclear transport  
PC3-18 (CODIM12)    
PC3-19 (CODIM18)    
PC3-20 Protein transport  
PC3-22 (CODIM11) Intracellular protein transport  
PC3-23 (CODIM19) Protein folding  
PC3-24 Polysaccharide biosynthetic process  
PC3-25 (CODIM20)    
PC3-26   DRD4;dopamine receptor D4
PC3-27 (CODIM21) Cell adhesion  
PC3-28 (CODIM22) Integrin-mediated signaling pathway  

HL60 - Drug-induced modules

Functional enrichment of genes Most prominent drug MOA
HL60-1 (CODIM13) Response to virus  
HL60-2 (CODIM5) Response to unfolded protein  
HL60-3 (CODIM1) Cell cycle phase P01; Antiprotozoals
HL60-4 Chemotaxis N05; Psycholeptics
HL60-5 (CODIM11) Translational elongation  
HL60-6 (CODIM10)    
HL60-7 (CODIM12)    
HL60-8 (CODIM9) mRNA metabolic process  
HL60-9 (CODIM14) Establishment of protein localization CYP1B1;cytochrome P450 1B1
HL60-11 Vesicle-mediated transport  
HL60-12 (CODIM6)    
HL60-13 (CODIM7) RNA splicing CYP1B1;cytochrome P450 1B1
HL60-15 (CODIM4) Response to bacterium D07; Corticosteroids, dermatological preparations
HL60-16 (CODIM8) Inner ear development  
HL60-17 Regulation of lymphocyte activation ADRB2;adrenergic, beta-2-, receptor, surface
HL60-18 (CODIM17)   C01; Cardiac therapy
HL60-19 (CODIM12) Nucleic acid transport  
HL60-20 (CODIM3) Nucleosome assembly  
HL60-22 Vesicle-mediated transport  
HL60-23 Regulation of transcription  
HL60-25 (CODIM2) Steroid biosynthetic process  
HL60-28 ER to Golgi vesicle-mediated transport  
HL60-29 Cell adhesion  

Rat Liver - Drug-induced modules

Functional enrichment of genes Most prominent drug MOA
RatLiver-1 Fatty acid metabolic process  
RatLiver-2 Apoptosis  
RatLiver-3 Regulation of protein kinase activity  
RatLiver-4 Response to xenobiotic stimulus  
RatLiver-5 Response to DNA damage stimulus  
RatLiver-8 Modification-dependent protein catabolic process  
RatLiver-9 Translational elongation  
RatLiver-10 Organic acid catabolic process  
RatLiver-12 Sterol biosynthetic process Statins
RatLiver-13 Phosphate metabolic process  
RatLiver-15 Energy derivation by oxidation of organic compounds  
RatLiver-16 Fatty acid metabolic process  
RatLiver-17 Metabolism of xenobiotics by cytochrome P450  
RatLiver-20 Aminoacyl-tRNA biosynthesis  
RatLiver-22 Chemotaxis  
RatLiver-23 Positive regulation of macromolecule biosynthetic process  
RatLiver-24 Translational elongation  
RatLiver-25 Fatty acid metabolic process  
RatLiver-27 RNA splicing, via transesterification reactions  
RatLiver-31 Positive regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolic process  
RatLiver-32 Regulation of protein transport  
RatLiver-33 Regulation of cell migration  
RatLiver-34 Muscle tissue development  
RatLiver-35 Nitrogen compound biosynthetic process  
RatLiver-36 Cell cycle  
RatLiver-38 Positive regulation of developmental process  
RatLiver-39 Protein transport  
RatLiver-40 Protein localization  
RatLiver-41 Protein transport  
RatLiver-42 Protein modification by small protein conjugation or removal  
RatLiver-43 Negative regulation of gene expression  

Supplementary Figure 1. Properties and preprocessing of input datasets. Row labeled ‘Genes’ contains the results of redundancy removal and orthology mapping between probe sets and charts the overlap as Venn diagram in the middle. Row labeled ‘Drugs’ charts the result of restricting the CMap data set to compounds tested in all three cell lines used here and how these overlap with compounds tested in rat liver (Venn diagram in the middle). Pie charts give the composition of the drug space profiled in each resource.

Supplementary Figure 2. Non-redundant drug-induced modules can be detected across a wide range of ISA parameter settings. The non-redundant final set of drug-induced modules from each cell line was compared to the whole set of raw drug-induced modules before redundancy removal (MCF7 n=2255, PC3 n=1661, HL60 n=1783) that were identified under a number of runs of the ISA algorithm that explored an extensive parameter grid (gene s.d. cutoff [from 5 to 2 in steps of 0.2], drug s.d. cutoff [from 4 to 1 in steps of 0.2]). The heatmap indicates the significance of gene overlap (Fisher’s exact test, negative log10 of the p-value, see color key) between any of the final non-redundant modules (columns) and the most similar module from a particular ISA parameter set (rows). In total, more than half of the final modules significantly overlapped (FDR-corrected p-value<0.01) with raw modules from more than 100 distinct parameter settings.

Supplementary Figure 3. Network of conserved drug-induced transcriptional modules (CODIM). Using a reciprocal best hit approach, drug-induced transcriptional modules across three cell lines and rat liver were linked to each other based on the significance of the overlap between their gene and drug members (with FDR corrected p-value < 0.01). Taken together, these modules covered a substantial fraction of the gene and drug space of the underlying expression data: on average 1,857 genes (20% of the input gene set) and 555 drugs (56%) per cell line, and for rat liver 1,587 genes (44%) and 311 drugs (90%) were represented within modules. All drug-induced transcriptional modules were labeled (number) in the order in which they were identified by the ISA procedure.

Supplementary Figure 4. Module conservation across model systems is significantly higher than expected for random independent biclusters. For each pair-wise comparison of 4 datasets (3 cancer cell lines and rat liver), the number of conserved modules was compared against the conservation between biclusters that were randomly generated, independently for each data set. Properties of 10,000 randomly generated biclusters resembled those of real ones (same distributions of gene members). In each pair-wise comparison between data sets, we observed significantly more conserved real drug-induced modules as assessed by their gene overlap (empirical p-value < 1E-4) than between independent randomly generated data sets.

Supplementary Figure 5. Assessing the impact of drug space on module conservation. The total number of treatments (990 drugs) in the CMap dataset were downsampled to smaller drug spaces comprising 10% to 100% (raw data) in increments of 10% (always using the same drug sample for all three cell lines). Five replicates were analyzed for each downsampled drug space. Drug-induced transcriptional modules were identified using exactly the same pipeline and parameters except for the number of seeds, which were reduced to 2,000 to reduce computation time. The total number of identified drug-induce modules and their conservation across cell lines and rat liver are displayed as stacked bar plots with error bars representing the standard error. These results indicate that a drug space of 200-400 drugs (~20-40%, assessed in three cell lines) appears to be a lower limit for robust conservation estimation.

Supplementary Figure 6. Functional coherence of transcriptional modules as benchmarked against STRING. Functionally related protein pairs were extracted from the STRING protein-protein association network for human and rat. Drug-induced transcriptional modules from human cell lines and rat liver were found to have a higher ratio of functionally associated gene pairs in comparison to random modules (1,000 gene sets with the same size and background). 30 out of 82 modules from human cell lines (35%), 12 out of 23 CODIMs (52%) and 13 out of 43 modules from rat liver (30%) clustered together more functionally related genes than expected by chance (at an empirical, permutation-based p-value cutoff < 0.05).

Supplementary Figure 7. CODIM2 is enriched in genes known to function in unfolded protein response (UPR) and cholesterol biosynthesis pathways. (A) The modulatory effect of CODIM2 on three branches of UPR signaling. Shown are genes known to be modulated by one of the key UPR sensors IRE-1, ATF6, PERK and how their expression changes in response to the drugs contained in CODIM2 (enriched for antipsychotics). Overall mean fold changes ±SEM in PC3 module 14, n=28, are in the following indicated in parentheses: for IRE1 targets HERPUD1 (1.75±0.11), EDEM1 (0.32±0.04), DNAJB9 (1.68±0.18), DNAJC3 (0.75±0.11), and HSPA5 (0.89±0.07); for the ATF6 branch ATF6 (0.2±0.03), GADD45 (0.68±0.07), HSP90B1 (0.48±0.06) and XBP1 (0.79±0.06); for PERK targets ATF4 (0.82±0.07), ATF3 (1.66±0.24), DDIT3 (1.92±0.14), DDIT4 (1.92±0.21), and TRIB3 (1.75±0.13) (for genes with multiple probe sets, the one with highest mean FC was selected). In contrary to reported effects of antipsychotics specifically on PERK signaling (Canfran-Duque et al, 2012), CODIM2 drugs were found to stimulate all branches of UPR signaling although to a lower extent in the ATF6 branch. Target genes were extracted from the literature (Canfran-Duque et al, 2012; Lee et al, 2003; Szegezdi et al, 2006). (IRE1, inositol required 1; PERK, PKR-like endoplasmic reticulum localized kinase; ATF6, activating transcription factor 6) (B) The curated pathway of cholesterol biosynthesis as illustrated in WikiPathways (Kelder et al, 2012). Genes in this pathway that were also contained in CODIM2 are labeled in green.

Supplementary Figure 8. Negative controls for assessing anti-proliferative effects of CODIM1 predictions. Three predicted candidate cell-cycle inhibitors were chosen for validation experiments: sulconazole (anti-fungal) and vinburnine (vasodilator) were successfully validated, whereas validation of mephentermine (cardiac stimulant) failed (see Fig. 3 in the main text). Butoconazole and moxisylate (both present in CMap) were chosen as negative controls for validation experiments with results shown here.

Supplementary Figure 9. Cell cycle analysis using propidium iodide staining and FACS. and assessing FL3-A signal to determine cell cycle phases of cells. (A) Treatment of HL-60 cells with reference compound nocodazole (200 ng/mL) resulted in G2/M arrest, most prominent after 24h of treatment. (B) Using “Flowing software”, the FL3-A signal indicating DNA content was used to determine cells in different cell cycle phases, subG1, G1, G2/M and endoR (endoreduplicated) (see Fig. 3 and Materials and Methods for details).

Supplementary Figure 10. Detailed results of experimentally tested drug-target interactions. Ten newly predicted drug-target interactions were experimentally evaluated with in vitro binding or cellular assays as indicated. Two drugs from the PC3-9 module, zaprinast (phosphodiesterase inhibitor) and raubasine (anti-hypertensive) were tested for PPARγ binding. Four drug repositioning candidates from MCF7-9 module were screened for binding affinity to estrogen receptor alpha (nitrendipine, antihypertensive; theobromine, alkaloid of the cacao plant, present at high levels in chocolate; bendroflumethiazide, thiazide diuretic) or an antagonistic effect in a cellular assay (dilazep, vasodilator). Predictions based on the HL60-17 module (hexetidine, anti-bacterial; (+)-chelidonine, alkaloid; vigabatrin, antiepileptic; podophyllotoxin, antiviral) were both tested against ADRA2C and ADRA1B, as it was not possible to a priori determine the exact target associated with the drugs from the HL60-17 module from cell line data. Ki values were obtained for predictions with 40% and more activity at 50 µM concentration (nh: Hill coefficient).

Supplementary Table 1. Robustness of drug-induced modules to redundancy removal parameters. Redundancies among modules were eliminated within each cell line using the recommended standard procedure ISAUnique (R ‘eisa’ package). We used a correlation (cor.) threshold of 0.3. Subsequently, modules were subjected to a second round of redundancy removal based on gene overlap between modules (Fisher test, p-value <10-5, using the same priority as pre-defined by ISA). Additionally, we explored a wide range of parameters here (cor. threshold: from 1 to 0.1 in decrements of 0.1 and overlap p-value cutoffs 10-2, 10-5, and 10-10) and assessed the effect on the resulting non-redundant set of modules (percentages in parentheses indicate the coverage of the non-redundant set of modules described in the main text). In conclusion, the combination of these redundancy removal routines yields similar results across a broad range of parameter settings and is thus quite robust to the actual parameters used (shaded in gray).

Supplementary Table 2. Robustness of constitutive co-expression modules to detection parameter. To obtain a final set of drug-induced modules, we discarded those ‘constitutive’ co-expression modules in which 10% or more of the gene pairs were correlated (Pearson coefficient > 0.6) in untreated control samples. Here we additionally explored a wider parameter grid by varying the percentage of correlated gene pairs between 1 and 50% (see Table). When we assessed how many of the non-redundant modules filtered with these parameters overlapped with the final modules presented in the main text, we found that for any cutoff, the majority of final drug-induced modules (in most cases the vast majority) were detected invariably.

Supplementary Table 3 (xlsx). Comparison of gene and drug members of drug-induced modules across human cell lines and organisms (separate file).

Supplementary Table 4. Characterization of gene and drug members of conserved drug-induced modules (CODIMs), CODIM-associated modules of rat liver and selected cell line-specific modules linked to a certain MOA. Functional enrichment of genes (biological process) and MOA enrichment of drugs were labeled in green if identified in the automated annotation process. For drug annotation, manually curated drug information was additionally gathered for CODIMs (colored in orange). Moreover, we checked the consistency of gene and drug annotations and categorized these associations in terms of their agreement with the current literature: green labels for previously reported associations, yellow for less obvious cases with poorly characterized drug mechanism of action, and grey for novel associations for which no clear understanding could be gathered from the literature (cases without detectable and interpretable enrichment in white).

Supplementary Table 5. Detailed information on the characterization of gene and drug members of unified conserved drug-induced modules (CODIMs), CODIM-associated modules of rat liver and selected cell line-specific modules linked to a certain MOA.

Supplementary Table 6. Experimental results of RNAi-screen for selected genes from CODIM2. We found that siRNA mediated knock-down of 10 out of 24 poorly characterized genes (not previously linked to cellular cholesterol regulation) have a significant effect (labeled in orange) in the filipin and/or the DiI-LDL uptake assays with an inhibitory (yellow) or stimulatory (blue) role. An effect was deemed significant if 2 or more siRNAs had a consistent effect with absolute (average) Z-score >1 and FDR corrected p-value <0.01.

Here is the bundled supplementary materials file with captions and references. (PDF document).

Iskar M, Zeller G, Blattmann P, Campillos M, Kuhn M, Kaminska KH, Runz H, Claude-Gavin A, Pepperkok R, van Noort V, Bork P. (2013) Characterization of drug-induced transcriptional modules: towards drug repositioning and functional understanding. Molecular Systems Biology 9. doi:10.1038/msb.2013.20.

PDF Pubmed ID: 23632384


Drug-induced transcriptional modules in mammalian biology: Implications for drug repositioning and resistance

PDF Slides

Dagstuhl Seminar, May 19 – 24, 2013 Wadern (Germany)
"Computational Methods Aiding Early-Stage Drug Design"


Drug-induced transcriptional modules in mammalian biology: implications for drug repositioning and resistance

The recent explosion of publicly available data on small molecules and their biological effects has had a profound impact on research in pharmacology and toxicology. It enables integrative analysis of various large-scale heterogeneous resources, which holds great promise for attaining a better understanding of drug action at various biological levels [1]. One type of these new resources catalogs genome-wide expression profiles of cells treated with drugs. They summarize the effects of drug perturbations at the transcriptome level and thus provide a rich database for computational methods that aim to bridge between molecular basis of drug action and higher-level phenotypic consequences such as therapeutic indications and side effects. In this study, an integrative framework was developed and applied to analyze large-scale microarray compendia of drug treatments to systematically identify and predict functional features related to both drug and gene action. As a prerequisite, we implemented a new methodology to overcome batch effect in drug-induced gene expression profiles, which substantially improved the integrity and reliability of the signals obtained for cell-type specific analysis [2]. By associating drug-target interactions with drug-induced gene expression profiles, the newly developed method revealed which drug targets are differentially regulated at the expression level upon modulation by drugs. This analysis provided a lower limit (8%) of drug targets that are subjected to feedback regulation and identified novel cases that might be implicated in drug tolerance and drug resistance [2].

Having established these advanced analysis tools, we set out to more comprehensively characterize and dissect the expression responses of mammalian cells to drug perturbations. To this end, we exploited an unsupervised biclustering approach to identify a comprehensive set of drug-induced transcriptional modules for three human cell lines and rat liver [3]. For the first time, we quantified the conservation of these modules between human cell lines (72%) and between human and the rat model organism (15%). In-depth analysis of drug-induced modules allowed us to gain insights into drug action and human molecular biology in various ways. Specifically, we predicted new therapeutic indications for a dozen of marketed drugs (drug repositioning/drug repurposing), and as a proof-of-principle, experimentally validated two of these drug repositioning candidates as novel cell cycle blockers that can potentially be used for chemotherapy. Additionally, drug-induced modules suggested functional roles for hypothetical genes. We predicted several poorly annotated genes to be novel regulators of cellular cholesterol levels and successfully validated a large proportion of these experimentally [3].

In summary, the integrative framework presented here has comprehensively delineated the expression responses of mammalian cells to drug perturbations into a collection of drug-induced transcriptional modules ( [3]. This valuable and versatile resource of drug-induced transcriptional modules not only led to a deeper understanding of drug action with implications for drug repositioning, but also provided new insights into biological networks of rat animal model and the human organism.

1. Iskar M, Zeller G, Zhao X-M, van Noort V, Bork P (2011) Drug discovery in the age of systems biology: the rise of computational approaches for data integration. Current Opinion Biotechnology 23(4):609-16.
Link to paper: PDF Pubmed: 22153034

2. Iskar M, Campillos M, Kuhn M, Jensen LJ, van Noort V, et al. (2010) Drug-induced regulation of target expression. PLoS Computational Biology 9;6(9). pii: e1000925.
Link to paper: PDF Pubmed: 20838579
3. Iskar M, Zeller G, Blattmann P, Campillos M, Kuhn M, et al. (2013) Characterization of drug-induced transcriptional modules: towards drug repositioning and functional understanding. Molecular Systems Biology 30;9:662.
Link to paper: PDF Pubmed: 23632384

6507 visitors since March 15th, 2013.