Kontrollmechanismen der Glycolyse

Control mechanisms of enzymic reactions are generally based on interactions between activators, inhibitors, substrates and products with enzyme proteins or on induction and repression of enzyme synthesis. All main types of control can be recognized in glycolysis. They are the basis of the network which controls the over-all glycolytic function and operates according to the feed-back principle. — Enzyme profiles may not be used for a functional definition of the metabolic state. Enzyme activities are governed by a variety of control mechanisms, which can best be recognized by steady-state and transient state analysis of metabolites and by an analysis of the system's response in titration experiments with pure enzymes under conditions whereby the system displays oscillatory behaviour of its over-all flux. The important parameter for the definition of the metabolic state is the net-flux through the system, since this parameter, along with the steady-state levels of the meabolites, gives the steady-state flux pattern, reveals the kinetic state of enzymic reactions and points to control points of metabolism. Continuous glycolytic oscillations in a cell-free extract ofSaccharomyces carlsbergensis have been observed over a period of 22 hours with a constant frequency of 0.17 min^?1 and a rate of 7.2 nMol ethanol per mg protein per min. Titration of such an extract by pure yeast enzymes reveals the gain (FDP, ADP) and damping components (ATP), which are fed back to the enzymes PFK and PK, respectively, PFK, PGK and PK operating as the control units. On the basis of the titration data as well as metabolite and enzyme activity phase relationship, the mechanism of this oscillation can be understood as a crossed feed-back interaction. Furthermore, it is discussed as the biochemical model of a physiological clock mechanism.     

Intensity of glycolysis with brain edema in rats and application of some pharmacological agents

In the in vitro experiments the state of glycolytic processes was studied in the rat brain tissues with development of edema induced by injury and administered nicotine and monoiodoacetate as well as with preliminary introduction of phentolamine, benactyzine, aminazine and trifluoperazine. The brain edema is shown to develop against a background of sharp inhibition of glycolysis independent of the initial content of lactic acid in the incubation medium. Phentolamine, benactyzine and aminazine which prevent from the edema contribute to the normal level of the glycolytic processes. Trifluoperazine ineffective with edema does not prevent from glycolysis inhibition. Inhibition of glycolysis resulting from a decrease in the activity of the corresponding enzymic systems is supposed to be one of the main pathogenetic factors of the brain edema development as it leads to deficit of energy necessary for normal water-electrolytic metabolism in cells.     

A theoretical treatment of damped oscillations in the transient state kinetics of single-enzyme reactions

An extension of the available kinetic theory for reactions in the transient state is presented which establishes that single-enzyme reactions may exhibit damped oscillations under the conditions of standard kinetic experiments performed by stopped-flow techniques. Such oscillations may occur for reasonable magnitudes of rate constants in the enzymic reaction mechanism and at physiological concentrations of enzyme and substrate. In the simplest reaction systems, the oscillations will be strongly damped and lead to progress curves resembling those of a reaction governed by standard exponential transients; statistical regression methods may then have to be applied for their detection and characterization. The observation that single-enzyme reactions may exhibit oscillatory behaviour points to a previously unrecognized possible source of the damped oscillations observed in metabolic systems such as the pathways of glycolysis or photosynthesis.     

Analysis of the oscillatory kinetics of glycolytic intermediates in a yeast extract by FT-IR spectroscopy

In the present work we demonstrate that FT-IR spectroscopy is a powerful tool for the time resolved and noninvasive measurement of multi-substrate/product interactions in complex metabolic networks as exemplified by the oscillating glycolysis in yeast extract. We found that many of the glycolytic intermediates can be identified with FT-IR spectroscopy. For this, we have constructed a spectral library of most of the glycolytic intermediates and obtained the kinetics of single components in spectra from glycolysing yeast extract by the use of mathematical fitting procedures. The results are in good agreement with the known phase relationships of oscillatory glycolysis. They provide the basis for future application of this method to investigate the energy metabolism of living cells.     

Oscillations of the glycolytic pathway and the purine nucleotide cycle.

In a treatment modeled after the oscillatory behavior of the glycolytic pathway and the purine nucleotide cycle observed in skeletal muscle extracts, it is shown that the basis of the oscillations is the AMP-dependent activation of phosphofructokinase by fructose diphosphate. Control of phosphofructokinase by the adenine nucleotides alone leads to the establishment of a steady state. Whether steady state or oscillatory behavior occurs depends in part on the activity of glyceraldehyde-3-phosphate dehydrogenase, which controls the rate of removal of fructose diphosphate. Under appropriate conditions oscillatory behavior can maintain a higher [ATP]/[ADP] ratio than steady state behavior. Viewed in the context of conditions that may be encountered in skeletal muscle in vivo, oscillatory behavior of glycolysis is shown to have additional advantages for maintaining a high [ATP]/[ADP] ratio.     

Glycolysis Is Governed by Growth Regime and Simple Enzyme Regulation in Adherent MDCK Cells

Due to its vital importance in the supply of cellular pathways with energy and precursors, glycolysis has been studied for several decades regarding its capacity and regulation. For a systems-level understanding of the Madin-Darby canine kidney (MDCK) cell metabolism, we couple a segregated cell growth model published earlier with a structured model of glycolysis, which is based on relatively simple kinetics for enzymatic reactions of glycolysis, to explain the pathway dynamics under various cultivation conditions. The structured model takes into account in vitro enzyme activities, and links glycolysis with pentose phosphate pathway and glycogenesis. Using a single parameterization, metabolite pool dynamics during cell cultivation, glucose limitation and glucose pulse experiments can be consistently reproduced by considering the cultivation history of the cells. Growth phase-dependent glucose uptake together with cell-specific volume changes generate high intracellular metabolite pools and flux rates to satisfy the cellular demand during growth. Under glucose limitation, the coordinated control of glycolytic enzymes re-adjusts the glycolytic flux to prevent the depletion of glycolytic intermediates. Finally, the model''s predictive power supports the design of more efficient bioprocesses.     

Kinetics of metabolic pathways. Transient response of the glycolytic system after phosphofructokinase reaction to ADP input

1. Transient response of the glycolytic system, after phosphofructokinase reaction, to ADP input was studied in a metabolic model in vitro, obtained from rat muscle. 2. Under ADP addition to the system, NADH concentration is largely enhanced to be decreased later, showing a very characteristic pulse whose parameters can be used to describe some features of the transition state. 3. NADH pulse was repeatedly obtained by means of successive ADP inputs to the system. 4. The effect of different variables on NADH pulse was studied, these results suggesting a possible role for different steps on this response. 5. The possible relation between the pulse of NADH and oscillatory behaviour of glycolysis is discussed.     

Oscillations and efficiency in glycolysis

We suggest that temporal oscillations of concentrations of intermediates in biochemical reaction systems may enhance the efficiency of free energy conversion (reduce dissipation) in those reactions. Experiments on glycolysis are used to estimate the Gibbs free energy changes along the glycolysis mechanism, and to postulate a construct for the glycolysis "machine" which involves: the PFK reaction as the primary oscillophor; the GAPDH reaction as a phase-shifting device; and the PK reaction with the property of intrinsic oscillatory response at resonance with the driving frequency. Analysis of a simple reaction mechanism with these postulated properties shows that the conversion of free energy from reactants to products is more efficient in an oscillatory than a steady state operation. The efficiency of free energy conversion in glycolysis from glucose + ADP to products + ATP is estimated to be increased by 5--10% due to oscillations. This may have been relevant for the evolutionary development of oscillations such as in glycolysis, especially in anaerobic cells.     

Identification of metabolic fluxes in hepatic cells from transient 13C-labeling experiments: Part II. Flux estimation

This contribution addresses the identification of metabolic fluxes and metabolite concentrations in mammalian cells from transient (13)C-labeling experiments. Whilst part I describes experimental set-up and acquisition of required metabolite and (13)C-labeling data, part II focuses on setting up network models and the estimation of intracellular fluxes. Metabolic fluxes were determined in glycolysis, pentose-phosphate pathway (PPP), and citric acid cycle (TCA) in a hepatoma cell line grown in aerobic batch cultures. In glycolytic and PPP metabolite pools isotopic stationarity was observed within 30 min, whereas in the TCA cycle the labeling redistribution did not reach isotopic steady state even within 180 min. In silico labeling dynamics were in accordance with in vivo (13)C-labeling data. Split ratio between glycolysis and PPP was 57%:43%; intracellular glucose concentration was estimated at 101.6 nmol per 10(6) cells. In contrast to isotopic stationary (13)C-flux analysis, transient (13)C-flux analysis can also be applied to industrially relevant mammalian cell fed-batch and batch cultures.     

Transduction of intracellular and intercellular dynamics in yeast glycolytic oscillations

Under certain well-defined conditions, a population of yeast cells exhibits glycolytic oscillations that synchronize through intercellular acetaldehyde. This implies that the dynamic phenomenon of the oscillation propagates within and between cells. We here develop a method to establish by which route dynamics propagate through a biological reaction network. Application of the method to yeast demonstrates how the oscillations and the synchronization signal can be transduced. That transduction is not so much through the backbone of glycolysis, as via the Gibbs energy and redox coenzyme couples (ATP/ADP, and NADH/NAD), and via both intra- and intercellular acetaldehyde.     

Expression of genes encoding F(1)-ATPase results in uncoupling of glycolysis from biomass production in Lactococcus lactis

We studied how the introduction of an additional ATP-consuming reaction affects the metabolic fluxes in Lactococcus lactis. Genes encoding the hydrolytic part of the F(1) domain of the membrane-bound (F(1)F(0)) H(+)-ATPase were expressed from a range of synthetic constitutive promoters. Expression of the genes encoding F(1)-ATPase was found to decrease the intracellular energy level and resulted in a decrease in the growth rate. The yield of biomass also decreased, which showed that the incorporated F(1)-ATPase activity caused glycolysis to be uncoupled from biomass production. The increase in ATPase activity did not shift metabolism from homolactic to mixed-acid fermentation, which indicated that a low energy state is not the signal for such a change. The effect of uncoupled ATPase activity on the glycolytic flux depended on the growth conditions. The uncoupling stimulated the glycolytic flux threefold in nongrowing cells resuspended in buffer, but in steadily growing cells no increase in flux was observed. The latter result shows that glycolysis occurs close to its maximal capacity and indicates that control of the glycolytic flux under these conditions resides in the glycolytic reactions or in sugar transport.     

Coherent and robust modulation of a metabolic network by cytoskeletal organization and dynamics

In order to investigate the influence of cytoskeletal organization and dynamics on cellular biochemistry, a mathematical model was formulated based on our own experimental evidence. The model couples microtubular protein (MTP) dynamics to the glycolytic pathway and its branches: the Krebs cycle, ethanolic fermentation, and the pentose phosphate (PP) pathway. Results show that the flux through glycolysis coherently and coordinately increases or decreases with increased or decreased levels of polymerized MTP, respectively. The rates of individual enzymatic steps and metabolite concentrations change with the polymeric status of MTP throughout the metabolic network. Negative control is exerted by the PP pathway on the glycolytic flux, and the extent of inhibition depends inversely on the polymerization state of MTP, i.e. a high degree of polymerization relieves the negative control. The stability of the model's steady state dynamics for a wide range of variation of metabolic parameters increased with the degree of polymerized MTP. The findings indicate that the organization of the cytoskeleton bestows coherence and robustness to the coordination of cellular metabolism.     

Mechanistic model of cardiac energy metabolism predicts localization of glycolysis to cytosolic subdomain during ischemia

A new multidomain mathematical model of cardiac cellular metabolism was developed to simulate metabolic responses to reduced myocardial blood flow. The model is based on mass balances and reaction kinetics that describe transport and metabolic processes of 31 key chemical species in cardiac tissue. The model has three distinct domains (blood, cytosol, and mitochondria) with interdomain transport of chemical species. In addition to distinguishing between cytosol and mitochondria, the model includes a subdomain in the cytosol to account for glycolytic metabolic channeling. Myocardial ischemia was induced by a 60% reduction in coronary blood flow, and model simulations were compared with experimental data from anesthetized pigs. Simulations with a previous model without compartmentation showed a slow activation of glycogen breakdown and delayed lactate production compared with experimental results. The addition of a subdomain for glycolysis resulted in simulations showing faster rates of glycogen breakdown and lactate production that closely matched in vivo experimental data. The dynamics of redox (NADH/NAD+) and phosphorylation (ADP/ATP) states were also simulated. These controllers are coupled to energy transfer reactions and play key regulatory roles in the cytosol and mitochondria. Simulations showed a similar dynamic response of the mitochondrial redox state and the rate of pyruvate oxidation during ischemia. In contrast, the cytosolic redox state displayed a time response similar to that of lactate production. In conclusion, this novel mechanistic model effectively predicted the rapid activation of glycogen breakdown and lactate production at the onset of ischemia and supports the concept of localization of glycolysis to a subdomain of the cytosol.     

Identification of flux control in metabolic networks using non-equilibrium thermodynamics

A method is presented to identify flux controlling reactions in metabolic networks using experimentally determined flux distributions. The method is based on the application of Ziegler's principle for the maximization of entropy production. According to this principle a metabolic network tends to maximize the entropy production rate while satisfying mass balances and maximal rate constraints. Experimental flux data corresponding to four different metabolic states of Saccharomyces cerevisiae were used to identify the corresponding flux controlling reactions. The bottleneck nature of several of the identified reactions was confirmed by earlier studies on over-expression of the identified target genes. The method also explains the failure of all the previous trials of increasing the glycolysis rate by direct over-expression of several glycolytic enzymes. These findings point to a wider use of the method for identification of novel targets for metabolic engineering of microorganisms used for sustainable production of fuels and chemicals.     

Design of glycolysis

The design of the glycolytic pathway resulting from the continuous refinement of evolution is discussed with regard to three aspects. 1. Functional and structural properties of individual enzymes. The catalytic constants of the glycolytic enzymes are remarkably optimized; the turnover numbers are within one order of magnitude. The same is true for the molarities of catalytic centres in the cytosol, as is noted for yeast. Functional properties of the enzymes are reflected in their tertiary and quaternary structures. 2. Regulatory mechanisms of single enzymes. A classification of the various types of enzymic control mechanisms operating in the glycolytic pathway is given. In addition to the usual Michaelis-Menten saturation kinetics and the various types of inhibition there is control by positive and negative effectors based on oligomeric structures (fast acting, fine control) as well as regulation by chemical interconversion structures (fast acting, fine control) as well as regulation by chemical based on enzymes cascades (slow acting, very effective). 3. Functional and regulatory mechanisms of the whole glycolytic reaction pathway. A prominent feature is the high enzyme:substrate ratio, which guarantees fast response times. However, a quantitative treatment of the overall kinetics is limited by an incomplete knowledge of the enzymes' dynamic and chemical compartmentation as well as some of their control properties. From an analysis of the oscillatory state, certain control points in the glycolytic chain can be located that coincide with major branching points to other metabolic pathways. These points are controlled by fast-acting cooperative enzymes that operate in a flip-flop mechanism together with the respective antagonistic enzymes, preventing futile cycles. The gating enzymes leading to the glycogen store and the citric acid cycle are of the slow-acting but very effective interconvertible type. The combination of all the complex and intricate features of design yields a glycolytic network that enables the cell to respond to its various metabolic needs quickly, effectively and economically.     

Expression of Genes Encoding F1-ATPase Results in Uncoupling of Glycolysis from Biomass Production in Lactococcus lactis

We studied how the introduction of an additional ATP-consuming reaction affects the metabolic fluxes in Lactococcus lactis. Genes encoding the hydrolytic part of the F₁ domain of the membrane-bound (F₁F₀) H⁺-ATPase were expressed from a range of synthetic constitutive promoters. Expression of the genes encoding F₁-ATPase was found to decrease the intracellular energy level and resulted in a decrease in the growth rate. The yield of biomass also decreased, which showed that the incorporated F₁-ATPase activity caused glycolysis to be uncoupled from biomass production. The increase in ATPase activity did not shift metabolism from homolactic to mixed-acid fermentation, which indicated that a low energy state is not the signal for such a change. The effect of uncoupled ATPase activity on the glycolytic flux depended on the growth conditions. The uncoupling stimulated the glycolytic flux threefold in nongrowing cells resuspended in buffer, but in steadily growing cells no increase in flux was observed. The latter result shows that glycolysis occurs close to its maximal capacity and indicates that control of the glycolytic flux under these conditions resides in the glycolytic reactions or in sugar transport.     

The fundamental organization of cardiac mitochondria as a network of coupled oscillators

Mitochondria can behave as individual oscillators whose dynamics may obey collective, network properties. We have shown that cardiomyocytes exhibit high-amplitude, self-sustained, and synchronous oscillations of bioenergetic parameters when the mitochondrial network is stressed to a critical state. Computational studies suggested that additional low-amplitude, high-frequency oscillations were also possible. Herein, employing power spectral analysis, we show that the temporal behavior of mitochondrial membrane potential (DeltaPsi(m)) in cardiomyocytes under physiological conditions is oscillatory and characterized by a broad frequency distribution that obeys a homogeneous power law (1/f(beta)) with a spectral exponent, beta = 1.74. Additionally, relative dispersional analysis shows that mitochondrial oscillatory dynamics exhibits long-term memory, characterized by an inverse power law that scales with a fractal dimension (D(f)) of 1.008, distinct from random behavior (D(f) = 1.5), over at least three orders of magnitude. Analysis of a computational model of the mitochondrial oscillator suggests that the mechanistic origin of the power law behavior is based on the inverse dependence of amplitude versus frequency of oscillation related to the balance between reactive oxygen species production and scavenging. The results demonstrate that cardiac mitochondria behave as a network of coupled oscillators under both physiological and pathophysiological conditions.     

Simple and complex spatiotemporal structures in a glycolytic allosteric enzyme model

Pattern formation in glycolysis is studied with a classical reaction-diffusion allosteric enzyme model. It is found that, similar to recent experimental reports in the yeast extracts, a small magnitude local perturbation can induce transient target waves in a two dimensional oscillatory medium. An above threshold stimulation generates target waves which eventually evolve into spatiotemporal chaos upon collisions with the boundary or other wave activities. Detailed simulation studies show that the studied simple glycolytic reaction-diffusion model can support three types of spatiotemporal behaviors which are independent of the boundary conditions: (1) a spatially uniform stable steady state, (2) periodic global oscillations and (3) spatiotemporal chaos.