Optimal protein production by a synthetic microbial consortium: coexistence, distribution of labor, and syntrophy.

The bacterium E. coli is widely used to produce recombinant proteins such as growth hormone and insulin. One inconvenience with E. coli cultures is the secretion of acetate through overflow metabolism. Acetate inhibits cell growth and represents a carbon diversion, which results in several negative effects on protein production. One way to overcome this problem is the use of a synthetic consortium of two different E. coli strains, one producing recombinant proteins and one reducing the acetate concentration. In this paper, we study a mathematical model of such a synthetic community in a chemostat where both strains are allowed to produce recombinant proteins. We give necessary and sufficient conditions for the existence of a coexistence equilibrium and show that it is unique. Based on this equilibrium, we define a multi-objective optimization problem for the maximization of two important bioprocess performance metrics, process yield and productivity. Solving numerically this problem, we find the best available trade-offs between the metrics. Under optimal operation of the mixed community, both strains must produce the protein of interest, and not only one (distribution instead of division of labor). Moreover, in this regime acetate secretion by one strain is necessary for the survival of the other (syntrophy). The results thus illustrate how complex multi-level dynamics shape the optimal production of recombinant proteins by synthetic microbial consortia.

Resource allocation accounts for the large variability of rate-yield phenotypes across bacterial strains.

Different strains of a microorganism growing in the same environment display a wide variety of growth rates and growth yields. We developed a coarse-grained model to test the hypothesis that different resource allocation strategies, corresponding to different compositions of the proteome, can account for the observed rate-yield variability. The model predictions were verified by means of a database of hundreds of published rate-yield and uptake-secretion phenotypes of Escherichia coli strains grown in standard laboratory conditions. We found a very good quantitative agreement between the range of predicted and observed growth rates, growth yields, and glucose uptake and acetate secretion rates. These results support the hypothesis that resource allocation is a major explanatory factor of the observed variability of growth rates and growth yields across different bacterial strains. An interesting prediction of our model, supported by the experimental data, is that high growth rates are not necessarily accompanied by low growth yields. The resource allocation strategies enabling high-rate, high-yield growth of E. coli lead to a higher saturation of enzymes and ribosomes, and thus to a more efficient utilization of proteomic resources. Our model thus contributes to a fundamental understanding of the quantitative relationship between rate and yield in E. coli and other microorganisms. It may also be useful for the rapid screening of strains in metabolic engineering and synthetic biology.

Optimal proteome allocation and the temperature dependence of microbial growth laws.

Although the effect of temperature on microbial growth has been widely studied, the role of proteome allocation in bringing about temperature-induced changes remains elusive. To tackle this problem, we propose a coarse-grained model of microbial growth, including the processes of temperature-sensitive protein unfolding and chaperone-assisted (re)folding. We determine the proteome sector allocation that maximizes balanced growth rate as a function of nutrient limitation and temperature. Calibrated with quantitative proteomic data for Escherichia coli, the model allows us to clarify general principles of temperature-dependent proteome allocation and formulate generalized growth laws. The same activation energy for metabolic enzymes and ribosomes leads to an Arrhenius increase in growth rate at constant proteome composition over a large range of temperatures, whereas at extreme temperatures resources are diverted away from growth to chaperone-mediated stress responses. Our approach points at risks and possible remedies for the use of ribosome content to characterize complex ecosystems with temperature variation.

Global dynamics of the chemostat with overflow metabolism.

Fast growing E. coli cells, in glucose-aerobic conditions, excrete fermentation by-products such as acetate. This phenomenon is known as overflow metabolism and has been observed in a diverse range of microorganisms. In this paper, we study a chemostat model subject to overflow metabolism: if the substrate uptake rate (or the specific growth rate) is above a threshold rate (different from zero), then secretion of a by-product happens. We assume that the presence of the by-product has an inhibitory effect on the growth of the microorganism. The model is described by a non-smooth differential system of dimension three. We prove the existence of at most one equilibrium (or steady-state) with presence of microorganism, which is globally stable. We use these results to discuss the performance of chemostat-type systems to produce biomass or recombinant proteins.

Optimal bacterial resource allocation: metabolite production in continuous bioreactors.

We show novel results addressing the problem of synthesizing a metabolite of interest in continuous bioreactors through resource allocation control. Our approach is based on a coarse-grained self-replicator dynamical model that accounts for microbial culture growth inside the bioreactor, and incorporates a synthetic growth switch that allows to externally modify the RNA polymerase concentration of the bacterial population, thus disrupting the natural process of allocation of available resources in bacteria. Further on, we study its asymptotic behavior using dynamical systems theory, and we provide conditions for the persistence of the bacterial population. We aim to maximize the synthesis of the metabolite of interest during a fixed interval of time in terms of the resource allocation decision. The latter is formulated as an Optimal Control Problem which is then explored using Pontryagin's Maximum Principle. We analyze the solution of the problem and propose a sub-optimal control strategy given by a constant allocation decision, which eventually takes the system to the optimal steady-state production regime. On this basis, we study and compare the two most significant steady-state production objectives in continuous bioreactors: biomass production and metabolite production. For this last purpose, and in addition to the allocation parameter, we control the dilution rate of the bioreactor, and we analyze the results through a numerical approach. The resulting two-dimensional optimization problem is defined in terms of Michaelis-Menten kinetics, and takes into account the constraints for the existence of the equilibrium of interest.

Enhanced production of heterologous proteins by a synthetic microbial community: Conditions and trade-offs.

Synthetic microbial consortia have been increasingly utilized in biotechnology and experimental evidence shows that suitably engineered consortia can outperform individual species in the synthesis of valuable products. Despite significant achievements, though, a quantitative understanding of the conditions that make this possible, and of the trade-offs due to the concurrent growth of multiple species, is still limited. In this work, we contribute to filling this gap by the investigation of a known prototypical synthetic consortium. A first E. coli strain, producing a heterologous protein, is sided by a second E. coli strain engineered to scavenge toxic byproducts, thus favoring the growth of the producer at the expense of diverting part of the resources to the growth of the cleaner. The simplicity of the consortium is ideal to perform an in depth-analysis and draw conclusions of more general interest. We develop a coarse-grained mathematical model that quantitatively accounts for literature data from different key growth phenotypes. Based on this, assuming growth in chemostat, we first investigate the conditions enabling stable coexistence of both strains and the effect of the metabolic load due to heterologous protein production. In these conditions, we establish when and to what extent the consortium outperforms the producer alone in terms of productivity. Finally, we show in chemostat as well as in a fed-batch scenario that gain in productivity comes at the price of a reduced yield, reflecting at the level of the consortium resource allocation trade-offs that are well-known for individual species.

Reducing a model of sugar metabolism in peach to catch different patterns among genotypes.

Several studies have been conducted to understand the dynamic of primary metabolisms in fruit by translating them into mathematics models. An ODE kinetic model of sugar metabolism has been developed by Desnoues et al. (2018) to simulate the accumulation of different sugars during peach fruit development. Two major drawbacks of this model are (a) the number of parameters to calibrate and (b) its integration time that can be long due to non-linearity and time-dependent input functions. Together, these issues hamper the use of the model for a large panel of genotypes, for which few data are available. In this paper, we present a model reduction scheme that explicitly addresses the specificity of genetic studies in that: (i) it yields a reduced model that is adapted to the whole expected genetic diversity (ii) it maintains network structure and variable identity, in order to facilitate biological interpretation. The proposed approach is based on the combination and the systematic evaluation of different reduction methods. Thus, we combined multivariate sensitivity analysis, structural simplification and timescale-based approaches to simplify the number and the structure of ordinary differential equations of the model. The original and reduced models were compared based on three criteria, namely the corrected Aikake Information Criterion (AICC), the calibration time and the expected error of the reduced model over a progeny of virtual genotypes. The resulting reduced model not only reproduces the predictions of the original one but presents many advantages including a reduced number of parameters to be estimated and shorter calibration time, opening new promising perspectives for genetic studies and virtual breeding. The validity of the reduced model was further evaluated by calibration on 30 additional genotypes of an inter-specific peach progeny for which few data were available.

Analysis of a genetic-metabolic oscillator with piecewise linear models.

Interactions between gene regulatory networks and metabolism produce a diversity of dynamics, including multistability and oscillations. Here, we characterize a regulatory mechanism that drives the emergence of periodic oscillations in metabolic networks subject to genetic feedback regulation by pathway intermediates. We employ a qualitative formalism based on piecewise linear models to systematically analyze the behavior of gene-regulated metabolic pathways. For a pathway with two metabolites and three enzymes, we prove the existence of two co-existing oscillatory behaviors: damped oscillations towards a fixed point or sustained oscillations along a periodic orbit. We show that this mechanism closely resembles the "metabolator", a genetic-metabolic circuit engineered to produce autonomous oscillations in vivo.

Optimal control of bacterial growth for the maximization of metabolite production.

Microorganisms have evolved complex strategies for controlling the distribution of available resources over cellular functions. Biotechnology aims at interfering with these strategies, so as to optimize the production of metabolites and other compounds of interest, by (re)engineering the underlying regulatory networks of the cell. The resulting reallocation of resources can be described by simple so-called self-replicator models and the maximization of the synthesis of a product of interest formulated as a dynamic optimal control problem. Motivated by recent experimental work, we are specifically interested in the maximization of metabolite production in cases where growth can be switched off through an external control signal. We study various optimal control problems for the corresponding self-replicator models by means of a combination of analytical and computational techniques. We show that the optimal solutions for biomass maximization and product maximization are very similar in the case of unlimited nutrient supply, but diverge when nutrients are limited. Moreover, external growth control overrides natural feedback growth control and leads to an optimal scheme consisting of a first phase of growth maximization followed by a second phase of product maximization. This two-phase scheme agrees with strategies that have been proposed in metabolic engineering. More generally, our work shows the potential of optimal control theory for better understanding and improving biotechnological production processes.

Modeling the bioconversion of polysaccharides in a continuous reactor: A case study of the production of oligogalacturonates by Dickeya dadantii.

In the quest for a sustainable economy of the Earth's resources and for renewable sources of energy, a promising avenue is to exploit the vast quantity of polysaccharide molecules contained in green wastes. To that end, the decomposition of pectin appears to be an interesting target because this polymeric carbohydrate is abundant in many fruit pulps and soft vegetables. To quantitatively study this degradation process, here we designed a bioreactor that is continuously fed with de-esterified pectin (PGA). Thanks to the pectate lyases produced by bacteria cultivated in the vessel, the PGA is depolymerized into oligogalacturonates (UGA), which are continuously extracted from the tank. A mathematical model of our system predicted that the conversion efficiency of PGA into UGA increases in a range of coefficients of dilution until reaching an upper limit where the fraction of UGA that is extracted from the bioreactor is maximized. Results from experiments with a continuous reactor hosting a strain of the plant pathogenic bacterium Dickeya dadantii and in which the dilution coefficients were varied quantitatively validated the predictions of our model. A further theoretical analysis of the system enabled an a priori comparison of the efficiency of eight other pectate lyase-producing microorganisms with that of D. dadantii Our findings suggest that D. dadantii is the most efficient microorganism and therefore the best candidate for a practical implementation of our scheme for the bioproduction of UGA from PGA.

Principal process analysis of biological models.

BACKGROUND: Understanding the dynamical behaviour of biological systems is challenged by their large number of components and interactions. While efforts have been made in this direction to reduce model complexity, they often prove insufficient to grasp which and when model processes play a crucial role. Answering these questions is fundamental to unravel the functioning of living organisms. RESULTS: We design a method for dealing with model complexity, based on the analysis of dynamical models by means of Principal Process Analysis. We apply the method to a well-known model of circadian rhythms in mammals. The knowledge of the system trajectories allows us to decompose the system dynamics into processes that are active or inactive with respect to a certain threshold value. Process activities are graphically represented by Boolean and Dynamical Process Maps. We detect model processes that are always inactive, or inactive on some time interval. Eliminating these processes reduces the complex dynamics of the original model to the much simpler dynamics of the core processes, in a succession of sub-models that are easier to analyse. We quantify by means of global relative errors the extent to which the simplified models reproduce the main features of the original system dynamics and apply global sensitivity analysis to test the influence of model parameters on the errors. CONCLUSION: The results obtained prove the robustness of the method. The analysis of the sub-model dynamics allows us to identify the source of circadian oscillations. We find that the negative feedback loop involving proteins PER, CRY, CLOCK-BMAL1 is the main oscillator, in agreement with previous modelling and experimental studies. In conclusion, Principal Process Analysis is a simple-to-use method, which constitutes an additional and useful tool for analysing the complex dynamical behaviour of biological systems.

Reduction and Stability Analysis of a Transcription-Translation Model of RNA Polymerase.

The aim of this paper is to analyze the dynamical behavior of biological models of gene transcription and translation. We focus on a particular positive feedback loop governing the synthesis of RNA polymerase, needed for transcribing its own gene. We write a high-dimension model based on mass action laws and reduce it to a two-variable model (RNA polymerase and its mRNA) by means of monotone system theory and timescale arguments. We show that the reduced model has either a single globally stable trivial equilibrium in (0, 0), or an unstable zero equilibrium and a globally stable positive one. We give generalizations of this model, notably with a variable growth rate. The dynamical behavior of this system can be related to biological observations on the bacterium Escherichia coli.

Mathematical modelling of microbes: metabolism, gene expression and growth.

The growth of microorganisms involves the conversion of nutrients in the environment into biomass, mostly proteins and other macromolecules. This conversion is accomplished by networks of biochemical reactions cutting across cellular functions, such as metabolism, gene expression, transport and signalling. Mathematical modelling is a powerful tool for gaining an understanding of the functioning of this large and complex system and the role played by individual constituents and mechanisms. This requires models of microbial growth that provide an integrated view of the reaction networks and bridge the scale from individual reactions to the growth of a population. In this review, we derive a general framework for the kinetic modelling of microbial growth from basic hypotheses about the underlying reaction systems. Moreover, we show that several families of approximate models presented in the literature, notably flux balance models and coarse-grained whole-cell models, can be derived with the help of additional simplifying hypotheses. This perspective clearly brings out how apparently quite different modelling approaches are related on a deeper level, and suggests directions for further research.

Dynamical Allocation of Cellular Resources as an Optimal Control Problem: Novel Insights into Microbial Growth Strategies.

Microbial physiology exhibits growth laws that relate the macromolecular composition of the cell to the growth rate. Recent work has shown that these empirical regularities can be derived from coarse-grained models of resource allocation. While these studies focus on steady-state growth, such conditions are rarely found in natural habitats, where microorganisms are continually challenged by environmental fluctuations. The aim of this paper is to extend the study of microbial growth strategies to dynamical environments, using a self-replicator model. We formulate dynamical growth maximization as an optimal control problem that can be solved using Pontryagin's Maximum Principle. We compare this theoretical gold standard with different possible implementations of growth control in bacterial cells. We find that simple control strategies enabling growth-rate maximization at steady state are suboptimal for transitions from one growth regime to another, for example when shifting bacterial cells to a medium supporting a higher growth rate. A near-optimal control strategy in dynamical conditions is shown to require information on several, rather than a single physiological variable. Interestingly, this strategy has structural analogies with the regulation of ribosomal protein synthesis by ppGpp in the enterobacterium Escherichia coli. It involves sensing a mismatch between precursor and ribosome concentrations, as well as the adjustment of ribosome synthesis in a switch-like manner. Our results show how the capability of regulatory systems to integrate information about several physiological variables is critical for optimizing growth in a changing environment.

Links between topology of the transition graph and limit cycles in a two-dimensional piecewise affine biological model.

A class of piecewise affine differential (PWA) models, initially proposed by Glass and Kauffman (in J Theor Biol 39:103-129, 1973), has been widely used for the modelling and the analysis of biological switch-like systems, such as genetic or neural networks. Its mathematical tractability facilitates the qualitative analysis of dynamical behaviors, in particular periodic phenomena which are of prime importance in biology. Notably, a discrete qualitative description of the dynamics, called the transition graph, can be directly associated to this class of PWA systems. Here we present a study of periodic behaviours (i.e. limit cycles) in a class of two-dimensional piecewise affine biological models. Using concavity and continuity properties of Poincaré maps, we derive structural principles linking the topology of the transition graph to the existence, number and stability of limit cycles. These results notably extend previous works on the investigation of structural principles to the case of unequal and regulated decay rates for the 2-dimensional case. Some numerical examples corresponding to minimal models of biological oscillators are treated to illustrate the use of these structural principles.

Global stability of enzymatic chains of full reversible Michaelis-Menten reactions.

We consider a chain of metabolic reactions catalyzed by enzymes, of reversible Michaelis-Menten type with full dynamics, i.e. not reduced with any quasi-steady state approximations. We study the corresponding dynamical system and show its global stability if the equilibrium exists. If the system is open, the equilibrium may not exist. The main tool is monotone systems theory. Finally we study the implications of these results for the study of coupled genetic-metabolic systems.

Global stability of reversible enzymatic metabolic chains.

We consider metabolic networks with reversible enzymatic reactions. The model is written as a system of ordinary differential equations, possibly with inputs and outputs. We prove the global stability of the equilibrium (if it exists), using techniques of monotone systems and compartmental matrices. We show that the equilibrium does not always exist. Finally, we consider a metabolic system coupled with a genetic network, and we study the dependence of the metabolic equilibrium (if it exists) with respect to concentrations of enzymes. We give some conclusions concerning the dynamical behavior of coupled genetic/metabolic systems.

Stabilizing effect of cannibalism in a two stages population model.

In this paper we build a prey-predator model with discrete weight structure for the predator. This model will conserve the number of individuals and the biomass and both growth and reproduction of the predator will depend on the food ingested. Moreover the model allows cannibalism which means that the predator can eat the prey but also other predators. We will focus on a simple version with two weight classes or stage (larvae and adults) and present some general mathematical results. In the last part, we will assume that the dynamics of the prey is fast compared to the predator's one to go further in the results and eventually conclude that under some conditions, cannibalism can stabilize the system: more precisely, an unstable equilibrium without cannibalism will become almost globally stable with some cannibalism. Some numerical simulations are done to illustrate this result.

Probabilistic approach for predicting periodic orbits in piecewise affine differential models.

Piecewise affine models provide a qualitative description of the dynamics of a system, and are often used to study genetic regulatory networks. The state space of a piecewise affine system is partitioned into hyperrectangles, which can be represented as nodes in a directed graph, so that the system's trajectories follow a path in a transition graph. This paper proposes and compares two definitions of probability of transition between two nodes A and B of the graph, based on the volume of the initial conditions on the hyperrectangle A whose trajectories cross to B. The parameters of the system can thus be compared to the observed transitions between two hyperrectangles. This property may become useful to identify sets of parameters for which the system yields a desired periodic orbit with a high probability, or to predict the most likely periodic orbit given a set of parameters, as illustrated by a gene regulatory system composed of two intertwined negative loops.