Physiological, biochemical, and mathematical studies of micro-aerobic continuous ethanol fermentation by Saccharomyces cerevisiae. II: intracellular metabolite and enzyme assays at steady state chemostat cultures

Intracellular metabolite concentration and enzyme activity measurements were made to explain the new metabolic and growth phenomena seen in the micro-aerobic, continuous yeast cultures described in Part I. The results of these assays suggested mechanisms for the observed maximum in the specific ethanol productivity as a function of the oxygen feed rate, changing ATP yields, the effects of antifoam, and the sharp changes in the biomass concentration with small changes in the oxygenation. Measured were the intracellular concentrations of ATP, NADH, glucose 6-phosphate, pyruvate, glycerol, and ethanol, and the activities of hexokinase and alcohol dehydrogenase. Rate-limiting steps were identified by the accumulation of metabolites upstream and the depletion of metabolites downstream of the step.A potential mechanism for the stimulation of fermentation with decreasing oxygenation was an activation of glucose transport by an accumulating intracellular ATP concentration. The inhibition of fermentation at yet lower oxygenation rates may have been caused by the continued accumulation of ATP to the point that the glycolytic kineses were inhibited. A mechanism for the changing ATP yields and intracellular ATP concentration proposed the existence of ATPases or ATP waste reactions stimulated by both oxygen and ATP. Antifoam had the effect of decreasing the resistance for glycerol transport out of the cell. The resulting stimulation of glycerol production and inhibition of ethanol production decreased the intracellular ATP content. Finally, intracellular ethanol was found not to accumulate to levels of higher than the extracellular concentration.     

Physiological, biochemical, and mathematical studies of micro-aerobic continuous ethanol fermentation by Saccharomyces cerevisiae. I: hysteresis, oscillations, and maximum specific ethanol productivities in chemostat culture

The growth and metabolism of Saccharomyces cerevisiae was studied in steady-state chemostat cultures under conditions of scarce oxygen and excess glucose. The specific ethanol productivity and specific glucose uptake rate were stimulated by 50% within a narrow range of air/nitrogen mixtures to the fermentor. Fermentation was inhibited at slightly higher and lower air/nitrogen ratios, confirming similar results by previous investigators. This stimulation could not be caused by obvious mechanisms, such as the Pasteur or Crabtree effects. Since this maximum in the fermentation rate occurred in a steady-state chemostat and at a constant dilution rate, the ATP yield of the culture necessarily attained a minimum. Thus, changes in the energetic efficiency of growth or the degree of wasting of ATP were surmised. The steady-state biomass concentration at various oxygenation rates exhibited hysteresis phenomena. Ignition and extinction of the biomass concentration occurred as critical oxygen feed rates were passed. The hysteresis was prevented by adding yeast extract to or reducing the antifoam concentration in the medium. These medium alterations had the simultaneous effect of stimulating the fermentation rate, suggesting that ATP has a critical role in dictating the biomass concentration in micro-aerobic culture. Silicone polymer antifoam was found to stimulate glycerol production at the expense of ethanol production, having consequences for the energy generation and the biomass concentration of the culture.     

Glycolytic flux is conditionally correlated with ATP concentration in Saccharomyces cerevisiae: a chemostat study under carbon- or nitrogen-limiting conditions.

Anaerobic and aerobic chemostat cultures of Saccharomyces cerevisiae were performed at a constant dilution rate of 0.10 h(-1). The glucose concentration was kept constant, whereas the nitrogen concentration was gradually decreasing; i.e., the conditions were changed from glucose and energy limitation to nitrogen limitation and energy excess. This experimental setup enabled the glycolytic rate to be separated from the growth rate. There was an extensive uncoupling between anabolic energy requirements and catabolic energy production when the energy source was present in excess both aerobically and anaerobically. To increase the catabolic activity even further, experiments were carried out in the presence of 5 mM acetic acid or benzoic acid. However, there was almost no effect with acetate addition, whereas both respiratory (aerobically) and fermentative activities were elevated in the presence of benzoic acid. There was a strong negative correlation between glycolytic flux and intracellular ATP content; i.e., the higher the ATP content, the lower the rate of glycolysis. No correlation could be found with the other nucleotides tested (ADP, GTP, and UTP) or with the ATP/ADP ratio. Furthermore, a higher rate of glycolysis was not accompanied by an increasing level of glycolytic enzymes. On the contrary, the glycolytic enzymes decreased with increasing flux. The most pronounced reduction was obtained for HXK2 and ENO1. There was also a correlation between the extent of carbohydrate accumulation and glycolytic flux. A high accumulation was obtained at low glycolytic rates under glucose limitation, whereas nitrogen limitation during conditions of excess carbon and energy resulted in more or less complete depletion of intracellular storage carbohydrates irrespective of anaerobic or aerobic conditions. However, there was one difference in that glycogen dominated anaerobically whereas under aerobic conditions, trehalose was the major carbohydrate accumulated. Possible mechanisms which may explain the strong correlation between glycolytic flux, storage carbohydrate accumulation, and ATP concentrations are discussed.     

A mathematical model of the nutrient dynamics of phytoplankton in a nitrate-limited environment

Insofar as saturation kinetics are applicable to the growth of phytoplankton in laboratory experiments and to growth in nature, the computer modeling of intracellular nutrient partitioning in populations of cells can lead to better understanding of the dynamics of natural populations. A three-compartment mathematical model was developed to represent a phytoplankton population having the capability to store nitrogen in a nitrate-limited environment. Parameters were estimated by fitting the model to data from two chemostat experiments reported by Caperon (1968). The model was used to simulate growth dynamics observed in chemostat and batch experiments. The model demonstrated the changes which may occur in the nitrogenous constituents of a phytoplankton population with time and environmental conditions. The model also demonstrates three phenomena which have been observed in field and laboratory experiments but which are not represented by the customary Monod model: (1) uptake rates may significantly exceed not growth rates, (2) high growth rates may be encountered at very low environmental nitrate concentrations, and (3) the ratio of internal nitrogen to population size may change significantly during a study period. It is suggested that the amount of nitorgen in storage may be used as an indicator of the physiological state of a monospecific population. Parameters for the one-compartment Monod model were estimated by customary methods form data generated by the three-compartment model. It was shown that difficulties encountered in estimating the yield coefficient and the decay coefficient may be attributed to the intracellular storage phenomenon. It was also demonstrated that the one-compartment Monod model was inadequate to accurately represent population growth in chemostat experiments when intracellular storage is a significant factor.     

On the control mechanisms of the nitrite level in Escherichia coli cells: the mathematical model

Background

Due to a high toxicity of nitrite and its metabolites, it is of high interest to study mechanisms underlying the low NO₂ level maintenance in the cell. During anaerobic growth of Escherichia coli the main nitrite-reducing enzymes are NrfA and NirB nitrite reductases. NrfA reductase is localized in the cell periplasm and uses NO₂ as an electron acceptor to create a proton gradient; NirB reductase is restricted to the cytoplasm and metabolizes excessive nitrite inside the cell, the uptake of which is mediated by the transporter protein NirC. While it is known that these three systems, periplasmic, cytoplasmic and transport, determine nitrite uptake and assimilation in the cell as well as its excretion, little is known about their co-ordination.

Results

Using a mathematical model describing the nitrite utilization in E. coli cells cultured in a flow chemostat, the role of enzymes involved in nitrite metabolism and transport in controlling nitrite intracellular levels was investigated. It was demonstrated that the model adapted to the experimental data on expression of nrfA and nirB genes encoding NrfA and NirB nitrite reductases, can describe nitrite accumulation kinetics in the chemostat in the millimolar range of added substrate concentrations without any additional assumptions. According to the model, in this range, low intracellular nitrite level, weakly dependent on its concentration in the growth media, is maintained (mcM). It is not sufficient to consider molecular-genetic mechanisms of NrfA reductase activity regulation to describe the nitrite accumulation dynamics in the chemostat in the micromolar range (≤1 mM) of added nitrite concentrations. Analysis of different hypotheses has shown that the mechanism of local enzyme concentration change due to membrane potential-induced diffusion from the cytoplasm to the periplasm at low nitrite levels is sufficient to explain the nitrite accumulation dynamics in the chemostat.

Conclusions

At nitrite concentrations in the media more than 2 mM, the model adapted to the experimental data on nitrite utilization dynamics in E. coli cells cultured in the flow chemostat demonstrates the largest contribution of genetic mechanisms involved in nrf and nir operons activity regulation to the control of nitrite intracellular levels. The model predicts a significant contribution of the membrane potential to the periplasmic NrfA nitrite reductase activity regulation and nitrite utilization dynamics at substrate concentrations ≤1 mM.

     

Dynamics of energy charge and adenine nucleotides during uncoupling of catabolism and anabolism in Penicillium ochrochloron

Filamentous fungi are able to spill energy when exposed to energy excess by uncoupling catabolism from anabolism, e.g. via overflow metabolism. In current study we tested the hypothesis that overflow metabolism is regulated via the energetic status of the hyphae (i.e. energy charge, ATP concentration). This hypothesis was studied in Penicillium ochrochloron during the steady state of glucose- or ammonium-limited chemostat cultures as well as during three transient states ((i) glucose pulse to a glucose-limited chemostat, (ii) shift from glucose-limited to ammonium-limited conditions in a chemostat, and (iii) ammonium exhaustion in batch culture). Organic acids were excreted under all conditions, even during exponential growth in batch culture as well as under glucose-limited conditions in a chemostat. Partial uncoupling of catabolism and anabolism via overflow metabolism was thus constitutively present. Under all tested conditions, overflow metabolism was independent of the energy charge or the ATP concentration of the hyphae. There was a reciprocal correlation between glucose uptake rate and intracellular adenine nucleotide content. During all transients states a rapid decrease in energy charge and the concentrations of nucleotides was observed shortly after a change in glycolytic flux (“ATP paradoxon”). A possible connection between the change in adenine nucleotide concentrations and the purine salvage pathway is discussed.     

Effects of glucose supply on myeloma growth and metabolism in chemostat culture

The effects of the glucose supply on growth and metabolism of an SP2/0 derived recombinant myeloma cell line were studied in chemostat culture during growth on IMDM medium at a fixed dilution rate of 0.032 h^?1. Lowering of the feed medium glucose concentration from 25.0 to 1.4 mmol/L resulted in a decrease of steady-state viable cell concentration from 1.9 × 10⁹ L^?1, whereas viability remained above 90%. Mass balances indicated that only a minor amount of glucose was utilized via the TCA cycle irrespective of the glucose concentration in the feed medium. The apparent biosynthetic yield of cells from ATP was independent of the ratio between the specific glucose and glutamine consumption rate. It is concluded that the primary role of glucose is the provision of intermediates for anabolic reactions. In addition, glucose may play an indirect catabolic role in the process of glutaminolysis by providing the pyruvate for the transamination of glutamate to alanine and α-ketoglutarate. At low glucose concentrations in the feed medium, glutamine is probably the sole energy source for this myeloma in chemostat culture. © 1995 Wiley-Liss, Inc.     

Maintenance requirements: energy supply from simultaneous endogenous respiration and substrate consumption

Abstract A model is presented that describes energy for maintenance purposes (ATP) as being obtained simultaneously from biomass degradation as well as from substrate degradation in excess of growth requirements. The ratio between both catabolic processes was taken to be growth rate dependent. As such, this approach is intermediate between established models; its significant features are negative growth and the absence of substrate consumption at zero substrate concentration, and the attainability of the maximum specific growth rate (the model parameter μ _max) at elevated substrate concentrations. As a simple case, the amounts of ATP obtained from direct substrate catabolism or from the degradation of an equivalent amount of biomass were taken as identical. Also, the maintenance demand in terms of ATP per unit time and biomass was taken to be constant. True growth rate dependency of maintenance can be implemented by relaxing either of these assumptions.     

A dynamic mathematical model of the chemostat

A number of experimental studies on the dynamic, behavior of the chemostat have shown that the specific growth rate does not, instantaneously adjust to changes in the concentration of limiting substrate in the chemostat following disturbances in the steady state input limiting substrate concentration or in the steady state dilution rate. Instead of an instantaneous response, as would be predicted by the Monod equation, experimental studies have shown that the specific growth rate experiences a dynamic lag in responding to the changes in the concentration of limiting substrate in the culture vessel. The observed dynamic lag has been recognized by researchers in such terms as an inertial phenomenon and as a hysteresis effect, but as yet a systems engineering approach has not been applied to the observed data. The present paper criticizes the use of the Monod equation as a dynamic relationship and offers as an alternative a dynamic equation relating specific growth rate to the limiting substrate concentration in the chemostat. Following the development of equations, experimental methods of evaluating parameters are discussed. Dynamic responses of analog simulations (incorporating the newly derived equations) are compared with the dynamic responses predicted by the Monod equation and with the dynamic responses of experimental chemostats.     

Chemostat cultures of yeasts, continuous culture fundamentals and simple unstructured mathematical models.

Fundamental aspects of chemostat cultures are reviewed. Using yeast cultures as examples, it is shown that steady states in chemostats may be predicted quantitatively by combining the correct number of unstructured kinetic models with expressions for existing stoichiometric constraints. The necessary number of such kinetic models corresponds to the number of limiting substrates and increases with the number of different metabolic pathways available to the strain. This is demonstrated by an experimental comparison of yeast growth limited by glucose alone for which metabolism is oxidative, and growth doubly limited by both glucose and oxygen, which occurs according to an oxido-reductive metabolism. The steady state data for such experiments can in principle be predicted based on a minimal amount of information by a simple stoichiometric model. It represents the overall stoichiometry of growth by a superposition of a fully oxidative and a fully reductive growth reaction and uses the concept of "aerobicity" to characterize the relative importance of the two reactions.     

On the control mechanisms of the nitrite level in <Emphasis Type="Italic">Escherichia coli</Emphasis> cells: the mathematical model

Background

Due to a high toxicity of nitrite and its metabolites, it is of high interest to study mechanisms underlying the low NO₂ level maintenance in the cell. During anaerobic growth of Escherichia coli the main nitrite-reducing enzymes are NrfA and NirB nitrite reductases. NrfA reductase is localized in the cell periplasm and uses NO₂ as an electron acceptor to create a proton gradient; NirB reductase is restricted to the cytoplasm and metabolizes excessive nitrite inside the cell, the uptake of which is mediated by the transporter protein NirC. While it is known that these three systems, periplasmic, cytoplasmic and transport, determine nitrite uptake and assimilation in the cell as well as its excretion, little is known about their co-ordination.

Results

Using a mathematical model describing the nitrite utilization in E. coli cells cultured in a flow chemostat, the role of enzymes involved in nitrite metabolism and transport in controlling nitrite intracellular levels was investigated. It was demonstrated that the model adapted to the experimental data on expression of nrfA and nirB genes encoding NrfA and NirB nitrite reductases, can describe nitrite accumulation kinetics in the chemostat in the millimolar range of added substrate concentrations without any additional assumptions. According to the model, in this range, low intracellular nitrite level, weakly dependent on its concentration in the growth media, is maintained (mcM). It is not sufficient to consider molecular-genetic mechanisms of NrfA reductase activity regulation to describe the nitrite accumulation dynamics in the chemostat in the micromolar range (≤1 mM) of added nitrite concentrations. Analysis of different hypotheses has shown that the mechanism of local enzyme concentration change due to membrane potential-induced diffusion from the cytoplasm to the periplasm at low nitrite levels is sufficient to explain the nitrite accumulation dynamics in the chemostat.

Conclusions

At nitrite concentrations in the media more than 2 mM, the model adapted to the experimental data on nitrite utilization dynamics in E. coli cells cultured in the flow chemostat demonstrates the largest contribution of genetic mechanisms involved in nrf and nir operons activity regulation to the control of nitrite intracellular levels. The model predicts a significant contribution of the membrane potential to the periplasmic NrfA nitrite reductase activity regulation and nitrite utilization dynamics at substrate concentrations ≤1 mM.

     

Physiological analysis of Methylobacterium extorquens AM1 grown in continuous and batch cultures

Chemostat cultures of Methylobacterium extorquens AM1 grown on methanol or succinate at a range of dilution rates were compared to batch cultures in terms of enzyme levels, poly-β-hydroxybutyrate content, and intracellular concentrations of adenine and pyridine nucleotides. In both chemostat and batch cultures, enzymes specific to C₁ metabolism were up-regulated during growth on methanol and down-regulated during growth on succinate, polyhydroxybutyrate levels were higher on succinate, intracellular ATP levels and the energy charge were higher during growth on methanol, while the pools of reducing equivalents were higher during growth on succinate. For most of the tested parameters, little alteration occurred in response to growth rate. Overall, we conclude that the chemostat cultivation conditions developed in this study roughly mimic the growth in batch cultures, but provide a better control over the culturing conditions and a better data reproducibility, which are important for integrative functional studies. This study provides baseline data for future work using chemostat cultures, defining key similarities and differences in the physiology compared to existing batch culture data.     

The pool of ADP and ATP regulates anaerobic product formation in resting cells of Lactococcus lactis

Lactococcus lactis grows homofermentatively on glucose, while its growth on maltose under anaerobic conditions results in mixed acid product formation in which formate, acetate, and ethanol are formed in addition to lactate. Maltose was used as a carbon source to study mixed acid product formation as a function of the growth rate. In batch and nitrogen-limited chemostat cultures mixed acid product formation was shown to be linked to the growth rate, and homolactic fermentation occurred only in resting cells. Two of the four lactococcal strains investigated with maltose, L. lactis 65.1 and MG1363, showed more pronounced mixed acid product formation during growth than L. lactis ATCC 19435 or IL-1403. In resting cell experiments all four strains exhibited homolactic fermentation. In resting cells the intracellular concentrations of ADP, ATP, and fructose 1,6-bisphosphate were increased and the concentration of P(i) was decreased compared with the concentrations in growing cells. Addition of an ionophore (monensin or valinomycin) to resting cultures of L. lactis 65.1 induced mixed acid product formation concomitant with decreases in the ADP, ATP, and fructose 1,6-bisphosphate concentrations. ADP and ATP were shown to inhibit glyceraldehyde-3-phosphate dehydrogenase, lactate dehydrogenase, and alcohol dehydrogenase in vitro. Alcohol dehydrogenase was the most sensitive enzyme and was totally inhibited at an adenine nucleotide concentration of 16 mM, which is close to the sum of the intracellular concentrations of ADP and ATP of resting cells. This inhibition of alcohol dehydrogenase might be partially responsible for the homolactic behavior of resting cells. A hypothesis regarding the level of the ATP-ADP pool as a regulating mechanism for the glycolytic flux and product formation in L. lactis is discussed.     

The energetics of embryonic growth and development. I. Oxygen consumption, biomass growth, and heat production.

A quantitative phenomenological model to describe the relationships between biomass growth rate, oxygen consumption, and heat production in developing embryos has been developed and tested using a wide range of experimental data. The model employs generalized material and energy balances, principles of enzyme kinetics, and an overall metabolic model scheme based on known biochemical principles. The phosphorylation concentration ratio of ATP and ADP occurs naturally and becomes a significant parameter in the analysis. The model is applied to the growth of Escherichia coli, Oryzias latipes, chick spinal cord, and whole chicken eggs. Excellent agreement between the model and the experimental data is obtained. In a succeeding paper (Part II) environmental effects and growth efficiency are discussed.     

Metabolic and energetic control of Pseudomonas mendocina growth during transitions from aerobic to oxygen-limited conditions in chemostat cultures.

Several metabolic fluxes were analyzed during gradual transitions from aerobic to oxygen-limited conditions in chemostat cultures of Pseudomonas mendocina growing in synthetic medium at a dilution rate of 0.25 h-1. P. mendocina growth was glucose limited at high oxygen partial pressures (70 and 20% pO2) and exhibited an oxidative type of metabolism characterized by respiratory quotient (RQ) values of 1.0. A similar RQ value was obtained at low pO2 (2%), and detectable levels of acetic, formic, and lactic acids were determined in the extracellular medium. RQs of 0.9 +/- 0.12 were found at 70% pO2 for growth rates ranging from 0.025 to 0.5 h-1. At high pO2, the control coefficients of oxygen on catabolic fluxes were 0.19 and 0.22 for O2 uptake and CO2 production, respectively. At low pO2 (2%), the catabolic and anabolic fluxes were highly controlled by oxygen. P. mendocina showed a mixed-type fermentative metabolism when nitrogen was flushed into chemostat cultures. Ethanol and acetic, lactic, and formic acids were excreted and represented 7.5% of the total carbon recovered. Approximately 50% of the carbon was found as uronic acids in the extracellular medium. Physiological studies were performed under microaerophilic conditions (nitrogen flushing) in continuous cultures for a wide range of growth rates (0.03 to 0.5 h-1). A cell population, able to exhibit a near-maximum theoretical yield of ATP (YmaxATP = 25 g/mol) with a number of ATP molecules formed during the transfer of an electron towards oxygen along the respiration chain (P/O ratio) of 3, appears to have adapted to microaerophilic conditions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamic analysis of a microbial process: A systems engineering approach

A general mathematical model of the chemostat system is developed in order to define an experimental program of dynamic testing. A glucose-limited culture ofSaccharomyces cerevisiae was grown in a chemostat using chemically defined medium. The chemostat was perturbed from an initial steady state by changes in input glucose concentration, dilution rate, pH, and temperature. Dynamic responses of cell mass, glucose, cell number, RNA, and protein concentrations were measured. A number of simulation techniques were used in developing a dynamic mathematical model and in comparing the developed model with experimental data as well as the Monod model. The resulting model was found to be quantitatively accurate and superior to the Monod model. The developed model was interpreted in the light of cell physiology. Adjustment of intracellular RNA fraction was found to be rate limiting in acceleration of cell specific growth rate.     

Metabolic and energetic control of Pseudomonas mendocina growth during transitions from aerobic to oxygen-limited conditions in chemostat cultures.

Several metabolic fluxes were analyzed during gradual transitions from aerobic to oxygen-limited conditions in chemostat cultures of Pseudomonas mendocina growing in synthetic medium at a dilution rate of 0.25 h-1. P. mendocina growth was glucose limited at high oxygen partial pressures (70 and 20% pO2) and exhibited an oxidative type of metabolism characterized by respiratory quotient (RQ) values of 1.0. A similar RQ value was obtained at low pO2 (2%), and detectable levels of acetic, formic, and lactic acids were determined in the extracellular medium. RQs of 0.9 +/- 0.12 were found at 70% pO2 for growth rates ranging from 0.025 to 0.5 h-1. At high pO2, the control coefficients of oxygen on catabolic fluxes were 0.19 and 0.22 for O2 uptake and CO2 production, respectively. At low pO2 (2%), the catabolic and anabolic fluxes were highly controlled by oxygen. P. mendocina showed a mixed-type fermentative metabolism when nitrogen was flushed into chemostat cultures. Ethanol and acetic, lactic, and formic acids were excreted and represented 7.5% of the total carbon recovered. Approximately 50% of the carbon was found as uronic acids in the extracellular medium. Physiological studies were performed under microaerophilic conditions (nitrogen flushing) in continuous cultures for a wide range of growth rates (0.03 to 0.5 h-1). A cell population, able to exhibit a near-maximum theoretical yield of ATP (YmaxATP = 25 g/mol) with a number of ATP molecules formed during the transfer of an electron towards oxygen along the respiration chain (P/O ratio) of 3, appears to have adapted to microaerophilic conditions.(ABSTRACT TRUNCATED AT 250 WORDS)     

A theoretical and empirical investigation of delayed growth response in the continuous culture of bacteria

When the growth of bacteria in a chemostat is controlled by limiting the supply of a single essential nutrient, the growth rate is affected both by the concentration of this nutrient in the culture medium and by the amount of time that it takes for the chemical and physiological processes that result in the production of new biomass. Thus, although the uptake of nutrient by cells is an essentially instantaneous process, the addition of new biomass is delayed by the amount of time that it takes to metabolize the nutrient. Mathematical models that incorporate this "delayed growth response" (DGR) phenomenon have been developed and analysed. However, because they are formulated in terms of parameters that are difficult to measure directly, these models are of limited value to experimentalists. In this paper, we introduce a DGR model that is formulated in terms of measurable parameters. In addition, we provide for this model a complete set of criteria for determining persistence versus extinction of the bacterial culture in the chemostat. Specifically, we show that DGR plays a role in determining persistence versus extinction only under certain ranges of chemostat operating parameters. It is also shown, however, that DGR plays a role in determining the steady-state nutrient and bacteria concentrations in all instances of persistence. The steady state and transient behavior of solutions of our model is found to be in agreement with data that we obtained in growing Escherichia coli 23716 in a chemostat with glucose as a limiting nutrient. One of the theoretical predictions of our model that does not occur in other DGR models is that under certain conditions a large delay in growth response might actually have a positive effect on the bacteria's ability to persist.