Targeting changes in cancer: assessing pathway stability by comparing pathway gene expression coherence levels in tumor and normal tissues

Computational algorithms for in silico profiling of activating mutations in cancer

Methods to catalog and computationally assess the mutational landscape of proteins in human cancers are desirable. One approach is to adapt evolutionary or data-driven methods developed for predicting whether a single-nucleotide polymorphism (SNP) is deleterious to protein structure and function. In cases where understanding the mechanism of protein activation and regulation is desired, an alternative approach is to employ structure-based computational approaches to predict the effects of point mutations. Through a case study of mutations in kinase domains of three proteins, namely, the anaplastic lymphoma kinase (ALK) in pediatric neuroblastoma patients, serine/threonine-protein kinase B-Raf (BRAF) in melanoma patients, and erythroblastic oncogene B 2 (ErbB2 or HER2) in breast cancer patients, we compare the two approaches above. We find that the structure-based method is most appropriate for developing a binary classification of several different mutations, especially infrequently occurring ones, concerning the activation status of the given target protein. This approach is especially useful if the effects of mutations on the interactions of inhibitors with the target proteins are being sought. However, many patients will present with mutations spread across different target proteins, making structure-based models computationally demanding to implement and execute. In this situation, data-driven methods-including those based on machine learning techniques and evolutionary methods-are most appropriate for recognizing and illuminate mutational patterns. We show, however, that, in the present status of the field, the two methods have very different accuracies and confidence values, and hence, the optimal choice of their deployment is context-dependent.

Blood vessel tortuosity selects against evolution of aggressive tumor cells in confined tissue environments: A modeling approach

Cancer is a disease of cellular regulation, often initiated by genetic mutation within cells, and leading to a heterogeneous cell population within tissues. In the competition for nutrients and growth space within the tumors the phenotype of each cell determines its success. Selection in this process is imposed by both the microenvironment (neighboring cells, extracellular matrix, and diffusing substances), and the whole of the organism through for example the blood supply. In this view, the development of tumor cells is in close interaction with their increasingly changing environment: the more cells can change, the more their environment will change. Furthermore, instabilities are also introduced on the organism level: blood supply can be blocked by increased tissue pressure or the tortuosity of the tumor-neovascular vessels. This coupling between cell, microenvironment, and organism results in behavior that is hard to predict. Here we introduce a cell-based computational model to study the effect of blood flow obstruction on the micro-evolution of cells within a cancerous tissue. We demonstrate that stages of tumor development emerge naturally, without the need for sequential mutation of specific genes. Secondly, we show that instabilities in blood supply can impact the overall development of tumors and lead to the extinction of the dominant aggressive phenotype, showing a clear distinction between the fitness at the cell level and survival of the population. This provides new insights into potential side effects of recent tumor vasculature normalization approaches.

From molecular interaction to acute promyelocytic leukemia: Calculating leukemogenesis and remission from endogenous molecular-cellular network

Acute promyelocytic leukemia (APL) remains the best example of a malignancy that can be cured clinically by differentiation therapy. We demonstrate that APL may emerge from a dynamical endogenous molecular-cellular network obtained from normal, non-cancerous molecular interactions such as signal transduction and translational regulation under physiological conditions. This unifying framework, which reproduces APL, normal progenitor, and differentiated granulocytic phenotypes as different robust states from the network dynamics, has the advantage to study transition between these states, i.e. critical drivers for leukemogenesis and targets for differentiation. The simulation results quantitatively reproduce microarray profiles of NB4 and HL60 cell lines in response to treatment and normal neutrophil differentiation, and lead to new findings such as biomarkers for APL and additional molecular targets for arsenic trioxide therapy. The modeling shows APL and normal states mutually suppress each other, both in "wiring" and in dynamical cooperation. Leukemogenesis and recovery under treatment may be a consequence of spontaneous or induced transitions between robust states, through "passes" or "dragging" by drug effects. Our approach rationalizes leukemic complexity and constructs a platform towards extending differentiation therapy by performing "dry" molecular biology experiments.

Systematic analysis of time-series gene expression data on tumor cell-selective apoptotic responses to HDAC inhibitors

SAHA (suberoylanilide hydroxamic acid or vorinostat) is the first nonselective histone deacetylase (HDAC) inhibitor approved by the US Food and Drug Administration (FDA). SAHA affects histone acetylation in chromatin and a variety of nonhistone substrates, thus influencing many cellular processes. In particularly, SAHA induces selective apoptosis of tumor cells, although the mechanism is not well understood. A series of microarray experiments was recently conducted to investigate tumor cell-selective proapoptotic transcriptional responses induced by SAHA. Based on that gene expression time series, we propose a novel framework for detailed analysis of the mechanism of tumor cell apoptosis selectively induced by SAHA. Our analyses indicated that SAHA selectively disrupted the DNA damage response, cell cycle, p53 expression, and mitochondrial integrity of tumor samples to induce selective tumor cell apoptosis. Our results suggest a possible regulation network. Our research extends the existing research.

A small molecule (pluripotin) as a tool for studying cancer stem cell biology: proof of concept

BACKGROUND

Cancer stem cells (CSC) are thought to be responsible for tumor maintenance and heterogeneity. Bona fide CSC purified from tumor biopsies are limited in supply and this hampers study of CSC biology. Furthermore, purified stem-like CSC subpopulations from existing tumor lines are unstable in culture. Finding a means to overcome these technical challenges would be a useful goal. In a first effort towards this, we examined whether a chemical probe that promotes survival of murine embryonic stem cells without added exogenous factors can alter functional characteristics in extant tumor lines in a fashion consistent with a CSC phenotype.

METHODOLOGY/PRINCIPAL FINDINGS

The seven tumor lines of the NCI60 colon subpanel were exposed to SC-1 (pluripotin), a dual kinase and GTPase inhibitor that promotes self-renewal, and then examined for tumorigenicity under limiting dilution conditions and clonogenic activity in soft agar. A statistically significant increase in tumor formation following SC-1 treatment was observed (p<0.04). Cloning efficiencies and expression of putative CSC surface antigens (CD133 and CD44) were also increased. SC-1 treatment led to sphere formation in some colon tumor lines. Finally, SC-1 inhibited in vitro kinase activity of RSK2, and another RSK2 inhibitor increased colony formation implicating a role for this kinase in eliciting a CSC phenotype.

CONCLUSIONS/SIGNIFICANCE

These findings validate a proof of concept study exposure of extant tumor lines to a small molecule may provide a tractable in vitro model for understanding CSC biology.

Integrating constitutive gene expression and chemoactivity: mining the NCI60 anticancer screen

Studies into the genetic origins of tumor cell chemoactivity pose significant challenges to bioinformatic mining efforts. Connections between measures of gene expression and chemoactivity have the potential to identify clinical biomarkers of compound response, cellular pathways important to efficacy and potential toxicities; all vital to anticancer drug development. An investigation has been conducted that jointly explores tumor-cell constitutive NCI60 gene expression profiles and small-molecule NCI60 growth inhibition chemoactivity profiles, viewed from novel applications of self-organizing maps (SOMs) and pathway-centric analyses of gene expressions, to identify subsets of over- and under-expressed pathway genes that discriminate chemo-sensitive and chemo-insensitive tumor cell types. Linear Discriminant Analysis (LDA) is used to quantify the accuracy of discriminating genes to predict tumor cell chemoactivity. LDA results find 15% higher prediction accuracies, using ∼30% fewer genes, for pathway-derived discriminating genes when compared to genes derived using conventional gene expression-chemoactivity correlations. The proposed pathway-centric data mining procedure was used to derive discriminating genes for ten well-known compounds. Discriminating genes were further evaluated using gene set enrichment analysis (GSEA) to reveal a cellular genetic landscape, comprised of small numbers of key over and under expressed on- and off-target pathway genes, as important for a compound’s tumor cell chemoactivity. Literature-based validations are provided as support for chemo-important pathways derived from this procedure. Qualitatively similar results are found when using gene expression measurements derived from different microarray platforms. The data used in this analysis is available at http://pubchem.ncbi.nlm.nih.gov/andhttp://www.ncbi.nlm.nih.gov/projects/geo (GPL96, GSE32474).

Identifying responsive functional modules from protein-protein interaction network

Proteins interact with each other within a cell, and those interactions give rise to the biological function and dynamical behavior of cellular systems. Generally, the protein interactions are temporal, spatial, or condition dependent in a specific cell, where only a small part of interactions usually take place under certain conditions. Recently, although a large amount of protein interaction data have been collected by high-throughput technologies, the interactions are recorded or summarized under various or different conditions and therefore cannot be directly used to identify signaling pathways or active networks, which are believed to work in specific cells under specific conditions. However, protein interactions activated under specific conditions may give hints to the biological process underlying corresponding phenotypes. In particular, responsive functional modules consist of protein interactions activated under specific conditions can provide insight into the mechanism underlying biological systems, e.g. protein interaction subnetworks found for certain diseases rather than normal conditions may help to discover potential biomarkers. From computational viewpoint, identifying responsive functional modules can be formulated as an optimization problem. Therefore, efficient computational methods for extracting responsive functional modules are strongly demanded due to the NP-hard nature of such a combinatorial problem. In this review, we first report recent advances in development of computational methods for extracting responsive functional modules or active pathways from protein interaction network and microarray data. Then from computational aspect, we discuss remaining obstacles and perspectives for this attractive and challenging topic in the area of systems biology.

Cytometric assessment of DNA damage induced by DNA topoisomerase inhibitors

Exposure of cells to inhibitors of DNA topoisomerase I (topo I) or topoisomerase II (topo II) leads to DNA damage that often involves formation of DNA double-strand breaks (DSBs). DNA damage, particularly induction of DSBs, manifests by phosphorylation of histone H2AX on Ser-139 which is mediated by one of the protein kinases of the phosphoinositide kinase family, namely ATM, ATR, and/or DNA-PK. The presence of Ser-139 phosphorylated H2AX (gammaH2AX) is thus a reporter of DNA damage. This protocol describes quantitative assessment of gammaH2AX detected immunocytochemically in individual cells combined with quantification of cellular DNA content by cytometry. The bivariate analysis of gammaH2AX expression versus DNA content allows one to correlate DNA damage with the cell cycle phase or DNA ploidy. The protocol can also be used to assess activation (Ser-1981 phosphorylation) of ATM; this event also revealing DNA damage induced by topo I or topo II inhibitors. Examples where DNA damage was induced by topotecan (topo I) and etoposide (topo II) inhibitors are provided.

Identifying differential correlation in gene/pathway combinations

BACKGROUND

An important emerging trend in the analysis of microarray data is to incorporate known pathway information a priori. Expression level "summaries" for pathways, obtained from the expression data for the genes constituting the pathway, permit the inclusion of pathway information, reduce the high dimensionality of microarray data, and have the power to elucidate gene-interaction dependencies which are not already accounted for through known pathway identification.

RESULTS

We present a novel method for the analysis of microarray data that identifies joint differential expression in gene-pathway pairs. This method takes advantage of known gene pathway memberships to compute a summary expression level for each pathway as a whole. Correlations between the pathway expression summary and the expression levels of genes not already known to be associated with the pathway provide clues to gene interaction dependencies that are not already accounted for through known pathway identification, and statistically significant differences between gene-pathway correlations in phenotypically different cells (e.g., where the expression level of a single gene and a given pathway summary correlate strongly in normal cells but weakly in tumor cells) may indicate biologically relevant gene-pathway interactions. Here, we detail the methodology and present the results of this method applied to two gene-expression datasets, identifying gene-pathway pairs which exhibit differential joint expression by phenotype.

CONCLUSION

The method described herein provides a means by which interactions between large numbers of genes may be identified by incorporating known pathway information to reduce the dimensionality of gene interactions. The method is efficient and easily applied to data sets of ~102 arrays. Application of this method to two publicly-available cancer data sets yields suggestive and promising results. This method has the potential to complement gene-at-a-time analysis techniques for microarray analysis by indicating relationships between pathways and genes that have not previously been identified and which may play a role in disease.