DNA looping in prokaryotes: experimental and theoretical approaches

The unreasonable effectiveness of equilibrium gene regulation through the cell cycle

Systems like the prototypical lac operon can reliably hold repression of transcription upon DNA replication across cell cycles with just 10 repressor molecules per cell and behave as if they were at equilibrium. The origin of this phenomenology is still an unresolved question. Here, we develop a general theory to analyze strong perturbations in quasi-equilibrium systems and use it to quantify the effects of DNA replication in gene regulation. We find a scaling law linking actual with predicted equilibrium transcription via a single kinetic parameter. We show that even the lac operon functions beyond the physical limits of naive regulation through compensatory mechanisms that suppress non-equilibrium effects. Synthetic systems without adjuvant activators, such as the cAMP receptor protein (CRP), lack this reliability. Our results provide a rationale for the function of CRP, beyond just being a tunable activator, as a mitigator of cell cycle perturbations.

Molecular insights into the fine-tuning of pH-dependent ArsR-mediated regulation of the SabA adhesin in Helicobacter pylori

Adaptation to variations in pH is crucial for the ability of Helicobacter pylori to persist in the human stomach. The acid responsive two-component system ArsRS, constitutes the global regulon that responds to acidic conditions, but molecular details of how transcription is affected by the ArsR response regulator remains poorly understood. Using a combination of DNA-binding studies, in vitro transcription assays, and H. pylori mutants, we demonstrate that phosphorylated ArsR (ArsR-P) forms an active protein complex that binds DNA with high specificity in order to affect transcription. Our data showed that DNA topology is key for DNA binding. We found that AT-rich DNA sequences direct ArsR-P to specific sites and that DNA-bending proteins are important for the effect of ArsR-P on transcription regulation. The repression of sabA transcription is mediated by ArsR-P with the support of Hup and is affected by simple sequence repeats located upstream of the sabA promoter. Here stochastic events clearly contribute to the fine-tuning of pH-dependent gene regulation. Our results reveal important molecular aspects for how ArsR-P acts to repress transcription in response to acidic conditions. Such transcriptional control likely mediates shifts in bacterial positioning in the gastric mucus layer.

Spatial proximity and gene function: a new dimension in prokaryotic gene association network analysis with 3D-GeneNet

Understanding the biological functions and processes of genes, particularly those not yet characterized, is crucial for advancing molecular biology and identifying therapeutic targets. The hypothesis guiding this study is that the 3D proximity of genes correlates with their functional interactions and relevance in prokaryotes. We introduced 3D-GeneNet, an innovative software tool that utilizes high-throughput sequencing data from chromosome conformation capture techniques and integrates topological metrics to construct gene association networks. Through a series of comparative analyses focused on spatial versus linear distances, we explored various dimensions such as topological structure, functional enrichment levels, distribution patterns of linear distances among gene pairs, and the area under the receiver operating characteristic curve by utilizing model organism Escherichia coli K-12. Furthermore, 3D-GeneNet was shown to maintain good accuracy compared to multiple algorithms (neighbourhood, co-occurrence, coexpression, and fusion) across multiple bacteria, including E. coli, Brucella abortus, and Vibrio cholerae. In addition, the accuracy of 3D-GeneNet's prediction of long-distance gene interactions was identified by bacterial two-hybrid assays on E. coli K-12 MG1655, where 3D-GeneNet not only increased the accuracy of linear genomic distance tripled but also achieved 60% accuracy by running alone. Finally, it can be concluded that the applicability of 3D-GeneNet will extend to various bacterial forms, including Gram-negative, Gram-positive, single-, and multi-chromosomal bacteria through Hi-C sequencing and analysis. Such findings highlight the broad applicability and significant promise of this method in the realm of gene association network. 3D-GeneNet is freely accessible at https://github.com/gaoyuanccc/3D-GeneNet.

Nanoscale Interaction of Endonuclease APE1 with DNA

Apurinic/apyrimidinic endonuclease 1 (APE1) is involved in DNA repair and transcriptional regulation mechanisms. This multifunctional activity of APE1 should be supported by specific structural properties of APE1 that have not yet been elucidated. Herein, we applied atomic force microscopy (AFM) to characterize the interactions of APE1 with DNA containing two well-separated G-rich segments. Complexes of APE1 with DNA containing G-rich segments were visualized, and analysis of the complexes revealed the affinity of APE1 to G-rich DNA sequences, and their yield was as high as 53%. Furthermore, APE1 is capable of binding two DNA segments leading to the formation of loops in the DNA-APE1 complexes. The analysis of looped APE1-DNA complexes revealed that APE1 can bridge G-rich segments of DNA. The yield of loops bridging two G-rich DNA segments was 41%. Analysis of protein size in various complexes was performed, and these data showed that loops are formed by APE1 monomer, suggesting that APE1 has two DNA binding sites. The data led us to a model for the interaction of APE1 with DNA and the search for the specific sites. The implication of these new APE1 properties in organizing DNA, by bringing two distant sites together, for facilitating the scanning for damage and coordinating repair and transcription is discussed.

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.

Global Effects of the Developmental Regulator BldB in Streptomyces venezuelae

In Streptomyces, the Bld (Bald) regulators control formation of the reproductive aerial hyphae. The functions of some of these regulators have been well characterized, but BldB has remained enigmatic. In addition to the bldB gene itself, Streptomyces venezuelae has 10 paralogs of bldB that sit next to paralogs of whiJ and abaA. Transcriptome sequencing (RNA-seq) revealed that loss of BldB function causes the dramatic transcriptional upregulation of the abaA paralogs and a novel inhibitor of sporulation, iosA, and that cooverexpression of just two of these genes, iosA and abaA6, was sufficient to recapitulate the bldB mutant phenotype. Further RNA-seq analysis showed that the transcription factor WhiJ9 is required for the activation of iosA seen in the bldB mutant, and biochemical studies showed that WhiJ9 mediates the activation of iosA expression by binding to direct repeats in the iosA-whiJ9 intergenic region. BldB and BldB9 hetero-oligomerize, providing a potential link between BldB and the iosA-whiJ9-bldB9 locus. This work greatly expands our overall understanding of the global effects of the BldB developmental regulator. IMPORTANCE To reproduce and disperse, the filamentous bacterium Streptomyces develops specialized reproductive structures called aerial hyphae. The formation of these structures is controlled by the bld (bald) genes, many of which encode transcription factors whose functions have been characterized. An exception is BldB, a protein whose biochemical function is unknown. In this study, we gain insight into the global effects of BldB function by examining the genome-wide transcriptional effects of deleting bldB. We identify a small set of genes that are dramatically upregulated in the absence of BldB. We show that their overexpression causes the bldB phenotype and characterize a transcription factor that mediates the upregulation of one of these target genes. Our results provide new insight into how BldB influences Streptomyces development.

A tetramerization domain in prokaryotic and eukaryotic transcription regulators homologous to p53

Transcriptional regulation usually requires the action of several proteins that either repress or activate a promotor of an open reading frame. These proteins can counteract each other, thus allowing tight regulation of the transcription of the corresponding genes, where tight repression is often linked to DNA looping or cross-linking. Here, the tetramerization domain of the bacterial gene repressor Rco from Bacillus subtilis plasmid pLS20 (RcopLS20) has been identified and its structure is shown to share high similarity to the tetramerization domain of the well known p53 family of human tumor suppressors, despite lacking clear sequence homology. In RcopLS20, this tetramerization domain is responsible for inducing DNA looping, a process that involves multiple tetramers. In accordance, it is shown that RcopLS20 can form octamers. This domain was named TetDloop and its occurrence was identified in other Bacillus species. The TetDloop fold was also found in the structure of a transcriptional repressor from Salmonella phage SPC32H. It is proposed that the TetDloop fold has evolved through divergent evolution and that the TetDloop originates from a common ancestor predating the occurrence of multicellular life.

Active and thermal fluctuations in multi-scale polymer structure and dynamics

The presence of athermal noise or biological fluctuations control and maintain crucial life-processes. In this work, we present an exact analytical treatment of the dynamic behavior of a flexible polymer chain that is subjected to both thermal and active forces. Our model for active forces incorporates temporal correlation associated with the characteristic time scale and processivity of enzymatic function (driven by ATP hydrolysis), leading to an active-force time scale that competes with relaxation processes within the polymer chain. We analyze the structure and dynamics of an active-Brownian polymer using our exact results for the dynamic structure factor and the looping time for the chain ends. The spectrum of relaxation times within a polymer chain implies two different behaviors at small and large length scales. Small length-scale relaxation is faster than the active-force time scale, and the dynamic and structural behavior at these scales are oblivious to active forces and, are thus governed by the true thermal temperature. Large length-scale behavior is governed by relaxation times that are much longer than the active-force time scale, resulting in an effective active-Brownian temperature that dramatically alters structural and dynamic behavior. These complex multi-scale effects imply a time-dependent temperature that governs living and non-equilibrium systems, serving as a unifying concept for interpreting and predicting their physical behavior.

RNA polymerase efficiently transcribes DNA-scaffolded, cooperative bacteriophage repressor complexes

DNA can act as a scaffold for the cooperative binding of protein oligomers. For example, the phage 186 CI repressor forms a wheel of seven dimers wrapped in DNA with specific binding sites, while phage λ CI repressor dimers bind to two well-separated sets of operators, forming a DNA loop. Atomic force microscopy was used to measure transcription elongation by E. coli RNA polymerase through these protein complexes. 186 CI, or λ CI, bound along unlooped DNA negligibly interfered with transcription by RNAP. Wrapped and looped topologies induced by these scaffolded, cooperatively bound repressor oligomers did not form significantly better roadblocks to transcription. Thus, despite binding with high affinity, these repressors are not effective roadblocks to transcription.

Molecular insights into the role of heme in the transcriptional regulatory system AppA/PpsR

Heme has been shown to have a crucial role in the signal transduction mechanism of the facultative photoheterotrophic bacterium Rhodobacter sphaeroides. It interacts with the transcriptional regulatory complex AppA/PpsR in which AppA and PpsR function as the antirepressor and repressor, respectively of photosynthesis gene expression. The mechanism, however of this interaction remains incompletely understood. In this study, we combined EPR spectroscopy and FRET to demonstrate the ligation of heme in PpsR with a proposed cysteine residue. We show that heme binding in AppA affects the fluorescent properties of the dark-adapted state of the protein, suggesting a less constrained flavin environment compared to the absence of heme and the light-adapted state. We performed ultrafast transient absorption measurements in order to reveal potential differences in the dynamic processes in the full-length AppA and its heme-binding domain alone. Comparison of the CO-binding dynamics demonstrates a more open heme pocket in the holo-protein, qualitatively similar to what has been observed in the CO sensor RcoM-2, and suggests a communication path between the BLUF and SCHIC domains of AppA. We have also examined quantitatively, the affinity of PpsR to bind to individual DNA fragments of the puc promoter using fluorescence anisotropy assays. We conclude that oligomerization of PpsR is initially triggered by binding of one of the two DNA fragments and observe a ∼10-fold increase in the dissociation constant K_d for DNA binding upon heme binding to PpsR. Our study provides significant new insight at the molecular level on the regulatory role of heme that modulates the complex transcriptional regulation in R. sphaeroides and supports the two levels of heme signaling, via its binding to AppA and PpsR and via the sensing of gases like oxygen.

A chromosomal loop anchor mediates bacterial genome organization

Nucleoprotein complexes play an integral role in genome organization of both eukaryotes and prokaryotes. Apart from their role in locally structuring and compacting DNA, several complexes are known to influence global organization by mediating long-range anchored chromosomal loop formation leading to spatial segregation of large sections of DNA. Such megabase-range interactions are ubiquitous in eukaryotes, but have not been demonstrated in prokaryotes. Here, using a genome-wide sedimentation-based approach, we found that a transcription factor, Rok, forms large nucleoprotein complexes in the bacterium Bacillus subtilis. Using chromosome conformation capture and live-imaging of DNA loci, we show that these complexes robustly interact with each other over large distances. Importantly, these Rok-dependent long-range interactions lead to anchored chromosomal loop formation, thereby spatially isolating large sections of DNA, as previously observed for insulator proteins in eukaryotes.

Single-molecule insights into torsion and roadblocks in bacterial transcript elongation

During transcription, RNA polymerase (RNAP) translocates along the helical template DNA while maintaining high transcriptional fidelity. However, all genomes are dynamically twisted, writhed, and decorated by bound proteins and motor enzymes. In prokaryotes, proteins bound to DNA, specifically or not, frequently compact DNA into conformations that may silence genes by obstructing RNAP. Collision of RNAPs with these architectural proteins, may result in RNAP stalling and/or displacement of the protein roadblock. It is important to understand how rapidly transcribing RNAPs operate under different levels of supercoiling or in the presence of roadblocks. Given the broad range of asynchronous dynamics exhibited by transcriptional complexes, single-molecule assays, such as atomic force microscopy, fluorescence detection, optical and magnetic tweezers, etc. are well suited for detecting and quantifying activity with adequate spatial and temporal resolution. Here, we summarize current understanding of the effects of torsion and roadblocks on prokaryotic transcription, with a focus on single-molecule assays that provide real-time detection and readout.

Anticipating response function in gene regulatory networks

The origin of an ordered genetic response of a complex and noisy biological cell is intimately related to the detailed mechanism of protein-DNA interactions present in a wide variety of gene regulatory (GR) systems. However, the quantitative prediction of genetic response and the correlation between the mechanism and the response curve is poorly understood. Here, we report in silico binding studies of GR systems to show that the transcription factor (TF) binds to multiple DNA sites with high cooperativity spreads from specific binding sites into adjacent non-specific DNA and bends the DNA. Our analysis is not limited only to the isolated model system but also can be applied to a system containing multiple interacting genes. The controlling role of TF oligomerization, TF-ligand interactions, and DNA looping for gene expression has been also characterized. The predictions are validated against detailed grand canonical Monte Carlo simulations and published data for the lac operon system. Overall, our study reveals that the expression of target genes can be quantitatively controlled by modulating TF-ligand interactions and the bending energy of DNA.

DNA looping by protamine follows a nonuniform spatial distribution

DNA looping plays an important role in cells in both regulating and protecting the genome. Often, studies of looping focus on looping by prokaryotic transcription factors like lac repressor or by structural maintenance of chromosomes proteins such as condensin. Here, however, we are interested in a different looping method whereby condensing agents (charge ≥+3) such as protamine proteins neutralize the DNA, causing it to form loops and toroids. We considered two previously proposed mechanisms for DNA looping by protamine. In the first mechanism, protamine stabilizes spontaneous DNA fluctuations, forming randomly distributed loops along the DNA. In the second mechanism, protamine binds and bends the DNA to form a loop, creating a distribution of loops that is biased by protamine binding. To differentiate between these mechanisms, we imaged both spontaneous and protamine-induced loops on short-length (≤1 μm) DNA fragments using atomic force microscopy. We then compared the spatial distribution of the loops to several model distributions. A random looping model, which describes the mechanism of spontaneous DNA folding, fit the distribution of spontaneous loops, but it did not fit the distribution of protamine-induced loops. Specifically, it failed to predict a peak in the spatial distribution of loops at an intermediate location along the DNA. An electrostatic multibinding model, which was created to mimic the bind-and-bend mechanism of protamine, was a better fit of the distribution of protamine-induced loops. In this model, multiple protamines bind to the DNA electrostatically within a particular region along the DNA to coordinate the formation of a loop. We speculate that these findings will impact our understanding of protamine’s in vivo role for looping DNA into toroids and the mechanism of DNA condensation by condensing agents more broadly.

DNA Looping Mediated by Site-Specific SfiI-DNA Interactions

Interactions between distant DNA segments play important roles in various biological processes, such as DNA recombination. Certain restriction enzymes create DNA loops when two sites are held together and then cleave the DNA. DNA looping is important during DNA synapsis. Here we investigated the mechanisms of DNA looping by restriction enzyme SfiI by measuring the properties of the system at various temperatures. Different sized loop complexes, mediated by SfiI-DNA interactions, were visualized with AFM. The experimental results revealed that small loops are more favorable compared to other loop sizes at all temperatures. Our theoretical model found that entropic cost dominates at all conditions, which explains the preference for short loops. Furthermore, specific loop sizes were predicted as favorable from an energetic point of view. These predictions were tested by experiments with transiently assembled SfiI loops on a substrate with a single SfiI site.

Lsr2 and Its Novel Paralogue Mediate the Adjustment of Mycobacterium smegmatis to Unfavorable Environmental Conditions

Nucleoid-associated proteins (NAPs) are the most abundant proteins involved in bacterial chromosome organization and global transcription regulation. The mycobacterial NAP family includes many diverse proteins; some are unique to actinobacteria, and many are crucial for survival under stress (e.g., HupB and Lsr2) and/or optimal growth conditions (e.g., mycobacterial integration host factor [mIHF]).

Lsr2 is a nucleoid-associated protein (NAP) that has been found strictly in actinobacteria, including mycobacteria. It is a functional homolog of histone-like nucleoid-structuring protein (H-NS); it acts as a DNA-bridging protein that plays a role in chromosomal organization and transcriptional regulation. To date, the studies on Lsr2 have focused mainly on Mycobacterium tuberculosis. In this study, we analyze the role of Lsr2 as a transcription factor in Mycobacterium smegmatis, a saprophytic bacterium whose natural habitat (soil and water) substantially differs from those of the obligatory mycobacterial pathogens. Our chromatin immunoprecipitation-sequencing (ChIP-seq) data revealed that Lsr2 binds preferentially to AT-rich regions of the M. smegmatis chromosome. We found that Lsr2 acts mainly as a repressor, controlling gene expression either directly by binding promoter regions or indirectly through DNA loop formation and DNA coating. One of the Lsr2-repressed genes encodes polyketide synthase (MSMEG_4727), which is involved in the synthesis of lipooligosaccharides (LOSs). An M. smegmatis strain deprived of Lsr2 produces more LOSs, which is mirrored by changes in the smoothness of cells and their susceptibilities to antibiotics. Unlike M. tuberculosis, M. smegmatis additionally encodes a paralogue of Lsr2, MSMEG_1060, which is a novel member of the mycobacterial NAP family. The Lsr2 and MSMEG_1060 proteins exhibit different DNA binding specificities and chromosomal localizations. Our results suggest that these proteins help M. smegmatis cells cope with stress conditions, including hypoxia and exposure to antibiotics. Thus, the present work provides novel insight into the role of Lsr2 paralogues in the ability of a saprophytic mycobacterial species to adjust to environmental changes.

IMPORTANCE Nucleoid-associated proteins (NAPs) are the most abundant proteins involved in bacterial chromosome organization and global transcription regulation. The mycobacterial NAP family includes many diverse proteins; some are unique to actinobacteria, and many are crucial for survival under stress (e.g., HupB and Lsr2) and/or optimal growth conditions (e.g., mycobacterial integration host factor [mIHF]). Here, we present a comprehensive study concerning two functional homologues of mycobacterial H-NS: Lsr2 and its paralogue from M. smegmatis, MSMEG_1060. We found that Lsr2 plays a role in transcriptional regulation, mainly by repressing gene expression via DNA loop formation and/or DNA-coating mechanisms. Intriguingly, the number of Lsr2-mediated genes was found to increase under hypoxia. Compared to Lsr2, MSMEG_1060 exhibits a different DNA binding specificity and chromosomal localization. Since tuberculosis remains a serious worldwide health problem, studies on stress response-mediating agents, such as Lsr2, may contribute to the development of novel antituberculosis drugs.

The loopometer: a quantitative in vivo assay for DNA-looping proteins

Proteins that can bring together separate DNA sites, either on the same or on different DNA molecules, are critical for a variety of DNA-based processes. However, there are no general and technically simple assays to detect proteins capable of DNA looping in vivo nor to quantitate their in vivo looping efficiency. Here, we develop a quantitative in vivo assay for DNA-looping proteins in Escherichia coli that requires only basic DNA cloning techniques and a LacZ assay. The assay is based on loop assistance, where two binding sites for the candidate looping protein are inserted internally to a pair of operators for the E. coli LacI repressor. DNA looping between the sites shortens the effective distance between the lac operators, increasing LacI looping and strengthening its repression of a lacZ reporter gene. Analysis based on a general model for loop assistance enables quantitation of the strength of looping conferred by the protein and its binding sites. We use this ‘loopometer’ assay to measure DNA looping for a variety of bacterial and phage proteins.

DNA Looping and DNA Conformational Fluctuations Can Accelerate Protein Target Search

Protein searching and binding to specific sites on DNA is a fundamentally important process that marks the beginning of all major cellular transformations. While the dynamics of protein-DNA interactions in in vitro settings is well investigated, the situation is much more complex for in vivo conditions because the DNA molecules in live cells are packed into chromosomal structures where they are undergoing strong dynamic and conformational fluctuations. In this work, we present a theoretical investigation on the role of DNA looping and DNA conformational fluctuations in the protein target search. It is based on a discrete-state stochastic analysis that allows for explicit calculations of dynamic properties, which is also supplemented by Monte Carlo computer simulations. It is found that for stronger nonspecific interactions between DNA and proteins the search occurs faster on the DNA looped conformation in comparison with the unlooped conformation, and the fastest search is observed when the loop is formed near the target site. It is also shown that DNA fluctuations between the looped and unlooped conformations influence the search dynamics, and this depends on the magnitude of conformational transition rates and on which conformation is more energetically stable. Physical-chemical arguments explaining these observations are presented. Our theoretical study suggests that the geometry and conformational changes in DNA are additional factors that might efficiently control the gene regulation processes.

The heme-binding protein PhuS transcriptionally regulates the Pseudomonas aeruginosa tandem sRNA prrF1,F2 locus

Pseudomonas aeruginosa is an opportunistic pathogen requiring iron for its survival and virulence. P. aeruginosa can acquire iron from heme via the nonredundant heme assimilation system and Pseudomonas heme uptake (Phu) systems. Heme transported by either the heme assimilation system or Phu system is sequestered by the cytoplasmic protein PhuS. Furthermore, PhuS has been shown to specifically transfer heme to the iron-regulated heme oxygenase HemO. As the PhuS homolog ShuS from Shigella dysenteriae was observed to bind DNA as a function of its heme status, we sought to further determine if PhuS, in addition to its role in regulating heme flux through HemO, functions as a DNA-binding protein. Herein, through a combination of chromatin immunoprecipitation–PCR, EMSA, and fluorescence anisotropy, we show that apo-PhuS but not holo-PhuS binds upstream of the tandem iron-responsive sRNAs prrF1,F2. Previous studies have shown the PrrF sRNAs are required for sparing iron for essential proteins during iron starvation. Furthermore, under certain conditions, a heme-dependent read through of the prrF1 terminator yields the longer PrrH transcript. Quantitative PCR analysis of P. aeruginosa WT and ΔphuS strains shows that loss of PhuS abrogates the heme-dependent regulation of PrrF and PrrH levels. Taken together, our data show that PhuS, in addition to its role in extracellular heme metabolism, also functions as a transcriptional regulator by modulating PrrF and PrrH levels in response to heme. This dual function of PhuS is central to integrating extracellular heme utilization into the PrrF/PrrH sRNA regulatory network that is critical for P. aeruginosa adaptation and virulence within the host.

Multiplexed characterization of rationally designed promoter architectures deconstructs combinatorial logic for IPTG-inducible systems

A crucial step towards engineering biological systems is the ability to precisely tune the genetic response to environmental stimuli. In the case of Escherichia coli inducible promoters, our incomplete understanding of the relationship between sequence composition and gene expression hinders our ability to predictably control transcriptional responses. Here, we profile the expression dynamics of 8269 rationally designed, IPTG-inducible promoters that collectively explore the individual and combinatorial effects of RNA polymerase and LacI repressor binding site strengths. We then fit a statistical mechanics model to measured expression that accurately models gene expression and reveals properties of theoretically optimal inducible promoters. Furthermore, we characterize three alternative promoter architectures and show that repositioning binding sites within promoters influences the types of combinatorial effects observed between promoter elements. In total, this approach enables us to deconstruct relationships between inducible promoter elements and discover practical insights for engineering inducible promoters with desirable characteristics.

Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.

Determinants of cyclization-decyclization kinetics of short DNA with sticky ends

Cyclization of DNA with sticky ends is commonly used to measure DNA bendability as a function of length and sequence, but how its kinetics depend on the rotational positioning of the sticky ends around the helical axis is less clear. Here, we measured cyclization (looping) and decyclization (unlooping) rates (k_loop and k_unloop) of DNA with sticky ends over three helical periods (100-130 bp) using single-molecule fluorescence resonance energy transfer (FRET). k_loop showed a nontrivial undulation as a function of DNA length whereas k_unloop showed a clear oscillation with a period close to the helical turn of DNA (∼10.5 bp). The oscillation of k_unloop was almost completely suppressed in the presence of gaps around the sticky ends. We explain these findings by modeling double-helical DNA as a twisted wormlike chain with a finite width, intrinsic curvature, and stacking interaction between the end base pairs. We also discuss technical issues for converting the FRET-based cyclization/decyclization rates to an equilibrium quantity known as the J factor that is widely used to characterize DNA bending mechanics.

Protamine loops DNA in multiple steps

Protamine proteins dramatically condense DNA in sperm to almost crystalline packing levels. Here, we measure the first step in the in vitro pathway, the folding of DNA into a single loop. Current models for DNA loop formation are one-step, all-or-nothing models with a looped state and an unlooped state. However, when we use a Tethered Particle Motion (TPM) assay to measure the dynamic, real-time looping of DNA by protamine, we observe the presence of multiple folded states that are long-lived (∼100 s) and reversible. In addition, we measure folding on DNA molecules that are too short to form loops. This suggests that protamine is using a multi-step process to loop the DNA rather than a one-step process. To visualize the DNA structures, we used an Atomic Force Microscopy (AFM) assay. We see that some folded DNA molecules are loops with a ∼10-nm radius and some of the folded molecules are partial loops—c-shapes or s-shapes—that have a radius of curvature of ∼10 nm. Further analysis of these structures suggest that protamine is bending the DNA to achieve this curvature rather than increasing the flexibility of the DNA. We therefore conclude that protamine loops DNA in multiple steps, bending it into a loop.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

The effect of obstacles in multi-site protein target search with DNA looping

Many fundamental biological processes are regulated by protein-DNA complexes called synaptosomes, which possess multiple interaction sites. Despite the critical importance of synaptosomes, the mechanisms of their formation are not well understood. Because of the multisite nature of participating proteins, it is widely believed that their search for specific sites on DNA involves the formation and breaking of DNA loops and sliding in the looped configurations. In reality, DNA in live cells is densely covered by other biological molecules that might interfere with the formation of synaptosomes. In this work, we developed a theoretical approach to evaluate the role of obstacles in the target search of multisite proteins when the formation of DNA loops and the sliding in looped configurations are possible. Our theoretical method is based on analysis of a discrete-state stochastic model that uses a master equations approach and extensive computer simulations. It is found that the obstacle slows down the search dynamics in the system when DNA loops are long-lived, but the effect is minimal for short-lived DNA loops. In addition, the relative positions of the target and the obstacle strongly influence the target search kinetics. Furthermore, the presence of the obstacle might increase the noise in the system. These observations are discussed using physical-chemical arguments. Our theoretical approach clarifies the molecular mechanisms of formation of protein-DNA complexes with multiple interactions sites.

Identification and characterization of sequence signatures in the Bacillus subtilis promoter Pylb for tuning promoter strength

Objective

To thoroughly characterize the P_ylb promoter and identify the elements that affect the promoter activity.

Result

The sequences flanking the − 35 and − 10 box of the P_ylb promoter were divided into six segments, and six random-scanning mutant promoter libraries fused to an enhanced green fluorescent protein EGFP were made and analyzed by flow cytometry. Our results showed that the four nucleotides flanking the − 35 box could mostly influence the promoter activity, and this influence was related to the GC content. The promoters mutated in these regions were successfully used for expressing the gene ophc2 encoding organophosphorus hydrolase (OPHC2) and the gene katA encoding catalase (KatA).

Conclusion

Our work identified and characterized the sequence signatures of the P_ylb promoter that could tune the promoter strength, providing further information for the potential application of this promoter. Meanwhile, the sequence signatures have the potential to be used for tuning gene expression in enzyme production, metabolic engineering, and synthetic biology.

Electronic supplementary material

The online version of this article (10.1007/s10529-019-02749-4) contains supplementary material, which is available to authorized users.

Twist-bend coupling, twist waves, and the shape of DNA loops

By combining analytical and numerical calculations, we investigate the minimal-energy shape of short DNA loops of approximately 100 base pairs (bp). We show that in these loops the excess twist density oscillates as a response to an imposed bending stress, as recently found in DNA minicircles and observed in nucleosomal DNA. These twist oscillations, here referred to as twist waves, are due to the coupling between twist and bending deformations, which in turn originates from the asymmetry between DNA major and minor grooves. We introduce a simple analytical variational shape that reproduces the exact loop energy up to the fourth significant digit and is in very good agreement with shapes obtained from coarse-grained simulations. We, finally, analyze the loop dynamics at room temperature, and show that the twist waves are robust against thermal fluctuations. They perform a normal diffusive motion, whose origin is briefly discussed.

Identification of a basal system for unwinding a bacterial chromosome origin

Genome duplication is essential for cell proliferation, and DNA synthesis is generally initiated by dedicated replication proteins at specific loci termed origins. In bacteria, the master initiator DnaA binds the chromosome origin (oriC) and unwinds the DNA duplex to permit helicase loading. However, despite decades of research it remained unclear how the information encoded within oriC guides DnaA-dependent strand separation. To address this fundamental question, we took a systematic genetic approach in vivo and identified the core set of essential sequence elements within the Bacillus subtilis chromosome origin unwinding region. Using this information, we then show in vitro that the minimal replication origin sequence elements are necessary and sufficient to promote the mechanical functions of DNA duplex unwinding by DnaA. Because the basal DNA unwinding system characterized here appears to be conserved throughout the bacterial domain, this discovery provides a framework for understanding oriC architecture, activity, regulation and diversity.

FoxH1 represses miR-430 during early embryonic development of zebrafish via non-canonical regulation

Background

FoxH1 is a forkhead transcription factor with conserved key functions in vertebrate mesoderm induction and left-right patterning downstream of the TGF-beta/Nodal signaling pathway. Binding of the forkhead domain (FHD) of FoxH1 to a highly conserved proximal sequence motif was shown to regulate target gene expression.

Results

We identify the conserved microRNA-430 family (miR-430) as a novel target of FoxH1. miR-430 levels are increased in foxH1 mutants, resulting in a reduced expression of transcripts that are targeted by miR-430 for degradation. To determine the underlying mechanism of miR-430 repression, we performed chromatin immunoprecipitation studies and overexpression experiments with mutant as well as constitutive active and repressive forms of FoxH1. Our studies reveal a molecular interaction of FoxH1 with miR-430 loci independent of the FHD. Furthermore, we show that previously described mutant forms of FoxH1 that disrupt DNA binding or that lack the C-terminal Smad Interaction Domain (SID) dominantly interfere with miR-430 repression, but not with the regulation of previously described FoxH1 targets.

Conclusions

We were able to identify the distinct roles of protein domains of FoxH1 in the regulation process of miR-430. We provide evidence that the indirect repression of miR-430 loci depends on the connection to a distal repressive chromosome environment via a non-canonical mode. The widespread distribution of such non-canonical binding sites of FoxH1, found not only in our study, argues against a function restricted to regulating miR-430 and for a more global role of FoxH1 in chromatin folding.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0683-z) contains supplementary material, which is available to authorized users.

Facilitation of DNA loop formation by protein-DNA non-specific interactions

Complex DNA topological structures, including polymer loops, are frequently observed in biological processes when protein molecules simultaneously bind to several distant sites on DNA. However, the molecular mechanisms of formation of these systems remain not well understood. Existing theoretical studies focus only on specific interactions between protein and DNA molecules at target sequences. However, the electrostatic origin of primary protein-DNA interactions suggests that interactions of proteins with all DNA segments should be considered. Here we theoretically investigate the role of non-specific interactions between protein and DNA molecules on the dynamics of loop formation. Our approach is based on analyzing a discrete-state stochastic model via a method of first-passage probabilities supplemented by Monte Carlo computer simulations. It is found that depending on a protein sliding length during the non-specific binding event three different dynamic regimes of the DNA loop formation might be observed. In addition, the loop formation time might be optimized by varying the protein sliding length, the size of the DNA molecule, and the position of the specific target sequences on DNA. Our results demonstrate the importance of non-specific protein-DNA interactions in the dynamics of DNA loop formations.

Bacterial gene control by DNA looping using engineered dimeric transcription activator like effector (TALE) proteins

Genetic switches must alternate between states whose probabilities are dependent on regulatory signals. Classical examples of transcriptional control in bacteria depend on repressive DNA loops anchored by proteins whose structures are sensitive to small molecule inducers or co-repressors. We are interested in exploiting these natural principles to engineer artificial switches for transcriptional control of bacterial genes. Here, we implement designed homodimeric DNA looping proteins ('Transcription Activator-Like Effector Dimers'; TALEDs) for this purpose in living bacteria. Using well-studied FKBP dimerization domains, we build switches that mimic regulatory characteristics of classical Escherichia coli lactose, galactose and tryptophan operon promoters, including induction or co-repression by small molecules. Engineered DNA looping using TALEDs is thus a new approach to tuning gene expression in bacteria. Similar principles may also be applicable for gene control in eukaryotes.

The biology and polymer physics underlying large-scale chromosome organization

Chromosome large-scale organization is a beautiful example of the interplay between physics and biology. DNA molecules are polymers and thus belong to the class of molecules for which physicists have developed models and formulated testable hypotheses to understand their arrangement and dynamic properties in solution, based on the principles of polymer physics. Biologists documented and discovered the biochemical basis for the structure, function and dynamic spatial organization of chromosomes in cells. The underlying principles of chromosome organization have recently been revealed in unprecedented detail using high-resolution chromosome capture technology that can simultaneously detect chromosome contact sites throughout the genome. These independent lines of investigation have now converged on a model in which DNA loops, generated by the loop extrusion mechanism, are the basic organizational and functional units of the chromosome.

Analytical expressions for the closure probability of a stiff wormlike chain for finite capture radius

Estimating the probability that two monomers of the same polymer chain are close together is a key ingredient to characterize intramolecular reactions and polymer looping. In the case of stiff wormlike polymers (rigid fluctuating elastic rods), for which end-to-end encounters are rare events, we derive an explicit analytical formula for the probability η(r_{c}) that the distance between the chain extremities is smaller than some capture radius r_{c}. The formula is asymptotically exact in the limit of stiff chains, and it leads to the identification of two distinct scaling regimes for the closure factor, originating from a strong variation of the fluctuations of the chain orientation at closure. Our theory is compatible with existing analytical results from the literature that cover the cases of a vanishing capture radius and of nearly fully extended chains.

Mechanisms of peptide sex pheromone regulation of conjugation in Enterococcus faecalis

In many gram positive bacteria, horizontal transfer and virulence are regulated by peptide-mediated cell-cell signaling. The heptapeptide cCF10 (C) activates conjugative transfer of the Enterococcus faecalis plasmid pCF10, whereas the iCF10 (I) peptide inhibits transfer. Both peptides bind to the same domain of the master transcription regulator PrgX, a repressor of transcription of the prgQ operon encoding conjugation genes. We show that repression of prgQ by PrgX tetramers requires formation of a pCF10 DNA loop where each of two PrgX DNA-binding sites is occupied by a dimer. I binding to PrgX enhances prgQ repression, while C binding has the opposite effect. Previous models suggested that differential effects of these two peptides on the PrgX oligomerization state accounted for their distinct functions. Our new results demonstrate that both peptides have similar, high-binding affinity for PrgX, and that both peptides actually promote formation of PrgX tetramers with higher DNA-binding affinity than Apo-PrgX. We propose that differences in repression ability of PrgX/peptide complexes result from subtle differences in the structures of DNA-bound PrgX/peptide complexes. Changes in the induction state of a donor cell likely results from replacement of one type of DNA-bound peptide/PrgX tetramer with the other.

Unexpected sequences and structures of mtDNA required for efficient transcription from the first heavy-strand promoter

Human mtDNA contains three promoters, suggesting a need for differential expression of the mitochondrial genome. Studies of mitochondrial transcription have used a reductionist approach, perhaps masking differential regulation. Here we evaluate transcription from light-strand (LSP) and heavy-strand (HSP1) promoters using templates that mimic their natural context. These studies reveal sequences upstream, hypervariable in the human population (HVR3), and downstream of the HSP1 transcription start site required for maximal yield. The carboxy-terminal tail of TFAM is essential for activation of HSP1 but not LSP. Images of the template obtained by atomic force microscopy show that TFAM creates loops in a discrete region, the formation of which correlates with activation of HSP1; looping is lost in tail-deleted TFAM. Identification of HVR3 as a transcriptional regulatory element may contribute to between-individual variability in mitochondrial gene expression. The unique requirement of HSP1 for the TFAM tail may enable its regulation by post-translational modifications.

DOI:http://dx.doi.org/10.7554/eLife.27283.001

Effect of transcription factor resource sharing on gene expression noise

Gene expression is intrinsically a stochastic (noisy) process with important implications for cellular functions. Deciphering the underlying mechanisms of gene expression noise remains one of the key challenges of regulatory biology. Theoretical models of transcription often incorporate the kinetics of how transcription factors (TFs) interact with a single promoter to impact gene expression noise. However, inside single cells multiple identical gene copies as well as additional binding sites can compete for a limiting pool of TFs. Here we develop a simple kinetic model of transcription, which explicitly incorporates this interplay between TF copy number and its binding sites. We show that TF sharing enhances noise in mRNA distribution across an isogenic population of cells. Moreover, when a single gene copy shares it’s TFs with multiple competitor sites, the mRNA variance as a function of the mean remains unaltered by their presence. Hence, all the data for variance as a function of mean expression collapse onto a single master curve independent of the strength and number of competitor sites. However, this result does not hold true when the competition stems from multiple copies of the same gene. Therefore, although previous studies showed that the mean expression follows a universal master curve, our findings suggest that different scenarios of competition bear distinct signatures at the level of variance. Intriguingly, the introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution. These results demonstrate the impact of limited availability of TF resource on the regulation of noise in gene expression.

Genetically identical cells, even when they are exposed to the same environmental conditions, display incredible diversity. Gene expression noise is attributed to be a key source of this phenotypic diversity. Transcriptional dynamics is a dominant source of expression noise. Although scores of theoretical and experimental studies have explored how noise is regulated at the level of transcription, most of them focus on the gene specific, cis regulatory elements, such as the number of transcription factor (TF) binding sites, their binding strength, etc. However, how the global properties of transcription, such as the limited availability of TFs impact noise in gene expression remains rather elusive. Here we build a theoretical model that incorporates the effect of limiting TF pool on gene expression noise. We find that competition between genes for TFs leads to enhanced variability in mRNA copy number across an isogenic population. Moreover, for gene copies sharing TFs with other competitor sites, mRNA variance as a function of the mean shows distinct imprints for one gene copy and multiple gene copies respectively. This stands in sharp contrast to the universal behavior found in mean expression irrespective of the different scenarios of competition. An interesting feature of competition is that introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution, which could lead to phenotypic variability.

Recent Advances in Helicobacter pylori Replication: Possible Implications in Adaptation to a Pathogenic Lifestyle and Perspectives for Drug Design

DNA replication is an important step in the life cycle of every cell that ensures the continuous flow of genetic information from one generation to the next. In all organisms, chromosome replication must be coordinated with overall cell growth. Helicobacter pylori growth strongly depends on its interaction with the host, particularly with the gastric epithelium. Moreover, H. pylori actively searches for an optimal microniche within a stomach, and it has been shown that not every microniche equally supports growth of this bacterium. We postulate that besides nutrients, H. pylori senses different, unknown signals, which presumably also affect chromosome replication to maintain H. pylori propagation at optimal ratio allowing H. pylori to establish a chronic, lifelong infection. Thus, H. pylori chromosome replication and particularly the regulation of this process might be considered important for bacterial pathogenesis. Here, we summarize our current knowledge of chromosome and plasmid replication in H. pylori and discuss the mechanisms responsible for regulating this key cellular process. The results of extensive studies conducted thus far allow us to propose common and unique traits in H. pylori chromosome replication. Interestingly, the repertoire of proteins involved in replication in H. pylori is significantly different to that in E. coli, strongly suggesting that novel factors are engaged in H. pylori chromosome replication and could represent attractive drug targets.

Random versus Cell Cycle-Regulated Replication Initiation in Bacteria: Insights from Studying Vibrio cholerae Chromosome 2

Bacterial chromosomes initiate replication at a fixed time in the cell cycle, whereas there is generally no particular time for plasmid replication initiation or chromosomal replication initiation from integrated plasmids. In bacteria with divided genomes, the replication system of one of the chromosomes typically resembles that of bacteria with undivided genomes, whereas the remaining chromosomes have plasmid-like replication systems. For example, in Vibrio cholerae, a bacterium with two chromosomes (chromosome 1 [Chr1] and Chr2), the Chr1 system resembles that of the Escherichia coli chromosome, and the Chr2 system resembles that of iteron-based plasmids. However, Chr2 still initiates replication at a fixed time in the cell cycle and thus offers an opportunity to understand the molecular basis for the difference between random and cell cycle-regulated modes of replication. Here we review studies of replication control in Chr2 and compare it to those of plasmids and chromosomes. We argue that although the Chr2 control mechanisms in many ways are reminiscent of those of plasmids, they also appear to combine more regulatory features than are found on a typical plasmid, including some that are more typical of chromosomes. One of the regulatory mechanisms is especially novel, the coordinated timing of replication initiation of Chr1 and Chr2, providing the first example of communication between chromosomes for replication initiation.

Pseudomonas aeruginosa AmrZ Binds to Four Sites in the algD Promoter, Inducing DNA-AmrZ Complex Formation and Transcriptional Activation

During late stages of cystic fibrosis pulmonary infections, Pseudomonas aeruginosa often overproduces the exopolysaccharide alginate, protecting the bacterial community from host immunity and antimicrobials. The transcription of the alginate biosynthesis operon is under tight control by a number of factors, including AmrZ, the focus of this study. Interestingly, multiple transcription factors interact with the far-upstream region of this promoter (PalgD), in which one AmrZ binding site has been identified previously. The mechanisms of AmrZ binding and subsequent activation remain unclear and require more-detailed investigation. In this study, in-depth examinations elucidated four AmrZ binding sites, and their disruption eliminated AmrZ binding and promoter activation. Furthermore, our in vitro fluorescence resonance energy transfer experiments suggest that AmrZ holds together multiple binding sites in PalgD and thereafter induces the formation of higher-order DNA-AmrZ complexes. To determine the importance of interactions between those AmrZ oligomers in the cell, a DNA phasing experiment was performed. PalgD transcription was significantly impaired when the relative phase between AmrZ binding sites was reversed (5 bp), while a full-DNA-turn insertion (10 bp) restored promoter activity. Taken together, the investigations presented here provide a deeper mechanistic understanding of AmrZ-mediated binding to PalgD IMPORTANCE: Overproduction of the exopolysaccharide alginate provides protection to Pseudomonas aeruginosa against antimicrobial treatments and is associated with chronic P. aeruginosa infections in the lungs of cystic fibrosis patients. In this study, we combined a variety of microbiological, genetic, biochemical, and biophysical approaches to investigate the activation of the alginate biosynthesis operon promoter by a key transcription factor named AmrZ. This study has provided important new information on the mechanism of activation of this extremely complex promoter.

Metal-responsive promoter DNA compaction by the ferric uptake regulator

Short-range DNA looping has been proposed to affect promoter activity in many bacterial species and operator configurations, but only few examples have been experimentally investigated in molecular detail. Here we present evidence for a metal-responsive DNA condensation mechanism controlled by the Helicobacter pylori ferric uptake regulator (Fur), an orthologue of the widespread Fur family of prokaryotic metal-dependent regulators. H. pylori Fur represses the transcription of the essential arsRS acid acclimation operon through iron-responsive oligomerization and DNA compaction, encasing the arsR transcriptional start site in a repressive macromolecular complex. A second metal-dependent regulator NikR functions as nickel-dependent anti-repressor at this promoter, antagonizing the binding of Fur to the operator elements responsible for the DNA condensation. The results allow unifying H. pylori metal ion homeostasis and acid acclimation in a mechanistically coherent model, and demonstrate, for the first time, the existence of a selective metal-responsive DNA compaction mechanism controlling bacterial transcriptional regulation.

The Fur protein regulates transcription of bacterial genes in response to metal ions. Here, the authors show that the Fur protein from Helicobacter pylori represses transcription by iron-responsive oligomerization and DNA compaction, encasing the transcriptional start site in a macromolecular complex.

Novel genomic rearrangements mediated by multiple genetic elements in Streptococcus pyogenes M23ND confer potential for evolutionary persistence

Symmetric genomic rearrangements around replication axes in genomes are commonly observed in prokaryotic genomes, including Group A Streptococcus (GAS). However, asymmetric rearrangements are rare. Our previous studies showed that the hypervirulent invasive GAS strain, M23ND, containing an inactivated transcriptional regulator system, covRS, exhibits unique extensive asymmetric rearrangements, which reconstructed a genomic structure distinct from other GAS genomes. In the current investigation, we identified the rearrangement events and examined the genetic consequences and evolutionary implications underlying the rearrangements. By comparison with a close phylogenetic relative, M18-MGAS8232, we propose a molecular model wherein a series of asymmetric rearrangements have occurred in M23ND, involving translocations, inversions and integrations mediated by multiple factors, viz., rRNA-comX (factor for late competence), transposons and phage-encoded gene segments. Assessments of the cumulative gene orientations and GC skews reveal that the asymmetric genomic rearrangements did not affect the general genomic integrity of the organism. However, functional distributions reveal re-clustering of a broad set of CovRS-regulated actively transcribed genes, including virulence factors and metabolic genes, to the same leading strand, with high confidence (p-value ~10-10). The re-clustering of the genes suggests a potential selection advantage for the spatial proximity to the transcription complexes, which may contain the global transcriptional regulator, CovRS, and other RNA polymerases. Their proximities allow for efficient transcription of the genes required for growth, virulence and persistence. A new paradigm of survival strategies of GAS strains is provided through multiple genomic rearrangements, while, at the same time, maintaining genomic integrity.

A part toolbox to tune genetic expression in Bacillus subtilis

Libraries of well-characterised components regulating gene expression levels are essential to many synthetic biology applications. While widely available for the Gram-negative model bacterium Escherichia coli, such libraries are lacking for the Gram-positive model Bacillus subtilis, a key organism for basic research and biotechnological applications. Here, we engineered a genetic toolbox comprising libraries of promoters, Ribosome Binding Sites (RBS), and protein degradation tags to precisely tune gene expression in B. subtilis. We first designed a modular Expression Operating Unit (EOU) facilitating parts assembly and modifications and providing a standard genetic context for gene circuits implementation. We then selected native, constitutive promoters of B. subtilis and efficient RBS sequences from which we engineered three promoters and three RBS sequence libraries exhibiting ∼14 000-fold dynamic range in gene expression levels. We also designed a collection of SsrA proteolysis tags of variable strength. Finally, by using fluorescence fluctuation methods coupled with two-photon microscopy, we quantified the absolute concentration of GFP in a subset of strains from the library. Our complete promoters and RBS sequences library comprising over 135 constructs enables tuning of GFP concentration over five orders of magnitude, from 0.05 to 700 μM. This toolbox of regulatory components will support many research and engineering applications in B. subtilis.

A Bayesian inference method for the analysis of transcriptional regulatory networks in metagenomic data

Background

Metagenomics enables the analysis of bacterial population composition and the study of emergent population features, such as shared metabolic pathways. Recently, we have shown that metagenomics datasets can be leveraged to characterize population-wide transcriptional regulatory networks, or meta-regulons, providing insights into how bacterial populations respond collectively to specific triggers. Here we formalize a Bayesian inference framework to analyze the composition of transcriptional regulatory networks in metagenomes by determining the probability of regulation of orthologous gene sequences. We assess the performance of this approach on synthetic datasets and we validate it by analyzing the copper-homeostasis network of Firmicutes species in the human gut microbiome.

Results

Assessment on synthetic datasets shows that our method provides a robust and interpretable metric for assessing putative regulation by a transcription factor on sets of promoter sequences mapping to an orthologous gene cluster. The inference framework integrates the regulatory contribution of secondary sites and can discern false positives arising from multiple instances of a clonal sequence. Posterior probabilities for orthologous gene clusters decline sharply when less than 20 % of mapped promoters have binding sites, but we introduce a sensitivity adjustment procedure to speed up computation that enhances regulation assessment in heterogeneous ortholog clusters. Analysis of the copper-homeostasis regulon governed by CsoR in the human gut microbiome Firmicutes reveals that CsoR controls itself and copper-translocating P-type ATPases, but not CopZ-type copper chaperones. Our analysis also indicates that CsoR frequently targets promoters with dual CsoR-binding sites, suggesting that it exploits higher-order binding conformations to fine-tune its activity.

Conclusions

We introduce and validate a method for the analysis of transcriptional regulatory networks from metagenomic data that enables inference of meta-regulons in a systematic and interpretable way. Validation of this method on the CsoR meta-regulon of gut microbiome Firmicutes illustrates the usefulness of the approach, revealing novel properties of the copper-homeostasis network in poorly characterized bacterial species and putting forward evidence of new mechanisms of DNA binding for this transcriptional regulator. Our approach will enable the comparative analysis of regulatory networks across metagenomes, yielding novel insights into the evolution of transcriptional regulatory networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-016-0082-8) contains supplementary material, which is available to authorized users.

Force distribution in a semiflexible loop

Loops undergoing thermal fluctuations are prevalent in nature. Ringlike or cross-linked polymers, cyclic macromolecules, and protein-mediated DNA loops all belong to this category. Stability of these molecules are generally described in terms of free energy, an average quantity, but it may also be impacted by local fluctuating forces acting within these systems. The full distribution of these forces can thus give us insights into mechanochemistry beyond the predictive capability of thermodynamics. In this paper, we study the force exerted by an inextensible semiflexible polymer constrained in a looped state. By using a simulation method termed "phase-space sampling," we generate the equilibrium distribution of chain conformations in both position and momentum space. We compute the constraint forces between the two ends of the loop in this chain ensemble using Lagrangian mechanics, and show that the mean of these forces is equal to the thermodynamic force. By analyzing kinetic and potential contributions to the forces, we find that the mean force acts in the direction of increasing extension not because of bending stress, but in spite of it. Furthermore, we obtain a distribution of constraint forces as a function of chain length, extension, and stiffness. Notably, increasing contour length decreases the average force, but the additional freedom allows fluctuations in the constraint force to increase. The force distribution is asymmetric and falls off less sharply than a Gaussian distribution. Our work exemplifies a system where large-amplitude fluctuations occur in a way unforeseen by a purely thermodynamic framework, and offers computational tools useful for efficient, unbiased simulation of a constrained system.

Genetic dissection of independent and cooperative transcriptional activation by the LysR-type activator ThnR at close divergent promoters

Regulation of tetralin biodegradation operons is one of the examples of unconventional LysR-type mediated transcriptional regulation. ThnR activates transcription from two divergent and closely located promoters PB and PC. Although ThnR activates each promoter independently, transcription from each one increases when both promoters are together. Mutational analysis of the intergenic region shows that cooperative transcription is achieved through formation of a ThnR complex when bound to its respective sites at each promoter, via formation of a DNA loop. Mutations also defined ThnR contact sites that are important for independent transcriptional activation at each promoter. A mutation at the PB promoter region, which abolishes its independent transcription, does not affect at all PB transcription in the presence of the divergent promoter PC, thus indicating that the complex formed via DNA loop can compensate for the deficiencies in the correct protein-DNA interaction at one of the promoters. Combination of mutations in both promoters identifies a region at PC that is not important for its independent transcription but it is essential for cooperative transcription from both promoters. This work provides new insights into the diversity and complexity of activation mechanisms used by the most abundant type of bacterial transcriptional regulators.

NagR Differentially Regulates the Expression of the glmS and nagAB Genes Required for Amino Sugar Metabolism by Streptococcus mutans

The ability of bacteria to metabolize glucosamine (GlcN) and N-acetyl-d-glucosamine (GlcNAc) is considered important for persistent colonization of the oral cavity. In the dental caries pathogen Streptococcus mutans, the NagR protein regulates the expression of glmS, which encodes a GlcN-6-P synthetase, and nagA (GlcNAc-6-P deacetylase) and nagB (GlcN-6-P deaminase), which are required for the catabolism of GlcNAc and GlcN. Two NagR-binding sites (dre) were identified in each of the promoter regions for nagB and glmS. Using promoter-reporter gene fusions, the role of each dre site was examined in the regulation of glmS and nagB promoter activities in cells grown with glucose, GlcNAc, or GlcN. A synergistic relationship between the two dre sites in the glmS promoter that required proper spacing was observed, but that was not the case for nagB. Binding of purified NagR to DNA fragments from both promoter regions, as well as to dre sites alone, was strongly influenced by particular sugar phosphates. Using a random mutagenesis approach that targeted the effector-binding domain of NagR, mutants that displayed aberrant regulation of both the glmS and nagAB genes were identified. Collectively, these findings provide evidence that NagR is essential for regulation of genes for both the synthesis and catabolism of GlcN and GlcNAc in S. mutans, and that NagR engages differently with the target promoter regions in response to specific metabolites interacting with the effector-binding domain of NagR.

IMPORTANCE

Glucosamine and N-acetylglucosamine are among the most abundant naturally occurring sugars on the planet, and they are catabolized by many bacterial species as sources of carbon and nitrogen. Representing a group called lactic acid bacteria (LAB), the human dental caries pathogen Streptococcus mutans is shown to differ from known paradigm organisms in that it possesses a GntR/HutC-type regulator, NagR, that is required for the regulation of both catabolism of GlcN and biosynthesis. Results reported here reveal a simple and elegant mechanism whereby NagR differentially regulates two opposing biological processes by surveying metabolic intermediates. This study provides insights that may contribute to the development of novel therapeutic tools to combat dental caries and other infectious diseases.

Growth phase dependent changes in the structure and protein composition of nucleoid in Escherichia coli

The genomic DNA of bacteria is highly compacted in a single or a few bodies known as nucleoids. Here, we have isolated Escherichia coli nucleoid by sucrose density gradient centrifugation. The sedimentation rates, structures as well as protein/ DNA composition of isolated nucleoids were then compared under various growth phases. The nucleoid structures were found to undergo changes during the cell growth; i. e., the nucleoid structure in the stationary phase was more tightly compacted than that in the exponential phase. In addition to factor for inversion stimulation (Fis), histone-like nucleoid structuring protein (H-NS), heat-unstable nucleoid protein (HU) and integration host factor (IHF) here we have identified, three new candidates of E. coli nucleoid, namely DNA-binding protein from starved cells (Dps), host factor for phage Qβ (Hfq) and suppressor of td(-) phenotype A (StpA). Our results reveal that the major components of exponential phase nucleoid are Fis, HU, H-NS, StpA and Hfq, while Dps occupies more than half of the stationary phase nucleoid. It has been known for a while that Dps is the main nucleoid-associated protein at stationary phase. From these results and the prevailing information, we propose a model for growth phase dependent changes in the structure and protein composition of nucleoid in E. coli.