Search

SUMMARY CellMiner-SCLC (https://discover.nci.nih.gov/SclcCellMinerCDB/) integrates drug sensitivity and genomic data, including high-resolution methylome and transcriptome from 118 patient-derived small cell lung cancer (SCLC) cell lines, providing a resource for research into this “recalcitrant cancer.” We demonstrate the reproducibility and stability of data from multiple sources and validate the SCLC consensus nomenclature on the basis of expression of master transcription factors NEUROD1, ASCL1, POU2F3, and YAP1. Our analyses reveal transcription networks linking SCLC subtypes with MYC and its paralogs and the NOTCH and HIPPO pathways. SCLC subsets express specific surface markers, providing potential opportunities for antibody-based targeted therapies. YAP1-driven SCLCs are notable for differential expression of the NOTCH pathway, epithelial-mesenchymal transition (EMT), and antigen-presenting machinery (APM) genes and sensitivity to mTOR and AKT inhibitors. These analyses provide insights into SCLC biology and a framework for future investigations into subtype-specific SCLC vulnerabilities., In Brief Tlemsani et al. provide a unique resource, SCLC-CellMiner, integrating drug sensitivity and multi-omics data from 118 small cell lung cancer (SCLC) cell lines. They demonstrate that SCLCs have differential transcriptional networks driven by lineage-specific transcription factors (NEUROD1, ASCL1, POU2F3, and YAP1). Furthermore, YAP1-driven SCLCs have distinct drug sensitivity profiles., Graphical Abstract

High-throughput and high-content databases are increasingly important resources in molecular medicine, systems biology, and pharmacology. However, the information usually resides in unwieldy databases, limiting ready data analysis and integration. One resource that offers substantial potential for improvement in this regard is the NCI-60 cell line database compiled by the U.S. National Cancer Institute, which has been extensively characterized across numerous genomic and pharmacologic response platforms. In this report, we introduce a CellMiner (http://discover.nci.nih.gov/cellminer/) web application designed to improve the use of this extensive database. CellMiner tools allowed rapid data retrieval of transcripts for 22,379 genes and 360 microRNAs along with activity reports for 20,503 chemical compounds including 102 drugs approved by the U.S. Food and Drug Administration. Converting these differential levels into quantitative patterns across the NCI-60 clarified data organization and cross-comparisons using a novel pattern match tool. Data queries for potential relationships among parameters can be conducted in an iterative manner specific to user interests and expertise. Examples of the in silico discovery process afforded by CellMiner were provided for multidrug resistance analyses and doxorubicin activity; identification of colon-specific genes, microRNAs, and drugs; microRNAs related to the miR-17-92 cluster; and drug identification patterns matched to erlotinib, gefitinib, afatinib, and lapatinib. CellMiner greatly broadens applications of the extensive NCI-60 database for discovery by creating web-based processes that are rapid, flexible, and readily applied by users without bioinformatics expertise. Cancer Res; 72(14); 3499–511. ©2012 AACR.

The “apoptotic ring” is characterized by the phosphorylation of histone H2AX at serine 139 (γ-H2AX) by DNA-dependent protein kinase (DNA-PK). The γ-H2AX apoptotic ring differs from the nuclear foci patterns observed in response to DNA-damaging agents. It contains phosphorylated DNA damage response proteins including activated Chk2, activated ATM, and activated DNA-PK itself but lacks MDC1 and 53BP1, which are required to initiate DNA repair. Because DNA-PK can phosphorylate heat shock protein 90α (HSP90α) in biochemical assays, we investigated whether HSP90α is involved in the apoptotic ring. Here we show that HSP90α is phosphorylated by DNA-PK on threonines 5 and 7 early during apoptosis and that both phosphorylated HSP90α and DNA-PK colocalize in the apoptotic ring. We also show that DNA-PK is a client of HSP90α and that HSP90α is required for full DNA-PK activation, γ-H2AX formation, DNA fragmentation, and apoptotic body formation. In contrast, HSP90 inhibition by geldanamycin markedly enhances TRAIL-induced DNA-PK and H2AX activation. Together, our results reveal that HSP90α is a substrate and chaperone of DNA-PK in the apoptotic response. The response of phosphorylated HSP90α to TRAIL and its localization to the γ-H2AX ring represent epigenetic features of apoptosis that offer insights for studying and monitoring nuclear apoptosis.

RNA splicing is required to remove introns from pre-mRNA, and alternative splicing generates protein diversity. Topoisomerase I (Top1) has been shown to be coupled with splicing by regulating serine/arginine-rich splicing proteins. Prior studies on isolated genes also showed that Top1 poisoning by camptothecin (CPT), which traps Top1 cleavage complexes (Top1cc), can alter RNA splicing. Here, we tested the effect of Top1 inhibition on splicing at the genome-wide level in human colon carcinoma HCT116 and breast carcinoma MCF7 cells. The RNA of HCT116 cells treated with CPT for various times was analyzed with ExonHit Human Splice Array. Unlike other exon array platforms, the ExonHit arrays include junction probes that allow the detection of splice variants with high sensitivity and specificity. We report that CPT treatment preferentially affects the splicing of splicing-related factors, such as RBM8A, and generates transcripts coding for inactive proteins lacking key functional domains. The splicing alterations induced by CPT are not observed with cisplatin or vinblastine and are not simply due to reduced Top1 activity, as Top1 downregulation by short interfering RNA did not alter splicing like CPT treatment. Inhibition of RNA polymerase II (Pol II) hyperphosphorylation by 5,6-dichloro-1-β-d-ribofuranosylbenzimidazole (DRB) blocked the splicing alteration induced by CPT, which suggests that the rapid Pol II hyperphosphorylation induced by CPT interferes with normal splicing. The preferential effect of CPT on genes encoding splicing factors may explain the abnormal splicing of a large number of genes in response to Top1cc. Cancer Res; 70(20); 8055–65. ©2010 AACR.

Aminoflavone (AF, NSC 686,288) is beginning clinical trials. It induces replication-mediated histone H2AX phosphorylation, DNA-protein crosslinks and activates p53. Here, we studied p21(CIP1/WAF1) and Mdm2 responses to AF. Although p53 stabilization and phosphorylation at serine 15 increased with dose and time of exposure, Mdm2 and p21(CIP1/WAF1) protein levels displayed a biphasic response, as they accumulated at submicromolar doses and then decreased with increasing AF. As both Mdm2 and p21(CIP1/WAF1) mRNA levels increased with AF concentration without reduction at higher concentrations, we measured the half-lives of Mdm2 and p21(CIP1/WAF1) proteins. Mdm2 and p21(CIP1/WAF1) half-lives were shortened with increasing AF concentrations. Proteasomal degradation appears responsible for the decrease of both Mdm2 and p21(CIP1/WAF1), as MG-132 prevented their degradation and revealed AF-induced Mdm2 polyubiquitylation. AF also induced protein kinase B (Akt) activation, which was reduced with increasing AF concentrations. Suppression of Akt by small interfering RNA was associated with downregulation of Mdm2 and p21(CIP1/WAF1) and with enhanced apoptosis. These results suggest that the cellular responses to AF are determined at least in part by Mdm2 and p21(CIP1/WAF1) protein levels, as well as by Akt activity, leading either to cell cycle arrest when Mdm2 and p21(CIP1/WAF1) are elevated, or to apoptosis when Mdm2 and p21(CIP1/WAF1) are degraded by the proteasome and Akt insufficiently activated to protect against apoptosis.

RecQ helicase BLM-deficient cells are characteristically hypersensitive to 4-nitroquinoline-1-oxide (4NQO). We recently reported that isogenic BLM-deficient cells (PNSG13) are more sensitive than BLM-complemented cells (PNSF5) to camptothecin, which specifically traps topoisomerase I cleavage complexes (Top1cc). We now report that PNSG13 are also 3.5-fold more sensitive to 4NQO compared with PNSF5 and that 4NQO induces higher levels of Top1cc and reduced histone γ-H2AX in PSNG13 than in PNSF5. Similarly, 4NQO induces more Top1cc in primary fibroblasts from a patient with Bloom syndrome than in normal human fibroblasts. 4NQO also induces Top1cc in colon cancer HCT116 and HT29 cells in a time- and concentration-dependent fashion. Of note, distinct from camptothecin, the Top1cc produced by 4NQO accumulate progressively after 4NQO addition and persist following 4NQO removal. The Top1cc induced by 4NQO are detectable by alkaline elution. To examine the functional relevance of the Top1cc induced by 4NQO, we used two stable topoisomerase I small interfering RNA (siRNA) cell lines derived from HCT116 and MCF7 cells. Both topoisomerase I siRNA cell lines are resistant to 4NQO, indicating that Top1cc contribute to the cellular activity of 4NQO. Collectively, these data show that 4NQO is an effective inducer of cellular Top1cc. Because 4NQO does not directly trap Top1cc in biochemical assays, we propose that active metabolites of 4NQO trap Top1cc by forming DNA adducts. Induction of Top1cc and histone γ-H2AX by 4NQO may contribute to the cellular effects of 4NQO, including its selective activity toward RecQ helicase BLM-deficient cells.(Cancer Res 2006; 66(13): 6540-5)

Camptothecin, originally discovered in 1957 as an antitumor activity in plant extracts, has recently become one of the most promising leads to new anticancer drugs. After lingering for many years, interest in camptothecin was revitalized in 1985 upon discovery of its specific action on topoisomerase I. Detailed elucidation of action mechanisms at the molecular, cellular, and pharmacologic levels has made camptothecin and its congeners perhaps the best understood among clinical anticancer drugs. Promising chemical variants of camptothecin, and recently other chemical categories of topoisomerase I-targeted drugs, provide unusually rich opportunities for rational drug selection and design. This is made possible by current concepts based, for the most part, on a sound experimental foundation, which points the way towards optimally effective therapy.

This study provides evidence for the importance of p21(CDKN1A) for the repair of replication-mediated DNA double-strand breaks (DSBs) induced by topoisomerase I. We report that defects of p21(CDKN1A) and p53 enhance camptothecin-induced histone H2AX phosphorylation (gammaH2AX), a marker for DNA DSBs. In human colon carcinoma HCT116 cells with wild-type (wt) p53, gammaH2AX reverses after camptothecin removal. By contrast, gammaH2AX increases after camptothecin removal in HCT116 cells deficient for p53 (p53-/-) or p21(CDKN1A) (p21-/-) as the cells reach the late-S and G2 phases. Since p21-/- cells exhibit similar S-phase arrest as wt cells in response to camptothecin and aphidicolin does not abrogate the enhanced gammaH2AX formation in p21-/- cells, we conclude that enhanced gammaH2AX formation in p21-/- cells is not due to re-replication. The cell cycle checkpoint abrogator and Chk1/Chk2 inhibitor 7-hydroxystaurosporine (UCN-01) also increases camptothecin-induced gammaH2AX formation and inhibits camptothecin-induced p21(CDKN1A) upregulation in HCT116 wt cells. TUNEL (terminal deoxynucleotidyl transferase (TdT)-mediated dUTP-biotin nick end labeling) assays demonstrate that gammaH2AX formation in late S and G2 cells following CPT treatment corresponds to DNA breaks. However, these breaks are not related to apoptotic DNA fragmentation. We propose that p21(CDKN1A) prevents the collapse of replication forks damaged by stabilized topoisomerase I cleavage complexes.

A standard for bioregulatory network diagrams is urgently needed in the same way that circuit diagrams are needed in electronics. Several graphical notations have been proposed, but none has become standard. We have prepared many detailed bioregulatory network diagrams using the molecular interaction map (MIM) notation, and we now feel confident that it is suitable as a standard. Here, we describe the MIM notation formally and discuss its merits relative to alternative proposals. We show by simple examples how to denote all of the molecular interactions commonly found in bioregulatory networks. There are two forms of MIM diagrams. “Heuristic” MIMs present the repertoire of interactions possible for molecules that are colocalized in time and place. “Explicit” MIMs define particular models (derived from heuristic MIMs) for computer simulation. We show also how pathways or processes can be highlighted on a canonical heuristic MIM. Drawing a MIM diagram, adhering to the rules of notation, imposes a logical discipline that sharpens one's understanding of the structure and function of a network.

Intrinsic (innate) and acquired (adaptive) resistance to chemotherapy critically limits the outcome of cancer treatments. For many years, it was assumed that the interaction of a drug with its molecular target would yield a lethal lesion, and that determinants of intrinsic drug resistance should therefore be sought either at the target level (quantitative changes or/and mutations) or upstream of this interaction, in drug metabolism or drug transport mechanisms. It is now apparent that independent of the factors above, cellular responses to a molecular lesion can determine the outcome of therapy. This review will focus on programmed cell death (apoptosis) and on survival pathways (Bcl-2, Apaf-1, AKT, NF-kappaB) involved in multidrug resistance. We will present our molecular interaction mapping conventions to summarize the AKT and IkappaB/NF-kappaB networks. They complement the p53, Chk2 and c-Abl maps published recently. We will also introduce the 'permissive apoptosis-resistance' model for the selection of multidrug-resistant cells.

Topoisomerase inhibitors are among the most efficient inducers of apoptosis. The main pathways leading from topoisomerase-mediated DNA damage to cell death involve activation of caspases in the cytoplasm by proapoptotic molecules released from mitochondria. In some cells, apoptotic response also involves the death receptor Fas (APO-1/CD95). The engagement of these apoptotic effector pathways is tightly controlled by upstream regulatory pathways that respond to DNA lesions-induced by topoisomerase inhibitors in cells undergoing apoptosis. These include the proapoptotic Chk2, c-Abl and SAPK/JNK pathways, the survival PI(3)kinase-Akt-dependent pathway and the transcription factors p53 and NF-kappaB. Initiation of cellular responses to DNA lesions-induced by topoisomerase inhibitors is ensured by the protein kinases DNA-PK, ATM and ATR, which bind to DNA breaks. These kinases commonly called "DNA sensors" mediate their effects (DNA repair, cell cycle arrest and/or apoptosis) by phosphorylating a large number of substrates, including several downstream kinases such as c-Abl and the checkpoint protein Chk2. c-Abl induces apoptosis by activating cell death pathways (e.g., SAPK, p53 and p73) and inhibiting cell survival pathways [e.g., PI(3)kinase]. The DNA-damage regulating kinase Chk2, in addition to its role in cell cycle arrest and/or DNA repair, can induce apoptosis by phosphorylation/activation of the promyelocytic leukemia (PML) protein and p53. Finally, we will review the recent observations that support a role for topoisomerases in chromatin fragmentation during the execution phase of apoptosis.

Recent results suggest that potentially lethal DNA lesions may result when replication forks encounter trapped topoisomerase-DNA complexes or some other types of DNA damage. Such events produce what are called replication-encounter lesions. These lesions have the characteristic that they may allow single stranded DNA-associated replication protein A (RPA) to become juxtaposed to dsDNA end-associated DNA-protein kinase. Our results suggest that DNA-protein kinases may then hyperphosphorylate the RPA2 subunit. We discuss a possible pathway by which hyperphosphorylation of RPA2 could lead to the release of active p53. This could constitute a pathway for signaling the presence of replication-encounter lesions to the p53-dependent cell cycle arrest and/or apoptosis initiator systems. *** DIRECT SUPPORT *** A01Q2011 00008

The present study was undertaken to examine structural features of L-chicoric acid (3) which are important for potency against purified HIV-1 integrase and for reported cytoprotective effects in cell-based systems. Through a progressive series of analogues, it was shown that enantiomeric D-chicoric acid (4) retains inhibitory potency against purified integrase equal to its L-counterpart and further that removal of either one or both carboxylic functionalities results in essentially no loss of inhibitory potency. Additionally, while two caffeoyl moieties are required, attachment of caffeoyl groups to the central linking structure can be achieved via amide or mixed amide/ester linkages. More remarkable is the finding that blockage of the catechol functionality through conversion to tetraacetate esters results in almost no loss of potency, contingent on the presence of at least one carboxyl group on the central linker. Taken as a whole, the work has resulted in the identification of new integrase inhibitors which may be regarded as bis-caffeoyl derivatives of glycidic acid and amino acids such as serine and beta-aminoalanine. The present study also examined the reported ability of chicoric acid to exert cytoprotective effects in HIV-infected cells. It was demonstrated in target and cell-based assays that the chicoric acids do not significantly inhibit other targets associated with HIV-1 replication, including reverse transcription, protease function, NCp7 zinc finger function, or replication of virus from latently infected cells. In CEM cells, for both the parent chicoric acid and selected analogues, antiviral activity was observable under specific assay conditions and with high dependence on the multiplicity of viral infection. However, against HIV-1- and HIV-2-infected MT-4 cells, the chicoric acids and their tetraacetylated esters exhibited antiviral activity (50% effective concentration (EC50) ranging from 1.7 to 20 microM and 50% inhibitory concentration (IC50) ranging from 40 to 60 microM).

An alkylating camptothecin (CPT) derivative, 7-chloromethyl-10,11-methylenedioxy-camptothecin (7-CM-MDO-CPT) was recently shown to produce irreversible topoisomerase I (top1) cleavage complexes by binding to the +1 base of the scissile strand of a top1 cleavage site. We demonstrate that 7-CM-EDO-CPT (7-chloromethyl-10,11-ethylenedioxy-camptothecin) also induces irreversible top1-DNA complexes. 7-CM-MDO-CPT, 7-CM-EDO-CPT, and the nonalkylating derivative 7-ethyl-10,11-methylenedioxy-camptothecin (7-E-MDO-CPT) also induced reversible top1 cleavable complexes, which were markedly more stable to salt-induced reversal than those induced by 7-ethyl-10-hyroxy-CPT, the active metabolite of CPT-11. This greater stability of the top1 cleavable complexes was contributed by the 7-alkyl and the 10,11-methylene- (or ethylene-) dioxy substitutions. Studies in SW620 cells showed that 7-E-MDO-CPT, 7-CM-MDO-CPT, and 7-CM-EDO-CPT are more potent inducers of cleavable complexes and more cytotoxic than CPT. The reversal of the cleavable complexes induced by 7-E-MDO-CPT, 7-CM-MDO-CPT, and 7-CM-EDO-CPT was markedly slower after drug removal than that for CPT, which is consistent with the data with purified top1. By contrast to CPT, 7-E-MDO-CPT, 7-CM-MDO-CPT, and 7-CM-EDO-CPT were cytotoxic irrespective of the presence of 10 microM aphidicolin. These results suggest that 7-E-MDO-CPT, 7-CM-MDO-CPT, and 7-CM-EDO-CPT are more potent top1 poisons than CPT and produce long lasting top1 cleavable complexes and greater cytotoxicity than CPT in cells.

Host site selection for full-site integration by human immunodeficiency virus type-1 (HIV-1) intergrase (IN) from nonionic detergent-lysed virions was investigated. Linear retrovirus-like DNA (469 bp) possessing 3' OH recessed long terminal repeat termini was efficiently inserted by a bimolecular donor reaction into a supercoiled DNA target (2867 bp), producing the HIV-1 5-bp host site duplication. Sequence data were analyzed from 193 donor-target recombinants obtained from the linear 3.8-kb DNA product. The selection of host target sites appeared randomly distributed and was independent of lysis and assay conditions. The fidelity of the 5-bp duplications in comparison to other size duplications was highest (94%) with high-salt (300 mM NaCl) lysis of the virions and 60 mM NaCl for strand transfer using Mg2+ as the divalent cation. Base sequence analysis demonstrated some biases in the 5-bp duplications at the sites of strand transfer and at the immediate host sequences surrounding the duplications. In addition to the observed duplications, approximately 30% of the recombinants isolated from the linear 3.8-kb DNA product contained specific and repetitive small-size deletions. No deletions smaller that 17 bp were observed and the distance between the deletion sets had a periodicity of approximately 10 bp. The mechanisms involved in how HIV-1 IN produces the 5-bp duplications and the repetitive host site deletions are discussed.

Topoisomerase I (Top1) relaxes DNA supercoiling by forming transient cleavage complexes (Top1cc) up- and downstream of transcription complexes. Top1cc can be trapped by carcinogenic and endogenous DNA lesions and by camptothecin, resulting in transcription blocks. Here, we undertook genome-wide analysis of camptothecin-treated cells at exon resolution. RNA samples from HCT116 and MCF7 cells were analyzed with the Affy Exon Array platform, allowing high-resolution mapping along 18,537 genes. Long genes that are highly expressed were the most susceptible to downregulation, whereas short genes were preferentially upregulated. Along the body of genes, downregulation was most important toward the 3′-end and increased with the number of exon–intron junctions. Ubiquitin and RNA degradation-related pathway genes were selectively downregulated. Parallel analysis of microRNA with the Agilent miRNA microarray platform revealed that miR-142-3p was highly induced by camptothecin. More than 10% of the downregulated genes were targets of this p53-dependent microRNA. Our study shows the profound impact of Top1cc on transcription elongation, especially at intron–exon junctions and on transcript stability by microRNA miR-142-3p upregulation. Cancer Res; 73(15); 4830–9. ©2013 AACR.

This personal account traces a series of studies that led from DNA physical chemistry to anticancer drug mechanisms. Chemical crosslinking as a basis for anticancer drug actions had been suspected since the time of the first clinical reports of the effectiveness of nitrogen mustard in 1946. After the elucidation of the DNA helix-coil transition, several nearly concurrent findings in the early 1960s established the paradigm of DNA interstrand crosslinking. The DNA filter elution phenomenon was discovered in the early 1970s, and lent itself to the development of practical assays for DNA crosslinks and other DNA lesions in mammalian cells. The assays allowed studies of the effects of DNA damaging agents at pharmacologically or toxicologically relevant doses, and have been widely applied in studies of mutagenic and chemotherapeutic agents. During the period 1979-1986, DNA filter elution studies led to the paradigm of DNA topoisomerases as targets of anticancer drug action, and this has become one of the most active areas of anticancer drug development.

Although there is extensive information on gene expression and molecular interactions in various cell types, integrating those data in a functionally coherent manner remains challenging. This study explores the premise that genes whose expression at the mRNA level is correlated over diverse cell lines are likely to function together in a network of molecular interactions. We previously derived expression-correlated gene clusters from the database of the NCI-60 human tumor cell lines and associated each cluster with function categories of the Gene Ontology (GO) database. From a cluster rich in genes associated with GO categories related to cell migration, we extracted 15 genes that were highly cross-correlated; prominent among them were RRAS, AXL, ADAM9, FN14, and integrin-beta1. We then used those 15 genes as bait to identify other correlated genes in the NCI-60 database. A survey of current literature disclosed, not only that many of the expression-correlated genes engaged in molecular interactions related to migration, invasion, and metastasis, but that highly cross-correlated subsets of those genes engaged in specific cell migration processes. We assembled this information in molecular interaction maps (MIMs) that depict networks governing 3 cell migration processes: degradation of extracellular matrix, production of transient focal complexes at the leading edge of the cell, and retraction of the rear part of the cell. Also depicted are interactions controlling the release and effects of calcium ions, which may regulate migration in a spaciotemporal manner in the cell. The MIMs and associated text comprise a detailed and integrated summary of what is currently known or surmised about the role of the expression cross-correlated genes in molecular networks governing those processes.

The present observations show that staurosporine can rapidly trigger both the morphological changes and intranucleosomal DNA fragmentation typical of apoptosis. This occurred in a number of cell lines from various origins regardless of the state of differentiation and cell cycle phase, suggesting the presence of a common inducible suicide pathway. The broad apoptotic activity of staurosporine appears to be unique among other protein kinase or phosphatase inhibitors we tested. Results obtained in a cell-free assay suggest that cytoplasmic proteins directly modulated by staurosporine may be involved in a ubiquitous signal for the induction of DNA fragmentation and apoptosis.

A remarkable overlap was observed between the gadd genes, a group of often coordinately expressed genes that are induced by genotoxic stress and certain other growth arrest signals, and the MyD genes, a set of myeloid differentiation primary response genes. The MyD116 gene was found to be the murine homolog of the hamster gadd34 gene, whereas MyD118 and gadd45 were found to represent two separate but closely related genes. Furthermore, gadd34/MyD116, gadd45, MyD118, and gadd153 encode acidic proteins with very similar and unusual charge characteristics; both this property and a similar pattern of induction are shared with mdm2, whic, like gadd45, has been shown previously to be regulated by the tumor suppressor p53. Expression analysis revealed that they are distinguished from other growth arrest genes in that they are DNA damage inducible and suggest a role for these genes in growth arrest and apoptosis either coupled with or uncoupled from terminal differentiation. Evidence is also presented for coordinate induction in vivo by stress. The use of a short-term transfection assay, in which expression vectors for one or a combination of these gadd/MyD genes were transfected with a selectable marker into several different human tumor cell lines, provided direct evidence for the growth-inhibitory functions of the products of these genes and their ability to synergistically suppress growth. Taken together, these observations indicate that these genes define a novel class of mammalian genes encoding acidic proteins involved in the control of cellular growth.

The National Cancer Institute currently tests approximately 400 compounds per week against a panel of human tumour cell lines in order to identify potential anti-cancer drugs. We describe several approaches, based on these in vitro data, to the problem of identifying the primary biochemical mechanism of action of a compound. Using linear and non-parametric discriminant procedures and cross-validation, we find that accurate identification of the mechanism of action is achieved for approximately 90 per cent of a diverse collection of 141 known compounds, representing six different mechanistic categories. We demonstrate that two-dimensional graphical displays of the compounds in terms of the initial three principal components (of the original data) result in suggestive visual clustering according to mechanism of action. Finally, we compare the classification accuracy of the statistical discrimination procedures with the accuracy obtained from a neural network approach and, for our example, we find that the results obtained from the various approaches are similar.

The National Cancer Institute's drug discovery program screens more than 20,000 chemical compounds and natural products a year for activity against a panel of 60 tumor cell lines in vitro. The result is an information-rich database of patterns that form the basis for what we term an "information-intensive" approach to the process of drug discovery. The first step was a demonstration, both by statistical methods (including the program COMPARE) and by neural networks, that patterns of activity in the screen can be used to predict a compound's mechanism of action. Given this finding, the overall plan has been to develop three large matrices of information: the first (designated A) gives the pattern of activity for each compound tested against each cell line in the screen; the second (S) encodes any of a number of types of 2-D or 3-D structural motifs for each compound; the third (T) indicates each cell's expression of molecular targets (e.g., from 2-dimensional protein gel electrophoresis). Construction and updating of these matrices is an ongoing process. The matrices can be concatenated in various ways to test a variety of specific hypotheses about compounds screened, as well as to "prioritize" candidate compounds for testing. To aid in these efforts, we have developed the DISCOVERY program package, which integrates the matrix data for visual pattern recognition. The "information-intensive" approach summarized here in some senses serves to bridge the perceived gap between screening and structure-based drug design.

The aim of this work was to identify nitrogen mustards that would react selectively with DNA, particularly in G-rich regions. A series of mono- and difunctional nitrogen mustards was synthesized in which the (2-chloroethyl)amino functions were connected to the N 9 of 9-aminoacridine by way of a spacer chain consisting of two to six methylene units. The length of the spacer chain connecting the alkylating and putative DNA-intercalating groups was found to affect the preference for the alkylation of different guanine-N 7 positions in a DNA sequence. All of the compounds reacted preferentially at G's that are followed by G as do most other types of nitrogen mustards, but the degree of selectivity was greater. The compounds reacted at much lower concentrations than were required for comparable reaction by mechlorethamine (HN2), consistent with initial noncovalent binding to DNA prior to guanine-N 7 alkylation. The degree of DNA-sequence selectivity increased as the spacer-chain length decreased below four methylene units. Most strikingly, long spacer compounds reacted strongly at 5'-GT-3' sequences, whereas this reaction was almost completely suppressed when the spacer length was reduced to two or three methylenes. Mono- and difunctional compounds of a given spacer length showed no consistent difference in DNA-sequence preference

Camptothecin-induced DNA photolesions were examined after UVA irradiation at 365 nm. DNA single-strand breaks were induced both in supercoiled and in relaxed SV40 DNA. In uniquely end-labeled human c-myc DNA, camptothecin-induced cleavage occurred exclusively at guanines and was markedly enhanced by hot piperidine treatment. Runs of polyguanines were the most cleaved, especially in their 5' flank. Primer extension experiments in the absence of piperidine treatment confirmed these results and did not show additional lesions. We found that synthetic single-stranded oligonucleotides were more reactive than duplex oligonucleotides. In addition, an excess of dideoxyguanosine triphosphates competed for camptothecin-induced DNA photolesions. Therefore, camptothecin stacking in DNA grooves is more likely than genuine drug intercalation. Groove shielding with sodium or magnesium reduced camptothecin-induced photodamage while minor groove occupancy with spermine extended damages. Photolesion mechanisms were investigated using scavengers. In aerobic conditions, the most effective scavengers were thiourea, sodium azide, and catalase. Protection by superoxide dismutase was weak, and mannitol was ineffective. In anaerobic conditions, lesions were more extensive. Taken together, these results show that photoactivated camptothecin interacts specifically and intimately with guanines. This finding is consistent with preferential stimulation of topoisomerase I cleavage at sites that bear a guanine at their 5'-DNA terminus [Jaxel, C., et al. (1991) J. Biol. Chem. 266, 1465-1469] and with the camptothecin stacking model at topoisomerase I DNA cleavage sites.

The effects of monocytic/macrophage and granulocytic differentiation induced by phorbol myristate acetate (TPA) and all-trans retinoic acid, respectively, were tested on the induction of apoptosis in human promyelocytic leukemia HL-60 cells treated with topoisomerase I and II inhibitors. Using a filter-binding assay, we observed a strong inhibition of DNA fragmentation induced by 3- and 24-hour continuous exposure to camptothecin, VP-16, VM-26, and m-AMSA in TPA- differentiated cells. The inhibition of the typical internucleosomal DNA fragmentation was confirmed by agarose gel electrophoresis. By contrast, drug-induced DNA fragmentation was not inhibited in retinoic acid-differentiated cells, and apoptosis occurred in these cells after 4 to 5 days in the absence of drug treatment. The TPA inhibitory effect was maximal after 24 hours of treatment and was correlated with differentiation, because phorbol dibutyrate ester was active, whereas 4- alpha-TPA, a nontumor promoter that does not induce differentiation, was not active. Using alkaline elution, we observed that TPA and retinoic acid differentiation were associated with changes in topoisomerase-mediated DNA breaks that were not correlated with their differential effects on drug-induced DNA fragmentation. Moreover, TPA also inhibited DNA fragmentation induced by vinblastine, cycloheximide, calphostin C, and x-rays. Using a cell-free system, we observed that DNA fragmentation was not inhibited in nuclei from TPA-differentiated cells. Rather, inhibition of apoptosis seemed to take place in the cytoplasm. We conclude that phenotypic changes associated with TPA- induced differentiation include inactivation of a cytoplasmic activity that can induce DNA fragmentation associated with apoptosis.

In this study, we further examined the sequence selectivity of camptothecin in mammalian topoisomerase I cDNA from human and Chinese hamster. In the absence of camptothecin, almost all the bases at the 3'-terminus of cleavage sites are T for calf thymus and wheat germ topoisomerase I. In addition, wheat germ topoisomerase I exhibits preference for C (or not T) at -3 and for T at -2 position. As for camptothecin-stimulated cleavage with topoisomerase I, G (or not T) at +1 is an additional strong preference. This sequence selectivity of camptothecin is similar to that previously found in SV40 DNA, suggesting that camptothecin preferentially interacts with topoisomerase I-mediated cleavage sites where G is the base at the 5'-terminus. These results support the stacking model of camptothecin (Jaxel et al. (1991) J. Biol. Chem. 266, 20418-20423). Comparison of calf thymus and wheat germ topoisomerase I-mediated cleavage sites in the presence of camptothecin shows that many major cleavage sites are similar. However, the relative intensities are often different. One of the differences was attributable to a bias at position -3 where calf thymus topoisomerase I prefers G and wheat germ topoisomerase I prefers C. This difference may explain the unique patterns of cleavage sites induced by the two enzymes. Sequencing analysis of camptothecin-stimulated cleavage sites in the surrounding regions of point mutations in topoisomerase I cDNA, which were found in camptothecin-resistant cell lines, reveals no direct relationship between DNA cleavage sites in vitro and mutation sites.

As part of the Spotlight on Molecular Profiling series, we present here new profiling studies of mRNA and microRNA expression for the 60 cell lines of the National Cancer Institute (NCI) Developmental Therapeutics program (DTP) drug screen (NCI-60) using the 41,000-probe Agilent Whole Human Genome Oligo Microarray and the 15,000-feature Agilent Human microRNA Microarray V2. The expression levels of ∼21,000 genes and 723 human microRNAs were measured. These profiling studies include quadruplicate technical replicates for six and eight cell lines for mRNA and microRNA, respectively, and duplicates for the remaining cell lines. The resulting data sets are freely available and searchable online in our CellMiner database. The result indicates high reproducibility for both platforms and an essential biological similarity across the various cell types. The mRNA and microRNA expression levels were integrated with our previously published 1,429-compound database of anticancer activity obtained from the NCI DTP drug screen. Large blocks of both mRNAs and microRNAs were identified with predominately unidirectional correlations to ∼1,300 drugs, including 121 drugs with known mechanisms of action. The data sets presented here will facilitate the identification of groups of mRNAs, microRNAs, and drugs that potentially affect and interact with one another. Mol Cancer Ther; 9(5); 1080–91. ©2010 AACR.

Described here are neural networks capable of predicting a drug's mechanism of action from its pattern of activity against a panel of 60 malignant cell lines in the National Cancer Institute's drug screening program. Given six possible classes of mechanism, the network misses the correct category for only 12 out of 141 agents (8.5 percent), whereas linear discriminant analysis, a standard statistical technique, misses 20 out of 141 (14.2 percent). The success of the neural net indicates several things. (i) The cell line response patterns are rich in information about mechanism. (ii) Appropriately designed neural networks can make effective use of that information. (iii) Trained networks can be used to classify prospectively the more than 10,000 agents per year tested by the screening program. Related networks, in combination with classical statistical tools, will help in a variety of ways to move new anticancer agents through the pipeline from in vitro studies to clinical application.

Ataxia telangiectasia mutated (ATM), the deficiency of which causes a severe neurodegenerative disease, is a crucial mediator for the DNA damage response (DDR). As neurons have high rates of transcription that require topoisomerase I (TOP1), we investigated whether TOP1 cleavage complexes (TOP1cc)-which are potent transcription-blocking lesions-also produce transcription-dependent DNA double-strand breaks (DSBs) with ATM activation. We show the induction of DSBs and DDR activation in post-mitotic primary neurons and lymphocytes treated with camptothecin, with the induction of nuclear DDR foci containing activated ATM, gamma-H2AX (phosphorylated histone H2AX), activated CHK2 (checkpoint kinase 2), MDC1 (mediator of DNA damage checkpoint 1) and 53BP1 (p53 binding protein 1). The DSB-ATM-DDR pathway was suppressed by inhibiting transcription and gamma-H2AX signals were reduced by RNase H1 transfection, which removes transcription-mediated R-loops. Thus, we propose that Top1cc produce transcription arrests with R-loop formation and generate DSBs that activate ATM in post-mitotic cells.

When replication is blocked by a template lesion or polymerase inhibitor while helicase continues unwinding the DNA, single stranded DNA (ssDNA) accumulates and becomes coated with RPA, which then initiates signals via PCNA mono-ubiquitination to activate trans-lesion polymerases and via ATR and Chk1 to inhibit Cdk2-dependent cell cycle progression. The signals are conveyed by way of a complex network of molecular interactions. To clarify those complexities, we have constructed a molecular interaction map (MIM) using a novel hierarchical assembly procedure. Molecules were arranged on the map in hierarchical levels according to interaction step distance from the DNA region of stalled replication. The hierarchical MIM allows us to disentangle the network's interlocking pathways and loops and to suggest functionally significant features of network architecture. The MIM shows how parallel pathways and multiple feedback loops can provide failsafe and robust switch-like responses to replication stress. Within the central level of hierarchy ATR and Claspin together appear to function as a nexus that conveys signals from many sources to many destinations. We noted a division of labor between those two molecules, separating enzymatic and structural roles. In addition, the network architecture disclosed by the hierarchical map, suggested a speculative model for how molecular crowding and the granular localization of network components in the cell nucleus can facilitate function.

Previous studies in rapidly proliferating rodent cells have suggested that the lethal effect of the DNA topoisomerase I inhibitor, camptothecin (CPT) is dependent upon the active participation of DNA replication (Holm et al. Cancer Res. 49:6365-6368; 1989). The purpose of the current study was to determine if this relationship applies to more slowly growing human cells. In our present study, we employed the human colon carcinoma cell line, HT-29 (45 hr doubling time). Flow cytometric determination of S-phase cells either by S-phase fit model or rectangle fit model analysis predicted that 21% of exponentially growing HT-29 cells were undergoing DNA replication. These findings were confirmed by immunofluorescence microscopy of bromodeoxyuridine labeled cells. Based on these findings, we would have expected only 20-30% of the cells to be susceptible to brief treatment (30 min) with CPT. Instead, 90-95% of HT-29 cells were killed. This apparent disparity was not due to prolonged cellular retention of drug after treatment because protein-linked DNA strand breaks reversed within 15 min of drug removal. Moreover, the DNA replication inhibitor, aphidicolin, fully protected HT-29 cells against CPT-induced killing but did not affect the production of CPT-induced protein-linked DNA strand breaks. Similar results were also obtained using the CPT-analog, 10,11-methylenedioxy-camptothecin, which was 5- to 10-fold more potent than camptothecin (O'Connor et al. Cancer Commun. 2:395-400; 1990).(ABSTRACT TRUNCATED AT 250 WORDS)

DNA filter elution assays have proved useful in studies of DNA strand breaks and crosslinks produced in mammalian cells or tissues by a wide variety of carcinogenic and cancer chemotherapeutic agents. The basic types of DNA lesions that can be measured include single and double-strand breaks, interstrand crosslinks and DNA-protein crosslinks. DNA filter elution has also been adapted to the assay of other lesions such as alkali-labile sites and the protein-associated strand breaks of topoisomerase DNA cleavage complexes. The essential concepts and theory of the technique are discussed and the applications of the technique to various types of studies are critically reviewed.

TRAIL is an endogenous death receptor ligand also used therapeutically because of its selective proapoptotic activity in cancer cells. In the present study, we examined chromatin alterations induced by TRAIL and show that TRAIL induces a rapid activation of DNA damage response (DDR) pathways with histone H2AX, Chk2, ATM, and DNA-PK phosphorylations. Within 1 h of TRAIL exposure, immunofluorescence confocal microscopy revealed gamma-H2AX peripheral nuclear staining (gamma-H2AX ring) colocalizing with phosphorylated/activated Chk2, ATM, and DNA-PK inside heterochromatin regions. The marginal distribution of DDR proteins in early apoptotic cells is remarkably different from the focal staining seen after DNA damage. TRAIL-induced DDR was suppressed upon caspase inhibition or Bax inactivation, demonstrating that the DDR activated by TRAIL is downstream from the mitochondrial death pathway. H2AX phosphorylation was dependent on DNA-PK, while Chk2 phosphorylation was dependent on both ATM and DNA-PK. Downregulation of Chk2 decreased TRAIL-induced cell detachment; delayed the activation of caspases 2, 3, 8, and 9; and reduced TRAIL-induced cell killing. Together, our findings suggest that nuclear activation of Chk2 by TRAIL acts as a positive feedback loop involving the mitochondrion-dependent activation of caspases, independently of p53.

We have assessed the ability of nucleosomes to influence the formation of mammalian topoisomerase II-DNA complexes by mapping the sites of cleavage induced by four unrelated topoisomerase II inhibitors in naked versus nucleosome-reconstituted SV40 DNA. DNA fragments were reconstituted with histone octamers from HeLa cells by the histone exchange method. Nucleosome positions were determined by comparing micrococcal nuclease cleavage patterns of nucleosome-reconstituted and naked DNA. Three types of DNA regions were defined: 1) regions with fixed nucleosome positioning; 2) regions lacking regular nucleosome phasing; and 3) a region around the replication origin (from position 5100 to 600) with no detectable nucleosomes. Topoisomerase II cleavage sites were suppressed in nucleosomes and persisted or were enhanced in linker DNA and in the nucleosome-free region around the replication origin. Incubation of reconstituted chromatin with topoisomerase II protected nucleosome-free regions from micrococcal nuclease cleavage without changing the overall micrococcal nuclease cleavage pattern. Thus, the present results indicate that topoisomerase II binds preferentially to nucleosome-free DNA and that the presence of nucleosomes at preferred DNA sequences influences drug-induced DNA breaks by topoisomerase II inhibitors.

The circadian clock is a set of regulatory steps that oscillate with a period of approximately 24 hours influencing many biological processes. These oscillations are robust to external stresses, and in the case of genotoxic stress (i.e. DNA damage), the circadian clock responds through phase shifting with primarily phase advancements. The effect of DNA damage on the circadian clock and the mechanism through which this effect operates remains to be thoroughly investigated. Here we build an in silico model to examine damage-induced circadian phase shifts by investigating a possible mechanism linking circadian rhythms to metabolism. The proposed model involves two DNA damage response proteins, SIRT1 and PARP1, that are each consumers of nicotinamide adenine dinucleotide (NAD), a metabolite involved in oxidation-reduction reactions and in ATP synthesis. This model builds on two key findings: 1) that SIRT1 (a protein deacetylase) is involved in both the positive (i.e. transcriptional activation) and negative (i.e. transcriptional repression) arms of the circadian regulation and 2) that PARP1 is a major consumer of NAD during the DNA damage response. In our simulations, we observe that increased PARP1 activity may be able to trigger SIRT1-induced circadian phase advancements by decreasing SIRT1 activity through competition for NAD supplies. We show how this competitive inhibition may operate through protein acetylation in conjunction with phosphorylation, consistent with reported observations. These findings suggest a possible mechanism through which multiple perturbations, each dominant during different points of the circadian cycle, may result in the phase advancement of the circadian clock seen during DNA damage., Author Summary Many physiological processes are regulated by the circadian clock, and we are continuing to learn about the role of the circadian clock in disease. Research in recent years has begun to shed light on the feedback mechanisms that exist between circadian regulation and other processes, including metabolism and the response to DNA damage. A challenge has been to understand the dynamic nature of the protein interactions of these processes, which often involve protein modification as a means of communicating cellular states, such as damaged DNA. Here we have devised a model that simulates an alteration of the circadian clock that is observed during DNA damage response. A novel aspect of this model is the inclusion of SIRT1, a protein that regulates core circadian proteins through modification and helps to repress gene expression. SIRT1 is dependent on a metabolite regulated by the circadian clock and is depleted during DNA damage. In conjunction with a second form of protein modification, our results suggest that multiple forms of protein modification may contribute to the experimentally observed alterations to circadian function.

Stability of the 'guardian of the genome' tumor suppressor protein p53 is regulated predominantly through its ubiquitination. The ubiquitin-specific protease HAUSP plays an important role in this process. Recent experiments showed that p53 demonstrates a differential response to changes in HAUSP which nature and significance are not understood yet. Here a data-driven mathematical model of the Mdm2-mediated p53 ubiquitination network is presented which offers an explanation for the cause of such a response. The model predicts existence of the HAUSP-regulated switch from auto- to p53 ubiquitination by Mdm2. This switch suggests a potential role of HAUSP as a downstream target of stress signals in cells. The model accounts for a significant amount of experimental data, makes predictions for some rate constants, and can serve as a building block for the larger model describing a complex dynamic response of p53 to cellular stresses.

To organize the rapidly accumulating information on bioregulatory networks related to the histone γ-H2AX-ATM-Chk2-p53-Mdm2 pathways in concise and unambiguous diagrams, we used the molecular interaction map notation (http://discover.nci.nih.gov/min). Molecular interaction maps are particularly useful for networks that include protein-protein binding and posttranslational modifications (e.g., phosphorylation). Both are important for nearly all of the proteins involved in DNA double-strand break signaling. Visualizing the regulatory circuits underlying cellular signaling may help identify key regulatory reactions and defects that can serve as targets for anticancer drugs.

Hereditary spinocerebellar ataxia with axonal neuropathy (SCAN1) is caused by an inactivating mutation (H493R) in the enzyme tyrosyl-DNA phosphodiesterase (Tdp1), which removes blocked 3'-termini at DNA strand breaks. Using SCAN1 cells treated with the specific topoisomerase I (Top1) inhibitor camptothecin, we find enhanced levels of Top1 cleavage complexes (Top1cc) and defective reversal of Top1cc in SCAN1 Tdp1-deficient cells, indicating a direct involvement of Tdp1 in the repair of Top1cc. Because the defective removal of Top1cc and the hypersensitivity of SCAN1 cells to camptothecin are not affected by aphidicolin, we propose that Tdp1 is involved in the repair of Top1cc associated with transcription damage in SCAN1 cells.

To help us understand how bioregulatory networks operate, we need a standard notation for diagrams analogous to electronic circuit diagrams. Such diagrams must surmount the difficulties posed by complex patterns of protein modifications and multiprotein complexes. To meet that challenge, we have designed the molecular interaction map (MIM) notation (http://discover.nci.nih.gov/mim/). Here we show the advantages of the MIM notation for three important types of diagrams: (1) explicit diagrams that define specific pathway models for computer simulation; (2) heuristic maps that organize the available information about molecular interactions and encompass the possible processes or pathways; and (3) diagrams of combinatorially complex models. We focus on signaling from the epidermal growth factor receptor family (EGFR, ErbB), a network that reflects the major challenges of representing in a compact manner the combinatorial complexity of multimolecular complexes. By comparing MIMs with other diagrams of this network that have recently been published, we show the utility of the MIM notation. These comparisons may help cell and systems biologists adopt a graphical language that is unambiguous and generally understood.

Purpose: The PathVisio-Faceted Search plugin helps users explore and understand complex pathways by overlaying experimental data and data from webservices, such as Ensembl BioMart, onto diagrams drawn using formalized notations in PathVisio. The plugin then provides a filtering mechanism, known as a faceted search, to find and highlight diagram nodes (e.g. genes and proteins) of interest based on imported data. The tool additionally provides a flexible scripting mechanism to handle complex queries. Availability: The PathVisio-Faceted Search plugin is compatible with PathVisio 3.0 and above. PathVisio is compatible with Windows, Mac OS X and Linux. The plugin, documentation, example diagrams and Groovy scripts are available at http://PathVisio.org/wiki/PathVisioFacetedSearchHelp. The plugin is free, open-source and licensed by the Apache 2.0 License. Contact:augustin@mail.nih.gov or jakeyfried@gmail.com

A complex bioregulatory network could be more easily comprehended if its essential function could be described by a small “core” subsystem, and if its response characteristics were switch-like. We tested this proposition by simulation studies of the hypoxia response control network. We hypothesized that a small subsystem governs the basics of the cellular response to hypoxia and that this response has a sharp oxygen-dependent transition. A molecular interaction map of the network was prepared, and an evolutionarily conserved core subsystem was extracted that could control the activity of hypoxia response promoter elements on the basis of oxygen concentration. The core subsystem included the hypoxia-inducible transcription factor (HIFα:ARNT heterodimer), proline hydroxylase, and the von Hippel-Lindau protein. Simulation studies showed that the same core subsystem can exhibit switch-like responses both to oxygen level and to HIFα synthesis rate, thus suggesting a mechanism for hypoxia response promoter element-dependent responses common to both hypoxia and growth factor signaling. The studies disclosed the mechanism responsible for the sharp transitions. We show how parameter sets giving switch-like behavior can be found and how this type of behavior provides a foundation for quantitative studies in cells.

Topoisomerase I (Top1) catalyzes two transesterification reactions: single-strand DNA cleavage and religation that are normally coupled for the relaxation of DNA supercoiling in transcribing and replicating chromatin. A variety of endogenous DNA modifications, potent anticancer drugs and carcinogens uncouple these two reactions, resulting in the accumulation of Top1 cleavage complexes. Top1 cleavage complexes damage DNA and kill cells by generating replication-mediated DNA double-strand breaks (DSBs) and by stalling transcription complexes. The repair of Top1-mediated DNA lesions involves integrated pathways that are conserved from yeasts to humans. Top1-mediated DNA damage and cell cycle checkpoint responses can be studied biochemically and genetically in yeast and human cells with known genetic defects. Defects in these repair/checkpoint pathways, which promote tumor development, explain, at least in part, the selectivity of camptothecins and other Top1 inhibitors for cancer cells.

Camptothecin (CPT) and its derivatives target mammalian DNA topoisomerase I (top1) and are among the most effective novel anticancer drugs. However, the activity of CPTs is limited by several factors, including drug inactivation by lactone ring opening, tumor drug resistance, and toxicity in patients. Novel top1 inhibitors are being searched to overcome such limitations and expand the anticancer spectrum of camptothecins. MJ-III-65 (NSC 706744) is among the most promising indenoisoquinolines to date. In this study, we show that MJ-III-65 enhances top1 cleavage complexes by both inhibiting their reversal (religation) more efficiently than CPT and by enhancing their formation. The top1 DNA cleavage complexes induced by MJ-III-65 exhibit a different distribution pattern compared with CPT and exhibit different base sequence preferences immediately around the top1 cleavage sites. Although CPTs have a preference for thymine at the (-1) position and guanine at the (+1) position of the top1-mediated DNA cleavage sites, MJ-III-65 can accommodate different base pairs at the (-1), (+1), or (+2) position, with a preference for a cytosine at the (-1) position on the scissile strand. Another difference with CPTs is the activity of MJ-III-65 against CPT-resistant top1 enzymes, implying that the amino acid residue interactions with top1 are different for MJ-III-65 and CPTs. As with CPT, MJ-III-65 is inactive against vaccinia top1. This study shows the specific molecular interactions of MJ-III-65 with top1 and demonstrates that MJ-III-65 is a potentially useful top1 inhibitor that enhances and traps top1 cleavage sites not sensitive to CPTs.

To uncover transcriptional stress responses related to p53, we used cDNA microarrays (National Cancer Institute Oncochips comprising 6500 different genes) to characterize the gene expression profiles of wild-type p53 HCT-116 cells and an isogenic p53 knockout counterpart after treatment with topotecan, a specific topoisomerase I inhibitor. The use of the p53 knockout cells had the advantage over p53-overexpressing systems in that p53 activation is mediated physiologically. RNA was extracted after low (0.1 microM)- and high (1 microM)-dose topotecan at multiple time points within the first 6 h of treatment. To facilitate simultaneous study of the p53 status and pharmacological effects on gene expression, we developed a novel "cross-referenced network" experimental design and used multiple linear least squares fitting to optimize estimates of relative transcript levels in the network of experimental conditions. Approximately 10% of the transcripts were up- or down-regulated in response to topotecan in the p53+/+ cells, whereas only 1% of the transcripts changed in the p53-/- cells, indicating that p53 has a broad effect on the transcriptional response to this stress. Individual transcripts and their relationships were analyzed using clustered image maps and by a novel two-dimensional analysis/visualization, gene expression map, in which each gene expression level is represented as a function of both the genotypic/phenotypic difference (i.e., p53 status) and the treatment effect (i.e., of topotecan dose and time of exposure). Overall, drug-induced p53 activation was associated with a coherent genetic program leading to cell cycle arrest and apoptosis. We identified novel p53-induced and DNA damage-induced genes (the proapoptotic SIVA gene and a set of transforming growth factor beta-related genes). Genes induced independently of p53 included the antiapoptotic cFLIP gene and known stress genes related to the mitogen-activated protein kinase pathway and the Fos/Jun pathway. Genes that were negatively regulated by p53 included members of the antiapoptotic protein chaperone heat shock protein 70 family. Finally, among the p53-dependent genes whose expression was independent of drug treatment was S100A4, a small Ca(2+)-binding protein that has recently been implicated in p53 binding and regulation. The new experimental design and gene expression map analysis introduced here are applicable to a wide range of studies that encompass both treatment effects and genotypic or phenotypic differences.

To study the molecular mechanisms by which drug resistance develops, we compared DU145 humanprostate cancer cells with a subline selected for resistance to camptothecin. Differences in gene expression level were assessed by hybridizing the two cell types against each other using quadruplicate "Oncochip" cDNA microarrays that included 1648 cancer-related genes. Expression levels differing by a factor of1.5 were detected for 181 of the genes. These differences were judged statistically reliable on the basis of a stratum-adjusted Kruskal-Wallis test, after taking into account a dye-dependent variable. The 181 expression-altered genes included a larger than expected number of the "apoptosis-related" genes (P = 0.04). To assess whether this observation reflected a generalized resistance of RCO.1 to apoptosis, we exposed the cells to a range of stresses (cisplatin, staurosporine, UV, ionizing radiation, and serum starvation) and found greatly reduced apoptotic responses for RC0.1 (relative to DU145) using flow cytometric Annexin V and terminal deoxynucleotidyl transferase-mediated nick end labeling assays. We next examined the apoptosis-related genes in the context of a molecular interaction map and found expression differences in the direction "expected" on the basis of the apoptosis-resistance of RC0.1 for BAD, caspase-6, and genes that signal via the Akt pathway. Exposure of the cells to wortmannin, an inhibitor of the Akt effector phosphatidylinositol 3-kinase, provided functional support for involvement of the Akt pathway. However, closer examination of the molecular interaction map revealed a paradox: many of the expression differences observed for apoptosis-related genes were in the direction "contrary" to that expected given the resistance of RC0.1. The map indicated that most of these unexpected expression differences were associated with genes involved in the nuclear factor kappa B and transforming growth factor beta pathways. Overall, the patterns that emerged suggested a two-step model for the selection process that led to resistance in RC0.1 cells. The first hypothesized step would involve a decrease in apoptotic susceptibility through changes in the apoptosis-control machinery associated with the Bcl-2 and caspase gene families, and also in antiapoptotic pathways operating through Akt/PKB. The second step would involve changes in multifunctional upstream genes (including some genes in the nuclear factor kappa B and transforming growth factor beta pathways) that can facilitate apoptosis but that would also tend to contribute to cell proliferation in the presence of drug. Thus, we propose that a downstream blockade of apoptosis was "permissive" for the selection of upstream pathway changes that would otherwise have induced apoptosis. This model is analogous to one suggested previously for the relationship between oncogene function and apoptosis in carcinogenesis.

UCN-01 (7-hydroxystaurosporine) is a cell-cycle checkpoint abrogator that sensitizes cells to ionizing radiation (IR) and chemotherapeutic agents. It has been shown previously that UCN-01 abrogates DNA-damage-induced G(2) checkpoint most selectively in p53-defective cells, by primarily targeting Chk1. Here we show that UCN-01 prevented IR-induced p53 up-regulation and p53 phosphorylation on serine 20, a site previously identified for Chk2 (or/and Chk1) kinase. We found that in human colon carcinoma HCT116 cells, IR treatment enhanced Chk2 kinase activity, whereas Chk1 activity remained unchanged, which suggested that UCN-01 may interrupt IR-induced p53 response by inhibiting Chk2 kinase. This conclusion is supported by in vitro kinase assays, showing that UCN-01 inhibits Chk2 immunoprecipitated from HCT116 cells (IC(50), approximately 10 nM). In addition, UCN-01 efficiently abrogated both the initiation and maintenance of IR-induced G(2) arrest in HCT116 cells and their isogenic p53 (-/-) derivative, indicating that G(2) checkpoint abrogation by UCN-01 is p53 independent. In the p53 (-/-) cells, there was no p21(Waf1/Cip1) induction nor UCN-01-induced apoptosis. Taken together, these observations indicate that UCN-01 can modulate both Chk1 and Chk2 in intact cells and enhance IR-induced apoptosis in p53-deficient, and consequently p21-deficient, cells.

cDNA microarray technology can be used to establish associations between characteristic gene expression patterns and molecular responses to drug therapy. In this study, we used cDNA microarrays of 1694 cancer-related genes to monitor the gene expression consequences of the treatment of HCT116 colon cancer cells with the topoisomerase I inhibitor camptothecin (CPT). To obtain a more homogeneous cellular response, we synchronized the cells in S-phase using aphidicolin (APH) before CPT treatment. Brief incubation with 20 and 1000 nM CPT caused reversible and irreversible G(2) arrest, respectively, and the patterns of gene expression change (with reference to untreated controls) were strikingly different at the two concentrations. Thirty-three genes, mainly divided into three groups, showed characteristic changes in the first 20 h as a consequence of treatment. Northern blots performed for five of these genes (each under eight experimental conditions) were quite consistent with the microarray results (average correlation coefficient, 0.86). Several p53-activated stress response genes were up-regulated after treatment with 1000 nM CPT or prolonged exposure to APH, but it seemed that the up-regulation did not directly cause cell cycle arrest because the up-regulation induced by prolonged treatment with APH did not prevent cell cycle progression after removal of APH. In contrast, cell cycle-dependent up-regulation of a group of mitosis-related genes was delayed or blocked after CPT treatments. The interrupted up-regulation of this group of genes was directly associated with G(2) arrest. In addition, we observed down-regulation of gene expression in cells that were recovering from cell cycle delay. The observations reported here suggest a fundamental difference at the gene expression level between the molecular mechanism of reversible G(2) delay that follows mild DNA damage and the mechanism of permanent G(2) arrest that follows more extensive DNA damage.

Exome sequencing provides unprecedented insights into cancer biology and pharmacological response. Here we assess these two parameters for the NCI-60, which is among the richest genomic and pharmacological publicly available cancer cell line databases. Homozygous genetic variants that putatively affect protein function were identified in 1,199 genes (approximately 6% of all genes). Variants that are either enriched or depleted compared to non-cancerous genomes, and thus may be influential in cancer progression and differential drug response were identified for 2,546 genes. Potential gene knockouts are made available. Assessment of cell line response to 19,940 compounds, including 110 FDA-approved drugs, reveals ≈80-fold range in resistance versus sensitivity response across cell lines. 103,422 gene variants were significantly correlated with at least one compound (at p

Using gene expression data to enhance our knowledge of control networks relevant to cancer biology and therapy is a challenging but urgent task. Based on the premise that genes that are expressed together in a variety of cell types are likely to functions together, we derived mutually correlated genes that function together in various processes in epithelial-like tumor cells. Expression-correlated genes were derived from data for the NCI-60 human tumor cell lines, as well as data from the Broad Institute's CCLE cell lines. NCI-60 cell lines that selectively expressed a mutually correlated subset of tight junction genes served as a signature for epithelial-like cancer cells. Those signature cell lines served as a seed to derive other correlated genes, many of which had various other epithelial-related functions. Literature survey yielded molecular interaction and function information about those genes, from which molecular interaction maps were assembled. Many of the genes had epithelial functions unrelated to tight junctions, demonstrating that new function categories were elicited. The most highly correlated genes were implicated in the following epithelial functions: interactions at tight junctions (CLDN7, CLDN4, CLDN3, MARVELD3, MARVELD2, TJP3, CGN, CRB3, LLGL2, EPCAM, LNX1); interactions at adherens junctions (CDH1, ADAP1, CAMSAP3); interactions at desmosomes (PPL, PKP3, JUP); transcription regulation of cell-cell junction complexes (GRHL1 and 2); epithelial RNA splicing regulators (ESRP1 and 2); epithelial vesicle traffic (RAB25, EPN3, GRHL2, EHF, ADAP1, MYO5B); epithelial Ca(+2) signaling (ATP2C2, S100A14, BSPRY); terminal differentiation of epithelial cells (OVOL1 and 2, ST14, PRSS8, SPINT1 and 2); maintenance of apico-basal polarity (RAB25, LLGL2, EPN3). The findings provide a foundation for future studies to elucidate the functions of regulatory networks specific to epithelial-like cancer cells and to probe for anti-cancer drug targets.

Topoisomerase 1 (top1) inhibitors are proving useful against a range of refractory tumors, and there is considerable interest in the development of additional top1 agents. Despite crystallographic studies, the binding site and ligand properties that lead to activity are poorly understood. Here we report a unique approach to quantitative structure-activity relationship (QSAR) analysis based on the National Cancer Institute's (NCI) drug databases. In 1990, the NCI established a drug discovery program in which compounds are tested for their ability to inhibit the growth of 60 different human cancer cell lines in culture. More than 70 000 compounds have been screened, and patterns of activity against the 60 cell lines have been found to encode rich information on mechanisms of drug action and drug resistance. Here, we use hierarchical clustering to define antitumor activity patterns in a data set of 167 tested camptothecins (CPTs) in the NCI drug database. The average pairwise Pearson correlation coefficient between activity patterns for the CPT set was 0.70. Coherence between chemical structures and their activity patterns was observed. QSAR studies were carried out using the mean 50% growth inhibitory concentrations (GI(50)) for 60 cell lines as the dependent variables. Different statistical methods, including stepwise linear regression, principal component regression (PCR), partial least-squares regression (PLS), and fully cross-validated genetic function approximation (GFA) were applied to construct quantitative structure-antitumor relationship models. For our data set, the GFA method performed better in terms of correlation coefficients and cross-validation analysis. A number of molecular descriptors were identified as being correlated with antitumor activity. Included were partial atomic charges and three interatomic distances that define the relative spatial dispositions of three significant atoms (the hydroxyl hydrogen of the E-ring, the lactone carbonyl oxygen of the E-ring, and the carbonyl oxygen of the D-ring). The cross-validated r(2) for the final GFA model was 0.783, indicating a predictive QSAR model.