Microbial growth on pollutant mixtures is an important aspect of bioremediation and wastewater treatment. However, efforts to develop mathematical models for mixed substrate kinetics have been limited. Nearly all models group either the microbial population (as "biomass") or the chemical species (e.g., as biological oxygen demand). When individual chemical species are considered, most models assume either no interaction or that the nature of the interaction is competition for the same rate-limiting enzyme. And when individual microbial species are considered, simple competition for the growth substrate is the only interaction included. Here, we present results using Pseudomonas putida F1 and Burkholderia sp. strain JS150 growing individually and together on benzene, toluene, phenol, and their mixtures and compare mathematical models to describe these results. We demonstrate that the simple models do not accurately predict the outcome of these biodegradation experiments, and we describe the development of a new model for substrate mixtures, the sum kinetics with interaction parameters (SKIP) model. In mixed-culture experiments, the interactions between species were substrate dependent and could not be predicted by simple competition models. Together, this set of experimental and modeling results presents our current state of work in this area and identifies challenges for future modeling efforts.

The biodegradation kinetics of toluene, phenol, and a mixture of toluene and phenol by Burkholderia species JS150 was measured and modeled. Both of these compounds can serve as the sole source of carbon and energy for this microorganism. The single‐substrate biodegradation kinetics was described well using the Monod model, with model constants of μmax,T = 0.39 h−1 and KS,T = 0.011 mM for growth on toluene and μmax,P = 0.309 h−1 and KS,P = 0.0054 mM for growth on phenol. Degradation of the mixture of toluene and phenol followed simultaneous utilization kinetics with toluene being the preferred substrate. Toluene was found to inhibit the rate of utilization of phenol while the presence of phenol had little effect on the rate of degradation of toluene. Of the kinetic models that were tested, one developed for microbial degradation of multiple substrates was able to describe substrate interactions and to model the mixture utilization by strain JS150. Simple competitive, noncompetitive, or uncompetitive substrate kinetics were not sufficient to describe the observed inhibitory interactions. © 2000 John Wiley & Sons, Inc. Biotechnol Bioeng 70: 428–435, 2000.

Interspecies interactions and changes in the rate and extent of biodegradation in mixed culture–mixed substrate studies were investigated. A binary mixed culture of Pseudomonas putida F1 and Burkholderia sp. JS150 degraded toluene, phenol, and their mixture. Both toluene and phenol can serve as sole sources of carbon and energy for both P. putida F1 and strain JS150. To investigate the population dynamics of this system, a fluorescent in‐situ hybridization method was chosen because of its ability to produce quantitative data, its low standard error, and the ease of use of this method. When the binary mixed culture was grown on toluene or phenol alone, significant interactions between the species were observed. These interactions could not be explained by a pure‐and‐simple competition model and were substrate dependent. Strain JS150 growth was slightly inhibited when grown with P. putida F1 on phenol, and P. putida F1 grew more rapidly than expected. Conversely, when the two species were grown together on toluene alone, P. putida F1 was inhibited while strain JS150 was unaffected. During growth of the mixed culture on a combination of toluene and phenol, the interactions were similar to that observed during growth on phenol alone; P. putida F1 growth was enhanced while strain JS150 was unaffected. Because of the observed interspecies interactions, monoculture kinetic parameters were not sufficient to describe the mixed culture kinetics in any experiment. This is one of the first reports of microbial population dynamics in which molecular microbial ecology and mathematical modeling have been combined. The use of the 16S‐rRNA‐based method allowed for observation and understanding of interspecies interactions that were not observable with standard culture‐based methods. These results suggest the need for more investigations that account for both substrate and microbial interactions when predicting the fate of organic pollutants in real systems. © 2000 John Wiley & Sons, Inc. Biotechnol Bioeng 70: 436–445, 2000.

Interspecies interactions and changes in the rate and extent of biodegradation in mixed culture-mixed substrate studies were investigated. A binary mixed culture of Pseudomonas putida F1 and Burkholderia sp. JS150 degraded toluene, phenol, and their mixture. Both toluene and phenol can serve as sole sources of carbon and energy for both P. putida F1 and strain JS150. To investigate the population dynamics of this system, a fluorescent in-situ hybridization method was chosen because of its ability to produce quantitative data, its low standard error, and the ease of use of this method. When the binary mixed culture was grown on toluene or phenol alone, significant interactions between the species were observed. These interactions could not be explained by a pure-and-simple competition model and were substrate dependent. Strain JS150 growth was slightly inhibited when grown with P. putida F1 on phenol, and P. putida F1 grew more rapidly than expected. Conversely, when the two species were grown together on toluene alone, P. putida F1 was inhibited while strain JS150 was unaffected. During growth of the mixed culture on a combination of toluene and phenol, the interactions were similar to that observed during growth on phenol alone; P. putida F1 growth was enhanced while strain JS150 was unaffected. Because of the observed interspecies interactions, monoculture kinetic parameters were not sufficient to describe the mixed culture kinetics in any experiment. This is one of the first reports of microbial population dynamics in which molecular microbial ecology and mathematical modeling have been combined. The use of the 16S-rRNA-based method allowed for observation and understanding of interspecies interactions that were not observable with standard culture-based methods. These results suggest the need for more investigations that account for both substrate and microbial interactions when predicting the fate of organic pollutants in real systems.