Mathematical models of microbial growth and competition in the chemostat regulated by cell-bound extracellular enzymes

A mathematical model of growth and competitive interaction of microorganisms in the chemostat is analyzed. The growth-limiting nutrient is not in a form that can be directly assimilated by the microorganisms, and must first be transformed into an intermediate product by cell-bound extracellular enzymes. General monotone functions, including Michaelis-Menten and sigmoidal response functions, are used to describe nutrient conversion and growth due to consumption of the intermediate product. It is shown that the initial concentration of the species is an important determining factor for survival or washout. When there are two species whose growth is limited by the same nutrient, three different modes of competition are described. Competitive coexistence steady states are shown to be possible in two of them, but they are always unstable. In all of our numerical simulations, the system approaches a steady state corresponding to the washout of one or both of the species from the chemostat.Research supported by NSF grant DMS-90-96279Research supported by NSERC grant A-9358     

Stability and stabilization for models of chemostats with multiple limiting substrates

We study chemostat models in which multiple species compete for two or more limiting nutrients. First, we consider the case where the nutrient flow and species removal rates and input nutrient concentrations are all given as positive constants. In that case, we use Brouwer degree theory to give conditions guaranteeing that the models admit globally asymptotically stable componentwise positive equilibrium points, from all componentwise positive initial states. Then we use the results to develop stabilization theory for a class of controlled chemostats with two or more limiting nutrients. For cases where the dilution rate and input nutrient concentrations can be selected as controls, we prove that many different componentwise positive equilibria can be made globally asymptotically stable. This extends the existing control results for chemostats with one limiting nutrient. We demonstrate our methods in simulations.     

Stability and stabilization for models of chemostats with multiple limiting substrates

We study chemostat models in which multiple species compete for two or more limiting nutrients. First, we consider the case where the nutrient flow and species removal rates and input nutrient concentrations are all given as positive constants. In that case, we use Brouwer degree theory to give conditions guaranteeing that the models admit globally asymptotically stable componentwise positive equilibrium points, from all componentwise positive initial states. Then we use the results to develop stabilization theory for a class of controlled chemostats with two or more limiting nutrients. For cases where the dilution rate and input nutrient concentrations can be selected as controls, we prove that many different componentwise positive equilibria can be made globally asymptotically stable. This extends the existing control results for chemostats with one limiting nutrient. We demonstrate our methods in simulations.     

Population dynamics and competition in chemostat models with adaptive nutrient uptake

 The standard Monod model for microbial population dynamics in the chemostat is modified to take into consideration that cells can adapt to the change of nutrient concentration in the chemostat by switching between fast and slow nutrient uptake and growing modes with asymmetric thresholds for transition from one mode to another. This is a generalization of a modified Monod model which considers adaptation by transition between active growing and quiescent cells. Global analysis of the model equations is obtained using the theory of asymptotically autonomous systems. Transient oscillatory population density and hysteresis growth pattern observed experimentally, which do not occur for the standard Monod model, can be explained by such adaptive mechanism of the cells. Competition between two species that can switch between fast and slow nutrient uptake and growing modes is also considered. It is shown that generically there is no coexistence steady state, and only one steady state, corresponding to the survival of at most one species in the chemostat, is a local attractor. Numerical simulations reproduce the qualitative feature of some experimental data which show that the population density of the winning species approaches a positive steady state via transient oscillations while that of the losing species approaches the zero steady state monotonically. Received 4 August 1995; received in revised form 15 December 1995     

Global dynamics of microbial competition for two resources with internal storage

We study a chemostat model that describes competition between two species for two essential resources based on storage. The model incorporates internal resource storage variables that serve the direct connection between species growth and external resource availability. Mathematical analysis for the global dynamics of the model is carried out by using the monotone dynamical system theory. It is shown that the limiting system of the model basically exhibits the familiar Lotka-Volterra alternatives: competitive exclusion, coexistence, and bi-stability, and most of these results can be carried over to the original model. B. Li’s research was partially supported by NSF grants DMS 0211614 and DMS 0616445. H. L. Smith’s research was supported in part by NSF grant DMS 0414270.     

Chaos in a periodically forced chemostat with algal mortality

We study the possibility of chaotic dynamics in the externally driven Droop model. This model describes a phytoplankton population in a chemostat under periodic nutrient supply. Previously, it has been proven under very general assumptions, that such systems are not able to exhibit chaotic dynamics. We show that the simple introduction of algal mortality may lead to chaotic oscillations of algal density in the forced chemostat. Our numerical simulations show that the existence of chaos is intimately related to plankton overshooting in the unforced model. We provide a simple measure, based on stability analysis, for estimating the amount of overshooting. These findings are not restricted to the Droop model but also hold for other chemostat models with mortality. Our results suggest periodically driven chemostats as a simple model system for the experimental verification of chaos in ecology.     

Asymptotic Behavior of a Chemostat Model with Constant Recycle Sludge Concentration

In this work, we study a several species aerobic chemostat model with constant recycle sludge concentration in continuous culture. We reduce the number of parameters by considering a dimensionless model. First, the existence of a global positive uniform attractor for the model with different removal rates is proved using the theory of dissipative dynamical systems. Hence, we investigate the asymptotic behavior of the model under small perturbations using methods of singular perturbation theory and we prove that, in the case of two species in competition, the unique equilibrium which is positive is globally asymptotically stable. Finally, we establish the link between the open problem of the chemostat with different removal rates and monotone functional responses, and our model when two species compete on the same nutrient. We give some numerical simulations to illustrate the results.     

Do spatial effects appear at low dilution rate in chemostat?

The chemostat theory on two species competition has shown that the dilution rate where transition of dominance occurs – transition-dilution rate – is independent of limiting-nutrient concentration. However, we obtained the experimental data indicating that the transition-dilution rate changed with variations in limiting-ammonium concentrations, using the chemostat mixed-culture of the cyanobacterium Microcystis novacekii and the green alga Scenedesmus quadricauda. The transition-dilution rate was dependent on the concentration of limiting ammonium in the influx culture medium. We tried to simulate the experimental results. Though the dilution rate has been considered independent of nutrient concentration, we introduce the effective dilution rate that depends on nutrient concentration (ammonium concentration in this study). A hyperbolic Monod-type function is used to represent the effective dilution rate for each species. The maximum dilution rate of the function is set to be the mechanical dilution rate (nominal dilution rate) of the chemostat culture. The calculation shows that the nominal transition-dilution rate where transition of dominance occur decreases with increased concentration. This simulation is well consistent with our experimental data. These results may suggest that the species-specificity of limiting nutrients, here nitrogen. Or they may imply that the depreciation of nitrogen becomes critical when both dilution rate and concentration are very low, especially for the green algae. In the latter case, spatial effects are induced internally in the ecosystem.     

Mathematical investigations of growth of microorganisms in the gradostat

Microbial ecosystems with spatial distribution of substrate (nutrient) supply in the form of gradients have been studied in a laboratory system called the gradostat, which is a series of coupled chemostats. We investigate analytically and numerically a mathematical model, similar to the one for a single chemostat based on Michaelis-Menten kinetics, of the growth of one species of microorganisms in the gradostat in the presence of one limiting substrate and two limiting complementary substrates. Our analysis predicts various patterns of spatial distribution of microorganisms at steady state and suggests further experiments to be performed with the gradostat.This work was supported by the Deutsche Forschungsgemeinschaft     

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.     

Dynamics of chemostat culture:the effect of a delay in cell response

The effect of introducing a delay between a change in concentration of limiting nutrient and the response in division rate of the population is investigated. It is shown that this might or might not introduce a delay in the control of the chemostat system, depending upon how the uptake rate is assumed to double during the cell cycle. An area in the parameter space where autonomous oscillations will occur for one specific model is found analytically. This result is confirmed by computer simulations.     

Allelopathy of plasmid-bearing and plasmid-free organisms competing for two complementary resources in a chemostat

We consider a model of competition between plasmid-bearing and plasmid-free organisms for two complementary nutrients in a chemostat. We assume that the plasmid-bearing organism produces an allelopathic agent at the cost of its reproductive abilities which is lethal to plasmid-free organism. Our analysis leads to different thresholds in terms of the model parameters acting as conditions under which the organisms associated with the system cannot thrive even in the absence of competition. Local stability of the system is obtained in the absence of one or both the organisms. Also, global stability of the system is obtained in the presence of both the organisms. Computer simulations have been carried out to illustrate various analytical results.     

Practical coexistence of two species in the chemostat - A slow-fast characterization

We show that the chemostat model with two species having different but close break-even concentrations exhibits a slow-fast dynamics. Considering small perturbations about the dilution rate for which break-even concentrations are identical, we use the Fenichel theory to show the coexistence of species for large times. Then we determine the reduced dynamics, which is non-trivial and characterized by the slopes of the growth functions about their break-even concentrations.     

The analysis and regulation for the dynamics of a temperate bacteriophage model

The purpose of this paper is to study the asymptotical behavior of a temperate bacteriophage model in chemostat, which was first proposed by Levin et al. [B.R. Levin, F.M. Stewart, L. Chao, Resource-limited growth, competition and predation: A model and experiment studies with bacteria and bacteriophage, Am. Nat. 125 (1977) 3]. Firstly, a classification for the equilibria of the model and their stability are obtained; secondly, sufficient conditions for uniform persistence are obtained; thirdly, sufficient conditions for the global asymptotic behavior are given, and simulations for the model are presented. The theoretical results show that there are more than eight cases for the classification of the model, and that the decrease (increase) of the nutrient concentration or average lytic time (flow rate) is beneficial to the survival of the sensitive cells. Both the simulated and theoretical results show that there is a possibility of switch phenomena or a periodical outburst of the phages and the lysogens, which is caused by the internal factors rather than by some external environment. Finally, the simulation and regulation of the dynamics of the model with experimental data are presented.     

Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110.

Flux balance models of metabolism use stoichiometry of metabolic pathways, metabolic demands of growth, and optimality principles to predict metabolic flux distribution and cellular growth under specified environmental conditions. These models have provided a mechanistic interpretation of systemic metabolic physiology, and they are also useful as a quantitative tool for metabolic pathway design. Quantitative predictions of cell growth and metabolic by-product secretion that are experimentally testable can be obtained from these models. In the present report, we used independent measurements to determine the model parameters for the wild-type Escherichia coli strain W3110. We experimentally determined the maximum oxygen utilization rate (15 mmol of O2 per g [dry weight] per h), the maximum aerobic glucose utilization rate (10.5 mmol of Glc per g [dry weight] per h), the maximum anaerobic glucose utilization rate (18.5 mmol of Glc per g [dry weight] per h), the non-growth-associated maintenance requirements (7.6 mmol of ATP per g [dry weight] per h), and the growth-associated maintenance requirements (13 mmol of ATP per g of biomass). The flux balance model specified by these parameters was found to quantitatively predict glucose and oxygen uptake rates as well as acetate secretion rates observed in chemostat experiments. We have formulated a predictive algorithm in order to apply the flux balance model to describe unsteady-state growth and by-product secretion in aerobic batch, fed-batch, and anaerobic batch cultures. In aerobic experiments we observed acetate secretion, accumulation in the culture medium, and reutilization from the culture medium. In fed-batch cultures acetate is cometabolized with glucose during the later part of the culture period.(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.     

Allelopathy of plasmid-bearing and plasmid-free organisms competing for two complementary resources in a chemostat

We consider a model of competition between plasmid-bearing and plasmid-free organisms for two complementary nutrients in a chemostat. We assume that the plasmid-bearing organism produces an allelopathic agent at the cost of its reproductive abilities which is lethal to plasmid-free organism. Our analysis leads to different thresholds in terms of the model parameters acting as conditions under which the organisms associated with the system cannot thrive even in the absence of competition. Local stability of the system is obtained in the absence of one or both the organisms. Also, global stability of the system is obtained in the presence of both the organisms. Computer simulations have been carried out to illustrate various analytical results.     

Stoichiometric interpretation of Escherichia coli glucose catabolism under various oxygenation rates.

Metabolic by-product secretion is commonly observed in oxygen-limited cultures. Oxygen limitations occur because of limits in the capacity of the respiratory system or because of the oxygenation limits of the cultivation method used. The latter restriction is of considerable practical importance since it results in a critical cell concentration above which oxygenation is insufficient, leading to by-product secretion. In this study we used a flux balance approach to determine optimal metabolic performance of Escherichia coli under variable oxygen limitations. This method uses linear optimization to find optimal metabolic flux patterns with respect to cell growth. Cell growth was defined as precursor requirements on the basis of a composition analysis. A growth-associated maintenance requirement of 23 mmol of ATP per g of biomass and a non-growth-associated maintenance value of 5.87 mmol at ATP per g (dry weight)-h were incorporated on the basis of a comparison with experimental data. From computations of optimal growth increased oxygen limitations were found to result in the secretion of acetate, formate, and ethanol in that order. Consistent with the experimental data in the literature, by-product secretion rates increased linearly with the growth rate. The computed optimal growth under increasing oxygen limitation revealed four critical growth rates at which changes in the by-product secretion pattern were observed. Concomitant with by-product secretion under oxygen limitations were changes in metabolic pathway utilization. The shifts in metabolism were characterized by changes in the metabolic values (computed as shadow prices) of the various redox carriers. The redox potential was thus identified as a likely trigger that leads to metabolic shifts.2+ ă