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

Further developments towards a genome-scale metabolic model of yeast

BACKGROUND: To date, several genome-scale network reconstructions have been used to describe the metabolism of the yeast Saccharomyces cerevisiae, each differing in scope and content. The recent community-driven reconstruction, while rigorously evidenced and well annotated, under-represented metabolite transport, lipid metabolism and other pathways, and was not amenable to constraint-based analyses because of lack of pathway connectivity. RESULTS: We have expanded the yeast network reconstruction to incorporate many new reactions from the literature and represented these in a well-annotated and standards-compliant manner. The new reconstruction comprises 1102 unique metabolic reactions involving 924 unique metabolites - significantly larger in scope than any previous reconstruction. The representation of lipid metabolism in particular has improved, with 234 out of 268 enzymes linked to lipid metabolism now present in at least one reaction. Connectivity is emphatically improved, with more than 90% of metabolites now reachable from the growth medium constituents. The present updates allow constraint-based analyses to be performed; viability predictions of single knockouts are comparable to results from in vivo experiments and to those of previous reconstructions. CONCLUSIONS: We report the development of the most complete reconstruction of yeast metabolism to date that is based upon reliable literature evidence and richly annotated according to MIRIAM standards. The reconstruction is available in the Systems Biology Markup Language (SBML) and via a publicly accessible database http://www.comp-sys-bio.org/yeastnet/

Integrating food webs with metabolic networks: modeling contaminant degradation in marine ecosystems

Georg Basler is supported by a Marie Curie Intra European Fellowship within the 7th European Community Framework Programme, ERC grant agreement number 329682. Evangelos Simeonidis is supported by the Luxembourg Centre for Systems Biomedicine.Peer reviewedPeer Reviewe

Optimisation of biochemical systems

EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Flux balance analysis: A geometric perspective

Advances in the field of bioinformatics have led to reconstruction of genome-scale networks for a number of key organisms. The application of physicochemical constraints to these stoichiometric networks allows researchers, through methods such as flux balance analysis, to highlight key sets of reactions necessary to achieve particular objectives. The key benefits of constraint-based analysis lie in the minimal knowledge required to infer systemic properties. However, network degeneracy leads to a large number of flux distributions that satisfy any objective; moreover, these distributions may be dominated by biologically irrelevant internal cycles. By examining the geometry underlying the problem, we define two methods for finding a unique solution within the space of all possible flux distributions; such a solution contains no internal cycles, and is representative of the space as a whole. The first method draws on typical geometric knowledge, but cannot be applied to large networks because of the high computational complexity of the problem. Thus a second method, an iteration of linear programs which scales easily to the genome scale, is defined. The algorithm is run on four recent genome-scale models, and unique flux solutions are found. The algorithm set out here will allow researchers in flux balance analysis to exchange typical solutions to their models in a reproducible format. Moreover, having found a single solution, statistical analyses such as correlations may be performed

Computing life: Add logos to biology and bios to physics

This paper discusses the interrelations between physics and biology. Particularly, we analyse the approaches for reconstructing the emergent properties of physical or biological systems. We propose approaches to scale emergence according to the degree of state-dependency of the system's component properties. Since the component properties of biological systems are state-dependent to a high extent, biological emergence should be considered as very strong emergence – i.e. its reconstruction would require a lot of information about state-dependency of its component properties. However, due to its complexity and volume, this information cannot be handled in the naked human brain, or on the back of an envelope. To solve this problem, biological emergence can be reconstructed in silico based on experimentally determined rate laws and parameter values of the living cell. According to some rough calculations, the silicon human might comprise the mathematical descriptions of around 105 interactions. This is not a small number, but taking into account the exponentially increase of computational power, it should not prove to be our principal limitation. The bigger challenges will be located in different areas. For example they may be related to the observer effect – the limitation to measuring a system's component properties without affecting the system. Another obstacle may be hidden in the tradition of "shaving away" all “unnecessary” assumptions (the so-called Occam's razor) that, in fact, reflects the intention to model the system as simply as possible and thus to deem the emergence to be less strong than it possibly is. We argue here that that Occam's razor should be replaced with the law of completeness