A synthetic Escherichia coli predator–prey ecosystem

We have constructed a synthetic ecosystem consisting of two Escherichia coli populations, which communicate bi-directionally through quorum sensing and regulate each other's gene expression and survival via engineered gene circuits. Our synthetic ecosystem resembles canonical predator–prey systems in terms of logic and dynamics. The predator cells kill the prey by inducing expression of a killer protein in the prey, while the prey rescue the predators by eliciting expression of an antidote protein in the predator. Extinction, coexistence and oscillatory dynamics of the predator and prey populations are possible depending on the operating conditions as experimentally validated by long-term culturing of the system in microchemostats. A simple mathematical model is developed to capture these system dynamics. Coherent interplay between experiments and mathematical analysis enables exploration of the dynamics of interacting populations in a predictable manner

Modelling and Analysis of Central Metabolism Operating Regulatory Interactions in Salt Stress Conditions in a L-Carnitine Overproducing E. coli Strain

Based on experimental data from E. coli cultures, we have devised a mathematical model in the GMA-power law formalism that describes the central and L-carnitine metabolism in and between two steady states, non-osmotic and hyperosmotic (0.3 M NaCl). A key feature of this model is the introduction of type of kinetic order, the osmotic stress kinetic orders (gOSn), derived from the power law general formalism, which represent the effect of osmotic stress in each metabolic process of the model

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators

Monitoring Parallel Robotic Cultivations with Online Multivariate Analysis

In conditional microbial screening, a limited number of candidate strains are tested at different conditions searching for the optimal operation strategy in production (e.g., temperature and pH shifts, media composition as well as feeding and induction strategies). To achieve this, cultivation volumes of >10 mL and advanced control schemes are required to allow appropriate sampling and analyses. Operations become even more complex when the analytical methods are integrated into the robot facility. Among other multivariate data analysis methods, principal component analysis (PCA) techniques have especially gained popularity in high throughput screening. However, an important issue specific to high throughput bioprocess development is the lack of so-called golden batches that could be used as a basis for multivariate analysis. In this study, we establish and present a program to monitor dynamic parallel cultivations in a high throughput facility. PCA was used for process monitoring and automated fault detection of 24 parallel running experiments using recombinant E. coli cells expressing three different fluorescence proteins as the model organism. This approach allowed for capturing events like stirrer failures and blockage of the aeration system and provided a good signal to noise ratio. The developed application can be easily integrated in existing data- and device-infrastructures, allowing automated and remote monitoring of parallel bioreactor systems.BMBF, 031L0018A, ERASysApp2 - Verbundprojekt: LEANPROT - Entwicklung einer Systembiologie-Plattform für die Entwicklung von lean-proteome-Escherichia coli-Stämmen - Deutsches Teilprojekt ADFG, 414044773, Open Access Publizieren 2019 - 2020 / Technische Universität Berli

A Coarse-Grained Biophysical Model of E. coli and Its Application to Perturbation of the rRNA Operon Copy Number

We propose a biophysical model of Escherichia coli that predicts growth rate and an effective cellular composition from an effective, coarse-grained representation of its genome. We assume that E. coli is in a state of balanced exponential steadystate growth, growing in a temporally and spatially constant environment, rich in resources. We apply this model to a series of past measurements, where the growth rate and rRNA-to-protein ratio have been measured for seven E. coli strains with an rRNA operon copy number ranging from one to seven (the wild-type copy number). These experiments show that growth rate markedly decreases for strains with fewer than six copies. Using the model, we were able to reproduce these measurements. We show that the model that best fits these data suggests that the volume fraction of macromolecules inside E. coli is not fixed when the rRNA operon copy number is varied. Moreover, the model predicts that increasing the copy number beyond seven results in a cytoplasm densely packed with ribosomes and proteins. Assuming that under such overcrowded conditions prolonged diffusion times tend to weaken binding affinities, the model predicts that growth rate will not increase substantially beyond the wild-type growth rate, as indicated by other experiments. Our model therefore suggests that changing the rRNA operon copy number of wild-type E. coli cells growing in a constant rich environment does not substantially increase their growth rate. Other observations regarding strains with an altered rRNA operon copy number, such as nucleoid compaction and the rRNA operon feedback response, appear to be qualitatively consistent with this model. In addition, we discuss possible design principles suggested by the model and propose further experiments to test its validity