Temporal controls of the asymmetric cell division cycle in Caulobacter crescentus

Modeling the temporal dynamics of master regulators and CtrA proteolysis in Caulobacter crescentus cell cycle

The cell cycle of Caulobacter crescentus involves the polar morphogenesis and an asymmetric cell division driven by precise interactions and regulations of proteins, which makes Caulobacter an ideal model organism for investigating bacterial cell development and differentiation. The abundance of molecular data accumulated on Caulobacter motivates system biologists to analyze the complex regulatory network of cell cycle via quantitative modeling. In this paper, We propose a comprehensive model to accurately characterize the underlying mechanisms of cell cycle regulation based on the study of: a) chromosome replication and methylation; b) interactive pathways of five master regulatory proteins including DnaA, GcrA, CcrM, CtrA, and SciP, as well as novel consideration of their corresponding mRNAs; c) cell cycle-dependent proteolysis of CtrA through hierarchical protease complexes. The temporal dynamics of our simulation results are able to closely replicate an extensive set of experimental observations and capture the main phenotype of seven mutant strains of Caulobacter crescentus. Collectively, the proposed model can be used to predict phenotypes of other mutant cases, especially for nonviable strains which are hard to cultivate and observe. Moreover, the module of cyclic proteolysis is an efficient tool to study the metabolism of proteins with similar mechanisms.

Timed cellular events in both eukaryotes and prokaryotes, such as chromosome replication, transcription, cell differentiation, cytokinesis, and cell division, are controlled by remarkably complex genetic regulations and protein-protein interactions. In this work, we investigate the cell cycle of Caulobacter crescentus, an alphaproteobacterium undergoing asymmetric cell divisions, to understand mechanisms underlying temporal regulations of complex cellular events. The asymmetric lifestyle makes Caulobacter crescentus easily synchronized and tracked, which is the foundation of molecular data accumulation. Here, we utilize the mathematical modeling together with experimental information to systematically integrate the complex gene-protein and protein-protein interactions in cell cycle progression. Using the mathematical model, we capture core features of cell cycle-dependent methylation, transcription, and proteolysis. In mutant cases, we found the complex and redundant regulatory network ensure the robustness of Caulobacter crescentus system because the change of most molecules does not cause immediate mortality, although they influence the time points of cell differentiation and division. The overall model and individual modules such as simulating transcriptional regulations and protease complexes can be further extended to the study of cell development in other bacterial species.

Computational modeling of unphosphorylated CtrA:Cori binding in the Caulobacter cell cycle

In the alphaproteobacterium, Caulobacter crescentus, phosphorylated CtrA (CtrA∼P), a master regulatory protein, binds directly to the chromosome origin (Cori) to inhibit DNA replication. Using a mathematical model of CtrA binding at Cori site [d], we provide computational evidence that CtrA_U can displace CtrA∼P from Cori at the G1-S transition. Investigation of this interaction within a detailed model of the C. crescentus cell cycle suggests that CckA phosphatase may clear Cori of CtrA∼P by altering the [CtrA_U]/[CtrA∼P] ratio rather than by completely depleting CtrA∼P. Model analysis reveals that the mechanism allows for a speedier transition into S phase, stabilizes the timing of chromosome replication under fluctuating rates of CtrA proteolysis, and may contribute to the viability of numerous mutant strains. Overall, these results suggest that CtrA_U enhances the robustness of chromosome replication. More generally, our proposed regulation of CtrA:Cori dynamics may represent a novel motif for molecular signaling in cell physiology.

High-Temperature Live-Cell Imaging of Cytokinesis, Cell Motility, and Cell-Cell Interactions in the Thermoacidophilic Crenarchaeon Sulfolobus acidocaldarius

Significant technical challenges have limited the study of extremophile cell biology. Here we describe a system for imaging samples at 75°C using high numerical aperture, oil-immersion lenses. With this system we observed and quantified the dynamics of cell division in the model thermoacidophilic crenarchaeon Sulfolobus acidocaldarius with unprecedented resolution. In addition, we observed previously undescribed dynamic cell shape changes, cell motility, and cell-cell interactions, shedding significant new light on the high-temperature lifestyle of this organism.

Cell cycle control and environmental response by second messengers in Caulobacter crescentus

Background

Second messengers, c-di-GMP and (p)ppGpp, are vital regulatory molecules in bacteria, influencing cellular processes such as biofilm formation, transcription, virulence, quorum sensing, and proliferation. While c-di-GMP and (p)ppGpp are both synthesized from GTP molecules, they play antagonistic roles in regulating the cell cycle. In C. crescentus, c-di-GMP works as a major regulator of pole morphogenesis and cell development. It inhibits cell motility and promotes S-phase entry by inhibiting the activity of the master regulator, CtrA. Intracellular (p)ppGpp accumulates under starvation, which helps bacteria to survive under stressful conditions through regulating nucleotide levels and halting proliferation. (p)ppGpp responds to nitrogen levels through RelA-SpoT homolog enzymes, detecting glutamine concentration using a nitrogen phosphotransferase system (PTS ^Ntr). This work relates the guanine nucleotide-based second messenger regulatory network with the bacterial PTS ^Ntr system and investigates how bacteria respond to nutrient availability.

Results

We propose a mathematical model for the dynamics of c-di-GMP and (p)ppGpp in C. crescentus and analyze how the guanine nucleotide-based second messenger system responds to certain environmental changes communicated through the PTS ^Ntr system. Our mathematical model consists of seven ODEs describing the dynamics of nucleotides and PTS ^Ntr enzymes. Our simulations are consistent with experimental observations and suggest, among other predictions, that SpoT can effectively decrease c-di-GMP levels in response to nitrogen starvation just as well as it increases (p)ppGpp levels. Thus, the activity of SpoT (or its homologues in other bacterial species) can likely influence the cell cycle by influencing both c-di-GMP and (p)ppGpp.

Conclusions

In this work, we integrate current knowledge and experimental observations from the literature to formulate a novel mathematical model. We analyze the model and demonstrate how the PTS ^Ntr system influences (p)ppGpp, c-di-GMP, GMP and GTP concentrations. While this model does not consider all aspects of PTS ^Ntr signaling, such as cross-talk with the carbon PTS system, here we present our first effort to develop a model of nutrient signaling in C. crescentus.

Stochastic modeling and simulation of reaction-diffusion system with Hill function dynamics

Background

Stochastic simulation of reaction-diffusion systems presents great challenges for spatiotemporal biological modeling and simulation. One widely used framework for stochastic simulation of reaction-diffusion systems is reaction diffusion master equation (RDME). Previous studies have discovered that for the RDME, when discretization size approaches zero, reaction time for bimolecular reactions in high dimensional domains tends to infinity.

Results

In this paper, we demonstrate that in the 1D domain, highly nonlinear reaction dynamics given by Hill function may also have dramatic change when discretization size is smaller than a critical value. Moreover, we discuss methods to avoid this problem: smoothing over space, fixed length smoothing over space and a hybrid method.

Conclusion

Our analysis reveals that the switch-like Hill dynamics reduces to a linear function of discretization size when the discretization size is small enough. The three proposed methods could correctly (under certain precision) simulate Hill function dynamics in the microscopic RDME system.

3D time-varying simulations of Ca2+ dynamics in arterial coupled cells: A massively parallel implementation

Preferential locations of atherosclerotic plaque are strongly associated with the areas of low wall shear stress and disturbed haemodynamic characteristics such as flow detachment, flow recirculation and oscillatory flow. The areas of low wall shear stress are also associated with the reduced production of adenosine triphosphate in the endothelial layer, as well as the resulting reduced production of inositol trisphosphate (IP³ ). The subsequent variation in Ca²⁺ signalling and nitric oxide synthesis could lead to the impairment of the atheroprotective function played by nitric oxide. In previous studies, it has been suggested that the reduced IP³ and Ca²⁺ signalling can explain the correlation of atherosclerosis with induced low WSS and disturbed flow characteristics. The massively parallel implementation described in this article provides insight into the dynamics of coupled smooth muscle cells and endothelial cells mapped onto the surface of an idealised arterial bifurcation. We show that variations in coupling parameters, which model normal and pathological conditions, provide vastly different smooth muscle cell Ca²⁺ dynamics and wave propagation profiles. The extensibility of the coupled cells model and scalability of the implementation provide a solid framework for in silico investigations of the interaction between complex cellular chemistry and the macro-scale processes determined by fluid dynamics. © 2016 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.

Dynamical Localization of DivL and PleC in the Asymmetric Division Cycle of Caulobacter crescentus: A Theoretical Investigation of Alternative Models

Cell-fate asymmetry in the predivisional cell of Caulobacter crescentus requires that the regulatory protein DivL localizes to the new pole of the cell where it up-regulates CckA kinase, resulting in a gradient of CtrA~P across the cell. In the preceding stage of the cell cycle (the "stalked" cell), DivL is localized uniformly along the cell membrane and maintained in an inactive form by DivK~P. It is unclear how DivL overcomes inhibition by DivK~P in the predivisional cell simply by changing its location to the new pole. It has been suggested that co-localization of DivL with PleC phosphatase at the new pole is essential to DivL's activity there. However, there are contrasting views on whether the bifunctional enzyme, PleC, acts as a kinase or phosphatase at the new pole. To explore these ambiguities, we formulated a mathematical model of the spatiotemporal distributions of DivL, PleC and associated proteins (DivJ, DivK, CckA, and CtrA) during the asymmetric division cycle of a Caulobacter cell. By varying localization profiles of DivL and PleC in our model, we show how the physiologically observed spatial distributions of these proteins are essential for the transition from a stalked cell to a predivisional cell. Our simulations suggest that PleC is a kinase in predivisional cells, and that, by sequestering DivK~P, the kinase form of PleC enables DivL to be reactivated at the new pole. Hence, co-localization of PleC kinase and DivL is essential to establishing cellular asymmetry. Our simulations reproduce the experimentally observed spatial distribution and phosphorylation status of CtrA in wild-type and mutant cells. Based on the model, we explore novel combinations of mutant alleles, making predictions that can be tested experimentally.

Dynamical modeling of the cell cycle and cell fate emergence in Caulobacter crescentus

The division of Caulobacter crescentus, a model organism for studying cell cycle and differentiation in bacteria, generates two cell types: swarmer and stalked. To complete its cycle, C. crescentus must first differentiate from the swarmer to the stalked phenotype. An important regulator involved in this process is CtrA, which operates in a gene regulatory network and coordinates many of the interactions associated to the generation of cellular asymmetry. Gaining insight into how such a differentiation phenomenon arises and how network components interact to bring about cellular behavior and function demands mathematical models and simulations. In this work, we present a dynamical model based on a generalization of the Boolean abstraction of gene expression for a minimal network controlling the cell cycle and asymmetric cell division in C. crescentus. This network was constructed from data obtained from an exhaustive search in the literature. The results of the simulations based on our model show a cyclic attractor whose configurations can be made to correspond with the current knowledge of the activity of the regulators participating in the gene network during the cell cycle. Additionally, we found two point attractors that can be interpreted in terms of the network configurations directing the two cell types. The entire network is shown to be operating close to the critical regime, which means that it is robust enough to perturbations on dynamics of the network, but adaptable to environmental changes.

Computational and genetic reduction of a cell cycle to its simplest, primordial components

What are the minimal requirements to sustain an asymmetric cell cycle? Here we use mathematical modelling and forward genetics to reduce an asymmetric cell cycle to its simplest, primordial components. In the Alphaproteobacterium Caulobacter crescentus, cell cycle progression is believed to be controlled by a cyclical genetic circuit comprising four essential master regulators. Unexpectedly, our in silico modelling predicted that one of these regulators, GcrA, is in fact dispensable. We confirmed this experimentally, finding that ΔgcrA cells are viable, but slow-growing and elongated, with the latter mostly due to an insufficiency of a key cell division protein. Furthermore, suppressor analysis showed that another cell cycle regulator, the methyltransferase CcrM, is similarly dispensable with simultaneous gcrA/ccrM disruption ameliorating the cytokinetic and growth defect of ΔgcrA cells. Within the Alphaproteobacteria, gcrA and ccrM are consistently present or absent together, rather than either gene being present alone, suggesting that gcrA/ccrM constitutes an independent, dispensable genetic module. Together our approaches unveil the essential elements of a primordial asymmetric cell cycle that should help illuminate more complex cell cycles.

Potential role of a bistable histidine kinase switch in the asymmetric division cycle of Caulobacter crescentus

The free-living aquatic bacterium, Caulobacter crescentus, exhibits two different morphologies during its life cycle. The morphological change from swarmer cell to stalked cell is a result of changes of function of two bi-functional histidine kinases, PleC and CckA. Here, we describe a detailed molecular mechanism by which the function of PleC changes between phosphatase and kinase state. By mathematical modeling of our proposed molecular interactions, we derive conditions under which PleC, CckA and its response regulators exhibit bistable behavior, thus providing a scenario for robust switching between swarmer and stalked states. Our simulations are in reasonable agreement with in vitro and in vivo experimental observations of wild type and mutant phenotypes. According to our model, the kinase form of PleC is essential for the swarmer-to-stalked transition and to prevent premature development of the swarmer pole. Based on our results, we reconcile some published experimental observations and suggest novel mutants to test our predictions.

Transcriptomic and phylogenetic analysis of a bacterial cell cycle reveals strong associations between gene co-expression and evolution

BACKGROUND

The genetic network involved in the bacterial cell cycle is poorly understood even though it underpins the remarkable ability of bacteria to proliferate. How such network evolves is even less clear. The major aims of this work were to identify and examine the genes and pathways that are differentially expressed during the Caulobacter crescentus cell cycle, and to analyze the evolutionary features of the cell cycle network.

RESULTS

We used deep RNA sequencing to obtain high coverage RNA-Seq data of five C. crescentus cell cycle stages, each with three biological replicates. We found that 1,586 genes (over a third of the genome) display significant differential expression between stages. This gene list, which contains many genes previously unknown for their cell cycle regulation, includes almost half of the genes involved in primary metabolism, suggesting that these "house-keeping" genes are not constitutively transcribed during the cell cycle, as often assumed. Gene and module co-expression clustering reveal co-regulated pathways and suggest functionally coupled genes. In addition, an evolutionary analysis of the cell cycle network shows a high correlation between co-expression and co-evolution. Most co-expression modules have strong phylogenetic signals, with broadly conserved genes and clade-specific genes predominating different substructures of the cell cycle co-expression network. We also found that conserved genes tend to determine the expression profile of their module.

CONCLUSION

We describe the first phylogenetic and single-nucleotide-resolution transcriptomic analysis of a bacterial cell cycle network. In addition, the study suggests how evolution has shaped this network and provides direct biological network support that selective pressure is not on individual genes but rather on the relationship between genes, which highlights the importance of integrating phylogenetic analysis into biological network studies.

Spatial perturbation with synthetic protein scaffold reveals robustness of asymmetric cell division

Asymmetric cell division is an important mechanism for creating diversity in a cellular population. Stem cells commonly perform asymmetric division to generate both a daughter stem cell for self-renewal and a more differentiated daughter cell to populate the tissue. During asymmetric cell division, protein cell fate determinants asymmetrically localize to the opposite poles of a dividing cell to cause distinct cell fate. However, it remains unclear whether cell fate determination is robust to fluctuations and noise during this spatial allocation process. To answer this question, we engineered Caulobacter, a bacterial model for asymmetric division, to express synthetic scaffolds with modular protein interaction domains. These scaffolds perturbed the spatial distribution of the PleC-DivJ-DivK phospho-signaling network without changing their endogenous expression levels. Surprisingly, enforcing symmetrical distribution of these cell fate determinants did not result in symmetric daughter fate or any morphological defects. Further computational analysis suggested that PleC and DivJ form a robust phospho-switch that can tolerate high amount of spatial variation. This insight may shed light on the presence of similar phospho-switches in stem cell asymmetric division regulation. Overall, our study demonstrates that synthetic protein scaffolds can provide a useful tool to probe biological systems for better understanding of their operating principles.

Mathematical modeling of fission yeast Schizosaccharomyces pombe cell cycle: exploring the role of multiple phosphatases

Cell cycle is the central process that regulates growth and division in all eukaryotes. Based on the environmental condition sensed, the cell lies in a resting phase G0 or proceeds through the cyclic cell division process (G1→S→G2→M). These series of events and phase transitions are governed mainly by the highly conserved Cyclin dependent kinases (Cdks) and its positive and negative regulators. The cell cycle regulation of fission yeast Schizosaccharomyces pombe is modeled in this study. The study exploits a detailed molecular interaction map compiled based on the published model and experimental data. There are accumulating evidences about the prominent regulatory role of specific phosphatases in cell cycle regulations. The current study emphasizes the possible role of multiple phosphatases that governs the cell cycle regulation in fission yeast S. pombe. The ability of the model to reproduce the reported regulatory profile for the wild-type and various mutants was verified though simulations.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11693-011-9090-7) contains supplementary material, which is available to authorized users.

Phase resetting reveals network dynamics underlying a bacterial cell cycle

Genomic and proteomic methods yield networks of biological regulatory interactions but do not provide direct insight into how those interactions are organized into functional modules, or how information flows from one module to another. In this work we introduce an approach that provides this complementary information and apply it to the bacterium Caulobacter crescentus, a paradigm for cell-cycle control. Operationally, we use an inducible promoter to express the essential transcriptional regulatory gene ctrA in a periodic, pulsed fashion. This chemical perturbation causes the population of cells to divide synchronously, and we use the resulting advance or delay of the division times of single cells to construct a phase resetting curve. We find that delay is strongly favored over advance. This finding is surprising since it does not follow from the temporal expression profile of CtrA and, in turn, simulations of existing network models. We propose a phenomenological model that suggests that the cell-cycle network comprises two distinct functional modules that oscillate autonomously and couple in a highly asymmetric fashion. These features collectively provide a new mechanism for tight temporal control of the cell cycle in C. crescentus. We discuss how the procedure can serve as the basis for a general approach for probing network dynamics, which we term chemical perturbation spectroscopy (CPS).

During the cell cycle, the cell progresses through a series of stages that are associated with various cell cycle events such as replication of genetic materials. Genetic and molecular dissections have revealed that the cell cycle is regulated by a network of interacting molecules that produces oscillatory dynamics. The major cell cycle regulators have been identified previously in different species and the activity of these regulators oscillates. However, the question of how cell cycle regulators coordinate different cell cycle events during the cell cycle remains controversial. Here, we investigate this question in a model bacterial system for cell cycle, Caulobacter crescentus. We perturb the expression of the master cell cycle regulator ctrA in a pulsatile fashion and quantify the response of the cell cycle to such perturbations. The measured response is contradictory to the existing mechanism of Caulobacter cell cycle control, which views the cell cycle progression as a sequential activation/inhibition process. We propose a new model that involves coupling of multiple oscillators and show the quantitative agreement between this new model and our measurements. We expect this procedure to be generalized and applied to a broad range of systems to obtain information that complements that obtained from other methods.

Computational modeling of signaling pathways mediating cell cycle checkpoint control and apoptotic responses to ionizing radiation-induced DNA damage

The shape of dose response of ionizing radiation (IR) induced cancer at low dose region, either linear non-threshold or J-shaped, has been a debate for a long time. This dose response relationship can be influenced by built-in capabilities of cells that minimize the fixation of IR-mediated DNA damage as pro-carcinogenic mutations. Key capabilities include sensing of damage, activation of cell cycle checkpoint arrests that provide time needed for repair of the damage as well as apoptosis. Here we describe computational modeling of the signaling pathways that link sensing of DNA damage and checkpoint arrest activation/apoptosis to investigate how these molecular-level interactions influence the dose response relationship for IR induced cancer. The model provides qualitatively accurate descriptions of the IR-mediated activation of cell cycle checkpoints and the apoptotic pathway, and of time-course activities and dose response of relevant regulatory proteins (e.g. p53 and p21). Linking to a two-stage clonal growth cancer model, the model described here successfully captured a monotonically increasing to a J-shaped dose response curve and identified one potential mechanism leading to the J-shape: the cell cycle checkpoint arrest time saturates with the increase of the dose.

Regulation of chromosomal replication in Caulobacter crescentus

The alpha-proteobacterium Caulobacter crescentus is characterized by its asymmetric cell division, which gives rise to a replicating stalked cell and a non-replicating swarmer cell. Thus, the initiation of chromosomal replication is tightly regulated, temporally and spatially, to ensure that it is coordinated with cell differentiation and cell cycle progression. Waves of DnaA and CtrA activities control when and where the initiation of DNA replication will take place in C. crescentus cells. The conserved DnaA protein initiates chromosomal replication by directly binding to sites within the chromosomal origin (Cori), ensuring that DNA replication starts once and only once per cell cycle. The CtrA response regulator represses the initiation of DNA replication in swarmer cells and in the swarmer compartment of pre-divisional cells, probably by competing with DnaA for binding to Cori. CtrA and DnaA are controlled by multiple redundant regulatory pathways that include DNA methylation-dependent transcriptional regulation, temporally regulated proteolysis and the targeting of regulators to specific locations within the cell. Besides being critical regulators of chromosomal replication, CtrA and DnaA are also master transcriptional regulators that control the expression of many genes, thus connecting DNA replication with other events of the C. crescentus cell cycle.

Single-gene tuning of Caulobacter cell cycle period and noise, swarming motility, and surface adhesion

We established that the sensor histidine kinase DivJ has an important role in the regulation of C. crescentus cell cycle period and noise. This was accomplished by designing and conducting single-cell experiments to probe the dependence of cell cycle noise on divJ expression and constructing a simplified cell cycle model that captures the dependence of cell cycle noise on DivJ with molecular details.

In addition to its role in regulating the cell cycle, DivJ also affects polar cell development in C. crescentus, regulating swarming motility and surface adhesion. We propose that pleiotropic control of polar cell development by the DivJ–DivK–PleC signaling pathway underlies divJ-dependent tuning of cell swarming and adhesion behaviors.

We have integrated the study of single-cell fluorescence dynamics with a kinetic model simulation to provide direct quantitative evidence that the DivJ histidine kinase is localized to the cell pole through a dynamic diffusion-and-capture mechanism during the C. crescentus cell cycle.

Temporally-coordinated localization of various structural and signaling proteins is critical for proper cell cycle regulation and polar cell development in the bacterium, Caulobacter crescentus. Included among these dynamically-localized regulatory proteins is the sensor histidine kinase, DivJ (Wheeler and Shapiro, 1999). Co-localized with DivJ in the early stalked phase is the phosphorylated response regulator DivK∼P (Jacobs et al, 2001), and the protease ClpXP (McGrath et al, 2006), which degrades the master cell cycle regulator, CtrA (Jenal and Fuchs, 1998). Recent single-cell measurements of surface attached C. crescentus cells have revealed an intriguing role for DivJ in the control of noise in cell division period (Siegal-Gaskins and Crosson, 2008). The noise of the cell cycle increases significantly upon disruption of the divJ gene, with a relatively small accompanying increase in the mean cell cycle time. The deterministic nature of the existing cell cycle models (Li et al, 2008, 2009; Shen et al, 2008) cannot explain the measured increase in cell cycle period and noise in a divJ null strain. Moreover, mechanistic descriptions of how DivJ and its signaling partners are localized and how these proteins underlie the control of polar cell development and cell adhesion in C. crescentus remain immature.

The single-cell experiments and analysis presented herein reveal that C. crescentus cell cycle period and noise can be tuned by DivJ (Figure 2). Specifically, in the case of low (or no) divJ expression the cell cycle is perturbed, and this is quantified by way of the (measured) noise in the cell cycle period. The level of noise is readily controlled through regulated expression of the divJ gene (Figure 2B). A simplified protein interaction network of stalked C. crescentus cell cycle regulation involving minimal components (CtrA, CtrA∼P, DivK, DivK∼P, and DivJ) was constructed to explore such tunability at the molecular level. The agreement of our model with our (and other) experiments suggests this simplified protein regulatory network is sufficient to explain the major features of the C. crescentus cell cycle. Indeed, stochastic simulations of this model using the Gillespie method (Gillespie, 1976) establish the importance of robust DivJ-mediated phosphorylation of its cognate receiver protein, DivK, in regulating the variance of cell cycle oscillations. Increased variability in the concentration of DivK∼P at the single cell level under divJ depletion subsequently leads to increased noise in the regulation of CtrA phosphorylation and degradation. Our experiments and simulations provide evidence that the steady state level of DivK∼P at the single-cell level (as maintained by DivJ) is essential in maintaining regular timing of the cell division period in C. crescentus.

In addition to its role in regulating cell cycle, divJ expression also affects polar cell development in C. crescentus. Specifically, the capacity of swarmer cells to adhere to a glass surface is suppressed at high levels of divJ expression. The effect of elevated divJ expression on the adhesive capacity of the cell is reflected in a reduced rate of two-dimensional biofilm formation. This effect is quantitatively captured by our mathematical model that relates single-cell surface adhesion physiology and biofilm formation dynamics. This result, and our observation that divJ expression tunes swarming motility in semi-solid growth medium, suggests a model in which increased DivJ concentration in the swarmer compartment (due to constitutive overexpression) ultimately results in improper development of polar organelles that are required for adhesion of swarming motility.

Despite the appreciated significance of protein localization for bacterial physiological functions, the molecular mechanism of how polar protein localization is achieved has only been tested in a few cases (Shapiro et al, 2002; Thanbichler and Shapiro, 2008). Mechanisms such as the polar insertion model and diffusion-and-capture have been proposed but the community's knowledge is limited to very few examples (Charles et al, 2001; Rudner et al, 2002). We provide direct evidence from experiments and simulations that the DivJ histidine kinase becomes localized to the cell pole through a dynamic diffusion-and-capture mechanism during the C. crescentus cell cycle (Figure 7). We show that a kinetic model based on a Langmuir adsorption/desorption relationship (Figure 7D) is sufficient to explain the time evolution of the single cell fluorescence time traces (Figure 7C and E) and allows establishing quantitative correspondences between the simulated dynamics and experimentally determined DivJ–EGFP dynamics. This localization mechanism is consistent with a diffusion-and-capture model. In short, the model posits that proteins are randomly distributed and are freely diffusing until they are captured at the site where they ultimately reside (Rudner et al, 2002; Shapiro et al, 2002; Bardy and Maddock, 2007). With a diffusion-and-capture pathway, it has been argued that proteins can be adsorbed either dynamically or statically (Shapiro et al, 2009). Our analysis of DivJ–EGFP in single cells supports a dynamic diffuse-and-capture mechanism for DivJ localization.

Sensor histidine kinases underlie the regulation of a range of physiological processes in bacterial cells, from chemotaxis to cell division. In the gram-negative bacterium Caulobacter crescentus, the membrane-bound histidine kinase, DivJ, is a polar-localized regulator of cell cycle progression and development. We show that DivJ localizes to the cell pole through a dynamic diffusion and capture mechanism rather than by active localization. Analysis of single C. crescentus cells in microfluidic culture demonstrates that controlled expression of divJ permits facile tuning of both the mean and noise of the cell division period. Simulations of the cell cycle that use a simplified protein interaction network capture previously measured oscillatory protein profiles, and recapitulate the experimental observation that deletion of divJ increases the cell cycle period and noise. We further demonstrate that surface adhesion and swarming motility of C. crescentus in semi-solid media can also be tuned by divJ expression. We propose a model in which pleiotropic control of polar cell development by the DivJ–DivK–PleC signaling pathway underlies divJ-dependent tuning of cell swarming and adhesion behaviors.

High-throughput, subpixel precision analysis of bacterial morphogenesis and intracellular spatio-temporal dynamics

Bacteria display various shapes and rely on complex spatial organization of their intracellular components for many cellular processes. This organization changes in response to internal and external cues. Quantitative, unbiased study of these spatio-temporal dynamics requires automated image analysis of large microscopy datasets. We have therefore developed MicrobeTracker, a versatile and high-throughput image analysis program that outlines and segments cells with subpixel precision, even in crowded images and mini-colonies, enabling cell lineage tracking. MicrobeTracker comes with an integrated accessory tool, SpotFinder, which precisely tracks foci of fluorescently labelled molecules inside cells. Using MicrobeTracker, we discover that the dynamics of the extensively studied Escherichia coli Min oscillator depends on Min protein concentration, unveiling critical limitations in robustness within the oscillator. We also find that the fraction of MinD proteins oscillating increases with cell length, indicating that the oscillator has evolved to be most effective when cells attain an appropriate length. MicrobeTracker was also used to uncover novel aspects of morphogenesis and cell cycle regulation in Caulobacter crescentus. By tracking filamentous cells, we show that the chromosomal origin at the old-pole is responsible for most replication/separation events while the others remain largely silent despite contiguous cytoplasm. This surprising position-dependent silencing is regulated by division.

Modeling Dose-response at Low Dose: A Systems Biology Approach for Ionization Radiation

For ionization radiation (IR) induced cancer, a linear non-threshold (LNT) model at very low doses is the default used by a number of national and international organizations and in regulatory law. This default denies any positive benefit from any level of exposure. However, experimental observations and theoretical biology have found that both linear and J-shaped IR dose-response curves can exist at those very low doses. We develop low dose J-shaped dose-response, based on systems biology, and thus justify its use regarding exposure to IR. This approach incorporates detailed, molecular and cellular descriptions of biological/toxicological mechanisms to develop a dose-response model through a set of nonlinear, differential equations describing the signaling pathways and biochemical mechanisms of cell cycle checkpoint, apoptosis, and tumor incidence due to IR. This approach yields a J-shaped dose response curve while showing where LNT behaviors are likely to occur. The results confirm the hypothesis of the J-shaped dose response curve: the main reason is that, at low-doses of IR, cells stimulate protective systems through a longer cell arrest time per unit of IR dose. We suggest that the policy implications of this approach are an increasingly correct way to deal with precautionary measures in public health.