An energetic model for oxygen-limited metabolism

Microbial production of 2,3-butanediol by Klebsiella oxytoca occurs under conditions of an oxygen limitation. The extent to which substrate is oxidized to 2,3-butanediol and its coproducts, (acetic acid, acetoin, and ethanol) and the relative flow rates of substrate to energetic and biosynthetic pathways are controlled by the degree of oxygen limitation. Two energetic relationships which describe the response to an oxygen limitation have been derived. The first relationship describes the coupling between growth and energy production observed under oxygen-limited conditions. This allows calculation of energetic parameters and modeling of the cell mass and substrate profiles in terms of the degree of oxygen limitation only. The second relationship describes the average degree of oxidation and the rate of the end-product flow. The model has been tested with both batch and continuous culture. During these kinetic studies, two phases of growth have been observed: energy-coupled growth, which was described above; and, energy-uncoupled growth, which arises when the degree of oxygen limitation reaches a critical value. Optimal culture performance with respect to 2,3-butanediol productivity occurs during energy-coupled growth. (c) 1993 John Wiley & Sons, Inc.     

A kinetic model for product formation of microbial and mammalian cells

Growth of microbial and mammalian cells can be classified into substrate-limited and substrate-sufficient growth according to the relative availability of the substrate (carbon and energy source) and other nutrients. It has been observed for a number of microbial and mammalian cells that the consumption rate of substrate and energy (ATP) is generally higher under substratesufficient conditions than under substrate limitation. Accordingly, the product formation under substrate excess often exhibits different patterns from those under substrate limitation. The extent of increase or decrease in product formation may depend not only on the nature of limitation and cell growth rate but also on the residual substrate concentration in a relatively wide range. The product formation kinetic models existing in literature cannot describe these effects. In this study, the Luedeking-Piret kinetic is extended to include a term describing the effect of residual substrate concentration. The extended model has a similar structure to the kinetic model for substrate and energy consumption rate recently proposed by Zeng and Deckwer. The applicability of the extended model is demonstrated with three microbial cultures for the production of primary metabolites and three hybridoma cell cultures for the production of ammonia and lactic acid over a wide range of substrate concentration. The model describes the product formation in all these cultures satisfactorily. Using this model, the range of residual substrate concentration, in which the product formation is affected, can be quantitatively assessed. (c) 1995 John Wiley & Sons, Inc.     

Fundamental Escherichia coli biochemical pathways for biomass and energy production: creation of overall flux states

We have previously shown that the metabolism for most efficient cell growth can be realized by a combination of two types of elementary modes. One mode produces biomass while the second mode generates only energy. The identity of the four most efficient biomass and energy pathway pairs changes, depending on the degree of oxygen limitation. The identification of such pathway pairs for different growth conditions offers a pathway-based explanation of maintenance energy generation. For a given growth rate, experimental aerobic glucose consumption rates can be used to estimate the contribution of each pathway type to the overall metabolic flux pattern. All metabolic fluxes are then completely determined by the stoichiometries of involved pathways defining all nutrient consumption and metabolite secretion rates. We present here equations that permit computation of network fluxes on the basis of unique pathways for the case of optimal, glucose-limited Escherichia coli growth under varying levels of oxygen stress. Predicted glucose and oxygen uptake rates and some metabolite secretion rates are in remarkable agreement with experimental observations supporting the validity of the presented approach. The entire most efficient, steady-state, metabolic rate structure is explicitly defined by the developed equations without need for additional computer simulations. The approach should be generally useful for analyzing and interpreting genomic data by predicting concise, pathway-based metabolic rate structures.     

A simulation model for the continuous production of acetoin and butanediol using Bacillus subtilis with integrated pervaporation separation

The potential for producing acetoin and butanediol with a Bacillus subtilis strain was investigated with continuous culture using molasses as carbon substrate. The steady-state results were influenced by both oxygen and undetermined limiting compounds. Employing the known metabolic pathways, four overall stoichiometry relations were used with an energetic assumption on the energy requirements for biomass formation to establish a linear relations were used with an energetic assumption on the energy requirements for biomass formation to establish a linear relation between the overall rates, whose parameters were determined by linear regression. This provided a relationship for the product formation rate. The chemostat culture data were described with a growth kinetics model, which included limitation by molasses and oxygen as well as diauxic effects and product inhibition. The biokinetics model was combined with an experimentally verified model for the membrane Pervaporation. From this combined model were determined the influence of the membrane characteristics (enrichment factors and membrane area) and the dilution rate on the performance of the integrated process. Simulations revealed that an increase of the enrichment factor, possible by membrane improvement, would have counteracting influences, owing to decreased product inhibition but with lower biomass concentration. (c) 1993 Wiley & Sons, Inc.     

Phosphate-limited culture of Azotobacter vinelandii.

Batch cultures of Azotobacter vinelandii grown in phosphate-deficient media were compared with control cultures grown in phosphate-sufficient media. Phosphate limitation was assessed by total cell yield and by growth kinetics. Although cell protein, nucleic acids, and early growth rate were unaffected by phosphate deficiency, cell wall structure, oxygen uptake, and cell viability were significantly affected. Also, phosphate-limited cells contained much larger amounts of poly-beta-hydroxybutyric acid but lower adenylate nucleotide energy charge than did control cells. The ratio of adenosine 5'-triphosphate to adenosine 5'-diphosphate was much lower in phosphate-deficient cells. The data indicate a substrate saving choice of three metabolic pathways available to this organism under different growth conditions.     

Respirometric evaluation and modeling of glucose utilization by Escherichia coli under aerobic and mesophilic cultivation conditions

The study presents a mechanistic model for the evaluation of glucose utilization by Escherichia coli under aerobic and mesophilic growth conditions. In the first step, the experimental data was derived from batch respirometric experiments conducted at 37 degrees C, using two different initial substrate to microorganism (S(0)/X(0)) ratios of 15.0 and 1.3 mgCOD/mgSS. Acetate generation, glycogen formation and oxygen uptake rate profile were monitored together with glucose uptake and biomass increase throughout the experiments. The oxygen uptake rate (OUR) exhibited a typical profile accounting for growth on glucose, acetate and glycogen. No acetate formation (overflow) was detected at low initial S(0)/X(0) ratio. In the second step, the effect of culture history developed under long-term growth limiting conditions on the kinetics of glucose utilization by the same culture was evaluated in a sequencing batch reactor (SBR). The system was operated at cyclic steady state with a constant mean cell residence time of 5 days. The kinetic response of E.coli culture was followed by similar measurements within a complete cycle. Model calibration for the SBR system showed that E. coli culture regulated its growth metabolism by decreasing the maximum growth rate (lower microH) together with an increase of substrate affinity (lower K(S)) as compared to uncontrolled growth conditions. The continuous low rate operation of SBR system induced a significant biochemical substrate storage capability as glycogen in parallel to growth, which persisted throughout the operation. The acetate overflow was observed again as an important mechanism to be accounted for in the evaluation of process kinetics.     

Dynamic flux balance analysis of the metabolism of Saccharomyces cerevisiae during the shift from fully respirative or respirofermentative metabolic states to anaerobiosis

Dynamic flux balance analysis was utilized to simulate the metabolic behaviour of initially fully respirative and respirofermentative steady-state cultures of Saccharomyces?cerevisiae during sudden oxygen depletion. The hybrid model for the dynamic flux balance analysis included a stoichiometric genome-scale metabolic model as a static part and dynamic equations for the uptake of glucose and the cessation of respirative metabolism. The yeast consensus genome-scale metabolic model [Herrg?rd MJ et?al. (2008) Nat Biotechnol26, 1155-1160; Dobson PD et?al. (2010) BMC Syst Biol4, 145] was refined with respect to oxygen-dependent energy metabolism and further modified to reflect S.?cerevisiae anabolism in the absence of oxygen. Dynamic flux balance analysis captured well the essential features of the dynamic metabolic behaviour of S.?cerevisiae during adaptation to anaerobiosis. Modelling and simulation enabled the identification of short time-scale flux distribution dynamics under the transition to anaerobic metabolism, during which the specific growth rate was reduced, as well as longer time-scale process dynamics when the specific growth rate recovered. Expression of the metabolic genes was set into the context of the identified dynamics. Metabolic gene expression responses associated with the specific growth rate and with the cessation of respirative metabolism were distinguished.     

Myo-inositol transport in Saccharomyces cerevisiae.

myo-Inositol uptake in Saccharomyces cerevisiae was dependent on temperature, time, and substrate concentration. The transport obeyed saturation kinetics with an apparent Km for myo-inositol of 0.1 mM, myo-Inositol analogs, such as scyllo-inositol, 2-inosose, mannitol, and 1,2-cyclohexanediol, had no effect on myo-inositol uptake, myo-Inositol uptake required metabolic energy. Removal of D-glucose resulted in a loss of activity, and azide and cyanide ions were inhibitory. In the presence of D-glucose, myo-inositol was accumulated in the cells against a concentration gradient. A myo-inositol transport mutant was isolated from UV-mutagenized S. cerevisiae cells using the replica-printing technique. The defect in myo-inositol uptake was due to a single nuclear gene mutation. The activities of L-serine and D-glucose transport were not affected by the mutation. Thus it was shown that S. cerevisiae grown under the present culture conditions possessed a single and specific myo-inositol transport system. myo-Inositol transport activity was reduced by the addition of myo-inositol to the culture medium. The activity was reversibly restored by the removal of myo-inositol from the medium. This restoration of activity was completely abolished by cycloheximide.     

Utilization of the tricarboxylic acid cycle, a reactor design criterion for the microaerobic production of 2,3-butanediol

A new parameter, the relative utilization of tricarboxylic acid (TCA) cycle beta, is introduced to quantitatively account for the involvement of fermentation pathways and TCA cycle in the utilization of oxygen under oxygen-limiting (microaerobic) conditions. With the facultative anaerobe Enterobacter aerogenes, which produces 2,3-butanediol, a method is proposed to calculate beta from measurement of metabolites and exhaust gas. In continuous culture beta was found to be small under oxygen limitation, indicating that the fermentation pathways were preferred over the TCA cycle and oxygen was almost entirely consumed through oxidation of reduced nicotinamide adenine dinucleotide (NADH(2)) released by fermentation under these conditions. The increase of beta at high oxygen supply revealed a saturation of oxygen utilization through fermentation pathways. It could be concluded that, for the optimal performance of a microaerobic culture, oxygen uptake rate must be kept at such a level that as much NADH(2) as possible from fermentation pathways is oxidized by oxygen, and at the same time the utilization of TCA cycle is kept at a minimum. As the dynamics of the microaerobic culture can be fast, a significant effect of reactor hydrodynamics, i.e., mixing, on the overall performance can be expected. This was confirmed experimentally, and the parameter beta proved to be a useful reactor design criterium for the microaerobic cultivation. (c) 1992 John Wiley & Sons, Inc.     

Effect of growth rate and substrate limitation on the composition and structure of the cell wall of Saccharomyces cerevisiae

1. A study was made of the composition and structure of walls isolated from yeast grown in continuous culture at different rates, under three conditions of glucose limitation in which the concentrations of glucose and ammonium sulphate in the medium and the oxygen-transfer rate in the culture were varied, and one condition of NH(4) (+) limitation. 2. The contents of total glucan and total mannan in the walls were relatively little affected by the growth rate under any of the four sets of conditions. The phosphorus and protein contents of walls from yeast grown under each of the four conditions increased as the growth rate was decreased. Walls from yeast grown under NH(4) (+) limitation contained only half as much protein as walls from cells grown under glucose limitation. The proportion of lipid was greatest in walls from yeast grown under NH(4) (+) limitation. 3. A procedure was devised for fractionating isolated walls, based on the ease with which the glucan and mannan were extracted with water and with hot and cold 6% (w/v) potassium hydroxide solution. The percentage of glucan, mannan, protein and phosphorus in each of the fractions was affected by the rate of growth and by the nature of the substrate limitation. 4. The beta-fructofuranosidase activities of yeast grown under glucose limitation increased as the growth rate was lowered, but decreased at very low growth rates. The effects at low growth rates were probably due to repression of enzyme synthesis by residual glucose in the culture filtrate. The beta-fructofuranosidase activities of yeast grown under NH(4) (+) limitation were much lower than those from yeast grown under any of the conditions of glucose limitation. 5. Yeast cells grown at any of the rates under NH(4) (+) limitation were longer and thinner than those grown at the same rate under any of the conditions of glucose limitation. Mean cell volumes were dependent on growth rate but not on the nature of the substrate limitation. 6. Electron micrographs of thin sections of isolated walls showed that cells grown under NH(4) (+) limitation had a more porous structure than those from cells grown under any of the conditions of glucose limitation.     

Inhibition kinetics of phenol degradation from unstable steady-state data

Multiplicity of steady states of a continuous culture with an inhibitory substrate was used to estimate kinetic parameters under steady-state conditions. A continuous culture of Pseudomonas cepacia G4, using phenol as the sole source of carbon and energy, was overloaded by increasing the dilution rate above the critical dilution rate. The culture was then stabilized in the inhibitory branch by a proportional controller using the carbon dioxide concentration in the reactor exhaust gas as the controlled variable and the dilution rate as the manipulated variable. By variation of the set point, several unstable steady states in the inhibitory branch were investigated and the specific phenol conversion rates calculated. In addition, phenol degradation was investigated under substrate limitation (chemostat operation).The results show that the phenol degradation by P. cepacia can be described by the same set of inhibition parameters under substrate limitation and under high substrate concentrations in the inhibitory branch. Biomass yield and maintenance coefficients were identical. Fitting of the data to various inhibition models resulted in the best fit for the Yano and Koga equation. The well-known Haldane model, which is most often used to describe substrate inhibition by phenol, gave the poorest fit. The described method allows a precise data estimation under steady-state conditions from the maximum of the biological reaction rate up to high substrate concentrations in the inhibitory branch. Inhibition parameter estimation by controlling unstable steady states may thus be useful in avoiding discrepancies between data generated by batch runs and their application to continuous cultures which have been often described in the literature. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 54: 567-576, 1997.     

Modeling threshold phenomena, metabolic pathways switches and signals in chemostat-cultivated cells: the Crabtree effect in Saccharomyces cerevisiae

The long-term Crabtree effect in Saccharomyces cerevisiae cultivated in aerobic chemostat at steady state has been studied for three different substrate concentrations in the feed of the bioreactor (data: J. Gen. Microbiol., 129 (1983) 653). We have shown that a model using two ways of transport/metabolization (T/M) of hyperbolic form, with high and low affinity for the substrate, allowed to represent correctly the main characteristics of the phenomenon. The model is based on an explicit form of the T/M kinetics when the bioreactor is considered as a polyphasic dispersed system (PDS). Mass balances analysis also allows to quantify the critical dilution rate value (threshold), Dc, of the transition between respiratory and respirofermentative mode, for which ethanol is produced. A good approximation for the threshold is Dc = V(S)0 Y(Xc, S) where Y(Xc,S) is the average yield coefficient before transition and V(S)0, the maximum specific rate of high affinity T/M pathway. The theoretical value is 0.3 h(-1), and is equal to the experimental value. We thus show in a quantitative way that the transition depends both on culture conditions (global characteristic of the system) and on strain properties (intrinsic characteristic of the microorganism as well). Using two different methods to calculate the residual substrate has carried out the comparison between the simulations end the experimental data. This allowed showing that the latter is not well represented by Monod's model and has confirmed that the affinity for the substrate varies according to the biomass. We have then shown how to calculate the most important specific rates (or metabolic flux) related to biomass, ethanol, oxygen, hydrogen, respiratory and fermentative CO(2) and H(2)O within the cellular phase. It has appeared that the oxygen uptake rate directly depends on high-affinity T/M pathway. This let us think that the regulation of the Crabtree effect in S. cerevisiae depends on the saturation of some glucose metabolization and transport pathways rather than on saturation of the respiratory chains. The specific rates analysis has also allowed us to show, at least in this case, that the metabolization rate (biosynthesis+fueling) had its maximum value on the whole dilution rates interval; metabolites excretion (ethanol and fermentative CO(2)) only intervenes to drain a "surplus" glucose flux. As a consequence, the transport capacity must be higher than the one of metabolization. Maximization of the metabolization specific rate could then be used as an optimization criterion in the stoichiometric calculation of metabolic flux (and not the specific growth rate maximization because growth is limited in a chemostat (mu = D)). We have also shown that the mass balances based on the T/M processes are in agreement with molar and elementary balances of the general stoichiometric equation for glucose respiration and fermentation under aerobic conditions. Thanks to the specific rates calculating the stoichiometric coefficients has done this. The total mass balance difference does not exceed 4%, which is compatible with the experimental carbon balance. Finally, we have emphasized that the ratio of biosynthesis flux and metabolization flux is constant before and after transition. This observation could be applied as soon as the free substrate concentration in the cellular phase is low. The paper succinctly describes the former theoretical results on which the model is built and sufficiently explains the algorithm for straightforward implementation.     

A general model for aerobic yeast growth: batch growth

A general model for aerobic yeast growth in batch culture is presented. It is based on the concept that the aerobic metabolism of all yeasts is determined by the relative sizes of the transport rate of sugar into the cell and the transport rate of respiratory intermediates into the mitochondrion. If the rate of sugar uptake rate exceeds the rate of transport of respiratory intermediates into the mitochondrion (as in Saccharomyces cerevisiae, S. uvarum, and S. pombe), the metabolism exhibits the features of ethanol excretion and limited specific oxygen uptake rate. If the rate of transport of respiratory intermediates into the mitochondrion is of the same order as the transport of sugar into the cell (as in Candida utilis), the metabolism is characterized by little or no ethanol excretion and a much higher specific oxygen uptake rate. Batch data from an extensive range of yeast and carbon sources is used to illustrate the use of this model. The ability of this model to fit such an extensive range of experimental data suggests that it can be used as a generalized model for aerobic yeast growth.     

Complete microbial degradation of both enantiomers of the chiral herbicide mecoprop [(RS)-2-(4-chloro-2-methylphenoxy)propionic acid] in an enantioselective manner by Sphingomonas herbicidovorans sp. nov.

Sphingomonas herbicidovorans MH (previously designated Flavobacterium sp. strain MH) was able to utilize the chiral herbicide (RS)-2-(4-chloro-2-methylphenoxy)propionic acid (mecoprop) as the sole carbon and energy source. When strain MH was offered racemic mecoprop as the growth substrate, it could degrade both the (R) and the (S) enantiomer to completion, as shown by biomass formation, substrate consumption, and stoichiometric chloride release. However, the (S) enantiomer disappeared much faster from the culture medium than the (R) enantiomer. These results suggest the involvement of specific enzymes for the degradation of each enantiomer. This view was substantiated by the fact that resting cells of strain MH grown on (S)-mecoprop were able to degrade the (S) but not the (R) enantiomer of mecoprop. Accordingly, resting cells of strain MH grown on (R)-mecoprop preferentially metabolized the (R) enantiomer. Nevertheless, such cells could transform (S)-mecoprop at low rates. Oxygen uptake rates with resting cells confirmed the above view, as oxygen consumption was strongly dependent on the growth substrate. Cells grown on (R)-mecoprop showed oxygen uptake rates more than two times higher upon incubation with the (R) than upon incubation with the (S) enantiomer and vice versa.     

Quantitative comparison of transient growth of Saccharomyces cerevisiae,Saccharomyces kluyveri,and Kluyveromyces lactis

A multitude of metabolic regulations occur in yeast, particularly under dynamic process conditions, such as under sudden glucose excess. However, quantification of regulations and classification of yeast strains under these conditions have yet to be elucidated, which requires high-frequency and consistent quantification of the metabolic response. The present study aimed at quantifying the dynamic regulation of the central metabolism of strains Saccharomyces cerevisiae, S. kluyveri, and Kluyveromyces lactis upon sudden glucose excess, accomplished by a shift-up in dilution rate inside of the oxidative region using a small metabolic flux model. It was found that, under transient growth conditions, S. kluyveri behaved like K. lactis, while classification using steady-state conditions would position S. kluyveri close to S. cerevisiae. For transient conditions and based on the observation whether excess glucose is initially used for catabolism (energy) or anabolism (carbon), we propose to classify strains into energy-driven, such as S. cerevisiae, and carbon-driven, such as S. kluyveri and K. lactis, strains. Furthermore, it was found that the delayed onset of fermentative catabolism in carbon-driven strains is a consequence of low catabolic flux and the initial shunt of glucose in non-nitrogen-containing biomass constituents. The MFA model suggests that energy limitation forced the cell to ultimately increase catabolic flux, while the capacity of oxidative catabolism is not sufficient to process this flux oxidatively. The combination of transient experiments and its exploitation with reconciled intrinsic rates using a small metabolic model could corroborate earlier findings of metabolic regulations, such as tight glucose control in carbon-driven strains and transient changes in biomass composition, as well as explore new regulations, such as assimilation of ethanol before glucose. The benefit from using small metabolic flux models is the richness of information and the enhanced insight into intrinsic metabolic pathways without a priori knowledge of adaptation kinetics. Used in an online context, this approach serves as an efficient tool for strain characterization and physiological studies.     

Carbon conversion efficiency and limits of productive bacterial degradation of methyl tert-butyl ether and related compounds

The utilization of the fuel oxygenate methyl tert-butyl ether (MTBE) and related compounds by microorganisms was investigated in a mainly theoretical study based on the Y(ATP) concept. Experiments were conducted to derive realistic maintenance coefficients and K(s) values needed to calculate substrate fluxes available for biomass production. Aerobic substrate conversion and biomass synthesis were calculated for different putative pathways. The results suggest that MTBE is an effective heterotrophic substrate that can sustain growth yields of up to 0.87 g g(-1), which contradicts previous calculation results (N. Fortin et al., Environ. Microbiol. 3:407-416, 2001). Sufficient energy equivalents were generated in several of the potential assimilatory routes to incorporate carbon into biomass without the necessity to dissimilate additional substrate, efficient energy transduction provided. However, when a growth-related kinetic model was included, the limits of productive degradation became obvious. Depending on the maintenance coefficient m(s) and its associated biomass decay term b, growth-associated carbon conversion became strongly dependent on substrate fluxes. Due to slow degradation kinetics, the calculations predicted relatively high threshold concentrations, S(min), below which growth would not further be supported. S(min) strongly depended on the maximum growth rate mu(ma)(x), and b and was directly correlated with the half maximum rate-associated substrate concentration K(s), meaning that any effect impacting this parameter would also change S(min). The primary metabolic step, catalyzing the cleavage of the ether bond in MTBE, is likely to control the substrate flux in various strains. In addition, deficits in oxygen as an external factor and in reduction equivalents as a cellular variable in this reaction should further increase K(s) and S(min) for MTBE.