Summary: PySCeSToolbox is an extension to the Python Simulator for Cellular Systems (PySCeS) that includes tools for performing generalized supply-demand analysis, symbolic metabolic control analysis, and a framework for investigating the kinetic and thermodynamic aspects of enzyme-catalyzed reactions. Each tool addresses a different aspect of metabolic behaviour, control, and regulation; the tools complement each other and can be used in conjunction to better understand higher level system behaviour. Availability and implementation: PySCeSToolbox is available on Linux, Mac OS X and Windows. It is licensed under the BSD 3-clause licence. Code, setup instructions and a link to documentation can be found at https://github.com/PySCeS/PyscesToolbox. Contact: [email protected]. Supplementary information: Supplementary data are available at Bioinformatics online.

UNLABELLED: The Python Simulator for Cellular Systems (PySCeS) is an extendable research tool for the numerical analysis and investigation of cellular systems. AVAILABILITY: PySCeS is distributed as Open Source Software under the GNU General Public Licence and is available for download from http://pysces.sourceforge.net CONTACT: [email protected].

MOTIVATION: Molecular biotechnology now makes it possible to build elaborate systems models, but the systems biology community needs information standards if models are to be shared, evaluated and developed cooperatively. RESULTS: We summarize the Systems Biology Markup Language (SBML) Level 1, a free, open, XML-based format for representing biochemical reaction networks. SBML is a software-independent language for describing models common to research in many areas of computational biology, including cell signaling pathways, metabolic pathways, gene regulation, and others. AVAILABILITY: The specification of SBML Level 1 is freely available from http://www.sbml.org/

MOTIVATION: Realistic simulation of the kinetic properties of metabolic pathways requires rate equations to be expressed in reversible form, because substrate and product elasticities are drastically different in reversible and irreversible reactions. This presents no special problem for reactions that follow reversible Michaelis-Menten kinetics, but for enzymes showing cooperative kinetics the full reversible rate equations are extremely complicated, and anyway in virtually all cases the full equations are unknown because sufficiently complete kinetic studies have not been carried out. There is a need, therefore, for approximate reversible equations that allow convenient simulation without violating thermodynamic constraints. RESULTS: We show how the irreversible Hill equation can be generalized to a reversible form, including effects of modifiers. The proposed equation leads to behaviour virtually indistinguishable from that predicted by a kinetic form of the Adair equation, despite the fact that the latter is a far more complicated equation. By contrast, a reversible form of the Monod-Wyman-Changeux equation that has sometimes been used leads to predictions for the effects of modifiers at high substrate concentration that differ qualitatively from those given by the Adair equation.

MetaModel is a user-friendly program for calculating steady-state fluxes and metabolite concentrations of metabolic systems on the IBM PC and compatible computers. For any steady state that is obtained, one can then calculate a matrix of elasticity coefficients at that steady state, or a matrix of control and response coefficients. It thus offers a simple way to calculate the control structure of a pathway: it provides not only an educational tool that allows the student to verify empirically the classic summation relationships of metabolic control analysis but also a research tool for addressing 'what if?' questions about the behaviour of metabolic systems. Results can not only be printed or stored in a file, but can also be written to a special file that can be read by popular spreadsheet programs, thereby giving access to rapid, flexible and powerful methods for subsequent analysis and plotting of these results.

Steady-state modelling and control analysis by means of computer simulation provides valuable insight into the behavior of metabolic pathways. This review, which is aimed at the newcomer to this field, discusses the objectives of steady-state modelling and the steady-state properties of the four basic metabolic structures, namely linear and branched chains, loops and cycles. It is shown how the model definition in terms of stoichiometric reactions and rate equations leads to a set of balance equations from which the conservation constraints and flux relationships can be deduced, either informally or through a rigorous analysis of the stoichiometric matrix. The initial analysis of a steady-state metabolic model is summarized in an algorithm. Key references to the literature on metabolic modelling are given.

METAMOD, a BBC microcomputer-based software package for steady-state modelling and control analysis of model metabolic pathways, is described, The package consists of two programs. METADEF allows the user to define the pathway in terms of reactions, rate equations and initial concentrations of metabolites. METACAL uses one of two algorithms to calculate the steady-state concentrations and fluxes. One algorithm uses the current ratio of production and consumption rates of variable metabolites to adjust iteratively their concentrations in such a way that they converge towards the steady state. The other algorithm solves the roots of the system equations by means of a quasi-Newtonian procedure. Control analysis allows the calculation of elasticity, control and response coefficients, by means of finite difference approximation. METAMOD is interactive and easy to use, and suitable for teaching and research purposes.