We describe the construction and analysis of a genome-scale metabolic model of Arabidopsis (Arabidopsis thaliana) primarily derived from the annotations in the Aracyc database. We used techniques based on linear programming to demonstrate the following: (1) that the model is capable of producing biomass components (amino acids, nucleotides, lipid, starch, and cellulose) in the proportions observed experimentally in a heterotrophic suspension culture; (2) that approximately only 15% of the available reactions are needed for this purpose and that the size of this network is comparable to estimates of minimal network size for other organisms; (3) that reactions may be grouped according to the changes in flux resulting from a hypothetical stimulus (in this case demand for ATP) and that this allows the identification of potential metabolic modules; and (4) that total ATP demand for growth and maintenance can be inferred and that this is consistent with previous estimates in prokaryotes and yeast.

We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

ScrumPy is a software package used for the definition and analysis of metabolic models. It is written using the Python programming language that is also used as a user interface. ScrumPy has features for both kinetic and structural modelling, but the emphasis is on structural modelling and those features of most relevance to analysis of large (genome-scale) models. The aim is at describing ScrumPy's functionality to readers with some knowledge of metabolic modelling, but implementation, programming and other computational details are omitted. ScrumPy is released under the Gnu Public Licence, and available for download from http://mudshark.brookes.ac.uk/ScrumPy.