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

Elucidation of Spatial Positioning of Ribosomes around Chromosome in Escherichia coli Cytoplasm via a Data-Informed Polymer-Based Model

The spatial arrangement of ribosomes and chromosome in Escherichia coli's cytoplasm challenges conventional wisdom. Contrary to the notion of ribosomes acting as inert crowders to the chromosome in the cytoplasm, here we propose a nuanced view by integrating a wide array of experimental data sets into a polymer-based computer model. A set of data-informed computer simulations determines that a delicate balance of attractive and repulsive interactions between ribosomes and the chromosome is required in order to reproduce experimentally obtained linear densities and brings forth the view that ribosomes are not mere inert crowders in the cytoplasm. The model finds that the ribosomes represent themselves as a poor solvent for the chromosome with a 50 nm mesh size, consistent with previous experimental analysis. Our multidimensional analysis of ribosome distribution, both free (30S and 50S) and bound (70S polysome), uncovers a relatively less pronounced segregation pattern than previously thought. Notably, we identify a ribosome-rich central region within the innermost core of the nucleoid. Moreover, our exploration of the chromosome mesh size and the conformation of bound ribosomes suggests that these ribosomes maintain elongated shapes, enabling them to navigate through the chromosome mesh and access the central core. This dynamic localization challenges the static segregation model and underscores the pivotal role of ribosome-chromosome interactions in cellular media.

Physical models of bacterial chromosomes

The interplay between bacterial chromosome organization and functions such as transcription and replication can be studied in increasing detail using novel experimental techniques. Interpreting the resulting quantitative data, however, can be theoretically challenging. In this minireview, we discuss how connecting experimental observations to biophysical theory and modeling can give rise to new insights on bacterial chromosome organization. We consider three flavors of models of increasing complexity: simple polymer models that explore how physical constraints, such as confinement or plectoneme branching, can affect bacterial chromosome organization; bottom-up mechanistic models that connect these constraints to their underlying causes, for instance, chromosome compaction to macromolecular crowding, or supercoiling to transcription; and finally, data-driven methods for inferring interpretable and quantitative models directly from complex experimental data. Using recent examples, we discuss how biophysical models can both deepen our understanding of how bacterial chromosomes are structured and give rise to novel predictions about bacterial chromosome organization.

Development of a Data-Driven Integrative Model of a Bacterial Chromosome

The chromosome of archetypal bacteria E. coli is known for a complex topology with a 4.6 × 106 base pairs (bp) long sequence of nucleotides packed within a micrometer-sized cellular confinement. The inherent organization underlying this chromosome eludes general consensus due to the lack of a high-resolution picture of its conformation. Here we present our development of an integrative model of E. coli at a 500 bp resolution (https://github.com/JMLab-tifrh/ecoli_finer), which optimally combines a set of multiresolution genome-wide experimentally measured data within a framework of polymer based architecture. In particular the model is informed with an intragenome contact probability map at 5000 bp resolution derived via the Hi-C experiment and RNA-sequencing data at 500 bp resolution. Via dynamical simulations, this data-driven polymer based model generates an appropriate conformational ensemble commensurate with chromosome architectures that E. coli adopts. As a key hallmark of the E. coli chromosome the model spontaneously self-organizes into a set of nonoverlapping macrodomains and suitably locates plectonemic loops near the cell membrane. As novel extensions, it predicts a contact probability map simulated at a higher resolution than precedent experiments and can demonstrate segregation of chromosomes in a partially replicating cell. Finally, the modular nature of the model helps us devise control simulations to quantify the individual role of key features in hierarchical organization of the bacterial chromosome.

Differential roles of positive and negative supercoiling in organizing the E. coli genome

This study aims to explore whether and how positive and negative supercoiling contribute to the three-dimensional (3D) organization of the bacterial genome. We used recently published Escherichia coli GapR ChIP-seq and TopoI ChIP-seq (also called EcTopoI-seq) data, which marks positive and negative supercoiling sites, respectively, to study how supercoiling correlates with the spatial contact maps obtained from chromosome conformation capture sequencing (Hi-C and 5C). We find that supercoiled chromosomal loci have overall higher Hi-C contact frequencies than sites that are not supercoiled. Surprisingly, positive supercoiling corresponds to higher spatial contact than negative supercoiling. Additionally, positive, but not negative, supercoiling could be identified from Hi-C data with high accuracy. We further find that the majority of positive and negative supercoils coincide with highly active transcription units, with a minor group likely associated with replication and other genomic processes. Our results show that both positive and negative supercoiling enhance spatial contact, with positive supercoiling playing a larger role in bringing genomic loci closer in space. Based on our results, we propose new physical models of how the E. coli chromosome is organized by positive and negative supercoils.

Chromosome structure modeling tools and their evaluation in bacteria

The three-dimensional (3D) structure of bacterial chromosomes is crucial for understanding chromosome function. With the growing availability of high-throughput chromosome conformation capture (3C/Hi-C) data, the 3D structure reconstruction algorithms have become powerful tools to study bacterial chromosome structure and function. It is highly desired to have a recommendation on the chromosome structure reconstruction tools to facilitate the prokaryotic 3D genomics. In this work, we review existing chromosome 3D structure reconstruction algorithms and classify them based on their underlying computational models into two categories: constraint-based modeling and thermodynamics-based modeling. We briefly compare these algorithms utilizing 3C/Hi-C datasets and fluorescence microscopy data obtained from Escherichia coli and Caulobacter crescentus, as well as simulated datasets. We discuss current challenges in the 3D reconstruction algorithms for bacterial chromosomes, primarily focusing on software usability. Finally, we briefly prospect future research directions for bacterial chromosome structure reconstruction algorithms.

Dynamics of chromosome organization in a minimal bacterial cell

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculate in silico chromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM.

High-Resolution 3D Genome Map of Brucella Chromosomes in Exponential and Stationary Phases

The three-dimensional (3D) genome structure of an organism or cell is highly relevant to its biological activities, but the availability of 3D genome information for bacteria, especially intracellular pathogens, is still limited. Here, we used Hi-C (high-throughput chromosome conformation capture) technology to determine the 3D chromosome structures of exponential- and stationary-phase Brucella melitensis at a 1-kb resolution. We observed that the contact heat maps of the two B. melitensis chromosomes contain a prominent diagonal and a secondary diagonal. Then, 79 chromatin interaction domains (CIDs) were detected at an optical density at 600 nm (OD600) of 0.4 (exponential phase), with the longest CID being 106 kb and the shortest being 12 kb. Moreover, we obtained 49,363 significant cis-interaction loci and 59,953 significant trans-interaction loci. Meanwhile, 82 CIDs of B. melitensis at an OD600 of 1.5 (stationary phase) were detected, with the longest CID being 94 kb and the shortest being 16 kb. In addition, 25,965 significant cis-interaction loci and 35,938 significant trans-interaction loci were obtained in this phase. Furthermore, we found that as the B. melitensis cells grew from the logarithmic to the plateau phase, the frequency of short-range interactions increased, while that of long-range interactions decreased. Finally, combined analysis of 3D genome and whole-genome transcriptome (RNA-seq) data revealed that the strength of short-range interactions in Chr1 is specifically and strongly correlated with gene expression. Overall, our study provides a global view of the chromatin interactions in the B. melitensis chromosomes, which will serve as a resource for further study of the spatial regulation of gene expression in Brucella. IMPORTANCE The spatial structure of chromatin plays important roles in normal cell functions and in the regulation of gene expression. Three-dimensional genome sequencing has been performed in many mammals and plants, but the availability of such data for bacteria, especially intracellular pathogens, is still limited. Approximately 10% of sequenced bacterial genomes contain more than one replicon. However, how multiple replicons are organized within bacterial cells, how they interact, and whether these interactions help to maintain or segregate these multipartite genomes are unresolved issues. Brucella is a Gram-negative, facultative intracellular, and zoonotic bacterium. Except for Brucella suis biovar 3, Brucella species have two chromosomes. Here, we applied Hi-C technology to determine the 3D genome structures of exponential- and stationary-phase Brucella melitensis chromosomes at a 1-kb resolution. Combined analysis of the 3D genome and RNA-seq data indicated that the strength of short-range interactions in B. melitensis Chr1 is specifically and strongly correlated with gene expression. Our study provides a resource to achieve a deeper understanding of the spatial regulation of gene expression in Brucella.

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

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

Mycobacterial IHF is a highly dynamic nucleoid-associated protein that assists HupB in organizing chromatin

Nucleoid-associated proteins (NAPs) crucially contribute to organizing bacterial chromatin and regulating gene expression. Among the most highly expressed NAPs are the HU and integration host factor (IHF) proteins, whose functional homologues, HupB and mycobacterial integration host factor (mIHF), are found in mycobacteria. Despite their importance for the pathogenicity and/or survival of tubercle bacilli, the role of these proteins in mycobacterial chromosome organization remains unknown. Here, we used various approaches, including super-resolution microscopy, to perform a comprehensive analysis of the roles of HupB and mIHF in chromosome organization. We report that HupB is a structural agent that maintains chromosome integrity on a local scale, and that the lack of this protein alters chromosome morphology. In contrast, mIHF is a highly dynamic protein that binds DNA only transiently, exhibits susceptibility to the chromosomal DNA topology changes and whose depletion leads to the growth arrest of tubercle bacilli. Additionally, we have shown that depletion of Mycobacterium smegmatis integration host factor (msIHF) leads to chromosome shrinkage and replication inhibition.

Topology-driven spatial organization of ring polymers under confinement

Entropic repulsion between DNA ring polymers under confinement is a key mechanism governing the spatial segregation of bacterial DNA before cell division. Here we establish that "internal" loops within a modified-ring polymer architecture enhance entropic repulsion between two overlapping polymers confined in a cylinder. Interestingly, they also induce entropy-driven spatial organization of polymer segments as seen in vivo. Here we design polymers of different architectures in our simulations by introducing a minimal number of cross-links between specific monomers along the ring-polymer contour. The cross-links are likely induced by various bridging proteins inside living cells. We investigate the segregation of two polymers with modified topologies confined in a cylinder, which initially had spatially overlapping configurations. This helps us to identify the architectures that lead to higher success rates of segregation. We also establish the mechanism that leads to localization of specific polymer segments. We use the blob model to provide a theoretical understanding of why certain architectures lead to enhanced entropic repulsive forces between the polymers. Lastly, we establish a correspondence between the organizational patterns of the chromosome of the C.crescentus bacterium and our results for a specifically designed polymer architecture. However, the principles outlined here pertaining to the organization of polymeric segments are applicable to both synthetic and biological polymers.

Polymer architecture orchestrates the segregation and spatial organization of replicating E. coli chromosomes in slow growth

The mechanism of chromosome segregation and organization in the bacterial cell cycle of E. coli is one of the least understood aspects in its life cycle. The E. coli chromosome is often modelled as a bead spring ring polymer. We introduce cross-links in the DNA-ring polymer, resulting in the formation of loops within each replicating bacterial chromosome. We use simulations to show that the chosen polymer-topology ensures its self-organization along the cell long-axis, such that various chromosomal loci get spatially localized as seen in vivo. The localization of loci arises due to entropic repulsion between polymer loops within each daughter DNA confined in a cylinder. The cellular addresses of the loci in our model are in fair agreement with those seen in experiments as given in J. A. Cass et al., Biophys. J., 2016, 110, 2597-2609. We also show that the adoption of such modified polymer architectures by the daughter DNAs leads to an enhanced propensity of their spatial segregation. Secondly, we match other experimentally reported results, including observation of the cohesion time and the ter-transition. Additionally, the contact map generated from our simulations reproduces the macro-domain like organization as seen in the experimentally obtained Hi-C map. Lastly, we have also proposed a plausible reconciliation of the 'Train Track' and the 'Replication Factory' models which provide conflicting descriptions of the spatial organization of the replication forks. Thus, we reconcile observations from complementary experimental techniques probing bacterial chromosome organization.

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

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

Relationship between the Chromosome Structural Dynamics and Gene Expression—A Chicken and Egg Dilemma?

Prokaryotic transcription was extensively studied over the last half-century. A great deal of data has been accumulated regarding the control of gene expression by transcription factors regulating their target genes by binding at specific DNA sites. However, there is a significant gap between the mechanistic description of transcriptional control obtained from in vitro biochemical studies and the complexity of transcriptional regulation in the context of the living cell. Indeed, recent studies provide ample evidence for additional levels of complexity pertaining to the regulation of transcription in vivo, such as, for example, the role of the subcellular localization and spatial organization of different molecular components involved in the transcriptional control and, especially, the role of chromosome configurational dynamics. The question as to how the chromosome is dynamically reorganized under the changing environmental conditions and how this reorganization is related to gene expression is still far from being clear. In this article, we focus on the relationships between the chromosome structural dynamics and modulation of gene expression during bacterial adaptation. We argue that spatial organization of the bacterial chromosome is of central importance in the adaptation of gene expression to changing environmental conditions and vice versa, that gene expression affects chromosome dynamics.