# A bacterial size law revealed by a coarse-grained model of cell physiology

François Bertaux, Julius von Kügelgen, Samuel Marguerat, Vahid Shahrezaei

► To cite this version:

François Bertaux, Julius von Kügelgen, Samuel Marguerat, Vahid Shahrezaei. A bacterial size law revealed by a coarse-grained model of cell physiology. PLoS Computational Biology, 2020, 16 (9), 10.1371/journal.pcbi.1008245. [https://doi.org/10.1371/journal.pcbi.1008245](https://doi.org/10.1371/journal.pcbi.1008245).

HAL Id: hal-03152985

https://inria.hal.science/hal-03152985

Submitted on 18 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# A bacterial size law revealed by a coarse-grained model of cell physiology

François Bertaux¹,²,³,⁴, Julius von Kügelgen¹, Samuel Marguerat²,³*, Vahid Shahrezaei¹*

¹ Department of Mathematics, Imperial College London, London, United Kingdom, ² MRC London Institute of Medical Sciences (LMS), London, United Kingdom, ³ Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, United Kingdom, ⁴ Institut Pasteur, USR 3756 IP CNRS, Paris, France

* samuel.marguerat@imperial.ac.uk (SM); vshahrezaei@imperial.ac.uk (VS)

<figure><img src="image_2.png" /></figure>

## Abstract

Universal observations in Biology are sometimes described as “laws”. In *E. coli*, experimental studies performed over the past six decades have revealed major growth laws relating ribosomal mass fraction and cell size to the growth rate. Because they formalize complex emerging principles in biology, growth laws have been instrumental in shaping our understanding of bacterial physiology. Here, we discovered a novel size law that connects cell size to the inverse of the metabolic proteome mass fraction and the active fraction of ribosomes. We used a simple whole-cell coarse-grained model of cell physiology that combines the proteome allocation theory and the structural model of cell division. This integrated model captures all available experimental data connecting the cell proteome composition, ribosome activity, division size and growth rate in response to nutrient quality, antibiotic treatment and increased protein burden. Finally, a stochastic extension of the model explains non-trivial correlations observed in single cell experiments including the adder principle. This work provides a simple and robust theoretical framework for studying the fundamental principles of cell size determination in unicellular organisms.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Bertaux F, von Kügelgen J, Marguerat S, Shahrezaei V (2020) A bacterial size law revealed by a coarse-grained model of cell physiology. PLoS Comput Biol 16(9): e1008245. https://doi.org/10.1371/journal.pcbi.1008245

**Editor:** James R. Faeder, University of Pittsburgh, UNITED STATES

**Received:** May 18, 2020

**Accepted:** August 13, 2020

**Published:** September 28, 2020

**Copyright:** © 2020 Bertaux et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** Leverhulme Trust:Samuel Marguerat, Vahid Shahrezaei RPG-2014-408; UK Research and Innovation | Engineering and Physical Sciences Research Council (EPSRC):Vahid Shahrezaei EP/ N014529/1. Medical Research Council (MRC): Samuel Marguerat.

**Competing interests:** The authors declare that they have no conflict of interest.

## Author summary

Bacteria respond to environmental changes by adjusting their molecular composition, cell size and growth rate. This plasticity is thought to result from years of evolution and to be at least in part optimal for bacterial physiology. Over the past decades, quantitative studies of bacterial growth have revealed simple phenomenological relationships, called “growth laws”, which link cell size and cell composition to the growth rate. Simplified mathematical models of cell physiology are useful tools to gain quantitative understanding of the molecular mechanisms that underlie growth laws. For instance, these models helped explaining how optimal allocation of cellular resource to physiological processes and pathways governs the cell molecular composition in response to specific environmental conditions. In this study, we have extended and integrated existing mathematical models and used experimental data from several recent studies to understand the co-regulation of cell
composition, cell size and the cellular growth rate. The model predictions uncovered a
novel “size law” that links cell size to the levels of metabolic proteins and the fraction of
active ribosomes present in the cell. This work provides a useful theoretical tool and a
quantitative basis for understanding mechanistically bacterial physiology as a function of
external conditions.

**Introduction**

The behavior of complex biological systems can be described by surprisingly simple rules or “laws”, which connect quantitatively aspects of the cell composition with its physiology. Over the last six decades, the discovery of several rules, called the bacterial growth laws, transformed the field of microbial physiology [1–8]. Growth laws describe relationships between the exponential growth rate and cellular parameters. One such law states that the proportion of the cellular content dedicated to ribosomes is linearly dependent on the growth rate (we call this hereafter the first growth law). A second growth law states that average cell size is an exponential function of the growth rate (referred to hereafter as the second growth law). These laws formalize emerging relations between cellular processes quantitatively and therefore generate testable hypothesis for mechanistic studies of the molecular circuitry that underpin cell physiology [9,10].

What are the biological principles underlying the first growth law? The law was first ratio-
nalized as a direct consequence of the catalytic role of ribosomes in protein synthesis [11].
Indeed, if every ribosome in a cell synthesises proteins at a constant condition-independent
rate, then the steady-state growth rate will be proportional to the fraction of cellular content
devoted to ribosomes. However, this does not explain what sets the amount of resources
invested into ribosomes in a given environment, and what prevents the cell to increase this
investment in order to grow faster. Theoretical and experimental studies identified cellular
resource allocation as a key concept to address these questions [12]. The theory of resource
allocation describes how different groups of proteins with similar function (proteome sectors)
are regulated reciprocally in conditions that affect growth rate. In particular, these studies
showed that cells invest higher amounts into ribosomes when nutrient quality increases
because the investment in metabolic proteins required to achieve a given growth rate is
reduced. This explains the positive correlation of the proteome fraction dedicated to ribosomal
proteins with the growth rate and provides a model for how the fraction of ribosomes is set.
Interestingly, Hwa and co-workers [13] also showed that the first growth law is not always
valid when growth rate is not modulated by nutrient quality. While the law remains valid
when useless proteins are over-expressed, it breaks down when translation is inhibited. Based
on these findings, they proposed a phenomenological model of proteome allocation whose
predictions go beyond the first growth law because it captures these orthogonal types of
growth rate modulation [13]. More mechanistic coarse-grained models predicting both the
growth rate and the coarse-grained proteome as a function of growth conditions have also
been proposed confirming and extending these findings [14–25].

In contrast, the mechanisms underlying the second growth law, which connects growth
rate and cell size, are less well understood. In 1968, Donachie [26] proposed that, in bacteria,
an exponential dependency of cell size on growth rate indicates that: 1) the timing of DNA rep-
lication is fully controlled by cell size, and is triggered at a constant volume per origin of repli-
cation; 2) cell division is enslaved to the DNA replication processes, and occurs after a
constant time-interval following replication initiation, which encompasses DNA replication
and cell division (C+D period). However, while recent work has confirmed that DNA replica-
tion is strongly size-controlled [27] and that the initiation size per origin is invariant across a
large range of growth rates and types of growth rate modulations [28], what determines the
duration of the C+D period is still unclear. Moreover, the relationship between cell size and
growth rate does not follow the second growth law, when useless proteins are over-expressed
on gene expression is increased or translation is inhibited [28,29]. This indicates that duration
of the C+D period is regulated in non-trivial ways outside of the canonical growth modula-
tions based on nutrient quality. To date, unlike for the first growth law, mechanistic under-
standing of this phenomenon in connection to the second growth law is lacking.

The two growth laws discussed so far are based on average cell population measurements. However, the variations at the single cell level from the population average also provides rich and complementary insights into the biology of cell growth and division [30,31]. Cells can control their size through a ‘sizer’ mechanism, where regardless of birth size they grow to an average size and then they divide. But, recent advances in single-cell techniques led to the discovery that bacterial cells achieve cell size homeostasis via an ‘adder’ principle, where cells add a constant size between birth and division independent of the birth size [32,30,33]. The adder phenomenon has been observed in several bacteria, yeast, archea and mammalian cells [34,30,35–42]. Despite its broad conservation, the mechanistic basis of this phenomenon remains elusive and hotly debated. It has been suggested that the adder phenomenon is related to second growth law assuming that replication controls cell division [43,44]. However, this model is challenged by recent studies that questioned the role of DNA replication in controlling division of single cells [45–47].

Here, we present a coarse-grained model of bacterial physiology that unifies the proteome allocation theory with the structural model of division control [29,30,48]. Building on existing works we have developed a simple, quantitative and actionable model able to integrate processes which so far had only been modelled in isolation. Specifically, this model relies on a minimal number of fitting parameters to capture experimental measurements of proteome allocation, ribosome activity and cell size across several growth modulations. It allowed to formalise a new bacterial “size law” which states that cell size depends on the metabolism proteome sector mass fraction and the fraction of active ribosomes. This law is simple and more general than the original second growth law. Finally, the model exhibits single cell properties that are consistent with experimental observations including the adder size control.

Results

A minimalistic whole-cell coarse-grained model to predict cell composition, growth rate and cell size

We developed a minimalistic whole-cell coarse-grained mathematical model to predict cell composition, growth rate and cell size for at least three types of growth rate modulation: 1) change of nutrient quality, 2) translation inhibition by chloramphenicol and 3) over-expression of useless proteins. These modulations were selected because they have been instrumental in uncovering key principles of proteome allocation underpinning the first growth law [13] and have abundant experimental measurements available [13,29,30,49,28]. Our model extends the existing whole-cell coarse-grained models [14–25] by including for the first time both cell composition and cell size. In addition, our model captures a wide-range of experimental measurements including cell size [29–30], proteome allocation [13,49] and ribosome activity [49] data.

As in previous work, the model (Fig 1, see Methods for a mathematical description) considers only two types of molecular components: protein precursors (noted A) and proteins
<figure><img src="image_4.png" /><figcaption>**Fig 1.** A simple whole-cell coarse-grained model of bacterial growth reproduces proteome allocation and ribosome activity data. (A) Schematic of model structure displaying model components (proteome sectors E, R, Q and X and protein precursor A) and reactions (metabolism, protein synthesis). How the three types of growth rate modulation (nutrient quality, expression of useless protein and ribosome inactivation) are modeled is also highlighted. The rate of precursor creation is proportional to the amount of metabolic enzymes E. The total protein synthesis rate is proportional to the number of active ribosomes, and the synthesis rate per ribosome (the ribosome efficiency) is dependent on the precursor concentration via a Michaelis-Menten relationship. A full description of the model is given in the Methods section. (B) Allocation of total protein synthesis capacity between proteome sectors. A fixed, condition-independent fraction $f_Q$ is allocated to housekeeping proteins. Expression of useless protein imposes a fixed allocation $f_U$. The allocation to R proteins is proportional to precursor concentration, and the remainder is allocated to E proteins (assuming that $f_X \ll 1$). (C-F) Model predictions (solid lines) agree with experimental data (circles-Scott et al. 2010 [13], squares-Dai et al. 2017 [49]). Parameters were all fixed from published data (see Methods). (C) Ribosomal proteome fraction data for nutrient and chloramphenicol growth rate modulations. Colors indicate different nutrient qualities. (D) Relationship between growth rate and useless proteome fraction for forced expression of useless protein. (E-F) Ribosome efficiency ($\frac{a}{a+d_{sat}}$) and active ribosome fraction ($\frac{R_a}{R}$).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1008245.g001

[15,16,20]. Proteins are split between four different classes (or sectors): transport and metabolism proteins (referred to as only metabolic sector for brevity) (E), ribosomal proteins (R), housekeeping proteins (Q), and division proteins (X) [29,30]. E proteins catalyse the import and transformation of external nutrients (not explicitly modelled) into precursors A. R proteins are involved in the synthesis of proteins from precursors. Q proteins are involved in cellular functions that do not directly contribute to growth. X proteins are regulating cell division (see below). Lastly, when growth rate is modulated by over-expression of useless proteins, we consider an additional class of “useless” proteins noted U (Fig 1). All proteins are stable in our model and only diluted by cell growth. The model enforces mass conservation—cell mass being defined as the sum of all molecular components. Finally, except for the import and transformation of nutrients into precursors, all reactions leave cell mass unchanged and the model assumes mass density homeostasis (cell volume and cell mass are therefore equivalent, and so are mass fractions and concentrations), as observed experimentally [29].

To fully specify the model, one should describe how the protein synthesis allocation fractions (how much of the total protein synthesis capacity is devoted to each protein class) are set
(Fig 1B). A fixed, condition-independent fraction $f_Q$ is allocated to Q proteins. We assume that the allocation fraction to division proteins is small ($f_X \ll 1$). The allocation fraction $f_U$ is imposed and represents the level of useless protein expression. The main allocation decision is therefore between $f_R$ and $f_E$. In this study, rather than assuming growth rate maximization to set these fractions, we assume that $f_R$ (and hence $f_E$) is dynamically regulated as detailed below. This choice is motivated by the recent observation [49] that the ribosome translation elongation rate exhibits a Michaelis-Menten dependence on the total ribosome proteome fraction across nutrient and chloramphenicol-mediated growth modulations. This implies that addition of chloramphenicol to cells growing in poor media results in an increase in the translation elongation rate. This observation is somewhat unintuitive and had not been captured by previous modelling studies [50,51]. Interestingly, this dependency can be simply reproduced in our model by assuming that the allocation fraction to ribosomal proteins $f_R$ is proportional to the precursor concentration $a$, $f_R = \delta a$ (Methods). We further assume that $\delta$ is large enough to ensure that $a \ll 1$, as the steady-state concentration of free amino-acids is known to remain small [10]. Simulations show that this approximation is already valid when $\delta = 5$ (Fig A in S1 Text). Then, the steady-state cell composition and growth rate can be predicted and depend only on: 1) four condition-independent parameters: the maximal ribosome speed $\sigma$, the ribosome regulation constant $K$, the housekeeping protein allocation fraction $f_Q$ and the unbinding constant $k_{off}^{cm}$ of chloramphenicol-ribosome complexes; and 2) three growth modulation parameters: the medium nutrient quality $k$, the chloramphenicol-imposed inactivation rate of ribosomes $k_{on}^{cm}$, and the useless protein allocation fraction $f_U$.

The model reproduces quantitatively proteome allocation data [13,49] and ribosome activity data [49] for the three types of growth modulation (Fig 1C–1F). This agreement is obtained without parameter fitting, because the parameters $\sigma$, $K$, $f_Q$ and $k_{off}^{cm}$ are directly constrained from experimental measurements [13,49]. To our knowledge, this is the first coarse-grained model of bacterial growth to reproduce both proteome allocation and ribosome activity data.

While it is often assumed that proteome fractions are optimally allocated to maximize growth rate [15], in our model proteome fractions are dynamically set by the simple regulation $f_R = \delta a$. Interestingly, this regulation achieves near-optimal growth rates across the three types of growth rate modulations considered in this study (Fig B in S1 Text). This is consistent with previous modelling work that shows at fast growth bacterial allocation of ribosome fractions is tuned to approximate growth-rate maximisation [18]. However, we note that the model assuming optimal proteome allocation does not reproduce the ribosome activity data of Fig 1E.

Our model also explicitly includes cell division and assumes that it is triggered by the accumulation of proteins belonging to the X sector to an absolute copy number threshold $X_{div}$ (Fig 2A), as proposed before in the structural model [48,30,29]. This design enables prediction of changes in both cell size and coarse-grained cell composition as a function of cell growth. Note that for simplicity and because it appears unnecessary, we do not assume that X proteins are destroyed immediately after division, in contrast to other studies [30,48]. Altogether, our model assumptions result in simple and intuitive equations characterizing model steady-states (Fig 2B). Importantly, steady-state concentrations of cell components are not dependent on the cell cycle as expected under balanced growth and are valid when cell division is not modelled (as done for the fitting of the proteome fraction data above).

In summary, we have developed a whole-cell model of cell physiology that captures changes in proteome fractions and ribosome activity as a function of growth conditions while uniquely incorporating direct modeling of cell division control.
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig 2.** Integration of the structural model enables prediction of both cell composition and cell size. (A) Simulation of model dynamics. After cell division (triggered when *X* number reaches *X*<sub>div</sub>), cell content is partitioned equally among daughter cells and only one is followed ('mother machine' setting, Wang et al., 2010). After a transient adaptation period, concentrations of cellular components are constant, cellular growth (size increase) is exponential and division occurs at a constant size. Here, nutrient quality is such that the growth rate is 1 hr<sup>-1</sup> and there is no chloramphenicol nor useless expression. (B) Equations characterizing model steady-state together with their intuitive interpretation (see Methods for their derivation). *r*<sub>a</sub> denotes the concentration of active ribosomes.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1008245.g002

Cell size can be predicted from coarse-grained cell composition across all conditions revealing a novel size law

Experimental data from previous studies shows a complex relationship between cell size and growth rate across the three types of growth modulation (Fig 3B, Fig C in S1 Text). While cell size increases with growth rate when nutrient quality is varied as posited by the second growth law, it decreases with the growth rate in response to overexpression of useless proteins. On the other hand, chloramphenicol-mediated translation inhibition leads to cell size increase or decrease depending on the nutrient quality of the medium. Basan and colleagues [29] proposed that a structural model of cell division requires the expression of the division protein *X* to follow the inverse of cell size assuming an invariant division threshold *X*<sub>div</sub> (Fig 2B, bottom). We therefore asked whether these complex variations in cell size with the growth rate could be explained by the regulation of division proteins *X* by different coarse-grained proteome sectors via the allocation fraction *f*<sub>*X*</sub> (Fig 3A).

We first considered a model whereby $f_X$ follows the concentration of a single coarse-grained proteome sector. Interestingly, size measurements for both nutrient quality and useless protein expression growth rate modulations are explained if $f_X$ is set to follow the concentration (noted *e*, not to confuse with the mathematical constant) of the metabolic sector *E* exclusively (Fig 3D). However, critically, this simple dependency cannot explain chloramphenicol-mediated translation inhibition data (Fig D in S1 Text). We therefore extended this analysis to dependencies of $f_X$ on two coarse-grained proteome quantities and tested all pairwise combinations of the quantities *e*, *r*, $r_a$ and $r_r^a$ (Fig E in S1 Text). We found that when $f_X$ depends on both the *E* sector concentration *e* and the fraction of active ribosomes $r_r^a$, the observed relationship between size and growth rate is quantitatively captured for all types of growth rate modulation, including translation inhibition (Fig 3C and 3E). Interestingly, the exponent for the *E* sector concentration remained close to 1 when fitting two quantities instead of one and
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Regulation of division proteins by two proteome sectors quantitatively explain cell size across growth modulations. (A) Hypothesis stating that X expression could depend on the concentration of one or two coarse-grained proteome sectors. Here we assume that $X_{div}$ is invariant. (B) Empirical relationship between cell size and growth rate for three types of growth rate modulation (nutrient quality, chloramphenicol-mediated translation inhibition and expression of useless protein). Data aggregated from three studies, Basan et al. (2015), Si et al. (2017) and Taheri-Araghi et al. (2015). A scaling factor for size was applied on Si et al. and Taheri-Araghi et al. data to make the nutrient modulation data of the three studies consistent (see Fig C in S1 Text). Different branches for useless expression and chloramphenicol indicate modulations at different nutrient quality. (C) Predicted relationship between cell size and growth rate when X expression depends on both E concentration and the fraction of active ribosomes. (D-F) Log-log plots comparing model predictions with experimental data. The natural log is used. (D) Regulation of X expression by E alone can explain size data for nutrient modulation and useless expression modulation (but not chloramphenicol-mediated translation inhibition, Fig D in S1 Text). (E) Among possible regulations of X by two coarse-grained quantities, regulation by E concentration and the fraction of active ribosomes explains size data for all three types of growth rate modulations (this regulation model was used for (C)). (F) Predicted C+D durations agree with measurements by Si et al. (2017). C+D is predicted by $C + D = \frac{\log[V_{div}/S_0]}{\alpha}$, where $\alpha$ is the growth rate, $V_{div}$ is the model-predicted size (3C and 3E) and $S_0$ is a constant measured in Si et al. (2017).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1008245.g003

including chloramphenicol modulation (Fig 3D and 3E). This validates a consistent contribution of this sector across the three growth modulations. We note that the fitted exponent values are very close to 1 and -2/3 for e and $r_a/r$ respectively. Because $f_X$ and cell size are inversely proportional, it then follows that cell size V scales with $e^{-1} \times (r_a/r)^{2/3}$. Based on this we propose a novel size law formalized as:

$$V \propto e^{-1} \times (\frac{r_a}{r})^{2/3}$$

Our model assumes that in the absence of translation inhibition all ribosomes remain active ($r_a/r = 1$). However, it has been observed that for nutrient-limited slow growth the active ribosome fraction is reduced [49]. Interestingly, using these experimental data instead of the model predictions improves the accuracy of the size law at slow growth rates (Fig F in S1 Text), providing further independent support for the proposed relationship between size and
active ribosome fraction. The duration of the C+D period has been observed to have non-trivial variations (Si et al., 2017). We note that this size law can predict the duration of the C+D period (Fig 3F) when assuming that DNA replication occurs at a fixed condition-independent cell size (i.e. using the relationship $V_{\text{div}} = S_0 \exp(\alpha(C+D))$, where $\alpha$ is the growth rate and $S_0$ is the constant initiation size).

In summary, our minimalistic model quantitatively captures how both size and coarse-grained cell composition change with the growth rate for all three types of growth modulation with only two parameters (the two exponents of the $f_X$ dependency on $e$ and $r_a/r$, those exponents being very close to 1 and -2/3 respectively). Therefore, *E. coli* cell size can be predicted from only two coarse-grained proteome quantities, metabolism sector concentration and the fraction of active ribosomes leading to the definition of a novel growth law.

## Emergence of 'adder' size homeostasis and cellular individuality in the presence of noise

So far, we have neglected cell-to-cell variability and focused on predicting average cell composition, size and growth rate at steady-state as a function of growth conditions. However, isogenic cells growing in a constant environment show significant phenotypic variability, notably for global traits such as growth rate or division size [30,52–54]. Moreover, non-trivial correlations between single-cell traits are often observed and contain rich information about regulatory mechanisms [55,33,56,30,31,57,23,25]. Since our model is based on biochemical reactions between coarse-grained molecular components, dynamic cell-to-cell variability naturally emerges when adopting a stochastic interpretation of all reactions and random partitioning of components at division (Fig 4A).

In the deterministic model, the two parameters $X_{\text{div}}$ (the number of division molecules $X$ needed to trigger cell division) and $f_X^{\text{scale}}$ (the constant factor in $f_X = f_X^{\text{scale}} \times e^{\left(\frac{r_a}{r}\right)^{-2/3}}$) are only setting the size scale via their ratio (Fig 3A). In the presence of reaction noise, those parameters will now impact the extent of cell-to-cell variability in growth rate and in size at division. Even with this simple model of gene expression noise, we could find values of those two parameters that lead to reasonable estimates of both noise in growth rate and size at birth (Methods).

Stochastic simulations of our model lead to a near 'adder' size homeostasis, as observed experimentally (Fig 4B, red circles). Interestingly, the deviation towards a 'sizer' size homeostasis observed experimentally at very slow growth, [27] could also be explained by additional model assumptions (Fig G in S1 Text). These include the constant degradation of $X$ during the cell cycle, which could dominate dilution at slow growth (Fig 4B, blue squares). We then explored the relationship between the size at birth and the growth rate of single cells for several conditions across the three types of growth rate modulation (Fig 4C, triangles). This relationship at the level of single cells is known to deviate from the one connecting average size at birth and average growth rate when conditions are varied [30,31]. These deviations were qualitatively captured by our stochastic model (Fig 4C, compare grey lines and green triangle groups).

In summary, our minimalistic coarse-grained model generates realistic cell-to-cell variability in growth rate and cell size, as well as realistic correlations between added size and birth size ('adder' size homeostasis) and between growth rate and birth size.

## Discussion

Our aim in this work was to provide a mechanistic understanding of the variations in cell size across growth conditions recently observed in *E. coli*. To this end, we have proposed a
A

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

B

<figure><img src="image_10.png" /></figure>

C

<figure><img src="image_11.png" /></figure>

**Fig 4. Emergence of 'adder' size homeostasis and cellular individuality in the presence of reaction noise.** (A) Simulation of model dynamics with molecular noise (top: total cell size, bottom: concentrations of cellular components). At division, cellular components are randomly split between daughter cells, so the tracked daughter has a probability 1/2 of getting each mother cell component. Same parameters as in Fig 2A. The parameters $X_{div}$ and $f_X^{scale}$ were chosen to obtain realistic cell-to-cell variability in size at birth and growth rate (Methods). (B) Model leads to near 'adder' size homeostasis. Average added size during one cell cycle as a function of birth cell size (via binning) is plotted. The very weak deviation towards 'sizer' can be explained by a residual correlation between size at birth and X count at birth due to the non-zero share of X in total size. A model variant where X is actively degraded at a constant rate and displaying a stronger deviation towards 'sizer' behavior is also shown. Other model variants exploring size homeostasis properties are discussed in Fig G in S1 Text. (C) Emergence of cellular individuality. Stochastic stimulations are performed for ten different growth conditions: seven nutrient qualities (green triangle groups), two useless expression strengths (red triangle groups) and one chloramphenicol (blue triangle groups). For each condition, cell cycles are binned by growth rate and the corresponding birth size is plotted. Continuous lines show the prediction of the deterministic model for the three growth rate modulations. Grey lines indicate experimental trends for different nutrient qualities extracted from mother machine data (Taheri-Araghi et al., 2015).

https://doi.org/10.1371/journal.pcbi.1008245.g004

minimalistic whole-cell coarse-grained model of E. coli physiology to predict the relationship between growth rate, cell composition, ribosome activity and cell size for three different types of growth rate modulation (Fig 1A and Fig 2). Our model builds on and unifies previous efforts to understand proteome allocation [15,13,10,16,20] and division control [48,29,30,58]. While other theoretical models [13,10,16,20] have also captured proteome allocation data, our model stands out because: 1) it includes ribosome activity and captures more data; 2) it is holistic with respect to cell composition and cell size and does not a priori neglect the mass fraction of protein precursors; and 3) it is not based on optimization of growth rate but rather a positive regulation of ribosome synthesis by precursors. Moreover, despite its simplicity and low number of free parameters (only the exponents in the size law were fitted and obtained to be close to -1 and 2/3 respectively; the other four parameters were constrained by previous studies), our model quantitatively reproduces experimental data on proteome allocation (Fig 1C and 1D), ribosome activity (Fig 1E and 1F) and average cell size (Fig 3C and 3E) for the three types of
growth rate modulation (change of nutrient quality, chloramphenicol-mediated ribosome inactivation and expression of useless proteins). Our model therefore compares favourably to existing coarse-grained models in terms of number of parameters (for instance, the pioneering model by Weisse and colleagues [16] has more than 10 parameters and is only focused on data from Scott and colleagues [13]). In addition, when reaction noise is included in the model, experimental observations such as the *adder* principle of size homeostasis and non-trivial correlations between single-cell cellular growth rate and cell size are predicted (Fig 4).

A remarkable result of our study is the emergence of a novel size law related to the second growth law. This law states that only one coarse-grained quantity characterizing cell composition is sufficient to predict quantitatively cell size changes in response to nutrient quality and to over-expression of a useless protein (Fig 3C and 3E, Fig D in S1 Text). A second coarse-grained quantity is required for prediction of cell size when translation is inhibited. Those quantities are the mass fraction (or equivalently, concentration) of metabolic proteins (e in the model) and the fraction of total ribosomes that are active ($r_a/r$), respectively. Strikingly, best fit is robustly obtained for exponent values close to ~-1 for the metabolic protein concentration and ~2/3 for the fraction of active ribosomes. This results in a 'size law' stating:

$$V \propto e^{-1} \times \left(\frac{r_a}{r}\right)^{2/3}.$$ We did not impose constraints on the values of exponents when fitting size data, the fact that fitted exponents are close to such particular values is therefore remarkable and suggests a fundamental underlying physical mechanism. Also, as reported in Fig F in S1 Text, the quality of the size predictions for nutrient modulation data at slow growth increases when the empirical active fractions of ribosomes are accounted for (below 1 at very slow growth, while our model assumes 1). This further validates the contribution of the active fraction of ribosomes to the law. In the absence of translation inhibition and at intermediate to fast growth the fraction of active ribosomes is close to 1 (Dai et al. 2016), therefore the law reduces to the simple form $V \propto e^{-1}$. This relationship has been observed before for modulation of growth in the absence of translation inhibition [30]. Basan and colleagues [29] also noted that if the growth rate dependence of X protein follows the one of constitutively expressed proteins for modulation by either nutrient quality or over-expression of useless proteins, then their size data could be explained. Because expression of the metabolic sector has the same growth rate dependence as constitutive expression upon nutrient modulation [13], our results agree with those previous observations. However, neither studies were able to explain cell size measurements upon translation inhibition with chloramphenicol. Basan and colleagues [29] opted for a conservative approach designed to uncover only one structural protein and which relied on condition-independent division thresholds. Our approach, on the other hand, is less restrictive and integrate proteome allocation with cell size data leveraging all predictions of our coarse-grained model, including different proteome sectors and ribosome activities. This allowed to uncover the role for ribosome activity in setting cell size and to capture experimental data from changing conditions including translation inhibition.

Our 'size law' is very different in nature from the second growth law stating that $V \propto \exp[\alpha(C+D)]$, where $\alpha$ is the growth rate, $C$ the duration of DNA replication and $D$ the duration between end of DNA replication and cell division. By construction, the second growth law links size at division to the control of DNA replication by size [26,27,44]. It predicts a simple exponential dependency with growth rate for nutrient quality modulation because those three quantities are largely independent of growth rate. However, while average size per origin at replication initiation is invariant across a wide range of growth conditions, both $C$ and $D$ durations change with growth conditions in complex ways [28]. Therefore, the second growth law is replication-initiation-centric and of limited use to make predictions on cell size at division outside of the canonical type of growth rate modulation (change in nutrient quality). In
contrast, the ‘size law’ which we propose here links coarse-grained cell composition to cell size
directly; therefore, by assuming the second growth law, we can also explain the observed varia-
tion in the C and D durations at least across three different types of growth rate modulation.
This law can generate directly testable predictions, such as cell size and C+D durations for
growth modulations involving both over-expression of useless proteins and translation inhibi-
tion by chloramphenicol.

In agreement with our results, recent studies contradict the replication-initiation-centric
view of cell size control [44]. Based on single-cell correlation statistics, Micali and colleagues
[46] propose that two concurrent processes, each controlling DNA replication and division, co-
determine the effective size at division. Furthermore, Si and colleagues [47] experimentally
demonstrated that two independent ‘adder’ homeostasis mechanisms are in place for DNA rep-
lication initiation and cell division. This does not mean that cellular constraints related to DNA
replication are not affecting average cell size at division. Yet it indicates that cell division is not
enslaved to the process controlling DNA replication initiation. Also, recent experiments
reported different amounts of DNA in newborn cells, suggesting the chromosome is not
involved in size homeostasis or responsible for the adder behavior [59]. Here we find that
homeostatic properties at the level of cell division appear consistent with the structural model of
division control. As shown in Fig 4, the structural model implemented at the single cell level
recovers ‘adder’ size control [58], deviation towards ‘sizer’ [27] in some parameter regimes (Fig
G in S1 Text) and non-trivial deviations from growth law at the single cell level [30,31]. We
note that as suggested by a recent model [46], there could be still a role for DNA-replication or
segregation in division control, for example when these processes are slowed down [28,35].

In our model, the quantitative relationship between cell composition and cell size is
explained via a dependence of the allocation fraction of division proteins *X* in *e* and r/a (Fig 3A),
but this relationship is also valid in itself as a phenomenological observation (Fig 3C). Thus, the
observed growth law is independent of the validity of the structural model and our whole-cell
coarse-grained model. In addition, although the total *X* protein fraction is likely to be negligible,
as it is commonly the case for cell cycle regulators, it should still follow the regulation of the
major proteome fractions in the cell. What could explain the quantitative relationship between
the coarse-grained quantities *e* and r/a and cell size? Under the structural model of division con-
trol, assuming division threshold *X*<sub>div</sub> is invariant across conditions (Fig 3A), it would imply
that the *X* allocation fraction *f*<sub>*X*</sub> scales with *e* and the inverse of (r/a<sup>2/3</sup>). In such case, the depen-
dency on the active ribosome fraction could for example reflect chloramphenicol-mediated
alteration of *X* expression, such as operon polarity via premature transcription termination
[60]. However, under the same assumption that *X*<sub>div</sub> is invariant, Basan and colleagues [29]
could not identify in proteomics data candidate proteins whose proteome fraction is behaving
quantitatively as you would expect for the factor *X*, that is as the inverse of their cell size data.
This result suggests that there may be no single protein triggering division when reaching a
fixed, condition-independent amount. An alternative possibility consistent with our results is
that *X* follows exactly the relative abundance of metabolic proteins *E*, while the division thresh-
old *X*<sub>div</sub> scales with the fraction of active ribosomes to the power 2/3. Other constraints could
explain a condition-dependent threshold *X*<sub>div</sub> such as the cell geometry [45]. A recent study
reported that rod-shaped bacteria like *E. coli* maintain an approximately constant aspect ratio
across various growth conditions, resulting in the cell surface area scaling as the cell volume to
the power 2/3 [61].

It is also possible that the resource allocation that maximizes growth in a given condition
imposes a ratio between surface and cytoplasmic proteins that constrains the surface-to-vol-
ume ratio (i.e. the cell width for rod-shaped organisms). In turn, depending on the mechanism
‘counting’ the absolute amount of X molecules, the threshold can depend on cell width. In the context of a width-dependent threshold, as suggested by several studies FtsZ is an attractive candidate as the sole X factor [62,45,63,47]. Another mechanism involving geometrical constraints is the recent observation that ribosomes are spatially organized via nucleoid exclusion and that nucleoid size scales with cell size [64].

In this work, we have tried to explain the mechanistic origin of size regulation across growth conditions. Whether this is an evolved trait that optimizes fitness is a more difficult question. While size changes due to expression of useless proteins or translation inhibition by drugs may not be an evolved trait, the increased size with nutrient modulation of growth rate seems to be a universal property of microbes, which makes it likely to be an evolved property of unicellular systems [65,66]. In a recent *in silico* study, we showed that increased size at fast growth prevents an increase in gene expression noise for proteins with decreasing concentrations [67]. Metabolic proteins *E* are such proteins, and noise in their expression is likely detrimental to fitness. Interestingly, the inverse relationship between cell size and metabolic protein concentration *e* found in this study implies that the total number of metabolic proteins is kept constant across growth conditions. This raises the possibility that the evolution of size regulation in response to growth conditions has been driven by selective pressures acting on metabolic protein expression noise rather than on size directly. Testing fitness benefit of size regulation is difficult [63,68] but should ultimately shed light on the constraints that shaped its evolution.

Modelling cellular systems by considering the cellular context within which they operate is fundamental to a full mechanistic understanding of their function. Modelling has shown that function of simple synthetic networks can be affected by cell cycle and cell growth, when cell physiology and resource allocation is taken into account [6,16,69,67]. The current study is a first attempt at integrating coarse-grained models of cell physiology with simple models of division control, which has revealed a novel size law providing a link between cell size and cellular resource allocation. Overall, our study highlights the importance of systematic integration of whole-cell physiological models in the study of natural and synthetic cellular systems.

# Methods

## Model reactions

Our model considers the following reactions between coarse-grained cell components:

$$ \emptyset \xrightarrow{k_E} A \text{ (nutrient import and transformation into protein precursors A)} $$

$$ A \xrightarrow{f_E \sigma R_{a+\frac{a}{dsal}}} E \text{ (synthesis of E proteins)} $$

$$ A \xrightarrow{f_R \sigma R_{a+\frac{a}{dsal}}} R_a \text{ (synthesis of active R proteins)} $$

$$ R_a \xrightarrow{k_{on}^{cm} R_a} R_i \text{ (inactivation of active R proteins via chloramphenicol binding)} $$

$$ R_i \xrightarrow{k_{off}^{cm} R_i} R_a \text{ (re-activation of chloramphenicol-inactivated R proteins via chloramphenicol unbinding)} $$

$$ A \xrightarrow{f_Q \sigma R_{a+\frac{a}{dsal}}} Q \text{ (synthesis of Q proteins)} $$

$$ A \xrightarrow{f_X \sigma R_{a+\frac{a}{dsal}}} X \text{ (synthesis of X proteins)} $$

$$ A \xrightarrow{f_U \sigma R_{a+\frac{a}{dsal}}} U \text{ (synthesis of U proteins)} $$

Where the cell volume is $V = E+R+Q+X+U+A$ (with $R = R_a+R_i$), the protein synthesis allocation parameters ($f_E+f_R+f_Q+f_X+f_U=1$, and $a = \frac{A}{V}$ (we will also note $e, r, q, x, u, \dots$ the concentrations $E/V, R/V, \dots$)). The parameter *k* is the medium nutrient quality, $k_{on}^{cm}$ is the chloramphenicol-imposed inactivation rate of ribosomes and $f_U$ is the imposed allocation fraction of useless proteins. Those three parameters represent the three types of growth
rate modulations that we consider. The parameter $\sigma$ is the maximal rate of protein synthesis by one ribosome and $a_{sat}$ is the enzymatic saturation constant of ribosomes. The parameter $k_{off}^{cm}$ is the dissociation rate of chloramphenicol-ribosome complexes. Note that all protein synthesis reactions conserve mass, so that volume growth equals nutrient import ($\frac{dV}{dt} = kE$ as seen below).

## Differential equations

A deterministic interpretation of these reactions leads to the following set of differential equations for absolute amounts of coarse-grained cell components (whose absolute mass can be evaluated to approximately 0.15 femtograms, given typical values of Xdiv and fX considered here and the typical mass of a cell):

$$ \frac{dA}{dt} = kE - \sigma R_a \frac{a}{a + a_{sat}} $$

$$ \frac{dE}{dt} = f_E \sigma R_a \frac{a}{a + a_{sat}} $$

$$ \frac{dR_a}{dt} = f_R \sigma R_a \frac{a}{a + a_{sat}} - k_{on}^{cm} R_a + k_{off}^{cm} R_i $$

$$ \frac{dQ}{dt} = f_Q \sigma R_a \frac{a}{a + a_{sat}} $$

$$ \frac{dX}{dt} = f_X \sigma R_a \frac{a}{a + a_{sat}} $$

$$ \frac{dU}{dt} = f_U \sigma R_a \frac{a}{a + a_{sat}} $$

$$ \frac{dR_i}{dt} = k_{on}^{cm} R_a - k_{off}^{cm} R_i $$

It can be seen that the right-hand side of those equations sum up to $kE$. The above differential equations can be also written in terms of concentrations, noticing that $\frac{dV}{dt} = kE = keV$:

$$ \frac{da}{dt} = ke - \sigma r_a \frac{a}{a + a_{sat}} - kea $$

$$ \frac{de}{dt} = f_E \sigma r_a \frac{a}{a + a_{sat}} - ke^2 $$

$$ \frac{dr_a}{dt} = f_R \sigma r_a \frac{a}{a + a_{sat}} - k_{on}^{cm} r_a + k_{off}^{cm} r_i - ker_a $$

$$ \frac{dq}{dt} = f_Q \sigma r_a \frac{a}{a + a_{sat}} - keq $$
$$ \frac{dx}{dt} = f_X \sigma r_a \frac{a}{a + a_{sat}} - kex $$

$$ \frac{du}{dt} = f_U \sigma r_a \frac{a}{a + a_{sat}} - keu $$

$$ \frac{dr_i}{dt} = f_R \sigma r_a \frac{a}{a + a_{sat}} + k_{on}^{cm} r_a - k_{off}^{cm} r_i - ker_i $$

In the above equations, the last terms represent dilution by cell growth (ke). Because concentrations are conserved at division, their dynamics and steady-state values are independent of the division threshold $X_{div}$.

## Steady-state growth

In this model, steady-state growth (or balanced growth, i.e. steady concentrations of coarse-grained cell components) is necessarily exponential at rate $\alpha = ke$, because $\frac{dV}{dt} = kE = keV$. By setting the time derivatives to zero, we obtain the following equations constraining the steady-state cell composition ($e, r_a, r_i, q, x, u, a$):

$$ ke(1 - a) = \sigma r_a \frac{a}{a + a_{sat}} $$

$$ ke^2 = f_E \sigma r_a \frac{a}{a + a_{sat}} $$

$$ ker_a = f_R \sigma r_a \frac{a}{a + a_{sat}} - k_{on}^{cm} r_a + k_{off}^{cm} r_i $$

$$ keq = f_Q \sigma r_a \frac{a}{a + a_{sat}} $$

$$ kex = f_X \sigma r_a \frac{a}{a + a_{sat}} $$

$$ keu = f_U \sigma r_a \frac{a}{a + a_{sat}} $$

$$ ker_i = k_{on}^{cm} r_a - k_{off}^{cm} r_i $$

The first equation gives another expression for the growth rate: $\alpha = ke = \sigma \frac{r_a}{1-a} \frac{a}{a+a_{sat}}$, where the factor $\frac{r_a}{1-a}$ is in fact the proteome fraction of active ribosomes $\phi_{R_a}$. It also reflects the balance between the rate of nutrient import and transformation into precursors and the rate of total protein synthesis in the context of dilution imposed by mass density homeostasis.

Combining the first equation with every other enables expressing protein concentrations as a function of $(1-a)$:

$$ e = f_E(1 - a), r = f_R(1 - a), q = f_Q(1 - a) \text{ and } x = f_X(1 - a) $$
This is intuitive: a proteome sector concentration is equal to its allocation fraction times the total protein concentration $1-a$. We also see that steady-state proteome fractions $\varphi$ are equal to allocation fractions $f$.

Noting that $\frac{r_i}{r_a} = \frac{k_{on}^{cm}}{k_{off}^{cm} + k_{E}(1-a)}$ (from the steady-state equation for $r_i$ together with $e = f_E(1-a)$) and that $f_R \sigma \frac{a}{a+a_{sat}} = k_{on}^{cm} + k_{E}(1-a) - k_{off}^{cm} \frac{r_i}{r_a}$ (from the steady-state equation for $r_a$ divided by $r_a$) leads to an equation for $a$ by substituting the first expression in the second:

$$f_R \sigma \frac{a}{a + a_{sat}} = k_{E}(1 - a) \frac{k_{E}(1 - a) + k_{on}^{cm} + k_{off}^{cm}}{k_{E}(1 - a) + k_{off}^{cm}}$$

Solving this equation for $a$ allows to compute the steady-state solution for fixed allocation fractions $f_E, f_R, f_Q, f_X$ and $f_U$.

## Dynamic regulation of allocation fractions

Dai and colleagues [49] found that the translation elongation rate displays a Michaelis-Menten like dependence with the RNA to protein ratio (which is proportional to the total ribosome proteome fraction) across nutrient conditions and chloramphenicol-mediated translation inhibition. A simple way to reproduce these findings with our model is to assume that the allocation fraction $f_R$ is regulated by the concentration of precursors via direct proportionality: $f_R = \delta \cdot a$. Indeed, we would then have $\frac{a}{a+a_{sat}} = \frac{\delta^{-1}f_R}{\delta^{-1}f_R + a_{sat}} = \frac{f_R}{f_R + \delta a_{sat}}$. Note that such regulation is a form of supply-driven activation [10].

Therefore, our model for how allocations fractions are dynamically set is as follows:

* $f_Q, f_U, f_X$ are constants ($f_Q$ is also invariant across conditions)

* $f_R = \delta \cdot a$ (or $1-f_Q-f_U-f_X$ if $\delta \cdot a > 1-f_Q-f_U-f_X$)

* $f_E = 1-f_R-f_Q-f_U-f_X$

The ordinary differential equations (ODEs) presented above with such dynamically changing allocations fractions are used to simulate Fig 2A, together with the halving of all amounts at division, which is triggered when $X = X_{div}$. The integration of the ODEs together with the handling of division events is natively supported by MATLAB ODEs solving suite. We used ode45 with default tolerances. A pseudo-code description of the simulation is therefore:

* Initialize amounts of A, E, RA, RI, Q, U, and X

* Integrate the ODEs (using parameters $\delta, a_{sat} = K/\delta, f_Q, f_X, k_{off}^{cm}, k, k_{on}^{cm}, f_U$) until X reaches the division threshold $X_{div}$

* Divide all amounts by 2 and repeat previous step

## Steady-state proteome fractions for the regulation model

By substituting $f_R = \delta \cdot a$ in the steady-state growth equations above, we obtain the following equation if we further assume that $a \ll 1$ (i.e. $\delta$ is large enough):

$$f_R \sigma \frac{f_R}{f_R + K} = k f_E \left( 1 + \frac{k_{on}^{cm}}{k_{off}^{cm} + k f_E} \right), \text{ where } K = \delta a_{sat}$$

Because $f_E = 1-f_R-f_Q-f_U-f_X$, this can be solved numerically by scalar optimization on $f_R$.

This steady-state approximation was used for the results presented in Figs 1 and 3. However, when simulating the dynamic model (Figs 2 and 4), a value for $\delta$ has to be chosen, and
the allocation fraction $f_X$ has to be non-zero. We chose $\delta = 5$ and we subtracted $f_X$ from $f_Q$ to ensure $f_Q+f_X = 0.5$. We show in Fig A in S1 Text the validity of the $a \ll 1$ approximation.

## Optimality of proteome allocation fractions

Above we considered the case where the allocation fractions are dictated by a supply-driven activation of ribosome synthesis. For given growth conditions, one could also search for the proteome allocation fractions that maximize growth rate (under the constraints of a fixed, condition-independent allocation fraction to housekeeping proteins $f_Q$, an imposed allocation fraction to useless proteins $f_U$, and under the assumption that $f_X \ll 1$).

This is achieved by enumerating all possible allocations by varying $f_R$ in between 0 and 1 $-f_Q-f_U$, and computing the growth rate of the corresponding steady-state (we use the bounded scalar optimization MATLAB function fminbnd).

Such optimality assumptions are sometimes made in previous work, such as in Molenaar et al. (2009).

## Model parameterization

Cell composition data is taken from (Scott et al. 2010) and (Dai et al. 2016) [13,49]. It consists of steady-state ribosome proteome fractions for a combination of nutrient conditions and chloramphenicol concentrations (measured via the RNA/protein ratio). The dataset from Dai et al. (2016) also provides translation elongation rate and estimated fraction of ribosomes that are active. Following Scott and colleagues [13], we used a conversion factor of 0.76 to convert RNA/protein ratios into proteome fractions of extended ribosomes. For consistency, we also used the notion of extended ribosomes (1.67 times larger than a single ribosome, Scott et al., 2010) to convert the maximal elongation rate of 22 aa/s measured by Dai and colleagues [49] into a maximal rate of protein synthesis ($\sigma$ in our model). Also, based on maximal ribosome proteome fractions observed by Scott and colleagues, we set the value of $f_Q$ to 0.5. Finally, we used the value of the unbinding constant $k_{off}^{cm}$ for chloramphenicol-ribosome complexes given in [49]. Altogether, we obtained the following parameterization from those studies:

$$ \sigma = \frac{1}{\frac{7336 \ aa*1.67}{22 \ aa.s^{-1}}} s^{-1} = 6.46 \ hr^{-1} $$

$$ K = 0.11 * 0.76 = 0.0836 $$

$$ f_Q = 0.5 $$

$$ k_{off}^{cm} = 5.04 \ hr^{-1} $$

## Finding regulation of $f_X$ expression explaining cell size data across types of growth rate modulation

The structural model of cell division links the steady-state cell composition and cell division size via the concentration $x$ of X proteins: $V_{div} = \frac{x_{div}}{x} = X_{div} \frac{1}{f_X(1-a)}$, which approximates to $\frac{x_{div}}{f_X}$ under the low precursor concentration assumption.

To search for regulation of $f_X$ explaining cell size data across the three types of growth rate modulations, we considered the following form of $f_X$ : $f_X \propto c_1^{\beta_1} c_2^{\beta_2}$ where $c_{1,2}$ are quantities depending on the coarse-grained cell composition and $\beta_{1,2}$ are exponents. This functional form was chosen for two reasons. First it the simplest form beyond regulation only by a single
factor. Second, regulation of gene expression by several types of molecules (such as transcription factors and co-factors) is likely to be multiplicative.

We considered the following coarse-grained quantities: coarse-grained concentrations *e*, *r*, *r*<sub>a</sub> and the fraction of active ribosomes *r*<sub>a</sub>/r (not to confuse with the mass fraction of active ribosomes *r*<sub>a</sub> or the proteome fraction of active ribosomes *f*<sub>RA</sub> = (*r*<sub>a</sub>−*r*<sub>f</sub>)/(*1*−*a*)).

Fitting was performed using multilinear regression on the logarithm of cell size data using MATLAB function *regress*. 95% confidence intervals on exponents were provided by the *regress* function. Cell size data between different datasets was normalized as detailed in Fig C in *S1 Text*.

To estimate the C+D duration directly from our size estimate, we use the value of ‘unit cell’ or DNA replication initiation volume from (Si et al. 2017) (S₀ = 0.28 µm³, converted to the normalized size scale described in Fig C in S1 Text).

**Stochastic model**

The stochastic version of the model was simulated using the *Gillespie* algorithm and in a
mother machine setting. When *X* reaches the division threshold *X*<sub>div</sub>, each coarse-grained mol-
ecule is kept in the daughter cell that we keep simulating with probability 1/2. Note that we are
modeling protein synthesis with a single reaction for each sector, hence abstracting away tran-
scription and mRNA degradation. Therefore, there are two sources of noise in the model: sto-
chasticity in reactions and partitioning of molecules at cell division.

Code availability

All the code for solving model steady-state, simulating model dynamics, parameter fitting, etc. as well as the scripts generating all figures is available on GitHub. It is written in MATLAB except the stochastic simulation code written in C++. All simulations and calculations were performed on a DELL laptop (Precision 5520, Windows 10). MATLAB version 2016b was used and only built-in MATLAB functions were used.

**Supporting information**

**S1 Text.** Supporting information for “a bacterial size law revealed by a coarse-grained model of cell physiology”.

(PDF)

Acknowledgments

We would like to thank Suckjoon Jun, Istvan Kleijn, Guillaume Le Treut, Peter Sarkies, Fang-wei Si, Peter Swain and Philipp Thomas for feedback on our manuscript and members of the Marguerat and Shahrezaei groups for discussions.

**Author Contributions**

**Conceptualization:** François Bertaux, Samuel Marguerat, Vahid Shahrezaei.

**Data curation:** François Bertaux.

**Formal analysis:** François Bertaux, Vahid Shahrezaei.

**Funding acquisition:** Samuel Marguerat, Vahid Shahrezaei.

**Investigation:** François Bertaux, Julius von Kügelgen, Samuel Marguerat, Vahid Shahrezaei.
**Methodology:** François Bertaux, Julius von Kügelgen, Vahid Shahrezaei.

**Project administration:** Samuel Marguerat, Vahid Shahrezaei.

**Software:** François Bertaux.

**Supervision:** Samuel Marguerat, Vahid Shahrezaei.

**Writing – original draft:** François Bertaux, Samuel Marguerat, Vahid Shahrezaei.

**Writing – review & editing:** François Bertaux, Samuel Marguerat, Vahid Shahrezaei.

## References

1. Schaechter M, MaalOe O, Kjeldgaard NO. Dependency on Medium and Temperature of Cell Size and Chemical Composition during Balanced Growth of Salmonella typhimurium. J Gen Microbiol. 1958; 19: 592–606. https://doi.org/10.1099/00221287-19-3-592 PMID: 13611202

2. Schleif R. Control of production of ribosomal protein. J Mol Biol. 1967; 27: 41–55. https://doi.org/10.1016/0022-2836(67)90350-6 PMID: 5340469

3. Cooper S, Helmstetter CE. Chromosome replication and the division cycle of Escherichia coli. J Mol Biol. 1968; 31: 519–540. https://doi.org/10.1016/0022-2836(68)90425-7 PMID: 4866337

4. Maaløe O. An analysis of bacterial growth. Dev Biol. 1969; 3 (Suppl): 33–58.

5. Churchward G, Bremer H, Young R. Macromolecular composition of bacteria. J Theor Biol. 1982; 94: 651–670. https://doi.org/10.1016/0022-5193(82)90305-8 PMID: 6176814

6. Klumpp S, Hwa T. Bacterial growth: global effects on gene expression, growth feedback and proteome partition. Curr Opin Biotechnol. 2014; 28: 96–102. https://doi.org/10.1016/j.copbio.2014.01.001 PMID: 24495512

7. Vadia S, Levin PA. Growth rate and cell size: a re-examination of the growth law. Curr Opin Microbiol. 2015; 24: 96–103. https://doi.org/10.1016/j.mib.2015.01.011 PMID: 25662920

8. Jun S, Si F, Pugatch R, Scott M. Fundamental principles in bacterial physiology—history, recent progress, and the future with focus on cell size control: a review. Rep Prog Phys. 2018; 81: 056601. https://doi.org/10.1088/1361-6633/aaa628 PMID: 29313526

9. Hill NS, Buske PJ, Shi Y, Levin PA. A Moonlighting Enzyme Links Escherichia coli Cell Size with Central Metabolism. Casadesús J, editor. PLoS Genet. 2013; 9: e1003663. https://doi.org/10.1371/journal.pgen.1003663 PMID: 23935518

10. Scott M, Klumpp S, Mateescu EM, Hwa T. Emergence of robust growth laws from optimal regulation of ribosome synthesis. Mol Syst Biol. 2014; 10: 747–747. https://doi.org/10.15252/msb.20145379 PMID: 25149558

11. Neidhardt FC, Magasanik B. Studies on the role of ribonucleic acid in the growth of bacteria. Biochim Biophys Acta. 1960; 42: 99–116. https://doi.org/10.1016/0006-3002(60)90757-5 PMID: 13728193

12. Basan M. Resource allocation and metabolism: the search for governing principles. Curr Opin Microbiol. 2018; 45: 77–83. https://doi.org/10.1016/j.mib.2018.02.008 PMID: 29544124

13. Scott M, Gunderson CW, Mateescu EM, Zhang Z, Hwa T. Interdependence of Cell Growth and Gene Expression: Origins and Consequences. Science. 2010; 330: 1099–1102. https://doi.org/10.1126/science.1192588 PMID: 21097934

14. Marr AG. Growth rate of Escherichia coli. Microbiol Rev. 1991; 55: 316–333. PMID: 1886524

15. Molenaar D, Berlo R van, Ridder D de, Teusink B. Shifts in growth strategies reflect tradeoffs in cellular economics. Mol Syst Biol. 2009; 5. https://doi.org/10.1038/msb.2009.82 PMID: 19888218

16. Weiße AY, Oyarzún DA, Danos V, Swain PS. Mechanistic links between cellular trade-offs, gene expression, and growth. Proc Natl Acad Sci. 2015; 112: E1038–E1047. https://doi.org/10.1073/pnas.1416533112 PMID: 25695966

17. Maitra A, Dill KA. Bacterial growth laws reflect the evolutionary importance of energy efficiency. Proc Natl Acad Sci U S A. 2015; 112: 406–411. https://doi.org/10.1073/pnas.1421138111 PMID: 25548180

18. Bosdriesz E, Molenaar D, Teusink B, Bruggeman FJ. How fast-growing bacteria robustly tune their ribosome concentration to approximate growth-rate maximization. FEBS J. 2015; 282: 2029–2044. https://doi.org/10.1111/febs.13258 PMID: 25754869

19. Maitra A, Dill KA. Modeling the Overproduction of Ribosomes when Antibacterial Drugs Act on Cells. Biophys J. 2016; 110: 743–748. https://doi.org/10.1016/j.bpj.2015.12.016 PMID: 26840738
20. Pandey PP, Jain S. Analytic derivation of bacterial growth laws from a simple model of intracellular chemical dynamics. Theory Biosci. 2016; 135: 121–130. https://doi.org/10.1007/s12064-016-0227-9 PMID: 27167220

21. Giordano N, Mairet F, Gouzé J-L, Geiselmann J, de Jong H. Dynamical Allocation of Cellular Resources as an Optimal Control Problem: Novel Insights into Microbial Growth Strategies. Igoshin OA, editor. PLOS Comput Biol. 2016; 12: e1004802. https://doi.org/10.1371/journal.pcbi.1004802 PMID: 26958858

22. Liao C, Blanchard AE, Lu T. An integrative circuit–host modelling framework for predicting synthetic gene network behaviours. Nat Microbiol. 2017; 2: 1658–1666. https://doi.org/10.1038/s41564-017-0022-5 PMID: 28947816

23. Thomas P, Terradot G, Danos V, Weiße AY. Sources, propagation and consequences of stochasticity in cellular growth. Nat Commun. 2018; 9: 4528. https://doi.org/10.1038/s41467-018-06912-9 PMID: 30375377

24. Sharma P, Pandey PP, Jain S. Modeling the cost and benefit of proteome regulation in a growing bacterial cell. Phys Biol. 2018; 15: 046005. https://doi.org/10.1088/1478-3975/aabe43 PMID: 29658492

25. Pandey PP, Singh H, Jain S. Exponential trajectories, cell size fluctuations and the adder property in bacteria follow from simple chemical dynamics and division control. Systems Biology; 2018 Dec. https://doi.org/10.1101/487504

26. Donachie WD. Relationship between Cell Size and Time of Initiation of DNA Replication. Nature. 1968; 219: 1077–1079. https://doi.org/10.1038/2191077a0 PMID: 4876941

27. Wallden M, Fange D, Lundius EG, Baltekin Ö, Elf J. The Synchronization of Replication and Division Cycles in Individual E. coli Cells. Cell. 2016; 166: 729–739. https://doi.org/10.1016/j.cell.2016.06.052 PMID: 27471967

28. Si F, Li D, Cox SE, Sauls JT, Azizi O, Sou C, et al. Invariance of Initiation Mass and Predictability of Cell Size in Escherichia coli. Curr Biol CB. 2017; 27: 1278–1287. https://doi.org/10.1016/j.cub.2017.03.022 PMID: 28416114

29. Basan M, Zhu M, Dai X, Warren M, Sévin D, Wang Y-P, et al. Inflating bacterial cells by increased protein synthesis. Mol Syst Biol. 2015;11. https://doi.org/10.15252/msb.20156178 PMID: 26519362

30. Taheri-Araghi S, Bradde S, Sauls JT, Hill NS, Levin PA, Paulsson J, et al. Cell-Size Control and Homeostasis in Bacteria. Curr Biol. 2015; 25: 385–391. https://doi.org/10.1016/j.cub.2014.12.009 PMID: 25544609

31. Kennard AS, Osella M, Javer A, Grilli J, Nghe P, Tans SJ, et al. Individuality and universality in the growth-division laws of single E. coli cells. Phys Rev E. 2016; 93: 012408. https://doi.org/10.1103/PhysRevE.93.012408 PMID: 26871102

32. Campos M, Surovtsev IV, Kato S, Paintdakhi A, Beltran B, Ebmeier SE, et al. A Constant Size Extension Drives Bacterial Cell Size Homeostasis. Cell. 2014; 159: 1433–1446. https://doi.org/10.1016/j.cell.2014.11.022 PMID: 25480302

33. Amir A. Cell Size Regulation in Bacteria. Phys Rev Lett. 2014; 112: 208102. https://doi.org/10.1103/PhysRevLett.112.208102

34. Iyer-Biswas S, Wright CS, Henry JT, Lo K, Burov S, Lin Y, et al. Scaling laws governing stochastic growth and division of single bacterial cells. Proc Natl Acad Sci U S A. 2014; 111: 15912–15917. https://doi.org/10.1073/pnas.1403232111 PMID: 25349411

35. Deforet M, van Ditmarsch D, Xavier JB. Cell-Size Homeostasis and the Incremental Rule in a Bacterial Pathogen. Biophys J. 2015; 109: 521–528. https://doi.org/10.1016/j.bpj.2015.07.002 PMID: 26244734

36. Soifer I, Robert L, Amir A. Single-Cell Analysis of Growth in Budding Yeast and Bacteria Reveals a Common Size Regulation Strategy. Curr Biol. 2016; 26: 356–361. https://doi.org/10.1016/j.cub.2015.11.067 PMID: 26776734

37. Chandler-Brown D, Schmoller KM, Winetraub Y, Skotheim JM. The Adder Phenomenon Emerges from Independent Control of Pre- and Post-Start Phases of the Budding Yeast Cell Cycle. Curr Biol. 2017; 27: 2774–2783.e3. https://doi.org/10.1016/j.cub.2017.08.015 PMID: 28889980

38. Logsdon MM, Ho P-Y, Papavinasasundaram K, Richardson K, Cokol M, Sassetti CM, et al. A Parallel Adder Coordinates Mycobacterial Cell-Cycle Progression and Cell-Size Homeostasis in the Context of Asymmetric Growth and Organization. Curr Biol. 2017; 27: 3367–3374.e7. https://doi.org/10.1016/j.cub.2017.09.046 PMID: 29107550

39. Priestman M, Thomas P, Robertson BD, Shahrezaei V. Mycobacteria Modify Their Cell Size Control under Sub-Optimal Carbon Sources. Front Cell Dev Biol. 2017; 5: 64. https://doi.org/10.3389/fcell.2017.00064 PMID: 28748182
40. Varsano G, Wang Y, Wu M. Probing Mammalian Cell Size Homeostasis by Channel-Assisted Cell Reshaping. Cell Rep. 2017; 20: 397–410. https://doi.org/10.1016/j.celrep.2017.06.057 PMID: 28700941

41. Eun Y-J, Ho P-Y, Kim M, LaRussa S, Robert L, Renner LD, et al. Archaeal cells share common size control with bacteria despite noisier growth and division. Nat Microbiol. 2018; 3: 148–154. https://doi.org/10.1038/s41564-017-0082-6 PMID: 29255255

42. Cadart C, Monnier S, Grilli J, Sáez PJ, Srivastava N, Attia R, et al. Size control in mammalian cells involves modulation of both growth rate and cell cycle duration. Nat Commun. 2018; 9: 3275. https://doi.org/10.1038/s41467-018-05393-0 PMID: 30115907

43. Ho P-Y, Amir A. Simultaneous regulation of cell size and chromosome replication in bacteria. Front Microbiol. 2015; 6: 662. https://doi.org/10.3389/fmicb.2015.00662 PMID: 26217311

44. Amir A. Is cell size a spandrel? eLife. 2017; 6: e22186. https://doi.org/10.7554/eLife.22186 PMID: 28102818

45. Harris LK, Theriot JA. Relative Rates of Surface and Volume Synthesis Set Bacterial Cell Size. Cell. 2016; 165: 1479–1492. https://doi.org/10.1016/j.cell.2016.05.045 PMID: 27259152

46. Micali G, Grilli J, Osella M, Cosentino Lagomarsino M. Concurrent processes set E. coli cell division. Sci Adv. 2018; 4: eaau3324. https://doi.org/10.1126/sciadv.aau3324 PMID: 30417095

47. Si F, Le Treut G, Sauls JT, Vadia S, Levin PA, Jun S. Mechanistic Origin of Cell-Size Control and Homeostasis in Bacteria. Curr Biol. 2019; 29: 1760–1770.e7. https://doi.org/10.1016/j.cub.2019.04.062 PMID: 31104932

48. Fantes PA, Grant WD, Pritchard RH, Sudbery PE, Wheals AE. The regulation of cell size and the control of mitosis. J Theor Biol. 1975; 50: 213–244. https://doi.org/10.1016/0022-5193(75)90034-x PMID: 1127959

49. Dai X, Zhu M, Warren M, Balakrishnan R, Patsalo V, Okano H, et al. Reduction of translating ribosomes enables Escherichia coli to maintain elongation rates during slow growth. Nat Microbiol. 2017; 2: 16231. https://doi.org/10.1038/nmicrobiol.2016.231 PMID: 27941827

50. Ehrenberg M, Kurland CG. Costs of accuracy determined by a maximal growth rate constraint. Q Rev Biophys. 1984; 17: 45–82. https://doi.org/10.1017/s0033583500005254 PMID: 6484121

51. Klumpp S, Scott M, Pedersen S, Hwa T. Molecular crowding limits translation and cell growth. Proc Natl Acad Sci. 2013; 110: 16754–16759. https://doi.org/10.1073/pnas.1310377110 PMID: 24082144

52. Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, et al. Robust growth of Escherichia coli. Curr Biol. 2010; 20: 1099–1103. https://doi.org/10.1016/j.cub.2010.04.045 PMID: 20537537

53. Elowitz MB. Stochastic Gene Expression in a Single Cell. Science. 2002; 297: 1183–1186. https://doi.org/10.1126/science.1070919 PMID: 12183631

54. Levine E, Hwa T. Stochastic fluctuations in metabolic pathways. Proc Natl Acad Sci U S A. 2007; 104: 9224–9229. https://doi.org/10.1073/pnas.0610987104 PMID: 17517669

55. Dunlop MJ, Cox RS, Levine JH, Murray RM, Elowitz MB. Regulatory activity revealed by dynamic correlations in gene expression noise. Nat Genet. 2008; 40: 1493–1498. https://doi.org/10.1038/ng.281 PMID: 19029898

56. Kiviet DJ, Nghe P, Walker N, Boulineau S, Sunderlikova V, Tans SJ. Stochasticity of metabolism and growth at the single-cell level. Nature. 2014; 514: 376–379. https://doi.org/10.1038/nature13582 PMID: 25186725

57. Kleijn IT, Krah LHJ, Hermsen R. Noise propagation in an integrated model of bacterial gene expression and growth. Morozov AV, editor. PLOS Comput Biol. 2018; 14: e1006386. https://doi.org/10.1371/journal.pcbi.1006386 PMID: 30289879

58. Ghusinga KR, Vargas-Garcia CA, Singh A. A mechanistic stochastic framework for regulating bacterial cell division. Sci Rep. 2016; 6: 30229. https://doi.org/10.1038/srep30229 PMID: 27456660

59. Huls PG, Vischer NOE, Woldringh CL. Different Amounts of DNA in Newborn Cells of Escherichia coli Preclude a Role for the Chromosome in Size Control According to the “Adder” Model. Front Microbiol. 2018; 9: 664. https://doi.org/10.3389/fmicb.2018.00664 PMID: 29675011

60. Zhu M, Mori M, Hwa T, Dai X. Disruption of transcription–translation coordination in Escherichia coli leads to premature transcriptional termination. Nat Microbiol. 2019 [cited 11 Oct 2019]. https://doi.org/10.1038/s41564-019-0543-1 PMID: 31451774

61. Ojkic N, Serbanescu D, Banerjee S. Surface-to-volume scaling and aspect ratio preservation in rod-shaped bacteria. eLife. 2019; 8: e47033. https://doi.org/10.7554/eLife.47033 PMID: 31456563

62. Palacios P, Vicente M, Sánchez M. Dependency of Escherichia coli cell-division size, and indepen-dency of nucleoid segregation on the mode and level of ftsZ expression. Mol Microbiol. 1996; 20: 1093–1098. https://doi.org/10.1111/j.1365-2958.1996.tb02549.x PMID: 8809761
63. Zheng H, Ho P-Y, Jiang M, Tang B, Liu W, Li D, et al. Interrogating the Escherichia coli cell cycle by cell dimension perturbations. Proc Natl Acad Sci. 2016; 113: 15000–15005. https://doi.org/10.1073/pnas.1617932114 PMID: 27956612

64. Gray WT, Govers SK, Xiang Y, Parry BR, Campos M, Kim S, et al. Nucleoid Size Scaling and Intracellular Organization of Translation across Bacteria. Cell. 2019; 177: 1632–1648.e20. https://doi.org/10.1016/j.cell.2019.05.017 PMID: 31150626

65. Chien A-C, Hill NS, Levin PA. Cell Size Control in Bacteria. Curr Biol. 2012; 22: R340–R349. https://doi.org/10.1016/j.cub.2012.02.032 PMID: 22575476

66. Turner JJ, Ewald JC, Skotheim JM. Cell Size Control in Yeast. Curr Biol. 2012; 22: R350–R359. https://doi.org/10.1016/j.cub.2012.02.041 PMID: 22575477

67. Bertaux F, Marguerat S, Shahrezaei V. Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits. R Soc Open Sci. 2018; 5: 172234. https://doi.org/10.1098/rsos.172234 PMID: 29657814

68. Monds RD, Lee TK, Colavin A, Ursell T, Quan S, Cooper TF, et al. Systematic Perturbation of Cytoskeletal Function Reveals a Linear Scaling Relationship between Cell Geometry and Fitness. Cell Rep. 2014; 9: 1528–1537. https://doi.org/10.1016/j.celrep.2014.10.040 PMID: 25456141

69. Borkowski O, Goelzer A, Schaffer M, Calabre M, Mäder U, Aymerich S, et al. Translation elicits a growth rate-dependent, genome-wide, differential protein production in Bacillus subtilis. Mol Syst Biol. 2016; 12: 870. https://doi.org/10.15252/msb.20156608 PMID: 27193784