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

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

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.

Translational Components Contribute to Acclimation Responses to High Light, Heat, and Cold in Arabidopsis

Plant metabolism is broadly reprogrammed during acclimation to abiotic changes. Most previous studies have focused on transitions from standard to single stressful conditions. Here, we systematically analyze acclimation processes to levels of light, heat, and cold stress that subtly alter physiological parameters and assess their reversibility during de-acclimation. Metabolome and transcriptome changes were monitored at 11 different time points. Unlike transcriptome changes, most alterations in metabolite levels did not readily return to baseline values, except in the case of cold acclimation. Similar regulatory networks operate during (de-)acclimation to high light and cold, whereas heat and high-light responses exhibit similar dynamics, as determined by surprisal and conditional network analyses. In all acclimation models tested here, super-hubs in conditional transcriptome networks are enriched for components involved in translation, particularly ribosomes. Hence, we suggest that the ribosome serves as a common central hub for the control of three different (de-)acclimation responses.

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

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.