Convergence of logic of cellular regulation in different premalignant cells by an information theoretic approach

Overexpressed kinetochore genes are used by cancer cells as genome destabilizers and transformation catalysts

Cancer cells have an altered transcriptome, which contributes to their abnormal behavior. Many tumors have high levels of kinetochore genes, which play important roles in genome stability. This overexpression could be utilized to destabilize cancer cell genomes, however this has not been proven specifically. We investigated the link between kinetochore gene overexpression, chromosomal number variations (CNVs) and genomic instability. Data on RNA expression and CNV from 12 different cancer types were evaluated using information theory. In all cancer types, we looked at the relationship between RNA expression and CNVs. Kinetochore gene expression was found to be substantially linked with CNV levels. In all cancer types, with the exception of thyroid cancer, highly expressed kinetochore genes were enriched in the most dominant cancer-specific co-expression subnetworks characterizing the largest patient subgroups. Except for thyroid cancer, kinetochore inner protein CENPA was among the transcripts most strongly associated with CNV values in all cancer types studied, with significantly higher expression levels in patients with high CNVs than in patients with low CNVs. CENPA function was investigated further in cell models by transfecting genomically stable (HCT116) and unstable (MCF7 and HT29) cancer cell lines using CENPA overexpression vectors. This overexpression increased the number of abnormal cell divisions in the stable cancer cell line HCT116 and, to a lesser extent, in the unstable cell lines MCF7 and HT29. Overexpression improved anchorage-independent growth properties of all cell lines. Our findings suggest that overexpression of kinetochore genes in general, and CENPA in particular, can cause genomic instability and cancer progression.

TMEA: A Thermodynamically Motivated Framework for Functional Characterization of Biological Responses to System Acclimation

The objective of gene set enrichment analysis (GSEA) in modern biological studies is to identify functional profiles in huge sets of biomolecules generated by high-throughput measurements of genes, transcripts, metabolites, and proteins. GSEA is based on a two-stage process using classical statistical analysis to score the input data and subsequent testing for overrepresentation of the enrichment score within a given functional coherent set. However, enrichment scores computed by different methods are merely statistically motivated and often elusive to direct biological interpretation. Here, we propose a novel approach, called Thermodynamically Motivated Enrichment Analysis (TMEA), to account for the energy investment in biological relevant processes. Therefore, TMEA is based on surprisal analysis, which offers a thermodynamic-free energy-based representation of the biological steady state and of the biological change. The contribution of each biomolecule underlying the changes in free energy is used in a Monte Carlo resampling procedure resulting in a functional characterization directly coupled to the thermodynamic characterization of biological responses to system perturbations. To illustrate the utility of our method on real experimental data, we benchmark our approach on plant acclimation to high light and compare the performance of TMEA with the most frequently used method for GSEA.

Exploring Alzheimer's Disease Molecular Variability via Calculation of Personalized Transcriptional Signatures

Despite huge investments and major efforts to develop remedies for Alzheimer’s disease (AD) in the past decades, AD remains incurable. While evidence for molecular and phenotypic variability in AD have been accumulating, AD research still heavily relies on the search for AD-specific genetic/protein biomarkers that are expected to exhibit repetitive patterns throughout all patients. Thus, the classification of AD patients to different categories is expected to set the basis for the development of therapies that will be beneficial for subpopulations of patients. Here we explore the molecular heterogeneity among a large cohort of AD and non-demented brain samples, aiming to address the question whether AD-specific molecular biomarkers can progress our understanding of the disease and advance the development of anti-AD therapeutics. We studied 951 brain samples, obtained from up to 17 brain regions of 85 AD patients and 22 non-demented subjects. Utilizing an information-theoretic approach, we deciphered the brain sample-specific structures of altered transcriptional networks. Our in-depth analysis revealed that 7 subnetworks were repetitive in the 737 diseased and 214 non-demented brain samples. Each sample was characterized by a subset consisting of ~1–3 subnetworks out of 7, generating 52 distinct altered transcriptional signatures that characterized the 951 samples. We show that 30 different altered transcriptional signatures characterized solely AD samples and were not found in any of the non-demented samples. In contrast, the rest of the signatures characterized different subsets of sample types, demonstrating the high molecular variability and complexity of gene expression in AD. Importantly, different AD patients exhibiting similar expression levels of AD biomarkers harbored distinct altered transcriptional networks. Our results emphasize the need to expand the biomarker-based stratification to patient-specific transcriptional signature identification for improved AD diagnosis and for the development of subclass-specific future treatment.

Metabolic, Physiological, and Transcriptomics Analysis of Batch Cultures of the Green Microalga Chlamydomonas Grown on Different Acetate Concentrations

Acetate can be efficiently metabolized by the green microalga Chlamydomonas reinhardtii. The regular concentration is 17 mM, although higher concentrations are reported to increase starch and fatty acid content. To understand the responses to higher acetate concentrations, Chlamydomonas cells were cultivated in batch mode in the light at 17, 31, 44, and 57 mM acetate. Metabolic analyses show that cells grown at 57 mM acetate possess increased contents of all components analyzed (starch, chlorophylls, fatty acids, and proteins), with a three-fold increased volumetric biomass yield compared to cells cultivated at 17 mM acetate at the entry of stationary phase. Physiological analyses highlight the importance of photosynthesis for the low-acetate and exponential-phase samples. The stationary phase is reached when acetate is depleted, except for the cells grown at 57 mM acetate, which still divide until ammonium exhaustion. Surprisal analysis of the transcriptomics data supports the biological significance of our experiments. This allows the establishment of a model for acetate assimilation, its transcriptional regulation and the identification of candidates for genetic engineering of this metabolic pathway. Altogether, our analyses suggest that growing at high-acetate concentrations could increase biomass productivities in low-light and CO₂-limiting air-bubbled medium for biotechnology.

Intracellular redox potential is correlated with miRNA expression in MCF7 cells under hypoxic conditions

Hypoxia is a ubiquitous feature of cancers, encouraging glycolytic metabolism, proliferation, and resistance to therapy. Nonetheless, hypoxia is a poorly defined term with confounding features described in the literature. Redox biology provides an important link between the external cellular microenvironment and the cell's response to changing oxygen pressures. In this paper, we demonstrate a correlation between intracellular redox potential (measured using optical nanosensors) and the concentrations of microRNAs (miRNAs) involved in the cell's response to changes in oxygen pressure. The correlations were established using surprisal analysis (an approach derived from thermodynamics and information theory). We found that measured redox potential changes reflect changes in the free energy computed by surprisal analysis of miRNAs. Furthermore, surprisal analysis identified groups of miRNAs, functionally related to changes in proliferation and metastatic potential that played the most significant role in the cell's response to changing oxygen pressure.

Personalized disease signatures through information-theoretic compaction of big cancer data

Every individual cancer develops and grows in its own specific way, giving rise to a recognized need for the development of personalized cancer diagnostics. This suggested that the identification of patient-specific oncogene markers would be an effective diagnostics approach. However, tumors that are classified as similar according to the expression levels of certain oncogenes can eventually demonstrate divergent responses to treatment. This implies that the information gained from the identification of tumor-specific biomarkers is still not sufficient. We present a method to quantitatively transform heterogeneous big cancer data to patient-specific transcription networks. These networks characterize the unbalanced molecular processes that deviate the tissue from the normal state. We study a number of datasets spanning five different cancer types, aiming to capture the extensive interpatient heterogeneity that exists within a specific cancer type as well as between cancers of different origins. We show that a relatively small number of altered molecular processes suffices to accurately characterize over 500 tumors, showing extreme compaction of the data. Every patient is characterized by a small specific subset of unbalanced processes. We validate the result by verifying that the processes identified characterize other cancer patients as well. We show that different patients may display similar oncogene expression levels, albeit carrying biologically distinct tumors that harbor different sets of unbalanced molecular processes. Thus, tumors may be inaccurately classified and addressed as similar. These findings highlight the need to expand the notion of tumor-specific oncogenic biomarkers to patient-specific, comprehensive transcriptional networks for improved patient-tailored diagnostics.

Surprisal analysis of genome-wide transcript profiling identifies differentially expressed genes and pathways associated with four growth conditions in the microalga Chlamydomonas

The usual cultivation mode of the green microalga Chlamydomonas is liquid medium and light. However, the microalga can also be grown on agar plates and in darkness. Our aim is to analyze and compare gene expression of cells cultivated in these different conditions. For that purpose, RNA-seq data are obtained from Chlamydomonas samples of two different labs grown in four environmental conditions (agar@light, agar@dark, liquid@light, liquid@dark). The RNA seq data are analyzed by surprisal analysis, which allows the simultaneous meta-analysis of all the samples. First we identify a balance state, which defines a state where the expression levels are similar in all the samples irrespectively of their growth conditions, or lab origin. In addition our analysis identifies additional constraints needed to quantify the deviation with respect to the balance state. The first constraint differentiates the agar samples versus the liquid ones; the second constraint the dark samples versus the light ones. The two constraints are almost of equal importance. Pathways involved in stress responses are found in the agar phenotype while the liquid phenotype comprises ATP and NADH production pathways. Remodeling of membrane is suggested in the dark phenotype while photosynthetic pathways characterize the light phenotype. The same trends are also present when performing purely statistical analysis such as K-means clustering and differentially expressed genes.

Intercellular signaling through secreted proteins induces free-energy gradient-directed cell movement

Controlling cell migration is important in tissue engineering and medicine. Cell motility depends on factors such as nutrient concentration gradients and soluble factor signaling. In particular, cell-cell signaling can depend on cell-cell separation distance and can influence cellular arrangements in bulk cultures. Here, we seek a physical-based approach, which identifies a potential governed by cell-cell signaling that induces a directed cell-cell motion. A single-cell barcode chip (SCBC) was used to experimentally interrogate secreted proteins in hundreds of isolated glioblastoma brain cancer cell pairs and to monitor their relative motions over time. We used these trajectories to identify a range of cell-cell separation distances where the signaling was most stable. We then used a thermodynamics-motivated analysis of secreted protein levels to characterize free-energy changes for different cell-cell distances. We show that glioblastoma cell-cell movement can be described as Brownian motion biased by cell-cell potential. To demonstrate that the free-energy potential as determined by the signaling is the driver of motion, we inhibited two proteins most involved in maintaining the free-energy gradient. Following inhibition, cell pairs showed an essentially random Brownian motion, similar to the case for untreated, isolated single cells.

A Thermodynamic-Based Interpretation of Protein Expression Heterogeneity in Different Glioblastoma Multiforme Tumors Identifies Tumor-Specific Unbalanced Processes

We describe a thermodynamic-motivated, information theoretic analysis of proteomic data collected from a series of 8 glioblastoma multiforme (GBM) tumors. GBMs are considered here as prototypes of heterogeneous cancers. That heterogeneity is viewed here as manifesting in different unbalanced biological processes that are associated with thermodynamic-like constraints. The analysis yields a molecular description of a stable steady state that is common across all tumors. It also resolves molecular descriptions of unbalanced processes that are shared by several tumors, such as hyperactivated phosphoprotein signaling networks. Further, it resolves unbalanced processes that provide unique classifiers of tumor subgroups. The results of the theoretical interpretation are compared against those of statistical multivariate methods and are shown to provide a superior level of resolution for identifying unbalanced processes in GBM tumors. The identification of specific constraints for each GBM tumor suggests tumor-specific combination therapies that may reverse this imbalance.

Critical Points in Tumorigenesis: A Carcinogen-Initiated Phase Transition Analyzed via Single-Cell Proteomics

A kinetic, single-cell proteomic study of chemically induced carcinogenesis is interpreted by treating the single-cell data as fluctuations of an open system transitioning between different steady states. In analogy to a first-order transition, phase coexistence and the loss of degrees of freedom are observed. The transition is detected well before the appearance of the traditional biomarker of the carcinogenic transformation.

Statistical thermodynamics of transcription profiles in normal development and tumorigeneses in cohorts of patients

Experimental biology is providing the distribution of numerous different biological molecules inside cells and in body fluids of patients. Statistical methods of analysis have very successfully examined these rather large databases. We seek to use a thermodynamic analysis to provide a physical understanding and quantitative characterization of human cancers and other pathologies within a molecule-centered approach. The key technical development is the introduction of a Lagrangian. By imposing constraints the minimal value of the Lagrangian defines a thermodynamically stable state of the cellular system. The minimization also allows using experimental data measured at a number of different conditions to evaluate the steady-state distribution of biomolecules such as messenger RNAs. Thereby the number of effectively accessible quantum states of biomolecules is determined from the experimentally measured expression levels. With the increased resolution provided by the minimization of the Lagrangian one can differentiate between normal and diseased patients and further between disease subtypes. Each such refinement corresponds to imposing an additional constraint of biological origin. The constraints are the unbalanced ongoing biological processes in the system. MicroRNA expression level data for control and diseased lung cancer patients are analyzed as an example.

Computational surprisal analysis speeds-up genomic characterization of cancer processes

Surprisal analysis is increasingly being applied for the examination of transcription levels in cellular processes, towards revealing inner network structures and predicting response. But to achieve its full potential, surprisal analysis should be integrated into a wider range computational tool. The purposes of this paper are to combine surprisal analysis with other important computation procedures, such as easy manipulation of the analysis results – e.g. to choose desirable result sub-sets for further inspection –, retrieval and comparison with relevant datasets from public databases, and flexible graphical displays for heuristic thinking. The whole set of computation procedures integrated into a single practical tool is what we call Computational Surprisal Analysis. This combined kind of analysis should facilitate significantly quantitative understanding of different cellular processes for researchers, including applications in proteomics and metabolomics. Beyond that, our vision is that Computational Surprisal Analysis has the potential to reach the status of a routine method of analysis for practitioners. The resolving power of Computational Surprisal Analysis is here demonstrated by its application to a variety of cellular cancer process transcription datasets, ours and from the literature. The results provide a compact biological picture of the thermodynamic significance of the leading gene expression phenotypes in every stage of the disease. For each transcript we characterize both its inherent steady state weight, its correlation with the other transcripts and its variation due to the disease. We present a dedicated website to facilitate the analysis for researchers and practitioners.

Surprisal analysis of Glioblastoma Multiform (GBM) microRNA dynamics unveils tumor specific phenotype

Gliomablastoma multiform (GBM) is the most fatal form of all brain cancers in humans. Currently there are limited diagnostic tools for GBM detection. Here, we applied surprisal analysis, a theory grounded in thermodynamics, to unveil how biomolecule energetics, specifically a redistribution of free energy amongst microRNAs (miRNAs), results in a system deviating from a non-cancer state to the GBM cancer -specific phenotypic state. Utilizing global miRNA microarray expression data of normal and GBM patients tumors, surprisal analysis characterizes a miRNA system response capable of distinguishing GBM samples from normal tissue biopsy samples. We indicate that the miRNAs contributing to this system behavior is a disease phenotypic state specific to GBM and is therefore a unique GBM-specific thermodynamic signature. MiRNAs implicated in the regulation of stochastic signaling processes crucial in the hallmarks of human cancer, dominate this GBM-cancer phenotypic state. With this theory, we were able to distinguish with high fidelity GBM patients solely by monitoring the dynamics of miRNAs present in patients' biopsy samples. We anticipate that the GBM-specific thermodynamic signature will provide a critical translational tool in better characterizing cancer types and in the development of future therapeutics for GBM.

Surprisal analysis characterizes the free energy time course of cancer cells undergoing epithelial-to-mesenchymal transition

The epithelial-to-mesenchymal transition (EMT) initiates the invasive and metastatic behavior of many epithelial cancers. Mechanisms underlying EMT are not fully known. Surprisal analysis of mRNA time course data from lung and pancreatic cancer cells stimulated to undergo TGF-β1-induced EMT identifies two phenotypes. Examination of the time course for these phenotypes reveals that EMT reprogramming is a multistep process characterized by initiation, maturation, and stabilization stages that correlate with changes in cell metabolism. Surprisal analysis characterizes the free energy time course of the expression levels throughout the transition in terms of two state variables. The landscape of the free energy changes during the EMT for the lung cancer cells shows a stable intermediate state. Existing data suggest this is the previously proposed maturation stage. Using a single-cell ATP assay, we demonstrate that the TGF-β1-induced EMT for lung cancer cells, particularly during the maturation stage, coincides with a metabolic shift resulting in increased cytosolic ATP levels. Surprisal analysis also characterizes the absolute expression levels of the mRNAs and thereby examines the homeostasis of the transcription system during EMT.

Glioblastoma cellular architectures are predicted through the characterization of two-cell interactions

To understand how pairwise cellular interactions influence cellular architectures, we measured the levels of functional proteins associated with EGF receptor (EGFR) signaling in pairs of U87EGFR variant III oncogene receptor cells (U87EGFRvIII) at varying cell separations. Using a thermodynamics-derived approach we analyzed the cell-separation dependence of the signaling stability, and identified that the stable steady state of EGFR signaling exists when two U87EGFRvIII cells are separated by 80-100 μm. This distance range was verified as the characteristic intercellular separation within bulk cell cultures. EGFR protein network signaling coordination for the U87EGFRvIII system was lowest at the stable state and most similar to isolated cell signaling. Measurements of cultures of less tumorigenic U87PTEN cells were then used to correctly predict that stable EGFR signaling occurs for those cells at smaller cell-cell separations. The intimate relationship between functional protein levels and cellular architectures explains the scattered nature of U87EGFRvIII cells relative to U87PTEN cells in glioblastoma multiforme tumors.

miRNA and mRNA cancer signatures determined by analysis of expression levels in large cohorts of patients

Toward identifying a cancer-specific gene signature we applied surprisal analysis to the RNAs expression behavior for a large cohort of breast, lung, ovarian, and prostate carcinoma patients. We characterize the cancer phenotypic state as a shared response of a set of mRNA or microRNAs (miRNAs) in cancer patients versus noncancer controls. The resulting signature is robust with respect to individual patient variability and distinguishes with high fidelity between cancer and noncancer patients. The mRNAs and miRNAs that are implicated in the signature are correlated and are known to contribute to the regulation of cancer-signaling pathways. The miRNA and mRNA networks are common to the noncancer and cancer patients, but the disease modulates the strength of the connectivities. Furthermore, we experimentally assessed the cancer-specific signatures as possible therapeutic targets. Specifically we restructured a single dominant connectivity in the cancer-specific gene network in vitro. We find a deflection from the cancer phenotype, significantly reducing cancer cell proliferation and altering cancer cellular physiology. Our approach is grounded in thermodynamics augmented by information theory. The thermodynamic reasoning is demonstrated to ensure that the derived signature is bias-free and shows that the most significant redistribution of free energy occurs in programming a system between the noncancer and cancer states. This paper introduces a platform that can elucidate miRNA and mRNA behavior on a systems level and provides a comprehensive systematic view of both the energetics of the expression levels of RNAs and of their changes during tumorigenicity.

Surprisal analysis of transcripts expression levels in the presence of noise: a reliable determination of the onset of a tumor phenotype

Towards a reliable identification of the onset in time of a cancer phenotype, changes in transcription levels in cell models were tested. Surprisal analysis, an information-theoretic approach grounded in thermodynamics, was used to characterize the expression level of mRNAs as time changed. Surprisal Analysis provides a very compact representation for the measured expression levels of many thousands of mRNAs in terms of very few - three, four - transcription patterns. The patterns, that are a collection of transcripts that respond together, can be assigned definite biological phenotypic role. We identify a transcription pattern that is a clear marker of eventual malignancy. The weight of each transcription pattern is determined by surprisal analysis. The weight of this pattern changes with time; it is never strictly zero but it is very low at early times and then rises rather suddenly. We suggest that the low weights at early time points are primarily due to experimental noise. We develop the necessary formalism to determine at what point in time the value of that pattern becomes reliable. Beyond the point in time when a pattern is deemed reliable the data shows that the pattern remain reliable. We suggest that this allows a determination of the presence of a cancer forewarning. We apply the same formalism to the weight of the transcription patterns that account for healthy cell pathways, such as apoptosis, that need to be switched off in cancer cells. We show that their weight eventually falls below the threshold. Lastly we discuss patient heterogeneity as an additional source of fluctuation and show how to incorporate it within the developed formalism.

Free energy rhythms in Saccharomyces cerevisiae: a dynamic perspective with implications for ribosomal biogenesis

To describe the time course of cellular systems, we integrate ideas from thermodynamics and information theory to discuss the work needed to change the state of the cell. The biological example analyzed is experimental microarray transcription level oscillations of yeast in the different phases as characterized by oxygen consumption. Surprisal analysis was applied to identify groups of transcripts that oscillate in concert and thereby to compute changes in free energy with time. Three dominant transcript groups were identified by surprisal analysis. The groups correspond to the respiratory, early, and late reductive phases. Genes involved in ribosome biogenesis peaked at the respiratory phase. The work to prepare the state is shown to be the sum of the contributions of these groups. We paid particular attention to work requirements during ribosomal building, and the correlation with ATP levels and dissolved oxygen. The suggestion that cells in the respiratory phase likely build ribosomes, an energy intensive process, in preparation for protein production during the S phase of the cell cycle is validated by an experiment. Surprisal analysis thereby provided a useful tool for determining the synchronization of transcription events and energetics in a cell in real time.

Perspectives in metabolic engineering: understanding cellular regulation towards the control of metabolic routes

Metabolic engineering seeks to redirect metabolic pathways through the modification of specific biochemical reactions or the introduction of new ones with the use of recombinant technology. Many of the chemicals synthesized via introduction of product-specific enzymes or the reconstruction of entire metabolic pathways into engineered hosts that can sustain production and can synthesize high yields of the desired product as yields of natural product-derived compounds are frequently low, and chemical processes can be both energy and material expensive; current endeavors have focused on using biologically derived processes as alternatives to chemical synthesis. Such economically favorable manufacturing processes pursue goals related to sustainable development and "green chemistry". Metabolic engineering is a multidisciplinary approach, involving chemical engineering, molecular biology, biochemistry, and analytical chemistry. Recent advances in molecular biology, genome-scale models, theoretical understanding, and kinetic modeling has increased interest in using metabolic engineering to redirect metabolic fluxes for industrial and therapeutic purposes. The use of metabolic engineering has increased the productivity of industrially pertinent small molecules, alcohol-based biofuels, and biodiesel. Here, we highlight developments in the practical and theoretical strategies and technologies available for the metabolic engineering of simple systems and address current limitations.

On a fundamental structure of gene networks in living cells

Computers are organized into hardware and software. Using a theoretical approach to identify patterns in gene expression in a variety of species, organs, and cell types, we found that biological systems similarly are comprised of a relatively unchanging hardware-like gene pattern. Orthogonal patterns of software-like transcripts vary greatly, even among tumors of the same type from different individuals. Two distinguishable classes could be identified within the hardware-like component: those transcripts that are highly expressed and stable and an adaptable subset with lower expression that respond to external stimuli. Importantly, we demonstrate that this structure is conserved across organisms. Deletions of transcripts from the highly stable core are predicted to result in cell mortality. The approach provides a conceptual thermodynamic-like framework for the analysis of gene-expression levels and networks and their variations in diseased cells.