Computational analysis of reciprocal association of metabolism and epigenetics in the budding yeast: a genome-scale metabolic model (GSMM) approach

Reconstruction of cell-specific models capturing the influence of metabolism on DNA methylation in cancer

The imbalance of epigenetic regulatory mechanisms such as DNA methylation, which can promote aberrant gene expression profiles without affecting the DNA sequence, may cause the deregulation of signaling, regulatory, and metabolic processes, contributing to a cancerous phenotype. Since some metabolites are substrates and cofactors of epigenetic regulators, their availability can be affected by characteristic cancer cell metabolic shifts, feeding cancer onset and progression through epigenetic deregulation. Hence, there is a need to study the influence of cancer metabolic reprogramming in DNA methylation to design new effective treatments. In this study, a generic Genome-Scale Metabolic Model (GSMM) of a human cell, integrating DNA methylation or demethylation reactions, was obtained and used for the reconstruction of Genome-Scale Metabolic Models enhanced with Enzymatic Constraints using Kinetic and Omics data (GECKOs) of 31 cancer cell lines. Furthermore, cell-line-specific DNA methylation levels were included in the models, as coefficients of a DNA composition pseudo-reaction, to depict the influence of metabolism over global DNA methylation in each of the cancer cell lines. Flux simulations demonstrated the ability of these models to provide simulated fluxes of exchange reactions similar to the equivalent experimentally measured uptake/secretion rates and to make good functional predictions. In addition, simulations found metabolic pathways, reactions and enzymes directly or inversely associated with the gene promoter methylation. Two potential candidates for targeted cancer epigenetic therapy were identified.

Genome-Wide Analysis of Yeast Metabolic Cycle through Metabolic Network Models Reveals Superiority of Integrated ATAC-seq Data over RNA-seq Data

Saccharomyces cerevisiae undergoes robust oscillations to regulate its physiology for adaptation and survival under nutrient-limited conditions. Environmental cues can induce rhythmic metabolic alterations in order to facilitate the coordination of dynamic metabolic behaviors. Of such metabolic processes, the yeast metabolic cycle enables adaptation of the cells to varying nutritional status through oscillations in gene expression and metabolite production levels. In this process, yeast metabolism is altered between diverse cellular states based on changing oxygen consumption levels: quiescent (reductive charging [RC]), growth (oxidative [OX]), and proliferation (reductive building [RB]) phases. We characterized metabolic alterations during the yeast metabolic cycle using a variety of approaches. Gene expression levels are widely used for condition-specific metabolic simulations, whereas the use of epigenetic information in metabolic modeling is still limited despite the clear relationship between epigenetics and metabolism. This prompted us to investigate the contribution of epigenomic information to metabolic predictions for progression of the yeast metabolic cycle. In this regard, we determined altered pathways through the prediction of regulated reactions and corresponding model genes relying on differential chromatin accessibility levels. The predicted metabolic alterations were confirmed via data analysis and literature. We subsequently utilized RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data sets in the contextualization of the yeast model. The use of ATAC-seq data considerably enhanced the predictive capability of the model. To the best of our knowledge, this is the first attempt to use genome-wide chromatin accessibility data in metabolic modeling. The preliminary results showed that epigenomic data sets can pave the way for more accurate metabolic simulations.

IMPORTANCE Dynamic chromatin organization mediates the emergence of condition-specific phenotypes in eukaryotic organisms. Saccharomyces cerevisiae can alter its metabolic profile via regulation of genome accessibility and robust transcriptional oscillations under nutrient-limited conditions. Thus, both epigenetic information and transcriptomic information are crucial in the understanding of condition-specific metabolic behavior in this organism. Based on genome-wide alterations in chromatin accessibility and transcription, we investigated the yeast metabolic cycle, which is a remarkable example of coordinated and dynamic yeast behavior. In this regard, we assessed the use of ATAC-seq and RNA-seq data sets in condition-specific metabolic modeling. To our knowledge, this is the first attempt to use chromatin accessibility data in the reconstruction of context-specific metabolic models, despite the extensive use of transcriptomic data. As a result of comparative analyses, we propose that the incorporation of epigenetic information is a promising approach in the accurate prediction of metabolic dynamics.

Next-Generation Genome-Scale Metabolic Modeling through Integration of Regulatory Mechanisms

Genome-scale metabolic models (GEMs) are powerful tools for understanding metabolism from a systems-level perspective. However, GEMs in their most basic form fail to account for cellular regulation. A diverse set of mechanisms regulate cellular metabolism, enabling organisms to respond to a wide range of conditions. This limitation of GEMs has prompted the development of new methods to integrate regulatory mechanisms, thereby enhancing the predictive capabilities and broadening the scope of GEMs. Here, we cover integrative models encompassing six types of regulatory mechanisms: transcriptional regulatory networks (TRNs), post-translational modifications (PTMs), epigenetics, protein-protein interactions and protein stability (PPIs/PS), allostery, and signaling networks. We discuss 22 integrative GEM modeling methods and how these have been used to simulate metabolic regulation during normal and pathological conditions. While these advances have been remarkable, there remains a need for comprehensive and widespread integration of regulatory constraints into GEMs. We conclude by discussing challenges in constructing GEMs with regulation and highlight areas that need to be addressed for the successful modeling of metabolic regulation. Next-generation integrative GEMs that incorporate multiple regulatory mechanisms and their crosstalk will be invaluable for discovering cell-type and disease-specific metabolic control mechanisms.

Tying Metabolic Branches With Histone Tails Using Systems Biology

Histone modifications represent an innate cellular mechanism to link nutritional status to gene expression. Metabolites such as acetyl-CoA and S-adenosyl methionine influence gene expression by serving as substrates for modification of histones. Yet, we lack a predictive model for determining histone modification levels based on cellular metabolic state. The numerous metabolic pathways that intersect with histone marks makes it highly challenging to understand their interdependencies. Here, we highlight new systems biology tools to unravel the impact of nutritional cues and metabolic fluxes on histone modifications.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

Genome-Scale Metabolic Modeling from Yeast to Human Cell Models of Complex Diseases: Latest Advances and Challenges

Genome-scale metabolic models (GEMs) are mathematical models that enable systematic analysis of metabolism. This modeling concept has been applied to study the metabolism of many organisms including the eukaryal model organism, the yeast Saccharomyces cerevisiae, that also serves as an important cell factory for production of fuels and chemicals. With the application of yeast GEMs, our knowledge of metabolism is increasing. Therefore, GEMs have also been used for modeling human cells to study metabolic diseases. Here we introduce the concept of GEMs and provide a protocol for reconstructing GEMs. Besides, we show the historic development of yeast GEMs and their applications. Also, we review human GEMs as well as their uses in the studies of complex diseases.

Perspectives on Systems Modeling of Human Peripheral Blood Mononuclear Cells

Human peripheral blood mononuclear cells (PBMCs) are the key drivers of the immune responses. These cells undergo activation, proliferation and differentiation into various subsets. During these processes they initiate metabolic reprogramming, which is coordinated by specific gene and protein activities. PBMCs as a model system have been widely used to study metabolic and autoimmune diseases. Herein we review various omics and systems-based approaches such as transcriptomics, epigenomics, proteomics, and metabolomics as applied to PBMCs, particularly T helper subsets, that unveiled disease markers and the underlying mechanisms. We also discuss and emphasize several aspects of T cell metabolic modeling in healthy and disease states using genome-scale metabolic models.

Genome scale models of yeast: towards standardized evaluation and consistent omic integration

Genome scale models (GEMs) have enabled remarkable advances in systems biology, acting as functional databases of metabolism, and as scaffolds for the contextualization of high-throughput data. In the case of Saccharomyces cerevisiae (budding yeast), several GEMs have been published and are currently used for metabolic engineering and elucidating biological interactions. Here we review the history of yeast's GEMs, focusing on recent developments. We study how these models are typically evaluated, using both descriptive and predictive metrics. Additionally, we analyze the different ways in which all levels of omics data (from gene expression to flux) have been integrated in yeast GEMs. Relevant conclusions and current challenges for both GEM evaluation and omic integration are highlighted.