Timing the start of division in E. coli: a single-cell study

Whole-cell modeling of E. coli confirms that in vitro tRNA aminoacylation measurements are insufficient to support cell growth and predicts a positive feedback mechanism regulating arginine biosynthesis

In Escherichia coli, inconsistencies between in vitro tRNA aminoacylation measurements and in vivo protein synthesis demands were postulated almost 40 years ago, but have proven difficult to confirm. Whole-cell modeling can test whether a cell behaves in a physiologically correct manner when parameterized with in vitro measurements by providing a holistic representation of cellular processes in vivo. Here, a mechanistic model of tRNA aminoacylation, codon-based polypeptide elongation, and N-terminal methionine cleavage was incorporated into a developing whole-cell model of E. coli. Subsequent analysis confirmed the insufficiency of aminoacyl-tRNA synthetase kinetic measurements for cellular proteome maintenance, and estimated aminoacyl-tRNA synthetase kcats that were on average 7.6-fold higher. Simulating cell growth with perturbed kcats demonstrated the global impact of these in vitro measurements on cellular phenotypes. For example, an insufficient kcat for HisRS caused protein synthesis to be less robust to the natural variability in aminoacyl-tRNA synthetase expression in single cells. More surprisingly, insufficient ArgRS activity led to catastrophic impacts on arginine biosynthesis due to underexpressed N-acetylglutamate synthase, where translation depends on repeated CGG codons. Overall, the expanded E. coli model deepens understanding of how translation operates in an in vivo context.

Interpretation of organizational role of proteins on E. coli nucleoid via Hi-C integrated model

Several proteins in Escherichia coli work together to maintain the complex organization of its chromosome. However, the individual roles of these so-called nucleoid-associated proteins (NAPs) in chromosome architectures are not well characterized. Here, we quantitatively dissect the organizational roles of Heat Unstable (HU), a ubiquitous protein in E. coli and MatP, an NAP specifically binding to the Ter macrodomain of the chromosome. Toward this end, we employ a polymer physics-based computer model of wild-type chromosome and their HU- and MatP-devoid counterparts by incorporating their respective experimentally derived Hi-C contact matrix, cell dimensions, and replication status of the chromosome commensurate with corresponding growth conditions. Specifically, our model for the HU-devoid chromosome corroborates well with the microscopy observation of compaction of chromosome at short genomic range but diminished long-range interactions, justifying precedent hypothesis of segregation defect upon HU removal. Control simulations point out that the change in cell dimension and chromosome content in the process of HU removal holds the key to the observed differences in chromosome architecture between wild-type and HU-devoid cells. On the other hand, simulation of MatP-devoid chromosome led to locally enhanced contacts between Ter and its flanking macrodomains, consistent with previous recombination assay experiments and MatP's role in insulation of the Ter macrodomain from the rest of the chromosome. However, the simulation indicated no change in matS sites' localization. Rather, a set of designed control simulations showed that insulation of Ter is not caused by bridging of distant matS sites, also lending credence to a recent mobility experiment on various loci of the E. coli chromosome. Together, the investigations highlight the ability of an integrative model of the bacterial genome in elucidating the role of NAPs and in reconciling multiple experimental observations.

Intra-colony channel morphology in Escherichia coli biofilms is governed by nutrient availability and substrate stiffness

Nutrient-transporting channels have been recently discovered in mature Escherichia coli biofilms, however the relationship between intra-colony channel structure and the surrounding environmental conditions is poorly understood. Using a combination of fluorescence mesoscopy and a purpose-designed open-source quantitative image analysis pipeline, we show that growth substrate composition and nutrient availability have a profound effect on the morphology of intra-colony channels in mature E. coli biofilms. Under all nutrient conditions, intra-colony channel width was observed to increase non-linearly with radial distance from the centre of the biofilm. Notably, the channels were around 25% wider at the centre of carbon-limited biofilms compared to nitrogen-limited biofilms. Channel density also differed in colonies grown on rich and minimal media, with the former creating a network of tightly packed channels and the latter leading to well-separated, wider channels with defined edges. Our approach paves the way for measurement of internal patterns in a wide range of biofilms, offering the potential for new insights into infection and pathogenicity.

Label-Free Optical Detection of Pathogenic Bacteria and Fungi at Extremely Low Cell Densities for Rapid Antibiotic Susceptibility Testing

Antibiotic resistance is a rapidly expanding public health problem across the globe leading to prolonged hospital admissions, increased morbidity and mortality, and associated high healthcare costs. Effective treatment of bacterial infections requires timely and correct antibiotic administration to the patients which relies on rapid phenotyping of disease-causing bacteria. Currently, antibiotic susceptibility tests can take several days and as a result, indiscriminate antibiotic use has exacerbated the evolution and spread of antibiotic resistance in clinical and community settings. In order to address this problem, we have developed a novel optical apparatus that we called RUSD (Rapid Ultra-Sensitive Detection). RUSD is built around a hollow silica fiber and utilizes bacterial cells as spatial light modulators. This generates a highly sensitive modulation transfer function due to the narrow reflectivity angle in the fiber-media interface. We leveraged the RUSD technology to allow for robust bacterial and fungal detection. RUSD can now detect pathogenic cell densities in a large dynamic window (OD600 from ∼10−7 to 10−1). Finally, we can generate dose response curves for various pathogens and antimicrobial compounds within one to three hours by using RUSD. Our antibiotic- susceptibility testing (AST) assay that we call iFAST (in-Fiber-Antibiotic-Susceptibility-Testing) is fast, highly sensitive, and does not change the existing workflow in clinical settings as it is compatible with FDA-approved AST. Thus, RUSD platform is a viable tool that will expedite decision-making process in the treatment of infectious diseases and positively impact the antibiotic resistance problem in the long term by minimizing the use of ineffective antibiotics.

Hi-C embedded polymer model of Escherichia coli reveals the origin of heterogeneous subdiffusion in chromosomal loci

Underneath its apparently simple architecture, the circular chromosome of Escherichia coli is known for displaying complex dynamics in its cytoplasm, with past investigations hinting at inherently diverse mobilities of chromosomal loci across the genome. To decipher its origin, we simulate the dynamics of genome-wide spectrum of E. coli chromosomal loci, via integrating its experimentally derived Hi-C interaction matrix within a polymer-based model. Our analysis demonstrates that, while the dynamics of the chromosome is subdiffusive in a viscoelastic media, the diffusion constants are strongly dependent of chromosomal loci coordinates and diffusive exponents (α) are widely heterogenous with α ≈ 0.36-0.60. The loci-dependent heterogeneous dynamics and mean first-passage times of interloci encounter were found to be modulated via genetically distant interloci communications and is robust even in the presence of active, ATP-dependent noises. Control investigations reveal that the absence of Hi-C-derived interactions in the model would have abolished the traits of heterogeneous loci diffusion, underscoring the key role of loci-specific genetically distant interaction in modulating the underlying heterogeneity of the loci diffusion.

Stochastic nucleoid segregation dynamics as a source of the phenotypic variability in E. coli

Segregation of the replicating chromosome from a single to two nucleoid bodies is one of the major processes in growing bacterial cells. The segregation dynamics is tuned by intricate interactions with other cellular processes such as growth and division, ensuring flexibility in a changing environment. We hypothesize that the internal stochasticity of the segregation process may be the source of cell-to-cell phenotypic variability, in addition to the well-established gene expression noise and uneven partitioning of low copy number components. We compare dividing cell lineages with filamentous cells, where the lack of the diffusion barriers is expected to reduce the impact of other factors on the variability of nucleoid segregation dynamics. The nucleoid segregation was monitored using time-lapse microscopy in live E. coli cells grown in linear grooves. The main characteristics of the segregation process, namely, the synchrony of partitioning, rates of separation, and final positions, as well as the variability of these characteristics, were determined for dividing and filamentous lineages growing under the same conditions. Indeed, the gene expression noise was considerably homogenized along filaments as determined from the distribution of CFP and YFP stochastically expressed from the chromosome. We find that 1) the synchrony of nucleoid partitioning is progressively decreasing during consecutive cell cycles, but to a significantly lesser degree in filamentous than in dividing cells; 2) the mean partitioning rate of nucleoids is essentially the same in dividing and filamentous cells, displaying a substantial variability in both; and 3) nucleoids segregate to the same distances in dividing and filamentous cells. Variability in distances is increasing during successive cell cycles, but to a much lesser extent in filamentous cells. Our findings indicate that the variability of the chromosome segregation dynamics is reduced upon removal of boundaries between nucleoids, whereas the remaining variability is essentially inherent to the nucleoid itself.

Hyperosmotic Shock Transiently Accelerates Constriction Rate in Escherichia coli

Bacterial cells in their natural environments encounter rapid and large changes in external osmolality. For instance, enteric bacteria such as Escherichia coli experience a rapid decrease when they exit from host intestines. Changes in osmolality alter the mechanical load on the cell envelope, and previous studies have shown that large osmotic shocks can slow down bacterial growth and impact cytoplasmic diffusion. However, it remains unclear how cells maintain envelope integrity and regulate envelope synthesis in response to osmotic shocks. In this study, we developed an agarose pad-based protocol to assay envelope stiffness by measuring population-averaged cell length before and after a hyperosmotic shock. Pad-based measurements exhibited an apparently larger length change compared with single-cell dynamics in a microfluidic device, which we found was quantitatively explained by a transient increase in division rate after the shock. Inhibiting cell division led to consistent measurements between agarose pad-based and microfluidic measurements. Directly after hyperosmotic shock, FtsZ concentration and Z-ring intensity increased, and the rate of septum constriction increased. These findings establish an agarose pad-based protocol for quantifying cell envelope stiffness, and demonstrate that mechanical perturbations can have profound effects on bacterial physiology.

A Hi-C data-integrated model elucidates E. coli chromosome's multiscale organization at various replication stages

The chromosome of Escherichia coli is riddled with multi-faceted complexity. The emergence of chromosome conformation capture techniques are providing newer ways to explore chromosome organization. Here we combine a beads-on-a-spring polymer-based framework with recently reported Hi-C data for E. coli chromosome, in rich growth condition, to develop a comprehensive model of its chromosome at 5 kb resolution. The investigation focuses on a range of diverse chromosome architectures of E. coli at various replication states corresponding to a collection of cells, individually present in different stages of cell cycle. The Hi-C data-integrated model captures the self-organization of E. coli chromosome into multiple macrodomains within a ring-like architecture. The model demonstrates that the position of oriC is dependent on architecture and replication state of chromosomes. The distance profiles extracted from the model reconcile fluorescence microscopy and DNA-recombination assay experiments. Investigations into writhe of the chromosome model reveal that it adopts helix-like conformation with no net chirality, earlier hypothesized in experiments. A genome-wide radius of gyration map captures multiple chromosomal interaction domains and identifies the precise locations of rrn operons in the chromosome. We show that a model devoid of Hi-C encoded information would fail to recapitulate most genomic features unique to E. coli.

Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation

The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle-and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Successive Paradigm Shifts in the Bacterial Cell Cycle and Related Subjects

A paradigm shift in one field can trigger paradigm shifts in other fields. This is illustrated by the paradigm shifts that have occurred in bacterial physiology following the discoveries that bacteria are not unstructured, that the bacterial cell cycle is not controlled by the dynamics of peptidoglycan, and that the growth rates of bacteria in the same steady-state population are not at all the same. These paradigm shifts are having an effect on longstanding hypotheses about the regulation of the bacterial cell cycle, which appear increasingly to be inadequate. I argue that, just as one earthquake can trigger others, an imminent paradigm shift in the regulation of the bacterial cell cycle will have repercussions or "paradigm quakes" on hypotheses about the origins of life and about the regulation of the eukaryotic cell cycle.

Constriction Rate Modulation Can Drive Cell Size Control and Homeostasis in C. crescentus

Rod-shaped bacteria typically grow first via sporadic and dispersed elongation along their lateral walls and then via a combination of zonal elongation and constriction at the division site to form the poles of daughter cells. Although constriction comprises up to half of the cell cycle, its impact on cell size control and homeostasis has rarely been considered. To reveal the roles of cell elongation and constriction in bacterial size regulation during cell division, we captured the shape dynamics of Caulobacter crescentus with time-lapse structured illumination microscopy and used molecular markers as cell-cycle landmarks. We perturbed the constriction rate using a hyperconstriction mutant or fosfomycin ([(2R,3S)-3-methyloxiran-2-yl]phosphonic acid) inhibition. We report that the constriction rate contributes to both size control and homeostasis, by determining elongation during constriction and by compensating for variation in pre-constriction elongation on a single-cell basis.

•

Perturbing constriction rate changes cell length

•

Faster constriction rate results in blunter cell poles

•

Early constriction rate modulation balances elongation before and during constriction

•

We propose that constriction rate is set by the accumulation of precursors during elongation

Surface Area to Volume Ratio: A Natural Variable for Bacterial Morphogenesis

An immediately observable feature of bacteria is that cell size and shape are remarkably constant and characteristic for a given species in a particular condition, but vary quantitatively with physiological parameters such as growth rate, indicating both genetic and environmental regulation. However, despite decades of research, the molecular mechanisms underlying bacterial morphogenesis have remained incompletely characterized. We recently demonstrated that a wide range of bacterial species exhibit a robust surface area to volume ratio (SA/V) homeostasis. Because cell size, shape, and SA/V are mathematically interconnected, if SA/V is indeed the natural variable that cells actively monitor, this finding has critical implications for our understanding of bacterial morphogenesis, placing fundamental constraints on the sizes and shapes that cells can adopt. In this Opinion article we discuss the broad implications that this novel perspective has for the field of bacterial growth and morphogenesis.

Division-Based, Growth Rate Diversity in Bacteria

To investigate the nature and origins of growth rate diversity in bacteria, we grew Escherichia coli and Bacillus subtilis in liquid minimal media and, after different periods of ¹⁵N-labeling, analyzed and imaged isotope distributions in individual cells with Secondary Ion Mass Spectrometry. We find a striking inter- and intra-cellular diversity, even in steady state growth. This is consistent with the strand-dependent, hyperstructure-based hypothesis that a major function of the cell cycle is to generate coherent, growth rate diversity via the semi-conservative pattern of inheritance of strands of DNA and associated macromolecular assemblies. We also propose quantitative, general, measures of growth rate diversity for studies of cell physiology that include antibiotic resistance.

Different Amounts of DNA in Newborn Cells of Escherichia coli Preclude a Role for the Chromosome in Size Control According to the "Adder" Model

According to the recently-revived adder model for cell size control, newborn cells of Escherichia coli will grow and divide after having added a constant size or length, ΔL, irrespective of their size at birth. Assuming exponential elongation, this implies that large newborns will divide earlier than small ones. The molecular basis for the constant size increment is still unknown. As DNA replication and cell growth are coordinated, the constant ΔL could be based on duplication of an equal amount of DNA, ΔG, present in newborn cells. To test this idea, we measured amounts of DNA and lengths of nucleoids in DAPI-stained cells growing in batch culture at slow and fast rates. Deeply-constricted cells were divided in two subpopulations of longer and shorter lengths than average; these were considered to represent large and small prospective daughter cells, respectively. While at slow growth, large and small prospective daughter cells contained similar amounts of DNA, fast growing cells with multiforked replicating chromosomes, showed a significantly higher amount of DNA (20%) in the larger cells. This observation precludes the hypothesis that ΔL is based on the synthesis of a constant ΔG. Growth curves were constructed for siblings generated by asymmetric division and growing according to the adder model. Under the assumption that all cells at the same growth rate exhibit the same time between initiation of DNA replication and cell division (i.e., constant C+D-period), the constructions predict that initiation occurs at different sizes (Li) and that, at fast growth, large newborn cells transiently contain more DNA than small newborns, in accordance with the observations. Because the state of segregation, measured as the distance between separated nucleoids, was found to be more advanced in larger deeply-constricted cells, we propose that in larger newborns nucleoid separation occurs faster and at a shorter length, allowing them to divide earlier. We propose a composite model in which both differential initiation and segregation leads to an adder-like behavior of large and small newborn cells.

Does the eclipse limit bacterial nucleoid complexity and cell width?

Cell size of bacteria M is related to 3 temporal parameters: chromosome replication time C, period from replication-termination to subsequent division D, and doubling time τ. Steady-state, bacillary cells grow exponentially by extending length L only, but their constant width W is larger at shorter τ's or longer C's, in proportion to the number of chromosome replication positions n (= C/τ), at least in Escherichia coli and Salmonella typhimurium. Extending C by thymine limitation of fast-growing thyA mutants result in continuous increase of M, associated with rising W, up to a limit before branching. A set of such puzzling observations is qualitatively consistent with the view that the actual cell mass (or volume) at the time of replication-initiation Mi (or Vi), usually relatively constant in growth at varying τ's, rises with time under thymine limitation of fast-growing, thymine-requiring E. coli strains. The hypothesis will be tested that presumes existence of a minimal distance l_min between successive moving replisomes, translated into the time needed for a replisome to reach l_min before a new replication-initiation at oriC is allowed, termed Eclipse E. Preliminary analysis of currently available data is inconsistent with a constant E under all conditions, hence other explanations and ways to test them are proposed in an attempt to elucidate these and other results. The complex hypothesis takes into account much of what is currently known about Bacterial Physiology: the relationships between cell dimensions, growth and cycle parameters, particularly nucleoid structure, replication and position, and the mode of peptidoglycan biosynthesis. Further experiments are mentioned that are necessary to test the discussed ideas and hypotheses.

Vibrio cholerae Combines Individual and Collective Sensing to Trigger Biofilm Dispersal

Bacteria can generate benefits for themselves and their kin by living in multicellular, matrix-enclosed communities, termed biofilms, which are fundamental to microbial ecology and the impact bacteria have on the environment, infections, and industry [1, 2, 3, 4, 5, 6]. The advantages of the biofilm mode of life include increased stress resistance and access to concentrated nutrient sources [3, 7, 8]. However, there are also costs associated with biofilm growth, including the metabolic burden of biofilm matrix production, increased resource competition, and limited mobility inside the community [9, 10, 11]. The decision-making strategies used by bacteria to weigh the costs between remaining in a biofilm or actively dispersing are largely unclear, even though the dispersal transition is a central aspect of the biofilm life cycle and critical for infection transmission [12, 13, 14]. Using a combination of genetic and novel single-cell imaging approaches, we show that Vibrio cholerae integrates dual sensory inputs to control the dispersal response: cells use the general stress response, which can be induced via starvation, and they also integrate information about the local cell density and molecular transport conditions in the environment via the quorum sensing apparatus. By combining information from individual (stress response) and collective (quorum sensing) avenues of sensory input, biofilm-dwelling bacteria can make robust decisions to disperse from large biofilms under distress, while preventing premature dispersal when biofilm populations are small. These insights into triggers and regulators of biofilm dispersal are a key step toward actively inducing biofilm dispersal for technological and medical applications, and for environmental control of biofilms.

•

Cells in V. cholerae biofilms decide to disperse by combining two sensory mechanisms

•

Quorum sensing and RpoS provide information on different environmental parameters

•

Integration of both sensory inputs yields robust and optimal dispersal decisions

HU content and dynamics in Escherichia coli during the cell cycle and at different growth rates

DNA-binding proteins play an important role in maintaining bacterial chromosome structure and functions. Heat-unstable (HU) histone-like protein is one of the most abundant of these proteins and participates in all major chromosome-related activities. Owing to its low sequence specificity, HU fusions with fluorescent proteins were used for general staining of the nucleoid, aiming to reveal its morphology and dynamics. We have exploited a single chromosomal copy of hupA-egfp fusion under the native promoter and used quantitative microscopy imaging to investigate the amount and dynamics of HUα in Escherichia coli cells. We found that in steady-state growing populations the cellular HUα content is proportional to the cell size, whereas its concentration is size independent. Single-cell live microscopy imaging confirmed that the amount of HUα exponentially increases during the cell cycle, but its concentration is maintained constant. This supports the existence of an auto-regulatory mechanism underlying the HUα cellular level, in addition to reflecting the gene copy number. Both the HUα amount and concentration strongly increase with the cell growth rate in different culture media. Unexpectedly, the HU/DNA stoichiometry also remarkably increases with the growth rate. This last finding may be attributed to a higher requirement for maintaining the chromosome structure in nucleoids with higher complexity.

Z-ring Structure and Constriction Dynamics in E. coli

The Z-ring plays a central role in bacterial division. It consists of FtsZ filaments, but the way these reorganize in the ring-like structure during septation remains largely unknown. Here, we measure the effective constriction dynamics of the ring. Using an oscillating optical trap, we can switch individual rod-shaped E. coli cells between horizontal and vertical orientations. In the vertical orientation, the fluorescent Z-ring image appears as a symmetric circular structure that renders itself to quantitative analysis. In the horizontal orientation, we use phase-contrast imaging to determine the extent of the cell constriction and obtain the effective time of division. We find evidence that the Z-ring constricts at a faster rate than the cell envelope such that its radial width (inwards from the cytoplasmic membrane) grows during septation. In this respect, our results differ from those recently obtained using photoactivated localization microscopy (PALM) where the radial width of the Z-ring was found to be approximately constant as the ring constricts. A possible reason for the different behavior of the constricting Z-rings could be the significant difference in the corresponding cell growth rates.

Threshold effect of growth rate on population variability of Escherichia coli cell lengths

A long-standing question in biology is the effect of growth on cell size. Here, we estimate the effect of Escherichia coli growth rate (r) on population cell size distributions by estimating the coefficient of variation of cell lengths (CV_L) from image analysis of fixed cells in DIC microscopy. We find that the CV_L is constant at growth rates less than one division per hour, whereas above this threshold, CV_L increases with an increase in the growth rate. We hypothesize that stochastic inhibition of cell division owing to replication stalling by a RecA-dependent mechanism, combined with the growth rate threshold of multi-fork replication (according to Cooper and Helmstetter), could form the basis of such a threshold effect. We proceed to test our hypothesis by increasing the frequency of stochastic stalling of replication forks with hydroxyurea (HU) treatment and find that cell length variability increases only when the growth rate exceeds this threshold. The population effect is also reproduced in single-cell studies using agar-pad cultures and 'mother machine'-based experiments to achieve synchrony. To test the role of RecA, critical for the repair of stalled replication forks, we examine the CV_L of E. coli ΔrecA cells. We find cell length variability in the mutant to be greater than wild-type, a phenotype that is rescued by plasmid-based RecA expression. Additionally, we find that RecA-GFP protein recruitment to nucleoids is more frequent at growth rates exceeding the growth rate threshold and is further enhanced on HU treatment. Thus, we find growth rates greater than a threshold result in increased E. coli cell lengths in the population, and this effect is, at least in part, mediated by RecA recruitment to the nucleoid and stochastic inhibition of division.

Threshold effect of growth rate on population variability of Escherichia coli cell lengths

A long-standing question in biology is the effect of growth on cell size. Here, we estimate the effect of Escherichia coli growth rate (r) on population cell size distributions by estimating the coefficient of variation of cell lengths (CV_L) from image analysis of fixed cells in DIC microscopy. We find that the CV_L is constant at growth rates less than one division per hour, whereas above this threshold, CV_L increases with an increase in the growth rate. We hypothesize that stochastic inhibition of cell division owing to replication stalling by a RecA-dependent mechanism, combined with the growth rate threshold of multi-fork replication (according to Cooper and Helmstetter), could form the basis of such a threshold effect. We proceed to test our hypothesis by increasing the frequency of stochastic stalling of replication forks with hydroxyurea (HU) treatment and find that cell length variability increases only when the growth rate exceeds this threshold. The population effect is also reproduced in single-cell studies using agar-pad cultures and 'mother machine'-based experiments to achieve synchrony. To test the role of RecA, critical for the repair of stalled replication forks, we examine the CV_L of E. coli ΔrecA cells. We find cell length variability in the mutant to be greater than wild-type, a phenotype that is rescued by plasmid-based RecA expression. Additionally, we find that RecA-GFP protein recruitment to nucleoids is more frequent at growth rates exceeding the growth rate threshold and is further enhanced on HU treatment. Thus, we find growth rates greater than a threshold result in increased E. coli cell lengths in the population, and this effect is, at least in part, mediated by RecA recruitment to the nucleoid and stochastic inhibition of division.

Effects of cell-cycle-dependent expression on random fluctuations in protein levels

Expression of many genes varies as a cell transitions through different cell-cycle stages. How coupling between stochastic expression and cell cycle impacts cell-to-cell variability (noise) in the level of protein is not well understood. We analyse a model where a stable protein is synthesized in random bursts, and the frequency with which bursts occur varies within the cell cycle. Formulae quantifying the extent of fluctuations in the protein copy number are derived and decomposed into components arising from the cell cycle and stochastic processes. The latter stochastic component represents contributions from bursty expression and errors incurred during partitioning of molecules between daughter cells. These formulae reveal an interesting trade-off: cell-cycle dependencies that amplify the noise contribution from bursty expression also attenuate the contribution from partitioning errors. We investigate the existence of optimum strategies for coupling expression to the cell cycle that minimize the stochastic component. Intriguingly, results show that a zero production rate throughout the cell cycle, with expression only occurring just before cell division, minimizes noise from bursty expression for a fixed mean protein level. By contrast, the optimal strategy in the case of partitioning errors is to make the protein just after cell division. We provide examples of regulatory proteins that are expressed only towards the end of the cell cycle, and argue that such strategies enhance robustness of cell-cycle decisions to the intrinsic stochasticity of gene expression.

Intercellular Variability in Protein Levels from Stochastic Expression and Noisy Cell Cycle Processes

Inside individual cells, expression of genes is inherently stochastic and manifests as cell-to-cell variability or noise in protein copy numbers. Since proteins half-lives can be comparable to the cell-cycle length, randomness in cell-division times generates additional intercellular variability in protein levels. Moreover, as many mRNA/protein species are expressed at low-copy numbers, errors incurred in partitioning of molecules between two daughter cells are significant. We derive analytical formulas for the total noise in protein levels when the cell-cycle duration follows a general class of probability distributions. Using a novel hybrid approach the total noise is decomposed into components arising from i) stochastic expression; ii) partitioning errors at the time of cell division and iii) random cell-division events. These formulas reveal that random cell-division times not only generate additional extrinsic noise, but also critically affect the mean protein copy numbers and intrinsic noise components. Counter intuitively, in some parameter regimes, noise in protein levels can decrease as cell-division times become more stochastic. Computations are extended to consider genome duplication, where transcription rate is increased at a random point in the cell cycle. We systematically investigate how the timing of genome duplication influences different protein noise components. Intriguingly, results show that noise contribution from stochastic expression is minimized at an optimal genome-duplication time. Our theoretical results motivate new experimental methods for decomposing protein noise levels from synchronized and asynchronized single-cell expression data. Characterizing the contributions of individual noise mechanisms will lead to precise estimates of gene expression parameters and techniques for altering stochasticity to change phenotype of individual cells.

Inside individual cells, gene products often occur at low molecular counts and are subject to considerable stochastic fluctuations (noise) in copy numbers over time. An important consequence of noisy expression is that the level of a protein can vary considerably even among genetically identical cells exposed to the same environment. Such non-genetic phenotypic heterogeneity is physiologically relevant and critically influences diverse cellular processes. In addition to noise sources inherent in gene product synthesis, recent experimental studies have uncovered additional noise mechanisms that critically effect expression. For example, the time within the cell cycle when a gene duplicates, and the time taken to complete cell cycle are governed by random processes. The key contribution of this work is development of novel mathematical results quantifying how cell cycle-related noise sources combine with stochastic expression to drive intercellular variability in protein molecular counts. Derived formulas lead to many counterintuitive results, such as increasing randomness in the timing of cell division can lower noise in the level of a protein. Finally, these results inform experimental strategies to systematically dissect the contributions of different noise sources in the expression of a gene of interest.

Phospholipids and glycolipids mediate proton containment and circulation along the surface of energy-transducing membranes

Proton bioenergetics provides the energy for growth and survival of most organisms in the biosphere ranging from unicellular marine phytoplankton to humans. Chloroplasts harvest light and generate a proton electrochemical gradient (proton motive force) that drives the production of ATP needed for carbon dioxide fixation and plant growth. Mitochondria, bacteria and archaea generate proton motive force to energize growth and other physiologies. Energy transducing membranes are at the heart of proton bioenergetics and are responsible for catalyzing the conversion of energy held in high-energy electrons→electron transport chain→proton motive force→ATP. Whereas the electron transport chain is understood in great detail there are major gaps in understanding mechanisms of proton transfer or circulation during proton bioenergetics. This paper is built on the proposition that phospho- and glyco-glycerolipids form proton transport circuitry at the membrane's surface. By this proposition, an emergent membrane property, termed the hyducton, confines active/unbound protons or hydronium ions to a region of low volume close to the membrane surface. In turn, a von Grotthuß mechanism rapidly moves proton substrate in accordance with nano-electrochemical poles on the membrane surface created by powerful proton pumps such as ATP synthase.

Noise and Epigenetic Inheritance of Single-Cell Division Times Influence Population Fitness

The fitness effect of biological noise remains unclear. For example, even within clonal microbial populations, individual cells grow at different speeds. Although it is known that the individuals' mean growth speed can affect population-level fitness, it is unclear how or whether growth speed heterogeneity itself is subject to natural selection. Here, we show that noisy single-cell division times can significantly affect population-level growth rate. Using time-lapse microscopy to measure the division times of thousands of individual S. cerevisiae cells across different genetic and environmental backgrounds, we find that the length of individual cells' division times can vary substantially between clonal individuals and that sublineages often show epigenetic inheritance of division times. By combining these experimental measurements with mathematical modeling, we find that, for a given mean division time, increasing heterogeneity and epigenetic inheritance of division times increases the population growth rate. Furthermore, we demonstrate that the heterogeneity and epigenetic inheritance of single-cell division times can be linked with variation in the expression of catabolic genes. Taken together, our results reveal how a change in noisy single-cell behaviors can directly influence fitness through dynamics that operate independently of effects caused by changes to the mean. These results not only allow a better understanding of microbial fitness but also help to more accurately predict fitness in other clonal populations, such as tumors.

Defining the rate-limiting processes of bacterial cytokinesis

Bacterial cytokinesis is accomplished by the essential 'divisome' machinery. The most widely conserved divisome component, FtsZ, is a tubulin homolog that polymerizes into the 'FtsZ-ring' ('Z-ring'). Previous in vitro studies suggest that Z-ring contraction serves as a major constrictive force generator to limit the progression of cytokinesis. Here, we applied quantitative superresolution imaging to examine whether and how Z-ring contraction limits the rate of septum closure during cytokinesis in Escherichia coli cells. Surprisingly, septum closure rate was robust to substantial changes in all Z-ring properties proposed to be coupled to force generation: FtsZ's GTPase activity, Z-ring density, and the timing of Z-ring assembly and disassembly. Instead, the rate was limited by the activity of an essential cell wall synthesis enzyme and further modulated by a physical divisome-chromosome coupling. These results challenge a Z-ring-centric view of bacterial cytokinesis and identify cell wall synthesis and chromosome segregation as limiting processes of cytokinesis.

Impact of the cell division cycle on gene circuits

In growing cells, protein synthesis and cell growth are typically not synchronous, and, thus, protein concentrations vary over the cell division cycle. We have developed a theoretical description of genetic regulatory systems in bacteria that explicitly considers the cell division cycle to investigate its impact on gene expression. We calculate the cell-to-cell variations arising from cells being at different stages in the division cycle for unregulated genes and for basic regulatory mechanisms. These variations contribute to the extrinsic noise observed in single-cell experiments, and are most significant for proteins with short lifetimes. Negative autoregulation buffers against variation of protein concentration over the division cycle, but the effect is found to be relatively weak. Stronger buffering is achieved by an increased protein lifetime. Positive autoregulation can strongly amplify such variation if the parameters are set to values that lead to resonance-like behaviour. For cooperative positive autoregulation, the concentration variation over the division cycle diminishes the parameter region of bistability and modulates the switching times between the two stable states. The same effects are seen for a two-gene mutual-repression toggle switch. By contrast, an oscillatory circuit, the repressilator, is only weakly affected by the division cycle.

Contributions of cell growth and biochemical reactions to nongenetic variability of cells

Cell-to-cell variability in the molecular composition of isogenic, steady-state growing cells arises spontaneously from the inherent stochasticity of intracellular biochemical reactions and cell growth. Here, we present a general decomposition of the total variance in the copy number per cell of a particular molecule. It quantifies the individual contributions made by processes associated with cell growth, biochemical reactions, and their control. We decompose the growth contribution further into variance contributions of random partitioning of molecules at cell division, mother-cell heterogeneity, and variation in cell-cycle progression. The contribution made by biochemical reactions is expressed in variance generated by molecule synthesis, degradation, and their regulation. We use this theory to study the influence of different growth and reaction-related processes, such as DNA replication, variable molecule-partitioning probability, and synthesis bursts, on stochastic cell-to-cell variability. Using simulations, we characterize the impact of noise in the generation-time on cell-to-cell variability. This article offers a widely-applicable theory on the influence of biochemical reactions and cellular growth on the phenotypic variability of growing, isogenic cells. The theory aids the design and interpretation of experiments involving single-molecule counting or real-time imaging of fluorescent reporter constructs.

Metabolism at evolutionary optimal States

Metabolism is generally required for cellular maintenance and for the generation of offspring under conditions that support growth. The rates, yields (efficiencies), adaptation time and robustness of metabolism are therefore key determinants of cellular fitness. For biotechnological applications and our understanding of the evolution of metabolism, it is necessary to figure out how the functional system properties of metabolism can be optimized, via adjustments of the kinetics and expression of enzymes, and by rewiring metabolism. The trade-offs that can occur during such optimizations then indicate fundamental limits to evolutionary innovations and bioengineering. In this paper, we review several theoretical and experimental findings about mechanisms for metabolic optimization.

A multi-layered protein network stabilizes the Escherichia coli FtsZ-ring and modulates constriction dynamics

The prokaryotic tubulin homolog, FtsZ, forms a ring-like structure (FtsZ-ring) at midcell. The FtsZ-ring establishes the division plane and enables the assembly of the macromolecular division machinery (divisome). Although many molecular components of the divisome have been identified and their interactions extensively characterized, the spatial organization of these proteins within the divisome is unclear. Consequently, the physical mechanisms that drive divisome assembly, maintenance, and constriction remain elusive. Here we applied single-molecule based superresolution imaging, combined with genetic and biophysical investigations, to reveal the spatial organization of cellular structures formed by four important divisome proteins in E. coli: FtsZ, ZapA, ZapB and MatP. We show that these interacting proteins are arranged into a multi-layered protein network extending from the cell membrane to the chromosome, each with unique structural and dynamic properties. Further, we find that this protein network stabilizes the FtsZ-ring, and unexpectedly, slows down cell constriction, suggesting a new, unrecognized role for this network in bacterial cell division. Our results provide new insight into the structure and function of the divisome, and highlight the importance of coordinated cell constriction and chromosome segregation.

Bacterial cell division is a highly regulated process that must be coordinated with other cellular processes (i.e. DNA replication and chromosome segregation) to promote faithful reproduction. In Escherichia coli, this regulation is most often mediated through the polymerization of the prokaryotic tubulin homolog, FtsZ, which forms a ring-like structure (FtsZ-ring) at midcell. The establishment of the FtsZ-ring marks the site of division and enables the assembly of the macromolecular division machinery (divisome). Here we applied single-molecule based superresolution imaging to reveal the three-dimensional structure of FtsZ in the context of its regulatory proteins: ZapA, ZapB and MatP. We found that these four proteins exist in a multi-layered network that extends from the cell membrane to the chromosome. This layered organization not only helps to stabilize the FtsZ-ring, but also serves to coordinate division with DNA status by influencing constriction rate. Our results not only provide a comprehensive view of the divisome, but also allow new insight to be garnered regarding the structure and function of the divisome.

Nanoscale imaging of the growth and division of bacterial cells on planar substrates with the atomic force microscope

With the use of the atomic force microscope (AFM), the Nanomicrobiology field has advanced drastically. Due to the complexity of imaging living bacterial processes in their natural growing environments, improvements have come to a standstill. Here we show the in situ nanoscale imaging of the growth and division of single bacterial cells on planar substrates with the atomic force microscope. To achieve this, we minimized the lateral shear forces responsible for the detachment of weakly adsorbed bacteria on planar substrates with the use of the so called dynamic jumping mode with very soft cantilever probes. With this approach, gentle imaging conditions can be maintained for long periods of time, enabling the continuous imaging of the bacterial cell growth and division, even on planar substrates. Present results offer the possibility to observe living processes of untrapped bacteria weakly attached to planar substrates.

How and why cells grow as rods

The rod is a ubiquitous shape adopted by walled cells from diverse organisms ranging from bacteria to fungi to plants. Although rod-like shapes are found in cells of vastly different sizes and are constructed by diverse mechanisms, the geometric similarities among these shapes across kingdoms suggest that there are common evolutionary advantages, which may result from simple physical principles in combination with chemical and physiological constraints. Here, we review mechanisms of constructing rod-shaped cells and the bases of different biophysical models of morphogenesis, comparing and contrasting model organisms in different kingdoms. We then speculate on possible advantages of the rod shape, and suggest strategies for elucidating the relative importance of each of these advantages.

Structure-related antibacterial activity of a titanium nanostructured surface fabricated by glancing angle sputter deposition

The aim of this study was to reproduce the physico-mechanical antibacterial effect of the nanocolumnar cicada wing surface for metallic biomaterials by fabrication of titanium (Ti) nanocolumnar surfaces using glancing angle sputter deposition (GLAD). Nanocolumnar Ti thin films were fabricated by GLAD on silicon substrates. S. aureus as well as E. coli were incubated with nanostructured or reference dense Ti thin film test samples for one or three hours at 37 °C. Bacterial adherence, morphology, and viability were analyzed by fluorescence staining and scanning electron microscopy and compared to human mesenchymal stem cells (hMSCs).Bacterial adherence was not significantly different after short (1 h) incubation on the dense or the nanostructured Ti surface. In contrast to S. aureus the viability of E. coli was significantly decreased after 3 h on the nanostructured film compared to the dense film and was accompanied by an irregular morphology and a cell wall deformation. Cell adherence, spreading and viability of hMSCs were not altered on the nanostructured surface. The results show that the selective antibacterial effect of the cicada wing could be transferred to a nanostructured metallic biomaterial by mimicking the natural nanocolumnar topography.

Transimulation - protein biosynthesis web service

Although translation is the key step during gene expression, it remains poorly characterized at the level of individual genes. For this reason, we developed Transimulation – a web service measuring translational activity of genes in three model organisms: Escherichia coli, Saccharomyces cerevisiae and Homo sapiens. The calculations are based on our previous computational model of translation and experimental data sets. Transimulation quantifies mean translation initiation and elongation time (expressed in SI units), and the number of proteins produced per transcript. It also approximates the number of ribosomes that typically occupy a transcript during translation, and simulates their propagation. The simulation of ribosomes’ movement is interactive and allows modifying the coding sequence on the fly. It also enables uploading any coding sequence and simulating its translation in one of three model organisms. In such a case, ribosomes propagate according to mean codon elongation times of the host organism, which may prove useful for heterologous expression. Transimulation was used to examine evolutionary conservation of translational parameters of orthologous genes. Transimulation may be accessed at http://nexus.ibb.waw.pl/Transimulation (requires Java version 1.7 or higher). Its manual and source code, distributed under the GPL-2.0 license, is freely available at the website.

Three-dimensional structure of the Z-ring as a random network of FtsZ filaments

The spatial organization of the Z-ring, the central element of the bacterial division machinery, is not yet fully understood. Using optical tweezers and subpixel image analysis, we have recently shown that the radial width of the Z-ring in unconstricted Escherichia coli is about 100 nm. The relatively large width is consistent with the observations of others. Moreover, simulation of the experimental FtsZ distribution using the theoretical three-dimensional (3D) point spread function was strongly in favour of a toroidal rather than a thin cylindrical model of the Z-ring. Here, we show that the low density of FtsZ filaments in the ring coincides within experimental uncertainty with the critical density of a 3D random network of cylindrical sticks. This suggests that the Z-ring may consist of a percolating network of FtsZ filaments. Several factors that are expected to affect the polymerization state and the extent of self-interaction of FtsZ within the Z-ring, as well as the functional implications of its sparse toroidal structure, are discussed in terms of percolation theory.

New approaches to the problem of generating coherent, reproducible phenotypes

Fundamental, unresolved questions in biology include how a bacterium generates coherent phenotypes, how a population of bacteria generates a coherent set of such phenotypes, how the cell cycle is regulated and how life arose. To try to help answer these questions, we have developed the concepts of hyperstructures, competitive coherence and life on the scales of equilibria. Hyperstructures are large assemblies of macromolecules that perform functions. Competitive coherence describes the way in which organisations such as cells select a subset of their constituents to be active in determining their behaviour; this selection results from a competition between a process that is responsible for a historical coherence and another process responsible for coherence with the current environment. Life on the scales of equilibria describes how bacteria depend on the cell cycle to negotiate phenotype space and, in particular, to satisfy the conflicting constraints of having to grow in favourable conditions so as to reproduce yet not grow in hostile conditions so as to survive. Both competitive coherence and life on the scales deal with the problem of reconciling conflicting constraints. Here, we bring together these concepts in the common framework of hyperstructures and make predictions that may be tested using a learning program, Coco, and secondary ion mass spectrometry.

Influence of polyelectrolyte film stiffness on bacterial growth

Photo-cross-linkable polyelectrolyte films, whose nanomechanical properties can be varied under UV light illumination, were prepared from poly(l-lysine) (PLL) and a hyaluronan derivative modified with photoreactive vinylbenzyl groups (HAVB). The adhesion and the growth of two model bacteria, namely Escherichia coli and Lactococcus lactis , were studied on non-cross-linked and cross-linked films to investigate how the film stiffness influences the bacterial behavior. While the Gram positive L. lactis was shown to grow slowly on both films, independently of their rigidity, the Gram negative E. coli exhibited a more rapid growth on non-cross-linked softer films compared to the stiffer ones. Experiments performed on photopatterned films showing both soft and stiff regions, confirmed a faster development of E. coli colonies on softer regions. Interestingly, this behavior is opposite to the one reported before for mammalian cells. Therefore, the photo-cross-linked (PLL/HAVB) films are interesting coatings for tissue engineering since they promote the growth of mammalian cells while limiting the bacterial colonization.

Dynamics of the quorum sensing switch: stochastic and non-stationary effects

Background

A wide range of bacteria species are known to communicate through the so called quorum sensing (QS) mechanism by means of which they produce a small molecule that can freely diffuse in the environment and in the cells. Upon reaching a threshold concentration, the signalling molecule activates the QS-controlled genes that promote phenotypic changes. This mechanism, for its simplicity, has become the model system for studying the emergence of a global response in prokaryotic cells. Yet, how cells precisely measure the signal concentration and act coordinately, despite the presence of fluctuations that unavoidably affects cell regulation and signalling, remains unclear.

Results

We propose a model for the QS signalling mechanism in Vibrio fischeri based on the synthetic strains lux01 and lux02. Our approach takes into account the key regulatory interactions between LuxR and LuxI, the autoinducer transport, the cellular growth and the division dynamics. By using both deterministic and stochastic models, we analyze the response and dynamics at the single-cell level and compare them to the global response at the population level. Our results show how fluctuations interfere with the synchronization of the cell activation and lead to a bimodal phenotypic distribution. In this context, we introduce the concept of precision in order to characterize the reliability of the QS communication process in the colony. We show that increasing the noise in the expression of LuxR helps cells to get activated at lower autoinducer concentrations but, at the same time, slows down the global response. The precision of the QS switch under non-stationary conditions decreases with noise, while at steady-state it is independent of the noise value.

Conclusions

Our in silico experiments show that the response of the LuxR/LuxI system depends on the interplay between non-stationary and stochastic effects and that the burst size of the transcription/translation noise at the level of LuxR controls the phenotypic variability of the population. These results, together with recent experimental evidences on LuxR regulation in wild-type species, suggest that bacteria have evolved mechanisms to regulate the intensity of those fluctuations.

Computational modeling of synthetic microbial biofilms

Microbial biofilms are complex, self-organized communities of bacteria, which employ physiological cooperation and spatial organization to increase both their metabolic efficiency and their resistance to changes in their local environment. These properties make biofilms an attractive target for engineering, particularly for the production of chemicals such as pharmaceutical ingredients or biofuels, with the potential to significantly improve yields and lower maintenance costs. Biofilms are also a major cause of persistent infection, and a better understanding of their organization could lead to new strategies for their disruption. Despite this potential, the design of synthetic biofilms remains a major challenge, due to the complex interplay between transcriptional regulation, intercellular signaling, and cell biophysics. Computational modeling could help to address this challenge by predicting the behavior of synthetic biofilms prior to their construction; however, multiscale modeling has so far not been achieved for realistic cell numbers. This paper presents a computational method for modeling synthetic microbial biofilms, which combines three-dimensional biophysical models of individual cells with models of genetic regulation and intercellular signaling. The method is implemented as a software tool (CellModeller), which uses parallel Graphics Processing Unit architectures to scale to more than 30,000 cells, typical of a 100 μm diameter colony, in 30 min of computation time.

Growth, metabolic partitioning, and the size of microorganisms

Population growth rate is a fundamental ecological and evolutionary characteristic of living organisms, but individuals must balance the metabolism devoted to biosynthesis and reproduction against the maintenance of existing structure and other functionality. Here we present a mathematical model that relates metabolic partitioning to the form of growth. The model captures the observed growth trajectory of single cells and individuals for a variety of species and taxa spanning prokaryotes, unicellular eukaryotes, and small multicellular eukaryotes. Our analysis suggests that the per-unit costs of biosynthesis and maintenance are conserved across prokaryotes and eukaryotes. However, the relative metabolic expenditure on growth and maintenance of whole organisms clearly differentiates taxa: prokaryotes spend an increasing fraction of their entire metabolism on growth with increasing cell size, whereas eukaryotes devote a diminishing fraction. These differences allow us to predict the minimum and maximum size for each taxonomic group, anticipating observed evolutionary life-history transitions. The framework provides energetic insights into taxonomic tradeoffs related to growth and metabolism and constrains traits that are important for size-structured modeling of microbial communities and their ecological and biogeochemical effects.

Timing of Z-ring localization in Escherichia coli

Bacterial cell division takes place in three phases: Z-ring formation at midcell, followed by divisome assembly and building of the septum per se. Using time-lapse microscopy of live bacteria and a high-precision cell edge detection method, we have previously found the true time for the onset of septation, τ(c), and the time between consecutive divisions, τ(g). Here, we combine the above method with measuring the dynamics of the FtsZ-GFP distribution in individual Escherichia coli cells to determine the Z-ring positioning time, τ(z). To analyze the FtsZ-GFP distribution along the cell, we used the integral fluorescence profile (IFP), which was obtained by integrating the fluorescence intensity across the cell width. We showed that the IFP may be approximated by an exponential peak and followed the peak evolution throughout the cell cycle, to find a quantitative criterion for the positioning of the Z-ring and hence the value of τ(z). We defined τ(z) as the transition from oscillatory to stable behavior of the mean IFP position. This criterion was corroborated by comparison of the experimental results to a theoretical model for the FtsZ dynamics, driven by Min oscillations. We found that τ(z) < τ(c) for all the cells that were analyzed. Moreover, our data suggested that τ(z) is independent of τ(c), τ(g) and the cell length at birth, L(0). These results are consistent with the current understanding of the Z-ring positioning and cell septation processes.

Condition-dependent cell volume and concentration of Escherichia coli to facilitate data conversion for systems biology modeling

Systems biology modeling typically requires quantitative experimental data such as intracellular concentrations or copy numbers per cell. In order to convert population-averaging omics measurement data to intracellular concentrations or cellular copy numbers, the total cell volume and number of cells in a sample need to be known. Unfortunately, even for the often studied model bacterium Escherichia coli this information is hardly available and furthermore, certain measures (e.g. cell volume) are also dependent on the growth condition. In this work, we have determined these basic data for E. coli cells when grown in 22 different conditions so that respective data conversions can be done correctly. First, we determine growth-rate dependent cell volumes. Second, we show that in a 1 ml E. coli sample at an optical density (600 nm) of 1 the total cell volume is around 3.6 µl for all conditions tested. Third, we demonstrate that the cell number in a sample can be determined on the basis of the sample's optical density and the cells' growth rate. The data presented will allow for conversion of E. coli measurement data normalized to optical density into volumetric cellular concentrations and copy numbers per cell--two important parameters for systems biology model development.

Self-organized partitioning of dynamically localized proteins in bacterial cell division

The Min proteins are equally partitioned between daughter cells at division.

The mechanism allowing this accurate distribution is intrinsic to the Min system.

Individual oscillations appear in each daughter cell before cytokinesis is completed.

Diffusion through the gradually constricting septum is key to this process.

One of the central problems of cell division is the proper distribution of all components to the progeny, which is essential to avoid the adverse effects that an unequal distribution—when not actively sought for differentiation purposes—would have on cell growth and regulation. Fast-growing bacterial cells are particularly exposed to this problem, as corrections of inequalities in protein distribution by biosynthesis could be too slow compared with the generation time. Moreover, bacterial proteins are usually stable and, therefore, their levels are not easily adjustable in one generation. Although for homogeneously distributed proteins an equal partitioning at division is readily achieved, dedicated mechanisms must exist to segregate proteins or cellular structures that possess a specific cellular location, but these mechanisms are largely unknown.

An extremely challenging case is represented by the Min proteins—MinC, MinD and MinE—that in Escherichia coli oscillate from pole to pole to inhibit the assembly of the cytokinetic ring anywhere except at mid-cell. The oscillations stem solely from local interactions among the proteins at the cytoplasmic membrane. In this work, we show that self-organization is also responsible for the distribution of Min proteins between daughter cells at division. Our combined experimental and computational results demonstrate that the equal protein partitioning stems from interplay between the self-organized oscillations and changes in the cell geometry during division, with no need for any additional regulatory network.

Using high-resolution time-lapse microscopy, we detected changes in the Min oscillatory regime that correlate with the amount of septal constriction (Figure 3A, B, E and F). When the cell is unconstricted, oscillations run from pole to pole (Figure 3A). When the constriction reaches a certain degree, typically corresponding to a septum of 600–500 nm, the oscillations change into a ‘half-cell to half-cell' mode during which the fluorescence covers, alternatively, the entire membrane of one daughter cell (Figure 3A, B and E). This mode persists for several minutes and, just before cell division when the septum is smaller than 200 nm, gives way to yet another oscillatory pattern wherein oscillations split and run independently in each daughter cell (Figure 3A, B and F).

Our 3D stochastic computer simulations revealed that these different regimes are an outcome of impaired diffusion through the closing septum and that oscillations finally split because protein exchange between the two future daughter cells becomes critically slow, so that independent oscillations on both sides of the septum become the stable solution (Figure 6A and E). FRAP experiments confirmed that the presence of the septum is enough to slow down the passage of molecules from one side of the cell to the other (Figure 6F). As oscillations become independent in each daughter cell before completion of cytokinesis, diffusion through the septum can still occur, which further equilibrates the levels of the Min proteins in the daughter cells (Figure 3C and D and Figure 6B, C and D).

In summary, our results suggest that E. coli cells have evolved a very simple and elegant way to ensure equal concentrations of the Min proteins in the progeny, based entirely on the intrinsic self-organizing properties of the Min system. This provides a clear example of self-organizing partitioning, which we expect to be a widely used strategy given its simplicity and low evolutionary cost.

How cells manage to get equal distribution of their structures and molecules at cell division is a crucial issue in biology. In principle, a feedback mechanism could always ensure equality by measuring and correcting the distribution in the progeny. However, an elegant alternative could be a mechanism relying on self-organization, with the interplay between system properties and cell geometry leading to the emergence of equal partitioning. The problem is exemplified by the bacterial Min system that defines the division site by oscillating from pole to pole. Unequal partitioning of Min proteins at division could negatively impact system performance and cell growth because of loss of Min oscillations and imprecise mid-cell determination. In this study, we combine live cell and computational analyses to show that known properties of the Min system together with the gradual reduction of protein exchange through the constricting septum are sufficient to explain the observed highly precise spontaneous protein partitioning. Our findings reveal a novel and effective mechanism of protein partitioning in dividing cells and emphasize the importance of self-organization in basic cellular processes.

Pre-dispositions and epigenetic inheritance in the Escherichia coli lactose operon bistable switch

The lactose operon regulation in Escherichia coli is a primary model of phenotypic switching, reminiscent of cell fate determination in higher organisms. Under conditions of bistability, an isogenic cell population partitions into two subpopulations, with the operon's genes turned on or remaining off. It is generally hypothesized that the final state of a cell depends solely on stochastic fluctuations of the network's protein concentrations, particularly on bursts of lactose permease expression. Nevertheless, the mechanisms underlying the cell switching decision are not fully understood. We designed a microfluidic system to follow the formation of a transiently bimodal population within growing microcolonies. The analysis of genealogy and cell history revealed the existence of pre-disposing factors for switching that are epigenetically inherited. Both the pre-induction expression stochasticity of the lactose operon repressor LacI and the cellular growth rate are predictive factors of the cell's response upon induction, with low LacI concentration and slow growth correlating with higher switching probability. Thus, stochasticity at the local level of the network and global physiology are synergistically involved in cell response determination.

Using buoyant mass to measure the growth of single cells

We used a suspended microchannel resonator (SMR) combined with picoliter-scale microfluidic control to measure buoyant mass and determine the ‘instantaneous’ growth rates of individual cells. The SMR measures mass with femtogram precision, allowing rapid determination of the growth rate in a fraction of a complete cell cycle. We found that for individual cells of Bacillus subtilis, Escherichia coli, Saccharomyces cerevisiae and mouse lymphoblasts, heavier cells grew faster than lighter cells.