Metabolic modeling of spatial heterogeneity of biofilms in microbial fuel cells reveals substrate limitations in electrical current generation

Model development of bioelectrochemical systems: A critical review from the perspective of physiochemical principles and mathematical methods

Bioelectrochemical systems (BESs) are promising devices for wastewater treatment and bio-energy production. Since various processes are interacted and affect the overall performance of the device, the development of theoretical modeling is an efficient approach to understand the fundamental mechanisms that govern the performance of the BES. This review aims to summarize the physiochemical principle and mathematical method in BES models, which is of great importance for the establishment of an accurate model while has received little attention in previous reviews. In this review, we begin with a classification of existing models including bioelectrochemical models, electronic models, and machine learning models. Subsequently, physiochemical principles and mathematical methods in models are discussed from two aspects: one is the description of methodology how to build a framework for models, and the other is to further review additional methods that can enrich model functions. Finally, the advantages/disadvantages, extended applications, and perspectives of models are discussed. It is expected that this review can provide a viewpoint from methodologies to understand BES models.

A Catalytic Effectiveness Factor for a Microbial Electrolysis Cell Biofilm Model

The aim of this work is to propose a methodology to obtain an effectiveness factor for biofilm in a microbial electrolysis cell (MEC) system and use it to reduce a partial differential equation (PDE) biofilm MEC model to an ordinary differential equation (ODE) MEC model. The biofilm mass balances of the different species are considered. In addition, it is considered that all the involved microorganisms are attached to the anodic biological film. Three effectiveness factors are obtained from partial differential equations describing the spatial distributions of potential and substrate in the biofilm. Then, a model reduction is carried out using the global mass balances of the different species in the system. The reduced model with three uncertain but bounded effectiveness factors is evaluated numerically and analyzed in the sense of stability and parametric sensibility to demonstrate its applicability. The reduced ODE model is compared with a validated model taken from the literature, and the results are in good agreement. The biofilm effectiveness factor in MEC systems can be extended to the reduction of PDE models to obtain ODE models that are commonly used in optimization and control problems.

Bio-Electrochemical System Depollution Capabilities and Monitoring Applications: Models, Applicability, Advanced Bio-Based Concept for Predicting Pollutant Degradation and Microbial Growth Kinetics via Gene Regulation Modelling

Microbial fuel cells (MFC) are an emerging technology for waste, wastewater and polluted soil treatment. In this manuscript, pollutants that can be treated using MFC systems producing energy are presented. Furthermore, the applicability of MFC in environmental monitoring is described. Common microbial species used, release of genome sequences, and gene regulation mechanisms, are discussed. However, although scaling-up is the key to improving MFC systems, it is still a difficult challenge. Mathematical models for MFCs are used for their design, control and optimization. Such models representing the system are presented here. In such comprehensive models, microbial growth kinetic approaches are essential to designing and predicting a biosystem. The empirical and unstructured Monod and Monod-type models, which are traditionally used, are also described here. Understanding and modelling of the gene regulatory network could be a solution for enhancing knowledge and designing more efficient MFC processes, useful for scaling it up. An advanced bio-based modelling concept connecting gene regulation modelling of specific metabolic pathways to microbial growth kinetic models is presented here; it enables a more accurate prediction and estimation of substrate biodegradation, microbial growth kinetics, and necessary gene and enzyme expression. The gene and enzyme expression prediction can also be used in synthetic and systems biology for process optimization. Moreover, various MFC applications as a bioreactor and bioremediator, and in soil pollutant removal and monitoring, are explored.

Spatiotemporal metabolic modeling of bacterial life in complex habitats

The combination of genome-scale metabolic networks with spatially explicit representation of microbial habitats (spatiotemporal metabolic network modeling) paves the way to predict complex metabolic landscapes to a hitherto unparalleled detail, thus providing new insights into trophic interactions occurring at different scales. Placing detailed bacterial metabolism in realistic physical environment highlights the roles of physical barriers and diffusional bottlenecks on bacterial community interactions, structure and stability. We review recent advances in spatiotemporal metabolic network modeling using a few illustrative examples that highlight the immense potential of these novel approaches to interpret and design metabolic mediated interactions in structures (natural and engineered) environments.

Electron transfer pathways and kinetic analysis of cathodic simultaneous nitrification and denitrification process in microbial fuel cell system

Microbial fuel cell (MFC) is an innovative bioconversion technology for wastewater treatment accompanied with electricity recovery. In this study, a kinetic model was developed base on Activated Sludge Model No.1 (ASM1) to describe electron transfer pathways during the simultaneous nitrification and denitrification (SND) process in the biocathode system of a dual-chamber MFC. The batch running of the dual-chamber MFC system showed that it produced a power density up to 2.96 W m^-3 within 48 h, the achieved SND efficiency and autotrophic denitrification ratio in the cathodic denitrification process were up to 87.3 ± 0.8% and 69.5 ± 6.6%, respectively. Meanwhile, by integrating nitrification, autotrophic denitrification, heterotrophic denitrification, organic carbon oxidation, and oxygen reduction in the cathode, the model was able to precisely fit the concentration variations of NH₃-N, dissolved oxygen (DO) and chemical oxygen demand (COD) during the cathodic SND process (R² ≥ 0.9876). The cathode electrons tended to be completely utilized with the increase of autotrophic denitrification ratio in the cathodic denitrification process. When the nitrification rate was enhanced, the autotrophic denitrification would prevail in the competition with the heterotrophic denitrification. In summary, the developed model was confirmed to be effective and reliable for describing the electron transfer pathways and predicting the performance of the nitrogen removal reactions during the cathodic SND process in a double-chamber MFC.

A Review of Control-Oriented Bioelectrochemical Mathematical Models of Microbial Fuel Cells

A microbial fuel cell (MFC) is a potentially viable renewable energy option which promises effective and commercial harvesting of electrical power by bacterial movement and at the same time also treats wastewater. Microbial fuel cells are complicated devices and therefore research in this field needs interdisciplinary knowledge and involves diverse areas such as biological, chemical, electrical, etc. In recent decades, rapid strides have taken place in fuel cell research and this technology has become more efficient. For effective usage, such devices need advanced control techniques for maintaining a balance between substrate supply, mass, charge, and external load. Most of the research work in this area focuses on experimental work and have been described from the design perspective. Recently, the development in mathematical modeling of such cells has taken place which has provided a few mathematical models. Mathematical modeling provides a better understanding of the operations and the dynamics of MFCs, which will help to develop control and optimization strategies. Control-oriented bio-electrochemical models with mass and charge balance of MFCs facilitate the development of advanced nonlinear controllers. This work reviews the different mathematical models of such cells available in the literature and then presents suitable parametrization to develop control-oriented bio-electrochemical models of three different types of cells with their uncertain parameters.

Challenges of biofilm control and utilization: lessons from mathematical modelling

This article reviews modern applications of mathematical descriptions of biofilm formation. The focus is on theoretically obtained results which have implications for areas including the medical sector, food industry and wastewater treatment. Examples are given as to how models have contributed to the overall knowledge on biofilms and how they are used to predict biofilm behaviour. We conclude that the use of mathematical models of biofilms has demonstrated over the years the ability to significantly contribute to the vast field of biofilm research. Among other things, they have been used to test various hypotheses on the nature of interspecies interactions, viability of biofilm treatment methods or forces behind observed biofilm pattern formations. Mathematical models can also play a key role in future biofilm research. Many models nowadays are analysed through computer simulations and continue to improve along with computational capabilities. We predict that models will keep on providing answers to important challenges involving biofilm formation. However, further strengthening of the ties between various disciplines is necessary to fully use the tools of collective knowledge in tackling the biofilm phenomenon.

In Silico Metabolic Design of Two-Strain Biofilm Systems Predicts Enhanced Biomass Production and Biochemical Synthesis

Engineered biofilm consortia have the potential to solve important biotechnological problems that have proved difficult for monoculture biofilms and planktonic consortia, such as conversion of lignocellulosic material to useful biochemicals. While considerable experimental progress has been reported for engineering and characterizing biofilm consortia, the field still lacks in silico tools for simulation, design, and optimization of stable, robust, and productive designed consortia. We developed biofilm consortia metabolic models for two coculture systems centered around the ecological design motif of a primary cell type that utilizes a supplied electron donor and secretes acetate as a byproduct and a secondary cell type that consumes the acetate, relieving byproduct inhibition on the primary cell type and enhancing overall system biomass. The models presented in this paper predict that distinct metabolic niches for the two cell types could be established by supplying electron donors and acceptors at opposite ends of the biofilm and that acetate consumption by the secondary cell type could increase total biomass accumulation and the synthesis of valuable biochemicals, such as isobutanol, by the primary cell type. System tunability is enhanced when each cell type is supplied with a unique terminal electron acceptor at opposite ends of the biofilm rather than competing for a common electron acceptor. Our model provides good qualitative agreement with data for a synthetic Escherichia coli coculture system, suggesting that the proposed design rules may have wide applicability to engineered biofilm consortia.

Competitive resource allocation to metabolic pathways contributes to overflow metabolisms and emergent properties in cross-feeding microbial consortia

Resource scarcity is a common stress in nature and has a major impact on microbial physiology. This review highlights microbial acclimations to resource scarcity, focusing on resource investment strategies for chemoheterotrophs from the molecular level to the pathway level. Competitive resource allocation strategies often lead to a phenotype known as overflow metabolism; the resulting overflow byproducts can stabilize cooperative interactions in microbial communities and can lead to cross-feeding consortia. These consortia can exhibit emergent properties such as enhanced resource usage and biomass productivity. The literature distilled here draws parallels between in silico and laboratory studies and ties them together with ecological theories to better understand microbial stress responses and mutualistic consortia functioning.

A Genome-Scale Modeling Approach to Quantify Biofilm Component Growth of Salmonella Typhimurium

Salmonella typhimurium (S. typhimurium) is an extremely dangerous foodborne bacterium that infects both animal and human subjects, causing fatal diseases around the world. Salmonella's robust virulence, antibiotic-resistant nature, and capacity to survive under harsh conditions are largely due to its ability to form resilient biofilms. Multiple genome-scale metabolic models have been developed to study the complex and diverse nature of this organism's metabolism; however, none of these models fully integrated the reactions and mechanisms required to study the influence of biofilm formation. This work developed a systems-level approach to study the adjustment of intracellular metabolism of S. typhimurium during biofilm formation. The most advanced metabolic reconstruction currently available, STM_v1.0, was 1st extended to include the formation of the extracellular biofilm matrix. Flux balance analysis was then employed to study the influence of biofilm formation on cellular growth rate and the production rates of biofilm components. With biofilm formation present, biomass growth was examined under nutrient rich and nutrient deficient conditions, resulting in overall growth rates of 0.8675 and 0.6238 h^-1 respectively. Investigation of intracellular flux variation during biofilm formation resulted in the elucidation of 32 crucial reactions, and associated genes, whose fluxes most significantly adapt during the physiological response. Experimental data were found in the literature to validate the importance of these genes for the biofilm formation of S. typhimurium. This preliminary investigation on the adjustment of intracellular metabolism of S. typhimurium during biofilm formation will serve as a platform to generate hypotheses for further experimental study on the biofilm formation of this virulent bacterium.

Metabolic modeling of a chronic wound biofilm consortium predicts spatial partitioning of bacterial species

Background

Chronic wounds are often colonized by consortia comprised of different bacterial species growing as biofilms on a complex mixture of wound exudate. Bacteria growing in biofilms exhibit phenotypes distinct from planktonic growth, often rendering the application of antibacterial compounds ineffective. Computational modeling represents a complementary tool to experimentation for generating fundamental knowledge and developing more effective treatment strategies for chronic wound biofilm consortia.

Results

We developed spatiotemporal models to investigate the multispecies metabolism of a biofilm consortium comprised of two common chronic wound isolates: the aerobe Pseudomonas aeruginosa and the facultative anaerobe Staphylococcus aureus. By combining genome-scale metabolic reconstructions with partial differential equations for metabolite diffusion, the models were able to provide both temporal and spatial predictions with genome-scale resolution. The models were used to analyze the metabolic differences between single species and two species biofilms and to demonstrate the tendency of the two bacteria to spatially partition in the multispecies biofilm as observed experimentally. Nutrient gradients imposed by supplying glucose at the bottom and oxygen at the top of the biofilm induced spatial partitioning of the two species, with S. aureus most concentrated in the anaerobic region and P. aeruginosa present only in the aerobic region. The two species system was predicted to support a maximum biofilm thickness much greater than P. aeruginosa alone but slightly less than S. aureus alone, suggesting an antagonistic metabolic effect of P. aeruginosa on S. aureus. When each species was allowed to enhance its growth through consumption of secreted metabolic byproducts assuming identical uptake kinetics, the competitiveness of P. aeruginosa was further reduced due primarily to the more efficient lactate metabolism of S. aureus. Lysis of S. aureus by a small molecule inhibitor secreted from P. aeruginosa and/or P. aeruginosa aerotaxis were predicted to substantially increase P. aeruginosa competitiveness in the aerobic region, consistent with in vitro experimental studies.

Conclusions

Our biofilm modeling approach allows the prediction of individual species metabolism and interspecies interactions in both time and space with genome-scale resolution. This study yielded new insights into the multispecies metabolism of a chronic wound biofilm, in particular metabolic factors that may lead to spatial partitioning of the two bacterial species. We believe that P. aeruginosa lysis of S. aureus combined with nutrient competition is a particularly relevant scenario for which model predictions could be tested experimentally.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0334-8) contains supplementary material, which is available to authorized users.

Carbon nanotube modification of microbial fuel cell electrodes

The use of carbon nanotubes (CNTs) for energy harvesting devices is preferable due to their unique mechanical, thermal, and electrical properties. On the other hand, microbial fuel cells (MFCs) are promising devices to recover carbon-neutral energy from the organic matters, and have been hindered with major setbacks towards commercialization. Nanoengineered CNT-based materials show remarkable electrochemical properties, and therefore have provided routes towards highly effective modification of MFC compartments to ultimately reach the theoretical limits of biomass energy recovery, low-cost power production, and thus the commercialization of MFCs. Moreover, these CNT-based composites offer significant flexibility in the design of MFCs that enable their use for a broad spectrum of applications ranging from scaled-up power generation to medically related devices. This article reviews the recent advances in the modification of MFCs using CNTs and CNT-based composites, and the extent to which each modification route impacts MFC power and current generation.

Spatiotemporal modeling of microbial metabolism

BACKGROUND

Microbial systems in which the extracellular environment varies both spatially and temporally are very common in nature and in engineering applications. While the use of genome-scale metabolic reconstructions for steady-state flux balance analysis (FBA) and extensions for dynamic FBA are common, the development of spatiotemporal metabolic models has received little attention.

RESULTS

We present a general methodology for spatiotemporal metabolic modeling based on combining genome-scale reconstructions with fundamental transport equations that govern the relevant convective and/or diffusional processes in time and spatially varying environments. Our solution procedure involves spatial discretization of the partial differential equation model followed by numerical integration of the resulting system of ordinary differential equations with embedded linear programs using DFBAlab, a MATLAB code that performs reliable and efficient dynamic FBA simulations. We demonstrate our methodology by solving spatiotemporal metabolic models for two systems of considerable practical interest: (1) a bubble column reactor with the syngas fermenting bacterium Clostridium ljungdahlii; and (2) a chronic wound biofilm with the human pathogen Pseudomonas aeruginosa. Despite the complexity of the discretized models which consist of 900 ODEs/600 LPs and 250 ODEs/250 LPs, respectively, we show that the proposed computational framework allows efficient and robust model solution.

CONCLUSIONS

Our study establishes a new paradigm for formulating and solving genome-scale metabolic models with both time and spatial variations and has wide applicability to natural and engineered microbial systems.

Genome-scale modelling of microbial metabolism with temporal and spatial resolution

Most natural microbial systems have evolved to function in environments with temporal and spatial variations. A major limitation to understanding such complex systems is the lack of mathematical modelling frameworks that connect the genomes of individual species and temporal and spatial variations in the environment to system behaviour. The goal of this review is to introduce the emerging field of spatiotemporal metabolic modelling based on genome-scale reconstructions of microbial metabolism. The extension of flux balance analysis (FBA) to account for both temporal and spatial variations in the environment is termed spatiotemporal FBA (SFBA). Following a brief overview of FBA and its established dynamic extension, the SFBA problem is introduced and recent progress is described. Three case studies are reviewed to illustrate the current state-of-the-art and possible future research directions are outlined. The author posits that SFBA is the next frontier for microbial metabolic modelling and a rapid increase in methods development and system applications is anticipated.

Metabolic modeling of synthesis gas fermentation in bubble column reactors

BACKGROUND

A promising route to renewable liquid fuels and chemicals is the fermentation of synthesis gas (syngas) streams to synthesize desired products such as ethanol and 2,3-butanediol. While commercial development of syngas fermentation technology is underway, an unmet need is the development of integrated metabolic and transport models for industrially relevant syngas bubble column reactors.

RESULTS

We developed and evaluated a spatiotemporal metabolic model for bubble column reactors with the syngas fermenting bacterium Clostridium ljungdahlii as the microbial catalyst. Our modeling approach involved combining a genome-scale reconstruction of C. ljungdahlii metabolism with multiphase transport equations that govern convective and dispersive processes within the spatially varying column. The reactor model was spatially discretized to yield a large set of ordinary differential equations (ODEs) in time with embedded linear programs (LPs) and solved using the MATLAB based code DFBAlab. Simulations were performed to analyze the effects of important process and cellular parameters on key measures of reactor performance including ethanol titer, ethanol-to-acetate ratio, and CO and H2 conversions.

CONCLUSIONS

Our computational study demonstrated that mathematical modeling provides a complementary tool to experimentation for understanding, predicting, and optimizing syngas fermentation reactors. These model predictions could guide future cellular and process engineering efforts aimed at alleviating bottlenecks to biochemical production in syngas bubble column reactors.