Histone acetylation plays a central role in gene regulation and is sensitive to the levels of metabolic intermediates. However, predicting the impact of metabolic alterations on acetylation in pathological conditions is a significant challenge. Here, we present a genome-scale network model that predicts the impact of nutritional environment and genetic alterations on histone acetylation. It identifies cell types that are sensitive to histone deacetylase inhibitors based on their metabolic state, and we validate metabolites that alter drug sensitivity. Our model provides a mechanistic framework for predicting how metabolic perturbations contribute to epigenetic changes and sensitivity to deacetylase inhibitors. Electronic supplementary material The online version of this article (10.1186/s13059-019-1661-z) contains supplementary material, which is available to authorized users.

No abstract

The ultimate goal of metabolic engineering is to produce desired compounds on an industrial scale in a cost effective manner. To address challenges in metabolic engineering, computational strain optimization algorithms based on genome-scale metabolic models have increasingly been used to aid in overproducing products of interest. However, most of these strain optimization algorithms utilize a metabolic network alone, with few approaches providing strategies that also include transcriptional regulation. Moreover previous integrated approaches generally require a pre-existing regulatory network. In this study, we developed a novel strain design algorithm, named OptRAM (Optimization of Regulatory And Metabolic Networks), which can identify combinatorial optimization strategies including overexpression, knockdown or knockout of both metabolic genes and transcription factors. OptRAM is based on our previous IDREAM integrated network framework, which makes it able to deduce a regulatory network from data. OptRAM uses simulated annealing with a novel objective function, which can ensure a favorable coupling between desired chemical and cell growth. The other advance we propose is a systematic evaluation metric of multiple solutions, by considering the essential genes, flux variation, and engineering manipulation cost. We applied OptRAM to generate strain designs for succinate, 2,3-butanediol, and ethanol overproduction in yeast, which predicted high minimum predicted target production rate compared with other methods and previous literature values. Moreover, most of the genes and TFs proposed to be altered by OptRAM in these scenarios have been validated by modification of the exact genes or the target genes regulated by the TFs, for overproduction of these desired compounds by in vivo experiments cataloged in the LASER database. Particularly, we successfully validated the predicted strain optimization strategy for ethanol production by fermentation experiment. In conclusion, OptRAM can provide a useful approach that leverages an integrated transcriptional regulatory network and metabolic network to guide metabolic engineering applications.

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

Multi-scale investigation from gene transcript level to metabolic activity is important to uncover plant response to environment perturbation. Here we integrated a genome-scale constraint-based metabolic model with transcriptome data to explore Arabidopsis thaliana response to both elevated and low CO2 conditions. The four condition-specific models from low to high CO2 concentrations show differences in active reaction sets, enriched pathways for increased/decreased fluxes, and putative post-transcriptional regulation, which indicates that condition-specific models are necessary to reflect physiological metabolic states. The simulated CO2 fixation flux at different CO2 concentrations is consistent with the measured Assimilation-CO2intercellular curve. Interestingly, we found that reactions in primary metabolism are affected most significantly by CO2 perturbation, whereas secondary metabolic reactions are not influenced a lot. The changes predicted in key pathways are consistent with existing knowledge. Another interesting point is that Arabidopsis is required to make stronger adjustment on metabolism to adapt to the more severe low CO2 stress than elevated CO2 . The challenges of identifying post-transcriptional regulation could also be addressed by the integrative model. In conclusion, this innovative application of multi-scale modeling in plants demonstrates potential to uncover the mechanisms of metabolic response to different conditions.

BackgroundRice (Oryza sativa) is one of the most important grain crops, which serves as food source for nearly half of the world population. The study of rice development process as well as related strategies for production has made significant progress. However, the comprehensive study on development of different rice tissues at both transcriptomic and metabolic flux level across different stages was lacked.ResultsIn this study, we performed RNA-Seq and characterized the expression profiles of differentiated tissues from Oryza sativa Zhonghua 11, including leaves, sheath, stamen, pistil, lemma and palea of the booting stage, and embryo, endosperm, lemma and palea of the mature grain stage. By integrating this set of transcriptome data of different rice tissues at different stages with a genome-scale rice metabolic model, we generated tissue-specific models and investigated the shift of metabolic patterns, and the discrepancy between transcriptomic and metabolic level. We found although the flux patterns are not very similar with the gene expression pattern, the tissues at booting stage and mature grain stage can be separately clustered by primary metabolism at either level. While the gene expression and flux distribution of secondary metabolism is more diverse across tissues and stages. The critical rate-limiting reactions and pathways were also identified. In addition, we compared the patterns of the same tissue at different stages and the different tissues at same stage. There are more altered pathways at gene expression level than metabolic level, which indicate the metabolism is more robust to reflect the phenotype, and might largely because of the complex post-transcriptional modification.ConclusionsThe tissue-specific models revealed more detail metabolic pattern shift among different tissues and stages, which is of great significance to uncover mechanism of rice grain development and further improve production and quality of rice.Electronic supplementary materialThe online version of this article (10.1186/s12918-018-0574-x) contains supplementary material, which is available to authorized users.

Cell cycle is a fundamental process for cell growth and proliferation, and its dysregulation leads to many diseases. How metabolic networks are regulated and rewired during the cell cycle is unknown. Here we apply a dynamic genome-scale metabolic modeling framework (DFA) to simulate a cell cycle of cytokine-activated murine pro-B cells. Phase-specific reaction activity predicted by DFA using time-course metabolomics were validated using matched time-course proteomics and phospho-proteomics data. Our model correctly predicted changes in methionine metabolism at the G1/S transition and the activation of lysine metabolism, nucleotides synthesis, fatty acid elongation and heme biosynthesis at the critical G0/G1 transition into cell growth and proliferation. Metabolic fluxes predicted from proteomics and phosphoproteomics constrained metabolic models were highly consistent with DFA fluxes and revealed that most reaction fluxes are regulated indirectly. Our model can help predict the impact of changes in nutrients, enzymes, or regulators on this critical cellular process.

8The metabolism of most organisms is controlled by a diverse cast of regulatory processes, including 9transcriptional regulation and post-translational modifications (PTMs). Yet how metabolic control is 10 distributed between these regulatory processes is unknown. Here we present Comparative Analysis of 11Regulators of Metabolism (CAROM), an approach that compares regulators based on network 12 connectivity, flux, and essentiality of their reaction targets. Using CAROM, we analyze transcriptome, 13proteome, acetylome and phospho-proteome dynamics during transition to stationary phase in E. coli 14 and S. cerevisiae. CAROM uncovered that the targets of each regulatory process shared unique 15 metabolic properties: growth-limiting reactions were regulated by acetylation, while isozymes and futile-16 cycles were preferentially regulated by phosphorylation. Reversibility, essentiality, and molecular-17weight further distinguished reactions controlled through diverse mechanisms. While every enzyme can 18 be potentially regulated by multiple mechanisms, analysis of context-specific datasets reveals a 19 conserved partitioning of metabolic regulation based on reaction attributes. 20Author summary 21There are several ways to regulate an enzyme's activity in a cell. Yet, the design principles that 22determine when an enzyme is regulated by transcription, translation or post-translational modifications 23 are unknown. Each control mechanism, such as transcription, comprises several regulators that control 24 a distinct set of targets. So far, it is unclear if similar partitioning of targets occurs at a higher level, 25between different control mechanisms. Here we systematically analyze patterns of metabolic regulation 26 in model microbes. We find that five key parameters can distinguish the targets of each mechanism. 27These key parameters provide insights on specific roles played by each mechanism in determining 28 overall metabolic activity. This approach may help define the basic regulatory architecture of metabolic 29 networks. 30