A living organism must not only organize itself from within; it must also maintain its organization in the face of changes in its environment and degradation of its components. We show here that a simple (M,R)-system consisting of three interlocking catalytic cycles, with every catalyst produced by the system itself, can both establish a non-trivial steady state and maintain this despite continuous loss of the catalysts by irreversible degradation. As long as at least one catalyst is present at a sufficient concentration in the initial state, the others can be produced and maintained. The system shows bistability, because if the amount of catalyst in the initial state is insufficient to reach the non-trivial steady state the system collapses to a trivial steady state in which all fluxes are zero. It is also robust, because if one catalyst is catastrophically lost when the system is in steady state it can recreate the same state. There are three elementary flux modes, but none of them is an enzyme-maintaining mode, the entire network being necessary to maintain the two catalysts.

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.