Ribosomes are responsible for the synthesis of proteins, the major component of cellular biomass. Classical experiments have established a linear relationship between the fraction of resources invested in ribosomal proteins and the rate of balanced growth of a microbial population. Very little is known, however, about how the investment in ribosomes varies over individual cells in a population. We therefore extended the study of ribosomal resource allocation from populations to single cells, using a combination of time-lapse fluorescence microscopy and statistical inference. We found a large variability of ribosome concentrations and growth rates in conditions of balanced growth of the model bacterium Escherichia coli in a given medium, which cannot be accounted for by the population-level growth law. A large variability in the allocation of resources to ribosomes was also found during the transition of the bacteria from a poor to a rich growth medium. While some cells immediately adapt their ribosome synthesis rate to the new environment, others do so only gradually. Our results thus reveal a range of strategies for investing resources in the molecular machines at the heart of cellular self-replication. This raises the fundamental question whether the observed variability is an intrinsic consequence of the stochastic nature of the underlying biochemical processes or whether it improves the fitness of Escherichia coli in its natural environment.

Using microorganisms for bioproduction requires the reorientation of metabolic fluxes from biomass synthesis to the production of compounds of interest. We previously engineered a synthetic growth switch in Escherichia coli based on inducible expression of the β- and β'-subunits of RNA polymerase. Depending on the level of induction, the cells stop growing or grow at a rate close to that of the wild-type strain. This strategy has been successful in transforming growth-arrested bacteria into biofactories with a high production yield, releasing cellular resources from growth towards biosynthesis. However, high selection pressure is placed on a growth-arrested population, favoring mutations that allow cells to escape from growth control. Accordingly, we made the design of the growth switch more robust by building in genetic redundancy. More specifically, we added the rpoA gene, encoding for the α-subunit of RNA polymerase, under the control of a copy of the same inducible promoter used for expression control of ββ'. The improved growth switch is much more stable (escape frequency <10-9), while preserving the capacity to improve production yields. Moreover, after a long period of growth inhibition the population can be regenerated within a few generations. This opens up the possibility to alternate biomass accumulation and product synthesis over a longer period of time and is an additional step towards the dynamical control of bioproduction.

Ribosomes are responsible for the synthesis of proteins, the major component of cellular biomass. Classical experiments have established a linear relationship between the fraction of resources invested in ribosomal proteins and the rate of balanced growth of a microbial population. We extended the study of ribosomal resource allocation from populations to single cells, using a combination of time-lapse fluorescence microscopy and statistical inference. We found a large variability of ribosome concentrations and growth rates in conditions of balanced growth of the model bacterium Escherichia coli. Moreover, the ribosome concentrations and growth rates of individual cells are uncorrelated, contrary to what would be expected from the population-level growth law. A similar large heterogeneity was found during the transition of the bacteria from a poor to a rich growth medium. Whereas some cells immediately adapt ribosomal resource allocation to the new environment, others do so only gradually. Our results thus reveal distinct strategies for investing resources in the molecular machines at the heart of cellular self-replication. This raises the interesting question whether the observed variability is an intrinsic consequence of the stochastic nature of the underlying biochemical processes or whether it improves the fitness of Escherichia coli in its natural environment. ### Competing Interest Statement The authors have declared no competing interest.

Fluorescent proteins (FPs) are a powerful tool to quantitatively monitor gene expression. The dynamics of a promoter and its regulation can be inferred from fluorescence data. The interpretation of fluorescent data, however, is strongly dependent on the maturation of FPs since different proteins mature in distinct ways. We propose a novel approach for analyzing fluorescent reporter data by incorporating maturation dynamics in the reconstruction of promoter activities. Our approach consists of developing and calibrating mechanistic maturation models for distinct FPs. These models are then used alongside a Bayesian approach to estimate promoter activities from fluorescence data. We demonstrate by means of targeted experiments in Escherichia coli that our approach provides robust estimates and that accounting for maturation is, in many cases, essential for the interpretation of gene expression data.

This is a document to group thoughts, ideas and results that come out of the meetings we discuss about the microfluidics data. Feel free to add your thoughts and remarks.

We investigated the scalability of a previously developed growth switch based on external control of RNA polymerase expression. Our results indicate that, in liter-scale bioreactors operating in fed-batch mode, growth-arrested Escherichia coli cells are able to convert glucose to glycerol at an increased yield. A multiomics quantification of the physiology of the cells shows that, apart from acetate production, few metabolic side effects occur. However, a number of specific responses to growth slow-down and growth arrest are launched at the transcriptional level. These notably include the downregulation of genes involved in growth-associated processes, such as amino acid and nucleotide metabolism and translation. Interestingly, the transcriptional responses are buffered at the proteome level, probably due to the strong decrease of the total mRNA concentration after the diminution of transcriptional activity and the absence of growth dilution of proteins. Growth arrest thus reduces the opportunities for dynamically adjusting the proteome composition, which poses constraints on the design of biotechnological production processes but may also avoid the initiation of deleterious stress responses.

During aerobic growth on glucose, Escherichia coli excretes acetate, a mechanism called "overflow metabolism." At high concentrations, the secreted acetate inhibits growth. Several mechanisms have been proposed for explaining this phenomenon, but a thorough analysis is hampered by the diversity of experimental conditions and strains used in these studies. Here, we describe the construction of a set of isogenic strains that remove different parts of the metabolic network involved in acetate metabolism. Analysis of these strains reveals that (i) high concentrations of acetate in the medium inhibit growth without significantly perturbing central metabolism; (ii) growth inhibition persists even when acetate assimilation is completely blocked; and (iii) regulatory interactions mediated by acetyl-phosphate play a small but significant role in growth inhibition by acetate. The major contribution to growth inhibition by acetate may originate in systemic effects like the uncoupling effect of organic acids or the perturbation of the anion composition of the cell, as previously proposed. Our data suggest, however, that under the conditions considered here, the uncoupling effect plays only a limited role.IMPORTANCE High concentrations of organic acids such as acetate inhibit growth of Escherichia coli and other bacteria. This phenomenon is of interest for understanding bacterial physiology but is also of practical relevance. Growth inhibition by organic acids underlies food preservation and causes problems during high-density fermentation in biotechnology. What causes this phenomenon? Classical explanations invoke the uncoupling effect of acetate and the establishment of an anion imbalance. Here, we propose and investigate an alternative hypothesis: the perturbation of acetate metabolism due to the inflow of excess acetate. We find that this perturbation accounts for 20% of the growth-inhibitory effect through a modification of the acetyl phosphate concentration. Moreover, we argue that our observations are not expected based on uncoupling alone.

BACKGROUND:Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolised. Although this system is one of the paradigms of regulation in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell - responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signalling, has motivated important modelling efforts over the past four decades, especially in the enterobacterium Escherichia coli.RESULTS:Starting from a simple core model with only four intracellular metabolites, we develop an ensemble of model variants, all showing diauxic growth behaviour during a batch process. The model variants fall into one of the four categories: flux balance models, kinetic models with growth dilution, kinetic models with regulation, and resource allocation models. The model variants differ from one another in only a single aspect, each breaking the symmetry between the two substrate assimilation pathways in a different manner, and can be quantitatively compared using a so-called diauxic growth index. For each of the model variants, we predict the behaviour in two new experimental conditions, namely a glucose pulse for a culture growing in minimal medium with lactose and a batch culture with different initial concentrations of the components of the transport systems. When qualitatively comparing these predictions with experimental data for these two conditions, a number of models can be excluded while other model variants are still not discriminable. The best-performing model variants are based on inducer inclusion and activation of enzymatic genes by a global transcription factor, but the other proposed factors may complement these well-known regulatory mechanisms.CONCLUSIONS:The model ensemble presented here offers a better understanding of the variety of mechanisms that have been proposed to play a role in CCR. In addition, it provides an educational resource for systems biology, as it gives an introduction to a broad range of modeling approaches in the context of a simple but biologically relevant example.

Bacteria have evolved complex regulatory networks to control the activity of transcription and translation, and thus the growth rate, over a range of environmental conditions. Reengineering RNA polymerase and ribosomes allows modifying naturally evolved regulatory networks and thereby profoundly reorganizing the manner in which bacteria allocate resources to different cellular functions. This opens new opportunities for our fundamental understanding of microbial physiology and for a variety of applications. Recent breakthroughs in genome engineering and the miniaturization and automation of culturing methods have offered new perspectives for the reengineering of the transcription and translation machinery in bacteria as well as the development of novel in vitro and in vivo gene expression systems. We review different examples from the unifying perspective of resource reallocation, and discuss the impact of these approaches for microbial systems biology and biotechnological applications.

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.

The recently published model of Bosdriesz et al. [1] provides a synthesis of the currently available knowledge of the ppGpp regulatory system. Through the mechanisms of ppGpp production and degradation, it describes regulation of the synthesis of ribosomal RNA. We explain below how we use the model to compare the action of the ppGpp system with the on-off control strategy. The denomination of variables and parameters follows the Supporting Information of [1] , and is reproduced in Table S4.1 in order to make the text self-contained. The evolution of the cellular concentration of ppGpp is described in [1] by

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.

Motivation: Time-series observations from reporter gene experiments are commonly used for inferring and analyzing dynamical models of regulatory networks. The robust estimation of promoter activities and protein concentrations from primary data is a difficult problem due to measurement noise and the indirect relation between the measurements and quantities of biological interest. Results: We propose a general approach based on regularized linear inversion to solve a range of estimation problems in the analysis of reporter gene data, notably the inference of growth rate, promoter activity, and protein concentration profiles. We evaluate the validity of the approach using in silico simulation studies, and observe that the methods are more robust and less biased than indirect approaches usually encountered in the experimental literature based on smoothing and subsequent processing of the primary data. We apply the methods to the analysis of fluorescent reporter gene data acquired in kinetic experiments with Escherichia coli. The methods are capable of reliably reconstructing time-course profiles of growth rate, promoter activity and protein concentration from weak and noisy signals at low population volumes. Moreover, they capture critical features of those profiles, notably rapid changes in gene expression during growth transitions. Availability and implementation: The methods described in this article are made available as a Python package (LGPL license) and also accessible through a web interface. For more information, see https://team.inria.fr/ibis/wellinverter. Contact: Hidde.de-Jong@inria.fr Supplementary information: Supplementary data are available at Bioinformatics online.

The ability to control growth is essential for fundamental studies of bacterial physiology and biotechnological applications. We have engineered an Escherichia coli strain in which the transcription of a key component of the gene expression machinery, RNA polymerase, is under the control of an inducible promoter. By changing the inducer concentration in the medium, we can adjust the RNA polymerase concentration and thereby switch bacterial growth between zero and the maximal growth rate supported by the medium. We show that our synthetic growth switch functions in a medium‐independent and reversible way, and we provide evidence that the switching phenotype arises from the ultrasensitive response of the growth rate to the concentration of RNA polymerase. We present an application of the growth switch in which both the wild‐type E. coli strain and our modified strain are endowed with the capacity to produce glycerol when growing on glucose. Cells in which growth has been switched off continue to be metabolically active and harness the energy gain to produce glycerol at a twofold higher yield than in cells with natural control of RNA polymerase expression. Remarkably, without any further optimization, the improved yield is close to the theoretical maximum computed from a flux balance model of E. coli metabolism. The proposed synthetic growth switch is a promising tool for gaining a better understanding of bacterial physiology and for applications in synthetic biology and biotechnology. Bacterial growth can be switched on/off by expressing the RNA polymerase large subunits under the control of an inducible promoter. When growth is switched off, cells reorient nutrient fluxes to increase the production of a metabolite of interest. Bacterial growth can be switched on/off by expressing the RNA polymerase large subunits under the control of an inducible promoter. When growth is switched off, cells reorient nutrient fluxes to increase the production of a metabolite of interest.

The inference of regulatory interactions and quantitative models of gene regulation from time-series transcriptomics data has been extensively studied and applied to a range of problems in drug discovery, cancer research, and biotechnology. The application of existing methods is commonly based on implicit assumptions on the biological processes under study. First, the measurements of mRNA abundance obtained in transcriptomics experiments are taken to be representative of protein concentrations. Second, the observed changes in gene expression are assumed to be solely due to transcription factors and other specific regulators, while changes in the activity of the gene expression machinery and other global physiological effects are neglected. While convenient in practice, these assumptions are often not valid and bias the reverse engineering process. Here we systematically investigate, using a combination of models and experiments, the importance of this bias and possible corrections. We measure in real time and in vivo the activity of genes involved in the FliA-FlgM module of the E. coli motility network. From these data, we estimate protein concentrations and global physiological effects by means of kinetic models of gene expression. Our results indicate that correcting for the bias of commonly-made assumptions improves the quality of the models inferred from the data. Moreover, we show by simulation that these improvements are expected to be even stronger for systems in which protein concentrations have longer half-lives and the activity of the gene expression machinery varies more strongly across conditions than in the FliA-FlgM module. The approach proposed in this study is broadly applicable when using time-series transcriptome data to learn about the structure and dynamics of regulatory networks. In the case of the FliA-FlgM module, our results demonstrate the importance of global physiological effects and the active regulation of FliA and FlgM half-lives for the dynamics of FliA-dependent promoters.

In bacteria, selective promoter recognition by RNA polymerase is achieved by its association with σ factors, accessory subunits able to direct RNA polymerase "core enzyme" (E) to different promoter sequences. Using Chromatin Immunoprecipitation-sequencing (ChIP-seq), we searched for promoters bound by the σ(S)-associated RNA polymerase form (Eσ(S)) during transition from exponential to stationary phase. We identified 63 binding sites for Eσ(S) overlapping known or putative promoters, often located upstream of genes (encoding either ORFs or non-coding RNAs) showing at least some degree of dependence on the σ(S)-encoding rpoS gene. Eσ(S) binding did not always correlate with an increase in transcription level, suggesting that, at some σ(S)-dependent promoters, Eσ(S) might remain poised in a pre-initiation state upon binding. A large fraction of Eσ(S)-binding sites corresponded to promoters recognized by RNA polymerase associated with σ(70) or other σ factors, suggesting a considerable overlap in promoter recognition between different forms of RNA polymerase. In particular, Eσ(S) appears to contribute significantly to transcription of genes encoding proteins involved in LPS biosynthesis and in cell surface composition. Finally, our results highlight a direct role of Eσ(S) in the regulation of non coding RNAs, such as OmrA/B, RyeA/B and SibC.

The understanding of the regulatory networks controlling the adaptation of bacteria to changes in their environment is critically dependent on the ability to monitor the dynamics of gene expression. Here, we review the use of fluorescent reporter genes for dynamically quantifying promoter activity and other quantities characterizing gene expression. We discuss critical physical and biological parameters in the design, development, and use of fluorescent reporter strains. Moreover, we review measurement models that have been proposed to interpret primary fluorescence data and inference methods for estimating gene expression profiles from these data. As an illustration of the use of fluorescent reporter strains for analyzing bacterial regulatory networks, we consider two applications in the model bacterium Escherichia coli in some detail: the joint control of gene expression by global physiological effects and specific regulatory interactions, and the importance of protein stability for the inference and analysis of transcriptional regulatory networks. We conclude by discussing some current trends in the use of fluorescent reporter genes.

Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolized. Although this system is one of the paradigms of the regulation of gene expression in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell that are responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signaling, and that are subject to global physical and physiological constraints, has motivated important modeling efforts over the past four decades, especially in the enterobacterium Escherichia coli. Different hypotheses concerning the dynamic functioning of the system have been explored by a variety of modeling approaches. We review these studies and summarize their contributions to the quantitative understanding of CCR, focusing on diauxic growth in E. coli. Moreover, we propose a highly simplified representation of diauxic growth that makes it possible to bring out the salient features of the models proposed in the literature and confront and compare the explanations they provide.