Interconnecting solvent quality, transcription, and chromosome folding in Escherichia coli.

All cells fold their genomes, including bacterial cells, where the chromosome is compacted into a domain-organized meshwork called the nucleoid. How compaction and domain organization arise is not fully understood. Here, we describe a method to estimate the average mesh size of the nucleoid in Escherichia coli. Using nucleoid mesh size and DNA concentration estimates, we find that the cytoplasm behaves as a poor solvent for the chromosome when the cell is considered as a simple semidilute polymer solution. Monte Carlo simulations suggest that a poor solvent leads to chromosome compaction and DNA density heterogeneity (i.e., domain formation) at physiological DNA concentration. Fluorescence microscopy reveals that the heterogeneous DNA density negatively correlates with ribosome density within the nucleoid, consistent with cryoelectron tomography data. Drug experiments, together with past observations, suggest the hypothesis that RNAs contribute to the poor solvent effects, connecting chromosome compaction and domain formation to transcription and intracellular organization.

A Well-Mixed E. coli Genome: Widespread Contacts Revealed by Tracking Mu Transposition.

The three-dimensional structures of chromosomes are increasingly being recognized as playing a major role in cellular regulatory states. The efficiency and promiscuity of phage Mu transposition was exploited to directly measure in vivo interactions between genomic loci in E. coli. Two global organizing principles have emerged: first, the chromosome is well-mixed and uncompartmentalized, with transpositions occurring freely between all measured loci; second, several gene families/regions show "clustering": strong three-dimensional co-localization regardless of linear genomic distance. The activities of the SMC/condensin protein MukB and nucleoid-compacting protein subunit HU-α are essential for the well-mixed state; HU-α is also needed for clustering of 6/7 ribosomal RNA-encoding loci. The data are explained by a model in which the chromosomal structure is driven by dynamic competition between DNA replication and chromosomal relaxation, providing a foundation for determining how region-specific properties contribute to both chromosomal structure and gene regulation.

Machine learning unravels inherent structural patterns in Escherichia coli Hi-C matrices and predicts chromosome dynamics.

High dimensional nature of the chromosomal conformation contact map ('Hi-C Map'), even for microscopically small bacterial cell, poses challenges for extracting meaningful information related to its complex organization. Here we first demonstrate that an artificial deep neural network-based machine-learnt (ML) low-dimensional representation of a recently reported Hi-C interaction map of archetypal bacteria Escherichia coli can decode crucial underlying structural pattern. The ML-derived representation of Hi-C map can automatically detect a set of spatially distinct domains across E. coli genome, sharing reminiscences of six putative macro-domains previously posited via recombination assay. Subsequently, a ML-generated model assimilates the intricate relationship between large array of Hi-C-derived chromosomal contact probabilities and respective diffusive dynamics of each individual chromosomal gene and identifies an optimal number of functionally important chromosomal contact-pairs that are majorly responsible for heterogenous, coordinate-dependent sub-diffusive motions of chromosomal loci. Finally, the ML models, trained on wild-type E. coli show-cased its predictive capabilities on mutant bacterial strains, shedding light on the structural and dynamic nuances of ΔMatP30MM and ΔMukBEF22MM chromosomes. Overall our results illuminate the power of ML techniques in unraveling the complex relationship between structure and dynamics of bacterial chromosomal loci, promising meaningful connections between ML-derived insights and biological phenomena.

Phase-separated ParB enforces diverse DNA compaction modes and stabilizes the parS-centered partition complex.

The tripartite ParABS system mediates chromosome segregation in the majority of bacterial species. Typically, DNA-bound ParB proteins around the parS sites condense the chromosomal DNA into a higher-order multimeric nucleoprotein complex for the ParA-driven partition. Despite extensive studies, the molecular mechanism underlying the dynamic assembly of the partition complex remains unclear. Herein, we demonstrate that Bacillus subtilis ParB (Spo0J), through the multimerization of its N-terminal domain, forms phase-separated condensates along a single DNA molecule, leading to the concurrent organization of DNA into a compact structure. Specifically, in addition to the co-condensation of ParB dimers with DNA, the engagement of well-established ParB condensates with DNA allows for the compression of adjacent DNA and the looping of distant DNA. Notably, the presence of CTP promotes the formation of condensates by a low amount of ParB at parS sites, triggering two-step DNA condensation. Remarkably, parS-centered ParB-DNA co-condensate constitutes a robust nucleoprotein architecture capable of withstanding disruptive forces of tens of piconewton. Overall, our findings unveil diverse modes of DNA compaction enabled by phase-separated ParB and offer new insights into the dynamic assembly and maintenance of the bacterial partition complex.

Transcription factor shapes chromosomal conformation and regulates gene expression in bacterial adaptation.

Genomic mutations allow bacteria to adapt rapidly to adverse stress environments. The three-dimensional conformation of the genome may also play an important role in transcriptional regulation and environmental adaptation. Here, using chromosome conformation capture, we investigate the high-order architecture of the Zymomonas mobilis chromosome in response to genomic mutation and ambient stimuli (acetic acid and furfural, derived from lignocellulosic hydrolysate). We find that genomic mutation only influences the local chromosome contacts, whereas stress of acetic acid and furfural restrict the long-range contacts and significantly change the chromosome organization at domain scales. Further deciphering the domain feature unveils the important transcription factors, Ferric uptake regulator (Fur) proteins, which act as nucleoid-associated proteins to promote long-range (>200 kb) chromosomal communications and regulate the expression of genes involved in stress response. Our work suggests that ubiquitous transcription factors in prokaryotes mediate chromosome organization and regulate stress-resistance genes in bacterial adaptation.

Insights into the molecular mechanism of ParABS system in chromosome partition by HpParA and HpParB.

The ParABS system, composed of ParA (an ATPase), ParB (a DNA binding protein), and parS (a centromere-like DNA), regulates bacterial chromosome partition. The ParB-parS partition complex interacts with the nucleoid-bound ParA to form the nucleoid-adaptor complex (NAC). In Helicobacter pylori, ParA and ParB homologs are encoded as HpSoj and HpSpo0J (HpParA and HpParB), respectively. We determined the crystal structures of the ATP hydrolysis deficient mutant, HpParAD41A, and the HpParAD41A-DNA complex. We assayed the CTPase activity of HpParB and identified two potential DNA binding modes of HpParB regulated by CTP, one is the specific DNA binding by the DNA binding domain and the other is the non-specific DNA binding through the C-terminal domain under the regulation of CTP. We observed an interaction between HpParAD41A and the N-terminus fragment of HpParB (residue 1-10, HpParBN10) and determined the crystal structure of the ternary complex, HpParAD41A-DNA-HpParBN10 complex which mimics the NAC formation. HpParBN10 binds near the HpParAD41A dimer interface and is clamped by flexible loops, L23 and L34, through a specific cation-π interaction between Arg9 of HpParBN10 and Phe52 of HpParAD41A. We propose a molecular mechanism model of the ParABS system providing insight into chromosome partition in bacteria.

Tuned SMC Arms Drive Chromosomal Loading of Prokaryotic Condensin.

SMC proteins support vital cellular processes in all domains of life by organizing chromosomal DNA. They are composed of ATPase "head" and "hinge" dimerization domains and a connecting coiled-coil "arm." Binding to a kleisin subunit creates a closed tripartite ring, whose ∼47-nm-long SMC arms act as barrier for DNA entrapment. Here, we uncover another, more active function of the bacterial Smc arm. Using high-throughput genetic engineering, we resized the arm in the range of 6-60 nm and found that it was functional only in specific length regimes following a periodic pattern. Natural SMC sequences reflect these length constraints. Mutants with improper arm length or peptide insertions in the arm efficiently target chromosomal loading sites and hydrolyze ATP but fail to use ATP hydrolysis for relocation onto flanking DNA. We propose that SMC arms implement force transmission upon nucleotide hydrolysis to mediate DNA capture or loop extrusion.

CRISPR-Cas Systems Optimize Their Immune Response by Specifying the Site of Spacer Integration.

CRISPR-Cas systems defend prokaryotes against viruses and plasmids. Short DNA segments of the invader, known as spacers, are stored in the CRISPR array as immunological memories. New spacers are added invariably to the 5' end of the array; therefore, the first spacer matches the latest foreign threat. Whether this highly polarized order of spacer insertion influences CRISPR-Cas immunity has not been explored. Here we show that a conserved sequence located immediately upstream of the CRISPR array specifies the site of new spacer integration. Mutation of this sequence results in erroneous incorporation of new spacers into the middle of the array. We show that spacers added through polarized acquisition give rise to more robust CRISPR-Cas immunity than spacers added to the middle of the array. This study demonstrates that the CRISPR-Cas system specifies the site of spacer integration to optimize the immune response against the most immediate threat to the host.

Colocalization of distant chromosomal loci in space in E. coli: a bacterial nucleolus.

The spatial organization of DNA within the bacterial nucleoid remains unclear. To investigate chromosome organization in Escherichia coli, we examined the relative positions of the ribosomal RNA (rRNA) operons in space. The seven rRNA operons are nearly identical and separated from each other by as much as 180° on the circular genetic map, a distance of ≥2 million base pairs. By inserting binding sites for fluorescent proteins adjacent to the rRNA operons and then examining their positions pairwise in live cells by epifluorescence microscopy, we found that all but rrnC are in close proximity. Colocalization of the rRNA operons required the rrn P1 promoter region but not the rrn P2 promoter or the rRNA structural genes and occurred with and without active transcription. Non-rRNA operon pairs did not colocalize, and the magnitude of their physical separation generally correlated with that of their genetic separation. Our results show that E. coli bacterial chromosome folding in three dimensions is not dictated entirely by genetic position but rather includes functionally related, genetically distant loci that come into close proximity, with rRNA operons forming a structure reminiscent of the eukaryotic nucleolus.

The MatP/matS site-specific system organizes the terminus region of the E. coli chromosome into a macrodomain.

The organization of the Escherichia coli chromosome into insulated macrodomains influences the segregation of sister chromatids and the mobility of chromosomal DNA. Here, we report that organization of the Terminus region (Ter) into a macrodomain relies on the presence of a 13 bp motif called matS repeated 23 times in the 800-kb-long domain. matS sites are the main targets in the E. coli chromosome of a newly identified protein designated MatP. MatP accumulates in the cell as a discrete focus that colocalizes with the Ter macrodomain. The effects of MatP inactivation reveal its role as main organizer of the Ter macrodomain: in the absence of MatP, DNA is less compacted, the mobility of markers is increased, and segregation of Ter macrodomain occurs early in the cell cycle. Our results indicate that a specific organizational system is required in the Terminus region for bacterial chromosome management during the cell cycle.

Pathogenicity island-directed transfer of unlinked chromosomal virulence genes.

In recent decades, the notorious pathogen Staphylococcus aureus has become progressively more contagious, more virulent, and more resistant to antibiotics. This implies a rather dynamic evolutionary capability, representing a remarkable level of genomic plasticity, most probably maintained by horizontal gene transfer. Here we report that the staphylococcal pathogenicity islands have a dual role in gene transfer: they not only mediate their own transfer, but they can independently direct the transfer of unlinked chromosomal segments containing virulence genes. While transfer of the island itself requires specific helper phages, transfer of unlinked chromosomal segments does not, so potentially any pac-type phage will serve. These results reveal that SaPIs can increase the horizontal exchange of accessory genes associated with disease and may shape pathogen genomes beyond the confines of their attachment sites.

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

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

The coordinated replication of Vibrio cholerae's two chromosomes required the acquisition of a unique domain by the RctB initiator.

Vibrio cholerae, the pathogenic bacterium that causes cholera, has two chromosomes (Chr1, Chr2) that replicate in a well-orchestrated sequence. Chr2 initiation is triggered only after the replication of the crtS site on Chr1. The initiator of Chr2 replication, RctB, displays activities corresponding with its different binding sites: initiator at the iteron sites, repressor at the 39m sites, and trigger at the crtS site. The mechanism by which RctB relays the signal to initiate Chr2 replication from crtS is not well-understood. In this study, we provide new insights into how Chr2 replication initiation is regulated by crtS via RctB. We show that crtS (on Chr1) acts as an anti-inhibitory site by preventing 39m sites (on Chr2) from repressing initiation. The competition between these two sites for RctB binding is explained by the fact that RctB interacts with crtS and 39m via the same DNA-binding surface. We further show that the extreme C-terminal tail of RctB, essential for RctB self-interaction, is crucial for the control exerted by crtS. This subregion of RctB is conserved in all Vibrio, but absent in other Rep-like initiators. Hence, the coordinated replication of both chromosomes likely results from the acquisition of this unique domain by RctB.

Genome replication dynamics of a bacteriophage and its satellite reveal strategies for parasitism and viral restriction.

Phage-inducible chromosomal island-like elements (PLEs) are bacteriophage satellites found in Vibrio cholerae. PLEs parasitize the lytic phage ICP1, excising from the bacterial chromosome, replicating, and mobilizing to new host cells following cell lysis. PLEs protect their host cell populations by completely restricting the production of ICP1 progeny. Previously, it was found that ICP1 replication was reduced during PLE(+) infection. Despite robust replication of the PLE genome, relatively few transducing units are produced. We investigated if PLE DNA replication itself is antagonistic to ICP1 replication. Here we identify key constituents of PLE replication and assess their role in interference of ICP1. PLE encodes a RepA_N initiation factor that is sufficient to drive replication from the PLE origin of replication during ICP1 infection. In contrast to previously characterized bacteriophage satellites, expression of the PLE initiation factor was not sufficient for PLE replication in the absence of phage. Replication of PLE was necessary for interference of ICP1 DNA replication, but replication of a minimalized PLE replicon was not sufficient for ICP1 DNA replication interference. Despite restoration of ICP1 DNA replication, non-replicating PLE remained broadly inhibitory against ICP1. These results suggest that PLE DNA replication is one of multiple mechanisms contributing to ICP1 restriction.

Biphasic unbinding of a metalloregulator from DNA for transcription (de)repression in Live Bacteria.

Microorganisms use zinc-sensing regulators to alter gene expression in response to changes in the availability of zinc, an essential micronutrient. Under zinc-replete conditions, the Fur-family metalloregulator Zur binds to DNA tightly in its metallated repressor form to Zur box operator sites, repressing the transcription of zinc uptake transporters. Derepression comes from unbinding of the regulator, which, under zinc-starvation conditions, exists in its metal-deficient non-repressor forms having no significant affinity with Zur box. While the mechanism of transcription repression by Zur is well-studied, little is known on how derepression by Zur could be facilitated. Using single-molecule/single-cell measurements, we find that in live Escherichia coli cells, Zur's unbinding rate from DNA is sensitive to Zur protein concentration in a first-of-its-kind biphasic manner, initially impeded and then facilitated with increasing Zur concentration. These results challenge conventional models of protein unbinding being unimolecular processes and independent of protein concentration. The facilitated unbinding component likely occurs via a ternary complex formation mechanism. The impeded unbinding component likely results from Zur oligomerization on chromosome involving inter-protein salt-bridges. Unexpectedly, a non-repressor form of Zur is found to bind chromosome tightly, likely at non-consensus sequence sites. These unusual behaviors could provide functional advantages in Zur's facile switching between repression and derepression.

A nucleotide-dependent oligomerization of the Escherichia coli replication initiator DnaA requires residue His136 for remodeling of the chromosomal origin.

Escherichia coli replication initiator protein DnaA binds ATP with high affinity but the amount of ATP required to initiate replication greatly exceeds the amount required for binding. Previously, we showed that ATP-DnaA, not ADP-DnaA, undergoes a conformational change at the higher nucleotide concentration, which allows DnaA oligomerization at the replication origin but the association state remains unclear. Here, we used Small Angle X-ray Scattering (SAXS) to investigate oligomerization of DnaA in solution. Whereas ADP-DnaA was predominantly monomeric, AMP-PNP-DnaA (a non-hydrolysable ATP-analog bound-DnaA) was oligomeric, primarily dimeric. Functional studies using DnaA mutants revealed that DnaA(H136Q) is defective in initiating replication in vivo. The mutant retains high-affinity ATP binding, but was defective in producing replication-competent initiation complexes. Docking of ATP on a structure of E. coli DnaA, modeled upon the crystallographic structure of Aquifex aeolicus DnaA, predicts a hydrogen bond between ATP and imidazole ring of His136, which is disrupted when Gln is present at position 136. SAXS performed on AMP-PNP-DnaA (H136Q) indicates that the protein has lost its ability to form oligomers. These results show the importance of high ATP in DnaA oligomerization and its dependence on the His136 residue.

Nonspecific DNA binding by P1 ParA determines the distribution of plasmid partition and repressor activities.

The faithful segregation, or "partition," of many low-copy number bacterial plasmids is driven by plasmid-encoded ATPases that are represented by the P1 plasmid ParA protein. ParA binds to the bacterial nucleoid via an ATP-dependent nonspecific DNA (nsDNA)-binding activity, which is essential for partition. ParA also has a site-specific DNA-binding activity to the par operator (parOP), which requires either ATP or ADP, and which is essential for it to act as a transcriptional repressor but is dispensable for partition. Here we examine how DNA binding by ParA contributes to the relative distribution of its plasmid partition and repressor activities, using a ParA with an alanine substitution at Arg351, a residue previously predicted to participate in site-specific DNA binding. In vivo, the parAR351A allele is compromised for partition, but its repressor activity is dramatically improved so that it behaves as a "super-repressor." In vitro, ParAR351A binds and hydrolyzes ATP, and undergoes a specific conformational change required for nsDNA binding, but its nsDNA-binding activity is significantly damaged. This defect in turn significantly reduces the assembly and stability of partition complexes formed by the interaction of ParA with ParB, the centromere-binding protein, and DNA. In contrast, the R351A change shows only a mild defect in site-specific DNA binding. We conclude that the partition defect is due to altered nsDNA binding kinetics and affinity for the bacterial chromosome. Furthermore, the super-repressor phenotype is explained by an increased pool of non-nucleoid bound ParA that is competent to bind parOP and repress transcription.

The role of Helicobacter pylori DnaA domain I in orisome assembly on a bipartite origin of chromosome replication.

The main roles of the DnaA protein are to bind the origin of chromosome replication (oriC), to unwind DNA and to provide a hub for the step-wise assembly of a replisome. DnaA is composed of four domains, with each playing a distinct functional role in the orisome assembly. Out of the four domains, the role of domain I is the least understood and appears to be the most species-specific. To better characterise Helicobacter pylori DnaA domain I, we have constructed a series of DnaA variants and studied their interactions with H. pylori bipartite oriC. We show that domain I is responsible for the stabilisation and organisation of DnaA-oriC complexes and provides cooperativity in DnaA-DNA interactions. Domain I mediates cross-interactions between oriC subcomplexes, which indicates that domain I is important for long-distance DnaA interactions and is essential for orisosme assembly on bipartite origins. HobA, which interacts with domain I, increases the DnaA binding to bipartite oriC; however, it does not stimulate but rather inhibits DNA unwinding. This suggests that HobA helps DnaA to bind oriC, but an unknown factor triggers DNA unwinding. Together, our results indicate that domain I self-interaction is important for the DnaA assembly on bipartite H. pylori oriC.

PresRAT: a server for identification of bacterial small-RNA sequences and their targets with probable binding region.

Bacterial small-RNA (sRNA) sequences are functional RNAs, which play an important role in regulating the expression of a diverse class of genes. It is thus critical to identify such sRNA sequences and their probable mRNA targets. Here, we discuss new procedures to identify and characterize sRNA and their targets via the introduction of an integrated online platform 'PresRAT'. PresRAT uses the primary and secondary structural attributes of sRNA sequences to predict sRNA from a given sequence or bacterial genome. PresRAT also finds probable target mRNAs of sRNA sequences from a given bacterial chromosome and further concentrates on the identification of the probable sRNA-mRNA binding regions. Using PresRAT, we have identified a total of 66,209 potential sRNA sequences from 292 bacterial genomes and 2247 potential targets from 13 bacterial genomes. We have also implemented a protocol to build and refine 3D models of sRNA and sRNA-mRNA duplex regions and generated 3D models of 50 known sRNAs and 81 sRNA-mRNA duplexes using this platform. Along with the server part, PresRAT also contains a database section, which enlists the predicted sRNA sequences, sRNA targets, and their corresponding 3D models with structural dynamics information.

Collapse transition of a heterogeneous polymer in a crowded medium.

Long chain molecules can be entropically compacted in a crowded medium. We study the compaction transition of a heterogeneous polymer with ring topology by crowding effects in a free or confined space. For this, we use molecular dynamics simulations in which the effects of crowders are taken into account through effective interactions between chain segments. Our parameter choices are inspired by the Escherichia coli chromosome. The polymer consists of small and big monomers; the big monomers dispersed along the backbone are to mimic the binding of RNA polymerases. Our results show that the compaction transition is a two-step process: initial compaction induced by the association (clustering) of big monomers followed by a gradual overall compaction. They also indicate that cylindrical confinement makes the initial transition more effective; for representative parameter choices, the initial compaction accounts for about 60% reduction in the chain size. Our simulation results support the view that crowding promotes clustering of active transcription units into transcription factories.