Population balance modeling is undergoing phenomenal growth in its applications, and this growth is accompanied by multifarious reviews. This review aims to fortify the model's fundamental base, as well as point to a variety of new applications, including modeling of crystal morphology, cell growth and differentiation, gene regulatory processes, and transfer of drug resistance. This is accomplished by presenting the many faces of population balance equations that arise in the foregoing applications.

Flux balance analysis (FBA) in combination with the decomposition of metabolic networks into elementary modes has provided a route to modeling cellular metabolism. It is dependent, however, on the availability of external fluxes such as substrate uptake or growth rate before estimates can become available of intracellular fluxes. The framework classically does not allow modeling of metabolic regulation or the formulation of dynamic models except through dynamic measurement of external fluxes. The cybernetic modeling approach of Ramkrishna and coworkers provides a dynamic framework for modeling metabolic systems because of its focus on describing regulatory processes based on cybernetic arguments and hence has the capacity to describe both external and internal fluxes. In this article, we explore the alternative of developing hybrid models combining cybernetic models for the external fluxes with the flux balance approach for estimation of the internal fluxes. The approach has the merit of the simplicity of the early cybernetic models and hence computationally facile while also providing detailed information on intracellular fluxes. The hybrid model of this article is based on elementary mode decomposition of the metabolic network. The uptake rates for the various elementary modes are combined using global cybernetic variables based on maximizing substrate uptake rates. Estimation of intracellular metabolism is based on its stoichiometric coupling with the external fluxes under the assumption of (pseudo‐) steady state conditions. The set of parameters of the hybrid model was estimated with the aid of nonlinear optimization routine, by fitting simulations with dynamic experimental data on concentrations of biomass, substrate, and fermentation products. The hybrid model estimations were tested with FBA (based on measured substrate uptake rate) for two different metabolic networks (one is a reduced network which fixes ATP contribution to the biomass and maintenance requirement of ATP, and the other network is a more complex network which has a separate reaction for maintenance.) for the same experiment involving anaerobic growth of E. coli GJT001. The hybrid model estimated glucose consumption and all fermentation byproducts to better than 10%. The FBA makes similar estimations of fermentation products, however, with the exception of succinate. The simulation results show that the global cybernetic variables alone can regulate the metabolic reactions obtaining a very satisfactory fit to the measured fermentation byproducts. In view of the hybrid model's ability to predict biomass growth and fermentation byproducts of anaerobic E. coli GJT001, this reduced order model offers a computationally efficient alternative to more detailed models of metabolism and hence useful for the simulation of bioreactors.

In recent times, stochastic treatments of gene regulatory processes have appeared in the literature in which a cell exposed to a signaling molecule in its environment triggers the synthesis of a specific protein through a network of intracellular reactions. The stochastic nature of this process leads to a distribution of protein levels in a population of cells as determined by a Fokker-Planck equation. Often instability occurs as a consequence of two (stable) steady state protein levels, one at the low end representing the “off” state, and the other at the high end representing the “on” state for a given concentration of the signaling molecule within a suitable range. A consequence of such bistability has been the appearance of bimodal distributions indicating two different populations, one in the “off” state and the other in the “on” state. The bimodal distribution can come about from stochastic analysis of a single cell. However, the concerted action of the population altering the extracellular concentration in the environment of individual cells and hence their behavior can only be accomplished by an appropriate population balance model which accounts for the reciprocal effects of interaction between the population and its environment. In this study, we show how to formulate a population balance model in which stochastic gene expression in individual cells is incorporated. Interestingly, the simulation of the model shows that bistability is neither sufficient nor necessary for bimodal distributions in a population. The original notion of linking bistability with bimodal distribution from single cell stochastic model is therefore only a special consequence of a population balance model. Traditionally cells in a population have been assumed to behave identically by using deterministic mathematical equations describing average cell behavior, thus ignoring its inherent randomness. A single cell stochastic model has therefore evolved in the literature to overcome this drawback. However, this single cell perspective does not account for interaction between the cell population and its environment. Since stochastic behavior leads to each cell acting differently, the cumulative impact of individual cells on their environment and consequent influence of the latter on each cell could constitute a behavior at variance. Thus in nature, cells are constantly under the influence of a highly dynamic environment which in turn is influenced by the dynamics of the cell population. A typical single cell stochastic model ignores such an interaction between the population and its environment, and uses probability distribution of a single cell to represent the entire population, which may lead to inappropriate predictions. In this study, we propose a population balance model coupled with stochastic gene regulation to demonstrate the behavior of a population in which its interactive behavior with its environment is considered. Our simulation results show that bistability is neither sufficient nor necessary for bimodal distributions in a population.