Modeling of an aerobic biofilm reactor with double-limiting substrate kinetics: bifurcational and dynamical analysis

A mathematical model of an aerobic biofilm reactor is presented to investigate the bifurcational patterns and the dynamical behavior of the reactor as a function of different key operating parameters. Suspended cells and biofilm are assumed to grow according to double limiting kinetics with phenol inhibition (carbon source) and oxygen limitation. The model presented by Russo et al. is extended to embody key features of the phenomenology of the granular‐supported biofilm: biofilm growth and detachment, gas–liquid oxygen transport, phenol, and oxygen uptake by both suspended and immobilized cells, and substrate diffusion into the biofilm. Steady‐state conditions and stability, and local dynamic behavior have been characterized. The multiplicity of steady states and their stability depend on key operating parameter values (dilution rate, gas–liquid mass transfer coefficient, biofilm detachment rate, and inlet substrate concentration). Small changes in the operating conditions may be coupled with a drastic change of the steady‐state scenario with transcritical and saddle‐node bifurcations. The relevance of concentration profiles establishing within the biofilm is also addressed. When the oxygen level in the liquid phase is <10% of the saturation level, the biofilm undergoes oxygen starvation and the active biofilm fraction becomes independent of the dilution rate. © 2011 American Institute of Chemical Engineers Biotechnol. Prog., 2011     

Dynamic effects related to steady-state multiplicity in continuous Saccharomyces cerevisiae cultivations

The behavioral differences between chemostat and productostat cultivation of aerobic glucose-limited Saccharomyces cerevisiae were investigated. Three types of experiments were conducted: a chemostat, where the dilution rate was shifted up or down in stepwise manner; and a productostat, with either stepwise changed or a rampwise increased ethanol setpoint, i.e., an accelero-productostat. The transient responses from chemostat and productostat experiments were interpreted using a simple metabolic flux model. In a productostat it was possible to obtain oxido-reductive steady states at dilution rates far below Dcrit due to a strong repression of the respiratory system. However, these steady states could not be obtained in a chemostat, since a dilution rate shift-down from an oxido-reductive steady state led to a derepression of the respiratory system. It can therefore be concluded that the range of dilution rates where steady-state multiplicity can be obtained differs depending on the operation mode and that this dilution rate multiplicity range may appear larger in a productostat than in a chemostat. A more narrow multiplicity range, however, was obtained when the productostat was operated as an accelero-productostat.     

Transient behavior of a continuous stirred tank biological reactor utilizing phenol as an inhibitory substrate

Transient experiments were conducted on a Pseudomomas utilizing phenol in a continuous culture by disturbing the influent substrate concentration and dilution rate. Two stable steady states existed for some ranges of the parameters. Highly damped oscillations were observed in approaching a new high conversion steady state or in returning to a new high conversion steady state following a small disturbance. When a large disturbance was applied there was a smooth (overdamped) approach to a new low conversion steady state. The observed oscillatory behavior for small disturbances was predicted by a modified Powell-Ierusalemskii bottleneck model, but could not be predicted by a Monod-Haldane model; neither model was accurate for predicting the effect of large disturbances. A constant wall growth factor was used to account for microbial film activity, and the existence of two stable states was directly due to the presence of the film.     

Multiplicity and stability analysis of microorganisms in continuous culture: effects of metabolic overflow and growth inhibition

Metabolic overflow (enhanced uptake of substrate and secretion of intermediates) is a phenomenon often observed for cells grown under substrate excess. Growth inhibition by substrate and/or product is also normally found for this kind of culture. An effort is made in this work to analyze the dynamic behavior of a continuous culture subject to metabolic overflow and growth inhibition by substrate and/or product. Analysis of a model system shows that in a certain range of operating conditions three nonwashout steady state solutions are possible. Local stability analysis indicates that only two of them are stable thus leading to multiplicity and hysteresis. Further analysis of the intrinsic effects of different terms describing the metabolic overflow and growth inhibitions reveals that for the model system and the parameters considered, the combined effects of product inhibition and an enhanced formation rate of product under substrate excess cause the multiplicity and hysteresis. Growth inhibition by substrate and/or an enhanced substrate uptake appear not to be necessary conditions. The combined effects of enhanced product formation and product inhibition can also lead to unusual dynamic behavior such as a prolonged time period to reach a steady state, oscillatory transition from one steady state to another, and sustained oscillations. Using the occurrence of multiplicity and oscillation as criteria, the operating regime of a continuous culture can be divided into four domains: one with multiplicity and oscillation, one with unique steady state but possible oscillatory behavior, the other two with unique and stable steady state. The model predictions are in accordance with recent experimental results. The results presented in this work may be used as guidelines for choosing proper operating conditions of similar culture systems to avoid undesired instability and multiplicity. Copyright 1998 John Wiley & Sons, Inc.     

Cybernetic modeling and regulation of metabolic pathways in multiple steady states of hybridoma cells

Hybridoma cells utilize a pair of complementary and partially substitutable substrates, glucose and glutamine, for growth. It has been shown that cellular metabolism shifts under different culture conditions. When those cultures at different metabolic states are switched to a continuous mode, they reach different steady states under the same operating conditions. A cybernetic model was constructed to describe the complementary and partial substitutable nature of substrate utilization. The model successfully predicted the metabolic shift and multiple steady-state behavior. The results are consistent with the experimental observation that the history of the culture affects the resulting steady state.     

Multiple stable states and hysteresis in continuous, oscillating cultures of budding yeast

The conditions that precede the onset of autonomous oscillations in continuous yeast cultures were studied in three different types of experiments. It was found that the final state of the culture depended on the protocol used to start up the reactor. Batch cultures, switched to continuous operation at different stages of the batch growth curve, all exhibited similar dynamics-ethanol depletion followed by autonomous oscillations. Small perturbations of the distribution of states in the reactor, achieved by addition of externally grown cells, were able to quench the oscillatory dynamics. Reaching the desired operating point by slow dilution rate changes gave rise to different final states, two oscillatory states and one steady state, depending on the rate of change in dilution rate. The multiplicity of stable states at a single operating point is not explained by any current distributed model and points toward a segregated mechanism of these oscillations.     

Chemostat-cultivated Escherichia coli at high dilution rate: multiple steady states and drift

The representation of metabolic network reaction kinetics in a scaled, polynomial form can allow for the prediction of multiple steady states. The polynomial formalism is used to study chemostat-cultured Escherichia coli which has been observed to exhibit two multiple steady states under ammonium ion-limited growth conditions: a high cell density-low ammonium ion concentration steady state and a low cell density-high ammonium ion concentration steady state. Additionally, the low-cell-density steady state has been observed to drift to the high-cell-density steady state. Inspection of the steady-state rate expressions for the ammonium ion transport/assimilation network (in polynomial form) suggests that at low ammonium ion concentrations, two steady states are possible. One corresponds to heavy use of the glutamine synthetase-glutamate synthase (GLNS-GS) branch and the second to heavy use of the glutamate dehydrogenase (GDH) branch. Realization of the predicted intracellular steady states is also found to be dependent on the parameters of the transport process. Moreover, the two steady states differ in where their energy intensity lies. To explain the drift, GLNS, which is inducible under low ammonium ion concentrations, is suggested to be a "memory element." A chemostat-based model is developed to illustrate that perturbations in dilution rate can lead to drift between the two steady states provided that the disturbance in dilution rate is sufficiently large and/or long in duration.     

Metabolic flux analysis of hybridoma continuous culture steady state multiplicity.

Physiological state multiplicity was observed in continuous cultures of the hybridoma cell line ATCC CRL-1606 cultivated in glutamine-limited steady state chemostats. At the same dilution rate (0.04 h-1), two physiologically different cultures were obtained which exhibited similar growth rates and viabilities but drastically different cell concentrations (7.36 x 10(5) and 1.36 x 10(6) cells/mL). Metabolic flux analysis conducted using metabolite and gas exchange rate measurements revealed a more efficient culture for the steady state with the higher cell concentration, as measured by the fraction of pyruvate carbon flux shuttled into the TCA cycle for energy generation. The low-efficiency steady state was achieved after innoculation by growing the cells in a nutrient rich environment, first in batch mode followed by a stepwise increase of the dilution rate to its set point at 0.04 h-1. The high-efficiency steady state was achieved by reducing the dilution rate to progressively lower values to 0.01 h-1 resulting in conditions of stricter nutrient limitation. The high energetic efficiency attained under such conditions was preserved upon increasing the chemostat dilution rate back to 0.04 h-1 with a higher nutrient consumption, resulting in approximate doubling of the steady state cell concentration. This metabolic adaptation is unlikely due to favorable genetic mutations and could be implemented for improving cell culture performance by inducing cellular metabolic shifts to more efficient flux distribution patterns.     

Metabolic modeling of Saccharomyces cerevisiae using the optimal control of homeostasis: a cybernetic model definition

A model is presented to describe the observed behavior of microorganisms that aim at metabolic homeostasis while growing and adapting to their environment in an optimal way. The cellular metabolism is seen as a network with a multiple controller system with both feedback and feedforward control, i.e., a model based on a dynamic optimal metabolic control. The dynamic network consists of aggregated pathways, each having a control setpoint for the metabolic states at a given growth rate. This set of strategies of the cell forms a true cybernetic model with a minimal number of assumptions. The cellular strategies and constraints were derived from metabolic flux analysis using an identified, biochemically relevant, stoichiometry matrix derived from experimental data on the cellular composition of continuous cultures of Saccharomyces cerevisiae. Based on these data a cybernetic model was developed to study its dynamic behavior. The growth rate of the cell is determined by the structural compounds and fluxes of compounds related to central metabolism. In contrast to many other cybernetic models, the minimal model does not consist of any assumed internal kinetic parameters or interactions. This necessitates the use of a stepwise integration with an optimization of the fluxes at every time interval. Some examples of the behavior of this model are given with respect to steady states and pulse responses. This model is very suitable for describing semiquantitatively dynamics of global cellular metabolism and may form a useful framework for including structured and more detailed kinetic models.     

Unveiling steady-state multiplicity in hybridoma cultures: the cybernetic approach

Mammalian cells grown in suspension produce waste metabolites such as lactate, alanine, and ammonia, which reduce the yield of cell mass and the desired product on the nutrients supplied. Previous studies (Cruz et al., 1999; Europa et al., 2000; Follstad et al., 1999) have shown that the cells can be made to alter their metabolism by starving them on their nutrients in continuous cultures at low dilution rates or starting the culture as a fed-batch. This leads to multiple steady states in continuous reactors, with some states being more favorable than others. Mathematical models that take into account the metabolic regulation that leads to these multiple steady states are invaluable tools for bioreactor control. In this article we present a cybernetic modeling strategy in which Metabolic Flux Analysis (MFA) is used to guide the cybernetic formulation. The hybridoma model presented as a result of this strategy considers the partially substitutable, partially complementary nature of glucose and glutamine. The choice of competitions within the network is guided by MFA and the model is successful in explaining the three multiple steady states observed. The cybernetic model though identified for the hybridoma experiments of Hu and others (Europa et al., 2000) seem generally applicable to mammalian systems as it captures the pathways that are common to mammalian cells grown in suspension. The model presented here could be used for start-up strategies for continuous reactors and model-based feedback control for maintaining high productivity of the reactor.     

Continuous cultivation of Saccharomyces cerevisiae. II. Growth on ethanol under transient-state conditions

Growth of Saccharomyces cerevisiae LBG H 1022 on ethanol under transient-state conditions was studied. As a cultivation device, an aerated Chemap fermentor combined with continuously working gas analyzers for oxygen and carbon dioxide was used. Yeast cell dry matter, substrate concentration, specific oxygen uptake, specific carbon dioxide release, and respiration quotient were measured during the different transient states. Depending on which range of the dilution rate the initial steady state was found, we obtain different responses to the shift experiment. For the lower range, up to D = 0.07, we deal with damped oscillations ranging above and below the steady-state values. For the higher specific growth rates, the rate of damping is strongly enhanced and the shape of the curves becomes an asymptotic approach to the final steady states.     

Dynamics and stability analysis of the growth and astaxanthin production system of Haematococcus pluvialis

This paper investigated high cell density cultivation of Haematococcus pluvialis for astaxanthin production in 3.7-L bioreactors. A biomass concentration of 2.74 g L⁻¹and an astaxanthin yield of 64.4 mg L⁻¹ were obtained. Based on the experimental results, a new and simple dynamic model is proposed, differing from Monod kinetics, to describe cell growth, product formation and substrate consumption. Good agreement was found between the model predictions and experimental data. The model revealed that there was cell growth inhibition on product formation and product feedback compensation for substrate consumption, but no substrate inhibition or product inhibition of cell growth. Stability analysis demonstrated that no multiplicity of steady states was observed; the unique positive steady state was locally asymptotically stable; and the effect of dilution rate on steady states was greater than that of the initial substrate concentration. Received 23 February 1999/ Accepted in revised form 08 June 1999     

A biochemically structured model for Saccharomyces cerevisiae.

A biochemically structured model for the aerobic growth of Saccharomyces cerevisiae on glucose and ethanol is presented. The model focuses on the pyruvate and acetaldehyde branch points where overflow metabolism occurs when the growth changes from oxidative to oxido-reductive. The model is designed to describe the onset of aerobic alcoholic fermentation during steady-state as well as under dynamical conditions, by triggering an increase in the glycolytic flux using a key signalling component which is assumed to be closely related to acetaldehyde. An investigation of the modelled process dynamics in a continuous cultivation revealed multiple steady states in a region of dilution rates around the transition between oxidative and oxido-reductive growth. A bifurcation analysis using the two external variables, the dilution rate, D, and the inlet concentration of glucose, S(f), as parameters, showed that a fold bifurcation occurs close to the critical dilution rate resulting in multiple steady-states. The region of dilution rates within which multiple steady states may occur depends strongly on the substrate feed concentration. Consequently a single steady state may prevail at low feed concentrations, whereas multiple steady states may occur over a relatively wide range of dilution rates at higher feed concentrations.     

Nonlinear dynamics of eucaryotic pyruvate dehydrogenase multienzyme complex: decarboxylation rate,oscillations, and multiplicity

Pyruvate conversion to acetyl-CoA by the pyruvate dehydrogenase (PDH) multienzyme complex is known as a key node in affecting the metabolic fluxes of animal cell culture. However, its possible role in causing possible nonlinear dynamic behavior such as oscillations and multiplicity of animal cells has received little attention. In this work, the kinetic and dynamic behavior of PDH of eucaryotic cells has been analyzed by using both in vitro and simplified in vivo models. With the in vitro model the overall reaction rate (nu(1)) of PDH is shown to be a nonlinear function of pyruvate concentration, leading to oscillations under certain conditions. All enzyme components affect nu(1) and the nonlinearity of PDH significantly, the protein X and the core enzyme dihydrolipoamide acyltransferase (E2) being mostly predominant. By considering the synthesis rates of pyruvate and PDH components the in vitro model is expanded to emulate in vivo conditions. Analysis using the in vivo model reveals another interesting kinetic feature of the PDH system, namely, multiple steady states. Depending on the pyruvate and enzyme levels or the operation mode, either a steady state with high pyruvate decarboxylation rate or a steady state with significantly lower decarboxylation rate can be achieved under otherwise identical conditions. In general, the more efficient steady state is associated with a lower pyruvate concentration. A possible time delay in the substrate supply and enzyme synthesis can also affect the steady state to be achieved and leads to oscillations under certain conditions. Overall, the predictions of multiplicity for the PDH system agree qualitatively well with recent experimental observations in animal cell cultures. The model analysis gives some hints for improving pyruvate metabolism in animal cell culture.     

Prediction of optimal biofilm thickness for membrane-attached biofilms growing in an extractive membrane bioreactor

This article presents a mathematical model of membrane-attached biofilm (MAB) behavior in a single-tube extractive membrane bioreactor (STEMB). MABs can be used for treatment of wastewaters containing VOCs, treatment of saline wastewaters, and nitrification processes. Extractive membrane bioreactors (EMBs) are employed to prevent the direct contact between a toxic volatile pollutant and the aerated gas by allowing counterdiffusion of substrates; i.e., pollutant diffuses from the tube side into the biofilm, whereas oxygen diffuses from the shell side into the biofilm. This reduces the air stripping problems usually found in conventional bioreactors. In this study, the biodegradation of a toxic VOC (1,2-dichloroethane, DCE) present in a synthetic wastewater has been investigated. An unstructured model is used to describe cell growth and cell decay in the MAB. The model has been verified by comparing model predicted trends with experimental data collected over 5 to 20-day periods, and has subsequently been used to model steady states in biofilm behavior over longer time scales. The model is capable of predicting the correct trends in system variables such as biofilm thickness, DCE flux across the membrane, carbon dioxide evolution, and suspended biomass. Steady states (constant biofilm thickness and DCE flux) are predicted, and factors that affect these steady states, i.e., cell endogeneous decay rate, and biofilm attrition, are investigated. Biofilm attrition does not have a great influence on biofilm behavior at low values of detachment coefficient close to those typically reported in the literature. Steady-state biofilm thickness is found to be an important variable; a thin biofilm results in a high DCE flux across the membrane, but with the penalty of a high loss of DCE via air stripping. The optimal biofilm thickness at steady state can be determined by trading off the decrease in air stripping (desirable) and the decrease in DCE flux (undesirable) which occur simultaneously as the thickness increases. (c) 1996 John Wiley & Sons, Inc.     

Multiplicity and stability analysis of continuous pre-fermentation of cheese culture

Unlike the usual multiplicity of steady states in continuous stirred-tank bioreactors caused by the non-monotonicity of substrate dependent rate expression(s), this paper described a different multiplicity caused by the effect of [H⁺]. The multiplicity of the bioreactor for cheese pre-fermentation is also described in this paper. Phase plane analysis shows clearly the separatrix that separates the zone of stable steady state from that of unstable steady state. The steady state becomes unique and elimination of multiplicity by supplying a seed stream, which also guarantees stability. Stability analysis of the multiple steady state shows that unlike the usual necessary condition (∂μ/∂S>0), this process stability demands the opposite necessary condition (∂μ/∂[H⁺]<0).     

Inhibition Kinetics of Phenol Degradation by Pseudomonas aeruginosa from Continuous Culture and Washout Data

Steady states of a continuous culture with an inhibitory substrate were used to estimate kinetic parameters under substrate limitation (chemostat operation). Pure cultures of an indigenous Pseudomonas aeruginosa were grown in continuous culture on phenol, the sole source of carbon and energy, at dilution rates of 0.010 to 0.20 h^{- 1}. Using different dilution rates, several steady states were investigated and the specific phenol consumption rates were calculated. In addition, phenol degradation was investigated by increasing the dilution rate above the critical dilution rate (washout cultivation). The results showed that the specific phenol consumption rate increased with increased dilution rate at steady state and that the degradation by Pseudomonas aeruginosa can be described by simple substrate inhibition kinetics under substrate limitation but cannot be described by simple substrate inhibition kinetics under washout cultivation. Fitting of the steady-state data from continuous cultivation to various inhibition models resulted in the best fit for the Yano and Koga kinetic inhibition model. The r_{s max} value of 0.278 mg/mg/h obtained from the Yano and Koga equation was comparable to the experimentally calculated r_{s max} value of 0.283 mg/mg/h obtained under washout cultivation.     

Influence of microbial growth kinetics on steady state multiplicity and stability of a two-step nitrification (SHARON) model

In this paper, the influence of microbial growth kinetics on the number and the stability of steady states for a nitrogen removal process is addressed. A two-step nitrification model is studied, in which the maximum growth rate of ammonium oxidizers is larger than the one of nitrite oxidizers. This model describes the behavior of a SHARON reactor for the treatment of wastewater streams with high ammonium concentrations. Steady states are identified through direct calculation using a canonical state space model representation, for several types of microbial kinetics. The stability of the steady states is assessed and the corresponding phase portraits are analyzed. Practical operation of a SHARON reactor aims at reaching ammonium conversion to nitrite while suppressing further conversion to nitrate. Regions in the input space are identified that result in this desired behavior, with only nitrite formation. It is demonstrated that not only the dilution rate plays a role, as is commonly known, but also the influent ammonium concentration. Besides, the type of microbial (inhibition) kinetics has a nonnegligible influence. While the results indicate that product inhibition does not affect the number of steady states of a (bio)reactor model, it is shown that substrate inhibition clearly yields additional steady states. Particular attention is devoted to the physical interpretation of these phenomena.     

Experimental investigations of multiple steady states in aerobic continuous cultivations of Saccharomyces cerevisiae

The steady-state behavior of a glucose-limited, aerobic, continuous cultivation of Saccharomyces cerevisiae CEN.PK113-7D was investigated around the critical dilution rate. Oxido-reductive steady states were obtained at dilution rates up to 0.09 h(-1) lower than the critical dilution rate by operating the bioreactor as a productostat, where the dilution rate was controlled on the basis of an ethanol measurement. Thus, the experimental investigations revealed that multiple steady states exist in a region of dilution rates below the critical dilution rate. The existence of multiple steady states was attributed to two distinct physiological effects occurring when growth changed from oxidative to oxido-reductive: (i) a decrease in the efficiency of ATP production and utilization (at ethanol concentrations below 3 g/L) and (ii) repression of the oxidative metabolism (at higher ethanol concentrations). The first effect was best observed at low ethanol concentrations, where multiple steady states were observed even when no repression of the oxidative metabolism was evident, i.e., the oxidative capacity was constant. However, at higher ethanol concentrations repression of the oxidative metabolism was observed (the oxidative capacity decreased), and this resulted in a broader range of dilution rates where multiple steady states could be found.     

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.