RESUMEN
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.
Asunto(s)
Redes y Vías Metabólicas/genética , Modelos Biológicos , Adipocitos/metabolismo , Algoritmos , Aminoácidos de Cadena Ramificada/metabolismo , Ciclo del Ácido Cítrico , Biología Computacional , Simulación por Computador , Ácidos Grasos/metabolismo , Genoma Humano , Humanos , Enfermedades Metabólicas/genética , Enfermedades Metabólicas/metabolismo , Análisis de Flujos Metabólicos/estadística & datos numéricos , Modelos Genéticos , Obesidad/genética , Obesidad/metabolismo , Análisis de Componente PrincipalRESUMEN
Elementary Flux Modes (EFMs) are a tool for constraint-based modeling and metabolic network analysis. However, systematic and automated visualization of EFMs, capable of integrating various data types is still a challenge. In this study, we developed an extension for the widely adopted COBRA Toolbox, EFMviz, for analysis and graphical visualization of EFMs as networks of reactions, metabolites and genes. The analysis workflow offers a platform for EFM visualization to improve EFM interpretability by connecting COBRA toolbox with the network analysis and visualization software Cytoscape. The biological applicability of EFMviz is demonstrated in two use cases on medium (Escherichia coli, iAF1260) and large (human, Recon 2.2) genome-scale metabolic models. EFMviz is open-source and integrated into COBRA Toolbox. The analysis workflows used for the two use cases are detailed in the two tutorials provided with EFMviz along with the data used in this study.
RESUMEN
The link between the experimental laboratory studies and bioinformatic approaches aims to find procedures to connect tools from both branches producing workflows that bring together different techniques that are capable of exploiting data at many levels. Thanks to the open access sources and the numerous tools available, it is possible to create various pipelines capable of solving specific problems. Nevertheless the lack of connectivity between them that interconnect different approaches complicates the exploitation of these workflows. Here, we present a detailed description of a workflow composed of different bioinformatics tools that exploits data from large-scale gene expression experiments, contextualizing them at many biological levels. To illustrate the relevance of our workflow for the vitamin D community we applied it to data from myeloid cell models treated with the hormonally active form of vitamin. From raw files of functional genomic studies it is possible to utilize the whole information to obtain a biological insight. Different software and algorithms are included to analyse at pathway, metabolic, ontology and molecular biology level the effects on gene expression. The usage of different databases to analyse gene expression data allows to perform a complete interpretation of functional genomic studies and the implementation of analysis and visualization software tools gives a better understanding of the biological meaning of the results. This review is an example of how to select and bring together several software modules to create one pipeline that processes and analyses genomic data at many biological levels making it open, reproducible and user friendly. Finally, the application of our bioinformatic pipeline revealed that vitamin D modulates crucial metabolic pathways in different myeloid cells that may play an important role in their immune function.