Search

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Escherichia coli adapts its lifestyle to the variations of environmental growth conditions, swapping between swimming motility or biofilm formation. The stationary-phase sigma factor RpoS is an important regulator of this switch, since it stimulates adhesion and represses flagellar biosynthesis. By measuring the dynamics of gene expression, we show that RpoS inhibits the transcription of the flagellar sigma factor, FliA, in exponential growth phase. RpoS also partially controls the expression of CsgD and CpxR, two transcription factors important for bacterial adhesion. We demonstrate that these two regulators repress the transcription of fliA, flgM, and tar and that this regulation is dependent on the growth medium. CsgD binds to the flgM and fliA promoters around their -10 promoter element, strongly suggesting direct repression. We show that CsgD and CpxR also affect the expression of other known modulators of cell motility. We propose an updated structure of the regulatory network controlling the choice between adhesion and motility. [ABSTRACT FROM AUTHOR]

Background: Fluorescent and luminescent reporter genes have become popular tools for the real-time monitoring of gene expression in living cells. However, mathematical models are necessary for extracting biologically meaningful quantities from the primary data. Results: We present a rigorous method for deriving relative protein synthesis rates (mRNA concentrations) and protein concentrations by means of kinetic models of gene expression. We experimentally and computationally validate this approach in the case of the protein Fis, a global regulator of transcription in Escherichia coli. We show that the mRNA and protein concentration profiles predicted from the models agree quite well with direct measurements obtained by Northern and Western blots, respectively. Moreover, we present computational procedures for taking into account systematic biases like the folding time of the fluorescent reporter protein and differences in the half-lives of reporter and host gene products. The results show that large differences in protein half-lives, more than mRNA half-lives, may be critical for the interpretation of reporter gene data in the analysis of the dynamics of regulatory systems. Conclusions: The paper contributes to the development of sound methods for the interpretation of reporter gene data, notably in the context of the reconstruction and validation of models of regulatory networks. The results have wide applicability for the analysis of gene expression in bacteria and may be extended to higher organisms. [ABSTRACT FROM AUTHOR]

Motivation: Fluorescent and luminescent reporter gene systems in combination with automated microplate readers allow real-time monitoring of gene expression on the population level at high precision and sampling density. This generates large amounts of data for the analysis of which computer tools are missing to date. Results: We have developed WellReader, a MATLAB program for the analysis of fluorescent and luminescent reporter gene data. WellReader allows the user to load the output files of microplate readers, remove outliers, correct for background effects and smooth and fit the data. Moreover, it computes biologically relevant quantities from the measured signals, notably promoter activities and protein concentrations, and compares the resulting expression profiles of different genes under different conditions. [ABSTRACT FROM AUTHOR]

Abstract: The alternative sigma factor RpoS is a central regulator of the stress response in many Proteobacteria, acting both during exponential growth and in stationary phase. The small protein Crl increases the interaction between RpoS and RNA polymerase and thereby activates certain RpoS-dependent promoters. However, the growth-phase dependence of the interaction of Crl with different forms of polymerase remains unknown. We use 41 GFP transcriptional fusions to study the dynamics of gene regulation by RpoS and Crl during growth transition from exponential to stationary phase in Escherichia coli. We confirm that RpoS can regulate gene expression in exponential phase, both positively and negatively. Crl slightly stimulates transcription by RpoS in exponential phase and controls a subset of RpoS-dependent genes in stationary phase. Growth temperature strongly affects induction of specific promoters by RpoS, whereas its impact on gene regulation by Crl is much less significant. In addition, we identify five new genes regulated by Crl (ada, cbpA, glgS, sodC and flgM) and demonstrate that Crl improves promoter binding and opening by RpoS-containing RNA polymerase at the hdeA promoter. Our study also shows that Crl is a cognate enhancer of RpoS activity under different growth conditions, since its deletion has no effect on genes transcribed by other sigma factors. [Copyright &y& Elsevier]

Abstract: Gene regulatory networks consist of direct interactions, but also include indirect interactions mediated by metabolism. We investigate to which extent these indirect interactions arising from metabolic coupling influence the dynamics of the system. To this end, we build a qualitative model of the gene regulatory network controlling carbon assimilation in Escherichia coli, and use this model to study the changes in gene expression following a diauxic shift from glucose to acetate. In particular, we compare the relative variation in the steady-state concentrations of enzymes and transcription regulators during growth on glucose and acetate, as well as the dynamic response of gene expression to the exhaustion of glucose and the subsequent assimilation of acetate. We find significant differences between the dynamics of the system in the absence and presence of metabolic coupling. This shows that interactions arising from metabolic coupling cannot be ignored when studying the dynamics of gene regulatory networks. [Copyright &y& Elsevier]

The alternative sigma factor σs (RpoS) of Escherichia coli RNA polymerase regulates the expression of stationary phase and stress-response genes, σs is also required for the transcription of the cryptic genes csgBA that encode the subunits of the curli proteins. The expression of the csgBA genes is regulated in response to a multitude of physiological signals. In stationary phase, these genes are transcribed by the σs factor, and expression of the operon is enhanced by the small protein Crl. It has been shown that Crl stimulates the activity of σs, leading to an increased transcription rate of a subset of genes of the rpoS regulon in stationary phase. However, the underlying molecular mechanism has remained elusive. We show here that Crl interacts directly with σs and that this interaction promotes binding of the σs holoenzyme (Eσs) to the csgBA promoter. Expression of Crl is increased during the transition from growing to stationary phase. Crl accumulates in stationary phase cells at low temperature (30 °C) but not at 37 °C. We therefore propose that Crl is a second thermosensor, besides DsrA, controlling σs activity. [ABSTRACT FROM AUTHOR]

Twelve populations of Escherichia coli evolved in and adapted to a glucose-limited environment from a common ancestor. We used two-dimensional protein electrophoresis to compare two evolved clones, isolated from independently derived populations after 20,000 generations. Exceptional parallelism was detected. We compared the observed changes in protein expression profiles with previously characterized global transcription profiles of the same clones; this is the first time such a comparison has been made in an evolutionary context where these changes are often quite subtle. The two methodologies exhibited some remarkable similarities that highlighted two different levels of parallel regulatory changes that were beneficial during the evolution experiment. First, at the higher level, both methods revealed extensive parallel changes in the same global regulatory network, reflecting the involvement of beneficial mutations in genes that control the ppGpp regulon. Second, both methods detected expression changes of identical gene sets that reflected parallel changes at a lower level of gene regulation. The protein profiles led to the discovery of beneficial mutations affecting the malT gene, with strong genetic parallelism across independently evolved populations. Functional and evolutionary analyses of these mutations revealed parallel phenotypic decreases in the maltose regulon expression and a high level of polymorphism at this locus in the evolved populations. [ABSTRACT FROM AUTHOR]