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 analysis of mammalian serine biosynthesis. Feedback inhibition on the final step.

The flux of serine biosynthesis in the liver of the normal rabbit, and of the rat on a low protein diet, is most sensitive to the activity of phosphoserine phosphatase (flux control coefficient up to 0.97), the last of the three enzymes in the pathway after it branches from glycolysis. The concentration of the pathway product, serine, has a strong controlling influence on the flux (response coefficient up to -0.64) through feedback inhibition at this step. The pathway is therefore controlled primarily by the demand for serine rather than the supply of the pathway precursor, 3-phosphoglycerate. Under conditions where there is a lower biosynthetic flux, the flux control coefficients of the first two enzymes of the pathway are increased, and are probably dominant in the rat on a normal diet. In rabbit liver, when ethanol is used to inhibit serine biosynthesis, control can be distributed between the three enzymes, even though the reactions catalysed by the first two remain close to equilibrium. Apart from their intrinsic value in aiding the understanding of the regulation of mammalian serine metabolism, our findings illustrate the danger of assuming that there are invariant design principles in the regulation of metabolic pathways, such as feedback control on the first step after a branch.     

Metabolism of hydroxybiphenyl and chloro-hydroxybiphenyl by biphenyl/chlorobiphenyl degradingPseudomonas testosteroni,strain B-356

Summary A biphenyl (BP) and chlorobiphenyl (CBP) metabolizingPseudomonas testosteroni, strain B-356 was also capable of utilizing 2-, 3-, and 4-hydroxybiphenyl. Data presented here suggest that utilization of biphenyl and mono-subtituted biphenyls involves the enzymes of the same pathway. Chloro-hydroxybiphenyls were also metabolized by strain B-356. The unsubstituted ring is first hydroxylated in position 2 and 3 and then cleaved in ameta 1, and 2, position to ultimately generate the benzoic acid derivatives. Since strain B-356 was capable of utilizing benzoic acid and mono-hydroxybenzoic acids, the utilization of biphenyl, 2-, 3-, and 4-hydroxybiphenyl is complete at non-toxic concentrations of the substrates. Chlorobenzoic acids and chloro-hydroxybenzoic acids were not metabolized further by this strain. Studies usingPseudomonas putida, strain KT2440 carrying cloned BP/CBP genes from strain B-356 provided further evidence for the presence of a common pathway for the metabolism of the above compounds inP. testosteroni, strain B-356. Suggestions are made on significance of the broad substrate specificity of the enzymes of biphenyl/chlorobiphenyl pathway in regard to their possible origin and in relation to PCB mixture degradation.     

Control and regulation of metabolic fluxes in microbes by substrates and enzymes

The control of enzymes and substrates on the flux through microbial metabolic pathways can be quantified in terms of flux control coefficients. In pathways involving group transfer, the summation theorem for flux control by the enzymes has to be modified: the sum of control by all enzymes is between 1 and 2. The phosphoenolpyruvate:glucose phosphotransferase system is such a pathway. Experimental determination of the control by the enzymes in this pathway is under way. The control of the enzymes on the glycolytic flux in yeast is low, with the possible exception of the uptake step. InKlebsiella pneumoniae potassium and ammonium ions can simultaneously be limiting, (i.e. have significant control on growth) at pH 6, but not at pH 8. This may be due to the fact that at pH 8 the high-affinity potassium uptake system is absent.     

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.     

Localization and organization of phenol degradation genes ofPseudomonas putida strain H

The genetic organization of the DNA region encoding the phenol degradation pathway ofPseudomonas putida H has been investigated. This strain can utilize phenol or some of its methylated derivatives as its sole source of carbon and energy. The first step in this process is the conversion of phenol into catechol. Catechol is then further metabolized via themeta-cleavage pathway into TCA cycle intermediates. Genes encoding these enzymes are clustered on the plasmid pPGH1. A region of contiguous DNA spanning about 16 kb contains all of the genetic information necessary for inducible phenol degradation. The analysis of mutants generated by insertion of transposons and cassettes indicates that all of the catabolic genes are contained in a single operon. This codes for a multicomponent phenol hydroxylase andmeta-cleavage pathway enzymes. Catabolic genes are subject to positive control by the gene product(s) of a second locus.     

Has the Bacterial Biphenyl Catabolic Pathway Evolved Primarily To Degrade Biphenyl? The Diphenylmethane Case

In this work, we have compared the ability of Pandoraea pnomenusa B356 and of Burkholderia xenovorans LB400 to metabolize diphenylmethane and benzophenone, two biphenyl analogs in which the phenyl rings are bonded to a single carbon. Both chemicals are of environmental concern. P. pnomenusa B356 grew well on diphenylmethane. On the basis of growth kinetics analyses, diphenylmethane and biphenyl were shown to induce the same catabolic pathway. The profile of metabolites produced during growth of strain B356 on diphenylmethane was the same as the one produced by isolated enzymes of the biphenyl catabolic pathway acting individually or in coupled reactions. The biphenyl dioxygenase oxidizes diphenylmethane to 3-benzylcyclohexa-3,5-diene-1,2-diol very efficiently, and ultimately this metabolite is transformed to phenylacetic acid, which is further metabolized by a lower pathway. Strain B356 was also able to cometabolize benzophenone through its biphenyl pathway, although in this case, this substrate was unable to induce the biphenyl catabolic pathway and the degradation was incomplete, with accumulation of 2-hydroxy-6,7-dioxo-7-phenylheptanoic acid. Unlike strain B356, B. xenovorans LB400 did not grow on diphenylmethane. Its biphenyl pathway enzymes metabolized diphenylmethane, but they poorly metabolize benzophenone. The fact that the biphenyl catabolic pathway of strain B356 metabolized diphenylmethane and benzophenone more efficiently than that of strain LB400 brings us to postulate that in strain B356, this pathway evolved divergently to serve other functions not related to biphenyl degradation.     

THE EVOLUTION OF GENES IN BRANCHED METABOLIC PATHWAYS

Simulation models of the evolution of genes in a branched metabolic pathway subject to stabilizing selection on flux are described and analyzed. The models are based either on metabolic control theory (MCT), with the assumption that enzymes are far from saturation, or on Michaelis–Menten kinetics, which allows for saturation and near saturation. Several predictions emerge from the models: (1) flux control evolves to be concentrated at pathway branch points, including the first enzyme in the pathway. (2) When flux is far from its optimum, adaptive substitutions occur disproportionately often in branching enzymes. (3) When flux is near its optimum, adaptive substitutions occur disproportionately often in nonbranching enzymes. (4) Slightly deleterious substitutions occur disproportionately often in nonbranching enzymes. (5) In terms of both flux control and patterns of substitution, pathway branches are similar to those predicted for linear pathways. These predictions provide null hypotheses for empirical examination of the evolution of genes in metabolic pathways.     

大肠杆菌苏氨酸合成途径动力学模型的构建与分析

动力学模型分析有利于理解生物系统的调控机制,从而为高效细胞工厂的理性设计提供指导。基于以往发表的相关途径动力学模型和测量的酶动力学数据,开发了大肠杆菌苏氨酸合成途径的动力学模型。模型包含从天冬氨酸至苏氨酸的合成途径及葡萄糖开始的为合成途径提供前体以及能量的代谢途径。与以往模型不同的是新模型中考虑了能量和还原力的平衡,从而使模型模拟的系统自身成为一个不需要从外界提供能量和还原力的自洽系统。模型稳态分析的结果表明PTS、G6PDH和HDH等反应对苏氨酸合成反应的通量控制系数较大,通过过表达这些反应的酶可以有效增加苏氨酸合成反应的通量。     

Regulation and control of compartmentalized glycolysis in bloodstream formTrypanosoma brucei

Unlike other eukaryotic cells, trypanosomes possess a compartmentalized glycolytic pathway. The conversion of glucose into 3-phosphoglycerate takes place in specialized peroxisomes, called glycosomes. Further conversion of this intermediate into pyruvate occurs in the cytosol. Due to this compartmentation, many regulatory mechanisms operating in other cell types cannot work in trypanosomes. This is reflected by the insensitivity of the glycosomal enzymes to compounds that act as activity regulators in other cell types. Several speculations have been raised about the function of compartmentation of glycolysis in trypanosomes. We calculate that even in a noncompartmentalized trypanosome the flux through glycolysis should not be limited by diffusion. Therefore, the sequestration of glycolytic enzymes in an organelle may not serve to overcome a diffusion limitation. We also search the available data for a possible relation between compartmentation and the distribution of control of the glycolytic flux among the glycolytic enzymes. Under physiological conditions, the rate of glycolytic ATP production in the bloodstream form of the parasite is possibly controlled by the oxygen tension, but not by the glucose concentration. Within the framework of Metabolic Control Analysis, we discuss evidence that glucose transport, although it does not qualify as the sole rate-limiting step, does have a high flux control coefficient. This, however, does not distinguish trypanosomes from other eukaryotic cell types without glycosomes.     

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.     

Common induction and regulation of biphenyl, xylene/toluene, and salicylate catabolism in Pseudomonas paucimobilis.

A strain of Pseudomonas paucimobilis (strain Q1) capable of utilizing biphenyl was isolated from soil. This strain grew not only on substituted biphenyls, but also on salicylate, xylene or toluene or both (xylene/toluene), and substituted benzoates. Evidence is presented that the catabolism of biphenyl, xylene/toluene, and salicylate is regulated by a common unit in this strain. The catabolism of biphenyl, xylene/toluene, and salicylate is interrelated, since benzoate and toluate are common metabolic intermediates of biphenyl and xylene/toluene, and salicylate is produced from 2-hydroxybiphenyl (o-phenylphenol). All the oxidative enzymes of the biphenyl, xylene/toluene, and salicylate degradative pathways were induced when the cells were grown on either biphenyl, xylene/toluene or salicylate. The P. paucimobilis Q1 cells showed induction of the meta-cleavage enzymes of both 2,3-dihydroxybiphenyl and catechol. Biphenyl-negative derivatives of strain Q1 were simultaneously rendered xylene/toluene and salicylate negative, whereas reversion to the biphenyl-positive character of such derivatives invariably led to a xylene/toluene- and salicylate-positive phenotype. Growth of the P. paucimobilis Q1 cells with benzoate as a sole carbon source allowed the induction of only the ortho pathway enzymes, suggesting that biphenyl, xylene/toluene, or salicylate specifically induced the meta pathway enzymes for the oxidative degradation of these compounds.     

Control of the flux in the arginine pathway of Neurospora crassa. Modulations of enzyme activity and concentration.

The influence of particular enzyme activities on the flux of metabolites in a pathway can be estimated by 'modulating' enzymes (i.e. changing turnover or concentration) and measuring the response in various parts of the system. By controlling the nuclear ration of two genetically different nuclear types in heterokaryons, the enzyme concentrations at four different steps in the arginine pathway were decreased over a range. This range was extended by the use of bradytrophs, mutant strains specifying enzymes with greatly diminished enzyme activities. Strains altered simultaneously at more than one step were also constructed by genetic recombination. By measuring the outputs of the pathway and the steady-state concentrations of intermediate pools, the fluxes in different parts of the pathway were calculated. This allowed the construction of flux/enzyme relationships, the slope of which is a measure of the sensitivity of a flux to the change in enzyme activity at that step. All fluxes were found to be considerably buffered for quite substantial decreases in the activities of all enzymes. Mass action plays an important part in this phenomenon, as do inhibition and repression. Because of the existence of expansion fluxes in growing systems, we find quantitatively different fluxes in different parts of the single pathway. For the same reason some enzyme modulations given decreased fluxes in one part and increased fluxes in another. The understanding of control in the pathway thus involves consideration of many mechanisms operating simultaneously and the estimation of changes in the whole system. The concept of a 'rate-limiting step' is found to be inadequate and is replaced by a quantitative measure, the Sensitivity Coefficient, which takes account of all the interactions. It is shown that control of the flux is shared among all the enzymes of the pathway. The results are discussed in terms of the theory of flux control.     

Hexokinase and not glycogen synthase controls the flux through the glycogen synthesis pathway in frog oocytes

Here we set out to evaluate the role of hexokinase and glycogen synthase in the control of glycogen synthesis in vivo. We used metabolic control analysis (MCA) to determine the flux control coefficient for each of the enzymes involved in the pathway. Acute microinjection experiments in frog oocytes were specifically designed to change the endogenous activities of the enzymes, either by directly injecting increasing amounts of a given enzyme (HK, PGM and UGPase) or by microinjection of a positive allosteric effector (glc-6P for GS). Values of 0.61 ± 0.07, 0.19 ± 0.03, 0.13 ± 0.03, and −0.06 ± 0.08 were obtained for the flux control coefficients of hexokinase EC 2.7.1.1 (HK), phosphoglucomutase EC 5.4.2.1 (PGM), UDPglucose pyrophosphorylase EC 2.7.7.9 (UGPase) and glycogen synthase EC 2.4.1.11 (GS), respectively. These values satisfy the summation theorem since the sum of the control coefficients for all the enzymes of the pathway is 0.87. The results show that, in frog oocytes, glycogen synthesis through the direct pathway is under the control of hexokinase. Phosphoglucomutase and UDPG-pyrophosphorylase have a modest influence, while the control exerted by glycogen synthase is null.     

In vitro control analysis of an enzyme system: Experimental and analytical developments

In this paper we describe a flow-through system for reconstituting parts of metabolism from purified enzymes. This involves pumping continuously into a reaction chamber, fresh enzymes and reagents so that metabolic reactions occur in the chamber. The waste products leave the chamber via the outflow so that a steady state can be setup. The system we chose consisted of a single enzyme, lactate dehydrogenase. This enzyme was chosen because it consumes NADH in the chamber which could be monitored spectrophotometrically. The aim of the work was to investigate whether a steady state could be achieved in the flow system and whether a metabolic control analysis could be done. We measured two control coefficients, C_LDH and C_pump for the enzyme flux and NADH concentration and confirmed that the summation theorem applied to this system. The advantage of a flow-through system is that the titrations necessary to estimate the control coefficients can be easily and precisely controlled; this means that accurate estimates for the control coefficients can be obtained. In the paper, we discuss some statistical aspects of the data analysis and some possible applications of the technique, including a method to determine the presence of metabolic channelling between two different enzymes.     

Flux control of sulphate assimilation in Arabidopsis thaliana: adenosine 5'-phosphosulphate reductase is more susceptible than ATP sulphurylase to negative control by thiols

The effect of externally applied L-cysteine and glutathione (GSH) on ATP sulphurylase and adenosine 5'-phosphosulphate reductase (APR), two key enzymes of assimilatory sulphate reduction, was examined in Arabidopsis thaliana root cultures. Addition of increasing L-cysteine to the nutrient solution increased internal cysteine, gamma-glutamylcysteine and GSH concentrations, and decreased APR mRNA, protein and extractable activity. An effect on APR could already be detected at 0.2 mm L-cysteine, whereas ATP sulphurylase was significantly affected only at 2 mm L-cysteine. APR mRNA, protein and activity were also decreased by GSH at 0.2 mm and higher concentrations. In the presence of L-buthionine-S, R-sulphoximine (BSO), an inhibitor of GSH synthesis, 0.2 mm L-cysteine had no effect on APR activity, indicating that GSH formed from cysteine was the regulating substance. Simultaneous addition of BSO and 0.5 mm GSH to the culture medium decreased APR mRNA, enzyme protein and activity. ATP sulphurylase activity was not affected by this treatment. Tracer experiments using (35)SO(4)(2-) in the presence of 0.5 mm L-cysteine or GSH showed that both thiols decreased sulphate uptake, APR activity and the flux of label into cysteine, GSH and protein, but had no effect on the activity of all other enzymes of assimilatory sulphate reduction and serine acetyltransferase. These results are consistent with the hypothesis that thiols regulate the flux through sulphate assimilation at the uptake and the APR step. Analysis of radioactive labelling indicates that the flux control coefficient of APR is more than 0.5 for the intracellular pathway of sulphate assimilation. This analysis also shows that the uptake of external sulphate is inhibited by GSH to a greater extent than the flux through the pathway, and that the flux control coefficient of APR for the pathway, including the transport step, is proportionately less, with a significant share of the control exerted by the transport step.     

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.     

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.     

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.