Organization of the bacterial nucleoid by DNA-bridging proteins and globular crowders.

The genomic DNA of bacteria occupies only a fraction of the cell called the nucleoid, although it is not bounded by any membrane and would occupy a volume hundreds of times larger than the cell in the absence of constraints. The two most important contributions to the compaction of the DNA coil are the cross-linking of the DNA by nucleoid proteins (like H-NS and StpA) and the demixing of DNA and other abundant globular macromolecules which do not bind to the DNA (like ribosomes). The present work deals with the interplay of DNA-bridging proteins and globular macromolecular crowders, with the goal of determining the extent to which they collaborate in organizing the nucleoid. In order to answer this question, a coarse-grained model was developed and its properties were investigated through Brownian dynamics simulations. These simulations reveal that the radius of gyration of the DNA coil decreases linearly with the effective volume ratio of globular crowders and the number of DNA bridges formed by nucleoid proteins in the whole range of physiological values. Moreover, simulations highlight the fact that the number of DNA bridges formed by nucleoid proteins depends crucially on their ability to self-associate (oligomerize). An explanation for this result is proposed in terms of the mean distance between DNA segments and the capacity of proteins to maintain DNA-bridging in spite of the thermal fluctuations of the DNA network. Finally, simulations indicate that non-associating proteins preserve a high mobility inside the nucleoid while contributing to its compaction, leading to a DNA/protein complex which looks like a liquid droplet. In contrast, self-associating proteins form a little deformable network which cross-links the DNA chain, with the consequence that the DNA/protein complex looks more like a gel.

Impact of Self-Association on the Architectural Properties of Bacterial Nucleoid Proteins.

The chromosomal DNA of bacteria is folded into a compact body called the nucleoid, which is composed essentially of DNA (∼80%), RNA (∼10%), and a number of different proteins (∼10%). These nucleoid proteins act as regulators of gene expression and influence the organization of the nucleoid by bridging, bending, or wrapping the DNA. These so-called architectural properties of nucleoid proteins are still poorly understood. For example, the reason why certain proteins compact the DNA coil in certain environments but make the DNA more rigid instead in other environments is the subject of ongoing debates. Here, we address the question of the impact of the self-association of nucleoid proteins on their architectural properties and try to determine whether differences in self-association are sufficient to induce large changes in the organization of the DNA coil. More specifically, we developed two coarse-grained models of proteins, which interact identically with the DNA but self-associate differently by forming either clusters or filaments in the absence of the DNA. We showed through Brownian dynamics simulations that self-association of the proteins dramatically increases their ability to shape the DNA coil. Moreover, we observed that cluster-forming proteins significantly compact the DNA coil (similar to the DNA-bridging mode of H-NS proteins), whereas filament-forming proteins significantly increase the stiffness of the DNA chain instead (similar to the DNA-stiffening mode of H-NS proteins). This work consequently suggests that the knowledge of the DNA-binding properties of the proteins is in itself not sufficient to understand their architectural properties. Rather, their self-association properties must also be investigated in detail because they might actually drive the formation of different DNA-protein complexes.

Requirements for DNA-Bridging Proteins to Act as Topological Barriers of the Bacterial Genome.

Bacterial genomes have been shown to be partitioned into several-kilobase-long chromosomal domains that are topologically independent from each other, meaning that change of DNA superhelicity in one domain does not propagate to neighbors. Both in vivo and in vitro experiments have been performed to question the nature of the topological barriers at play, leading to several predictions on possible molecular actors. Here, we address the question of topological barriers using polymer models of supercoiled DNA chains that are constrained such as to mimic the action of predicted molecular actors. More specifically, we determine under which conditions DNA-bridging proteins may act as topological barriers. To this end, we developed a coarse-grained bead-and-spring model and investigated its properties through Brownian dynamics simulations. As a result, we find that DNA-bridging proteins must exert rather strong constraints on their binding sites; they must block the diffusion of the excess of twist through the two binding sites on the DNA molecule and, simultaneously, prevent the rotation of one DNA segment relative to the other one. Importantly, not all DNA-bridging proteins satisfy this second condition. For example, single bridges formed by proteins that bind DNA nonspecifically, like H-NS dimers, are expected to fail with this respect. Our findings might also explain, in the case of specific DNA-bridging proteins like LacI, why multiple bridges are required to create stable independent topological domains. Strikingly, when the relative rotation of the DNA segments is not prevented, relaxation results in complex intrication of the two domains. Moreover, although the value of the torsional stress in each domain may vary, their differential is preserved. Our work also predicts that nucleoid-associated proteins known to wrap DNA must form higher protein-DNA complexes to efficiently work as topological barriers.

Bacterial Nucleoid: Interplay of DNA Demixing and Supercoiling.

This work addresses the question of the interplay of DNA demixing and supercoiling in bacterial cells. Demixing of DNA from other globular macromolecules results from the overall repulsion between all components of the system and leads to the formation of the nucleoid, which is the region of the cell that contains the genomic DNA in a rather compact form. Supercoiling describes the coiling of the axis of the DNA double helix to accommodate the torsional stress injected in the molecule by topoisomerases. Supercoiling is able to induce some compaction of the bacterial DNA, although to a lesser extent than demixing. In this work, we investigate the interplay of these two mechanisms with the goal of determining whether the total compaction ratio of the DNA is the mere sum or some more complex function of the compaction ratios due to each mechanism. To this end, we developed a coarse-grained bead-and-spring model and investigated its properties through Brownian dynamics simulations. This work reveals that there actually exist different regimes, depending on the crowder volume ratio and the DNA superhelical density. In particular, a regime in which the effects of DNA demixing and supercoiling on the compaction of the DNA coil simply add up is shown to exist up to moderate values of the superhelical density. In contrast, the mean radius of the DNA coil no longer decreases above this threshold and may even increase again for sufficiently large crowder concentrations. Finally, the model predicts that the DNA coil may depart from the spherical geometry very close to the jamming threshold as a trade-off between the need to minimize both the bending energy of the stiff plectonemes and the volume of the DNA coil to accommodate demixing.

Preferential Localization of the Bacterial Nucleoid.

Prokaryotes do not make use of a nucleus membrane to segregate their genetic material from the cytoplasm, so that their nucleoid is potentially free to explore the whole volume of the cell. Nonetheless, high resolution images of bacteria with very compact nucleoids show that such spherical nucleoids are invariably positioned at the center of mononucleoid cells. The present work aims to determine whether such preferential localization results from generic (entropic) interactions between the nucleoid and the cell membrane or instead requires some specific mechanism, like the tethering of DNA at mid-cell or periodic fluctuations of the concentration gradient of given chemical species. To this end, we performed numerical simulations using a coarse-grained model based on the assumption that the formation of the nucleoid results from a segregative phase separation mechanism driven by the de-mixing of the DNA and non-binding globular macromolecules. These simulations show that the abrupt compaction of the DNA coil, which takes place at large crowder density, close to the jamming threshold, is accompanied by the re-localization of the DNA coil close to the regions of the bounding wall with the largest curvature, like the hemispherical caps of rod-like cells, as if the DNA coil were suddenly acquiring the localization properties of a solid sphere. This work therefore supports the hypothesis that the localization of compact nucleoids at regular cell positions involves either some anchoring of the DNA to the cell membrane or some dynamical localization mechanism.

A segregative phase separation scenario of the formation of the bacterial nucleoid.

The mechanism responsible for the compaction of the genomic DNA of bacteria inside a structure called the nucleoid is a longstanding but still lively debated question. Most puzzling is the fact that the nucleoid occupies only a small fraction of the cell, although it is not separated from the rest of the cytoplasm by any membrane and would occupy a volume about a thousand times larger outside the cell. Here, by performing numerical simulations using coarse-grained models, we elaborate on the conjecture that the formation of the nucleoid may result from a segregative phase separation mechanism driven by the demixing of the DNA coil and non-binding globular macromolecules present in the cytoplasm, presumably functional ribosomes. Simulations performed with crowders having a spherical, dumbbell or octahedral geometry highlight the sensitive dependence of the level of DNA compaction on the dissymmetry of DNA/DNA, DNA/crowder, and crowder/crowder repulsive interactions, thereby supporting the segregative phase separation scenario. Simulations also consistently predict a much stronger DNA compaction close to the jamming threshold. Moreover, simulations performed with crowders of different sizes suggest that the final density distribution of each species results from the competition between thermodynamic forces and steric hindrance, so that bigger crowders are expelled selectively from the nucleoid only at moderate total crowder concentrations. This work leads to several predictions, which may eventually be tested experimentally.

Role of Salt Valency in the Switch of H-NS Proteins between DNA-Bridging and DNA-Stiffening Modes.

This work investigates the interactions of H-NS proteins and bacterial genomic DNA through computer simulations performed with a coarse-grained model. The model was developed specifically to study the switch of H-NS proteins from the DNA-stiffening to the DNA-bridging mode, which has been observed repeatedly upon addition of multivalent cations to the buffer but is still not understood. Unraveling the corresponding mechanism is all the more crucial, as the regulation properties of H-NS proteins, as well as other nucleoid proteins, are linked to their DNA-binding properties. The simulations reported here support a mechanism, according to which the primary role of multivalent cations consists in decreasing the strength of H-NS/DNA interactions compared to H-NS/H-NS interactions, with the latter ones becoming energetically favored with respect to the former ones above a certain threshold of the effective valency of the cations of the buffer. Below the threshold, H-NS dimers form filaments, which stretch along the DNA molecule but are quite inefficient in bridging genomically distant DNA sites (DNA-stiffening mode). In contrast, just above the threshold, H-NS dimers form three-dimensional clusters, which are able to connect DNA sites that are distant from the genomic point of view (DNA-bridging mode). The model provides clear rationales for the experimental observations that the switch between the two modes is a threshold effect and that the ability of H-NS dimers to form higher order oligomers is crucial for their bridging capabilities.

Recovery of mechanical pressure in a gas of underdamped active dumbbells with Brownian noise.

In contrast with a gas at thermodynamic equilibrium, the mean force exerted on a wall by a gas of active particles usually depends on the confining potential, thereby preventing a proper definition of mechanical pressure. In this paper, we investigate numerically the properties of a gas of underdamped self-propelled dumbbells subject to Brownian noise of increasing intensity, in order to understand how the notion of pressure is recovered as noise progressively masks the effects of self-propulsion and the system approaches thermodynamic equilibrium. The simulations performed for a mobile asymmetric wall separating two chambers containing an equal number of active dumbbells highlight some subtle and unexpected properties of the system. First, Brownian noise of moderate intensity is sufficient to let mean forces equilibrate for small values of the damping coefficient, while much stronger noise is required for larger values of the damping coefficient. Moreover, the displacement of the mean position of the wall upon increase of the intensity of the noise is not necessarily monotonous and may instead display changes of direction. Both facts actually reflect the existence of several mechanisms leading to the rupture of force balance, which tend to displace the mean position of the wall towards different directions and display different robustness against an increase of the intensity of Brownian noise. This work therefore provides a clear illustration of the fact that driving an autonomous system towards (or away from) thermodynamic equilibrium may not be a straightforward process, but may instead proceed through the variations of the relative weights of several conflicting mechanisms.

Coarse-Grained Model of the Demixing of DNA and Non-Binding Globular Macromolecules.

The volume occupied by the unconstrained genomic DNA of prokaryotes in saline solutions is thousand times larger than the cell. Moreover, it is not separated from the rest of the cell by a membrane. Nevertheless, it occupies only a small fraction of the cell called the nucleoid. The mechanisms leading to such compaction are the matter of ongoing debates. The present work aims at exploring a newly proposed mechanism, according to which the formation of the nucleoid would result from the demixing of the DNA and nonbinding globular macromolecules of the cytoplasm, like ribosomes. To this end, a coarse-grained model of prokaryotic cells was developed, and demixing was analyzed as a function of the size and number of crowders. The model suggests that compaction of the DNA is actually governed by the volume occupancy ratio of the crowders and remains weak almost up to the jamming critical density. Strong compaction is however observed just before jamming, suggesting that crowding and electrostatic repulsion work synergetically in this limit. Finally, simulations performed with crowders with different sizes indicate that the DNA and the largest crowders demix preferentially. Together with the recent observation of the gradual compaction of long DNA molecules upon increase of the concentration of bovine serum albumin proteins and silica nanoparticles, this work supports the demixing mechanism as a key player for the formation of the nucleoid.

Pressure of a gas of underdamped active dumbbells.

The pressure exerted on a wall by a gas at equilibrium does not depend on the shape of the confining potential defining the walls. In contrast, it has been shown recently [A. P. Solon et al., Nat. Phys. 11, 673 (2015)] that a gas of overdamped active particles exerts on a wall a force that depends on the confining potential, resulting in a net force on an asymmetric wall between two chambers at equal densities. Here, considering a model of underdamped self-propelled dumbbells in two dimensions, we study how the behavior of the pressure depends on the damping coefficient of the dumbbells, thus exploring inertial effects. We find in particular that the force exerted on a moving wall between two chambers at equal density continuously vanishes at low damping coefficient, and exhibits a complex dependence on the damping coefficient at low density, when collisions are scarce. We further show that this behavior of the pressure can to a significant extent be understood in terms of the trajectories of individual particles close to and in contact with the wall.

Compaction of bacterial genomic DNA: clarifying the concepts.

The unconstrained genomic DNA of bacteria forms a coil, whose volume exceeds 1000 times the volume of the cell. Since prokaryotes lack a membrane-bound nucleus, in sharp contrast with eukaryotes, the DNA may consequently be expected to occupy the whole available volume when constrained to fit in the cell. Still, it has been known for more than half a century that the DNA is localized in a well-defined region of the cell, called the nucleoid, which occupies only 15% to 25% of the total volume. Although this problem has focused the attention of many scientists in recent decades, there is still no certainty concerning the mechanism that enables such a dramatic compaction. The goal of this Topical Review is to take stock of our knowledge on this question by listing all possible compaction mechanisms with the proclaimed desire to clarify the physical principles they are based upon and discuss them in the light of experimental results and the results of simulations based on coarse-grained models. In particular, the fundamental differences between ψ-condensation and segregative phase separation and between the condensation by small and long polycations are highlighted. This review suggests that the importance of certain mechanisms, like supercoiling and the architectural properties of DNA-bridging and DNA-bending nucleoid proteins, may have been overestimated, whereas other mechanisms, like segregative phase separation and the self-association of nucleoid proteins, as well as the possible role of the synergy of two or more mechanisms, may conversely deserve more attention.

Free-energy landscape and characteristic forces for the initiation of DNA unzipping.

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.

Equilibration of complexes of DNA and H-NS proteins on charged surfaces: a coarse-grained model point of view.

The Histone-like Nucleoid Structuring protein (H-NS) is a nucleoid-associated protein, which is involved in both gene regulation and DNA compaction. Although it is a key player in genome organization by forming bridges between DNA duplexes, the precise structure of complexes of DNA and H-NS proteins is still not well understood. In particular, it is not clear whether the structure of DNA/H-NS complexes in the living cell is similar to that of complexes deposited on mica surfaces, which may be observed by AFM microscopy. A coarse-grained model, which helps getting more insight into this question, is described and analyzed in the present paper. This model is able of describing both the bridging of bacterial DNA by H-NS in the bulk and the deposition and equilibration of the complex on a charged surface. Simulations performed with the model reveal that a slight attraction between DNA and the charged surface is sufficient to let DNA/H-NS complexes reorganize from 3D coils to planar plasmids bridged by H-NS proteins similar to those observed by AFM microscopy. They furthermore highlight the antagonistic effects of the interactions between DNA and the surface. Indeed, increasing these interactions slows down the equilibration of naked plasmids on the surface but, on the other hand, enables a faster equilibration of DNA/H-NS complexes. Based on the distribution of the lifetimes of H-NS bridges and the time evolution of the number of trans-binding protein dimers during equilibration of the complexes on the surface, it is argued that the decrease of the equilibration time of the complex upon increase of the interaction strength between DNA and the surface is ascribable to the associated decrease of the probability to form new bridges between DNA and the proteins.

Elastic models of the fast traps of carnivorous Dionaea and Aldrovanda.

The carnivorous aquatic waterwheel plant (Aldrovanda vesiculosa L.) and the closely related terrestrial Venus flytrap (Dionaea muscipula Sol. ex J. Ellis) both feature elaborate snap-traps, which shut after reception of an external mechanical stimulus by prey animals. Although Aldrovanda is usually considered a miniature aquatic Dionaea, the shutting mechanisms of the two plants differ quite markedly. The fast shutting of Aldrovanda is indeed based on a simple swelling or shrinking mechanism, while the movement of Dionaea's traps is accelerated by the snap-buckling of the lobes. The purpose of this Brief Report is to describe several key improvements to the elastic models that have recently been introduced to elucidate these movements [Poppinga and Joyeux, Phys. Rev. E 84, 041928 (2011)]. In particular, a precise mechanism for the action of the motor cells of Aldrovanda is proposed, the facts that the opening of the leaves of Dionaea is an irreversible mechanism based on growth and that the strain field is anisotropic and much smaller on the inner than on the outer surface of the leaves during shutting are taken properly into account, and a more accurate formula for calculating mean curvatures is used. The improvements brought to the model are described in detail and the physical consequences of these improvements are discussed.

A model of H-NS mediated compaction of bacterial DNA.

The histone-like nucleoid structuring protein (H-NS) is a nucleoid-associated protein, which is involved in both gene regulation and DNA compaction. H-NS can bind to DNA in two different ways: in trans, by binding to two separate DNA duplexes, or in cis, by binding to different sites on the same duplex. Based on scanning force microscopy imaging and optical trap-driven unzipping assays, it has recently been suggested that DNA compaction may result from the antagonistic effects of H-NS binding to DNA in trans and cis configurations. To get more insight into the compaction mechanism, we constructed a coarse-grained model description of the compaction of bacterial DNA by H-NS. These simulations highlight the fact that DNA compaction indeed results from the subtle equilibrium between several competing factors, which include the deformation dynamics of the plasmid and the several binding modes of protein dimers to DNA, i.e., dangling configurations, cis- and trans-binding. In particular, the degree of compaction is extremely sensitive to the difference in binding energies of the cis and trans configurations. Our simulations also point out that the conformations of the DNA-protein complexes are significantly different in bulk and in planar conditions, suggesting that conformations observed on mica surfaces may differ significantly from those that prevail in living cells.

Different mechanics of snap-trapping in the two closely related carnivorous plants Dionaea muscipula and Aldrovanda vesiculosa.

The carnivorous aquatic waterwheel plant (Aldrovanda vesiculosa L.) and the closely related terrestrial venus flytrap (Dionaea muscipula Sol. ex J. Ellis) both feature elaborate snap-traps, which shut after reception of an external mechanical stimulus by prey animals. Traditionally, Aldrovanda is considered as a miniature, aquatic Dionaea, an assumption which was already established by Charles Darwin. However, videos of snapping traps from both species suggest completely different closure mechanisms. Indeed, the well-described snapping mechanism in Dionaea comprises abrupt curvature inversion of the two trap lobes, while the closing movement in Aldrovanda involves deformation of the trap midrib but not of the lobes, which do not change curvature. In this paper, we present detailed mechanical models for these plants, which are based on the theory of thin solid membranes and explain this difference by showing that the fast snapping of Aldrovanda is due to kinematic amplification of the bending deformation of the midrib, while that of Dionaea unambiguously relies on the buckling instability that affects the two lobes.

Thermal and mechanical denaturation properties of a DNA model with three sites per nucleotide.

In this paper, we show that the coarse grain model for DNA, which has been proposed recently by Knotts et al. [J. Chem. Phys. 126, 084901 (2007)], can be adapted to describe the thermal and mechanical denaturation of long DNA sequences by adjusting slightly the base pairing contribution. The adjusted model leads to (i) critical temperatures for long homogeneous sequences that are in good agreement with both experimental ones and those obtained from statistical models, (ii) a realistic step-like denaturation behaviour for long inhomogeneous sequences, and (iii) critical forces at ambient temperature of the order of 10 pN, close to measured values. The adjusted model furthermore supports the conclusion that the thermal denaturation of long homogeneous sequences corresponds to a first-order phase transition and yields a critical exponent for the critical force equal to σ = 0.70. This model is both geometrically and energetically realistic, in the sense that the helical structure and the grooves, where most proteins bind, are satisfactorily reproduced, while the energy and the force required to break a base pair lie in the expected range. It therefore represents a promising tool for studying the dynamics of DNA-protein specific interactions at an unprecedented detail level.

Mechanical model of the ultrafast underwater trap of Utricularia.

The underwater traps of the carnivorous plants of the Utricularia species catch their prey through the repetition of an "active slow deflation followed by passive fast suction" sequence. In this paper, we propose a mechanical model that describes both phases and strongly supports the hypothesis that the trap door acts as a flexible valve that buckles under the combined effects of pressure forces and the mechanical stimulation of trigger hairs, and not as a panel articulated on hinges. This model combines two different approaches, namely (i) the description of thin membranes as triangle meshes with strain and curvature energy, and (ii) the molecular dynamics approach, which consists of computing the time evolution of the position of each vertex of the mesh according to Langevin equations. The only free parameter in the expression of the elastic energy is the Young's modulus E of the membranes. The values for this parameter are unequivocally obtained by requiring that the trap model fires, like real traps, when the pressure difference between the outside and the inside of the trap reaches about 15 kPa. Among other results, our simulations show that, for a pressure difference slightly larger than the critical one, the door buckles, slides on the threshold, and finally swings wide open, in excellent agreement with the sequence observed in high-speed videos.

Ultra-fast underwater suction traps.

Carnivorous aquatic Utricularia species catch small prey animals using millimetre-sized underwater suction traps, which have fascinated scientists since Darwin's early work on carnivorous plants. Suction takes place after mechanical triggering and is owing to a release of stored elastic energy in the trap body accompanied by a very fast opening and closing of a trapdoor, which otherwise closes the trap entrance watertight. The exceptional trapping speed--far above human visual perception--impeded profound investigations until now. Using high-speed video imaging and special microscopy techniques, we obtained fully time-resolved recordings of the door movement. We found that this unique trapping mechanism conducts suction in less than a millisecond and therefore ranks among the fastest plant movements known. Fluid acceleration reaches very high values, leaving little chance for prey animals to escape. We discovered that the door deformation is morphologically predetermined, and actually performs a buckling/unbuckling process, including a complete trapdoor curvature inversion. This process, which we predict using dynamical simulations and simple theoretical models, is highly reproducible: the traps are autonomously repetitive as they fire spontaneously after 5-20 h and reset actively to their ready-to-catch condition.

Comparison of kinetic and dynamical models of DNA-protein interaction and facilitated diffusion.

It has long been asserted that proteins such as transcription factors may locate their target in DNA sequences at rates that surpass by several orders of magnitude the three-dimensional diffusion limit thanks to facilitated diffusion, that is, the combination of one-dimensional (sliding along the DNA) and three-dimensional diffusion. This claim has been supported throughout the years by several mass action kinetic models, while the dynamical model we proposed recently (J. Chem. Phys. 2009, 130, 015103) suggests that acceleration of targeting due to facilitated diffusion cannot be large. In order to solve this apparent contradiction, we performed additional simulations to compare the results obtained with our model to those obtained with the kinetic model of Klenin et al. (Phys. Rev. Lett. 2006, 96, 018104). We show in this paper that the two models actually support each other and agree in predicting a low efficiency for facilitated diffusion. Extrapolation of these results to real systems even indicates that facilitated diffusion necessarily slows down the targeting process compared to three-dimensional diffusion.