Augmenting Biogas Process Modeling by Resolving Intracellular Metabolic Activity

The process of anaerobic digestion in which waste biomass is transformed to methane by complex microbial communities has been modeled for more than 16 years by parametric gray box approaches that simplify process biology and do not resolve intracellular microbial activity. Information on such activity, however, has become available in unprecedented detail by recent experimental advances in metatranscriptomics and metaproteomics. The inclusion of such data could lead to more powerful process models of anaerobic digestion that more faithfully represent the activity of microbial communities. We augmented the Anaerobic Digestion Model No. 1 (ADM1) as the standard kinetic model of anaerobic digestion by coupling it to Flux-Balance-Analysis (FBA) models of methanogenic species. Steady-state results of coupled models are comparable to standard ADM1 simulations if the energy demand for non-growth associated maintenance (NGAM) is chosen adequately. When changing a constant feed of maize silage from continuous to pulsed feeding, the final average methane production remains very similar for both standard and coupled models, while both the initial response of the methanogenic population at the onset of pulsed feeding as well as its dynamics between pulses deviates considerably. In contrast to ADM1, the coupled models deliver predictions of up to 1,000s of intracellular metabolic fluxes per species, describing intracellular metabolic pathway activity in much higher detail. Furthermore, yield coefficients which need to be specified in ADM1 are no longer required as they are implicitly encoded in the topology of the species’ metabolic network. We show the feasibility of augmenting ADM1, an ordinary differential equation-based model for simulating biogas production, by FBA models implementing individual steps of anaerobic digestion. While cellular maintenance is introduced as a new parameter, the total number of parameters is reduced as yield coefficients no longer need to be specified. The coupled models provide detailed predictions on intracellular activity of microbial species which are compatible with experimental data on enzyme synthesis activity or abundance as obtained by metatranscriptomics or metaproteomics. By providing predictions of intracellular fluxes of individual community members, the presented approach advances the simulation of microbial community driven processes and provides a direct link to validation by state-of-the-art experimental techniques

Characterizing steady states of genome-scale metabolic networks in continuous cell cultures

We present a model for continuous cell culture coupling intra-cellular metabolism to extracellular variables describing the state of the bioreactor, taking into account the growth capacity of the cell and the impact of toxic byproduct accumulation. We provide a method to determine the steady states of this system that is tractable for metabolic networks of arbitrary complexity. We demonstrate our approach in a toy model first, and then in a genome-scale metabolic network of the Chinese hamster ovary cell line, obtaining results that are in qualitative agreement with experimental observations. More importantly, we derive a number of consequences from the model that are independent of parameter values. First, that the ratio between cell density and dilution rate is an ideal control parameter to fix a steady state with desired metabolic properties invariant across perfusion systems. This conclusion is robust even in the presence of multi-stability, which is explained in our model by the negative feedback loop on cell growth due to toxic byproduct accumulation. Moreover, a complex landscape of steady states in continuous cell culture emerges from our simulations, including multiple metabolic switches, which also explain why cell-line and media benchmarks carried out in batch culture cannot be extrapolated to perfusion. On the other hand, we predict invariance laws between continuous cell cultures with different parameters. A practical consequence is that the chemostat is an ideal experimental model for large-scale high-density perfusion cultures, where the complex landscape of metabolic transitions is faithfully reproduced. Thus, in order to actually reflect the expected behavior in perfusion, performance benchmarks of cell-lines and culture media should be carried out in a chemostat

Mechanistic understanding of mixed-culture fermentations by metabolic modelling

Biorefineries are set to become an important agent in the shift towards a circular economy due to their potential to valorise organic wastes into marketable products. Anaerobic fermentations yielding volatile fatty acids are a key process in this production scheme as their products act as intermediates between the organic wastes and the final biorefinery products. However, their product selectivity is highly influenced by the environmental conditions and the mechanisms governing the process remain unknown. In this thesis, predictive tools were developed with the objective of understanding the mechanisms governing anaerobic fermentations and of designing processes targeting specific volatile fatty acids with high productivity

Macromolecular crowding explains overflow metabolism in cells

Overflow metabolism is a metabolic phenotype of cells characterized by mixed oxidative phosphorylation (OxPhos) and fermentative glycolysis in the presence of oxygen. Recently, it was proposed that a combination of a protein allocation constraint and a higher proteome fraction cost of energy generation by OxPhos relative to fermentation form the basis of overflow metabolism in the bacterium, Escherichia coli. However, we argue that the existence of a maximum or optimal macromolecular density is another essential requirement. Here we re-evaluate our previous theory of overflow metabolism based on molecular crowding following the proteomic fractions formulation. We show that molecular crowding is a key factor in explaining the switch from OxPhos to overflow metabolism

Constrained Allocation Flux Balance Analysis

New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an "ensemble averaging" procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws.Comment: 21 pages, 6 figures (main) + 33 pages, various figures and tables (supporting); for the supplementary MatLab code, see http://tinyurl.com/h763es

Modeling human gut microbiota: from steady states to dynamic systems

Human gut microbes are an essential part of human sub-microscopic systems and involved in many critical biological processes such as Type 2 diabetes (T2D) and osteoporosis. However, the underlying mechanisms are unclear. Several mathematical modeling approaches, such as genome-scale metabolic models (GEMs) and ordinary differential equation (ODE) based models, have been used to simulate the dynamics of human gut microbiota. This thesis aims to explore, simulate, and predict the behavior of gut microbial ecosystems and the relationships between gut microbes and humans by modeling.The importance of the gut microbiome for bone metabolism and T2D has been demonstrated in mice and human cohorts. We first reconstructed a GEM for Limosilactobacillus reuteri ATCC PTA 6475, which is a probiotic that significantly reduces bone loss in older women with lower bone mineral density. To investigate the associations between T2D and the gut microbiota, GEMs for 827 gut microbial species and 1,779 community-level GEMs for T2D cohorts have also been constructed. With these GEMs, we investigated metabolic potentials such as short-chain fatty acids, amino acids, and vitamins that play vital roles in the host metabolism regulation. Furthermore, the integration of the models with machine learning method provides potential insights into the possible roles of gut microbiota in T2D.Cybernetic models, which simulate metabolic rates by integrating the control of enzyme synthesis and enzyme activities, have been applied to explore the dynamic behaviors of small-size metabolic networks. However, only a few studies have applied cybernetic theory to the microbial community so far. The remaining part of this thesis focuses on the use of cybernetic models to explore human gut microbiota\u27s interactions and population dynamics. Considering the high computing burden of the current cybernetic modeling approach for processing the full-size GEMs, we have developed a computing-efficient strategy for model reconstruction and simulation to reveal the metabolic dynamics of human gut microbiota.In this thesis, we explore the human gut microbiota from single L. reuteri species to microbial gut communities, from simple steady state systems by GEMs to complex dynamic systems by cybernetic model

A physical model of cell metabolism

Cell metabolism is characterized by three fundamental energy demands: to sustain cell maintenance, to trigger aerobic fermentation and to achieve maximum metabolic rate. The transition to aerobic fermentation and the maximum metabolic rate are currently understood based on enzymatic cost constraints. Yet, we are lacking a theory explaining the maintenance energy demand. Here we report a physical model of cell metabolism that explains the origin of these three energy scales. Our key hypothesis is that the maintenance energy demand is rooted on the energy expended by molecular motors to fluidize the cytoplasm and counteract molecular crowding. Using this model and independent parameter estimates we make predictions for the three energy scales that are in quantitative agreement with experimental values. The model also recapitulates the dependencies of cell growth with extracellular osmolarity and temperature. This theory brings together biophysics and cell biology in a tractable model that can be applied to understand key principles of cell metabolism