Dynamic Network Modeling of Stem Cell Metabolism

Protocol for CAROM: A machine learning tool to predict post-translational regulation from metabolic signatures

This protocol describes CAROM, a computational tool that combines genome-scale metabolic networks (GEMs) and machine learning to identify enzyme targets of post-translational modifications (PTMs). Condition-specific enzyme and reaction properties are used to predict targets of phosphorylation and acetylation in multiple organisms. CAROM is influenced by the accuracy of GEMs and associated flux-balance analysis (FBA), which generate the inputs of the model. We demonstrate the protocol using multi-omics data from E. coli. For complete details on the use and execution of this protocol, please refer to Smith et al. (2022).

A flux-based machine learning model to simulate the impact of pathogen metabolic heterogeneity on drug interactions

Drug combinations are a promising strategy to counter antibiotic resistance. However, current experimental and computational approaches do not account for the entire complexity involved in combination therapy design, such as the effect of pathogen metabolic heterogeneity, changes in the growth environment, drug treatment order, and time interval. To address these limitations, we present a comprehensive approach that uses genome-scale metabolic modeling and machine learning to guide combination therapy design. Our mechanistic approach (a) accommodates diverse data types, (b) accounts for time- and order-specific interactions, and (c) accurately predicts drug interactions in various growth conditions and their robustness to pathogen metabolic heterogeneity. Our approach achieved high accuracy (area under the receiver operating curve (AUROC) = 0.83 for synergy, AUROC = 0.98 for antagonism) in predicting drug interactions for Escherichia coli cultured in 57 metabolic conditions based on experimental validation. The entropy in bacterial metabolic response was predictive of combination therapy outcomes across time scales and growth conditions. Simulation of metabolic heterogeneity using population FBA identified two subpopulations of E. coli cells defined by the levels of three proteins (eno, fadB, and fabD) in glycolysis and lipid metabolism that influence cell tolerance to a broad range of antibiotic combinations. Analysis of the vast landscape of condition-specific drug interactions revealed a set of 24 robustly synergistic drug combinations with potential for clinical use.

Acetyl-CoA metabolism drives epigenome change and contributes to carcinogenesis risk in fatty liver disease

BACKGROUND

The incidence of non-alcoholic fatty liver disease (NAFLD)-associated hepatocellular carcinoma (HCC) is increasing worldwide, but the steps in precancerous hepatocytes which lead to HCC driver mutations are not well understood. Here we provide evidence that metabolically driven histone hyperacetylation in steatotic hepatocytes can increase DNA damage to initiate carcinogenesis.

METHODS

Global epigenetic state was assessed in liver samples from high-fat diet or high-fructose diet rodent models, as well as in cultured immortalized human hepatocytes (IHH cells). The mechanisms linking steatosis, histone acetylation and DNA damage were investigated by computational metabolic modelling as well as through manipulation of IHH cells with metabolic and epigenetic inhibitors. Chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) and transcriptome (RNA-seq) analyses were performed on IHH cells. Mutation locations and patterns were compared between the IHH cell model and genome sequence data from preneoplastic fatty liver samples from patients with alcohol-related liver disease and NAFLD.

RESULTS

Genome-wide histone acetylation was increased in steatotic livers of rodents fed high-fructose or high-fat diet. In vitro, steatosis relaxed chromatin and increased DNA damage marker γH2AX, which was reversed by inhibiting acetyl-CoA production. Steatosis-associated acetylation and γH2AX were enriched at gene clusters in telomere-proximal regions which contained HCC tumour suppressors in hepatocytes and human fatty livers. Regions of metabolically driven epigenetic change also had increased levels of DNA mutation in non-cancerous tissue from NAFLD and alcohol-related liver disease patients. Finally, genome-scale network modelling indicated that redox balance could be a key contributor to this mechanism.

CONCLUSIONS

Abnormal histone hyperacetylation facilitates DNA damage in steatotic hepatocytes and is a potential initiating event in hepatocellular carcinogenesis.

Lipid exposure activates gene expression changes associated with estrogen receptor negative breast cancer

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A and MCF-12A cells, and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility, and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

Metabolic signatures of regulation by phosphorylation and acetylation

Acetylation and phosphorylation are highly conserved posttranslational modifications (PTMs) that regulate cellular metabolism, yet how metabolic control is shared between these PTMs is unknown. Here we analyze transcriptome, proteome, acetylome, and phosphoproteome datasets in E. coli, S. cerevisiae, and mammalian cells across diverse conditions using CAROM, a new approach that uses genome-scale metabolic networks and machine learning to classify targets of PTMs. We built a single machine learning model that predicted targets of each PTM in a condition across all three organisms based on reaction attributes (AUC>0.8). Our model predicted phosphorylated enzymes during a mammalian cell-cycle, which we validate using phosphoproteomics. Interpreting the machine learning model using game theory uncovered enzyme properties including network connectivity, essentiality, and condition-specific factors such as maximum flux that differentiate targets of phosphorylation from acetylation. The conserved and predictable partitioning of metabolic regulation identified here between these PTMs may enable rational rewiring of regulatory circuits.

Metabolic remodelling during early mouse embryo development

During early mammalian embryogenesis, changes in cell growth and proliferation depend on strict genetic and metabolic instructions. However, our understanding of metabolic reprogramming and its influence on epigenetic regulation in early embryo development remains elusive. Here we show a comprehensive metabolomics profiling of key stages in mouse early development and the two-cell and blastocyst embryos, and we reconstructed the metabolic landscape through the transition from totipotency to pluripotency. Our integrated metabolomics and transcriptomics analysis shows that while two-cell embryos favour methionine, polyamine and glutathione metabolism and stay in a more reductive state, blastocyst embryos have higher metabolites related to the mitochondrial tricarboxylic acid cycle, and present a more oxidative state. Moreover, we identify a reciprocal relationship between α-ketoglutarate (α-KG) and the competitive inhibitor of α-KG-dependent dioxygenases, L-2-hydroxyglutarate (L-2-HG), where two-cell embryos inherited from oocytes and one-cell zygotes display higher L-2-HG, whereas blastocysts show higher α-KG. Lastly, increasing 2-HG availability impedes erasure of global histone methylation markers after fertilization. Together, our data demonstrate dynamic and interconnected metabolic, transcriptional and epigenetic network remodelling during early mouse embryo development.

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.

Dissecting Murine Muscle Stem Cell Aging through Regeneration Using Integrative Genomic Analysis

During aging, there is a progressive loss of volume and function in skeletal muscle that impacts mobility and quality of life. The repair of skeletal muscle is regulated by tissue-resident stem cells called satellite cells (or muscle stem cells [MuSCs]), but in aging, MuSCs decrease in numbers and regenerative capacity. The transcriptional networks and epigenetic changes that confer diminished regenerative function in MuSCs as a result of natural aging are only partially understood. Herein, we use an integrative genomics approach to profile MuSCs from young and aged animals before and after injury. Integration of these datasets reveals aging impacts multiple regulatory changes through significant differences in gene expression, metabolic flux, chromatin accessibility, and patterns of transcription factor (TF) binding activities. Collectively, these datasets facilitate a deeper understanding of the regulation tissue-resident stem cells use during aging and healing.

Nutrient Sensing by Histone Marks: Reading the Metabolic Histone Code Using Tracing, Omics, and Modeling

Several metabolites serve as substrates for histone modifications and communicate changes in the metabolic environment to the epigenome. Technologies such as metabolomics and proteomics have allowed us to reconstruct the interactions between metabolic pathways and histones. These technologies have shed light on how nutrient availability can have a dramatic effect on various histone modifications. This metabolism-epigenome cross talk plays a fundamental role in development, immune function, and diseases like cancer. Yet, major challenges remain in understanding the interactions between cellular metabolism and the epigenome. How the levels and fluxes of various metabolites impact epigenetic marks is still unclear. Discussed herein are recent applications and the potential of systems biology methods such as flux tracing and metabolic modeling to address these challenges and to uncover new metabolic-epigenetic interactions. These systems approaches can ultimately help elucidate how nutrients shape the epigenome of microbes and mammalian cells.

Current Status and Future Prospects of Genome-Scale Metabolic Modeling to Optimize the Use of Mesenchymal Stem Cells in Regenerative Medicine

Mesenchymal stem cells are a promising source for externally grown tissue replacements and patient-specific immunomodulatory treatments. This promise has not yet been fulfilled in part due to production scaling issues and the need to maintain the correct phenotype after re-implantation. One aspect of extracorporeal growth that may be manipulated to optimize cell growth and differentiation is metabolism. The metabolism of MSCs changes during and in response to differentiation and immunomodulatory changes. MSC metabolism may be linked to functional differences but how this occurs and influences MSC function remains unclear. Understanding how MSC metabolism relates to cell function is however important as metabolite availability and environmental circumstances in the body may affect the success of implantation. Genome-scale constraint based metabolic modeling can be used as a tool to fill gaps in knowledge of MSC metabolism, acting as a framework to integrate and understand various data types (e.g., genomic, transcriptomic and metabolomic). These approaches have long been used to optimize the growth and productivity of bacterial production systems and are being increasingly used to provide insights into human health research. Production of tissue for implantation using MSCs requires both optimized production of cell mass and the understanding of the patient and phenotype specific metabolic situation. This review considers the current knowledge of MSC metabolism and how it may be optimized along with the current and future uses of genome scale constraint based metabolic modeling to further this aim.

Tissue of origin dictates GOT1 dependence and confers synthetic lethality to radiotherapy

BACKGROUND

Metabolic programs in cancer cells are influenced by genotype and the tissue of origin. We have previously shown that central carbon metabolism is rewired in pancreatic ductal adenocarcinoma (PDA) to support proliferation through a glutamate oxaloacetate transaminase 1 (GOT1)-dependent pathway.

METHODS

We utilized a doxycycline-inducible shRNA-mediated strategy to knockdown GOT1 in PDA and colorectal cancer (CRC) cell lines and tumor models of similar genotype. These cells were analyzed for the ability to form colonies and tumors to test if tissue type impacted GOT1 dependence. Additionally, the ability of GOT1 to impact the response to chemo- and radiotherapy was assessed. Mechanistically, the associated specimens were examined using a combination of steady-state and stable isotope tracing metabolomics strategies and computational modeling. Statistics were calculated using GraphPad Prism 7. One-way ANOVA was performed for experiments comparing multiple groups with one changing variable. Student's t test (unpaired, two-tailed) was performed when comparing two groups to each other. Metabolomics data comparing three PDA and three CRC cell lines were analyzed by performing Student's t test (unpaired, two-tailed) between all PDA metabolites and CRC metabolites.

RESULTS

While PDA exhibits profound growth inhibition upon GOT1 knockdown, we found CRC to be insensitive. In PDA, but not CRC, GOT1 inhibition disrupted glycolysis, nucleotide metabolism, and redox homeostasis. These insights were leveraged in PDA, where we demonstrate that radiotherapy potently enhanced the effect of GOT1 inhibition on tumor growth.

CONCLUSIONS

Taken together, these results illustrate the role of tissue type in dictating metabolic dependencies and provide new insights for targeting metabolism to treat PDA.

The metabolic network model of primed/naive human embryonic stem cells underlines the importance of oxidation-reduction potential and tryptophan metabolism in primed pluripotency

Background

Pluripotency is proposed to exist in two different stages: Naive and Primed. Conventional human pluripotent cells are essentially in the primed stage. In recent years, several protocols have claimed to generate naive human embryonic stem cells (hESCs). To the best of our knowledge, none of these protocols is currently recognized as the gold standard method. Furthermore, the consistency of the resulting cells from these diverse protocols at the molecular level is yet to be shown. Additionally, little is known about the principles that govern the metabolic differences between naive and primed pluripotency. In this work, using a computational approach, we tried to shed light on these basic issues.

Results

We showed that, after batch effect removal, the transcriptome data of eight different protocols which supposedly produce naive hESCs are clustered consistently when compared to the primed ones. Next, by integrating transcriptomes of all hESCs obtained by these protocols, we reconstructed p-hESCNet and n-hESCNet, the first metabolic network models representing hESCs. By exploiting reporter metabolite analysis we showed that the status of NAD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{+}$$\end{document}+ and the metabolites involved in the TCA cycle are significantly altered between naive and primed hESCs. Furthermore, using flux variability analysis (FVA), the models showed that the kynurenine-mediated metabolism of tryptophan is remarkably downregulated in naive human pluripotent cells.

Conclusion

The aim of the present paper is twofold. Firstly, our findings confirm the applicability of all these protocols for generating naive hESCs, due to their consistency at the transcriptome level. Secondly, we showed that in silico metabolic models of hESCs can be used to simulate the metabolic states of naive and primed pluripotency. Our models confirmed the OXPHOS activation in naive cells and showed that oxidation-reduction potential vary between naive and primed cells. Tryptophan metabolism is also outlined as a key pathway in primed pluripotency and the models suggest that decrements in the activity of this pathway might be an appropriate marker for naive pluripotency.