Matrix method for determining steps most rate-limiting to metabolic fluxes in biotechnological processes

The metabolic control theory developed by Kacser, Burns, Heinrich, and Rapoport is briefly outlined, extended, and transformed so as optimally to address some biotechnological questions. The extensions include (i) a new theorem that relates the control of metabolite concentrations by enzyme activities to flux ratios at branches in metabolic pathways; (ii) a new theorem that does the same for the control of the distribution of the flux over two branches; (iii) a method that expresses these controls into properties (the so-called elasticity coefficients) of the enzymes in the pathway; and (iv) a theorem that relates the effects of changes in metabolite concentrations on reaction rates to the effects of changes in enzyme properties on the same rates. Matrix equations relating the flux control and concentration control coefficients to the elasticity coefficients of enzymes in simple linear and branched pathways incorporating feedback are given, together with their general solutions and a numerical example. These equations allow one to develop rigorous criteria by which to decide the optimal strategy for the improvement of a microbial process. We show how this could be used in deciding which property of which enzyme should be changed in order to obtain the maximal concentration of a metabolite or the maximal metabolic flux.     

Defining control coefficients in non-ideal metabolic pathways

The extent to which an enzyme controls a flux has been defined as the effect on that flux of a small modulation of the activity of that enzyme divided by the magnitude of the modulation. We here show that in pathways with metabolic channelling or high enzyme concentrations and conserved moieties involving both enzymic and non-enzymic species, this definition is ambiguous; the magnitude of the corresponding flux control coefficient depends on how the enzyme activity is modulated. This is illustrated with two models of biochemically relevant pathways, one in which dynamic metabolite channelling plays a role, and one with a moiety-conserved cycle. To avoid such ambiguity, we view biochemical pathways in a more detailed manner, i.e., as a network of elemental steps. We define 'elemental control coefficients' in terms of the effect on a flux of an equal modulation of the forward and reverse rate constant of any such elemental step (which may correspond to transitions between enzyme states). This elemental control coefficient is independent of the method of modulation. We show how metabolic control analysis can proceed when formulated in terms of the elemental control coefficients and how the traditional control coefficients are related to these elemental control coefficients. An 'impact' control coefficient is defined which quantifies the effect of an activation of all elemental processes in which an enzyme is involved. It equals the sum of the corresponding elemental control coefficients. In ideal metabolic pathways this impact control coefficient reduces to the traditional flux control coefficient. Differences between the traditional control coefficients are indicative of non-ideality of a metabolic pathway, i.e. of channelling or high enzyme concentrations.     

Getting to the inside of cells using metabolic control analysis

Metabolic control analysis can relate control properties of an intact system to kinetic properties (elasticity coefficients) of the enzymes within that system. The method formulating the former as matrix inverse of the latter is elaborated here for the general case and founded in standard metabolic control theory. Then a method is developed that accomplishes the reverse: it is shown that a matrix containing all elasticity coefficients and information concerning the pathway structure equals the inverse of a matrix containing flux and concentration control coefficients. As a consequence, by measuring the control properties of an intact system, one is able to deduce its in situ pathway structure and enzyme kinetic properties: This solves the ever-present question of whether the kinetic properties of enzymes in their isolated state differ from those under the conditions prevailing in the cell.     

Metabolic control and its analysis. Additional relationships between elasticities and control coefficients

Existing theorems from the analysis of metabolic control have been taken and embedded in a simple matrix algebra procedure for calculating the flux control coefficients of enzymes (formerly known as sensitivities) in a metabolic pathway from their kinetic properties (their elasticities). New theorems governing the flux control coefficients of branched pathways and substrate cycles have been derived to allow the procedure to be applied to complex pathway configurations. Modifications to the elasticity terms used in the equations have been theoretically justified so that the method remains valid for pathways with conserved metabolites (for example, the adenine nucleotide pool or the intermediates of a catalytic cycle such as the tricarboxylic acid cycle) or with pools of metabolites kept very near to equilibrium by very rapid reactions. The matrix equations generated using these theorems and relationships may be solved algebraically or numerically. Algebraic solutions have been used to determine the factors responsible for the degree of amplification of flux control coefficients by substrate cycles and to show that it is possible to derive expressions for the elasticities of a group of enzymes.     

Control of molecular transformations in polyenzyme systems: quantitative theory of the regulation of metabolism

An attempt of a comprehensive treatment of the theory of metabolic control is presented. The introductory section giving an outline of the early development of the theory, is followed by definitions quantifying the control in the metabolic system. By means of the perturbation method the complete system of equations is obtained which allows one to express all the enzyme control coefficients ("global" coefficients) through the elasticity coefficients characterizing kinetic properties of individual enzymes ("local" coefficients) and through the steady-state values of metabolic fluxes and concentrations. It is shown how connectivity relations between global and local coefficients should be modified when conserved sums of intermediates are present in the system. A new theorem is derived, it allows one to express the global response of the system to any change in the external parameter (such as external effector concentration, or temperature, pH, ionic strength, ets.) through the control coefficients and local responses of individual reaction steps. Explicit formulas are derived for response coefficients of the fluxes and concentrations to changes in the conserved sums of intermediates, which express the values of these global coefficients through the control and elasticity coefficients of enzymes and steady-state pools. The results obtained comprise as a special case all the results published so far in the literature.     

'Slave' metabolites and enzymes. A rapid way of delineating metabolic control.

Although control of fluxes and concentrations tends to be distributed rather than confined to a single rate-limiting enzyme, the extent of control can differ widely between enzymes in a metabolic network. In some cases, there are enzymes that lack control completely. This paper identifies one surprising origin of such lack of control: If, in a metabolic system, there is a metabolite that affects the catalytic rate of only one enzyme, the corresponding enzyme cannot control any metabolic variable other than the concentration of that metabolite. We call such enzymes 'slave enzymes', and the corresponding metabolites 'slave metabolites'. Implications of the existence of slave enzymes for the control properties of enzymes further down the metabolic pathway are discussed and examined for the glycolytic pathway of yeast. Inadvertent assumptions in metabolic models may cause the latter incorrectly to calculate absence of metabolic control. The phenomenon of slave enzymes may well be important in enhancing metabolic signal transduction.     

Control of gluconeogenesis in rat liver cells. Flux control coefficients of the enzymes in the gluconeogenic pathway in the absence and presence of glucagon.

We have used control analysis to quantify the distribution of control in the gluconeogenic pathway in liver cells from starved rats. Lactate and pyruvate were used as gluconeogenic substrates. The flux control coefficients of the various enzymes in the gluconeogenic pathway were calculated from the elasticity coefficients of the enzymes towards their substrates and products and the fluxes through the different branches in the pathway. The elasticity coefficients were either calculated from gamma/Keq. ratios (where gamma is the mass-action ratio and Keq. is the equilibrium constant) and enzyme-kinetic data or measured experimentally. It is concluded that the gluconeogenic enzyme pyruvate carboxylase and the glycolytic enzyme pyruvate kinase play a central role in control of gluconeogenesis. If pyruvate kinase is inactive, gluconeogenic flux from lactate is largely controlled by pyruvate carboxylase. The low elasticity coefficient of pyruvate carboxylase towards its product oxaloacetate minimizes control by steps in the gluconeogenic pathway located after pyruvate carboxylase. This situation occurs when maximal gluconeogenic flux is required, i.e. in the presence of glucagon. In the absence of the hormone, when pyruvate kinase is active, control of gluconeogenesis is distributed among many steps, including pyruvate carboxylase, pyruvate kinase, fructose-1,6-bisphosphatase and also steps outside the classic gluconeogenic pathway such as the adenine-nucleotide translocator.     

Total cell protein concentration as an evolutionary constraint on the metabolic control distribution in cells.

Cell protein occupies 15-35% of cell volume. This level is argued to be the maximum compatible with cell function. Because of this constraint, selection pressure during evolution is likely to have maximized pathway fluxes for minimum total protein level. Pathways optimized in this way are shown to have the following characteristics: (1) the "simple" flux control coefficients of all enzymes are equal, (2) the normal flux control coefficients depend on the relative kinetic constants of the enzymes, such that enzymes with low specific activity are present at relatively high levels and have high flux control, (3) the normal flux control coefficients are proportional to enzyme levels. A single rate limiting step located at the first step in a pathway is likely to be inefficient in terms of protein levels, and the major metabolic pathways are therefore expected to have control distributed throughout the pathway. This has important implications for metabolic control.     

CONTROL: software for the analysis of the control of metabolic networks

The program CONTROL is based on metabolic control theory anduses the method developed by Reder (1988). In this theory, twosets of parameters are defined in the vicinity of a steady-state:the elasticity coefficients which describe the local behaviourof the isolated enzymes, and the control coefficients whichexpress the response of the whole metabolic network to perturbationsat a given step. The theory shows that relationships exist betweenthe control coefficients (summation relationships or structuralrelationships) and also between the two types of coefficients(control and elasticity coefficients: connectivity relationships).The program CONTROL is divided into two parts (sub-menus). Thefirst one calculates all the control coefficients (flux andconcentrations) of a metabolic network from the elasticity coefficients.Using the second menu, the symbolic relationships are obtainedbetween the control coefficients (summation relationships) andbetween the control coefficients and the elasticity coefficients(connectivity relationships). These two sub-menus can be appliedindependently to any metabolic network (to date limited to 19steps and 19 metabolites).     

Subtleties in control by metabolic channelling and enzyme organization

Because of its importance to cell function, the free-energy metabolism of the living cell is subtly and homeostatically controlled. Metabolic control analysis enables a quantitative determination of what controls the relevant fluxes. However, the original metabolic control analysis was developed for idealized metabolic systems, which were assumed to lack enzyme-enzyme association and direct metabolite transfer between enzymes (channelling). We here review the recently developed molecular control analysis, which makes it possible to study non-ideal (channelled, organized) systems quantitatively in terms of what controls the fluxes, concentrations, and transit times. We show that in real, non-ideal pathways, the central control laws, such as the summation theorem for flux control, are richer than in ideal systems: the sum of the control of the enzymes participating in a non-ideal pathway may well exceed one (the number expected in the ideal pathways), but may also drop to values below one. Precise expressions indicate how total control is determined by non-ideal phenomena such as ternary complex formation (two enzymes, one metabolite), and enzyme sequestration. The bacterial phosphotransferase system (PTS), which catalyses the uptake and concomitant phosphorylation of glucose (and also regulates catabolite repression) is analyzed as an experimental example of a non-ideal pathway. Here, the phosphoryl group is channelled between enzymes, which could increase the sum of the enzyme control coefficients to two, whereas the formation of ternary complexes could decrease the sum of the enzyme control coefficients to below one. Experimental studies have recently confirmed this identification, as well as theoretically predicted values for the total control. Macromolecular crowding was shown to be a major candidate for the factor that modulates the non-ideal behaviour of the PTS pathway and the sum of the enzyme control coefficients.     

Modern theories of metabolic control and their applications

Existing, qualitative notions with respect to the way in which enzyme properties control metabolism are discussed in the light of the control analysis developed by H. Kacser and J. A. Burns ((1973) in: Rate Control of Biological Processes, Davies DD, ed., Cambridge University Press, pp. 63–104) and R. Heinrich and T. A. Rapoport ((1974) Eur. 3. Biochem.42, 89–95), and recent experimental data. Points at which the existing notions should be adjusted are: (i) Metabolic control is shared by enzymes rather than confined to one rate-limiting enzyme per pathway. (if) Whether an enzyme exercises strong control on a flux cannot be deduced solely from its own properties, nor is it directly related to its distance from equilibrium. With respect to metabolic control, enzymes should be classified into four groups, rather than two (reversible versus irreversible). (iii) The distribution of control among the enzymes depends on the metabolic conditions. (iv) Control structures of metabolic pathways probably differ with the function of that pathway.     

A 'top-down' approach to the determination of control coefficients in metabolic control theory

A new approach to the determination of flux and concentration control coefficients in metabolic pathways is outlined. Linear pathways are conceptually divided in two around an intermediate metabolite (or group or metabolites) and the control coefficients of the two parts are derived from the elasticity coefficients of the two parts to the intermediate. Branched pathways are treated similarly, the control coefficients of the branches being derived either from the elasticities of the branches to their common intermediate or from the relative flux changes of the branches. Repeating this analysis around other intermediates in the pathway allows the control coefficients of smaller and smaller groups of enzymes to be determined. In complex systems this approach to describing control may have several advantages over determining the control coefficients of individual enzymes and is a potentially useful complementary approach.     

Fluxes and metabolic pools as model traits for quantitative genetics. I. The L-shaped distribution of gene effects

The fluxes through metabolic pathways can be considered as model quantitative traits, whose QTL are the polymorphic loci controlling the activity or quantity of the enzymes. Relying on metabolic control theory, we investigated the relationships between the variations of enzyme activity along metabolic pathways and the variations of the flux in a population with biallelic QTL. Two kinds of variations were taken into account, the variation of the average enzyme activity across the loci, and the variation of the activity of each enzyme of the pathway among the individuals of the population. We proposed analytical approximations for the flux mean and variance in the population as well as for the additive and dominance variances of the individual QTL. Monte Carlo simulations based on these approximations showed that an L-shaped distribution of the contributions of individual QTL to the flux variance (R(2)) is consistently expected in an F(2) progeny. This result could partly account for the classically observed L-shaped distribution of QTL effects for quantitative traits. The high correlation we found between R(2) value and flux control coefficients variance suggests that such a distribution is an intrinsic property of metabolic pathways due to the summation property of control coefficients.     

Dominance is not inevitable

In a diploid organism, a mutant gene that results in elimination of an enzyme activity in the homozygote is almost universally found to be recessive, so that the heterozygote phenotype is virtually indistinguishable from the wild type. It has been argued (H. Kacser & J. A. Burns, Genetics 97, 639-666 (1981)) that there is no need to look to evolution for an explanation of this phenomenon, as it is an inevitable consequence of the low control coefficients for metabolic flux possessed by nearly all enzymes. However, it is possible to envisage pathways in which every enzyme is more than half-saturated, so that moderate changes in the concentration of any enzyme result in substantial changes in metabolic flux. Such behaviour can occur, for example, if the limiting rates of the enzymes decrease as one proceeds along the pathway and the precursor concentration is large compared with the Michaelis constants of all the enzymes. Consequently one does require an explanation in terms of natural selection of why such pathways are apparently not observed in nature.     

Metabolic control analysis of Aspergillus niger L-arabinose catabolism

A mathematical model of the L-arabinose/D-xylose catabolic pathway of Aspergillus niger was constructed based on the kinetic properties of the enzymes. For this purpose L-arabinose reductase, L-arabitol dehydrogenase and D-xylose reductase were purified using dye-affinity chromatography, and their kinetic properties were characterized. For the other enzymes of the pathway the kinetic data were available from the literature. The metabolic model was used to analyze flux and metabolite concentration control of the L-arabinose catabolic pathway. The model demonstrated that flux control does not reside at the enzyme following the intermediate with the highest concentration, L-arabitol, but is distributed over the first three steps in the pathway, preceding and following L-arabitol. Flux control appeared to be strongly dependent on the intracellular L-arabinose concentration. At 5 mM intracellular L-arabinose, a level that resulted in realistic intermediate concentrations in the model, flux control coefficients for L-arabinose reductase, L-arabitol dehydrogenase and L-xylulose reductase were 0.68, 0.17 and 0.14, respectively. The analysis can be used as a guide to identify targets for metabolic engineering aiming at either flux or metabolite level optimization of the L-arabinose catabolic pathway of A. niger. Faster L-arabinose utilization may enhance utilization of readily available organic waste containing hemicelluloses to be converted into industrially interesting metabolites or valuable enzymes or proteins.     

Theory of metabolism regulation: a complete system of equations for regulation coefficients

Basic quantitative parameters of control in a metabolic system are considered: control coefficients of enzymes with respect to metabolic fluxes and concentrations, and in the case when there are conservation laws, the response coefficients of metabolic fluxes and concentrations to changes in the conserved sums of metabolite concentrations (e. g. conserved moieties). Relationships are obtained which generalize the well known connectivity relations for the case of metabolites binding by conservation laws. Additional relationships are obtained which complement the set of connectivity relations up to the complete system of equations for determining all the control coefficients. The control coefficients are expressed through the enzyme elasticity coefficients, steady state metabolic fluxes and concentrations. Formulas are derived which express response coefficients of flux and concentrations through the enzyme control and elasticity coefficients and metabolite concentrations.     

Topological analysis of metabolic control

A topological approach is presented for the analysis of control and regulation in metabolic pathways. In this approach, the control structure of a metabolic pathway is represented by a weighted directed graph. From an inspection of the topology of the graph, the control coefficients of the enzymes are evaluated in a heuristic manner in terms of the enzyme elasticities. The major advantage of the topological approach is that it provides a visual framework for (1) calculating the control coefficients of the enzymes, (2) analyzing the cause-effect relationships of the individual enzymes, (3) assessing the relative importance of the enzymes in metabolic regulation, and (4) simplifying the structure of a given pathway, from a regulatory viewpoint. Results are obtained for (a) an unbranched pathway in the absence of feedback the feedforward regulation and (b) an unbranched pathway with feedback inhibition. Our formulation is based on the metabolic control theory of Kacser and Burns (1973) and Heinrich and Rapoport (1974).     

Steady-state analysis of metabolic pathways: comparing the double modulation method and the lin-log approach

The increasing interest in studying enzyme kinetics under in vivo conditions requires practical methods to estimate control parameters from experimental data. In contrast to currently established approaches of dynamic modelling, this paper addresses the steady-state analysis of metabolic pathways. Within the framework of metabolic control analysis (MCA), elasticity coefficients are used to describe the control properties of a local enzyme reaction. The double modulation method is one of the first experimental approaches to estimate elasticity coefficients from measurements of steady-state flux rates and metabolite concentrations. We propose a generalized form of the double modulation method and compare it to the recently developed linear-logarithmic approach.