Contribution of Network Connectivity in Determining the Relationship between Gene Expression and Metabolite Concentration Changes

Abstract

One of the primary mechanisms through which a cell exerts control over its metabolic state is by modulating expression levels of its enzyme-coding genes. However, the changes at the level of enzyme expression allow only indirect control over metabolite levels, for two main reasons. First, at the level of individual reactions, metabolite levels are non-linearly dependent on enzyme abundances as per the reaction kinetics mechanisms. Secondly, specific metabolite pools are tightly interlinked with the rest of the metabolic network through their production and consumption reactions. While the role of reaction kinetics in metabolite concentration control is well studied at the level of individual reactions, the contribution of network connectivity has remained relatively unclear. Here we report a modeling framework that integrates both reaction kinetics and network connectivity constraints for describing the interplay between metabolite concentrations and mRNA levels. We used this framework to investigate correlations between the gene expression and the metabolite concentration changes in Saccharomyces cerevisiae during its metabolic cycle, as well as in response to three fundamentally different biological perturbations, namely gene knockout, nutrient shock and nutrient change. While the kinetic constraints applied at the level of individual reactions were found to be poor descriptors of the mRNA-metabolite relationship, their use in the context of the network enabled us to correlate changes in the expression of enzyme-coding genes to the alterations in metabolite levels. Our results highlight the key contribution of metabolic network connectivity in mediating cellular control over metabolite levels, and have implications towards bridging the gap between genotype and metabolic phenotype.

The Human Blood Metabolome-Transcriptome Interface

“…We were able to show that strong associations between enzymes (represented by their respective transcripts) and metabolites are significantly accumulated at shorter pathway distances (Fig 3B and 3C), which is consistent with previous studies [60,78,79]. Further functional characterization identified transport, energy, lipid and amino acid subsystems to be predominately present at short pathway distances (Fig 3D and S5 Fig).…”

Gluconeogenesis in Leishmania mexicana

“…Furthermore, PPDK has a main role in the entry of alanine and lactate into this pathway in amastigotes; however, its participation in gluconeogenesis, if needed, can be circumvented in promastigotes. It has been shown in other organisms (40,41) that during gene expression perturbations, like the genetic ablation of PEPCK and PPDK reported here, compensatory routes could take place to prevent the accumulation of metabolites inside the cells. It could explain, in part, the minimal participation of PPDK and PEPCK in gluconeogenesis observed in promastigotes and amastigotes, respectively.…”

The role of tumor metabolism as a driver of prostate cancer progression and lethal disease: results from a nested case-control study

“…In this study, we were unable to measure metabolite concentrations and inferred the metabolic state based on expression of metabolic pathways. However, the regulatory mechanisms determining the relationship between gene expression and metabolite levels are complex and not yet fully understood [38]. Transcript levels and metabolite abundances do not overlap directly with the underlying biochemical pathways [39].…”