One lump or two?

We investigate methods for modelling metabolism within populations of cells. Typically one represents the interaction of a cloned population of cells with their environment as though it were one large cell. The question is as to whether any dynamics are lost by this assumption, and as to whether it might be more appropriate to instead model each cell individually. We show that it is sufficient to model at an intermediate level of granularity, representing the population as two interacting lumps of tissue

Characterisation of multiple substrate-specific (d)ITP/(d)XTPase and modelling of deaminated purine nucleotide metabolism

To be viable, organisms possess a number of (deoxy)nucleotide phosphohydrolases, which hydrolyze these nucleotides removing them from the active NTP and dNTP pools. Deamination of purine bases can result in accumulation of such nucleotides as ITP, dITP, XTP and dXTP. E. coli RdgB has been characterised as a deoxyribonucleoside triphosphate pyrophosphohydrolase that can act on these nucleotides. S. cerevisiae homologue encoded by YJR069C was purified and its (d)NTPase activity was assayed using fifteen nucleotide substrates. ITP, dITP, and XTP were identified as major substrates and kinetic parameters measured. Inhibition by ATP, dATP and GTP were established. On the basis of experimental and published data, modelling and simulation of ITP, dITP, XTP and dXTP metabolism was performed. (d)ITP/(d)XTPase is a new example of enzyme with multiple substrate-specificity demonstrating that multispecificity is not a rare phenomenon

Towards a genome-scale kinetic model of cellular metabolism

<p>Abstract</p> <p>Background</p> <p>Advances in bioinformatic techniques and analyses have led to the availability of genome-scale metabolic reconstructions. The size and complexity of such networks often means that their potential behaviour can only be analysed with constraint-based methods. Whilst requiring minimal experimental data, such methods are unable to give insight into cellular substrate concentrations. Instead, the long-term goal of systems biology is to use kinetic modelling to characterize fully the mechanics of each enzymatic reaction, and to combine such knowledge to predict system behaviour.</p> <p>Results</p> <p>We describe a method for building a parameterized genome-scale kinetic model of a metabolic network. Simplified linlog kinetics are used and the parameters are extracted from a kinetic model repository. We demonstrate our methodology by applying it to yeast metabolism. The resultant model has 956 metabolic reactions involving 820 metabolites, and, whilst approximative, has considerably broader remit than any existing models of its type. Control analysis is used to identify key steps within the system.</p> <p>Conclusions</p> <p>Our modelling framework may be considered a stepping-stone toward the long-term goal of a fully-parameterized model of yeast metabolism. The model is available in SBML format from the BioModels database (BioModels ID: MODEL1001200000) and at <url>http://www.mcisb.org/resources/genomescale/</url>.</p

A robust and efficient method for estimating enzyme complex abundance and metabolic flux from expression data

A major theme in constraint-based modeling is unifying experimental data, such as biochemical information about the reactions that can occur in a system or the composition and localization of enzyme complexes, with highthroughput data including expression data, metabolomics, or DNA sequencing. The desired result is to increase predictive capability resulting in improved understanding of metabolism. The approach typically employed when only gene (or protein) intensities are available is the creation of tissue-specific models, which reduces the available reactions in an organism model, and does not provide an objective function for the estimation of fluxes, which is an important limitation in many modeling applications. We develop a method, flux assignment with LAD (least absolute deviation) convex objectives and normalization (FALCON), that employs metabolic network reconstructions along with expression data to estimate fluxes. In order to use such a method, accurate measures of enzyme complex abundance are needed, so we first present a new algorithm that addresses quantification of complex abundance. Our extensions to prior techniques include the capability to work with large models and significantly improved run-time performance even for smaller models, an improved analysis of enzyme complex formation logic, the ability to handle very large enzyme complex rules that may incorporate multiple isoforms, and depending on the model constraints, either maintained or significantly improved correlation with experimentally measured fluxes. FALCON has been implemented in MATLAB and ATS, and can be downloaded from: https://github.com/bbarker/FALCON. ATS is not required to compile the software, as intermediate C source code is available, and binaries are provided for Linux x86-64 systems. FALCON requires use of the COBRA Toolbox, also implemented in MATLAB.Comment: 30 pages, 12 figures, 4 table

Quiescience as a mechanism for cyclical hypoxia and acidosis

Tumour tissue characteristically experiences fluctuations in substrate supply. This unstable microenvironment drives constitutive metabolic changes within cellular populations and, ultimately, leads to a more aggressive phenotype. Previously, variations in substrate levels were assumed to occur through oscillations in the hæmodynamics of nearby and distant blood vessels. In this paper we examine an alternative hypothesis, that cycles of metabolite concentrations are also driven by cycles of cellular quiescence and proliferation. Using a mathematical modelling approach, we show that the interdependence between cell cycle and the microenvironment will induce typical cycles with the period of order hours in tumour acidity and oxygenation. As a corollary, this means that the standard assumption of metabolites entering diffusive equilibrium around the tumour is not valid; instead temporal dynamics must be considered