|
|
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Division of Molecular Pharmacology, Cancer Chemotherapy Center, Japanese Foundation for Cancer Research, Koto-ku, Tokyo, Japan (N.N., T.N., K.Y., T.Y.); Division of Cellular and Molecular Toxicology, Biological Safety Research Center, National Institute of Health Sciences, Setagaya-ku, Tokyo, Japan (N.N., J.K.); and Second Department of Surgery, Tokai University School of Medicine, Boseidai, Isehara-City, Kanagawa, Japan (T.N., S.S., H.M.)
Received June 6, 2007; accepted August 16, 2007
| Abstract |
|---|
|
|
|---|
60%) satisfied our criteria for the further analysis and were classified by cluster analysis of the fingerprints of these chemicals and several standard anticancer drugs into the following three clusters: 1) anticancer drugs, 2) chemicals that shared similar action mechanisms (for example, ouabain and digoxin), and 3) chemicals whose action mechanisms were unknown. These results suggested that chemicals belonging to a cluster (i.e., a cluster of toxic chemicals, a cluster of anticancer drugs, etc.) shared similar action mechanism. In summary, the JFCR39 system can classify chemicals based on their fingerprints, even when their action mechanisms are unknown, and it is highly probable that the chemicals within a cluster share common action mechanisms.
A number of screening methods are currently available for discovering new anticancer drugs. One very powerful and unique approach using multiple cancer cell lines was developed at NCI (Paull et al., 1989
; Weinstein et al., 1992
, 1997
) and in our laboratory (Yamori et al., 1999
; Dan et al., 2002
, 2003
; Yamori, 2003
; Nakatsu et al., 2005
; Akashi and Yamori, 2007
; Akashi et al., 2007
; Nakamura et al., 2007
). This bioinformatics-based approach enables mechanism-oriented evaluation of anticancer drugs. For example, we can evaluate the cell toxicity in vitro by determining the 50% growth inhibition (GI50), total growth inhibition, and 50% lethal concentration across a panel of 39 human cancer cell lines (JFCR39). We can also predict the molecular targets or evaluate the action mechanisms of the test compounds by comparing the cell growth inhibition profiles (termed "fingerprints") across the panel for these compounds with those of the standard anticancer drugs using the COMPARE algorithm (Yamori et al., 1999
). We have used this system successfully and demonstrated that the molecular targets of the novel chemicals MS-274, FJ5002, and ZSTK474 were topoisomerases I and II (Yamori et al., 1999
), telomerase (Naasani et al., 1999
), and phosphatidylinositol 3-kinase (Yaguchi et al., 2006
), respectively. Several other interesting studies, based on a panel of cancer cells, classified anticancer drugs according to their action mechanisms or molecular targets by cluster analysis of their GI50 values (Weinstein et al., 1992
, 1997
; Dan et al., 2002
). Correlation analysis has also been used to explore the genes associated with the sensitivity of the cells in the panel to anticancer drugs (Scherf et al., 2000
; Okutsu et al., 2002
; Zembutsu et al., 2002
; Nakatsu et al., 2005
).
|
|
| Materials and Methods |
|---|
|
|
|---|
Determination of Cell Growth Inhibition Profiles. Growth inhibition experiments were performed to assess the sensitivity of the cells to various chemicals as described before (Yamori et al., 1999
; Dan et al., 2002
). Growth inhibition was measured by determining the changes in the amounts of total cellular protein after 48 h of chemical treatment using a sulforhodamine B assay. For each chemical, the growth assay was performed using a total of five different concentrations of the chemical (for example, 10–4,10–5, 10–6,10–7, and 10–8 M) and one negative control. All assays were performed in duplicate. This GI50 calculation method is well established and reliable through anticancer drug screen using NCI60 as well as JFCR39 (Paull et al., 1989
; Yamori et al., 1999
; Yamori, 2003
). At each test concentration, the percentage growth was calculated using the following seven absorbance measurements: growth at time 0 (T0), growth of the control cells (C), and test growth in the presence of five different concentrations (T) of a drug. The percentage growth inhibition was calculated as: % growth = 100 x [(T – T0)/(C – T0)] when T
T0, and % growth = 100 x [(T – T0)/T] when T < T0. The GI50 values, which represent 50% growth inhibition concentration, were calculated as 100 x [(T – T0)/(C – T0)] = 50. When the GI50 of a chemical could not be calculated, the highest used concentration was assigned as its GI50 value. Absolute values of GI50 were then log transformed for further analysis. We certified the accuracy of measured GI50 data by using reference control chemicals, such as mitomycin-C, paclitaxel, and SN-38, in every experiment and by checking the dose response curves.
Chemicals. Spironolactone, para-aminoazobenzene, para-cresidine, neostigmine bromide, para-dichlorobenzene, phenytoin, ortho-toluidine, imipramine, cobalt chloride, atrazine, propylthiouracil, (D,L)-thalidomide, carbon tetrachloride, hydroquinone, monocrotaline, vinyl chloride, tributyl-tin chloride, valproic acid, benzene, acrylamide, pentachlorophenol, aniline, 1,3-diphenylguanidine, polypropylene glycol, 10,10'-oxy-bis(phenoxyarsine), testosterone propionate, carbaryl, acephate, bisphenol A, 17-
-estradiol, diethylstilbestrol, and
-bungarotoxin were purchased from Wako (Tokyo, Japan). Snake venoms from Agkistrodon halys blomhoffii, Trimeresurus flavoviridis, Crotalus atrox, Naja nigricollis, and Naja naja kaouthia were purchased from Latoxan (Valence, France). 2-Aminomethylpyridine, 1H-1,2,4-triazole, 1H-1,2,3-triazole, 3,4,4'-trichlorocarbanilide, edifenphos, dichlorvos, O-ethyl O-4-nitrophenyl phenylphosphonothioate, 2,4-dinitrophenol, N-methylaniline, 1,2-dichloro-3-nitrobenzene, 4-ethylnitrobenzene, 2-vinylpyridine, 3-amino-1H-1,2,4-triazole, N-ethyl-N-nitrosourea, 5-aza-2'-deoxycytidine, ethynyl estradiol, 3-methylcholanthrene, phenobarbital, acetaminophen, isoniazid, capsaicin, N-deacetyl-N-methylcolchicine (Colcemid), 2.4-dinitrochlorobenzene, and dexamethasone were from Sigma Chemicals (St. Louis, MO). Methoprene acid, methoprene, all-trans retinoic acid, and 9-cis retinoic acid were from BIOMOL International L.P. (Plymouth Meeting, PA). Levothyroxine was from MP Biomedicals (Irvine, CA). 3-Iodo-2-propynyl butylcarbamate was from Olin Japan Inc. (Tokyo, Japan), p-chlorophenyl-3-iodopropargylformal was from Nagase ChemteX (Osaka, Japan), and 2,3,3,3–2',3',3',3'-octachlorodipropylether was from Sankyo Chemical Industries, Ltd. (Tokyo, Japan). 1,2-Benzisothiazolin-3-one was from Riverson (Osaka, Japan), zinc butylxanthate was from Ouchishinko Chemical Industrial Co., Ltd. (Tokyo, Japan), and 4-amino-2,6-dichlorophenol was from Tokyo Kasei Kogyo Co. Ltd. (Tokyo, Japan).
Hierarchical Clustering. Hierarchical clustering analysis was carried out using the average linkage method and the "GeneSpring" software (Silicon Genetics, Inc., Redwood, CA). Pearson correlation coefficients were used to determine the degree of similarity.
| Results |
|---|
|
|
|---|
|
|
Classification of the Chemicals by Hierarchical Clustering. Sixty-nine chemicals were selected for further analysis based on the following criteria: 1) GI50 values for the test chemical can be determined for at least 10 cell lines in the JFCR39 panel, and 2) the range of log GI50 for the test chemical is more than 0.6, suggesting differential growth inhibition. We analyzed the GI50 values of these 69 chemicals and 20 anticancer drugs by hierarchical clustering analysis (Fig. 3). We found approximately 12 clusters (threshold: r = 0, Fig. 3, clusters A–L), which were further divided into 49 subclusters (threshold: r = 0.408, Fig. 3, clusters A1–L6).
|
| Discussion |
|---|
|
|
|---|
We determined the growth inhibition of cells in the JFCR39 panel by 130 chemicals and calculated their GI50 values. Some of the chemicals were assessed twice or more to confirm the reproducibility of the assay. We had to exclude 61 chemicals from further analysis because they did not inhibit the cells in the JFCR39 panel significantly. The rest of the chemicals (69 of 130,
60%) met our selection criteria and were evaluated by cluster analysis.
First, we found that the chemicals tested in duplicate formed tight clusters, showing high reproducibility. Next, we investigated the difference between these 69 test chemicals and the anticancer drugs. Sixty-nine chemicals, which are not anticancer drugs, formed several clusters, which were different from the anticancer drug clusters. These results suggest that the action mechanisms of these chemicals are different from the action mechanisms of the anticancer drugs. However, we found that cisplatin did not belong to the cluster A, which consisted of DNA-targeting anticancer drugs. We do not understand the reason at present. However, there is a possibility that cisplatin has other action mechanisms, which may have made the fingerprint of cisplatin different from those of other DNA-targeting drugs. Indeed, cisplatin is known to form DNA-protein cross-links (Zwelling et al., 1979
; Chválová et al., 2007
).
Our analysis also identified several interesting clusters. For example, the cluster F3 consisted of cardiac glycosides digoxin and ouabain, both of which inhibit Na-K ATPase (Reuter et al., 2002
). The cluster D1 consisted of 9-cis retinoic acid, 13-cis retinoic acid, and 4-[E-2-(5,6,7,8-tetrahydro-5,5,8,8-tetra-methyl-2-naphthalenyl)-1-propenyl]benzoic acid, which are RAR agonists. These results suggest that chemicals other than the anticancer drugs also form clusters when they share the same action mechanisms. It is noteworthy that SV-NNK and SV-NN, from snakes that belonged to the elapidae family, formed one cluster (cluster D2). In contrast, the snake venoms from the C. atrox and T. flavoviridis, which belonged to the viperidae family, formed a cluster (cluster D3) different from the elapidae cluster. These results are reasonable because snake venoms from different snake families are known to differ not only in composition but also in levels of toxicity and mechanisms of action.
The agricultural chemicals paraquat, ziram, and thiram were also classified into a single cluster (cluster F1). Among these agricultural chemicals, the action mechanism of ziram is not known. However, both paraquat and thiram are known to induce oxidative stress (Cereser et al., 2001
; Suntres, 2002
). Therefore, based on our observations, we could suggest that ziram also acted by inducing oxidative stress. The agricultural chemicals methoprene (insect growth regulator) and carbaryl (chorine esterase inhibitor) formed cluster L3, although their common mechanism is unknown. Cluster D4 and D5 consist of the antibacterial agents or fungicides. 3-Iodo-2-propynyl-butylcarbamate and p-chlorophenyl-3'-iodopropargylformal, belonging to cluster D4, are the iodotype antibacterial agents.
Thus, cluster analysis of GI50 values of various chemicals, determined using the JFCR39 cell panel, suggests that the JFCR39 system could, at least in part, allow classification of chemical compounds on the basis of their action mechanisms. Our analysis also suggests that the chemicals belonging the same cluster share a common action mechanism. We are going to develop a larger library of reference chemicals with known action mechanisms (i.e., various inhibitors of biological pathways), and expand our database by integrating their GI50 measurements, which will make the cluster analysis as well as the COMPARE analysis more informative for predicting the mechanism of test chemicals.
In conclusion, to evaluate the potential of the JFCR39 system in predicting the action mechanisms of toxic chemicals, we investigated the fingerprints of 130 different types of chemical compounds including toxic chemicals, pesticides, drugs, and synthetic intermediates. Using the hierarchical clustering analysis, we classified 69 chemicals, at least in part, based on their action mechanisms. Thus, this approach using the JFCR39 cell panel is useful not only in predicting the action mechanisms of toxic chemicals but also in evaluating their toxicity.
| Acknowledgements |
|---|
| Footnotes |
|---|
N. N. and T. N. equally contributed to this study.
Article, publication date, and citation information can be found at http://molpharm.aspetjournals.org.
ABBREVIATIONS: GI50, 50% growth inhibition concentration; GI50, 50% growth inhibition; SN-38, 7-ethyl-10-hydroxycamptothecin; SV-NN, snake venom from N. nigricollis; SV-NNK; snake venom from N. naja kaouthia.
The online version of this article (available at http://molpharm.aspetjournals.org) contains supplemental material. ![]()
Address correspondence to: Takao Yamori, Division of Molecular Pharmacology, Cancer Chemotherapy Center, Japanese Foundation for Cancer Research, 3-10-6, Ariake, Koto-ku, Tokyo 135-8550, Japan. E-mail: yamori{at}jfcr.or.jp, 07a$sl
| References |
|---|
|
|
|---|
Akashi T and Yamori T (2007) A novel method for analyzing phosphoproteins using SELDI-TOF MS in combination with a series of recombinant proteins. Proteomics 7: 2350–2354[CrossRef][Medline]
Aström A, Pettersson U, Krust A, Chambon P, and Voorhees JJ (1990) Retinoic acid and synthetic analogs differentially activate retinoic acid receptor dependent transcription. Biochem Biophys Res Commun 173: 339–345.[CrossRef][Medline]
Cereser C, Boget S, Parvaz P, and Revol A (2001) An evaluation of thiram toxicity on cultured human skin fibroblasts. Toxicology 162: 89–101.[CrossRef][Medline]
Chválová K, Brabec V, and Kasparkova J (2007) Mechanism of the formation of DNA-protein cross-links by antitumor cisplatin. Nucleic Acids Res 35: 1812–1821.
Dan S, Shirakawa M, Mukai Y, Yoshida Y, Yamazaki K, Kawaguchi T, Matsuura M, Nakamura Y, and Yamori T (2003) Identification of candidate predictive markers of anticancer drug sensitivity using a panel of human cancer cell lines. Cancer Sci 94: 1074–1082.[CrossRef][Medline]
Dan S, Tsunoda T, Kitahara O, Yanagawa R, Zembutsu H, Katagiri T, Yamazaki K, Nakamura Y, and Yamori T (2002) An integrated database of chemosensitivity to 55 anticancer drugs and gene expression profiles of 39 human cancer cell lines. Cancer Res 62: 1139–1147.
Hershberger LG, Shipley EG, and Meyer RK (1953) Myotrophic activity of 19-nortestosterone and other steroids determined by modified levator ani muscle method. Proc Soc Exp Biol Med 83: 175–180.[CrossRef][Medline]
Jamieson ER and Lippard SJ (1999) Structure, Recognition, and Processing of Cisplatin-DNA Adducts. Chem Rev 99: 2467–2498.[CrossRef][Medline]
Kanno J, Kato H, Iwata T, and Inoue T (2002) Phytoestrogen-low diet for endocrine disruptor studies. J Agric Food Chem 50: 3883–3885.[CrossRef][Medline]
Naasani I, Seimiya H, Yamori T, and Tsuruo T (1999) FJ5002: a potent telomerase inhibitor identified by exploiting the disease-oriented screening program with COMPARE analysis. Cancer Res 59: 4004–4011.
Nakamura H, Dan S, Akashi T, Unno M, and Yamori T (2007) Absolute quantification of four isoforms of the class I phosphoinositide-3-kinase catalytic subunit by real-time RT-PCR. Biol Pharm Bull 30: 1181–1184.[CrossRef][Medline]
Nakatsu N, Yoshida Y, Yamazaki K, Nakamura T, Dan S, Fukui Y, and Yamori T (2005) Chemosensitivity profile of cancer cell lines and identification of genes determining chemosensitivity by an integrated bioinformatical approach using cDNA arrays. Mol Cancer Ther 4: 399–412.
Okutsu J, Tsunoda T, Kaneta Y, Katagiri T, Kitahara O, Zembutsu H, Yanagawa R, Miyawaki S, Kuriyama K, Kubota N, et al. (2002) Prediction of chemosensitivity for patients with acute myeloid leukemia, according to expression levels of 28 genes selected by genome-wide complementary DNA microarray analysis. Mol Cancer Ther 1: 1035–1042.
Paull KD, Shoemaker RH, Hodes L, Monks A, Scudiero DA, Rubinstein L, Plowman J, and Boyd MR (1989) Display and analysis of patterns of differential activity of drugs against human tumor cell lines: development of mean graph and COMPARE algorithm. J Natl Cancer Inst 81: 1088–1092.
Reuter H, Henderson SA, Han T, Ross RS, Goldhaber JI, and Philipson KD (2002) The Na+-Ca2+ exchanger is essential for the action of cardiac glycosides. Circ Res 90: 305–308.
Scherf U, Ross DT, Waltham M, Smith LH, Lee JK, Tanabe L, Kohn KW, Reinhold WC, Myers TG, Andrews DT, et al. (2000) A gene expression database for the molecular pharmacology of cancer. Nat Genet 24: 236–244.[CrossRef][Medline]
Suntres ZE (2002) Role of antioxidants in paraquat toxicity. Toxicology 180: 65–77.[CrossRef][Medline]
Weinstein JN, Kohn KW, Grever MR, Viswanadhan VN, Rubinstein LV, Monks AP, Scudiero DA, Welch L, Koutsoukos AD, and Chiausa AJ (1992) Neural computing in cancer drug development: predicting mechanism of action. Science 258: 447–451.
Weinstein JN, Myers TG, O'Connor PM, Friend SH, Fornace AJ Jr, Kohn KW, Fojo T, Bates SE, Rubinstein LV, Anderson NL, et al. (1997) An information-intensive approach to the molecular pharmacology of cancer. Science 275: 343–349.
Wong E and Giandomenico CM (1999) Current status of platinum-based antitumor drugs. Chem Rev 99: 2451–2466.[CrossRef][Medline]
Yaguchi S, Fukui Y, Koshimizu I, Matsuno T, Gouda H, Hirono S, Yamazaki K, and Yamori T (2006) Antitumor activity of ZSTK474, a new phosphatidyl inositol 3-kinase inhibitor. J Natl Cancer Inst 98: 545–556.
Yamori T (2003) Panel of human cancer cell lines provides valuable database for drug discovery and bioinformatics. Cancer Chemother Pharmacol 52 (Suppl 1): S74–S79[CrossRef][Medline]
Yamori T, Matsunaga A, Sato S, Yamazaki K, Komi A, Ishizu K, Mita I, Edatsugi H, Matsuba Y, Takezawa K, et al. (1999) Potent antitumor activity of MS-247, a novel DNA minor groove binder, evaluated by an in vitro and in vivo human cancer cell line panel. Cancer Res 59: 4042–4049,.
Zembutsu H, Ohnishi Y, Tsunoda T, Furukawa Y, Katagiri T, Ueyama Y, Tamaoki N, Nomura T, Kitahara O, Yanagawa R, et al. (2002) Genome-wide cDNA microarray screening to correlate gene expression profiles with sensitivity of 85 human cancer xenografts to anticancer drugs. Cancer Res 62: 518–527.
Zwelling LA, Anderson T, and Kohn KW (1979) DNA-protein and DNA interstrand cross-linking by cis- and trans-platinum(II) diamminedichloride in L1210 mouse leukemia cells and relation to cytotoxicity. Cancer Res 39: 365–369.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||