Randomizing genome-scale metabolic networks
Areejit Samal and Olivier Martin
Contact the author: Please use for correspondence this email.
Submission date: 07. Dec. 2010 (revised version: July 2011)
published in: PLOS ONE, 6 (2011) 7, art-no. 22295
DOI number (of the published article): 10.1371/journal.pone.0022295
Keywords and phrases: complex networks, Flux balance analysis, MCMC
Download full preprint: PDF (740 kB)
Networks coming from protein-protein interactions, transcriptional regulation, signaling, or metabolism may appear to have “unusual” properties. To quantify this, it is appropriate to randomize the network and test the hypothesis that the network is not statistically different from expected in a motivated ensemble. However, when dealing with metabolic networks, the randomization of the network using edge exchange generates fictitious reactions that are biochemically meaningless. Here we provide several natural ensembles of randomized metabolic networks. A first constraint is to use valid biochemical reactions. Further constraints correspond to imposing appropriate functional constraints. We explain how to perform these randomizations with the help of Markov Chain Monte Carlo (MCMC) and show that they allow one to approach the properties of biological metabolic networks. The implication of the present work is that the observed global structural properties of real metabolic networks are likely to be the consequence of simple biochemical and functional constraints.