DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling.

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Diversification of Ubiquinone Biosynthesis via Gene Duplications, Transfers, Losses, and Parallel Evolution.

The availability of an ever-increasing diversity of prokaryotic genomes and metagenomes represents a major opportunity to understand and decipher the mechanisms behind the functional diversification of microbial biosynthetic pathways. However, it remains unclear to what extent a pathway producing a specific molecule from a specific precursor can diversify. In this study, we focus on the biosynthesis of ubiquinone (UQ), a crucial coenzyme that is central to the bioenergetics and to the functioning of a wide variety of enzymes in Eukarya and Pseudomonadota (a subgroup of the formerly named Proteobacteria). UQ biosynthesis involves three hydroxylation reactions on contiguous carbon atoms. We and others have previously shown that these reactions are catalyzed by different sets of UQ-hydroxylases that belong either to the iron-dependent Coq7 family or to the more widespread flavin monooxygenase (FMO) family. Here, we combine an experimental approach with comparative genomics and phylogenetics to reveal how UQ-hydroxylases evolved different selectivities within the constrained framework of the UQ pathway. It is shown that the UQ-FMOs diversified via at least three duplication events associated with two cases of neofunctionalization and one case of subfunctionalization, leading to six subfamilies with distinct hydroxylation selectivity. We also demonstrate multiple transfers of the UbiM enzyme and the convergent evolution of UQ-FMOs toward the same function, which resulted in two independent losses of the Coq7 ancestral enzyme. Diversification of this crucial biosynthetic pathway has therefore occurred via a combination of parallel evolution, gene duplications, transfers, and losses.

Assessing in vivo the impact of gene context on transcription through DNA supercoiling.

Gene context can have significant impact on gene expression but is currently not integrated in quantitative models of gene regulation despite known biophysical principles and quantitative in vitro measurements. Conceptually, the simplest gene context consists of a single gene framed by two topological barriers, known as the twin transcriptional-loop model, which illustrates the interplay between transcription and DNA supercoiling. In vivo, DNA supercoiling is additionally modulated by topoisomerases, whose modus operandi remains to be quantified. Here, we bridge the gap between theory and in vivo properties by realizing in Escherichia coli the twin transcriptional-loop model and by measuring how gene expression varies with promoters and distances to the topological barriers. We find that gene expression depends on the distance to the upstream barrier but not to the downstream barrier, with a promoter-dependent intensity. We rationalize these findings with a first-principle biophysical model of DNA transcription. Our results are explained if TopoI and gyrase both act specifically, respectively upstream and downstream of the gene, with antagonistic effects of TopoI, which can repress initiation while facilitating elongation. Altogether, our work sets the foundations for a systematic and quantitative description of the impact of gene context on gene regulation.

The context of genes, particularly the arrangement of neighboring genes along the DNA, exerts an important impact on their expression. However, predicting this impact remains challenging due to the complex interplay of concurrent mechanisms. To gain a quantitative understanding, we experimentally implemented the simplest possible theoretical model, isolating a gene from its neighboring genes. This allowed us to investigate the role of DNA's mechanical and topological properties, along with the enzymes that shape these properties, including RNA polymerases and topoisomerases. Comparison of the experimental results to a mathematical model based on physical principles allowed us to parametrize the operating mode of topoisomerases. Our work paves the way towards a systematic understanding of the role of gene context in gene expression.

Computational Tools for the Multiscale Analysis of Hi-C Data in Bacterial Chromosomes.

Just as in eukaryotes, high-throughput chromosome conformation capture (Hi-C) data have revealed nested organizations of bacterial chromosomes into overlapping interaction domains. In this chapter, we present a multiscale analysis framework aiming at capturing and quantifying these properties. These include both standard tools (e.g., contact laws) and novel ones such as an index that allows identifying loci involved in domain formation independently of the structuring scale at play. Our objective is twofold. On the one hand, we aim at providing a full, understandable Python/Jupyter-based code which can be used by both computer scientists and biologists with no advanced computational background. On the other hand, we discuss statistical issues inherent to Hi-C data analysis, focusing more particularly on how to properly assess the statistical significance of results. As a pedagogical example, we analyze data produced in Pseudomonas aeruginosa, a model pathogenetic bacterium. All files (codes and input data) can be found on a GitHub repository. We have also embedded the files into a Binder package so that the full analysis can be run on any machine through Internet.

Dynamics of the compartmentalized Streptomyces chromosome during metabolic differentiation.

Bacteria of the genus Streptomyces are prolific producers of specialized metabolites, including antibiotics. The linear chromosome includes a central region harboring core genes, as well as extremities enriched in specialized metabolite biosynthetic gene clusters. Here, we show that chromosome structure in Streptomyces ambofaciens correlates with genetic compartmentalization during exponential phase. Conserved, large and highly transcribed genes form boundaries that segment the central part of the chromosome into domains, whereas the terminal ends tend to be transcriptionally quiescent compartments with different structural features. The onset of metabolic differentiation is accompanied by a rearrangement of chromosome architecture, from a rather 'open' to a 'closed' conformation, in which highly expressed specialized metabolite biosynthetic genes form new boundaries. Thus, our results indicate that the linear chromosome of S. ambofaciens is partitioned into structurally distinct entities, suggesting a link between chromosome folding, gene expression and genome evolution.

Multiscale Dynamic Structuring of Bacterial Chromosomes.

Since the nucleoid was isolated from bacteria in the 1970s, two fundamental questions emerged and are still in the spotlight: how bacteria organize their chromosomes to fit inside the cell and how nucleoid organization enables essential biological processes. During the last decades, knowledge of bacterial chromosome organization has advanced considerably, and today, such chromosomes are considered to be highly organized and dynamic structures that are shaped by multiple factors in a multiscale manner. Here we review not only the classical well-known factors involved in chromosome organization but also novel components that have recently been shown to dynamically shape the 3D structuring of the bacterial genome. We focus on the different functional elements that control short-range organization and describe how they collaborate in the establishment of the higher-order folding and disposition of the chromosome. Recent advances have opened new avenues for a deeper understanding of the principles and mechanisms of chromosome organization in bacteria.

Supercoiled DNA and non-equilibrium formation of protein complexes: A quantitative model of the nucleoprotein ParBS partition complex.

ParABS, the most widespread bacterial DNA segregation system, is composed of a centromeric sequence, parS, and two proteins, the ParA ATPase and the ParB DNA binding proteins. Hundreds of ParB proteins assemble dynamically to form nucleoprotein parS-anchored complexes that serve as substrates for ParA molecules to catalyze positioning and segregation events. The exact nature of this ParBS complex has remained elusive, what we address here by revisiting the Stochastic Binding model (SBM) introduced to explain the non-specific binding profile of ParB in the vicinity of parS. In the SBM, DNA loops stochastically bring loci inside a sharp cluster of ParB. However, previous SBM versions did not include the negative supercoiling of bacterial DNA, leading to use unphysically small DNA persistences to explain the ParB binding profiles. In addition, recent super-resolution microscopy experiments have revealed a ParB cluster that is significantly smaller than previous estimations and suggest that it results from a liquid-liquid like phase separation. Here, by simulating the folding of long (≥ 30 kb) supercoiled DNA molecules calibrated with realistic DNA parameters and by considering different possibilities for the physics of the ParB cluster assembly, we show that the SBM can quantitatively explain the ChIP-seq ParB binding profiles without any fitting parameter, aside from the supercoiling density of DNA, which, remarkably, is in accord with independent measurements. We also predict that ParB assembly results from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., ParB clusters behave like liquid-like protein condensates with unconventional "leaky" boundaries.

Distinct Activities of Bacterial Condensins for Chromosome Management in Pseudomonas aeruginosa.

Three types of structurally related structural maintenance of chromosomes (SMC) complexes, referred to as condensins, have been identified in bacteria. Smc-ScpAB is present in most bacteria, whereas MukBEF is found in enterobacteria and MksBEF is scattered over the phylogenic tree. The contributions of these condensins to chromosome management were characterized in Pseudomonas aeruginosa, which carries both Smc-ScpAB and MksBEF. In this bacterium, SMC-ScpAB controls chromosome disposition by juxtaposing chromosome arms. In contrast, MksBEF is critical for chromosome segregation in the absence of the main segregation system, and it affects the higher-order architecture of the chromosome by promoting DNA contacts in the megabase range. Strikingly, our results reveal a prevalence of Smc-ScpAB over MksBEF involving a coordination of their activities with chromosome replication. They also show that E. coli MukBEF can substitute for MksBEF in P. aeruginosa while prevailing over Smc-ScpAB. Our results reveal a hierarchy between activities of bacterial condensins on the same chromosome.

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.

Intracellular Positioning Systems Limit the Entropic Eviction of Secondary Replicons Toward the Nucleoid Edges in Bacterial Cells.

Bacterial genomes, organized intracellularly as nucleoids, are composed of the main chromosome coexisting with different types of secondary replicons. Secondary replicons are major drivers of bacterial adaptation by gene exchange. They are highly diverse in type and size, ranging from less than 2 to more than 1000 kb, and must integrate with bacterial physiology, including to the nucleoid dynamics, to limit detrimental costs leading to their counter-selection. We show that large DNA circles, whether from a natural plasmid or excised from the chromosome tend to localize in a dynamic manner in a zone separating the nucleoid from the cytoplasm at the edge of the nucleoid. This localization is in good agreement with silico simulations of DNA circles in the nucleoid volume. Subcellular positioning systems counteract this tendency, allowing replicons to enter the nucleoid space. In enterobacteria, these systems are found in replicons above 25 kb, defining the limit with small randomly segregated plasmids. Larger replicons carry at least one of the three described family of systems, ParAB, ParRM, and StbA. Replicons above 180 kb all carry a ParAB system, suggesting this system is specifically required in the cases of large replicons. Simulations demonstrated that replicon size profoundly affects localization, compaction, and dynamics of DNA circles in the nucleoid volume. The present work suggests that presence of partition systems on the larger plasmids or chromids is not only due to selection for accurate segregation but also to counteract their unmixing with the chromosome and consequent exclusion from the nucleoid.

A Two-Enzyme Adaptive Unit within Bacterial Folate Metabolism.

Enzyme function and evolution are influenced by the larger context of a metabolic pathway. Deleterious mutations or perturbations in one enzyme can often be compensated by mutations to others. We used comparative genomics and experiments to examine evolutionary interactions with the essential metabolic enzyme dihydrofolate reductase (DHFR). Analyses of synteny and co-occurrence across bacterial species indicate that DHFR is coupled to thymidylate synthase (TYMS) but relatively independent from the rest of folate metabolism. Using quantitative growth rate measurements and forward evolution in Escherichia coli, we demonstrate that the two enzymes adapt as a relatively independent unit in response to antibiotic stress. Metabolomic profiling revealed that TYMS activity must not exceed DHFR activity to prevent the depletion of reduced folates and the accumulation of the intermediate dihydrofolate. Comparative genomics analyses identified >200 gene pairs with similar statistical signatures of modular co-evolution, suggesting that cellular pathways may be decomposable into small adaptive units.

Universal and idiosyncratic characteristic lengths in bacterial genomes.

In condensed matter physics, simplified descriptions are obtained by coarse-graining the features of a system at a certain characteristic length, defined as the typical length beyond which some properties are no longer correlated. From a physics standpoint, in vitro DNA has thus a characteristic length of 300 base pairs (bp), the Kuhn length of the molecule beyond which correlations in its orientations are typically lost. From a biology standpoint, in vivo DNA has a characteristic length of 1000 bp, the typical length of genes. Since bacteria live in very different physico-chemical conditions and since their genomes lack translational invariance, whether larger, universal characteristic lengths exist is a non-trivial question. Here, we examine this problem by leveraging the large number of fully sequenced genomes available in public databases. By analyzing GC content correlations and the evolutionary conservation of gene contexts (synteny) in hundreds of bacterial chromosomes, we conclude that a fundamental characteristic length around 10-20 kb can be defined. This characteristic length reflects elementary structures involved in the coordination of gene expression, which are present all along the genome of nearly all bacteria. Technically, reaching this conclusion required us to implement methods that are insensitive to the presence of large idiosyncratic genomic features, which may co-exist along these fundamental universal structures.

Modeling Bacterial DNA: Simulation of Self-Avoiding Supercoiled Worm-Like Chains Including Structural Transitions of the Helix.

Under supercoiling constraints, naked DNA, such as a large part of bacterial DNA, folds into braided structures called plectonemes. The double-helix can also undergo local structural transitions, leading to the formation of denaturation bubbles and other alternative structures. Various polymer models have been developed to capture these properties, with Monte-Carlo (MC) approaches dedicated to the inference of thermodynamic properties. In this chapter, we explain how to perform such Monte-Carlo simulations, following two objectives. On one hand, we present the self-avoiding supercoiled Worm-Like Chain (ssWLC) model, which is known to capture the folding properties of supercoiled DNA, and provide a detailed explanation of a standard MC simulation method. On the other hand, we explain how to extend this ssWLC model to include structural transitions of the helix.

Evolution of Ubiquinone Biosynthesis: Multiple Proteobacterial Enzymes with Various Regioselectivities To Catalyze Three Contiguous Aromatic Hydroxylation Reactions.

The ubiquitous ATP synthase uses an electrochemical gradient to synthesize cellular energy in the form of ATP. The production of this electrochemical gradient relies on liposoluble proton carriers like ubiquinone (UQ), which is used in the respiratory chains of eukaryotes and proteobacteria. The biosynthesis of UQ requires three hydroxylation reactions on contiguous positions of an aromatic ring. In Escherichia coli, each of three UQ flavin monooxygenases (FMOs), called UbiF, UbiH, and UbiI, modifies a single position of the aromatic ring. This pattern of three hydroxylation reactions/three proteins has been accepted as a paradigm in UQ biology. Using a phylogenetic analysis, we found that UbiF, UbiH, and UbiI are detected only in a small fraction of proteobacteria, and we identified two new types of UQ FMOs: UbiM, which is distributed in members of the alpha, beta, and gamma classes of proteobacteria, and UbiL, which is restricted to members of the alphaproteobacteria. Remarkably, the ubiL and ubiM genes were found in genomes with fewer than three UQ hydroxylase-encoding genes. We demonstrated, using biochemical approaches, that UbiL from Rhodospirillum rubrum and UbiM from Neisseria meningitidis hydroxylate, respectively, two and three positions of the aromatic ring during UQ biosynthesis. We conclude that bacteria have evolved a large repertoire of hydroxylase combinations for UQ biosynthesis, including pathways with either three specialist enzymes or pathways with one or two generalist enzymes of broader regioselectivity. The emergence of the latter is potentially related to genome reduction events. IMPORTANCE UQ, a key molecule for cellular bioenergetics that is conserved from proteobacteria to humans, appeared in an ancestral proteobacterium more than 2 billion years ago. UQ biosynthesis has been studied only in a few model organisms, and thus, the diversity of UQ biosynthesis pathways is largely unknown. In the work reported here, we conducted a phylogenomic analysis of hydroxylases involved in UQ biosynthesis. Our results support the existence of at least two UQ hydroxylases in the proteobacterial ancestor, and yet, we show that their number varies from one to four in extant proteobacterial species. Our biochemical experiments demonstrated that bacteria containing only one or two UQ hydroxylases have developed generalist enzymes that are able to catalyze several steps of UQ biosynthesis. Our study documents a rare case where evolution favored the broadening of an enzyme's regioselectivity, which resulted in gene loss in several proteobacterial species with small genomes.

Insights into the Mechanisms of Basal Coordination of Transcription Using a Genome-Reduced Bacterium.

Coordination of transcription in bacteria occurs at supra-operonic scales, but the extent, specificity, and mechanisms of such regulation are poorly understood. Here, we tackle this problem by profiling the transcriptome of the model organism Mycoplasma pneumoniae across 115 growth conditions. We identify three qualitatively different levels of co-expression corresponding to distinct relative orientations and intergenic properties of adjacent genes. We reveal that the degree of co-expression between co-directional adjacent operons, and more generally between genes, is tightly related to their capacity to be transcribed en bloc into the same mRNA. We further show that this genome-wide pervasive transcription of adjacent genes and operons is specifically repressed by DNA regions preferentially bound by RNA polymerases, by intrinsic terminators, and by large intergenic distances. Taken together, our findings suggest that the basal coordination of transcription is mediated by the physical entities and mechanical properties of the transcription process itself, and that operon-like behaviors may strongly vary from condition to condition.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation.

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations.

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

From structure to function of bacterial chromosomes: Evolutionary perspectives and ideas for new experiments.

The link between chromosome structure and function is a challenging open question because chromosomes in vivo are highly dynamic and arduous to manipulate. Here, we examine several promising approaches to tackle this question specifically in bacteria, by integrating knowledge from different sources. Toward this end, we first provide a brief overview of experimental tools that have provided insights into the description of the bacterial chromosome, including genetic, biochemical and fluorescence microscopy techniques. We then explore the possibility of using comparative genomics to isolate functionally important features of chromosome organization, exploiting the fact that features shared between phylogenetically distant bacterial species reflect functional significance. Finally, we discuss possible future perspectives from the field of experimental evolution. Specifically, we propose novel experiments in which bacteria could be screened and selected on the basis of the structural properties of their chromosomes.

On the demultiplexing of chromosome capture conformation data.

How to describe the multiple chromosome structures that underlie interactions among genome loci and how to quantify the occurrence of these structures in a cell population remain important challenges to solve, which can be addressed via a proper demultiplexing of chromosome capture conformation related data. Here, we first aim to review two main methodologies that have been proposed to tackle this problem: restrained-based methods, in which the resulting chromosome structures stem from the multiple solutions of a distance satisfaction problem; and thermodynamic-based methods, in which the structures stem from the simulation of polymer models. Next, we propose a novel demultiplexing method based on a matrix decomposition of contact maps. To this end, we extend the notion of topologically associated domains (TADs) by introducing that of statistical interaction domains (SIDs). SIDs can overlap and occur in a cell population at certain frequencies, and we propose a simple method to estimate these frequency values. As an application, we show that SIDs that measure 100kb to tens of Mb long occur both frequently and specifically in the human genome.