Comparison of metabolic states using genome-scale metabolic models

Abstract

Genome-scale metabolic models (GEMs) are comprehensive knowledge bases of cellular metabolism and serve as mathematical tools for studying biological phenotypes and metabolic states or conditions in various organisms and cell types. Given the sheer size and complexity of human metabolism, selecting parameters for existing analysis methods such as metabolic objective functions and model constraints is not straightforward in human GEMs. In particular, comparing several conditions in large GEMs to identify condition- or disease-specific metabolic features is challenging. In this study, we showcase a scalable, model-driven approach for an in-depth investigation and comparison of metabolic states in large GEMs which enables identifying the underlying functional differences. Using a combination of flux space sampling and network analysis, our approach enables extraction and visualisation of metabolically distinct network modules. Importantly, it does not rely on known or assumed objective functions. We apply this novel approach to extract the biochemical differences in adipocytes arising due to unlimited vs blocked uptake of branched-chain amino acids (BCAAs, considered as biomarkers in obesity) using a human adipocyte GEM (iAdipocytes1809). The biological significance of our approach is corroborated by literature reports confirming our identified metabolic processes (TCA cycle and Fatty acid metabolism) to be functionally related to BCAA metabolism. Additionally, our analysis predicts a specific altered uptake and secretion profile indicating a compensation for the unavailability of BCAAs. Taken together, our approach facilitates determining functional differences between any metabolic conditions of interest by offering a versatile platform for analysing and comparing flux spaces of large metabolic networks. Author summary Cellular metabolism is a highly complex and interconnected system. As many lifestyle diseases in humans have a strong metabolic component, it is important to understand metabolic differences between healthy and diseased states. In systems biology, metabolic behaviours are investigated using genome-scale metabolic models. In addition to the sheer size and complexity of the genome-scale metabolic models of human systems, using existing analysis methods is challenging and the parameter selection is not straightforward. Therefore, novel methodological frameworks are necessary for analysing metabolic conditions despite the challenges posed by human models. Particularly, an ongoing challenge has been that of comparing several phenotypes for identifying condition- or disease-specific metabolic signatures. We address this significant challenge by developing a scalable and model-driven approach, ComMet (Comparison of Metabolic states). ComMet enables an in-depth investigation and comparison of metabolic phenotypes in large models while also identifying the underlying functional differences. Novel hypotheses can be generated using ComMet for not only understanding known metabolic phenotypes better but also for guiding the design of new experiments to validate the processes predicted by ComMet.