SCLC-CellMiner: A Resource for Small Cell Lung Cancer Cell Line Genomics and Pharmacology Based on Genomic Signatures

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.

CellMinerCDB for Integrative Cross-Database Genomics and Pharmacogenomics Analyses of Cancer Cell Lines

CellMinerCDB provides a web-based resource (https://discover.nci.nih.gov/cellminercdb/) for integrating multiple forms of pharmacological and genomic analyses, and unifying the richest cancer cell line datasets (the NCI-60, NCI-SCLC, Sanger/MGH GDSC, and Broad CCLE/CTRP). CellMinerCDB enables data queries for genomics and gene regulatory network analyses, and exploration of pharmacogenomic determinants and drug signatures. It leverages overlaps of cell lines and drugs across databases to examine reproducibility and expand pathway analyses. We illustrate the value of CellMinerCDB for elucidating gene expression determinants, such as DNA methylation and copy number variations, and highlight complexities in assessing mutational burden. We demonstrate the value of CellMinerCDB in selecting drugs with reproducible activity, expand on the dominant role of SLFN11 for drug response, and present novel response determinants and genomic signatures for topoisomerase inhibitors and schweinfurthins. We also introduce LIX1L as a gene associated with mesenchymal signature and regulation of cellular migration and invasiveness.

The NCI-60 Methylome and Its Integration into CellMiner

A unique resource for systems pharmacology and genomic studies is the NCI-60 cancer cell line panel, which provides data for the largest publicly available library of compounds with cytotoxic activity (∼21,000 compounds), including 108 FDA-approved and 70 clinical trial drugs as well as genomic data, including whole-exome sequencing, gene and miRNA transcripts, DNA copy number, and protein levels. Here, we provide the first readily usable genome-wide DNA methylation database for the NCI-60, including 485,577 probes from the Infinium HumanMethylation450k BeadChip array, which yielded DNA methylation signatures for 17,559 genes integrated into our open access CellMiner version 2.0 (https://discover.nci.nih.gov/cellminer). Among new insights, transcript versus DNA methylation correlations revealed the epithelial/mesenchymal gene functional category as being influenced most heavily by methylation. DNA methylation and copy number integration with transcript levels yielded an assessment of their relative influence for 15,798 genes, including tumor suppressor, mitochondrial, and mismatch repair genes. Four forms of molecular data were combined, providing rationale for microsatellite instability for 8 of the 9 cell lines in which it occurred. Individual cell line analyses showed global methylome patterns with overall methylation levels ranging from 17% to 84%. A six-gene model, including PARP1, EP300, KDM5C, SMARCB1, and UHRF1 matched this pattern. In addition, promoter methylation of two translationally relevant genes, Schlafen 11 (SLFN11) and methylguanine methyltransferase (MGMT), served as indicators of therapeutic resistance or susceptibility, respectively. Overall, our database provides a resource of pharmacologic data that can reinforce known therapeutic strategies and identify novel drugs and drug targets across multiple cancer types. Cancer Res; 77(3); 601-12. ©2016 AACR.

Predicted Role of NAD Utilization in the Control of Circadian Rhythms during DNA Damage Response

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.

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.

Using drug response data to identify molecular effectors, and molecular "omic" data to identify candidate drugs in cancer

The current convergence of molecular and pharmacological data provides unprecedented opportunities to gain insights into the relationships between the two types of data. Multiple forms of large-scale molecular data, including but not limited to gene and microRNA transcript expression, DNA somatic and germline variations from next-generation DNA and RNA sequencing, and DNA copy number from array comparative genomic hybridization are all potentially informative when one attempts to recognize the panoply of potentially influential events both for cancer progression and therapeutic outcome. Concurrently, there has also been a substantial expansion of the pharmacological data being accrued in a systematic fashion. For cancer cell lines, the National Cancer Institute cell line panel (NCI-60), the Cancer Cell Line Encyclopedia (CCLE), and the collaborative Genomics of Drug Sensitivity in Cancer (GDSC) databases all provide subsets of these forms of data. For the patient-derived data, The Cancer Genome Atlas (TCGA) provides analogous forms of genomic information along with treatment histories. Integration of these data in turn relies on the fields of statistics and statistical learning. Multiple algorithmic approaches may be chosen, depending on the data being considered, and the nature of the question being asked. Combining these algorithms with prior biological knowledge, the results of molecular biological studies, and the consideration of genes as pathways or functional groups provides both the challenge and the potential of the field. The ultimate goal is to provide a paradigm shift in the way that drugs are selected to provide a more targeted and efficacious outcome for the patient.

Gene expression correlations in human cancer cell lines define molecular interaction networks for epithelial phenotype

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.

SIRT1/PARP1 crosstalk: connecting DNA damage and metabolism

An intricate network regulates the activities of SIRT1 and PARP1 proteins and continues to be uncovered. Both SIRT1 and PARP1 share a common co-factor nicotinamide adenine dinucleotide (NAD+) and several common substrates, including regulators of DNA damage response and circadian rhythms. We review this complex network using an interactive Molecular Interaction Map (MIM) to explore the interplay between these two proteins. Here we discuss how NAD + competition and post-transcriptional/translational feedback mechanisms create a regulatory network sensitive to environmental cues, such as genotoxic stress and metabolic states, and examine the role of those interactions in DNA repair and ultimately, cell fate decisions.

Transcription poisoning by Topoisomerase I is controlled by gene length, splice sites, and miR-142-3p

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.

Abstract 5219: Molecular phenotype of an epithelial-like subset of the NCI-60 human tumor cell lines and relevance to gene expression patterns in TCGA normal breast and breast cancer tissue samples

Better cancer therapy is likely to come from molecular characterization of individual patients’ tumors that will permit selection and development of tailored therapies. To that end, many bioinformatic studies have been carried out on human tumor cell lines and tissue samples. Large gene expression databases have become available for the National Cancer institute's 60 human tumor cell lines (NCI-60) and for normal breast and breast cancer tissue samples from the Cancer Genome Atlas (TCGA) Research Network. We aimed to infer functional relationships among expression-correlated genes based on these data and on published information about molecular interactions. In addition, we asked how or to what extent inferences based on cell line data may apply to human tissue samples.

Among the best-defined cell phenotypes are those of epithelia. Epithelial cells have polarity that distinguishes an apical region, directed towards a surface or lumen, from a basolateral region resting on a connective tissue basement membrane. Polarity is created by transport of appropriate molecular components as cargo in vesicles that are moved by motor proteins along cytoskeletal tracks to the relevant cell regions. Epithelial polarity is maintained by molecular structures, such as tight junctions, that prevent back-flow between the apical and basolateral surfaces of the cell.

We initially defined an epithelial-like phenotype on the basis of expression of genes coding for tight junction and adherence junction proteins and their family members. Among the NCI-60, we found 11 cell lines that express a mutually correlated subset of those genes. The relevant cell lines were breast cancer MCF7 and T47D; colon cancer COLO205, HCC_2998, HCT_116, HCT_15, HT29, and KM12; lung cancer NCI_H322M; and ovary cancer OVCAR_3 and OVCAR_4; we call these the NCI-60 epithelial consensus (NEC) cell lines. The following tight junction genes were selectively expressed in the NEC cell lines: TJP3, CLDN3, CLDN4, CLDN7, OCLN, MARVELD2, and MARVELD3. Highly correlated with these genes was the adherence junction gene CDH1/E-cadherin (as expected).

We then assembled an expanded list of 75 genes that were selectively expressed in the NEC cell lines, and found that many of them function in the epithelia-specific processes summarized above.

With some notable exceptions, the tight junction and adherence junction genes co-expressed in the NEC cell lines were also co-expressed in the normal breast and breast cancer tissue samples. An unexpected finding was that some of the normal breast tissue samples expressed a different set of tight junction family genes.

This work delineates similarities and differences between epithelial-like NCI-60 cell lines and TCGA breast tissue samples in regard to the expression of functionally defined genes.

Citation Format: Kurt W. Kohn, Barry R. Zeeberg, William C. Reinhold, Ari Kahn, Yves Pommier. Molecular phenotype of an epithelial-like subset of the NCI-60 human tumor cell lines and relevance to gene expression patterns in TCGA normal breast and breast cancer tissue samples. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5219. doi:10.1158/1538-7445.AM2013-5219

PathVisio-Faceted Search: an exploration tool for multi-dimensional navigation of large pathways

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.

CellMiner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set

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.

Gene expression profiles of the NCI-60 human tumor cell lines define molecular interaction networks governing cell migration processes

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.

Functional categories associated with clusters of genes that are co-expressed across the NCI-60 cancer cell lines

Background

The NCI-60 is a panel of 60 diverse human cancer cell lines used by the U.S. National Cancer Institute to screen compounds for anticancer activity. In the current study, gene expression levels from five platforms were integrated to yield a single composite transcriptome profile. The comprehensive and reliable nature of that dataset allows us to study gene co-expression across cancer cell lines.

Methodology/Principal Findings

Hierarchical clustering revealed numerous clusters of genes in which the genes co-vary across the NCI-60. To determine functional categorization associated with each cluster, we used the Gene Ontology (GO) Consortium database and the GoMiner tool. GO maps genes to hierarchically-organized biological process categories. GoMiner can leverage GO to perform ontological analyses of gene expression studies, generating a list of significant functional categories.

Conclusions/Significance

GoMiner analysis revealed many clusters of coregulated genes that are associated with functional groupings of GO biological process categories. Notably, those categories arising from coherent co-expression groupings reflect cancer-related themes such as adhesion, cell migration, RNA splicing, immune response and signal transduction. Thus, these clusters demonstrate transcriptional coregulation of functionally-related genes.

PathVisio-Validator: a rule-based validation plugin for graphical pathway notations

PURPOSE

The PathVisio-Validator plugin aims to simplify the task of producing biological pathway diagrams that follow graphical standardized notations, such as Molecular Interaction Maps or the Systems Biology Graphical Notation. This plugin assists in the creation of pathway diagrams by ensuring correct usage of a notation, and thereby reducing ambiguity when diagrams are shared among biologists. Rulesets, needed in the validation process, can be generated for any graphical notation that a developer desires, using either Schematron or Groovy. The plugin also provides support for filtering validation results, validating on a subset of rules, and distinguishing errors and warnings.

PathVisio-MIM: PathVisio plugin for creating and editing Molecular Interaction Maps (MIMs)

MOTIVATION

A plugin for the Java-based PathVisio pathway editor has been developed to help users draw diagrams of bioregulatory networks according to the Molecular Interaction Map (MIM) notation. Together with the core PathVisio application, this plugin presents a simple to use and cross-platform application for the construction of complex MIM diagrams with the ability to annotate diagram elements with comments, literature references and links to external databases. This tool extends the capabilities of the PathVisio pathway editor by providing both MIM-specific glyphs and support for a MIM-specific markup language file format for exchange with other MIM-compatible tools and diagram validation.

AVAILABILITY

The PathVisio-MIM plugin is freely available and works with versions of PathVisio 2.0.11 and later on Windows, Mac OS X and Linux. Information about MIM notation and the MIMML format is available at http://discover.nci.nih.gov/mim. The plugin, along with diagram examples, instructions and Java source code, may be downloaded at http://discover.nci.nih.gov/mim/mim_pathvisio.html.

A formal MIM specification and tools for the common exchange of MIM diagrams: an XML-Based format, an API, and a validation method

Background

The Molecular Interaction Map (MIM) notation offers a standard set of symbols and rules on their usage for the depiction of cellular signaling network diagrams. Such diagrams are essential for disseminating biological information in a concise manner. A lack of software tools for the notation restricts wider usage of the notation. Development of software is facilitated by a more detailed specification regarding software requirements than has previously existed for the MIM notation.

Results

A formal implementation of the MIM notation was developed based on a core set of previously defined glyphs. This implementation provides a detailed specification of the properties of the elements of the MIM notation. Building upon this specification, a machine-readable format is provided as a standardized mechanism for the storage and exchange of MIM diagrams. This new format is accompanied by a Java-based application programming interface to help software developers to integrate MIM support into software projects. A validation mechanism is also provided to determine whether MIM datasets are in accordance with syntax rules provided by the new specification.

Conclusions

The work presented here provides key foundational components to promote software development for the MIM notation. These components will speed up the development of interoperable tools supporting the MIM notation and will aid in the translation of data stored in MIM diagrams to other standardized formats. Several projects utilizing this implementation of the notation are outlined herein. The MIM specification is available as an additional file to this publication. Source code, libraries, documentation, and examples are available at http://discover.nci.nih.gov/mim.

Genome-wide analysis of novel splice variants induced by topoisomerase I poisoning shows preferential occurrence in genes encoding splicing factors

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.

mRNA and microRNA expression profiles of the NCI-60 integrated with drug activities

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 approximately 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 approximately 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.

Abstract 1103: Linkage between endogenous Chk2 activation and p14ARF expression in p53-mutated cells revealed by proteomic and genomic analyses of the NCI60

Chk2 is a multi-substrate kinase involved in the DNA damage response and genomic stability. Chk2 inhibitors are currently being evaluated in cancer models, with the rationale of enhancing the activity of DNA damaging agents in p53-deficient tumors while protecting normal tissues. Part of the ATM-Chk2-p53 axis, Chk2 is overexpressed or activated in various malignant tumors. Here we determined Chk2 copy number, total and single-exon transcriptional levels, and total and activated (pT68) protein across the 60 cell lines of the NCI DTP screen. We correlated these parameters with p53 mutational status, miRNA expression (Agilent miRNA Microarray V2), and performed supervised transcript (Agilent Whole Human Genome Oligo Microarray) class comparison analysis between high and low pT68 Chk2 cell lines using BRB Array Tools 3.8. All the lines used for this purpose were known to be p53-mutated. Our phosphoproteomic analyses demonstrate that a number of cell lines have high endogenous Chk2 activation (ovarian OVCAR3 and OVCAR4, renal RFX-393, lung NCI-H322M, and EKVX, melanoma UACC-62 and −257, leukemia RPMI-8226 and breast BT-549) whereas the wild-type p53 tumors almost invariably show undetectable Chk2 activation. Expression of total Chk2 protein and mRNA were generally positively correlated to each other. Exon array analyses across the NCI-60 suggest exon 3 splicing in almost all Chk2 high expresser cell lines and premature truncation between exons 8 and 10 in the ovarian SKOV3 cells, which could be related to the low level of Chk2 expression in that cell line. One microRNA, hsa-miR-455-3p, showed a highly significant negative correlation with both Chk2 transcript levels and with total protein levels (Spearman's r =-0.48, p < 10-3 and r=-0.35, p < 0.01 respectively) and was predicted to hybridize to Chk2 mRNA by miRBase (www.mirbase.org). Among the 181 probes differentially expressed in cells with highly activated Chk2, 123 mapped to known gene IDs, and 58 were eligible for network analysis with Ingenuity Pathways Analysis 8.0 (Ingenuity Systems Inc.). Functions associated with the networks enriched in genes differentially expressed by high pT68 cells included: cell cycle, cell-to-cell signaling, cell growth and proliferation, cancer. Remarkably, the tumor suppressor CDKN2A (p14ARF) showed a 132-fold difference between the two groups considered (unpaired t test, two-tailed p < 10-7; false discovery rate < 10-7), being highly overexpressed in high pT68 Chk2 cell lines. Several EGFR-Ras-MAPK related genes were downregulated in high pT68 cell lines. Together, these data suggest that Chk2 is part of a multifaceted regulation network in cancer. Although traditionally considered as a low-penetrance tumor suppressor, it clearly plays more complex roles in tumor biology, possibly behaving as an activated oncogene in specific genomic contexts.

Note: This abstract was not presented at the AACR 101st Annual Meeting 2010 because the presenter was unable to attend.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 1103.

Predicted functions of MdmX in fine-tuning the response of p53 to DNA damage

Tumor suppressor protein p53 is regulated by two structurally homologous proteins, Mdm2 and MdmX. In contrast to Mdm2, MdmX lacks ubiquitin ligase activity. Although the essential interactions of MdmX are known, it is not clear how they function to regulate p53. The regulation of tumor suppressor p53 by Mdm2 and MdmX in response to DNA damage was investigated by mathematical modeling of a simplified network. The simplified network model was derived from a detailed molecular interaction map (MIM) that exhibited four coherent DNA damage response pathways. The results suggest that MdmX may amplify or stabilize DNA damage-induced p53 responses via non-enzymatic interactions. Transient effects of MdmX are mediated by reservoirs of p53∶MdmX and Mdm2∶MdmX heterodimers, with MdmX buffering the concentrations of p53 and/or Mdm2. A survey of kinetic parameter space disclosed regions of switch-like behavior stemming from such reservoir-based transients. During an early response to DNA damage, MdmX positively or negatively regulated p53 activity, depending on the level of Mdm2; this led to amplification of p53 activity and switch-like response. During a late response to DNA damage, MdmX could dampen oscillations of p53 activity. A possible role of MdmX may be to dampen such oscillations that otherwise could produce erratic cell behavior. Our study suggests how MdmX may participate in the response of p53 to DNA damage either by increasing dependency of p53 on Mdm2 or by dampening oscillations of p53 activity and presents a model for experimental investigation.

A Molecular Interaction Map (MIM) akin to a circuit diagram of an electric device can give a comprehensive view of cellular processes and help understand complex protein functions in cells. To this end, we generated a MIM focused on the p53-Mdm2-MdmX network proteins and performed computer simulations to help understand how Mdm2 and MdmX may regulate p53. Proper regulation of p53 is important for cell survival: elevated levels of p53 can lead to cell death, and decreased levels of p53 can lead to cancer. Mdm2 and MdmX are structurally homologous proteins that regulate p53. Mdm2 negatively regulates p53 by degradation, but MdmX regulation of p53 is not well understood. Recently, Mdm2 and MdmX have been recognized as potential cancer therapeutic targets. In an effort to better understand how MdmX can alter the p53 activity under various conditions, we used mathematical models based on the MIM network to generate hypotheses that can be tested by experiments. Our simulations suggest that MdmX may increase the dependency of p53 on Mdm2 or dampen p53 oscillations during DNA damage response.

Network architecture of signaling from uncoupled helicase-polymerase to cell cycle checkpoints and trans-lesion DNA synthesis

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.

HAUSP-regulated switch from auto- to p53 ubiquitination by Mdm2 (in silico discovery)

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.

Chromatin challenges during DNA replication: a systems representation

In a recent review, A. Groth and coworkers presented a comprehensive account of nucleosome disassembly in front of a DNA replication fork, assembly behind the replication fork, and the copying of epigenetic information onto the replicated chromatin. Understanding those processes however would be enhanced by a comprehensive graphical depiction analogous to a circuit diagram. Accordingly, we have constructed a molecular interaction map (MIM) that preserves in essentially complete detail the processes described by Groth et al. The MIM organizes and elucidates the information presented by Groth et al. on the complexities of chromatin replication, thereby providing a tool for system-level comprehension of the effects of genetic mutations, altered gene expression, and pharmacologic intervention.

Dose-response transition from cell cycle arrest to apoptosis with selective degradation of Mdm2 and p21WAF1/CIP1 in response to the novel anticancer agent, aminoflavone (NSC 686,288)

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.

Hereditary ataxia SCAN1 cells are defective for the repair of transcription-dependent topoisomerase I cleavage complexes

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.

Depicting combinatorial complexity with the molecular interaction map notation

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.

4-nitroquinoline-1-oxide induces the formation of cellular topoisomerase I-DNA cleavage complexes

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 gamma-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 gamma-H2AX by 4NQO may contribute to the cellular effects of 4NQO, including its selective activity toward RecQ helicase BLM-deficient cells.

Chk2 molecular interaction map and rationale for Chk2 inhibitors

To organize the rapidly accumulating information on bioregulatory networks related to the histone gamma-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.

Cancers as wounds that do not heal: differences and similarities between renal regeneration/repair and renal cell carcinoma

Cancers have been described as wounds that do not heal, suggesting that the two share common features. By comparing microarray data from a model of renal regeneration and repair (RRR) with reported gene expression in renal cell carcinoma (RCC), we asked whether those two processes do, in fact, share molecular features and regulatory mechanisms. The majority (77%) of the genes expressed in RRR and RCC were concordantly regulated, whereas only 23% were discordant (i.e., changed in opposite directions). The orchestrated processes of regeneration, involving cell proliferation and immune response, were reflected in the concordant genes. The discordant gene signature revealed processes (e.g., morphogenesis and glycolysis) and pathways (e.g., hypoxia-inducible factor and insulin-like growth factor-I) that reflect the intrinsic pathologic nature of RCC. This is the first study that compares gene expression patterns in RCC and RRR. It does so, in particular, with relation to the hypothesis that RCC resembles the wound healing processes seen in RRR. However, careful attention to the genes that are regulated in the discordant direction provides new insights into the critical differences between renal carcinogenesis and wound healing. The observations reported here provide a conceptual framework for further efforts to understand the biology and to develop more effective diagnostic biomarkers and therapeutic strategies for renal tumors and renal ischemia.

A novel norindenoisoquinoline structure reveals a common interfacial inhibitor paradigm for ternary trapping of the topoisomerase I-DNA covalent complex

We show that five topoisomerase I inhibitors (two indenoisoquinolines, two camptothecins, and one indolocarbazole) each intercalate between the base pairs flanking the cleavage site generated during the topoisomerase I catalytic cycle and are further stabilized by a network of hydrogen bonds with topoisomerase I. The interfacial inhibition paradigm described for topoisomerase I inhibitors can be generalized to a variety of natural products that trap macromolecular complexes as they undergo catalytic conformational changes that create hotspots for drug binding. Stabilization of such conformational states results in uncompetitive inhibition and exemplifies the relevance of screening for ligands and drugs that stabilize ("trap") these macromolecular complexes.

p21CDKN1A allows the repair of replication-mediated DNA double-strand breaks induced by topoisomerase I and is inactivated by the checkpoint kinase inhibitor 7-hydroxystaurosporine

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.

Molecular interaction maps of bioregulatory networks: a general rubric for systems biology

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.

Targeting Chk2 Kinase: Molecular Interaction Maps and Therapeutic Rationale

Most anticancer drugs presently used clinically target genomic DNA. The selectivity of these anticancer drugs for tumor tissues is probably due to tumor-specific defects suppressing cell cycle checkpoints and DNA repair, and enhancing apoptotic response in the tumor. We will review the molecular interactions within the ATM-Chk2 pathway implicating the DNA damage sensor kinases (ATM, ATR and DNA-PK), the adaptor BRCT proteins (Nbs1, Brca1, 53BP1, MDC1) and the effector kinases (Chk2, Chk1, Plk3, JNK, p38). The molecular interaction map convention (MIM) will be used for presenting this molecular network (http://discover.nci.nih.gov/mim/). A characteristic of the ATM-Chk2 pathway is its redundancy. First, ATM and Chk2 phosphorylate common substrates including p53, E2F1, BRCA1, and Chk2 itself, which suggests that Chk2 (also known as CHECK2, Cds1 in fission yeast, and Dmchk2 or Dmnk or Loki in the fruit fly) acts as a relay for ATM and/or as a salvage pathway when ATM is inactivated. Secondly, redundancy is apparent for the substrates, which can be phosphorylated/activated at similar residues by Chk2, Chk1, and the polo kinases (Plk's). Functionally, Chk2 can activate both apoptosis (via p53, E2F1 and PML) and cell cycle checkpoint (via Cdc25A and Cdc25C, p53, and BRCA1). We will review the short list of published Chk2 inhibitors. We will also propose a novel paradigm for screening interfacial inhibitors of Chk2. Chk2 inhibitors might be used to enhance the tumor selectivity of DNA targeted agents in p53-deficient tumors, and for the treatment of tumors whose growth depends on enhanced Chk2 activity.

Molecular interaction map of the p53 and Mdm2 logic elements, which control the Off-On switch of p53 in response to DNA damage

The molecular network that controls responses to genotoxic stress is centered at p53 and Mdm2. Recent findings have shown this network to be more complex than previously envisioned. Using a notation specifically designed for circuit diagram-like representations of bioregulatory networks, we have prepared an updated molecular interaction map of the immediate connections of p53 and Mdm2, which are described as logic elements of the network. We use the map as the basis for a comprehensive review of current concepts of signal processing by these logic elements (an interactive version of the maps-eMIMs can be examined at ). We also used molecular interaction maps to propose a p53 Off-On switch in response to DNA damage.

Properties of switch-like bioregulatory networks studied by simulation of the hypoxia response control system

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 (HIFalpha: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 HIFalpha 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.

Apoptosis defects and chemotherapy resistance: molecular interaction maps and networks

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.

Molecular interaction maps--a diagrammatic graphical language for bioregulatory networks

Molecular interaction maps (MIMs) use a clear, accurate, and versatile graphical language to depict complex biological processes. Here, we discuss the main features of the MIM language and its potential uses. MIMs can be used as database resources and simulation guides, and can serve to generate new hypotheses regarding the roles of specific molecules in the bioregulatory networks that control progression through the cell cycle, differentiation, and cell death.

Repair of and checkpoint response to topoisomerase I-mediated DNA damage

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.

Differential induction of topoisomerase I-DNA cleavage complexes by the indenoisoquinoline MJ-III-65 (NSC 706744) and camptothecin: base sequence analysis and activity against camptothecin-resistant topoisomerases I

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.

Apoptosis induced by topoisomerase inhibitors

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.

Impact of p53 knockout and topotecan treatment on gene expression profiles in human colon carcinoma cells: a pharmacogenomic study

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.

Apoptotic susceptibility of cancer cells selected for camptothecin resistance: gene expression profiling, functional analysis, and molecular interaction mapping

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 of >1.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.

[Cell cycle and checkpoints in oncology: new therapeutic targets]

Leland H. Hartwell, Paul M. Nurse et R. Timothy Hunt just received the Nobel price for their discovery of the molecular components of the cell cycle and cell cycle checkpoints. This review is an update of the molecular networks driving the cell cycle and its regulation, and of the importance of this knowledge for understanding the mechanisms driving oncogenesis and therapeutic developments.

UCN-01 inhibits p53 up-regulation and abrogates gamma-radiation-induced G(2)-M checkpoint independently of p53 by targeting both of the checkpoint kinases, Chk2 and Chk1

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.

Importance of the fourth alpha-helix within the CAP homology domain of type II topoisomerase for DNA cleavage site recognition and quinolone action

We report that point mutations causing alteration of the fourth alpha-helix (alpha4-helix) of the CAP homology domain of eukaryotic (Saccharomyces cerevisiae) type II topoisomerases (Ser(740)Trp, Gln(743)Pro, and Thr(744)Pro) change the selection of type II topoisomerase-mediated DNA cleavage sites promoted by Ca(2+) or produced by etoposide, the fluoroquinolone CP-115,953, or mitoxantrone. By contrast, Thr(744)Ala substitution had minimal effect on Ca(2+)- and drug-stimulated DNA cleavage sites, indicating the selectivity of single amino acid substitutions within the alpha4-helix on type II topoisomerase-mediated DNA cleavage. The equivalent mutation in the gene for Escherichia coli gyrase causing Ser(83)Trp also changed the DNA cleavage pattern generated by Ca(2+) or quinolones. Finally, Thr(744)Pro substitution in the yeast type II topoisomerase rendered the enzyme sensitive to antibacterial quinolones. This study shows that the alpha4-helix within the conserved CAP homology domain of type II topoisomerases is critical for selecting the sites of DNA cleavage. It also demonstrates that selective amino acid residues in the alpha4-helix are important in determining the activity and possibly the binding of quinolones to the topoisomerase II-DNA complexes.

Changes in deoxyribonucleic acid linking number due to treatment of mammalian cells with the intercalating agent 4'-(9-acridinylamino)methanesulfon-m-anisidide

Treatment of mammalian cells with DNA intercalating agents produces protein-associated DNA strand breaks. These breaks have been proposed to represent the action of a topoisomerase, which would alter the DNA linking number. Changes in DNA linking number in cells treated with the intercalating agent 4'-(9-acridinylamino)methanesulfon-m-anisidide (m-AMSA) were studied by ethidium titration of nucleoid sedimentation. m-AMSA treatment was found to produce an increase in DNA linking number. Previously, we had proposed that intercalator-induced protein-associated DNA breaks act to reduce DNA torsional strain that results from the intercalator-induced decrease in DNA twist. In such a model, linking number would be expected to decrease. The finding that the DNA linking number increased following m-AMSA treatment suggests that intercalators may block enzymes that normally decrease linking number. Such enzymes would have DNA gyrase like properties. Consistent with this possibility, a DNA gyrase inhibitor, novobiocin, inhibited the restoration of normal linking number and, to a lesser degree, the reversal of protein-associated strand breaks after removal of intercalator.

Transcriptional regulation of mitotic genes by camptothecin-induced DNA damage: microarray analysis of dose- and time-dependent effects

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.

Quantitative structure-antitumor activity relationships of camptothecin analogues: cluster analysis and genetic algorithm-based studies

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.