DnaA and the timing of chromosome replication in Escherichia coli as a function of growth rate

The design of unit cells by combining the self-reproduction systems and metabolic cushioning loads

Recently, we published a comprehensive theoretical analysis of the self-reproduction processes in proto-cells (doubling of their components) composed of different combinations of cellular subsystems. In this paper, we extend the detailed analysis of structural and functional peculiarities of self-reproduction processes to unit cells of the Cooper-Helmstetter-Donachie cell cycle theory. We show that: 1. Our modelling framework allows to calculate physiological parameters (numbers of cell components, flux patterns, cellular composition, etc.) of unit cells, including also unit cell mass that determines the DNA replication initiation conditions. 2. Unit cells might have additional cell (cushioning) components that are responsible not only for carrying out various special functions, but also for regulating cell size and stabilizing the growth of cells. 3. The optimal productivity of the synthesis of cushioning components (useful cellular load) is observed at doubling time approximately two times longer than the minimal doubling time of the unit cells.

Regulatory elements coordinating initiation of chromosome replication to the Escherichia coli cell cycle

Escherichia coli coordinates replication and division cycles by initiating replication at a narrow range of cell sizes. By tracking replisomes in individual cells through thousands of division cycles in wild-type and mutant strains, we were able to compare the relative importance of previously described control systems. We found that accurate triggering of initiation does not require synthesis of new DnaA. The initiation size increased only marginally as DnaA was diluted by growth after dnaA expression had been turned off. This suggests that the conversion of DnaA between its active ATP- and inactive ADP-bound states is more important for initiation size control than the total free concentration of DnaA. In addition, we found that the known ATP/ADP converters DARS and datA compensate for each other, although the removal of them makes the initiation size more sensitive to the concentration of DnaA. Only disruption of the regulatory inactivation of DnaA mechanism had a radical impact on replication initiation. This result was corroborated by the finding that termination of one round of replication correlates with the next initiation at intermediate growth rates, as would be the case if RIDA-mediated conversion from DnaA-ATP to DnaA-ADP abruptly stops at termination and DnaA-ATP starts accumulating.

Beyond the average: An updated framework for understanding the relationship between cell growth, DNA replication, and division in a bacterial system

Our understanding of the bacterial cell cycle is framed largely by population-based experiments that focus on the behavior of idealized average cells. Most famously, the contributions of Cooper and Helmstetter help to contextualize the phenomenon of overlapping replication cycles observed in rapidly growing bacteria. Despite the undeniable value of these approaches, their necessary reliance on the behavior of idealized average cells masks the stochasticity inherent in single-cell growth and physiology and limits their mechanistic value. To bridge this gap, we propose an updated and agnostic framework, informed by extant single-cell data, that quantitatively accounts for stochastic variations in single-cell dynamics and the impact of medium composition on cell growth and cell cycle progression. In this framework, stochastic timers sensitive to medium composition impact the relationship between cell cycle events, accounting for observed differences in the relationship between cell cycle events in slow- and fast-growing cells. We conclude with a roadmap for potential application of this framework to longstanding open questions in the bacterial cell cycle field.

Robust replication initiation from coupled homeostatic mechanisms

The bacterium Escherichia coli initiates replication once per cell cycle at a precise volume per origin and adds an on average constant volume between successive initiation events, independent of the initiation size. Yet, a molecular model that can explain these observations has been lacking. Experiments indicate that E. coli controls replication initiation via titration and activation of the initiator protein DnaA. Here, we study by mathematical modelling how these two mechanisms interact to generate robust replication-initiation cycles. We first show that a mechanism solely based on titration generates stable replication cycles at low growth rates, but inevitably causes premature reinitiation events at higher growth rates. In this regime, the DnaA activation switch becomes essential for stable replication initiation. Conversely, while the activation switch alone yields robust rhythms at high growth rates, titration can strongly enhance the stability of the switch at low growth rates. Our analysis thus predicts that both mechanisms together drive robust replication cycles at all growth rates. In addition, it reveals how an origin-density sensor yields adder correlations.

Exponential trajectories, cell size fluctuations, and the adder property in bacteria follow from simple chemical dynamics and division control

Experiments on steady-state bacterial cultures have uncovered several quantitative regularities at the system level. These include, first, the exponential growth of cell size with time and the balanced growth of intracellular chemicals between cell birth and division, which are puzzling given the nonlinear and decentralized chemical dynamics in the cell. We model a cell as a set of chemical populations undergoing nonlinear mass action kinetics in a container whose volume is a linear function of the chemical populations. This turns out to be a special class of dynamical systems that generically has attractors in which all populations grow exponentially with time at the same rate. This explains exponential balanced growth of bacterial cells without invoking any regulatory mechanisms and suggests that this could be a robust property of protocells as well. Second, we consider the hypothesis that cells commit themselves to division when a certain internal chemical population reaches a threshold of N molecules. We show that this hypothesis leads to a simple explanation of some of the variability observed across cells in a bacterial culture. In particular, it reproduces the adder property of cell size fluctuations observed recently in E. coli; the observed correlations among interdivision time, birth volume, and added volume in a generation; and the observed scale of the fluctuations (CV ≈ 10-30%) when N is between 10 and 100. Third, upon including a suitable regulatory mechanism that optimizes the growth rate of the cell, the model reproduces the observed bacterial growth laws including the dependence of the growth rate and ribosomal protein fraction on the medium. Thus, the models provide a framework for unifying diverse aspects of bacterial growth physiology under one roof. They also suggest new questions for experimental and theoretical enquiry.

Cell Size Is Coordinated with Cell Cycle by Regulating Initiator Protein DnaA in E. coli

Sixty years ago, bacterial cell size was found to be an exponential function of growth rate. Fifty years ago, a more general relationship was proposed, in which cell mass was equal to the initiation mass multiplied by 2 to the power of the ratio of the total time of C and D periods to the doubling time. This relationship has recently been experimentally confirmed by perturbing doubling time, C period, D period, or initiation mass. However, the underlying molecular mechanism remains unclear. Here, we developed a theoretical model for initiator protein DnaA mediating DNA replication initiation in Escherichia coli. We introduced an initiation probability function for competitive binding of DnaA-ATP and DnaA-ADP at oriC. We established a kinetic description of regulatory processes (e.g., expression regulation, titration, inactivation, and reactivation) of DnaA. Cell size as a spatial constraint also participates in the regulation of DnaA. By simulating DnaA kinetics, we obtained a regular DnaA oscillation coordinated with cell cycle and a converged cell size that matches replication initiation frequency to the growth rate. The relationship between the simulated cell size and growth rate, C period, D period, or initiation mass reproduces experimental results. The model also predicts how DnaA number and initiation mass vary with perturbation parameters, comparable with experimental data. The results suggest that 1) when growth rate, C period, or D period changes, the regulation of DnaA determines the invariance of initiation mass; 2) ppGpp inhibition of replication initiation may be important for the growth rate independence of initiation mass because three possible mechanisms therein produce different DnaA dynamics, which is experimentally verifiable; and 3) perturbation of some DnaA regulatory process causes a changing initiation mass or even an abnormal cell cycle. This study may provide clues for concerted control of cell size and cell cycle in synthetic biology.

Low Affinity DnaA-ATP Recognition Sites in E. coli oriC Make Non-equivalent and Growth Rate-Dependent Contributions to the Regulated Timing of Chromosome Replication

Although the mechanisms that precisely time initiation of chromosome replication in bacteria remain unclear, most clock models are based on accumulation of the active initiator protein, DnaA-ATP. During each cell division cycle, sufficient DnaA-ATP must become available to interact with a distinct set of low affinity recognition sites in the unique chromosomal replication origin, oriC, and assemble the pre-replicative complex (orisome) that unwinds origin DNA and helps load the replicative helicase. The low affinity oriC-DnaA-ATP interactions are required for the orisome's mechanical functions, and may also play a role in timing of new rounds of DNA synthesis. To further examine this possibility, we constructed chromosomal oriCs with equal preference for DnaA-ADP or DnaA-ATP at one or more low affinity recognition sites, thereby lowering the DnaA-ATP requirement for orisome assembly, and measured the effect of the mutations on cell cycle timing of DNA synthesis. Under slow growth conditions, mutation of any one of the six low affinity DnaA-ATP sites in chromosomal oriC resulted in initiation earlier in the cell cycle, but the shift was not equivalent for every recognition site. Mutation of τ2 caused a greater change in initiation age, suggesting its occupation by DnaA-ATP is a temporal bottleneck during orisome assembly. In contrast, during rapid growth, all origins with a single mutated site displayed wild-type initiation timing. Based on these observations, we propose that E. coli uses two different, DnaA-ATP-dependent initiation timing mechanisms; a slow growth timer that is directly coupled to individual site occupation, and a fast growth timer comprising DnaA-ATP and additional factors that regulate DnaA access to oriC. Analysis of origins with paired mutated sites suggests that Fis is an important component of the fast growth timing mechanism.

Stochastic timing in gene expression for simple regulatory strategies

Timing is essential for many cellular processes, from cellular responses to external stimuli to the cell cycle and circadian clocks. Many of these processes are based on gene expression. For example, an activated gene may be required to reach in a precise time a threshold level of expression that triggers a specific downstream process. However, gene expression is subject to stochastic fluctuations, naturally inducing an uncertainty in this threshold-crossing time with potential consequences on biological functions and phenotypes. Here, we consider such ‘timing fluctuations’ and we ask how they can be controlled. Our analytical estimates and simulations show that, for an induced gene, timing variability is minimal if the threshold level of expression is approximately half of the steady-state level. Timing fluctuations can be reduced by increasing the transcription rate, while they are insensitive to the translation rate. In presence of self-regulatory strategies, we show that self-repression reduces timing noise for threshold levels that have to be reached quickly, while self-activation is optimal at long times. These results lay a framework for understanding stochasticity of endogenous systems such as the cell cycle, as well as for the design of synthetic trigger circuits.

Mathematical modeling of bacterial cell cycle: the problem of coordinating genome replication with cell growth

In this paper, we perform an analysis of bacterial cell-cycle models implementing different strategies to coordinately regulate genome replication and cell growth dynamics. It has been shown that the problem of coupling these processes does not depend directly on the dynamics of cell volume expansion, but does depend on the type of cell growth law. Our analysis has distinguished two types of cell growth laws, "exponential" and "linear", each of which may include both exponential and linear patterns of cell growth. If a cell grows following a law of the "exponential" type, including the exponential V(t) = V(0) exp (kt) and linear V(t) = V(0)(1 + kt) dynamic patterns, then the cell encounters the problem of coupling growth rates and replication. It has been demonstrated that to solve the problem, it is sufficient for a cell to have a repressor mechanism to regulate DNA replication initiation. For a cell expanding its volume by a law of the "linear" type, including exponential V(t) = V(0) + V(1) exp (kt) and linear V(t) = V(0) + kt dynamic patterns, the problem of coupling growth rates and replication does not exist. In other words, in the context of the coupling problem, a repressor mechanism to regulate DNA replication, and cell growth laws of the "linear" type displays the attributes of universality. The repressor-type mechanism allows a cell to follow any growth dynamic pattern, while the "linear" type growth law allows a cell to use any mechanism to regulate DNA replication.

Concerted control of Escherichia coli cell division

The coordination of cell growth and division is a long-standing problem in biology. Focusing on Escherichia coli in steady growth, we quantify cell division control using a stochastic model, by inferring the division rate as a function of the observable parameters from large empirical datasets of dividing cells. We find that (i) cells have mechanisms to control their size, (ii) size control is effected by changes in the doubling time, rather than in the single-cell elongation rate, (iii) the division rate increases steeply with cell size for small cells, and saturates for larger cells. Importantly, (iv) the current size is not the only variable controlling cell division, but the time spent in the cell cycle appears to play a role, and (v) common tests of cell size control may fail when such concerted control is in place. Our analysis illustrates the mechanisms of cell division control in E. coli. The phenomenological framework presented is sufficiently general to be widely applicable and opens the way for rigorous tests of molecular cell-cycle models.

Temperature-dependence of the DnaA-DNA interaction and its effect on the autoregulation of dnaA expression

The DnaA protein is a key factor for the regulation of the timing and synchrony of initiation of bacterial DNA replication. The transcription of the dnaA gene in Escherichia coli is regulated by two promoters, dnaAP1 and dnaAP2. The region between these two promoters contains several DnaA-binding sites that have been shown to play an important role in the negative auto-regulation of dnaA expression. The results obtained in the present study using an in vitro and in vivo quantitative analysis of the effect of mutations to the high-affinity DnaA sites reveal an additional effect of positive autoregulation. We investigated the role of transcription autoregulation in the change of dnaA expression as a function of temperature. While negative auto-regulation is lost at dnaAP1, the effects of both positive and negative autoregulation are maintained at the dnaAP2 promoter upon lowering the growth temperature. These observations can be explained by the results obtained in vitro showing a difference in the temperature-dependence of DnaA-ATP binding to its high- and low-affinity sites, resulting in a decrease in DnaA-ATP oligomerization at lower temperatures. The results of the present study underline the importance of the role for autoregulation of gene expression in the cellular adaptation to different growth temperatures.

Growth-rate-dependent dynamics of a bacterial genetic oscillator

Gene networks exhibiting oscillatory dynamics are widespread in biology. The minimal regulatory designs giving rise to oscillations have been implemented synthetically and studied by mathematical modeling. However, most of the available analyses generally neglect the coupling of regulatory circuits with the cellular "chassis" in which the circuits are embedded. For example, the intracellular macromolecular composition of fast-growing bacteria changes with growth rate. As a consequence, important parameters of gene expression, such as ribosome concentration or cell volume, are growth-rate dependent, ultimately coupling the dynamics of genetic circuits with cell physiology. This work addresses the effects of growth rate on the dynamics of a paradigmatic example of genetic oscillator, the repressilator. Making use of empirical growth-rate dependencies of parameters in bacteria, we show that the repressilator dynamics can switch between oscillations and convergence to a fixed point depending on the cellular state of growth, and thus on the nutrients it is fed. The physical support of the circuit (type of plasmid or gene positions on the chromosome) also plays an important role in determining the oscillation stability and the growth-rate dependence of period and amplitude. This analysis has potential application in the field of synthetic biology, and suggests that the coupling between endogenous genetic oscillators and cell physiology can have substantial consequences for their functionality.