Analysing algorithms and data sources for the tissue-specific reconstruction of liver healthy and cancer cells

Genome-Scale Metabolic Models (GSMMs), mathematical representations of the cell metabolism in different organisms including humans, are resourceful tools to simulate metabolic phenotypes and understand associated diseases, such as obesity, diabetes and cancer. In the last years, different algorithms have been developed to generate tissue-specific metabolic models that simulate different phenotypes for distinct cell types. Hepatocyte cells are one of the main sites of metabolic conversions, mainly due to their diverse physiological functions. Most of the liver's tissue is formed by hepatocytes, being one of the largest and most important organs regarding its biological functions. Hepatocellular carcinoma is, also, one of the most important human cancers with high mortality rates. In this study, we will analyze four different algorithms (MBA, mCADRE, tINIT and FASTCORE) for tissue-specific model reconstruction, based on a template model and two types of data sources: transcriptomics and proteomics. These methods will be applied to the reconstruction of metabolic models for hepatocyte cells and HepG2 cancer cell line. The models will be analyzed and compared under different perspectives, emphasizing their functional analysis considering a set of metabolic liver tasks. The results show that there is no ``ideal'' algorithm. However, with the current analysis, we were able to retrieve knowledge about the metabolism of the liver.This study was supported by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of UID/BIO/04469/2013 unit and COMPETE 2020 (POCI-01-0145-FEDER-006684), BioTecNorte operation (NORTE01-0145-FEDER-000004) and Search-ON2: Revitalization of HPC infrastructure of UMinho, (NORTE-07-0162-FEDER-000086), all funded by European Regional Development Fund under the scope of Norte2020—Programa Operacional Regional do Norte.info:eu-repo/semantics/publishedVersio

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution

A predictive toxicogenomics signature to classify genotoxic versus non-genotoxic chemicals in human TK6 cells

Genotoxicity testing is a critical component of chemical assessment. The use of integrated approaches in genetic toxicology, including the incorporation of gene expression data to determine the DNA damage response pathways involved in response, is becoming more common. In companion papers previously published in Environmental and Molecular Mutagenesis, Li et al. (2015) [6] developed a dose optimization protocol that was based on evaluating expression changes in several well-characterized stress-response genes using quantitative real-time PCR in human lymphoblastoid TK6 cells in culture. This optimization approach was applied to the analysis of TK6 cells exposed to one of 14 genotoxic or 14 non-genotoxic agents, with sampling 4 h post-exposure. Microarray-based transcriptomic analyses were then used to develop a classifier for genotoxicity using the nearest shrunken centroids method. A panel of 65 genes was identified that could accurately classify toxicants as genotoxic or non-genotoxic. In Buick et al. (2015) [1], the utility of the biomarker for chemicals that require metabolic activation was evaluated. In this study, TK6 cells were exposed to increasing doses of four chemicals (two genotoxic that require metabolic activation and two non-genotoxic chemicals) in the presence of rat liver S9 to demonstrate that S9 does not impair the ability to classify genotoxicity using this genomic biomarker in TK6cells

Integrated network analysis identifies nitric oxide response networks and dihydroxyacid dehydratase as a crucial target in Escherichia coli

Nitric oxide (NO) is used by mammalian immune systems to counter microbial invasions and is produced by bacteria during denitrification. As a defense, microorganisms possess a complex network to cope with NO. Here we report a combined transcriptomic, chemical, and phenotypic approach to identify direct NO targets and construct the biochemical response network. In particular, network component analysis was used to identify transcription factors that are perturbed by NO. Such information was screened with potential NO reaction mechanisms and phenotypic data from genetic knockouts to identify active chemistry and direct NO targets in Escherichia coli. This approach identified the comprehensive E. coli NO response network and evinced that NO halts bacterial growth via inhibition of the branched-chain amino acid biosynthesis enzyme dihydroxyacid dehydratase. Because mammals do not synthesize branched-chain amino acids, inhibition of dihydroxyacid dehydratase may have served to foster the role of NO in the immune arsenal

Nitric oxide reaction with red blood cells and hemoglobin under heterogeneous conditions

Understanding the interaction of nitric oxide (NO) with red blood cells (RBCs) is vital to elucidating the metabolic fate of NO in the vasculature. Because hemoglobin (Hb) is the most abundant intraerythrocytic protein and reacts rapidly with NO, the interaction of NO with Hb has been studied extensively. We and others have shown the NO reaction with RBCs is nearly 1,000-fold slower than the reaction with cell-free Hb. Because the reaction rate of NO with cell-free Hb and RBCs is quite different, we hypothesize that different reaction products evolve under locally high NO concentrations, which can be generated by bolus NO addition or NO inhalation. Here we use electron paramagnetic resonance to show that bolus NO addition to cell-free Hb solutions results in nitrosylhemoglobin [HbFe(II)NO] formation as a minor product through a MetHb-dependent pathway. Further, the RBC is shown to be more prone to form HbFe(II)NO under this heterogeneous condition compared with an equivalent free-Hb solution. In both cases, trapping MetHb with cyanide blocked the formation of HbFe(II)NO. We conclude that the formation of HbFe(II)NO is a heterogeneous phenomenon involving three successive reactions of NO with the same heme molecule. These results were supported further by mathematically modeling NO–Hb reactions and diffusion