BayFlux: ABayesian method to quantify metabolicFluxes and their uncertainty at the genome scale

Author:

Backman Tyler W. H.ORCID,Schenk Christina,Radivojevic Tijana,Ando David,Singh Janavi,Czajka Jeffrey J.,Costello Zak,Keasling Jay D.,Tang Yinjie,Akhmatskaya Elena,Martin Hector GarciaORCID

Abstract

AbstractMetabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function.13C Metabolic Flux Analysis (13C MFA) is considered to be the gold standard for measuring metabolic fluxes.13C MFA typically works by leveraging extracellular exchange fluxes as well as data from13C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the13C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in13C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from13C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P-13C MOMA and ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods. We provide an open source Python implementation of BayFlux athttps://github.com/JBEI/bayflux.Author summary13C MFA practitioners know that modeling results can be sensitive to minor modifications of the metabolic model. Certain parts of the metabolic model that are not well mapped to a molecular mechanism (e.g.drains to biomass or ATP maintenance) can have an inordinate impact on the final fluxes. The only way to ascertain the validity of the model is by checking that the result does not significantly differ from previously observed flux profiles. However, that approach diminishes the possibility of discovering truly novel flux profiles. Because of this strong dependence on metabolic model details, it would be very useful to have a systematic and repeatable way to produce these metabolic models. And indeed there is one: genome-scale metabolic models can be systematically obtained from genomic sequences, and represent all the known genomically encoded metabolic information. However, these models are much larger than the traditionally used central carbon metabolism models. Hence, the number of degrees of freedom of the model (fluxes) significantly exceeds the number of measurements (metabolite labeling profiles and exchange fluxes). As a result, one expects many flux profiles compatible with the experimental data. The best way to represent these is by identifying all fluxes compatible with the experimental data. Our novel method BayFlux, based on Bayesian inference and Markov Chain Monte Carlo sampling, provides this capability. Interestingly, this approach leads to the observation that traditional optimization approaches can significantly overestimate flux uncertainty, and that genome-scale models of metabolism produce narrower flux distributions than the small core metabolic models that are traditionally used in13C MFA. Furthermore, we show that the extra information provided by this approach allows us to improve knockout predictions, compared to traditional methods. Although the method scales well with more reactions, improvements will be needed to tackle the large metabolic models found in microbiomes and human metabolism.

Publisher

Cold Spring Harbor Laboratory

Reference54 articles.

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3