{'fraction_non_alphanumeric': 0.038259459078246766, 'fraction_numerical': 0.030431256391269087, 'mean_word_length': 5.158523344191097, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 6, 'lorem ipsum': 0, 'www.': 5, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Systematic polypharmacology and drug repurposing via an integrated L1000-based Connectivity Map database mining

Drug repurposing aims to find novel indications of clinically used or experimental drugs. Because drug data already exist, drug repurposing may save time and cost, and bypass safety concerns. Polypharmacology, one drug with multiple targets, serves as a basis for drug repurposing. Large-scale databases have been accumulated in recent years, and utilization and integration of these databases would be highly helpful for polypharmacology and drug repurposing. The Connectivity Map (CMap) is a database collecting gene-expression profiles of drug-treated human cancer cells, which has been widely used for investigation of polypharmacology and drug repurposing. In this study, we integrated the next-generation L1000-based CMap and an analytic Web tool, the L1000FWD, for systematic analyses of polypharmacology and drug repurposing. Two different types of anti-cancer drugs were used as proof-of-concept examples, including histone deacetylase (HDAC) inhibitors and topoisomerase inhibitors. We identified KM-00927 and BRD-K75081836 as novel HDAC inhibitors and mitomycin C as a topoisomerase IIB inhibitor. Our study provides a prime example of utilization and integration of the freely available public resources for systematic polypharmacology analysis and drug repurposing.


Introduction
Drug repurposing, also known as drug repositioning, aims to identify novel indication(s) of an approved or experimental drug [1]. Polypharmacology describes the ability of a drug to affect more than one molecular target, which has been viewed as a basic property of many therapeutic small molecules and serves as a principle for drug repurposing [2]. Because the safety and efficacy of an approved drug has already been optimized for its original indication, the time for a new indication approval may be significantly shorter and the cost is likely to be far less than those for a new drug. The Connectivity Map (CMap) is a database that collects microarray-based gene-expression profiles from cultured human cancer cell lines treated with various experimentally and clinically used small molecules, and provides a pattern-matching Web-based software to mine these data [3,4]. The current version (build 02; https://portals.broadinstitute.org/cmap/) of CMap collects more than 7000 gene-expression profiles representing 1309 compounds. Because most CMap compounds are the United States Food and Drug Administration (FDA)-approved drugs, this database has become a powerful tool for drug repurposing. By inputting a gene-expression profile of interest and querying it against the CMap data, a list of ranked CMap drugs is obtained. Drugs with positive scores indicate that these CMap drugs and the queried drug have similar gene-expressing signatures, and may share similar mechanism-ofaction (MOA). Thus, the queried drug may be able to treat the same diseases as the original indications of these CMap drugs. Moreover, CMap drugs can be repurposed for novel indications if the gene-expression signatures of established clinical drugs are queried.
Conventional microarray-based gene-expression profiling is expensive and not suitable for high-throughput small-molecule screening. The L1000 assay platform is a new low-cost and highthroughput gene-expression profiling method, which contains 1000 landmark transcripts used to estimate the genome-wide gene-expression signature generated by microarrays [5,6]. L1000 assay combines ligation-mediated amplification, optically addressed and barcoded microspheres (beads), and flow cytometric detection system for gene-expression signature analysis [6]. Based on L1000 assay platform, the Library of Integrated Cellular Signatures (LINCS) is then developed as the nextgeneration version of CMap [5,7]. Compared to CMap, the LINCS collects the gene-expression profiles of not only small-molecule treatments but also genetic manipulation for knocking-down genes by shRNA or over-expressing genes by cDNA. In addition, the LINCS has higher numbers of geneexpression signatures (8870 perturbagens) and cell lines (nine cancer cell lines). Previously, LINCS was accessed via a Web-based interface, the Lincscloud (http://www.lincscloud.org/). Recently, the Lincscloud has been deprecated and replaced by the CLUE (https://clue.io/).
There are several useful resources and databases that are derived from L1000-based LINCS data. For example, the L1000 Characteristic Direction Signature Search Engine (L1000CDS2; http://amp.pharm. mssm.edu/L1000CDS2/) is designed to query gene-expression signatures against the LINCS data to discover and prioritize small molecules that reverse or mimic the inputted gene-expression profile [8]. The CRowd Extracted Expression of Differential Signatures (CREEDS; http://amp.pharm.mssm.edu/ CREEDS/) combines gene-expression profiles from the gene-expression omnibus (GEO), a repository for omics data [9], with LINCS data to identify potential novel drug mimickers [10]. The L1000 Fireworks Display (L1000FWD; http://amp.pharm.mssm.edu/L1000FWD/) provides visualization of drug-induced gene-expression signatures and their similarity in MOA [11]. The integrative LINCS (iLINCS; http://www.ilincs.org/ilincs/) is an integrative Web platform for analysing transcriptomics and proteomics LINCS data. Given the fact that increasingly large-scale databases and analytic tools are developed, the strategies of utilization and integration of these resources have become more and more important.
In this study, we attempt to establish a systematic analysis strategy for polypharmacology and drug repurposing by the integration of the next-generation CMap and L1000FWD. Two types of anti-cancer drugs were used as proof-of-concept examples, including histone deacetylase (HDAC) inhibitors and topoisomerase inhibitors. Our results indicated that the L1000-based next-generation CMap had a convenient and reliable prediction for drug repurposing and polypharmacology. Accompanied analysis of candidate drugs by the L1000FWD further excluded the drugs with dissimilar rsos.royalsocietypublishing.org R. Soc. open sci. 5: 181321 gene signatures. Two novel HDAC inhibitors (KM-00927 and BRD-K75081836) and a novel topoisomerase IIB inhibitor (mitomycin C) were identified. Taken together, this study gives a prime example of systematic polypharmacology analysis and drug repurposing by integrating the freely available resources.

Prediction of novel histone deacetylase inhibitors by the Connectivity Map
Histone acetylation by the histone acetyltransferases (HATs) and deacetylation by the HDACs is essential for regulating gene transcription in eukaryotic cells [12]. The addition of acetyl groups to lysine residues of histones promotes a more-relaxed chromatin structure, allowing transcriptional activation. By contrast, the removal of acetyl group makes the chromatin structure more compact, resulting in gene silencing [13,14]. Because over-expression of HDACs was frequently found in tumours, which prevents tumour suppressor gene expressions, inhibition of HDACs is considered a potential anti-cancer strategy. Currently, several HDAC inhibitors have been approved by the FDA for clinical uses. Vorinostat (suberoylanilide hydroxamic acid; SAHA) and romidepsin (depsipeptide or FK228), have been approved for treating relapsed cutaneous T-cell lymphoma in 2006 and 2009, respectively [15,16]. Panobinostat has been approved for treating multiple myeloma in 2015 [17]. Belinostat has been approved for treating peripheral T-cell lymphoma in 2014 [18].
Transcriptional reprogramming is believed to contribute largely to the therapeutic benefits of HDAC inhibitors [19]. Because CMap is established based on drug transcriptome profiles, we proposed that CMap is a suitable tool for the discovery of novel HDAC inhibitors by comparing the similarity of gene-expression profiling (gene signatures) between the established HDAC inhibitors and CMap drugs. Differentially expressed genes (DEGs, listed in electronic supplementary material, table S1) in HDAC inhibitor-treated cancer cells, including SAHA-treated HCT116 (GSE22061 [20]), SAHA-treated MDA-MB-231 (GSE60125) and FK228-treated HCT116 cells (GSE22061 [20]), were prepared using the GEO online tool GEO2R [9]. These DEGs were queried using CMap (build 02; https://portals. broadinstitute.org/cmap/), and the results were shown in table 1 ranked according to the average scores. CMap drugs with positive mean scores share similar gene-expression patterns with the queried HDAC inhibitors. We found that several HDAC inhibitors (SAHA/vorinostat, trichostatin A, scriptaid, HC toxin and MS-275) were similar to rifabutin (table 1; its chemical structure as shown in figure 1a). Rifabutin is an anti-mycobacterial agent primarily used in the treatment of tuberculosis. Because rifabutin is a semisynthetic ansamycin that has been found to be a heat shock protein 90 (HSP90) inhibitor [21], the mean scores of common HSP90 inhibitors (geldanamycin, tanespimycin, alvespimycin) were also shown as a comparison (table 1). However, these HSP90 inhibitors were not similar to HDAC inhibitors.
To determine the anti-cancer activity of rifabutin against HCT116 cells, MTT assay was performed. The doses of rifabutin less than 10 mM slightly increased the cell viability whereas 50 mM rifabutin reduced the cell viability to 76% (figure 1b). Inhibition of HDAC activity can induce the acetylation of histones. To examine the effect of rifabutin on HDAC activity, whole cell lysates from HCT116 cells treated with rifabutin and SAHA (a positive control) were analysed by Western blotting. As shown in figure 1c, the expression of acetylated histones H3 and H4 increased in response to SAHA, but not rifabutin treatment. Therefore, rifabutin is not an HDAC inhibitor.

Prediction of novel histone deacetylase inhibitors by the Library of Integrated Cellular Signatures
Because the next-generation CMap, the LINCS, has higher numbers of drug-gene signatures than CMap, we also used the LINCS to predict potential HDAC inhibitors. By querying LINCS data with the same DEGs as described above (electronic supplementary material, table S1), we found that two drugs (KM-00927 and BRD-K75081836; their chemical structures are shown in figure 1d) were similar to HDAC inhibitors (table 2). Because only KM-00927 is commercially available, its effects on cell viability and histone H3/H4 acetylation were examined. As shown in figure 1e, KM-00927 dose-dependently killed HCT116 cells. In addition, the acetylation of histones H3/H4 was induced by KM-00927, as SAHA did (figure 1f ), suggesting that KM-00927 is a potential HDAC inhibitor. To investigate whether KM-00927 directly inhibited HDAC, an in vitro HDAC activity assay was performed. We found that rsos.royalsocietypublishing.org R. Soc. open sci. 5: 181321 Table 1. Prediction of potential HDAC inhibitors by the CMap. Differentially expressed genes (DEGs) induced by HDAC inhibitors were obtained from NCBI GEO database, including GSE22061 (SAHA-and FK228-treated HCT116 cells) and GSE60125 (SAHA-treated MDA-MB-231 cells). The full list of DEGs is shown in electronic supplementary material, table S1. These DEGs were inputted to the CMap build 02 and similar drugs were obtained. The mean scores of common HDAC and HSP90 inhibitors are listed in this table, which is ranked by average score that is calculated by adding up the mean scores and dividing the total by the numbers of scores.
CMap drug  Recently, an online Web tool, L1000 Fireworks Display (L1000FWD; http://amp.pharm.mssm.edu/ L1000FWD), has been developed, which can provide visualization of drug-induced gene-expression signatures and their similarity in MOA [11]. As shown in figure 2a (the upper part), the advantage of L1000FWD is that gene signatures from different drug-treat cell types while sharing common MOA can be clustered together. Because CMap and LINCS analyses only gave predictions ranked according to the gene signature similarities, we thought that L1000FWD analysis may provide an additional improvement for CMap/LINCS predictions based on drugs' similarities in MOA. By querying L1000FWD for the predicted CMap and LINCS drugs (figure 1a,d), we found that gene signatures from cells treated with KM-00927 and BRD-K75081836 were clustered in the region of HDAC inhibitors, which was similar to the patterns of vorinostat/SAHA and panobinostat (highlighted in red boxes; figure 2a, the lower part) By contrast, gene signatures from cells treated with rifabutin or geldanamycin were not (figure 2a, the lower part). These results further support our results showing that KM-00927 is an HDAC inhibitor. The above analyses also suggested that BRD-K75081836 may be a potent HDAC inhibitor. Actually, we found that the structure of BRD-K75081836 was very similar to hydroxamic acid-based HDAC inhibitors such as varinostat/SAHA, trichostatin A and scriptaid (figure 2b). Three structural moieties are typically found in hydroxamic acid-based HDAC inhibitors, including a surface recognition cap group that sits outside the active site of an HDAC enzyme and can tolerate extraordinary variability, a zinc-binding group (ZBG) that binds to the buried zinc ion within the active site, and between both, a linker fitting the tunnel of the active site [22]. To quantify the structural similarity between BRD-K75081836 and hydroxamic acid-based HDAC inhibitors, hierarchical clustering of drugs was performed using the ChemBioServer (http://chembioserver.vi-seem.eu/) [23]. Indeed, BRD-K75081836, but not other drugs (KM-00927, rifabutin and geldanamycin), was structurally similar to hydroxamic acid-based HDAC inhibitors including SAHA, trichostatin A, scriptaid and panobinostat. Among them, SAHA was the most similar one to BRD-K75081836 (electronic supplementary material, figure S1). Therefore, we hypothesized that BRD-K75081836 is a direct HDAC inhibitor. To support Table 2. Prediction of potential HDAC inhibitors by the LINCS. Differentially expressed genes (DEGs) induced by HDAC inhibitors were obtained from NCBI GEO database, including GSE22061 (SAHA-and FK228-treated HCT116 cells) and GSE60125 (SAHAtreated MDA-MB-231 cells). The full list of DEGs is shown in electronic supplementary material, table S1. These DEGs were inputted to the Lincscloud website and similar drugs were obtained. this notion, molecular docking analysis was performed to analyse the binding of BRD-K75081836 to the catalytic site of HDAC2. As shown in figure 2c, the binding modes of BRD-K75081836 and SAHA were similar, indicated by the overlapping hydrogen-bonding (green colour) and van der Waals (grey colour) interactions. In addition, the total fitness energy of BRD-K75081836 was higher than that of SAHA (table 3), suggesting the high binding affinity of BRD-K75081836 to HDAC2. These results support that BRD-K75081836 is a hydroxamic acid-based HDAC inhibitor, which warrants further investigation. Furthermore, we performed molecular docking of rifabutin and KM-00927 to HDAC2 in parallel. The results found that the total fitness energies of rifabutin and KM-00927 were lower than those of SAHA and BRD-K75081836 (table 3), further supporting the results that rifabutin and KM-00927 did not interact and inhibit HDAC activity in vitro (figure 1g).
There are several types of HDAC inhibitors with efficacies in the nanomolar to millimolar ranges. The most potent types include cyclic peptides (such as depsipeptide/FK228, apicidin and trapoxin), hydroxomic acids (such as trichostatin A and SAHA) and benzamides (such as entinostat/MS-275 and   mocetinostat), which inhibit HDAC activity in the nanomolar and low micromolar ranges and kill cancer cells in the sub-micromolar ranges. The ketones (such as a-ketomides and trifluoromethyl ketone) and short-chain fatty acids (such as valproic acid, phenylbutyrate and sodium butyrate) exhibit HDAC inhibition and cancer cell killing in the micromolar and millimolar ranges, respectively [24 -26].
Despite the least potency of sort-chain fatty acid to inhibition of HDAC, some of them are used in clinic. For example, valproic acid is used to treat epilepsy and bipolar disorder [27,28]. In addition, valproic acid is still identified as a potent anti-cancer agent and encouraging results are found in clinical trials [29]. KM-00927 is a synthetic drug with unknown function. Our results found that KM-00927 was a potent and indirect HDAC inhibitor, which induced histone acetylation and killed cancer cells in micromolar ranges (10 -50 mM). We thought that KM-00927 could be developed as an anticancer drug in the future. However, whether inhibition of HDAC activity contributes to its anti-cancer activity warrants further investigations.

Prediction of novel topoisomerase inhibitors by integrating the next-generation CMap and L1000FWD
Topoisomerases are enzymes that catalysed the winding and unwinding of DNA double helix through breaking the backbone of either one or both the DNA strands. There are two types of DNA topoisomerases, type I (subtypes: IA, IB and IC) and type II (subtypes: IIA and IIB), which break one or two strands of a DNA helix, respectively [30]. Topoisomerases are proved therapeutic targets of anti-cancer and antibacterial drugs. Clinically successful topoisomerase-targeting anti-cancer drugs act through topoisomerase poisoning, which leads to replication fork arrest and double-strand breaks [30].

Common topoisomerase I inhibitors include camptothecin and its derivatives such as topotecan, irinotecan and SN-38 (active metabolite of irinotecan). Examples of topoisomerase II inhibitors include etoposide, teniposide, doxorubicin and amsacrine.
To further support our drug repurposing strategy, analyses using the next-generation CMap and L1000FWD were performed to predict novel topoisomerase inhibitors. As shown in figure 3, 32 candidate drugs (highlighted in red) were obtained from top 50 perturbations that showed connections with topoisomerase I/II inhibitors. After L1000FWD analysis, we found that gene signatures of topoisomerase I/II inhibitors displayed three major clusters (figure 4a; electronic supplementary material, figure S2; highlighted in green, red and blue boxes). We chose mycophenolate mofetil, mycophenolic acid, mitomycin C and mepacrine for further experiments because of their higher ranking in the CLUE prediction (figure 3). Other candidate drugs with lower ranking were excluded because most of their gene signatures were overlapped with those of CDK inhibitors (electronic supplementary material, figure S3). Some drugs, such as SIB-1893, pitavastatin, gemcitabine and cladribine, have similar patterns with mycophenolate mofetil and mycophenolic acid (electronic supplementary material, figure S3; figure 4b). At the current stage, we only use mycophenolate mofetil/mycophenolic acid as a representative for further analysis.
Mycophenolate mofetil is a prodrug of mycophenolic acid that is an inhibitor of inosine monophosphate dehydrogenase (IMPDH). These two drugs are immunosuppressant drugs used to prevent rejection in organ transplantation [31]. Because only some gene signatures of mycophenolic acid and no gene signatures of mycophenolate mofetil were enriched in the clusters of topoisomerase inhibitors (figure 4b), we speculated that both mycophenolate mofetil and mycophenolic acid may not be topoisomerase inhibitors. Thus, mycophenolic acid was used as a negative example for validation. Mitomycin C, a DNA alkylating agent, has been widely used as chemotherapy treatment for various cancer types [32]. Mepacrine, also called quinacrine, is an anti-malarial drug. Both mitomycin C and mepacrine had gene signatures enriched in the clusters of topoisomerase inhibitors (figure 4b), indicating that they may be potential topoisomerase inhibitors. Indeed, mepacrine has been shown to display anti-cancer activity through inhibition of topoisomerase activity [33]. Therefore, we further examined and compared the effects of mitomycin C and mycophenolic acid (figure 5a). Topoisomerase activity was measured by a band-depletion assay that detects the presence of catalytic topoisomerase-DNA cleavage complexes trapped by topoisomerase inhibitors. Topoisomerase-DNA complexes have lower mobility during SDS-PAGE compared to free enzymes [34]. As shown in figure 5b, mitomycin C, but not mycophenolic acid, depleted the band of topoisomerase IIB (TOP2B), but not topoisomerase 1 (TOP1) and IIA (TOP2A). As positive controls, SN-38 and etoposide (VP-16) depleted TOP1 and TOP2B, respectively. Therefore, mitomycin C acts like etoposide and is a topoisomerase IIB inhibitor. It has been shown that intrinsic sensitivity of cancer cells to topoisomerase II inhibitors is positively and highly correlated with the topoisomerase II expression [35][36][37]. To investigate whether the cell-killing activity of mitomycin C was associated with topoisomerase II expression, drug sensitivity and gene-expression data were obtained by querying the Cancer Therapeutics Response Portal (CTRP [38][39][40]) using the CellMinerCDB (https://discover.nci.nih.gov/cellminercdb/) [41]. As shown in  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  47  48  49  50   name  function  topotecan  camptothecin   SN-38  irinotecan  amonafide  amsacrine  mitoxantrone  etoposide  doxorubicin  teniposide  ellipticine  pirarubicin  pidorubicine   figure 5c, drug activity of mitomycin C was positively correlated with both TOP2A and TOP2B, but not TOP1 mRNA expression. These results supported that topoisomerase II expression could be used as a therapeutic biomarker for mitomycin C.
To investigate the interaction between mitomycin C and topoisomerase IIB, molecular docking was performed. As shown in electronic supplementary material, table S2, the binding affinity of mitomycin C to topoisomerase IIB was lower than that of etoposide, suggesting that additional  mechanism(s) might be involved in the action of mitomycin C. An enzymatic bioreduction of mitomycin C is required for its activation and cytotoxic effect. Upon reduction, mitomycin C is converted into a highly bis-electrophilic intermediate that alkylates and crosslinks DNA [42]. Mitomycin C alkylations are specific for the guanine nucleosides in the CpG sequences [43]. It has been shown that psorospermin, a plantderived anti-tumour agent, alkylates the guanine at the topoisomerase II cleavage site to trap the topoisomerase II-DNA cleavage complex [44]. In addition, clerocidin, a topoisomerase II poison, alkylates and forms covalent adducts with guanines, then trapping topoisomerase II enzymes [45]. These studies support the action of mitomycin C to promote the formation of topoisomerase II-DNA cleavage complex (figure 5b). Therefore, our results provide an additional mechanistic insight into the mode of action of the anti-cancer agent mitomycin C, which may extend its clinical application for treating cancers or other diseases. However, the exact mechanism warrants further investigation.

Conclusion
More and more large-scale biological databases and user-friendly analytic tools are developed, allowing scientists to obtain comprehensive and systematic insights into the scientific problems, even those without a computational background. This study is a prime example of utilization and integration of these freely available public resources. We found that combined analyses using L1000-based CMap/ LINCS and L1000FWD provide reliable prediction for drug repurposing and investigation of drug MOA. According to our analysis, a workflow is therefore proposed ( topoisomerase IIB inhibitor (mitomycin C). The anti-cancer activity and MOA of these drugs will be investigated in detail in the future.

Connectivity Map analysis
For the CMap build 02 analyses, the DEGs (electronic supplementary material, File S1) from HDAC inhibitor-treated cancer cell lines, including SAHA-treated HCT116 (GSE22061 [20]), SAHA-treated MDA-MB-231 (GSE60125) and FK228-treated HCT116 cells (GSE22061 [20]), were prepared using the GEO online tool GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) [9]. GEO2R performs comparisons on original submitter-supplied processed data tables using the GEOquery and limma R packages from the Bioconductor project [9]. The cut-off criteria for DEGs were set at adjusted p-value ,0.01 and jlogFCj. 2. Gene symbols were converted into probe IDs of the Affymetrix GeneChip Human Genome U133A Array (HG-U133A) by an online tool, the NetAffx (https://www.affymetrix. com/analysis/index.affx; Affymetrix). Then, these probe sets were queried according to the instruction. For the LINCS analyses, theses DEGs (gene symbols) were directly inputted in the Lincscloud website (http://www.lincscloud.org/). It should be noted that the Lincscloud website has been deprecated while this study was in progress. For the next-generation CMap (the CLUE) analyses, the online 'Touchstone' tool was used to directly explore the connections of gene signatures of drugs. For the prediction of novel topoisomerase inhibitors, the 'perturbation classes' and 'perturbation type' were set as 'CMap class' and 'compound', respectively. Drugs with description 'topoisomerase inhibitor' were selected for generation of a heat map ( figure 3). The drug similarity was ranked  according to the CMap connectivity score (ranging from 2100 to 100) that corresponds to the fraction of reference gene sets with a greater similarity to the perturbagen than the current query. For example, a score of 95 indicates that only 5% of reference gene sets showed stronger connectivity than the current query to the perturbagen. In general, connectivity scores higher than 95 and lower than 295 were considered as strong scores [7]. For L1000FWD analysis, drug data can be queried and accessed by inputting drug name on the website (http://amp.pharm.mssm.edu/L1000FWD/). Drugs sharing similar MOAs will be clustered together [11].

Chemical structure similarity analysis by the ChemBioServer
The ChemBioServer is a Web application for filtering, clustering and visualizing of chemical compounds [23]. For the hierarchical clustering of compounds, a SDF file containing the simplified molecular-input line-entry system (SMILES) of drugs was uploaded to the ChemBioServer. The parameters were set as follows: distance ¼ Soergel (Tanimoto coefficient); clustering linkage ¼ Ward; and clustering threshold ¼ 0.4.

Molecular docking
Docking simulations were performed with iGEMDOCK [46]. The structures of HDAC2 (4LXZ) and topoisomerase IIB (3QX3) were obtained from the Protein Data Bank. The 3D drug structures in MOL2 format were prepared using the software, Open Babel v. 2.4.1 [47]. The binding site of the target was prepared, and energy-minimized compounds were imported. The docking protocol is set to 80 generations per ligand and a population size of 300 random individuals. All docking conformations were performed 10 times using a genetic evolutionary algorithm, and the fitness of the docked structures was calculated. The binding site for docking was defined by the bound ligand included in the downloaded HDAC2 protein structure, and the binding site radius was set to 8 Å . The binding affinity of a drug to HDAC2 was estimated by the fitness value calculated via the following formula: Fitness (kcal mol 21 ) ¼ VDW þ Hbond þ Elect, where the VDW term is van der Waal energy, Hbond term is hydrogen-bonding energy and Elect term is electro statistic energy.

Cell culture
HeLa human cervical cancer cells were cultured in RPMI-1640 medium. HCT116 human colorectal cancer cells were cultured in McCoy's 5A medium. The culture medium was supplemented with 10% FBS, 2 mM L-glutamine, 1 mM sodium pyruvate, 1% antibiotic-antimycotic. Cells were cultured at 378C and 5% CO 2 in a humidified incubator and passaged every 2 or 3 days.

Cell viability assay
Cell viability was measured by an MTT assay. Cells (3-5 Â 10 3 /well) were plated in 96-well plates and treated with drugs for 72 h. Then, 0.5 mg ml 21 MTT was added. After 4 h, the blue MTT formazan precipitates were dissolved in 200 ml DMSO. The absorbance at 550 nm was measured on a multi-well plate reader.