Modular metabolic control analysis of large responses

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 metabolic control analysis approach to introduce the study of systems in biochemistry: the glycolytic pathway in the red blood cell

Metabolic control analysis (MCA) is a promising approach in biochemistry aimed at understanding processes in a quantitative fashion. Here the contribution of enzymes and transporters to the control of a given pathway flux and metabolite concentrations is determined and expressed quantitatively by means of numerical coefficients. Metabolic flux can be influenced by a wide variety of modulators acting on one or more metabolic steps along the pathway. We describe a laboratory exercise to study metabolic regulation of human erythrocytes (RBCs). Within the framework of MCA, students use these cells to determine the sensitivity of the glycolytic flux to two inhibitors (iodoacetic acid: IA, and iodoacetamide: IAA) known to act on the enzyme glyceraldehyde-3-phosphate-dehydrogenase. Glycolytic flux was estimated by determining the concentration of extracellular lactate, the end product of RBC glycolysis. A low-cost colorimetric assay was implemented, that takes advantage of the straightforward quantification of the absorbance signal from the photographic image of the multi-well plate taken with a standard digital camera. Students estimate flux response coefficients for each inhibitor by fitting an empirical function to the experimental data, followed by analytical derivation of this function. IA and IAA exhibit qualitatively different patterns, which are thoroughly analyzed in terms of the physicochemical properties influencing their action on the target enzyme. IA causes highest glycolytic flux inhibition at lower concentration than IAA. This work illustrates the feasibility of using the MCA approach to study key variables of a simple metabolic system, in the context of an upper level biochemistry course. © 2018 International Union of Biochemistry and Molecular Biology, 46(5):502-515, 2018.

Control and regulation of pathways via negative feedback

The biochemical networks found in living organisms include a huge variety of control mechanisms at multiple levels of organization. While the mechanistic and molecular details of many of these control mechanisms are understood, their exact role in driving cellular behaviour is not. For example, yeast glycolysis has been studied for almost 80 years but it is only recently that we have come to understand the systemic role of the multitude of feedback and feed-forward controls that exist in this pathway. In this article, control theory is discussed as an approach to dissect the control logic of complex pathways. One of the key issues is distinguishing between the terms control and regulation and how these concepts are applied to regulated enzymes such as phosphofructokinase. In doing so, one of the paradoxes in metabolic regulation can be resolved where enzymes such as phosphofructokinase have little control but, nevertheless, possess significant regulatory influence.

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.

Optimization-based metabolic control analysis

In this work, a novel optimization-based metabolic control analysis (OMCA) method is introduced for reducing data requirement for metabolic control analysis (MCA). It is postulated that using the optimal control approach, the fluxes in a metabolic network are correlated to metabolite concentrations and enzyme activities as a state-feedback control system that is optimal with respect to a homeostasis objective. It is then shown that the optimal feedback gains are directly related to the elasticity coefficients (ECs) of MCA. This approach requires determination of the relative "importance" of metabolites and fluxes for the system, which is possible with significantly reduced experimental data, as compared with typical MCA requirements. The OMCA approach is applied to a top-down control model of glycolysis in hepatocytes. It is statistically demonstrated that the OMCA model is capable of predicting the ECs observed experimentally with few exceptions. Further, an OMCA-based model reconciliation study shows that the modification of four assumed stoichiometric coefficients in the model can explain most of the discrepancies, with the exception of elasticities with respect to the NADH/NAD ratio.

In silico strategy to rationally engineer metabolite production: A case study for threonine in Escherichia coli

Genetic engineering of metabolic pathways is a standard strategy to increase the production of metabolites of economic interest. However, such flux increases could very likely lead to undesirable changes in metabolite concentrations, producing deleterious perturbations on other cellular processes. These negative effects could be avoided by implementing a balanced increase of enzyme concentrations according to the Universal Method [Kacser and Acerenza (1993) Eur J Biochem 216:361-367]. Exact application of the method usually requires modification of many reactions, which is difficult to achieve in practice. Here, improvement of threonine production via pyruvate kinase deletion in Escherichia coli is used as a case study to demonstrate a partial application of the Universal Method, which includes performing sensitivity analysis. Our analysis predicts that manipulating a few reactions is sufficient to obtain an important increase in threonine production without major perturbations of metabolite concentrations.

Harnessing nature's toolbox: regulatory elements for synthetic biology

Synthetic biologists seek to engineer complex biological systems composed of modular elements. Achieving higher complexity in engineered biological organisms will require manipulating numerous systems of biological regulation: transcription; RNA interactions; protein signalling; and metabolic fluxes, among others. Exploiting the natural modularity at each level of biological regulation will promote the development of standardized tools for designing biological systems.

Stochastic control of metabolic pathways

We study the effect of extrinsic noise in metabolic networks. We introduce external random fluctuations at the kinetic level, and show how these lead to a stochastic generalization of standard metabolic control analysis. While summation and connectivity theorems hold true in the presence of extrinsic noise, control coefficients incorporate its effect through an explicit dependency on the noise intensity. New elasticities and response coefficients are also defined. Accordingly, the concept of control by noise is introduced as a way of tuning the systemic behaviour of metabolisms. We argue that this framework holds for intrinsic noise too, when time-scale separation is present in the system.

Kinetic properties required for sustained or paradoxical control of metabolic fluxes under large changes in enzyme activities

The effect that an increase in the activity of an enzyme has on its flux normally decreases with activity increase. To achieve a large increase in flux by manipulating a single step would therefore require a high initial effect that maintains or increases when the activity is increased, what has been called sustained or paradoxical control. Using metabolic control analysis for large responses, we derive conditions for sustained or paradoxical control in terms of elasticity coefficients. These are used to characterise types of rate laws contributing to this behaviour. The result that simple pathways, with normal kinetics, subject to large activity changes can lead to paradoxical control behaviour suggests that this type of pattern may be much more ubiquitous than could have, in principle, been suspected.