Modular metabolic control analysis of large responses in branched systems - application to aspartate metabolism

Control analysis in the identification of key enzymes driving metabolic adaptations: Towards drug target discovery

Metabolic Control Analysis (MCA) marked a turning point in understanding the design principles of metabolic network control by establishing control coefficients as a means to quantify the degree of control that an enzyme exerts on flux or metabolite concentrations. MCA has demonstrated that control of metabolic pathways is distributed among many enzymes rather than depending on a single rate-limiting step. MCA also proved that this distribution depends not only on the stoichiometric structure of the network but also on other kinetic determinants, such as the degree of saturation of the enzyme active site, the distance to thermodynamic equilibrium, and metabolite feedback regulatory loops. Consequently, predicting the alterations that occur during metabolic adaptation in response to strong changes involving a redistribution in such control distribution can be challenging. Here, using the framework provided by MCA, we illustrate how control distribution in a metabolic pathway/network depends on enzyme kinetic determinants and to what extent the redistribution of control affects our predictions on candidate enzymes suitable as targets for small molecule inhibition in the drug discovery process. Our results uncover that kinetic determinants can lead to unexpected control distribution and outcomes that cannot be predicted solely from stoichiometric determinants. We also unveil that the inference of key enzyme-drivers of an observed metabolic adaptation can be dramatically improved using mean control coefficients and ruling out those enzyme activities that are associated with low control coefficients. As the use of constraint-based stoichiometric genome-scale metabolic models (GSMMs) becomes increasingly prevalent for identifying genes/enzymes that could be potential drug targets, we anticipate that incorporating kinetic determinants and ruling out enzymes with low control coefficients into GSMM workflows will facilitate more accurate predictions and reveal novel therapeutic targets.

A modular modulation method for achieving increases in metabolite production

Increasing the production of overproducing strains represents a great challenge. Here, we develop a modular modulation method to determine the key steps for genetic manipulation to increase metabolite production. The method consists of three steps: (i) modularization of the metabolic network into two modules connected by linking metabolites, (ii) change in the activity of the modules using auxiliary rates producing or consuming the linking metabolites in appropriate proportions and (iii) determination of the key modules and steps to increase production. The mathematical formulation of the method in matrix form shows that it may be applied to metabolic networks of any structure and size, with reactions showing any kind of rate laws. The results are valid for any type of conservation relationships in the metabolite concentrations or interactions between modules. The activity of the module may, in principle, be changed by any large factor. The method may be applied recursively or combined with other methods devised to perform fine searches in smaller regions. In practice, it is implemented by integrating to the producer strain heterologous reactions or synthetic pathways producing or consuming the linking metabolites. The new procedure may contribute to develop metabolic engineering into a more systematic practice.

Regulatory design in a simple system integrating membrane potential generation and metabolic ATP consumption. Robustness and the role of energy dissipating processes

Bacterial physiological responses integrate energy-coupling processes at the membrane level with metabolic energy demand. The regulatory design behind these responses remains largely unexplored. Propionigenium modestum is an adequate organism to study these responses because it presents the simplest scheme known integrating membrane potential generation and metabolic ATP consumption. A hypothetical sodium leak is added to the scheme as the sole regulatory site. Allosteric regulation is assumed to be absent. Information of the rate equations is not available. However, relevant features of the patterns of responses may be obtained using Metabolic Control Analysis (MCA) and Metabolic Control Design (MCD). With these tools, we show that membrane potential disturbances can be compensated by adjusting the leak flux, without significant perturbations of ATP consumption. Perturbations of membrane potential by ATP demand are inevitable and also require compensatory changes in the leak. Numerical simulations were performed with a kinetic model exhibiting the responses for small changes obtained with MCA and MCD. A modest leak (10% of input) was assumed for the reference state. We found that disturbances in membrane potential and ATP consumption, produced by environmental perturbations of the cation concentration, may be reverted to the reference state adjusting the leak. Leak changes can also compensate for undesirable effects on membrane potential produced by changes in nutrient availability or ATP demand, in a wide range of values. The system is highly robust to parameter fluctuations. The regulatory role of energy dissipating processes and the trade-off between energetic efficiency and regulatory capacity are discussed.