Delving deeper: Relating the behaviour of a metabolic system to the properties of its components using symbolic metabolic control analysis

High-level behaviour of metabolic systems results from the properties of, and interactions between, numerous molecular components. Reaching a complete understanding of metabolic behaviour based on the system’s components is therefore a difficult task. This problem can be tackled by constructing and subsequently analysing kinetic models of metabolic pathways since such models aim to capture all the relevant properties of the system components and their interactions. Symbolic control analysis is a framework for analysing pathway models in order to reach a mechanistic understanding of their behaviour. By providing algebraic expressions for the sensitivities of system properties, such as metabolic fluxor steady-state concentrations, in terms of the properties of individual reactions it allows one to trace the high level behaviour back to these low level components. Here we apply this method to a model of pyruvate branch metabolism in Lactococcus lactis in order to explain a previously observed negative flux response towards an increase in substrate concentration. With this method we are able to show, first, that the sensitivity of flux towards changes in reaction rates (represented by flux control coefficients) is determined by the individual metabolic branches of the pathway, and second, how the sensitivities of individual reaction rates towards their substrates (represented by elasticity coefficients) contribute to this flux control. We also quantify the contributions of enzyme binding and mass-action to enzyme elasticity separately, which allows for an even finer-grained understanding of flux control. These analytical tools allow us to analyse the control properties of a metabolic model and to arrive at a mechanistic understanding of the quantitative contributions of each of the enzymes to this control. Our analysis provides an example of the descriptive power of the general principles of symbolic control analysis. Author summary Metabolic networks are complex systems consisting of numerous individual molecular components. The properties of these components, together with their non-linear interactions, give rise to high-level observed behaviour of the system in which they reside. Therefore, in order to fully understand the behaviour of a metabolic system, one has to consider the properties of all of its components. The analysis of computer models that capture and represent these systems and their components simplifies this task by allowing for an easy way to isolate the effects of each individual component. In this paper we use the framework of symbolic control analysis to investigate the sensitivity of the rate of flow of matter through one of the branches in a particular metabolic pathway towards changes in the rates of individual reactions. Here we are able to quantify how certain chains of reactions, individual reactions, and even thermodynamic and kinetic aspects of individual reactions contribute to the overall sensitivity of the rate of matter-flow. Thus, we are able to trace the behaviour of the system as a whole in a mechanistic way to the properties of the individual molecular components.