A steady-state model of microbial acclimation to substrate limitation

John R. Casey; Michael J. Follows

doi:10.1371/journal.pcbi.1008140

Abstract

Microbes acclimate to changes in substrate availability by altering the number of transporters on the cell surface, however there is some disagreement on just how. We revisit the physics of substrate uptake and consider the steady-state scenario whereby cells have acclimated to maximize fitness. Flux balance analysis of a stoichiometric model of Escherichia coli was used in conjunction with quantitative proteomics data and molecular modeling of membrane transporters to reconcile these opposing views. An emergent feature of the proposed model is a critical substrate concentration S*, which delineates two rate limits. At concentrations above S*, transporter abundance can be regulated to maintain uptake rates as demanded by maximal growth rates, whereas below S*, uptake rates are strictly diffusion limited. In certain scenarios, the proposed model can take on a qualitatively different shape from the familiar hyperbolic kinetics curves, instead resembling the long-forgotten Blackman kinetics.

Author summary

The mechanics of resource-limited microbial growth is a fundamental focus in cell biology and biophysics. Physiological acclimation plays a key role in microbial growth rate dependence on the availability of a limiting resource, but progress has been mostly rooted in theoretical studies due to a lack of relevant experimental data. In light of new quantitative proteomics data which disagree with current models, we revisited the physics of substrate transport and propose a model, based on a different set of assumptions, which applies to the steady-state scenario. Depending on the design of the transport system, the proposed model predicts that microbial growth rate dependence on substrate availability can take on a familiar hyperbolic shape (e.g., Monod) or a piecewise linear shape, not unlike the Blackman kinetics model which fell out of favor long ago. The sharp transition appears from a discontinuity which marks a critical substrate concentration, above which physiological acclimation can sustain maximal growth rates, and below which diffusion is strictly limiting. We implemented the model in Escherichia coli using a combination of flux balance analysis, molecular modeling, and quantitative proteomics data. Predicted kinetics closely matched experimental observations across a range of carbon substrates, and the model can be easily implemented for other systems.

Citation: Casey JR, Follows MJ (2020) A steady-state model of microbial acclimation to substrate limitation. PLoS Comput Biol 16(8): e1008140. https://doi.org/10.1371/journal.pcbi.1008140

Editor: William Cannon, Pacific Northwest National Laboratory, UNITED STATES

Received: March 18, 2020; Accepted: July 9, 2020; Published: August 26, 2020

Copyright: © 2020 Casey, Follows. 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: All relevant data are within the paper and its Supporting Information files. Accompanying code and documentation is available at https://github.com/jrcasey/NutrientUptake.

Funding: This work was supported by the Simons Foundation (https://www.simonsfoundation.org) as part of the Simons Collaboration on Computational Biogeochemical Modeling of Marine Ecosystems (https://cbiomes.org; Simons Foundation grant no. 549931 to M.J.F. and grant no. 549894 to J.R.C). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Microbial growth rates may be limited by catalytic rates or by the rate of diffusion of resources, which are often supplied with some irregularity in the natural environment. As the availability of nutrients and organic substrates change, various metabolic and physiological acclimation strategies can be leveraged to maximize growth rate and remain competitive. For instance, flexibility in the molecular composition of biomass can be used to compensate for resource deficits by adjusting demands through holistic changes in elemental stoichiometry. Given the dynamic nature of microbial ecosystems, these acclimated phenotypes are probably the norm rather than the exception.

An important physiological acclimation strategy to maximize growth on a limiting resource is to alter the number of transporters on the cell surface; however, two conflicting views of how microbes implement this strategy have emerged (Fig 1). One view (“Model A”) predicts that, at very low substrate concentrations, the diffusive flux is low and therefore very few transporters would be required to match uptake rates to encounter rates. As substrate concentrations increase, uptake rates become limited by the number of transporters, so more are synthesized [1–8]. Another view (“Model B”) predicts the opposite; at low substrate concentrations the encounter rate is limiting, so high transporter abundance is favored, while at higher concentrations resources are allocated for use elsewhere in the proteome [9–13]. Considering the evolution of these two divergent views, it is interesting that both evolved under different assumptions from the same underlying physics [14].

Download:

Fig 1. Observed glucose transporter abundances across a range of glucose concentrations ([15]; data plotted with permission from the authors).

Qualitative predictions of Models A and B are shown to illustrate their expected behavior.

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

The relevant data (i.e., the abundance of transporters across a range of substrate concentrations) which could resolve, or even reconcile these opposing views has, until relatively recently, been lacking; however, progress in proteomics has enabled the direct quantitation of the majority of expressed proteins in model organisms like Escherichia coli K12 [15], including its transporters. In glucose fed batch cultures and chemostats across a range of dilution rates, transporter abundances qualitatively matched neither of the two previous model predictions (Fig 1). This finding suggests perhaps that another interpretation of the underlying physics be posed which explains these observations. For instance, prior approaches in line with Model B have leaned heavily on resource allocation, whereby a tradeoff is established by the fraction of the proteome dedicated to substrate transport [8–12]. This fraction was assumed to modulate both the affinity and the maximum uptake rate, often requiring large proteome investments in transport (e.g., up to 77% of total protein; [8]), however the total contribution of all 90 transporter proteins on both the inner and outer membranes and their associated periplasmic binding proteins to the Escherichia coli K12 proteome was 11±3% across 22 different growth conditions, with any particular transporter complex contributing much less [15]. Clearly the allocation of resources to transport is not without cost, as we will explore, but may offer confounding predictions when modeled as a linear constraint on kinetic parameters.

We consider the simplest case of a single cell growing on a single limiting substrate S (Fig 2). The rate of change of the cell quota Q of element X (mol X cell^-1) can be expressed in terms of S as (1) where v is the substrate uptake rate (mol S cell^-1 s^-1, r is the cell radius (m), n is the number of transporters (cell^-1), D is the molecular diffusivity of the substrate (m² s^-1), Y_X/S is the yield (mol X [mol S]^-1), and M is a non-growth associated maintenance requirement (s^-1). Over the time-average in a continuous culture, , so we define the steady-state growth rate μ (s^-1) as (2)

Download:

Fig 2. Conceptual diagram of substrate uptake.

A substrate molecule diffuses down a gradient in concentration from S_∞ toward the cell surface S₀. The molecule encounters an unoccupied transporter at a rate no higher than the catalytic rate, where it is subsequently transported to the interior S_I and metabolized to generate reducing power or to synthesize new biomass. Coupling of a mechanistic model of substrate transport with flux balance analysis of a stoichiometric model of the metabolic network allows computation of yields, optimal uptake rates, and optimal transporter abundances.

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

The experimental data of Shmidt et al., [15] provide the opportunity to test and calibrate a model of nutrient uptake and transporter allocation since r, n, S, D, and μ are simultaneously known. Here we base the functional form of v on the model of Armstrong [4], an approximation to the model of Pasciak and Gavis [2] (described in detail in the Model section). To close the model according to Eq (2) we must also constrain the yield Y_X/S and maintenance M; we do so by employing a genome-scale metabolic model of Escherichia coli K12 [16].

In short, we will demonstrate that the model is compatible with the quantitative proteomic data of Schmidt et al., [15] under the assumption that allocation to transporters is optimized to two, theoretically anticipated regimes: (i) a regime where internal synthesis rates are the limiting factor and (ii) a regime where diffusive encounter with the resource is the limiting factor.

Results

Conceptual model

For simplicity, we consider a single cell of radius r suspended in an infinitely large volume, containing a single nutrient at bulk concentration S_∞ (mol m^-3; Fig 2; Table 1). The subscript ∞ denotes that the concentration is that of a volume an infinite distance from the cell surface. A crude description of transport of an extracellular substrate at a concentration at the surface of the cell S₀ (mol m^-3) across the inner membrane via a transporter E to the intracellular space at a concentration S_i can be expressed as the two-step reaction (3) where k₁ is the molecular encounter rate with a transporter (mol m^-3 s^-1) and k_cat is the catalytic constant or turnover number, the rate of dissociation of the transporter-substrate complex ES in the transport process (mol transporter^-1 s^-1). It should be noted that k₁ represents not only the encounter rate, but specifically the encounter rate of molecules colliding with the transporter with sufficient energy to overcome a reaction energy barrier. Thus, the back reaction, often denoted k₋₁, is implicit. Also note that k_cat should be interpreted to be the maximum apparent catalytic rate in vivo, rather than the constitutive in vitro k_cat [17].

Download:

Table 1. Parameters used in this article.

https://doi.org/10.1371/journal.pcbi.1008140.t001

At high substrate concentrations (S₀ → ∞) when the molecular encounter rate with each transporter far exceeds their maximum catalytic rate (k₁ ≫ k_cat), effectively all transporters are occupied, that is, in the process of carrying out a transport reaction. For consistency with previous studies, we refer to this limit as the porter limit, denoted by the superscript P. At the porter limit, the maximum uptake rate is set by the catalytic rate of a number of transporters n, (4)

However, another limit is reached as transporter abundances increase further; we refer to this limit condition as the growth limit , denoted here and elsewhere by the superscript G. In this state, increased uptake capacity would exceed metabolic demands and thus the excess transporters would be proportionately inhibited, so as not to violate . The maximum growth rate under optimal conditions can be set by a variety of internal bottlenecks (e.g., protein translation rates, oxygen diffusion, metabolic choke-points, molecular crowding, temperature), which are manifested by differences in metabolic and physiological designs in a particular environment. The growth limit represents a specific instance of v_max where cells have acclimated to growth at high nutrient concentrations by adjusting the number of transporters to some optimal abundance n* = n^G which supplies substrate at the rate required to satisfy the maximum growth rate. In this instance, the porter limit and growth limit converge.

Kinetic model

A model of the dependence of uptake rate kinetics on changes in cell physiology was previously developed [2]. In subsequent work, Armstrong [4] provided a convenient approximation of Pasciak and Gavis’ quadratic model, resembling an expansion of the hyperbolic Michaelis-Menten model . In this approximation, the apparent half-saturation concentration k_s = k_cat/k₁ can be described by the sum of two limits: a porter limit and a diffusion limit . Armstrong’s model incorporates similar dynamics to the quadratic model, but uses the simplifying assumption that cell radii are much smaller than, or much larger than a characteristic length scale of the reaction-diffusion process . A discussion on the validity of this assumption can be found elsewhere [6]. Thus, we apply Armstrong’s approximation (5) where is the effective half-saturation concentration at the porter limit, which is independent of the number of transporters, (6) Critically, α represents a dimensionless probability that a substrate molecule which enters the vicinity of the transporter catchment area A (m²) is captured and transported. This parameter indirectly accounts for the fraction of collisions which exceed the activation energy of the ligand binding reaction, without any explicit knowledge of the magnitude of this barrier or the energy of the collision, (7) We derive this probability α by equating the velocity (m s^-1) of molecular transport across the membrane when all transporters are saturated (ν^P at the porter limit; [4]) to the velocity of the nutrient molecule diffusing towards the cell surface (ν^D; [14]), which is dependent on the hydrated molecular diffusivity of the substrate D (m² s^-1), the cell size, and the advective velocity of the cell u (m s^-1), (8)

Thus, motility is inversely proportional to the effective half-saturation concentration at the porter limit, by increasing the encounter rate k₁ relative to the catalytic rate, as can be seen from Eqs (6), (7) and (8). The effective half-saturation concentration at the diffusion limit is not only cell size-dependent, but is dependent on the number of transporters, such that (9) where ϕ is a dimensionless cell shape factor and Sh is the dimensionless Sherwood number, which relates mass transfer by advective shear forces to those of the viscous forces. We include ϕ and Sh for completeness, but have excluded any explicit treatment of their effects on nutrient transport in this article by assigning both of their values to be 1. Both turbulence and cell shape influence the encounter rate, and a detailed discussion on their effects on nutrient transport can be found elsewhere [18, 19].

Acclimation

Afforded sufficient time, microbial isolates grown under viable conditions will acclimate both metabolically and physiologically so as to maximize growth rate (e.g., [20]). Given this observation, it follows that uptake rates should be matched as closely as possible to by regulating transporter abundances optimally. Maintenance of a number of transporters in excess of this optimal value would result in the intracellular accumulation of the substrate and allosteric feedback inhibition of metabolism and sub-optimal growth. Conversely, maintenance of an insufficient number of transporters would result in porter limitation, which is also a sub-optimal growth state. If both of these statements are true, then an optimal number of transporters n* should be predicted over two concentration intervals; a zero-order interval over which uptake rates may be maintained at , and a first-order interval over which uptake is limited by the substrate concentration dependent diffusive flux. A schematic of these transitions is given in panels A and C in Fig 3. At zero-order, (10)

Download:

Fig 3. Schematic of steady-state nutrient acclimation.

A)—Optimal transporter abundances n* (solid green line) lie at the porter limitation (blue shaded region) boundary. All points above and below this boundary fail to maximize growth rate; for all points below, the uptake rate is sub-maximal; for all points above, some transporters are unoccupied. Over an interval of bulk substrate concentrations S_∞ ≤ S*, this boundary is met with diffusion limitation (orange shaded region), where the catalytic rate exceeds the encounter rate; as concentrations exceed S*, the porter limitation boundary transitions to an internal growth rate limit (yellow shaded region), where the catalytic rate exceeds the rate of some downstream reaction. B)—For smaller cells, or for transporters with slower k_cat, membrane surface area limitation (a special case of porter limitation) may be encountered. In the interval bounded by and , n* is constrained to n_max (black dashed line). C)—Growth rates of optimally acclimated cells (solid green line) follow the intersection of two limits: the diffusive limit (purple dashed line) and the internal growth limit (black dashed line). These rate limits are, again, bisected by S* with a sharp transition. D)—If surface area limitation is encountered, growth rates in the interval follow a more gradual, hyperbolic transition from diffusion limitation to growth limitation.

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

By demanding that the uptake rate match , we find the optimal number of transporters to be (11)

In this expression, behaves proportionally to 1/S_∞ and it is plainly seen that a discontinuity is possible at a particular S_∞. This concentration represents the lower limit for which Eq 11 holds, (12)

For all concentrations , no abundance of transporters is sufficient to match uptake to the growth limit. Instead, the maximal uptake rate is set by the diffusive flux v^D = 4πϕShrDS_∞. Accordingly, the optimal number of transporters in this diffusion limited interval is (13)

In keeping with the maximum growth rate assumption, it follows that maintaining unoccupied transporters incurs some non-zero cost; we therefore consider the minimum n* over the full domain of S_∞. The intersection of the two transporter abundance optima S* represents the steady-state transition from diffusion limitation to growth limitation. Since is a linear function of S_∞, and since is symmetric about the discontinuity , the intersection S* (in the positive domain of S_∞) is the positive root of a quadratic, (14)

By assuming optimal transporter abundance within the cellular growth context of Eq 2, we arrive at the proposed steady-state acclimation model (Fig 3), which is piecewise linear with a transition at the discontinuity S*, (15)

Cell size and surface area limitation

It is worth noting that although cell size is clearly an important determinant of growth rate in the diffusion limited regime, the range of cell sizes observed within the range of phenotypes from batch-acclimated to the lowest dilution rate glucose-limited chemostat spanned only a 12% change in radius [15]. Thus we expect acclimation to nutrient availability in this case to be predominantly regulated by transporter abundance. It might also be noted that despite this small observed change in cell sizes over the substrate concentration interval tested, the most effective strategy to maximizing growth rate as S_∞ → 0 shifts from regulating n to regulating r. Perhaps Escherichia coli K12 is simply not capable of this adaptation, but should not be ignored in the case of other taxa which are.

As concentrations approach in the nutrient-replete regime, n* → ∞, but clearly this is not physiologically possible since there is finite membrane space afforded to transporters. In their conceptual model of transport through a sphere of known dimension and with a number n of pores of known dimension, Berg and Purcell [1] proposed that the addition of a number of transporters , where s is the radius of the transporter, does not appreciably increase the uptake rate [1]. Given their selection of 10 Å for transporter radius and for a cell of 1 μm radius, this corresponds to an areal coverage of less than 0.1%. Although this is a convenient parameterization for the maximum number of transporters n_max, we find that in the case of the glucose transporter PtsG, given experimental maximum n = 9402 cell^-1 [15] and corresponding cell size, s would need to be no larger than 2 Å which is roughly half the van der Waals radius of the glucose molecule alone, and much smaller than the barrel radius from our structural analysis (see Methods). On the contrary, Asknes and Egge [3] predicted a site coverage of 8.5%, corresponding to n_max = 21, 058 cell^-1, which is comfortably above the maximum observed n given our best estimates of the corresponding catchment area for the PtsG protein. Therefore, we required that all n ≤ n_max, where and the maximal areal coverage f_max = 0.085. Nevertheless, the wide discrepancy between f_max estimates between these studies is unclear, and it is certainly possible that f_max will vary between cell membrane designs, and is therefore a somewhat unsatisfying and arbitrary boundary which should be better constrained.

Variation in cell size introduces a nuance to our proposed model which provides a mechanistic basis for more gradual transitions observed between the diffusion-limited regime and the growth-limited regime. In some scenarios, n* exceeds n_max over a substrate concentration interval wherein uptake rates become surface area (SA) limited, a special case of porter limitation. This interval can be constrained to a lower and upper bound , both evaluated at the intersections of Eqs 11 and 13 with n_max, (16) (17) Within the interval , the optimal number of transporters is constrained to n_max, resulting in a hyperbolic transition in uptake rate from diffusion limitation to growth limitation (panels B and D in Fig 3). Although it was not relevant for our present case, to accommodate surface area limitation under such a scenario, a third conditional is appended to Eq 15, (18) where v(n_max) is the uptake rate defined in Eq 5, evaluated at n = n_max.

Model validation

Glucose transport in aerobic cultures of the model bacterium Escherichia coli K12 was used as a system to validate our acclimation model. Although the physiological response to glucose limitation in E. coli is complex, with several complementary transporters with broad substrate specificity induced by cAMP, the primary transporter under both glucose excess and glucose limited growth conditions is the phosphotransferase system (PTS) which relies on the glucose-specific permease PtsG [27]. Indeed, PtsG was the most abundant permease in both glucose-limited chemostats and glucose excess batch cultures (S1 Fig; [15]).

Whereas k_cat values are rarely reported for transporters, we were able to directly compare k_s values from a collection of 11 substrates (Table 2). The modeled effective half-saturation concentrations () for all substrates spanned 4 orders of magnitude between substrates. k_s values from experimental studies were compiled from the same Escherichia coli strain K12 (but not necessarily from the same sub-strain), harvested in exponential phase from substrate replete batch cultures. In a direct comparison, our model predictions agreed closely with these reported values (Fig 4; Model 2 linear regression, R² = 0.89, dF = 11), and the slope was not different from parity (regression slope = 0.98, standard error = 0.11). Note that the literature value reported for zinc was determined by a method which gives a value closer to the dissociation constant rather than the half-saturation constant [28], and should therefore be considered an upper bound. The deviation between predicted and observed k_s values for acetate transport is less obvious, but may be due to either the mechanism of the transporter ActP, which is the only symporter in our set, or the presence of a second, lower-affinity symporter YaaH which would introduce biphasic kinetics that our model does not currently resolve.

Download:

Table 2. Summary of parameter values corresponding to batch growth conditions.

Sulfate, phosphate, and zinc parameters are reported only for glucose-replete batch cultures. “Exp.” refers to experimentally determined k_s values, except in the case of zinc, which should be considered a dissociation concentration and is thus an upper bound.

https://doi.org/10.1371/journal.pcbi.1008140.t002

Download:

Fig 4. Comparison of predicted and experimental values of the effective half-saturation concentrations of 11 substrates in nutrient-replete batch acclimated cultures of Escherichia coli K12.

Parity is indicated as a solid line.

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

Our core assumption, that unoccupied transporters would not be maintained by acclimated phenotypes which have optimized their physiology to maximize growth rates, can be tested by comparing transporter abundances for non-limiting substrates in cells grown across a range of growth conditions, with the expectation that abundances scale with maximal uptake rates. Additionally, at constant temperature, pH, and salinity, the property k_cat would be expected to be invariant. Both the absolute abundance (cell^-1) and the aerial density (m^-2 of inner membrane surface area) of transporters of the inorganic nutrient ions phosphate, sulfate, and zinc, which were supplied in great excess in all media formulations [15], were linearly related to their uptake rates, as quantified by FBA, across all growth conditions (P = 9 * 10⁻⁵, 5 * 10⁻³, and 6 * 10⁻³, respectively). The coefficient of variation between growth conditions for k_cat was less than 25% for this subset, and no growth rate or cell size dependence of the residuals could be detected.

Acclimation to limiting nutrients

Escherichia coli K12 BW25113 cells maintained an optimal number of PtsG transporters corresponding to Eq 11 for all S_∞ > S*, and corresponding to Eq (13) for all S_∞ ≤ S*. Fig 5 shows contours of uptake rates in the plane of transporter abundance and glucose concentrations. The instantaneous uptake rate kinetics for an arbitrary acclimated phenotype may be described by the hyperbolic curve generated by following a horizontal section in the top panel, corresponding to a constant number of transporters. Accordingly, the uptake rate kinetics for all optimally acclimated phenotypes may be described by following n*, which is a linear function proportional to S_∞ in the diffusion-limited regime, and a reciprocal function of S_∞ in the growth limited regime. Predicted n* optima were in close agreement with measured PtsG glucose transporter abundances for Escherichia coli cultures grown in glucose-limited chemostats spanning a range of dilution rates from 0.12 h^-1 to 0.50 h^-1 and in nutrient replete batch acclimated cultures growing in exponential phase at 0.60 h^-1. In this particular scenario, n_max was not intersected by n*, so no intermediate surface area limitation transition was encountered (Fig 5—Top panel; [15]). Uptake rates for predicted optimally acclimated phenotypes also closely matched those derived from data (Fig 5—Bottom panel; [15]). The uptake rate kinetics for all optimally acclimated phenotypes was linear in the diffusion limited regime, and transitioned at S* to maintain in the growth limited regime.

Download:

Fig 5. Model predictions and observations of Escherichia coli acclimation to growth on glucose.

Top panel—Modeled and experimental abundances of the glucose transporter PtsG across steady-state concentrations of glucose (log scale). Contours indicate the corresponding uptake rates. The vertical green dashed line represents the critical substrate concentration S*. n_max is above the plotted range. Bottom panel—Modeled and experimental uptake rates over an interval of glucose concentrations. The uptake rate profile for a glucose-replete batch acclimated culture is shown in cyan. Uptake rates for all acclimated phenotypes over the concentration range are shown in red, with model predictions shown as a continuous line and experimentally derived values from published data [15] shown as markers. For guidance, the maximum diffusive flux is shown as a dashed blue line, S* is shown as a vertical dashed green line, and is indicated by the horizontal dashed orange line.

https://doi.org/10.1371/journal.pcbi.1008140.g005

Discussion

A model of substrate uptake which accounts for optimal physiological acclimation to sustained growth at a range of constant substrate concentrations was developed and implemented with a pseudo-mechanistic approach. Critically, the proposed model assumes that an optimal phenotype would maintain only the minimal number of transporters required to sustain either the maximum growth rate or the maximum diffusive flux. Maintaining more than this optimal abundance would incur a cost by exceeding substrate demand, and thus a growth rate penalty; maintaining fewer than this optimal abundance would result in sub-maximal substrate uptake rates, also a growth rate penalty. The resulting profile of both observed and optimal transporter abundances (Fig 5) shows an intermediate maximum, which is unexpected by conventional wisdom of, at least, the transcriptional response to nutrient limitation. The gap in observations of transporter abundances between the highest dilution rate glucose-limited chemostat and the glucose-replete batch acclimated culture is expected for practical reasons, but is somewhat unsatisfying and should be revisited to further validate the proposed model as quantitative proteomics datasets become more commonplace. Nevertheless, the steady-state kinetics associated with this optimal set of phenotypes (Eq 15) resembles those of the piecewise linear Blackman kinetics, and for similar reasons; a transition from one limitation to another. The critical substrate concentration S* therefore represents the transition between diffusion limitation and growth limitation, barring the scenario where membrane surface area limitation mediates that transition. Variations in cell size or k_cat introduce some nuance to the shape of this transition, as membrane surface area limitation presents another, intermediate limiting state. In those scenarios, n* is held at n_max, resulting in hyperbolic uptake rates in the interval between and , thus providing a mechanism to reconcile observations of a more gradual transition from first-order diffusion limitation to zero-order growth limitation (Eq 18). Although we have considered the simple case of a single limiting substrate, we speculate that the competition between transporters for membrane space may give rise to apparent co-limitation when multiple substrates are present at concentrations near S*, as may very well be the case for natural environments where S* would be expected in the vicinity of the minimum resource limit for net population growth (Tilman’s R*) [38].

Steady-state substrate-limited microbial growth has been routinely modeled as a hyperbolic function of the concentration of a limiting resource either internally to the cell, as is the case for the Droop model [39], or in the external medium, as is the case for the Monod model [40]. Discussion of the well known Monod and Droop models can be found elsewhere (e.g., [41]), but the Blackman model [42] is a somewhat forgotten model of growth. The Blackman model describes steady-state microbial growth kinetics by a linear dependence on a limiting factor, similar in form to (Eq 15), which persists over a discrete interval until another limiting factor is encountered. Despite their differences, critical evaluation of these foundational microbial growth models and their later modifications (e.g., [43, 44]) have shown support for each [45]. Intriguingly, Blackman kinetics, which were presumably known to Monod some 40 years later, would have provided a better fit to his chemostat data from glucose-limited Escherichia coli than did his own empirical model. This realization was the motivation behind several subsequent studies which expanded on Blackman kinetics [46, 47]. Piecewise linear growth dependence on substrate concentration is readily observed in natural microbial communities (e.g., Fig 4 in: [41]), though it is rarely attributed as such (c.f., [48]) and has been largely abandoned in favor of the more mathematically convenient hyperbolic models of Monod and Droop. The competition for membrane space in our proposed model provides a possible mechanism to reconcile these two fundamentally different descriptions of microbial growth, one which could conceivably be put to the test.

The distinction between the assumptions that frame our model from those of the two previously mentioned Models A and B on the dependence of transporter abundance on substrate concentration is fairly subtle, but can predict qualitatively different shapes to the growth dependence on substrate availability. As in vivo enzyme turnover number predictions improve and become more widely available (e.g., using machine learning algorithms [43]), the requirement for technically challenging and costly quantitative proteomics data may be obviated. This should enable more comprehensive validation of the model and more broad application to steady-state simulations of microbial growth, requiring only that a genome-scale stoichiometric model of metabolism be available, and that the cell size and the maximum growth rate be additionally known.

Methods

Determination of transporter k_cat by quantitative proteomics and flux balance analysis

k_cat values were determined by Eq 4 using flux balance analysis (FBA) quantitation of and corresponding protein abundances n^G reported in a quantitative proteomics study [15]. A genome-scale stoichiometric model of metabolism for Escherichia coli K12 MG1655 (iML1515; [16]) was used to quantify transport fluxes in simulations of each of the cultivation conditions employed by Schmidt et al. [15] in their quantitative proteomics study.

In their study, Escherichia coli K12 BW25113 was grown in batch culture on each of 11 sole carbon sources and harvested in logarithmic phase, or in continuous culture on glucose supplied at 4 dilution rates (0.12, 0.20, 0.35, and 0.50 h^-1). FBA was implemented by allowing for unconstrained transport of only those substrates which were supplied in each defined medium, and the “biomass reaction” was constrained to the experimental growth rate. An L1-norm minimization was implemented to remove loops. Extensive documentation and guides to the implementation of FBA with iML1515 and other stoichiometric models of metabolism within several programming environments is available at opencobra.github.io. FBA was implemented within the Matlab (The Mathworks, Inc.) Toolbox COBRA (Version 3.0; [21]), the Python (Python Software Foundation) package COBRApy [22], and optimizations were performed using Mosek (Version 9; Mosek ApS).

Determination of transporter capture area

Ligand capture area A for each transporter was interpreted to be the area of the transporter membrane domain. Protein sequences for each transporter were used to model quaternary structures using RaptorX [23]. The predicted structures were then rendered and transmembrane domain dimensions were measured using PyMOL (The PyMOL Molecular Graphics System, Version 2.3 Schrodinger, LLC).

Estimation of molecular diffusivities

Aqueous diffusion coefficients were determined as an empirical function of hydrated molecular volumes and water viscosity [24]. Hydrated molecular volumes were calculated using the LeBas incremental method [25]. Dynamic viscosity was calculated as a function of temperature and salinity [26] of the cultivation conditions [15].

Physiological data

To their credit, Schmidt and co-authors [15] had the foresight to supply detailed and complementary physiological measurements with their proteomics dataset. Electron micrographs provided measurements of inner and outer cell membrane dimensions, and cell concentrations and growth rates were reported with each medium formulation. We commend the authors on providing these data in an accessible format, as it has enabled our model validation and development. It is worth noting that although a different Escherichia coli K12 substrain (BW25113) was used for the majority of growth conditions to quantify protein abundances from the metabolic model substrain (MG1655) we used in FBA simulations, Schmidt et al., [15] showed only minor differences in absolute protein quantitation between these two substrains in glucose-replete batch cultures.

Supporting information

S1 Fig. Comparison of the abundances of permease domains of known glucose transport systems in Escherichia coli K12 grown in glucose-replete batch cultures and glucose-limited chemostats, as determined by quantitative proteomics ([15]; Data plotted with permission from the authors).

MalF—maltose permease; ManY—mannose permease; MglC—galactose permease; PtsG—glucose permease.

https://doi.org/10.1371/journal.pcbi.1008140.s001

(TIF)

References

1. Berg HC, Purcell EM. Physics of chemoreception. Biophysical Journal. 1977;20:193–219. pmid:911982
- View Article
- PubMed/NCBI
- Google Scholar
2. Pasciak WJ, Gavis J. Transport limitation of nutrient uptake in phytoplankton. Limnology and Oceanography. 1974;19(6):881–888.
- View Article
- Google Scholar
3. Aksnes D, Egge J. A theoretical model for nutrient uptake in phytoplankton. Marine Ecology Progress Series. 1991;70:65–72.
- View Article
- Google Scholar
4. Armstrong RA. Nutrient uptake rate as a function of cell size and surface transporter density: A Michaelis-like approximation to the model of Pasciak and Gavis. Deep Sea Research Part I: Oceanographic Research Papers. 2008;55(10):1311–1317.
- View Article
- Google Scholar
5. Aksnes D, Cao F. Inherent and apparent traits in microbial nutrient uptake. Marine Ecology Progress Series. 2011;440:41–51.
- View Article
- Google Scholar
6. Fiksen Ø, Follows MJ, Aksnes DL. Trait-based models of nutrient uptake in microbes extend the Michaelis-Menten framework. Limnology and Oceanography. 2013;58(1):193–202.
- View Article
- Google Scholar
7. Lindemann C,Fiksen Ø, Andersen KH, Aksnes DL. Scaling laws in phytoplankton nutrient uptake affinity. Frontiers in Marine Science. 2016;3.
- View Article
- Google Scholar
8. Zeng H, Yang A. Bridging substrate intake kinetics and bacterial growth phenotypes with flux balance analysis incorporating proteome allocation. Scientific Reports. 2020;10(1).
- View Article
- Google Scholar
9. Pahlow M. Linking chlorophyll-nutrient dynamics to the Redfield N:C ratio with a model of optimal phytoplankton growth. Marine Ecology Progress Series. 2005;287:33–43.
- View Article
- Google Scholar
10. Smith SL, Yamanaka Y. Optimization-based model of multinutrient uptake kinetics. Limnology and Oceanography. 2007;53(4):1545–1558.
- View Article
- Google Scholar
11. Smith S, Yamanaka Y, Pahlow M, Oschlies A. Optimal uptake kinetics: physiological acclimation explains the pattern of nitrate uptake by phytoplankton in the ocean. Marine Ecology Progress Series. 2009;384:1–12.
- View Article
- Google Scholar
12. Smith SL, Pahlow M, Merico A, Wirtz KW. Optimality-based modeling of planktonic organisms. Limnology and Oceanography. 2011;56(6):2080–2094.
- View Article
- Google Scholar
13. Bonachela JA, Raghib M, Levin SA. Dynamic model of flexible phytoplankton nutrient uptake. Proceedings of the National Academy of Sciences of the United States of America. 2011;108(51):20633–20638. pmid:22143781
- View Article
- PubMed/NCBI
- Google Scholar
14. Munk WH, Riley GA. Absorption of nutrients by aquatic plants. Journal of Marine Research. 1952;11:215–240.
- View Article
- Google Scholar
15. Schmidt A, Kochanowski K, Vedelaar S, Ahrné E, Volkmer B, Callipo L, et al. The quantitative and condition-dependent Escherichia coli proteome. Nature Biotechnology. 2015;34(1):104–110. pmid:26641532
- View Article
- PubMed/NCBI
- Google Scholar
16. Monk JM, Charusanti P, Aziz RK, Lerman JA, Premyodhin N, Orth JD, et al. Genome-scale metabolic reconstructions of multiple Escherichia coli strains highlight strain-specific adaptations to nutritional environments. Proceedings of the National Academy of Sciences. 2013;110(50):20338–20343.
- View Article
- Google Scholar
17. Davidi D, Noor E, Liebermeister W, Bar-Even A, Flamholz A, Tummler K, et al. Global characterization of in vivo enzyme catalytic rates and their correspondence to in vitro k_cat measurements. Proceedings of the National Academy of Sciences. 2016;113(12):3401–3406.
- View Article
- Google Scholar
18. Gavis J. Munk and Riley revisited: nutrient diffusion transport and rates of phytoplankton growth. Journal of Plankton Research. 1976;34(2):161–179.
- View Article
- Google Scholar
19. Lazier JRN, Mann KH. Turbulence and the diffusive layers around small organisms. Deep-Sea Research. 1989;36(11):1721–1733.
- View Article
- Google Scholar
20. Lenski RE, Rose MR, Simpson SC, Tadler SC. Long-term experimental evolution in Escherichia coli. I. Adaptation and divergence during 2,000 generations. American Naturalist. 1991;138(6):1315–1341.
- View Article
- Google Scholar
21. Heirendt L, Arreckx S, Pfau T, Mendoza SN, Richelle A, Heinken A, et al. Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nature Protocols. 2019;14(3):639–702. pmid:30787451
- View Article
- PubMed/NCBI
- Google Scholar
22. Ebrahim A, Lerman JA, Palsson BØ, Hyduke DR. COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Systems Biology. 2013;7:74. pmid:23927696
- View Article
- PubMed/NCBI
- Google Scholar
23. Källberg M, Wang H, Wang S, Peng J, Wang Z, Lu H, et al. Template-based protein structure modeling using the RaptorX web server. Nature Protocols. 2012;7(8):1511–1522. pmid:22814390
- View Article
- PubMed/NCBI
- Google Scholar
24. Hayduk W, Laudie H. Prediction of diffusion coefficients for nonelectrolytes in dilute aqueous solutions. AIChE Journal. 1974;20:611–615.
- View Article
- Google Scholar
25. Arnold JH. Studies in diffusion, II Kinetic theory of diffusion in liquid systems. Journal of the American Chemical Society. 1930;52:3937–3948.
- View Article
- Google Scholar
26. El-Dessouky HT, Ettouney HM. Fundamentals of Salt Water Desalination. 1st ed. Elsevier Science; 2002.
27. Steinsiek S, Bettenbrock K. Glucose transport in Escherichia coli mutant strains with defects in sugar transport systems. Journal of Bacteriology. 2012;194(21):5897–5908. pmid:22923596
- View Article
- PubMed/NCBI
- Google Scholar
28. Graham AI, Hunt S, Stokes SL, Bramall N, Bunch J, Cox AG, et al. Severe zinc depletion of Escherichia coli. Journal of Biological Chemistry. 2009;284(27):18377–18389. pmid:19377097
- View Article
- PubMed/NCBI
- Google Scholar
29. Lindblow-Kull C, Kull FJ, Shrift A. Single transporter for sulfate, selenate and selenite in Escherichia coli K12. Journal of Bacteriology. 1985;163(3):1267–1269.
- View Article
- Google Scholar
30. Willsky GR, Bennett RL, Malamy MH. Inorganic phosphate transport in Escherichia coli: Involvement of two genes which play a role in alkaline phosphatase regulation. Journal of Bacteriology. 1973;113(2):529–539. pmid:4570598
- View Article
- PubMed/NCBI
- Google Scholar
31. Gimenez R, Nunez MF, Badia J, Aguilar J, Baldoma L. The Gene yjcG, cotranscribed with the gene acs, encodes an acetate permease in Escherichia coli. Journal of Bacteriology. 2003;185(21):6448–6455. pmid:14563880
- View Article
- PubMed/NCBI
- Google Scholar
32. Charbit A, Reizer J, Saier MH Jr. Function of the duplicated IIB domain and oligomeric structure of the fructose permease of Escherichia coli. Journal of Biological Chemistry. 1996;271(17):9997–10003. pmid:8626640
- View Article
- PubMed/NCBI
- Google Scholar
33. Engel P, Kramer R, Unden G. Anaerobic fumarate transport in Escherichia coli by an fnr-dependent dicarboxylate uptake system which is different from the aerobic dicarboxylate uptake system. Journal of Bacteriology. 1992;174(17):5533–5539. pmid:1512189
- View Article
- PubMed/NCBI
- Google Scholar
34. McDonald TP, Walmsley AR, Henderson PJF. Asparagine 394 in putative helix 11 of the galactose-H+ symport protein (GalP) from Escherichia coli is associated with the internal binding site for cytochalasin B and sugar. Journal of Biological Chemistry. 1997;272(24):15189–15199. pmid:9182541
- View Article
- PubMed/NCBI
- Google Scholar
35. Hunter IS, Kornberg HL. Glucose transport of Escherichia coli growing in glucose-limited continuous culture. Biochemical Journal. 1979;178:97–101. pmid:373752
- View Article
- PubMed/NCBI
- Google Scholar
36. Truniger V, Boos W. Glycerol uptake in Escherichia coli is sensitive to membrane lipid composition. Research in Microbiology. 1993;144:565–574. pmid:8310182
- View Article
- PubMed/NCBI
- Google Scholar
37. Shamanna DK, Sanderson KE. Uptake and catabolism of D-Xylose in Salmonella typhimurium LT2. Journal of Bacteriology. 1979;139(1):64–70. pmid:222731
- View Article
- PubMed/NCBI
- Google Scholar
38. Tilman D. The Resource-Ratio Hypothesis of plant succession. The American Naturalist. 1985;125(6):827–852.
- View Article
- Google Scholar
39. Droop MR. Some thoughts on nutrient limitation in algae. Journal of Phycology. 1973;9:264–272.
- View Article
- Google Scholar
40. Monod J. The growth of bacterial cultures. Annual Reviews Microbiology. 1949;3:371–394.
- View Article
- Google Scholar
41. Sommer U. A comparison of the Droop and the Monod models of nutrient limited growth applied to natural populations of phytoplankton. Functional Ecology. 1991;5(4):535–544.
- View Article
- Google Scholar
42. Blackman FF. Optima and limiting factors. Annals of Botany. 1905;19(74):281–295.
- View Article
- Google Scholar
43. Powell EO. The growth rate of microorganisms as a function of substrate concentration. In: Powell PO, Evans CGT, Strange RE, Tempest DW, editors. Microbial physiology and continuous culture. H.M.S.O. London; 1967. p. 34–56.
44. Westerhoff HV, Lolkema LS, Otto R, Hellingwerf KJ. Thermodynamics of growth: Non-equilibrium thermodynamics of bacterial growth: The phenomenological and the mosaic approach. Bioehimica et Biophysica Acta. 1982;683:181–220.
- View Article
- Google Scholar
45. Jannasch HW, Egli T. Microbial growth kinetics: a historical perspective. Antonie van Leeuwnhoek. 1993;63:213–224.
- View Article
- Google Scholar
46. Condrey RE. The chemostat and Blackman kinetics. Biotechnology and Bioengineering. 1982;24:1705–1709. pmid:18546467
- View Article
- PubMed/NCBI
- Google Scholar
47. Dabes JN, Finn RK, Wilke CR. Equations of substrate-limited growth: The case for Blackman kinetics. Biotechnology and Bioengineering. 1973;15:1159–1177. pmid:4203180
- View Article
- PubMed/NCBI
- Google Scholar
48. Riley GA. Factors controlling phytoplankton populations on Georges Bank. Journal of Marine Research. 1946;6:54–73.
- View Article
- Google Scholar
49. Heckmann D, Lloyd CJ, Mih N, Ha Y, Zielinski DC, Haiman ZB, et al. Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models. Nature Communications. 2018;9(1):1–10.
- View Article
- Google Scholar

[ref1] 1. Berg HC, Purcell EM. Physics of chemoreception. Biophysical Journal. 1977;20:193–219. pmid:911982
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Pasciak WJ, Gavis J. Transport limitation of nutrient uptake in phytoplankton. Limnology and Oceanography. 1974;19(6):881–888.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref3] 3. Aksnes D, Egge J. A theoretical model for nutrient uptake in phytoplankton. Marine Ecology Progress Series. 1991;70:65–72.
View Article
Google Scholar

[9] View Article

[10] Google Scholar

[ref4] 4. Armstrong RA. Nutrient uptake rate as a function of cell size and surface transporter density: A Michaelis-like approximation to the model of Pasciak and Gavis. Deep Sea Research Part I: Oceanographic Research Papers. 2008;55(10):1311–1317.
View Article
Google Scholar

[12] View Article

[13] Google Scholar

[ref5] 5. Aksnes D, Cao F. Inherent and apparent traits in microbial nutrient uptake. Marine Ecology Progress Series. 2011;440:41–51.
View Article
Google Scholar

[15] View Article

[16] Google Scholar

[ref6] 6. Fiksen Ø, Follows MJ, Aksnes DL. Trait-based models of nutrient uptake in microbes extend the Michaelis-Menten framework. Limnology and Oceanography. 2013;58(1):193–202.
View Article
Google Scholar

[18] View Article

[19] Google Scholar

[ref7] 7. Lindemann C,Fiksen Ø, Andersen KH, Aksnes DL. Scaling laws in phytoplankton nutrient uptake affinity. Frontiers in Marine Science. 2016;3.
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref8] 8. Zeng H, Yang A. Bridging substrate intake kinetics and bacterial growth phenotypes with flux balance analysis incorporating proteome allocation. Scientific Reports. 2020;10(1).
View Article
Google Scholar

[24] View Article

[25] Google Scholar

[ref9] 9. Pahlow M. Linking chlorophyll-nutrient dynamics to the Redfield N:C ratio with a model of optimal phytoplankton growth. Marine Ecology Progress Series. 2005;287:33–43.
View Article
Google Scholar

[27] View Article

[28] Google Scholar

[ref10] 10. Smith SL, Yamanaka Y. Optimization-based model of multinutrient uptake kinetics. Limnology and Oceanography. 2007;53(4):1545–1558.
View Article
Google Scholar

[30] View Article

[31] Google Scholar

[ref11] 11. Smith S, Yamanaka Y, Pahlow M, Oschlies A. Optimal uptake kinetics: physiological acclimation explains the pattern of nitrate uptake by phytoplankton in the ocean. Marine Ecology Progress Series. 2009;384:1–12.
View Article
Google Scholar

[33] View Article

[34] Google Scholar

[ref12] 12. Smith SL, Pahlow M, Merico A, Wirtz KW. Optimality-based modeling of planktonic organisms. Limnology and Oceanography. 2011;56(6):2080–2094.
View Article
Google Scholar

[36] View Article

[37] Google Scholar

[ref13] 13. Bonachela JA, Raghib M, Levin SA. Dynamic model of flexible phytoplankton nutrient uptake. Proceedings of the National Academy of Sciences of the United States of America. 2011;108(51):20633–20638. pmid:22143781
View Article
PubMed/NCBI
Google Scholar

[39] View Article

[40] PubMed/NCBI

[41] Google Scholar

[ref14] 14. Munk WH, Riley GA. Absorption of nutrients by aquatic plants. Journal of Marine Research. 1952;11:215–240.
View Article
Google Scholar

[43] View Article

[44] Google Scholar

[ref15] 15. Schmidt A, Kochanowski K, Vedelaar S, Ahrné E, Volkmer B, Callipo L, et al. The quantitative and condition-dependent Escherichia coli proteome. Nature Biotechnology. 2015;34(1):104–110. pmid:26641532
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref16] 16. Monk JM, Charusanti P, Aziz RK, Lerman JA, Premyodhin N, Orth JD, et al. Genome-scale metabolic reconstructions of multiple Escherichia coli strains highlight strain-specific adaptations to nutritional environments. Proceedings of the National Academy of Sciences. 2013;110(50):20338–20343.
View Article
Google Scholar

[50] View Article

[51] Google Scholar

[ref17] 17. Davidi D, Noor E, Liebermeister W, Bar-Even A, Flamholz A, Tummler K, et al. Global characterization of in vivo enzyme catalytic rates and their correspondence to in vitro k_cat measurements. Proceedings of the National Academy of Sciences. 2016;113(12):3401–3406.
View Article
Google Scholar

[53] View Article

[54] Google Scholar

[ref18] 18. Gavis J. Munk and Riley revisited: nutrient diffusion transport and rates of phytoplankton growth. Journal of Plankton Research. 1976;34(2):161–179.
View Article
Google Scholar

[56] View Article

[57] Google Scholar

[ref19] 19. Lazier JRN, Mann KH. Turbulence and the diffusive layers around small organisms. Deep-Sea Research. 1989;36(11):1721–1733.
View Article
Google Scholar

[59] View Article

[60] Google Scholar

[ref20] 20. Lenski RE, Rose MR, Simpson SC, Tadler SC. Long-term experimental evolution in Escherichia coli. I. Adaptation and divergence during 2,000 generations. American Naturalist. 1991;138(6):1315–1341.
View Article
Google Scholar

[62] View Article

[63] Google Scholar

[ref21] 21. Heirendt L, Arreckx S, Pfau T, Mendoza SN, Richelle A, Heinken A, et al. Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nature Protocols. 2019;14(3):639–702. pmid:30787451
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref22] 22. Ebrahim A, Lerman JA, Palsson BØ, Hyduke DR. COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Systems Biology. 2013;7:74. pmid:23927696
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref23] 23. Källberg M, Wang H, Wang S, Peng J, Wang Z, Lu H, et al. Template-based protein structure modeling using the RaptorX web server. Nature Protocols. 2012;7(8):1511–1522. pmid:22814390
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref24] 24. Hayduk W, Laudie H. Prediction of diffusion coefficients for nonelectrolytes in dilute aqueous solutions. AIChE Journal. 1974;20:611–615.
View Article
Google Scholar

[77] View Article

[78] Google Scholar

[ref25] 25. Arnold JH. Studies in diffusion, II Kinetic theory of diffusion in liquid systems. Journal of the American Chemical Society. 1930;52:3937–3948.
View Article
Google Scholar

[80] View Article

[81] Google Scholar

[ref26] 26. El-Dessouky HT, Ettouney HM. Fundamentals of Salt Water Desalination. 1st ed. Elsevier Science; 2002.

[ref27] 27. Steinsiek S, Bettenbrock K. Glucose transport in Escherichia coli mutant strains with defects in sugar transport systems. Journal of Bacteriology. 2012;194(21):5897–5908. pmid:22923596
View Article
PubMed/NCBI
Google Scholar

[84] View Article

[85] PubMed/NCBI

[86] Google Scholar

[ref28] 28. Graham AI, Hunt S, Stokes SL, Bramall N, Bunch J, Cox AG, et al. Severe zinc depletion of Escherichia coli. Journal of Biological Chemistry. 2009;284(27):18377–18389. pmid:19377097
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref29] 29. Lindblow-Kull C, Kull FJ, Shrift A. Single transporter for sulfate, selenate and selenite in Escherichia coli K12. Journal of Bacteriology. 1985;163(3):1267–1269.
View Article
Google Scholar

[92] View Article

[93] Google Scholar

[ref30] 30. Willsky GR, Bennett RL, Malamy MH. Inorganic phosphate transport in Escherichia coli: Involvement of two genes which play a role in alkaline phosphatase regulation. Journal of Bacteriology. 1973;113(2):529–539. pmid:4570598
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref31] 31. Gimenez R, Nunez MF, Badia J, Aguilar J, Baldoma L. The Gene yjcG, cotranscribed with the gene acs, encodes an acetate permease in Escherichia coli. Journal of Bacteriology. 2003;185(21):6448–6455. pmid:14563880
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref32] 32. Charbit A, Reizer J, Saier MH Jr. Function of the duplicated IIB domain and oligomeric structure of the fructose permease of Escherichia coli. Journal of Biological Chemistry. 1996;271(17):9997–10003. pmid:8626640
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref33] 33. Engel P, Kramer R, Unden G. Anaerobic fumarate transport in Escherichia coli by an fnr-dependent dicarboxylate uptake system which is different from the aerobic dicarboxylate uptake system. Journal of Bacteriology. 1992;174(17):5533–5539. pmid:1512189
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref34] 34. McDonald TP, Walmsley AR, Henderson PJF. Asparagine 394 in putative helix 11 of the galactose-H+ symport protein (GalP) from Escherichia coli is associated with the internal binding site for cytochalasin B and sugar. Journal of Biological Chemistry. 1997;272(24):15189–15199. pmid:9182541
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref35] 35. Hunter IS, Kornberg HL. Glucose transport of Escherichia coli growing in glucose-limited continuous culture. Biochemical Journal. 1979;178:97–101. pmid:373752
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref36] 36. Truniger V, Boos W. Glycerol uptake in Escherichia coli is sensitive to membrane lipid composition. Research in Microbiology. 1993;144:565–574. pmid:8310182
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref37] 37. Shamanna DK, Sanderson KE. Uptake and catabolism of D-Xylose in Salmonella typhimurium LT2. Journal of Bacteriology. 1979;139(1):64–70. pmid:222731
View Article
PubMed/NCBI
Google Scholar

[123] View Article

[124] PubMed/NCBI

[125] Google Scholar

[ref38] 38. Tilman D. The Resource-Ratio Hypothesis of plant succession. The American Naturalist. 1985;125(6):827–852.
View Article
Google Scholar

[127] View Article

[128] Google Scholar

[ref39] 39. Droop MR. Some thoughts on nutrient limitation in algae. Journal of Phycology. 1973;9:264–272.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref40] 40. Monod J. The growth of bacterial cultures. Annual Reviews Microbiology. 1949;3:371–394.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref41] 41. Sommer U. A comparison of the Droop and the Monod models of nutrient limited growth applied to natural populations of phytoplankton. Functional Ecology. 1991;5(4):535–544.
View Article
Google Scholar

[136] View Article

[137] Google Scholar

[ref42] 42. Blackman FF. Optima and limiting factors. Annals of Botany. 1905;19(74):281–295.
View Article
Google Scholar

[139] View Article

[140] Google Scholar

[ref43] 43. Powell EO. The growth rate of microorganisms as a function of substrate concentration. In: Powell PO, Evans CGT, Strange RE, Tempest DW, editors. Microbial physiology and continuous culture. H.M.S.O. London; 1967. p. 34–56.

[ref44] 44. Westerhoff HV, Lolkema LS, Otto R, Hellingwerf KJ. Thermodynamics of growth: Non-equilibrium thermodynamics of bacterial growth: The phenomenological and the mosaic approach. Bioehimica et Biophysica Acta. 1982;683:181–220.
View Article
Google Scholar

[143] View Article

[144] Google Scholar

[ref45] 45. Jannasch HW, Egli T. Microbial growth kinetics: a historical perspective. Antonie van Leeuwnhoek. 1993;63:213–224.
View Article
Google Scholar

[146] View Article

[147] Google Scholar

[ref46] 46. Condrey RE. The chemostat and Blackman kinetics. Biotechnology and Bioengineering. 1982;24:1705–1709. pmid:18546467
View Article
PubMed/NCBI
Google Scholar

[149] View Article

[150] PubMed/NCBI

[151] Google Scholar

[ref47] 47. Dabes JN, Finn RK, Wilke CR. Equations of substrate-limited growth: The case for Blackman kinetics. Biotechnology and Bioengineering. 1973;15:1159–1177. pmid:4203180
View Article
PubMed/NCBI
Google Scholar

[153] View Article

[154] PubMed/NCBI

[155] Google Scholar

[ref48] 48. Riley GA. Factors controlling phytoplankton populations on Georges Bank. Journal of Marine Research. 1946;6:54–73.
View Article
Google Scholar

[157] View Article

[158] Google Scholar

[ref49] 49. Heckmann D, Lloyd CJ, Mih N, Ha Y, Zielinski DC, Haiman ZB, et al. Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models. Nature Communications. 2018;9(1):1–10.
View Article
Google Scholar

[160] View Article

[161] Google Scholar

Figures

Abstract

Author summary

Introduction

Results

Conceptual model

Kinetic model

Acclimation

Cell size and surface area limitation

Model validation

Acclimation to limiting nutrients

Discussion

Methods

Determination of transporter kcat by quantitative proteomics and flux balance analysis

Determination of transporter capture area

Estimation of molecular diffusivities

Physiological data

Supporting information

S1 Fig. Comparison of the abundances of permease domains of known glucose transport systems in Escherichia coli K12 grown in glucose-replete batch cultures and glucose-limited chemostats, as determined by quantitative proteomics ([15]; Data plotted with permission from the authors).

References

Determination of transporter k_cat by quantitative proteomics and flux balance analysis