Physics of Intracellular Organization in Bacteria

Contribution of different macromolecules to the diffusion of a 40 nm particle in Escherichia coli

Due to the high concentration of proteins, nucleic acids, and other macromolecules, the bacterial cytoplasm is typically described as a crowded environment. However, the extent to which each of these macromolecules individually affects the mobility of macromolecular complexes, and how this depends on growth conditions, is presently unclear. In this study, we sought to quantify the crowding experienced by an exogenous 40 nm fluorescent particle in the cytoplasm of E. coli under different growth conditions. By performing single-particle tracking measurements in cells selectively depleted of DNA and/or mRNA, we determined the contribution to crowding of mRNA, DNA, and remaining cellular components, i.e., mostly proteins and ribosomes. To estimate this contribution to crowding, we quantified the difference of the particle's diffusion coefficient in conditions with and without those macromolecules. We found that the contributions of the three classes of components were of comparable magnitude, being largest in the case of proteins and ribosomes. We further found that the contributions of mRNA and DNA to crowding were significantly larger than expected based on their volumetric fractions alone. Finally, we found that the crowding contributions change only slightly with the growth conditions. These results reveal how various cellular components partake in crowding of the cytoplasm and the consequences this has for the mobility of large macromolecular complexes.

The quantitative basis for the redistribution of immobile bacterial lipoproteins to division septa

The spatial localisation of proteins is critical for most cellular function. In bacteria, this is typically achieved through capture by established landmark proteins. However, this requires that the protein is diffusive on the appropriate timescale. It is therefore unknown how the localisation of effectively immobile proteins is achieved. Here, we investigate the localisation to the division site of the slowly diffusing lipoprotein Pal, which anchors the outer membrane to the cell wall of Gram-negative bacteria. While the proton motive force-linked TolQRAB system is known to be required for this repositioning, the underlying mechanism is unresolved, especially given the very low mobility of Pal. We present a quantitative, mathematical model for Pal relocalisation in which dissociation of TolB-Pal complexes, powered by the proton motive force across the inner membrane, leads to the net transport of Pal along the outer membrane and its deposition at the division septum. We fit the model to experimental measurements of protein mobility and successfully test its predictions experimentally against mutant phenotypes. Our model not only explains a key aspect of cell division in Gram-negative bacteria, but also presents a physical mechanism for the transport of low-mobility proteins that may be applicable to multi-membrane organelles, such as mitochondria and chloroplasts.

In order for bacteria to successfully survive it is vital that they are able to concentrate proteins at precise sub-cellular locations. This usually occurs by capturing a freely diffusive protein, with the requirement that the protein is sufficiently mobile to reach the target location within the appropriate timeframe. Currently, it is unclear how immobile proteins are localised. Here, we examine how a very slowly diffusing protein Pal is able to localise to the centre of dividing cells to fulfil its role in membrane constriction. We present a mathematical model in which Pal is made mobile through the binding of another protein, TolB, and then deposited at the septum via active dissociation. This method is similar to a conveyor belt where Pal is collected everywhere but only deposited at the centre of the cell. We are able to fit this model to measurements of protein mobility and also test its predictions against mutant phenotypes. Our study is not only able to explain how this key protein is relocalised but also presents a general mechanism for the transport of slowly diffusing proteins that may be applicable to other systems.

Emergent Myxobacterial Behaviors Arise from Reversal Suppression Induced by Kin Contacts

A wide range of biological systems, from microbial swarms to bird flocks, display emergent behaviors driven by coordinated movement of individuals. To this end, individual organisms interact by recognizing their kin and adjusting their motility based on others around them. However, even in the best-studied systems, the mechanistic basis of the interplay between kin recognition and motility coordination is not understood. Here, using a combination of experiments and mathematical modeling, we uncover the mechanism of an emergent social behavior in Myxococcus xanthus. By overexpressing the cell surface adhesins TraA and TraB, which are involved in kin recognition, large numbers of cells adhere to one another and form organized macroscopic circular aggregates that spin clockwise or counterclockwise. Mechanistically, TraAB adhesion results in sustained cell-cell contacts that trigger cells to suppress cell reversals, and circular aggregates form as the result of cells’ ability to follow their own cellular slime trails. Furthermore, our in silico simulations demonstrate a remarkable ability to predict self-organization patterns when phenotypically distinct strains are mixed. For example, defying naive expectations, both models and experiments found that strains engineered to overexpress different and incompatible TraAB adhesins nevertheless form mixed circular aggregates. Therefore, this work provides key mechanistic insights into M. xanthus social interactions and demonstrates how local cell contacts induce emergent collective behaviors by millions of cells.

IMPORTANCE In many species, large populations exhibit emergent behaviors whereby all related individuals move in unison. For example, fish in schools can all dart in one direction simultaneously to avoid a predator. Currently, it is impossible to explain how such animals recognize kin through brain cognition and elicit such behaviors at a molecular level. However, microbes also recognize kin and exhibit emergent collective behaviors that are experimentally tractable. Here, using a model social bacterium, we engineer dispersed individuals to organize into synchronized collectives that create emergent patterns. With experimental and mathematical approaches, we explain how this occurs at both molecular and population levels. The results demonstrate how the combination of local physical interactions triggers intracellular signaling, which in turn leads to emergent behaviors on a population scale.

Myxococcus xanthus as a Model Organism for Peptidoglycan Assembly and Bacterial Morphogenesis

A fundamental question in biology is how cell shapes are genetically encoded and enzymatically generated. Prevalent shapes among walled bacteria include spheres and rods. These shapes are chiefly determined by the peptidoglycan (PG) cell wall. Bacterial division results in two daughter cells, whose shapes are predetermined by the mother. This makes it difficult to explore the origin of cell shapes in healthy bacteria. In this review, we argue that the Gram-negative bacterium Myxococcus xanthus is an ideal model for understanding PG assembly and bacterial morphogenesis, because it forms rods and spheres at different life stages. Rod-shaped vegetative cells of M. xanthus can thoroughly degrade their PG and form spherical spores. As these spores germinate, cells rebuild their PG and reestablish rod shape without preexisting templates. Such a unique sphere-to-rod transition provides a rare opportunity to visualize de novo PG assembly and rod-like morphogenesis in a well-established model organism.

Chiral-filament self-assembly on curved manifolds

Rod-like and banana-shaped proteins, like BAR-domain proteins and MreB proteins, adsorb on membranes and regulate the membrane curvature. The formation of large filamentous complexes of these proteins plays an important role in cellular processes like membrane trafficking, cytokinesis and cell motion. We propose a simplified model to investigate such curvature-dependent self-assembly processes. Anisotropic building blocks, modeled as trimer molecules, which have a preferred binding site, interact via pair-wise Lennard-Jones potentials. When several trimers assemble, they form an elastic ribbon with an intrinsic curvature and twist, controlled by bending and torsional rigidity. For trimer self-assembly on the curved surface of a cylindrical membrane, this leads to a preferred spatial orientation of the ribbon. We show that these interactions can lead to the formation of helices with several windings around the cylinder. The emerging helix angle and pitch depend on the rigidities and the intrinsic curvature and twist values. In particular, a well-defined and controllable helix angle emerges in the case of equal bending and torsional rigidity. The dynamics of filament growth is characterized by three regimes, in which filament length increases with the power laws t^z in time, with z≃ 3/4, z = 1/2, and z = 0 for short, intermediate, and long times, respectively. A comparison with the solutions of the Smoluchowski aggregation equation allows the identification of the underlying mechanism in the short-time regime as a crossover from size-independent to diffusion-limited aggregation. Thus, helical structures, as often observed in biology, can arise by self-assembly of anisotropic and chiral proteins.

Effects of spatial heterogeneity on bacterial genetic circuits

Intracellular spatial heterogeneity is frequently observed in bacteria, where the chromosome occupies part of the cell's volume and a circuit's DNA often localizes within the cell. How this heterogeneity affects core processes and genetic circuits is still poorly understood. In fact, commonly used ordinary differential equation (ODE) models of genetic circuits assume a well-mixed ensemble of molecules and, as such, do not capture spatial aspects. Reaction-diffusion partial differential equation (PDE) models have been only occasionally used since they are difficult to integrate and do not provide mechanistic understanding of the effects of spatial heterogeneity. In this paper, we derive a reduced ODE model that captures spatial effects, yet has the same dimension as commonly used well-mixed models. In particular, the only difference with respect to a well-mixed ODE model is that the association rate constant of binding reactions is multiplied by a coefficient, which we refer to as the binding correction factor (BCF). The BCF depends on the size of interacting molecules and on their location when fixed in space and it is equal to unity in a well-mixed ODE model. The BCF can be used to investigate how spatial heterogeneity affects the behavior of core processes and genetic circuits. Specifically, our reduced model indicates that transcription and its regulation are more effective for genes located at the cell poles than for genes located on the chromosome. The extent of these effects depends on the value of the BCF, which we found to be close to unity. For translation, the value of the BCF is always greater than unity, it increases with mRNA size, and, with biologically relevant parameters, is substantially larger than unity. Our model has broad validity, has the same dimension as a well-mixed model, yet it incorporates spatial heterogeneity. This simple-to-use model can be used to both analyze and design genetic circuits while accounting for spatial intracellular effects.

A Skeptic's Guide to Bacterial Mechanosensing

Surface sensing in bacteria is a precursor to the colonization of biotic and abiotic surfaces, and an important cause of drug resistance and virulence. As a motile bacterium approaches and adheres to a surface from the bulk fluid, the mechanical forces that act on it change. Bacteria are able to sense these changes in the mechanical load through a process termed mechanosensing. Bacterial mechanosensing has featured prominently in recent literature as playing a key role in surface sensing. However, the changes in mechanical loads on different parts of the cell at a surface vary in magnitudes as well as in signs. This confounds the determination of a causal relationship between the activation of specific mechanosensors and surface sensing. Here, we explain how contrasting mechanical stimuli arise on a surface-adherent cell and how known mechanosensors respond to these stimuli. The evidence for mechanosensing in select bacterial species is reinterpreted, with a focus on mechanosensitive molecular motors. We conclude with proposed criteria that bacterial mechanosensors must satisfy to successfully mediate surface sensing.

Physical Views on ParABS-Mediated DNA Segregation

In this chapter, we will focus on ParABS: an apparently simple, three-component system, required for the segregation of bacterial chromosomes and plasmids. We will specifically describe how biophysical measurements combined with physical modeling advanced our understanding of the mechanism of ParABS-mediated complex assembly, segregation and positioning.

Quantitative imaging of gene-expressing liposomes reveals rare favorable phenotypes

The biosynthesis of proteins from genomic DNA is a universal process in every living organism. Building a synthetic cell using separate biological parts hence implies to reconstitute a minimal gene expression apparatus and to compartmentalize it in a cell-mimicking environment. Previous studies have demonstrated that the PURE (Protein synthesis Using Recombinant Elements) system could be functionally encapsulated inside lipid vesicles. However, quantitative insights on functional consequences of spatial confinement of PURE system reactions remain scarce, which has hampered the full exploitation of gene-expressing liposomes as the fundamental unit to build an artificial cell. We report on direct imaging of tens of thousands of gene-expressing liposomes per sample allowing us to assess sub-population features in a statistically relevant manner. Both the vesicle size (diameter  <10 μm) and lipid composition (mixture of phospholipids with zwitterionic and negatively charged headgroups, including cardiolipin) are compatible with the properties of bacterial cells. Therefore, our liposomes provide a suitable chassis to host the Escherichia coli-derived PURE translation machinery and other bacterial processes in future developments. The potential of high-content imaging to identify rare phenotypes is demonstrated by the fact that a subset of the liposome population exhibits a remarkably high yield of synthesized protein or a prolonged expression lifespan that surpasses the performance of ensemble liposome-averaged and bulk reactions. Among the three commercial PURE systems tested, PUREfrex2.0 offers the most favorable phenotypes displaying both high yield and long protein synthesis lifespan. Moreover, probing membrane permeability reveals a large heterogeneity amongst liposomes. In situ expression and membrane embedding of the pore-forming connexin leads to a characteristic permeability time profile, while increasing the fraction of permeable liposomes in the population. We see diversity in gene expression dynamics and membrane permeability as an opportunity to complement a rational design approach aiming at further implementing biological functions in liposome-based synthetic cells.

Cell geometry and leaflet bilayer asymmetry regulate domain formation in plasma membranes

We model how pattern formation in a multicomponent lipid bilayer pinned to an elastic substrate is governed by the interplay between lipid phase separation and the tendency of domains of high intrinsic curvature lipids to deform the membrane away from a stiff substrate such as the cell wall. The emergent patterns, which include compact and striped lipid microdomains, are anticorrelated across the two leaflets and depend on leaflet asymmetry, the ability of lipids to flip between leaflets, and the global geometry. We characterize analytically the dependence of stripe width on lipid parameters, and consider the implications of interleaflet patterning for curvature-dependent lipid localization.

Protein polarization driven by nucleoid exclusion of DnaK(HSP70)-substrate complexes

Many bacterial proteins require specific subcellular localization for function. How Escherichia coli proteins localize at one pole, however, is still not understood. Here, we show that the DnaK (HSP70) chaperone controls unipolar localization of the Shigella IpaC type III secretion substrate. While preventing the formation of lethal IpaC aggregates, DnaK promoted the incorporation of IpaC into large and dynamic complexes (LDCs) restricted at the bacterial pole through nucleoid occlusion. Unlike stable polymers and aggregates, LDCs show dynamic behavior indicating that nucleoid occlusion also applies to complexes formed through transient interactions. Fluorescence recovery after photobleaching analysis shows DnaK-IpaC exchanges between opposite poles and DnaKJE-mediated incorporation of immature substrates in LDCs. These findings reveal a key role for LDCs as reservoirs of functional DnaK-substrates that can be rapidly mobilized for secretion triggered upon bacterial contact with host cells.

Many bacterial proteins exhibit spatially defined localization important for function. Here the authors show that the polar localization of Shigella IpaC type III secretion substrate is mediated by its interaction with the DnaK chaperone and occlusion by the bacterial nucleoid.

The Enigmatic Genome of an Obligate Ancient Spiroplasma Symbiont in a Hadal Holothurian

Protective symbiosis has been reported in many organisms, but the molecular mechanisms of the mutualistic interactions between the symbionts and their hosts are unclear. Here, we sequenced the 424-kbp genome of "Candidatus Spiroplasma holothuricola," which dominated the hindgut microbiome of a sea cucumber, a major scavenger captured in the Mariana Trench (6,140 m depth). Phylogenetic relationships indicated that the dominant bacterium in the hindgut was derived from a basal group of Spiroplasma species. In this organism, the genes responsible for the biosynthesis of amino acids, glycolysis, and sugar transporters were lost, strongly suggesting endosymbiosis. The highly decayed genome consists of two chromosomes and harbors genes coding for proteolysis, microbial toxin, restriction-methylation systems, and clustered regularly interspaced short palindromic repeats (CRISPRs), composed of three cas genes and 76 CRISPR spacers. The holothurian host is probably protected against invading viruses from sediments by the CRISPRs/Cas and restriction systems of the endosymbiotic spiroplasma. The protective endosymbiosis indicates the important ecological role of the ancient Spiroplasma symbiont in the maintenance of hadal ecosystems.IMPORTANCE Sea cucumbers are major inhabitants in hadal trenches. They collect microbes in surface sediment and remain tolerant against potential pathogenic bacteria and viruses. This study presents the genome of endosymbiotic spiroplasmas in the gut of a sea cucumber captured in the Mariana Trench. The extreme reduction of the genome and loss of essential metabolic pathways strongly support its endosymbiotic lifestyle. Moreover, a considerable part of the genome was occupied by a CRISPR/Cas system to provide immunity against viruses and antimicrobial toxin-encoding genes for the degradation of microbes. This novel species of Spiroplasma is probably an important protective symbiont for the sea cucumbers in the hadal zone.

Transmembrane region of bacterial chemoreceptor is capable of promoting protein clustering

Many membrane proteins are known to form higher-order oligomers, but the degree to which membrane regions could facilitate protein complex assembly remains largely unclear. Clusters of chemotaxis receptors are among the most prominent structures in the bacterial cell membrane, and they play important functions in processing of chemotactic signals. Although much work has been done to elucidate mechanisms of cluster formation, it almost exclusively focused on cytoplasmic interactions among receptors and other chemotaxis proteins, whereas involvement of membrane-mediated interactions was only hypothesized. Here we used imaging of constructs composed of only a fluorescent protein and the TM helices of Tar to demonstrate that interactions between the lipid bilayer and transmembrane (TM) helices of Escherichia coli chemoreceptors alone are sufficient to mediate clustering. We found that the ability to cluster depends on the sequence or length of the TM helices, implying that certain conformations of these helices facilitate clustering, whereas others do not. Notably, observed sequence specificity was apparently consistent with differences in clustering between native E. coli receptors, with the TM sequence of better-clustering high-abundance receptors being more efficient in promoting membrane-mediated complex formation. These results indicate that being more than just membrane anchors, TM helices could play an important role in the clustering and organization of membrane proteins in bacteria.

Physical model of protein cluster positioning in growing bacteria

Chemotaxic receptors in bacteria form clusters at cell poles and also laterally, and this clustering plays an important role in signal transduction. These clusters were found to be periodically arranged on the surface of the bacterium Escherichia coli, independent of any known positioning mechanism. In this work we extend a model based on diffusion and aggregation to more realistic geometries and present a means based on "bursty" protein production to distinguish spontaneous positioning from an independently existing positioning mechanism. We also consider the case of isotropic cellular growth and characterize the degree of order arising spontaneously. Our model could also be relevant for other examples of periodically positioned protein clusters in bacteria.

Nucleoid and cytoplasmic localization of small RNAs in Escherichia coli

Bacterial small RNAs (sRNAs) regulate protein production by binding to mRNAs and altering their translation and degradation. sRNAs are smaller than most mRNAs but larger than many proteins. Therefore it is uncertain whether sRNAs can enter the nucleoid to target nascent mRNAs. Here, we investigate the intracellular localization of sRNAs transcribed from plasmids in Escherichia coli using RNA fluorescent in-situ hybridization. We found that sRNAs (GlmZ, OxyS, RyhB and SgrS) have equal preference for the nucleoid and cytoplasm, and no preferential localization at the cell membrane. We show using the gfp mRNA (encoding green fluorescent protein) that non-sRNAs can be engineered to have different proportions of nucleoid and cytoplasmic localization by altering their length and/or translation. The same localization as sRNAs was achieved by decreasing gfp mRNA length and translation, which suggests that sRNAs and other RNAs may enter the densely packed DNA of the nucleoid if they are sufficiently small. We also found that the Hfq protein, which binds sRNAs, minimally affects sRNA localization. Important implications of our findings for engineering synthetic circuits are: (i) sRNAs can potentially bind nascent mRNAs in the nucleoid, and (ii) localization patterns and distribution volumes of sRNAs can differ from some larger RNAs.

The Slow Mobility of the ParA Partitioning Protein Underlies Its Steady-State Patterning in Caulobacter

In bacteria, ParABS systems mediate intracellular transport of various cargos, including chromosomal regions in Caulobacter crescentus. Transport of the ParB/parS partition complex requires the DNA-binding activity of ParA, which transiently tethers the partition complex during translocation. In C. crescentus, the directionality of the transport is set up by a gradient of ParA whose concentration gradually increases from one end of the cell (old pole) to the other (new pole). Importantly, this ParA gradient is already observed before DNA replication and segregation are initiated when the partition complex is anchored at the old pole. How such micron-scale ParA pattern is established and maintained before the initiation of chromosome segregation has not been experimentally established. Although the stimulation of ParA ATPase activity by the localized ParB/parS partition complex is thought to be involved, this activity alone cannot quantitatively describe the ParA pattern observed inside cells. Instead, our experimental and theoretical study shows that the missing key component for achieving the experimentally observed steady-state ParA patterning is the slow mobility of ParA dimers (D ∼10(-3)μm(2)/s) due to intermittent DNA binding. Our model recapitulates the entire steady-state ParA distribution observed experimentally, including the shape of the gradient as well as ParA accumulation at the location of the partition complex. Stochastic simulations suggest that cell-to-cell variability in ParA pattern is due to the low ParA copy number in C. crescentus cells. The model also accounts for an apparent exclusion of ParA from regions with small spacing between partition complexes observed in filamentous cells. Collectively, our work demonstrates that in addition to its function in mediating transport, the conserved DNA-binding property of ParA has a critical function before DNA segregation by setting up a ParA pattern required for transport directionality.

DNA-relay mechanism is sufficient to explain ParA-dependent intracellular transport and patterning of single and multiple cargos

Spatial ordering of macromolecular components inside cells is important for cellular physiology and replication. In bacteria, ParA/B systems are known to generate various intracellular patterns that underlie the transport and partitioning of low-copy-number cargos such as plasmids. ParA/B systems consist of ParA, an ATPase that dimerizes and binds DNA upon ATP binding, and ParB, a protein that binds the cargo and stimulates ParA ATPase activity. Inside cells, ParA is asymmetrically distributed, forming a propagating wave that is followed by the ParB-rich cargo. These correlated dynamics lead to cargo oscillation or equidistant spacing over the nucleoid depending on whether the cargo is in single or multiple copies. Currently, there is no model that explains how these different spatial patterns arise and relate to each other. Here, we test a simple DNA-relay model that has no imposed asymmetry and that only considers the ParA/ParB biochemistry and the known fluctuating and elastic dynamics of chromosomal loci. Stochastic simulations with experimentally derived parameters demonstrate that this model is sufficient to reproduce the signature patterns of ParA/B systems: the propagating ParA gradient correlated with the cargo dynamics, the single-cargo oscillatory motion, and the multicargo equidistant patterning. Stochasticity of ATP hydrolysis breaks the initial symmetry in ParA distribution, resulting in imbalance of elastic force acting on the cargo. Our results may apply beyond ParA/B systems as they reveal how a minimal system of two players, one binding to DNA and the other modulating this binding, can transform directionally random DNA fluctuations into directed motion and intracellular patterning.