Human red cell glycolysis in high altitude chronic hypoxia

We have found that glycolysis in human red blood cells under the hypoxic conditions found at high altitudes is connected with changes in enzyme activities and levels of various metabolic intermediates. The sensitivity of the four kinases to hypoxia results in 1) glycolytic hyperactivity leading to a higher intracellular energy state, and 2) accumulation of 2–3 DPG, whose role in the adaptation of red blood cell respiration to high altitude has been shown by previous research. PEP, 3PG, and G6P appear to be the main regulating intermediates in glycolysis in this system. The reason for the very large increase in G1-6DP is still not clear.     

Experimental determination of flux control distribution in biochemical systems: in vitro model to analyze transient metabolite concentrations

Determination of the control coefficients allows the identification of rate-controlling steps in a reaction system. However, the measurement of the flux control coefficients in a biochemical system is not a trivial task, except for some special cases. We have developed a theoretical basis for the direct determination of these coefficients from dynamic responses. In order to show the validity of this methodology experimentally, the dynamic approach is applied to an in vitro reconstituted partial glycolytic pathway to determine the flux control coefficients of hexokinase and phosphofructokinase. It is shown that the dynamic approach gives consistent results, which agree well with values obtained by the direct enzyme titration method. The detailed procedure and potential applications to other systems, such as immobilized enzyme or cell reactors, are discussed. (c) 1993 Wiley & Sons, Inc.     

Model reduction and a priori kinetic parameter identifiability analysis using metabolome time series for metabolic reaction networks with linlog kinetics

In this work, we present a time-scale analysis based model reduction and parameter identifiability analysis method for metabolic reaction networks. The method uses the information obtained from short term chemostat perturbation experiments. We approximate the time constant of each metabolite pool by their turn-over time and classify the pools accordingly into two groups: fast and slow pools. We performed a priori model reduction, neglecting the dynamic term of the fast pools. By making use of the linlog approximative kinetics, we obtained a general explicit solution for the fast pools in terms of the slow pools by elaborating the degenerate algebraic system resulting from model reduction. The obtained relations yielded also analytical relations between a subset of kinetic parameters. These relations also allow to realize an analytical model reduction using lumped reaction kinetics. After solving these theoretical identifiability problems and performing model reduction, we carried out a Monte Carlo approach to study the practical identifiability problems. We illustrated the methodology on model reduction and theoretical/practical identifiability analysis on an example system representing the glycolysis in Saccharomyces cerevisiae cells.     

Aerobic Glycolysis by Proliferating Cells: Protection against Oxidative Stress at the Expense of Energy Yield

Primary cultures of mitogen-activated rat thymocytes were used to study energy metabolism, gene expression of glycolytic enzymes, and production of reactive oxygen species during cell cycle progression. During transition from the resting to the proliferating state a 7- to 10-fold increase of glycolytic enzyme induction occurs which enables the cells to meet the enhanced energy demand by increased aerobic glycolysis. Cellular redox changes have been found to regulate gene expression of glycolytic enzymes by reversible oxidative inactivation of Spl-binding to the cognate DNA-binding sites in the promoter region. In contrast to nonproliferating cells, production of phorbol 12-myristate 13-acetate (PMA)-primed reactive oxygen species (ROS) in proliferating rat thymocytes and HL-60 cells is nearly abolished. Pyruvate, a product of aerobic glycolysis, is an effective scavenger of ROS, which could be shown to be generated mainly at the site of complex III of the mitochondrial respiratory chain. Aerobic glycolysis by proliferating cells is discussed as a means to minimize oxidative stress during the phases of the cell cycle when maximally enhanced biosynthesis and cell division do occur.     

Pyruvate kinase-catalyzed ATP-formation in human red blood cell membranes

Previous studies on the linkage between enzymatically catalyzed ATP-generating reactions in the red blood cell membrane and the sodium and potassium transport in the control of overall glycolysis of human erythrocytes were controversial. In this study a significant amount of pyruvate kinase activity is shown to be localized within the membrane. Membrane fragments produce 20.5 mumol of ATP per 10(10) membranes per hour from phosphoenolpyruvate and ADP. The kinetics of the membrane-localized pyruvate kinase do not differ from those of the enzyme from hemolysates. The results clearly document the presence of the second ATP-generating enzyme of glycolysis, pyruvate kinase, in human red blood cell membranes. The main fraction of the enzyme is deeply hidden in the lipid layers of the membrane. It can be demasked by mechanical desintegration of membranes at high levels of activity. It is suggested that the amount of the membrane-localized fraction of pyruvate kinase is related to the clinical severity of the hemolytic process in pyruvate kinase deficiency.     

A computer model of pancreatic islet glycolysis

We have modeled an experiment with perifused pancreatic islet cells using our BIOSSIM language. The experiment and the resulting model are concerned with glucose uptake and glycolysis by the beta-cells of pancreatic islets. Although glycolysis appears to be involved in insulin release, we do not have enough information to represent insulin release in detail. The rapid entry of glucose into the beta-cell is promoted by a carrier having a very high tissue capacity. Phosphorylation of glucose by the low affinity enzyme glucokinase appears to be limiting for glycolysis. The effects of several hexose diphosphate activators of phosphofructokinase are modeled. Model behavior is described. The kinetic parameters of the enzyme submodels are given. Because of the difficulties of preparing large amounts of experimental material, information on pancreatic islet metabolism is limited. This model is a plausible explanation of the experimental results. Recent work on the genetically engineered glucose transporter and glucokinase is discussed.     

The activities of phosphorylase, hexokinase, phosphofructokinase, lactate dehydrogenase and the glycerol 3-phosphate dehydrogenase in muscles from vertebrates and invertebrates

1. The maximum activities of hexokinase, phosphorylase and phosphofructokinase have been measured in extracts from a variety of muscles and they have been used to estimate the maximum rates of operation of glycolysis in muscle. These estimated rates of glycolysis are compared with those calculated for the intact muscle from such information as oxygen uptake, glycogen degradation and lactate formation. Reasonable agreement between these determinations is observed, and this suggests that such enzyme activity measurements may provide a useful method for comparative investigations into quantitative aspects of maximum glycolytic flux in muscle. 2. The enzyme activities from insect flight muscle confirm and extend much of the earlier work and indicate the type of fuel that can support insect flight. The maximum activity of hexokinase in some insect flight muscles is about tenfold higher than that in vertebrate muscles. The activity of phosphorylase is greater, in general, in vertebrate muscle (particularly white muscle) than in insect flight muscle. This is probably related to the role of glycogen breakdown in vertebrate muscle (particularly white muscle) for the provision of ATP from anaerobic glycolysis and not from complete oxidation of the glucose residues. The activity of hexokinase was found to be higher in red than in white vertebrate muscle, thus confirming and extending earlier reports. 3. The maximum activity of the mitochondrial glycerophosphate dehydrogenase was always much lower than that of the cytoplasmic enzyme, indicating that the former enzyme is rate-limiting for the glycerol 3-phosphate cycle. From the maximum activity of the mitochondrial enzyme it can be calculated that the operation of this cycle would account for the reoxidation of all the glycolytically produced NADH in insect flight muscle but it could account for only a small amount in vertebrate muscle. Other mechanisms for this NADH reoxidation in vertebrate muscle are discussed briefly.     

Relations among variations in human red cell volume,density, membrane area,hemoglobin content and cation content

It is shown how variations in different properties of red cells can be inter-related provided relations exist among these properties at the single cell level. On the basis of the cell density dependence on cell volume and hemoglobin content, and the assumed volume dependence on red cell cation and hemoglobin content, nine relations among the variations in red cell volume, density, membrane area, hemoglobin content and cation content, and their correlations are derived. Values of seven correlation coefficients are theoretically predicted and are shown to be consistent with the experiments performed by density fractionated red blood cells. The cell volume dependence on cation and hemoglobin content obtained from relations among variations is compared with the predictions obtained by the existing model about the osmotic behavior of the red blood cell. Furthermore, it is shown that data on the variations of the red cell properties indicate the existence of the relation among cation content, hemoglobin content, and membrane area at the level of a single cell.     

Optimal stoichiometric designs of ATP-producing systems as determined by an evolutionary algorithm.

The design of metabolic pathways is thought to be the result of an optimization process such that the structure of contemporary metabolic routes maximizes a particular objective function. Recently, it has been shown that some essential stoichiometric properties of glycolysis can be explained on the basis of the requirement for a high ATP production rate. Because the number of stoichiometrically feasible designs increases strongly with the number of reactions involved, a systematic analysis of all the possibilities turns out to be inaccessible beyond a certain system size. We present, therefore, an alternative approach to compute in a more efficient way the optimal design of glycolysis interacting with an external ATP-consuming reaction. The algorithm is based on the laws of evolution by natural selection, and may be viewed as a particular version of evolutionary algorithms. The following conclusions are derived: (a) evolutionary algorithms are very useful search strategies in determining optimal stoichiometries of metabolic pathways. (b) Essential topological features of the glycolytic network may be explained on the basis of flux optimization. (c) There is a strong interrelation between the optimal stoichiometries and the thermodynamic and kinetic properties of the participating reactions. (d) Some subsequences of reactions in optimal pathways are strongly conserved at variation of system parameters, which may be understood by applying principles of metabolic control analysis.     

Nonequilibrium self-assembly of linear fibers: microscopic treatment of growth, decay, catastrophe and rescue

Many of the large structures of cells are constructed from fibers. These fibers self-assemble from individual proteins in a far-from-equilibrium fashion. Nonequilibrium self-assembly results in a highly dynamic process at the subcellular level that can be regulated and tuned to carry out many of the biological functions of the cell: growth, division and locomotion. We construct and analyze a nonequilibrium model of the dynamic end of a biological fiber that possesses site-resolved resolution. We solve for the steady states of this nonequilibrium system using a variational method. The results are compared to exact numerical solutions for systems with modest size. Using an effective reaction coordinate, we construct an effective potential from the steady-state distribution. The stochastic transitions of the system can be analyzed in this representation. We then apply this method to model microtubule systems. Predictions for macroscopic catastrophe, rescue and dynamic instability in the steady states are analyzed. We find that the length of the cap of the microtubule is small. The relations between the catastrophe/rescue rate and the growth rate are also discussed.     

Carbohydrate metabolism in transforming lymphocytes from the aged

There is an age-related decline in immune capacity which has been linked to a decreased response of lymphocytes to mitogens in vitro. During transformation, lymphocytes require a marked increase in energy production and biosynthesis which is supplied primarily by glycolysis. In the elderly, the glycolytic enzymes increase significantly in transforming lymphocytes at least 24 hr later than in the young and then at significantly reduced levels. Glucose utilization is also impaired in stimulated lymphocytes from the elderly but follows the impairment of glycolysis. In stimulated cells from the young, increases in glycolytic enzyme activity levels accompany sharp increases in blastogenesis while a delayed increase in glycolytic enzyme activity in the elderly is accompanied by a delay in blastogenesis. Maximal glycolytic enzyme activity levels are significantly reduced in transformed lymphocytes from the elderly though the number of transformed cells is also significantly reduced. However, glycolytic enzyme activity levels are significantly lower in the elderly than in the young even on a per transformed cell basis. Thus, this reduction cannot be attributed to the lower number of transformed cells that are present in the elderly. This defect in the increase of glycolysis in stimulated cells from the elderly suggests an intracellular mechanism which could be related to the impaired lymphocyte stimulation in vitro in the aged.     

Factors influencing Na+ transport in dog red cells

Na⁺ movements in dog red cells have been measured in a study of the relationship between cell volume, Na⁺ permeability and glycolysis. When dog red cells are shrunken by 20% at 38 °C the apparent Na⁺ influx increases by a factor of about fifty, and the effect remains when cells are deprived of glucose for –2.5 h. Flux returns to normal when the cells are restored to their initial volume. Glycolysis is required for the volume effect and we have studied the effect of glycolytic modifiers such as fluoride, sulfate, bisulfite and pyruvate on these glucose depleted dog red cells. The results indicate that the volume effect is associated with a change in the concentration of 3-phosphoglycerate and may be mediated by phosphoglycerate kinase, the membrame-associated enzyme which forms 3-phosphoglycerate from 1,3-diphosphoglycerate. The state of high Na⁺ permeability persists for several hours in the absence of glucose and it appears that shrinking the cells has opened a Na⁺-specific channel through which this cation can exchange easily.     

Metabolic pathway analysis of enzyme-deficient human red blood cells

Five enzymopathies (G6PDH, TPI, PGI, DPGM and PGK deficiencies) in the human red blood cells are investigated using a stoichiometric modeling approach, i.e., metabolic pathway analysis. Elementary flux modes (EFMs) corresponding to each enzyme deficiency case are analyzed in terms of functional capabilities. When available, experimental findings reported in literature related to metabolic behavior of the human red blood cells are compared with the results of EFM analysis. Control-effective flux (CEF) calculation, a novel approach which allows quantification and interpretation of determined EFMs, is performed for further analysis of enzymopathies. Glutathione reductase reaction is found to be the most effective reaction in terms of its CEF value in all enzymopathies in parallel with its known essential role for red blood cells. Efficiency profiles of the enzymatic reactions upon the degree of enzyme deficiency are obtained by the help of the CEF approach, as a basis for future experimental studies. CEF analysis, which is found to be promising in the analysis of erythrocyte enzymopathies, has the potential to be used in modeling efforts of human metabolism.     

Glycolysis of red cells suspended in solutions of impermeable solutes. Intracellular pH and glycolysis.

The glycolytic rate human red cells suspended in a sucrose medium of low or physiological pH was higher than that of the cells suspended in Ringer's medium of the same. pH. The medium pHP-glycolytic rate curve of red cells suspended in soucrose media shifted to the acidic side by about one unit compared with that of cells suspended in Ringer's medium. Similarly, the pattern of glycolytic intermediates in red cells suspended in a sucrose medium resembled that in cells suspended in Ringer's solution of about one unit higher pH. These phenomena could be ascribed to the change of intracellular pH, which was measured by the 5,5'-dimethyl-oxazolidine-2,4-dione method. A similar stimulation of glycolysis was observed when sodium citrate was added to red cells suspended in Ringer's solution at constant pH. These observations indicate that membrane-impermeable non-electrolytes or anions stimulate glycolysis of red cells by elevation ofthe intracellular pH. Red cell glycolysis is influenced mainly by the intracellular pH rather than by the pH of the suspending medium.     

Cellular energetics as a target for tumor cell elimination

Investigation of cancer cell metabolism has revealed variability of the metabolic profiles among different types of tumors. According to the most classical model of cancer bioenergetics, malignant cells primarily use glycolysis as the major metabolic pathway and produce large quantities of lactate with suppressed oxidative phosphorylation even in the presence of ample oxygen. This is referred to as aerobic glycolysis, or the Warburg effect. However, a growing number of recent studies provide evidence that not all cancer cells depend on glycolysis, and, moreover, oxidative phosphorylation is essential for tumorigenesis. Thus, it is necessary to consider distinctive patterns of cancer metabolism in each specific case. Chemoresistance of cancer cells is associated with decreased sensitivity to different types of antitumor agents. Stimulation of apoptosis is a major strategy for elimination of cancer cells, and therefore activation of mitochondrial functions with direct impact on mitochondria to destabilize them appears to be an important approach to the induction of cell death. Consequently, the design of combination therapies using acclaimed cytotoxic agents directed to induction of apoptosis and metabolic agents affecting cancer cell bioenergetics are prospective strategies for antineoplastic therapy.     

A novel tetrazolium method for peroxidase histochemistry and immunohistochemistry.

We developed a new method for the histochemical demonstration of peroxidase. This method, which has a novel reaction mechanism, is based on the oxidation of phenol by peroxidase and coupling of this reaction to the reduction of a tetrazolium salt, with the deposition of an insoluble formazan at sites of enzyme activity. This new method was compared with an established diaminobenzidine (DAB) technique for peroxidase histochemistry and immunohistochemistry. Although both methods identified peroxidase activity in myeloid cells of bone marrow biopsy specimens, there was no interference from red cell pseudoperoxidase activity with the phenol-tetrazolium method, in contrast to the diaminobenzidine method. The detection of cytokeratin using an indirect immunoperoxidase technique was compared with both methods for demonstrating peroxidase activity. The phenol-tetrazolium method gave results similar to that obtained with DAB and appeared to be at least as sensitive as DAB in detecting low amounts of antigen. In addition, the production of a formazan as the final reaction product means that the phenol-tetrazolium method is ideally suited for quantitative peroxidase histochemistry. Therefore, the phenol-tetrazolium method represents a useful alternative method to DAB and for certain applications offers significant advantages over DAB.     

Simplified modelling of metabolic pathways for flux prediction and optimization: lessons from an in vitro reconstruction of the upper part of glycolysis

Explicit modelling of metabolic networks relies on well-known mathematical tools and specialized computer programs. However, identifying and estimating the values of the very numerous enzyme parameters inherent to the models remain a tedious and difficult task, and the rate equations of the reactions are usually not known in sufficient detail. A way to circumvent this problem is to use 'non-mechanistic' models, which may account for the behaviour of the systems with a limited number of parameters. Working on the first part of glycolysis reconstituted in vitro, we showed how to derive, from titration experiments, values of effective enzyme activity parameters that do not include explicitly any of the classical kinetic constants. With a maximum of only two parameters per enzyme, this approach produced very good estimates for the flux values, and enabled us to determine the optimization conditions of the system, i.e. to calculate the set of enzyme concentrations that maximizes the flux. This fast and easy method should be valuable in the context of integrative biology or for metabolic engineering, where the challenge is to deal with the dramatic increase in the number of parameters when the systems become complex.     

Rheological behaviors of bovine blood forming artificial rouleaux

A purpose of the present study is to make an artificial rouleau of bovine red blood cells which is not capable of rouleau formation under physiological condition. Rheological behaviors of bovine blood forming artificial rouleaux were examined. The modification of cell surface by enzyme trypsin induced rouleau formation, whereas the modification of cell surface by neuraminidase did not cause any aggregate formation. The drastic elevation of the fibrinogen content in bovine red blood cells suspension also brought about the formation of rouleau. The value of dynamic rigidity modulus G' of bovine red blood cells in saline solution containing high concentration of fibrinogen is somewhat smaller than that of trypsin treated bovine red blood cells in plasma. The value of G' of trypsin treated bovine red blood cells in plasma first increased to a maximum value and then decreased with the time. It is supposed that the removal of macro-molecules from the cell surface facilitates the mutual approach of cells and causes the formation of rouleau which seems to be the same as that of human and horse bloods.