Analytic derivation of bacterial growth laws from a simple model of intracellular chemical dynamics.

Experiments have found that the growth rate and certain other macroscopic properties of bacterial cells in steady-state cultures depend upon the medium in a surprisingly simple manner; these dependencies are referred to as 'growth laws'. Here we construct a dynamical model of interacting intracellular populations to understand some of the growth laws. The model has only three population variables: an amino acid pool, a pool of enzymes that transport an external nutrient and produce the amino acids, and ribosomes that catalyze their own and the enzymes' production from the amino acids. We assume that the cell allocates its resources between the enzyme sector and the ribosomal sector to maximize its growth rate. We show that the empirical growth laws follow from this assumption and derive analytic expressions for the phenomenological parameters in terms of the more basic model parameters. Interestingly, the maximization of the growth rate of the cell as a whole implies that the cell allocates resources to the enzyme and ribosomal sectors in inverse proportion to their respective 'efficiencies'. The work introduces a mathematical scheme in which the cellular growth rate can be explicitly determined and shows that two large parameters, the number of amino acid residues per enzyme and per ribosome, are useful for making approximations.

Introduction

Bacterial cells contain thousands of molecular species and are exceedingly complex, yet they exhibit certain remarkable regularities at the system level which have been quantified experimentally. The regularities of concern in this paper are a subset of the so-called ‘bacterial growth laws’ (Monod 1949; Schaechter et al. 1958; Maaloe and Kjeldgaard 1966; Maaløe 1979; Bremer and Dennis 1996; Scott et al. 2010) which highlight the relationships between macroscopically measured quantities such as cell composition, size, growth rate and the environment or medium in which the cell grows. The empirical relationships are summarized in terms of phenomenological equations. In this paper we attempt to deduce these phenomenological relationships from a mathematical model of a cell containing a few interacting (pools of) molecular species. The population dynamics of these molecular species based on standard chemical kinetics, together with an optimization principle, gives rise to the growth laws.

When genetically identical bacterial cells drawn from an overnight culture are introduced into a vessel containing a medium with a certain concentration of nutrients, temperature, etc., they exhibit several phases of growth (Monod 1949). These include, in sequence, a lag phase where there is very little growth in the number of cells, an acceleration phase where growth picks up, an exponential phase in which the population of cells grows exponentially with time (at a constant growth rate), a deceleration phase with declining growth rate that sets in when the food begins to run out and a stationary phase where the population is constant, followed by an eventual population decline. Regularities are most apparent in the exponential phase which is often referred to as a steady state. In this phase the averages and distributions (across the population of cells) of cell doubling time, cell size at birth, intracellular concentration of ribosome, total protein and metabolites become constant in time (for as long as the exponential phase lasts). These constant average values depend upon the strain of bacteria and on the medium (its concentration of nutrients, temperature, etc.). Repeated experiments with the same strain and medium but with different initial conditions (corresponding to different overnight cultures) yield the same growth rate in the steady state and the same values of these averages. The growth laws are statements of how the growth rate and these averages depend upon the environment and cellular parameters. The first of these, due to Monod (1949), is the hyperbolic dependence of the steady-state growth rate upon the concentration [F] of a growth-limiting nutrient (or food molecule) in the medium: is the maximum value of the growth rate possible in the medium and the value of [F] at which the growth rate is half its maximum value.

In the cell, the ribosome which assembles amino acids to produce proteins from a messenger RNA template is an important catalyst of cell growth. The amount of cellular investment in ribosomes is found to depend upon the growth rate in a characteristic manner. In particular, the ratio of ribosomal protein in the cell to total protein in the cell (by weight), referred to as the ‘ribosomal fraction’ , is found to be a linear increasing function of when is increased by improving the nutritional quality of the medium (Schaechter et al. 1958; Maaløe 1979; Bremer and Dennis 1996):where and are constants. However, when is altered by changing the catalytic efficiency of ribosomes (e.g., by producing mutants with different catalytic efficiencies or by adding antibiotics in the medium that particularly affect the catalytic efficiency) keeping the nutritional quality of the medium the same, then is found to be a linear decreasing function of (Scott et al. 2010):where and are constants. The above three equations can be considered to be phenomenological equations describing bacterial growth steady states, with the six constants as phenomenological constants (Scott et al. 2010). The simplicity and universality of these phenomenological laws are surprising given the complexity and diversity of bacteria. In addition to the above growth laws, the size of bacterial cells also exhibits remarkable properties which are not the subject of this paper.

There have been several recent works which have attempted to understand the growth laws theoretically, through mathematical modeling (Molenaar et al. 2009; Scott et al. 2010, 2014; Maitra and Dill 2015; Weiße et al. 2015; Bosdriesz et al. 2015). Scott et al. (2010, 2014) have related the phenomenological constants to molecular parameters of the cell. Taking forward an idea due to Maaløe (1979), they have argued that the growth laws reflect regulatory mechanisms in the cell that optimize its growth rate in any given medium. They and other authors (Maitra and Dill 2015; Weiße et al. 2015; Bosdriesz et al. 2015) have constructed models for the molecular regulatory mechanisms inside the cell that can produce the above growth laws.

In this paper we adopt a different approach that is closer in spirit to the work of Molenaar et al. (2009). Molenaar et al. considered a nonlinear dynamical model of a cell with a few classes of metabolites and enzymes as well as ribosomes and showed through computer simulations that maximization of the cellular growth rate qualitatively reproduced some of the growth laws and other observed properties of cells. Here we consider a simpler nonlinear dynamical model of the cell containing only three molecular populations: one metabolite pool, one enzyme pool and ribosomes. We are able to obtain an explicit formula for the growth rate of the cell as a function of cellular and medium parameters, which has so far been lacking in existing models. Maximizing the growth rate with respect to one of the parameters, the fraction of ribosomes making ribosomes, we derive all the three growth laws analytically. The method produces analytic expressions for the phenomenological parameters in terms of the molecular parameters in the model. These expressions are generalizations of the ones obtained by Scott et al. and reduce to their results when certain processes are ignored. We show that the optimization of growth rate leads to a simple principle of cellular economy. The work provides a direct connection between growth rate optimization and the growth laws.

At a methodological level we identify natural large parameters in the cell that are useful in making approximations. This might prove useful in more complex cellular models and in modeling other cellular phenomena as well.

Precursor-Transporter-Ribosome (PTR) cell: a coarse grained model

Consider a simple mathematical model of a growing cell consisting of three types of molecules; precursors, transporters and ribosomes. We refer to this model as the Precursor-Transporter-Ribosome (PTR) model. The system has the following three reactions (Fig. 1):

Fig. 1

The PTR cell. Precursor molecules (P) are produced by the catalytic action of the metabolic proteins (T) on the external food molecules (F). Metabolic proteins and ribosomal proteins (R) are synthesized from the P molecules in reactions catalysed by R

, where external food molecules (F) are transported into the cell by the action of transporter proteins (T) and converted into precursor molecules (P) representing amino acids;

, where P molecules are converted into T by the catalytic action of ribosomes (R), and

, where R catalyses the production of itself using P.

All the molecules are produced in the interior of the cell. The membrane consists solely of transporter molecules, which are assumed to migrate immediately upon formation to the cell boundary. The interior of the cell consists of precursor molecules and ribosomes. The model is described by the following set of differential equations: where P represents the number of precursor molecules in the cell (amino acid pool), T the number of all metabolic protein molecules that transport food into the cell and convert it into precursor and R is number of ribosomes in the cell. The rate constant represents the efficiency of metabolism in making P from external food. It is an increasing function of the external food concentration [F] (explicit forms to be discussed later) and can also encapsulate the quality of the food source (e.g., the number of P molecules produced per food molecule transported in). The other production rate constants are parametrized as follows:where k represents ribosomal catalytic efficiency and is the rate at which a single ribosome consumes P molecules, per unit concentration of P, for the production of proteins. This accounts for the term kRP/V in the equation, the total rate of consumption of P. A fraction of the ribosomes makes the T proteins and the remaining fraction the ribosomal proteins. Thus, of the P consumption flux a part goes to produce T and the remaining part goes to produce R. Each T molecule (ribosome) contains () amino acid residues; hence the rate of production of T is and that of R is . This explains the assumed forms of and . and are the degradation rates of T and R, respectively, into a waste product; we assume a negligible degradation rate for P.

V is the instantaneous volume of the interior of the cell and we assume that it is a linear function of the molecular populations. Since molecular populations in the bulk are P and R, we can take it to be proportional to . Our results do not depend upon this particular choice and for generality we assumewhere are constants. Note that Eqs. (4a)–(4c) do not contain a term proportional to on the right-hand side because they refer to populations instead of concentrations.

Steady-state solution of the PTR cell

The steady state of a bacterial culture corresponds to cells growing exponentially with a constant rate. We look for an exponential solution for the chemical populations: , where , a constant, is the growth rate of the PTR cell. Substituting this ansatz into Eq. (4), we get where . Henceforth we drop the subscript 0 as the equations are valid for the time-dependent quantities P(t), T(t), R(t) as well. The last of these equations immediately givesSubstituting (8) in (7b) gives the ratio T/R:and substituting (8) and (9) in (7a) gives the ratio P/R:Thus the ratios of the populations and the concentrations of the three chemicals at steady state can be expressed in terms of and the parameters of the model. In order to solve the problem fully, we need to find in terms of the parameters.

Growth rate   The Eq. (8) gives . Note that V can be written as . Thus P/V is completely expressed in terms of the ratios T/R and P/R which are known as functions of and the parameters [Eqs. (9) and (10)]. Therefore, the equation becomes an equation that contains only and the parameters. Simplifying it, we get a quadratic equation in with coefficients depending on the parameters:with

Equation (11) has two solutions:The solution is the physically relevant one, in which the square-root is always taken with the negative sign. There are several ways to see this:We remark here that it has been possible to obtain an explicit solution for because we have expressed the cell volume as a function of the populations and further assumed that it is a linear function of the populations, (6). This assumption (a) makes the exponential ansatz a solution of (4), and (b) causes the absolute populations to be eliminated from (8), leaving an equation connecting and the parameters. In our view the volume assumption is a crucial one that has been missing from previous models.

The Eqs. (4a)–(4c) can be simulated numerically for a fixed set of parameter values and initial conditions. This was done for several parameter sets and initial conditions. We found that at large times P, T and R always grew exponentially with time and their rate of exponential growth was given by and not . Further, the observed asymptotic ratios were given by Eqs. (9), (10) with . (Parameter values had to be chosen such that . When parameter values were such that , an exponential decline of populations was observed instead of growth.)

One can examine the two limits and . When , in both these limits must go to zero. Physically, when , then and Eq. (4c) implies that ribosomes are not produced; hence R is a constant, or . When , then , and T is not produced; hence again . It is easy to see that goes to zero in both these limits and not .

We have verified analytically from Eq. (7) that when and , gives rise to negative populations while gives rise to positive populations.

Ribosomal fraction ()  The ratio of ribosomal protein to total protein (by weight) is given by . Using Eq. (9) becomesNotice that this expression for is a nonlinear function of if and a constant independent of if . This is quite different from the observed linear growth laws (2) and (3). Thus the PTR model does not reproduce the observed growth laws. The model as it stands is missing an important ingredient—regulation—that we now turn to.

The PTR model with ‘regulation’ and bacterial growth laws

Upto now we have treated and , the fraction of ribosomes catalysing the production of transporters and ribosomal protein, respectively, as constant parameters of the model. However, it is a well-known fact that regulatory mechanisms exist in bacteria that regulate how much ribosome is engaged in producing ribosomal protein and how much in producing metabolic protein. In the context of the PTR model these mechanisms would modulate the value of the parameter (and hence ). The absence of this mechanism in the PTR model as described above is the reason that it does not reproduce the observed growth laws.

Trade-off between metabolic and ribosomal protein production   Since is a function of the cellular and medium parameters [Eq. (13)], we first ask how it varies as is increased keeping the medium and all other cellular parameters fixed. Numerical analysis of the steady-state of the PTR model shows that when all other parameters are fixed, is a non-monotonic function of as shown in Fig. 2a. This reflects a trade-off between production of metabolic proteins and ribosomal proteins in the model. There is a distinct value of () where is a maximum (). , depend upon the other parameters and in particular, increases as is increased (keeping the others constant). For convenience we here write where q equals the number of P molecules produced per food molecule consumed (quality of the medium), and depends upon external food concentration. We observe in Fig. 2a that as the quality of medium is increased, increases. These two properties, namely non-monotonicity of with respect to and the increase of with medium quality have also been noted in Scott et al. (2014) using a different approach.

Fig. 2

The PTR model in the optimized steady state qualitatively reproduces the observed growth laws. a Trade-off between production of ribosomal and metabolic proteins: as a function of for different values of q (, , , , , , , ). b as a function of q. Other parameters same as in (a). c versus for different values of q and k (other parameters same as before). The coloured lines (positive slope) correspond to Eq. (2) (changing medium quality at fixed translational efficiency). The grey lines (negative slope) correspond to (3) (changing translational efficiency at fixed medium quality)

Incorporating the effect of regulatory mechanisms through an optimization assumption  In order to bring in regulatory mechanisms we can make the rate constants dependent on molecular concentrations reflecting feedback mechanisms or introduce other molecular species (the regulators) into the model (Scott et al. 2014; Maitra and Dill 2015; Weiße et al. 2015; Bosdriesz et al. 2015). However, in the interest of mathematical simplicity we take an alternative approach involving optimization, employed earlier by Molenaar et al. (2009) for a different model. We assume that for any fixed medium and other cellular parameters additional regulatory mechanisms existing in the cell act to modify (e.g., by changing the proportion of messenger RNA molecules corresponding to R and T) such that the cellular growth rate is optimized, i.e., for a given medium and other cellular parameters, the regulation adjusts to . This is in spirit similar to the optimality assumption made in flux balance analysis of metabolic networks (Orth et al. 2010). In other words, we assume that the steady state reached when these other (unspecified) regulatory dynamics are included is approximated by the steady state of the PTR model withwhere is the value of that maximizes [Eq. (13)] keeping all other parameters fixed. We call this steady state (when is set equal to ) as the optimized steady-state. A change in medium, in general, leads to a different since is a function of all the medium-dependent parameters parameters and .

Optimized steady state of the PTR cell reproduces qualitative features of observed growth laws  Figure 2b, c shows that the optimized steady state of the PTR cell qualitatively satisfies the growth laws summarized in Eqs. (1)–(3). In Fig. 2b to increase the growth rate for the PTR cell we only increase the medium quality q (keeping constant). For each medium quality q we numerically obtain using Eq. (13), i.e., the value of that gives the largest for the given q. We denote this optimized as since it depends on q. We plot the dependence of on q and find the qualitative behaviour similar to the Monod curve (1).

Next we show the dependence of ribosome fraction on in the optimized steady state when is increased by increasing q. For each medium quality q, we already have and . To obtain we use the relation Eq. (14) with and . Figure 2c shows the plot of versus as the quality of the medium is increased. The lines with positive slope in Fig. 2c correspond to this variation. Notice the linear behaviour of the curves as has been observed in experiments, Eq. (2).

For a smaller value of k (smaller ribosomal efficiency), the versus curve remains linear but with a larger slope (coloured lines in 2c) as has been observed in experiments (Scott et al. 2010). Figure 2c qualitatively reproduces the observed behaviour of [Eqs. (2), (3)] when the growth rate is varied, both by increasing medium quality and by decreasing ribosomal efficiency.

Analytic derivation of the growth laws for the PTR model—the large  approximation The above results obtained numerically and illustrated in Fig. 2 can be derived analytically. The expressions turn out to be very simple when and are much greater than unity, which we assume in the following ( being the number of amino acid molecules needed to make an enzyme is , and , the number of amino acids in all ribosomal protein per ribosome is Bremer and Dennis 1996). We also need to assume that the parameters and defined in (12) are independent of and ; in other words, , the rate of P production per unit T molecule in the cell, and , the rate at which a ribosome adds amino acids to a protein, are sufficiently large, in fact, respectively, of order . and will turn out to be the two natural time scales that determine the system level properties of the cell. The time scales , and the volume will also be assumed to be independent of . may be independent or only weakly dependent on , respectively (sublinear dependence). With these assumptions, and .

Then, as shown in the “Appendix”,and the optimized steady-state growth rate of the PTR cell is given byThis leads to the Monod curve as will be discussed later.

Substituting and in Eqs. (9) and (14) givesThis expresses the ribosomal fraction at the optimized steady state of the PTR cell in terms of medium and cellular parameters. The growth laws in the standard form (2), (3) follow from (19) and (17). For example, to understand the dependence of on when the medium quality is varied, one can eliminate in favour of in Eq. (17) and substitute that in Eq. (19). This yields (2) withSimilarly, one can eliminate in favour of from Eq.  (17) and substitute in Eq. (19) to get Eq. (3), withThis reproduces the equations of the growth laws and expresses the constants appearing in those equations in terms of the medium and cellular parameters. Equation (2) with parameters given by (20) describes the positive slope lines in Fig. 2c and Eq. (3) with parameters given in (21) describes the negative slope lines.

Discussion

Nutritional and ribosomal efficiency  We now discuss the meaning of the formulae obtained. The formulae are expressed in terms of two quantities and and it is useful to interpret these quantities first. We follow Scott et al. (2014) in calling the ‘nutritional efficiency’ of the PTR cell in the given medium. Note that the production term in is . Since is the total number P molecules locked up in T, is the number of amino acid molecules produced in the cell per unit time per amino acid residue locked up in the metabolic enzymes. , being the rate of P production per unit P invested in metabolic enzymes, is appropriately the ‘metabolic efficiency’ or ‘nutritional efficiency’ of the cell in the given environment. In order to see the meaning of it is convenient to consider the situation where the concentration of P is high enough so that its availability is no longer a limiting factor for ribosomal activity. In the model the largest value of P/V is , which arises when the contribution of P to the volume dominates over the contribution from T and R, i.e., . Then (4c) becomes . Then R by itself forms an autocatalytic set (ACS) with growth rate . is the maximal growth rate of this ACS (when and ), or the rate at which R can make copies of itself if it was solely focused on doing that (that is, if ). , being the maximal rate of R production per unit R present, will be referred to as the ‘ribosomal efficiency’ of the cell. The factor in is the rate at which a ribosome can add an amino acid to a protein when there is no limitation of P and the factor of accounts for the number of P required to make a ribosome. In Scott et al. (2014) is referred to as the ‘translational efficiency’ of the cell.

Optimization as a principle of cellular economy  As mentioned earlier, the growth laws (2) and (3) follow from (19). The latter is a more basic equation as it expresses directly in terms of the parameters without reference to the growth rate, and it encapsulates the consequence of growth rate maximization when . (19) or equivalently (18) can be recast asWe can interpret as the allocation or investment of the cell in the metabolic sector (measured in units of P) and as the investment in the ribosome sector. We define the ‘output’ of each sector as the ‘investment’ times ‘efficiency’ of the sector. Then the investment strategy of the cell, namely (22), can be stated asEquivalently, (22) can be stated as the following principle of cellular economy: the resources allocated to the enzyme and ribosomal sectors are inversely proportional to their respective efficiencies. In other words, the PTR cell follows the dictum: From each sector according to its ability, to each sector according to its need. Here ‘ability’ of a sector is the same as its ‘efficiency’, defined earlier, and ‘need’ is the allocation or investment in the sector that would make its ‘output’ equal to that of the other sector. This principle follows from the optimization of the growth rate of the PTR cell as a whole in the large approximation. Note that efficiency is hardwired into the cellular and medium parameters while the allocation, in the context of the present model, is a matter of cellular ‘choice’ (though, of course, in practice, even that is hardwired into the regulatory mechanisms that dynamically implement the ‘choice’.)

We remark that (22) is not a requirement for the system to have a steady state. Indeed, steady states are achieved in the model even when is not at its optimal value given by (16), as discussed earlier. When , we can still have a steady state with constant concentrations satisfying the Eqs. (8)–(14), but (22) does not hold. (22) is the condition that the steady state has the maximal possible value of given that all parameters other than are fixed.

The Monod curve   We turn to a discussion of the analytic expression for , Eq. (17). First we discuss the situation when . Then from (20), (21), and , and our results for and all the other quantities reproduce exactly the results of Scott et al. (2010, 2014). The growth rate reduces toThis is the same as the expression derived in Scott et al. (2010, 2014), when (20) and  (21) are used to set . To make contact with the Monod equation (1), one has to say how depends upon the concentration [F] of the external nutrient. As mentioned below (4) and hence is an increasing function of [F]. If one substitutes the simplest function , where is a constant, into (24), one obtains (1) with and . Alternatively, if the transport limited Michelis-Menten form of food uptake , where and K are constants, is substituted in (24), one gets (1) with and (Scott et al. 2014).

The difference between our derivation of (24) and that of Scott et al. is that the latter uses the growth laws (2), (3) as the starting point and obtains the above mentioned expression for . It does not require any further assumption of growth rate optimality in deriving that expression as (2), (3) already incorporate optimality. On the other hand, our derivation starts with equations (4) describing the dynamics of the three pools, obtains in the steady state before optimization and then uses the optimality assumption to derive (2), (3) as well as the optimized . This crisply establishes the relationship between optimality and the growth laws.

It may be helpful to make a few remarks about (24). The right-hand side is a symmetric function of and , which define the two natural time scales in the problem. (1) For fixed as a function of , it saturates at a maximum value . The saturation is not a consequence of a Michelis-Menten type saturation kinetics assumed in the model [Eq. (4) has no Michelis-Menten or Hill type terms], but is a consequence of the existence of these two time scales in cellular dynamics. When , in (16) approaches 1; thus the core autocatalytic set that drives the PTR cell—ribosome producing more ribosome—is focused largely on producing itself. Even then, we know that the maximal rate of R self-reproduction production can only be , which explains the saturation. (2) Interestingly, not only is the saturation value of equal to , the value of at which is half its maximum value is also . This has a simple explanation. If, to achieve the maximum growth rate, the ribosome pool is focused solely on making ribosome, then at half the maximal rate only half the pool is focused on making ribosome. The other half is then focused on making T, and this equal investment in both sectors means . But from the principle of cellular economy the two sectorial outputs are equal; therefore, must be equal to . An alternative way of saying this is to observe from (16) that at . (3) The symmetry between and implies that if is held fixed and is increased, will saturate at a value , and the value of at half-saturation is also .

Dissipation terms   Equation (17) is a generalization of (24) when are nonzero. Note that even with the additional terms there is a symmetry between the two sectors: is unchanged under the simultaneous interchange , . When , the factor in the numerator reflects that the metabolic efficiency has to be to sustain a nonzero growth rate. This is because for every molecules of P produced by each molecule of T per unit time, a number is lost through the term. Similarly, when , , the factor in the numerator means that the ribosomal efficiency has to be greater than for the ribosomal ACS to get off the ground. A nonzero requires a greater fraction of ribosomes to be making ribosomes, and a nonzero requires a greater fraction to be making T [see Eq. (16)]. However, the relative investment by the cell in the two sectors as measured by T/R or is independent of [see Eqs. (18), (19)]. The equality of and has been commented upon by Scott et al. (2010). They have considered models in which the degradation terms are zero. In the present model also when . But when , the two are not equal.

We note that in the model the phenomenological parameter is zero if , and the if [Eqs. (20), (21)]. In bacterial cells may be of the order of 0.1 h (Dressaire et al. 2009; Maitra and Dill 2015), while may be much lower (Zundel et al. 2009). This predicts a value of about 2–3 times smaller than the observed value given in Scott et al. (2010). Equation (20) predicts that when , as a function of is linear with a positive intercept. This feature is seen in the data (Scott et al. 2010). However, again the value of the intercept predicted by (20) is smaller than the value from the data. This suggests that other contributions to and , not described by the present model, are significant.

The ‘constant fraction’ sector  Scott et al. (2010) introduced another sector of proteins Q in addition to T and R which takes up a fixed fraction of the protein mass , to account for the fact that was observed in experiments to be less than unity. In the present model this sector can be added as follows (we consider the case ): to (4), add another equation , where . The other changes are in (5), where we now have , and in (6), where a term is added to the definition of V. is now defined by . In the optimization, is treated as a fixed number; can range between 0 and and is chosen to maximize the growth rate. Doing the analysis as for the PTR model, one reproduces the growth laws (2) and (3) in which are the same as for the PTR model, and , . Further, , and, as before, , .

A simpler derivation of the results—from a linear model  Above, we have presented a detailed derivation of , and from the PTR model assuming . It is worth mentioning that the same results follow from a much simpler heuristic argument. Supposing we assume that the dominant contribution to V is , i.e., we ignore the contribution of T and R to V. (This does not mean that the contribution of T and R to the mass of the cell is much smaller than that of P. If , the contribution of T and R to the mass of the cell could be large, even larger than the contribution of P, while their contribution to the volume is much smaller than that of P, as long as are independent of (or sufficiently weakly dependent on) ). Then (4) reduces to a set of linear equations withThe largest eigenvalue of A is b; hence the growth rate of the cell is . The eigenvector corresponding to b has , . Since we have or . We now ask the following: what is the largest value of possible, and for what value of does that occur? Since , one may naively think that the largest possible value of , namely will give the largest . However, we also have the inequality ; therefore, the largest value of occurs when . This is the same conclusion as reached in the “Appendix” for the full PTR model under the large approximation. The condition immediately yields with given by (16), and given by (17). Further the above eigenvector also reproduces (18) for T/R. In the linearized equation , we can recognize the two terms as the outputs of the metabolic and ribosomal sectors.

The above approximation is reasonable for . It is meaningless for because P turns negative in that regime under this approximation, though the full model has a perfectly reasonable behaviour even for . As seen earlier, this approximation is also good for deducing and T/R as these tend to finite limits when approaches from below. It is not useful for estimating P/R near which approaches zero in this approximation. One can see from the full model that P/R receives corrections in a small range of of size around , in which range it goes from a value to a smaller value. In the full model P/R does not go to zero at .

Conclusion

In this paper we have constructed a simple dynamical system describing a cell in terms of its three coarse-grained molecular pools and shown that the optimization of the steady-state growth rate of the cell with respect to a parameter that can be tuned by intracellular regulation leads to the growth laws (1), (2) and (3). We have reproduced and extended existing formulae for the growth rate and other physiological parameters. This deepens our understanding of the macroscopic physiological variables in terms of microscopic parameters. We expect that this kind of model can be extended to include other molecular sectors in the cell (Hui et al. 2015).

At a methodological level we have introduced a scheme that allows an explicit computation of the steady-state growth rate of the cell in terms of the cellular and medium parameters. In this scheme a key assumption is that the volume of the cell is determined by its molecular populations. We have also put to use two natural large parameters in the cell, and , to set the scale of certain other parameters and to make approximations. This has allowed us to get analytic results for the nonlinear system level dynamics.

Our model uses an optimization principle to fix an internal parameter, , the fraction of ribosomes making ribosomes. The model is silent on the dynamical mechanisms inside the cell that implement this optimization. These mechanisms have been the subject of several recent works (Scott et al. 2014; Maitra and Dill 2015; Weiße et al. 2015; Bosdriesz et al. 2015). We hope that combining some of the methods introduced here with the mechanisms discussed in these works will produce models that are more satisfactory than the present one.