DNA looping

Polymer dynamics under tension: Mean first passage time for looping

This study deals with polymer looping, an important process in many chemical and biological systems. We investigate basic questions on the looping dynamics of a polymer under tension using the freely jointed chain (FJC) model. Previous theoretical approaches to polymer looping under tension have relied on barrier escape methods, which assume local equilibrium, an assumption that may not always hold. As a starting point, we use an analytical expression for the equilibrium looping probability as a function of the number of monomers and applied force, predicting an inverse relationship between looping time and looping probability. Using molecular dynamics simulations, the predictions of this theoretical approach are validated within the numerical precision achieved. We compare our predictions to those of the barrier escape approach by way of a calculation of the mean first passage time (MFPT) for the ends of a polymer to cross. For this purpose, we derive the exact free energy landscape, but the resulting temporal predictions do not agree with the observed inverse scaling. We conclude that the traditional barrier escape approach does not provide satisfactory predictions for polymer looping dynamics and that the inverse scaling with looping probability offers a more reliable alternative.

Single cell variant to enhancer to gene map for coronary artery disease

Although genome wide association studies (GWAS) in large populations have identified hundreds of variants associated with common diseases such as coronary artery disease (CAD), most disease-associated variants lie within non-coding regions of the genome, rendering it difficult to determine the downstream causal gene and cell type. Here, we performed paired single nucleus gene expression and chromatin accessibility profiling from 44 human coronary arteries. To link disease variants to molecular traits, we developed a meta-map of 88 samples and discovered 11,182 single-cell chromatin accessibility quantitative trait loci (caQTLs). Heritability enrichment analysis and disease variant mapping demonstrated that smooth muscle cells (SMCs) harbor the greatest genetic risk for CAD. To capture the continuum of SMC cell states in disease, we used dynamic single cell caQTL modeling for the first time in tissue to uncover QTLs whose effects are modified by cell state and expand our insight into genetic regulation of heterogenous cell populations. Notably, we identified a variant in the COL4A1/COL4A2 CAD GWAS locus which becomes a caQTL as SMCs de-differentiate by changing a transcription factor binding site for EGR1/2. To unbiasedly prioritize functional candidate genes, we built a genome-wide single cell variant to enhancer to gene (scV2E2G) map for human CAD to link disease variants to causal genes in cell types. Using this approach, we found several hundred genes predicted to be linked to disease variants in different cell types. Next, we performed genome-wide Hi-C in 16 human coronary arteries to build tissue specific maps of chromatin conformation and link disease variants to integrated chromatin hubs and distal target genes. Using this approach, we show that rs4887091 within the ADAMTS7 CAD GWAS locus modulates function of a super chromatin interactome through a change in a CTCF binding site. Finally, we used CRISPR interference to validate a distal gene, AMOTL2, liked to a CAD GWAS locus. Collectively we provide a disease-agnostic framework to translate human genetic findings to identify pathologic cell states and genes driving disease, producing a comprehensive scV2E2G map with genetic and tissue level convergence for future mechanistic and therapeutic studies.

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.

Fructose-1-kinase has pleiotropic roles in Escherichia coli

In Escherichia coli, the master transcription regulator catabolite repressor activator (Cra) regulates >100 genes in central metabolism. Cra binding to DNA is allosterically regulated by binding to fructose-1-phosphate (F-1-P), but the only documented source of F-1-P is from the concurrent import and phosphorylation of exogenous fructose. Thus, many have proposed that fructose-1,6-bisphosphate (F-1,6-BP) is also a physiological regulatory ligand. However, the role of F-1,6-BP has been widely debated. Here, we report that the E. coli enzyme fructose-1-kinase (FruK) can carry out its "reverse" reaction under physiological substrate concentrations to generate F-1-P from F-1,6-BP. We further show that FruK directly binds Cra with nanomolar affinity and forms higher order, heterocomplexes. Growth assays with a ΔfruK strain and fruK complementation show that FruK has a broader role in metabolism than fructose catabolism. Since fruK itself is repressed by Cra, these newly-reported events add layers to the dynamic regulation of E. coli's central metabolism that occur in response to changing nutrients. These findings might have wide-spread relevance to other γ-proteobacteria, which conserve both Cra and FruK.

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.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

Quantitative Determination of DNA Bridging Efficiency of Chromatin Proteins

DNA looping is important for genome organization in all domains of life. The basis of DNA loop formation is the bridging of two separate DNA double helices. Detecting DNA bridge formation generally involves the use of complex single-molecule techniques (atomic force microscopy, magnetic or optical tweezers). Although DNA bridging can be qualitatively described, quantification of DNA bridging and bridging dynamics using these techniques is challenging. Here we describe a biochemical assay capable of not only detecting DNA bridge formation but also allowing for quantification of DNA bridging efficiency and the quantification of the effects of physicochemical conditions or protein interaction partners on DNA bridge formation.

Fructose-1-kinase has pleiotropic roles in Escherichia coli

In Escherichia coli, the master transcription regulator Catabolite Repressor Activator (Cra) regulates >100 genes in central metabolism. Cra binding to DNA is allosterically regulated by binding to fructose-1-phosphate (F-1-P), but the only documented source of F-1-P is from the concurrent import and phosphorylation of exogenous fructose. Thus, many have proposed that fructose-1,6-bisphosphate (F-1,6-BP) is also a physiological regulatory ligand. However, the role of F-1,6-BP has been widely debated. Here, we report that the E. coli enzyme fructose-1-kinase (FruK) can carry out its "reverse" reaction under physiological substrate concentrations to generate F-1-P from F-1,6-BP. We further show that FruK directly binds Cra with nanomolar affinity and forms higher order, heterocomplexes. Growth assays with a ΔfruK strain and fruK complementation show that FruK has a broader role in metabolism than fructose catabolism. The ΔfruK strain also alters biofilm formation. Since fruK itself is repressed by Cra, these newly-reported events add layers to the dynamic regulation of E. coli central metabolism that occur in response to changing nutrients. These findings might have wide-spread relevance to other γ-proteobacteria, which conserve both Cra and FruK.

Multidimensional scaling methods can reconstruct genomic DNA loops using Hi-C data properties

This paper proposes multidimensional scaling (MDS) applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Currently, the mechanisms underlying the regulation of gene expression are poorly understood, and where and when DNA loops are formed remains undetermined. Previous studies have focused on reproducing the entire three-dimensional structure of chromatin; however, identifying DNA loops using these data is time-consuming and difficult. MDS is an unsupervised method for reconstructing the original coordinates from a distance matrix. Here, MDS was applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Hi-C data were converted to distances by taking the inverse to reproduce loops via MDS, and the missing values were set to zero. Using the converted data, MDS was applied to the log-transformed genomic coordinate distances and this process successfully reproduced the DNA loops in the given structure. Consequently, the reconstructed DNA loops revealed significantly more DNA-transcription factor interactions involved in DNA loop formation than those obtained from previously applied methods. Furthermore, the reconstructed DNA loops were significantly consistent with chromatin immunoprecipitation followed by sequencing (ChIP-seq) peak positions. In conclusion, the proposed method is an improvement over previous methods for identifying DNA loops.

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.

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.

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.

Sophisticated CRISPR/Cas tools for fine-tuning plant performance

Over the last years, the discovery of various natural and the development of a row of engineered CRISPR/Cas nucleases have made almost every site of plant genomes accessible for the induction of specific changes. Newly developed tools open up a wide range of possibilities for the induction of genetic variability, from changing a single bp to Mbps, and thus to fine-tune plant performance. Whereas early approaches focused on targeted mutagenesis, recently developed tools enable the induction of precise and predefined genomic modifications. The use of base editors allows the substitution of single nucleotides, whereas the use of prime editors and gene targeting methods enables the induction of larger sequence modifications from a few bases to several kbp. Recently, through CRISPR/Cas-mediated chromosome engineering, it became possible to induce heritable inversions and translocations in the Mbp range. Thus, a novel way of breaking and fixing genetic linkages has come into reach for breeders. In addition, sequence-specific recruitment of various factors involved in transcriptional and post-transcriptional regulation has been shown to provide an additional class of methods for the fine tuning of plant performance. In this review, we provide an overview of the most recent progress in the field of CRISPR/Cas-based tool development for plant genome engineering and try to evaluate the importance of these developments for breeding and biotechnological applications.

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.

CRISPR/dCas9 platforms in plants: strategies and applications beyond genome editing

Clustered regularly interspaced short palindromic repeat (CRISPR) and Cas9-associated protein systems provide a powerful genetic manipulation tool that can drive plant research forward. Nuclease-dead Cas9 (dCas9) is an enzymatically inactive mutant of Cas9 in which its endonuclease activity is non-functional. The applications of CRISPR/dCas9 have expanded and diversified in recent years. Originally, dCas9 was used as a CRISPR/Cas9 re-engineering tool that enables targeted expression of any gene or multiple genes through recruitment of transcriptional effector domains without introducing irreversible DNA-damaging mutations. Subsequent applications have made use of its ability to recruit modifying enzymes and reporter proteins to DNA target sites. In this paper, the most recent progress in the applications of CRISPR/dCas9 in plants, which include gene activation and repression, epigenome editing, modulation of chromatin topology, live-cell chromatin imaging and DNA-free genetic modification, will be reviewed. The associated strategies for exploiting the CRISPR/dCas9 system for crop improvement with a dimer of the future of the CRISPR/dCas9 system in the functional genomics of crops and the development of traits will be briefly discussed.

Homolog comparisons further reconcile in vitro and in vivo correlations of protein activities by revealing over-looked physiological factors

To bridge biological and biochemical disciplines, the relationship between in vitro protein biochemical function and in vivo activity must be established. Such studies can (a) help determine whether properties measured in simple, dilute solutions extrapolate to the complex in vivo conditions and (b) illuminate cryptic biological factors that are new avenues for study. We have explored the in vivo-in vitro relationship for chimeras built from LacI/GalR transcription regulators. In prior studies of individual chimeras, amino acid changes that altered in vitro DNA binding affinity exhibited correlated changes in in vivo transcription repression. However, discrepancies arose when the two datasets were compared to each other: Although their DNA binding domains were identical and their in vitro binding affinities spanned the same range, their in vivo repression ranges differed by >50-fold. Here, we determined that the presence of endogenous ligand for one chimera further exacerbated the offset, but that different abilities to simultaneously bind and "loop" two DNA operators resolves the discrepancy. Indeed, results suggest that the lac operon can be looped by even weakly interacting repressor dimers. For looping-competent repressors, we measured in vitro binding to the secondary operator. Surprisingly, this was largely insensitive to amino acid changes in the repressor protein, which reflects altered specificity; this supports the emerging view that the locations of specificity determining positions can be unique to each protein homolog. In aggregate, this work illustrates how a comparative approach can enrich understanding of the in vivo-in vitro relationship and suggest unexpected avenues for future study.

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.

YY1 Is a Structural Regulator of Enhancer-Promoter Loops

There is considerable evidence that chromosome structure plays important roles in gene control, but we have limited understanding of the proteins that contribute to structural interactions between gene promoters and their enhancer elements. Large DNA loops that encompass genes and their regulatory elements depend on CTCF-CTCF interactions, but most enhancer-promoter interactions do not employ this structural protein. Here, we show that the ubiquitously expressed transcription factor Yin Yang 1 (YY1) contributes to enhancer-promoter structural interactions in a manner analogous to DNA interactions mediated by CTCF. YY1 binds to active enhancers and promoter-proximal elements and forms dimers that facilitate the interaction of these DNA elements. Deletion of YY1 binding sites or depletion of YY1 protein disrupts enhancer-promoter looping and gene expression. We propose that YY1-mediated enhancer-promoter interactions are a general feature of mammalian gene control.

A cautionary note on the use of chromosome conformation capture in plants

BACKGROUND

The chromosome conformation capture (3C) technique is a method to study chromatin interactions at specific genomic loci. Initially established for yeast the 3C technique has been adapted to plants in recent years in order to study chromatin interactions and their role in transcriptional gene regulation. As the plant scientific community continues to implement this technology, a discussion on critical controls, validations steps and interpretation of 3C data is essential to fully benefit from 3C in plants.

RESULTS

Here we assess the reliability and robustness of the 3C technique for the detection of chromatin interactions in Arabidopsis. As a case study, we applied this methodology to the genomic locus of a floral integrator gene SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1), and demonstrate the need of several controls and standard validation steps to allow a meaningful interpretation of 3C data. The intricacies of this promising but challenging technique are discussed in depth.

CONCLUSIONS

The 3C technique offers an interesting opportunity to study chromatin interactions at a resolution infeasible by microscopy. However, for interpretation of 3C interaction data and identification of true interactions, 3C technology demands a stringent experimental setup and extreme caution.

Efficient chromosomal-scale DNA looping in Escherichia coli using multiple DNA-looping elements

Genes are frequently regulated by interactions between proteins that bind to the DNA near the gene and proteins that bind to DNA sites located far away, with the intervening DNA looped out. But it is not understood how efficient looping can occur when the sites are very far apart. We develop a simple theoretical framework that relates looping efficiency to the energetic cost and benefit of looping, allowing prediction of the efficiency of single or multiple nested loops at different distances. Measurements of absolute loop efficiencies for Lac repressor and λ CI using gene expression reporters in Escherichia coli cells show that, as predicted by the model, long-range DNA looping between a pair of sites can be strongly enhanced by the use of nested DNA loops or by the use of additional protein-binding sequences. A combination of these approaches was able to generate efficient DNA looping at a 200 kb distance.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

Proteins mediating DNA loops effectively block transcription

Loops are ubiquitous topological elements formed when proteins simultaneously bind to two noncontiguous DNA sites. While a loop‐mediating protein may regulate initiation at a promoter, the presence of the protein at the other site may be an obstacle for RNA polymerases (RNAP) transcribing a different gene. To test whether a DNA loop alters the extent to which a protein blocks transcription, the lac repressor (LacI) was used. The outcome of in vitro transcription along templates containing two LacI operators separated by 400 bp in the presence of LacI concentrations that produced both looped and unlooped molecules was visualized with scanning force microscopy (SFM). An analysis of transcription elongation complexes, moving for 60 s at an average of 10 nt/s on unlooped DNA templates, revealed that they more often surpassed LacI bound to the lower affinity O2 operator than to the highest affinity Os operator. However, this difference was abrogated in looped DNA molecules where LacI became a strong roadblock independently of the affinity of the operator. Recordings of transcription elongation complexes, using magnetic tweezers, confirmed that they halted for several minutes upon encountering a LacI bound to a single operator. The average pause lifetime is compatible with RNAP waiting for LacI dissociation, however, the LacI open conformation visualized in the SFM images also suggests that LacI could straddle RNAP to let it pass. Independently of the mechanism by which RNAP bypasses the LacI roadblock, the data indicate that an obstacle with looped topology more effectively interferes with transcription.

Biasing Simulations of DNA Base Pair Parameters with Application to Propellor Twisting in AT/AT, AA/TT, and AC/GT Steps and Their Uracil Analogs

An accurate and efficient implementation of the six DNA base pair parameters as order parameters for enhanced sampling simulations is presented. The parameter definitions are defined by vector algebra operations on a reduced atomic set of the base pair, and correlate very well with standard definitions. Application of the model is illustrated by umbrella sampling simulations of propeller twisting within AT/AT, AA/TT, and AC/GT steps and their uracil analogs. Strong correlations are found between propeller twisting and a number of conformational parameters, including buckle, opening, BI/BII backbone configuration, and sugar puckering. The thymine methyl group is observed to notably alter the local conformational free energy landscape, with effects within and directly upstream of the thymine containing base pair.

Supercoiling Effects on Short-Range DNA Looping in E. coli

DNA-protein loops can be essential for gene regulation. The Escherichia coli lactose (lac) operon is controlled by DNA-protein loops that have been studied for decades. Here we adapt this model to test the hypothesis that negative superhelical strain facilitates the formation of short-range (6-8 DNA turns) repression loops in E. coli. The natural negative superhelicity of E. coli DNA is regulated by the interplay of gyrase and topoisomerase enzymes, adding or removing negative supercoils, respectively. Here, we measured quantitatively DNA looping in three different E. coli strains characterized by different levels of global supercoiling: wild type, gyrase mutant (gyrB226), and topoisomerase mutant (ΔtopA10). DNA looping in each strain was measured by assaying repression of the endogenous lac operon, and repression of ten reporter constructs with DNA loop sizes between 70-85 base pairs. Our data are most simply interpreted as supporting the hypothesis that negative supercoiling facilitates gene repression by small DNA-protein loops in living bacteria.

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.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

Organization of DNA in a bacterial nucleoid

BACKGROUND

It is unclear how DNA is packaged in a bacterial cell in the absence of nucleosomes. To investigate the initial level of DNA condensation in bacterial nucleoid we used in vivo DNA digestion coupled with high-throughput sequencing of the digestion-resistant fragments. To this end, we transformed E. coli cells with a plasmid expressing micrococcal nuclease. The nuclease expression was under the control of AraC repressor, which enabled us to perform an inducible digestion of bacterial nucleoid inside a living cell.

RESULTS

Analysis of the genomic localization of the digestion-resistant fragments revealed their non-random distribution. The patterns observed in the distribution of the sequenced fragments indicate the presence of short DNA segments protected from the enzyme digestion, possibly because of interaction with DNA-binding proteins. The average length of such digestion-resistant segments is about 50 bp and the characteristic repeat in their distribution is about 90 bp. The gene starts are depleted of the digestion-resistant fragments, suggesting that these genomic regions are more exposed than genomic sequences on average. Sequence analysis of the digestion-resistant segments showed that while the GC-content of such sequences is close to the genome-wide value, they are depleted of A-tracts as compared to the bulk genomic DNA or to the randomized sequence of the same nucleotide composition.

CONCLUSIONS

Our results suggest that DNA is packaged in the bacterial nucleoid in a non-random way that facilitates interaction of the DNA binding factors with regulatory regions of the genome.

Simulation Assisted Analysis of the Intrinsic Stiffness for Short DNA Molecules Imaged with Scanning Atomic Force Microscopy

Studying the mechanical properties of short segments of dsDNA can provide insight into various biophysical phenomena, from DNA looping to the organization of nucleosomes. Scanning atomic force microscopy (AFM) is able to acquire images of single DNA molecules with near-basepair resolution. From many images, one may use equilibrium statistical mechanics to quantify the intrinsic stiffness (or persistence length) of the DNA. However, this approach is highly dependent upon both the correct microscopic polymer model and a correct image analysis of DNA contours. These complications have led to significant debate over the flexibility of dsDNA at short length scales. We first show how to extract accurate measures of DNA contour lengths by calibrating to DNA traces of simulated AFM data. After this calibration, we show that DNA adsorbed on an aminopropyl-mica surface behaves as a worm-like chain (WLC) for contour lengths as small as ~20 nm. We also show that a DNA binding protein can modify the mechanics of the DNA from that of a WLC.

Chromatin and Genomic determinants of alternative splicing

Alternative splicing significantly contributes to proteomic diversity and mis-regulation of splicing can cause diseases in human. Although both genomic and chromatin features have been shown to associate with splicing, the mechanisms by which various chromatin marks influence splicing is not clear for the most part. Moreover, it is not known whether the influence of specific genomic features on splicing is potentially modulated by the chromatin context. Here we report a deep neural network (DNN) model for predicting exon inclusion based on comprehensive genomic and chromatin features. Our analysis in three cell lines shows that, while both genomic and chromatin features can predict splicing to varying degrees, genomic features are the primary drivers of splicing, and the predictive power of chromatin features can largely be explained by their correlation with genomic features; chromatin features do not yield substantial independent contribution to splicing predictability. However, our model identified specific interactions between chromatin and genomic features suggesting that the effect of genomic elements may be modulated by chromatin context.

The dam replacing gene product enhances Neisseria gonorrhoeae FA1090 viability and biofilm formation

Many Neisseriaceae do not exhibit Dam methyltransferase activity and, instead of the dam gene, possess drg (dam replacing gene) inserted in the leuS/dam locus. The drg locus in Neisseria gonorrhoeae FA1090 has a lower GC-pairs content (40.5%) compared to the whole genome of N. gonorrhoeae FA1090 (52%). The gonococcal drg gene encodes a DNA endonuclease Drg, with GmeATC specificity. Disruption of drg or insertion of the dam gene in gonococcal genome changes the level of expression of genes as shown by transcriptome analysis. For the drg-deficient N. gonorrhoeae mutant, a total of 195 (8.94% of the total gene pool) genes exhibited an altered expression compared to the wt strain by at least 1.5 fold. In dam-expressing N. gonorrhoeae mutant, the expression of 240 genes (11% of total genes) was deregulated. Most of these deregulated genes were involved in translation, DNA repair, membrane biogenesis and energy production as shown by cluster of orthologous group analysis. In vivo, the inactivation of drg gene causes the decrease of the number of live neisserial cells and long lag phase of growth. The insertion of dam gene instead of drg locus restores cell viability. We have also shown that presence of the drg gene product is important for N. gonorrhoeae FA1090 in adhesion, including human epithelial cells, and biofilm formation. Biofilm produced by drg-deficient strain is formed by more dispersed cells, compared to this one formed by parental strain as shown by scanning electron and confocal microscopy. Also adherence assays show a significantly smaller biomass of formed biofilm (OD₅₇₀ = 0.242 ± 0.038) for drg-deficient strain, compared to wild-type strain (OD₅₇₀ = 0.378 ± 0.057). Dam-expressing gonococcal cells produce slightly weaker biofilm with cells embedded in an extracellular matrix. This strain has also a five times reduced ability for adhesion to human epithelial cells. In this context, the presence of Drg is more advantageous for N. gonorrhoeae biology than Dam presence.

A complex genetic switch involving overlapping divergent promoters and DNA looping regulates expression of conjugation genes of a gram-positive plasmid

Plasmid conjugation plays a significant role in the dissemination of antibiotic resistance and pathogenicity determinants. Understanding how conjugation is regulated is important to gain insights into these features. Little is known about regulation of conjugation systems present on plasmids from Gram-positive bacteria. pLS20 is a native conjugative plasmid from the Gram-positive bacterium Bacillus subtilis. Recently the key players that repress and activate pLS20 conjugation have been identified. Here we studied in detail the molecular mechanism regulating the pLS20 conjugation genes using both in vivo and in vitro approaches. Our results show that conjugation is subject to the control of a complex genetic switch where at least three levels of regulation are integrated. The first of the three layers involves overlapping divergent promoters of different strengths regulating expression of the conjugation genes and the key transcriptional regulator Rco_LS20. The second layer involves a triple function of Rco_LS20 being a repressor of the main conjugation promoter and an activator and repressor of its own promoter at low and high concentrations, respectively. The third level of regulation concerns formation of a DNA loop mediated by simultaneous binding of tetrameric Rco_LS20 to two operators, one of which overlaps with the divergent promoters. The combination of these three layers of regulation in the same switch allows the main conjugation promoter to be tightly repressed during conditions unfavorable to conjugation while maintaining the sensitivity to accurately switch on the conjugation genes when appropriate conditions occur. The implications of the regulatory switch and comparison with other genetic switches involving DNA looping are discussed.

Plasmids are extrachromosomal, autonomously replicating units that are harbored by many bacteria. Many plasmids encode transfer function allowing them to be transferred into plasmid-free bacteria by a process named conjugation. Since many of them also carry antibiotic resistance genes, plasmid-mediated conjugation is a major mechanism in the dissemination of antibiotic resistance. In depth knowledge on the regulation of conjugation genes is a prerequisite to design measures interfering with the spread of antibiotic resistance. pLS20 is a conjugative plasmid of the soil bacterium Bacillus subtilis, which is also a gut commensal in animals and humans. Here we describe in detail the molecular mechanism by which the key transcriptional regulator tightly represses the conjugation genes during conditions unfavorable to conjugation without compromising the ability to switch on accurately the conjugation genes when appropriate. We found that conjugation is subject to the control of a unique genetic switch where at least three levels of regulation are integrated. The first level involves overlapping divergent promoters of different strengths. The second layer involves a triple function of the transcriptional regulator. And the third level of regulation concerns formation of a DNA loop mediated by the transcriptional regulator.

Multiple LacI-mediated loops revealed by Bayesian statistics and tethered particle motion

The bacterial transcription factor LacI loops DNA by binding to two separate locations on the DNA simultaneously. Despite being one of the best-studied model systems for transcriptional regulation, the number and conformations of loop structures accessible to LacI remain unclear, though the importance of multiple coexisting loops has been implicated in interactions between LacI and other cellular regulators of gene expression. To probe this issue, we have developed a new analysis method for tethered particle motion, a versatile and commonly used in vitro single-molecule technique. Our method, vbTPM, performs variational Bayesian inference in hidden Markov models. It learns the number of distinct states (i.e. DNA–protein conformations) directly from tethered particle motion data with better resolution than existing methods, while easily correcting for common experimental artifacts. Studying short (roughly 100 bp) LacI-mediated loops, we provide evidence for three distinct loop structures, more than previously reported in single-molecule studies. Moreover, our results confirm that changes in LacI conformation and DNA-binding topology both contribute to the repertoire of LacI-mediated loops formed in vitro, and provide qualitatively new input for models of looping and transcriptional regulation. We expect vbTPM to be broadly useful for probing complex protein–nucleic acid interactions.

Regulation of human CYP2C9 expression by electrophilic stress involves activator protein 1 activation and DNA looping

Cytochrome P450 (CYP)2C9 and CYP2C19 are important human enzymes that metabolize therapeutic drugs, environmental chemicals, and physiologically important endogenous compounds. Initial studies using primary human hepatocytes showed induction of both the CYP2C9 and CYP2C19 genes by tert-butylhydroquinone (tBHQ). As a pro-oxidant, tBHQ regulates the expression of cytoprotective genes by activation of redox-sensing transcription factors, such as the nuclear factor E2-related factor 2 (Nrf2) and members of the activator protein 1 (AP-1) family of proteins. The promoter region of CYP2C9 contains two putative AP-1 sites (TGAGTCA) at positions -2201 and -1930, which are also highly conserved in CYP2C19. The CYP2C9 promoter is activated by ectopic expression of cFos and JunD, whereas Nrf2 had no effect. Using specific kinase inhibitors for mitogen-activated protein kinase, we showed that extracellular signal-regulated kinase and Jun N-terminal kinase are essential for tBHQ-induced expression of CYP2C9. Electrophoretic mobility shift assays demonstrate that cFos distinctly interacts with the distal AP-1 site and JunD with the proximal site. Because cFos regulates target genes as heterodimers with Jun proteins, we hypothesized that DNA looping might be required to bring the distal and proximal AP-1 sites together to activate the CYP2C9 promoter. Chromosome conformation capture analyses confirmed the formation of a DNA loop in the CYP2C9 promoter, possibly allowing interaction between cFos at the distal site and JunD at the proximal site to activate CYP2C9 transcription in response to electrophiles. These results indicate that oxidative stress generated by exposure to electrophilic xenobiotics and metabolites induces the expression of CYP2C9 and CYP2C19 in human hepatocytes.

Lac repressor mediated DNA looping: Monte Carlo simulation of constrained DNA molecules complemented with current experimental results

Tethered particle motion (TPM) experiments can be used to detect time-resolved loop formation in a single DNA molecule by measuring changes in the length of a DNA tether. Interpretation of such experiments is greatly aided by computer simulations of DNA looping which allow one to analyze the structure of the looped DNA and estimate DNA-protein binding constants specific for the loop formation process. We here present a new Monte Carlo scheme for accurate simulation of DNA configurations subject to geometric constraints and apply this method to Lac repressor mediated DNA looping, comparing the simulation results with new experimental data obtained by the TPM technique. Our simulations, taking into account the details of attachment of DNA ends and fluctuations of the looped subsegment of the DNA, reveal the origin of the double-peaked distribution of RMS values observed by TPM experiments by showing that the average RMS value for anti-parallel loop types is smaller than that of parallel loop types. The simulations also reveal that the looping probabilities for the anti-parallel loop types are significantly higher than those of the parallel loop types, even for loops of length 600 and 900 base pairs, and that the correct proportion between the heights of the peaks in the distribution can only be attained when loops with flexible Lac repressor conformation are taken into account. Comparison of the in silico and in vitro results yields estimates for the dissociation constants characterizing the binding affinity between O1 and Oid DNA operators and the dimeric arms of the Lac repressor.

Quantitation of the DNA tethering effect in long-range DNA looping in vivo and in vitro using the Lac and λ repressors

Efficient and specific interactions between proteins bound to the same DNA molecule can be dependent on the length of the DNA tether that connects them. Measurement of the strength of this DNA tethering effect has been largely confined to short separations between sites, and it is not clear how it contributes to long-range DNA looping interactions, such as occur over separations of tens to hundreds of kilobase pairs in vivo. Here, gene regulation experiments using the LacI and λ CI repressors, combined with mathematical modeling, were used to quantitate DNA tethering inside Escherichia coli cells over the 250- to 10,000-bp range. Although LacI and CI loop DNA in distinct ways, measurements of the tethering effect were very similar for both proteins. Tethering strength decreased with increasing separation, but even at 5- to 10-kb distances, was able to increase contact probability 10- to 20-fold and drive efficient looping. Tethering in vitro with the Lac repressor was measured for the same 600-to 3,200-bp DNAs using tethered particle motion, a single molecule technique, and was 5- to 45-fold weaker than in vivo over this range. Thus, the enhancement of looping seen previously in vivo at separations below 500 bp extends to large separations, underlining the need to understand how in vivo factors aid DNA looping. Our analysis also suggests how efficient and specific looping could be achieved over very long DNA separations, such as what occurs between enhancers and promoters in eukaryotic cells.

Mechanosensing of DNA bending in a single specific protein-DNA complex

Many crucial biological processes are regulated by mechanical stimuli. Here, we report new findings that pico-Newton forces can drastically affect the stability of the site-specific DNA binding of a single transcription factor, the E. coli integration host factor (IHF), by stretching a short ~150 nm DNA containing a single IHF binding site. Dynamic binding and unbinding of single IHF were recorded and analyzed for the force-dependent stability of the IHF-DNA complex. Our results demonstrate that the IHF-DNA interaction is fine tuned by force in different salt concentration and temperature over physiological ranges, indicating that, besides other physiological factors, force may play equally important role in transcription regulation. These findings have broad implications with regard to general mechanosensitivity of site-specific DNA bending proteins.

Reliable prediction of complex phenotypes from a modular design in free energy space: an extensive exploration of the lac operon

The basic methodology for designing, altering, and constructing biological systems is increasingly relying on well-established engineering principles to move forward from trial and error approaches to reliably predicting the system behavior from the properties of the components and their interactions. The inherent complexity of even the simplest biological systems, however, often precludes achieving such predictive power. A prototypical example is the lac operon, one of the best-characterized genetic systems, which still poses serious challenges for understanding the results of combining its parts into novel setups. The reason is the pervasive complex hierarchy of events involved in gene regulation that extend from specific protein-DNA interactions to the combinatorial assembly of nucleoprotein complexes. Here, we integrate such complexity into a few-parameter model to accurately predict gene expression from a few simple rules to connect the parts. The model accurately reproduces the observed transcriptional activity of the lac operon over a 10,000-fold range for 21 different operator setups, different repressor concentrations, and tetrameric and dimeric forms of the repressor. Incorporation of the calibrated model into more complex scenarios accurately captures the induction curves for key operator configurations and the temporal evolution of the β-galactosidase activity of cell populations.

Quantitative influence of macromolecular crowding on gene regulation kinetics

We introduce macromolecular crowding quantitatively into the model for kinetics of gene regulation in Escherichia coli. We analyse and compute the specific-site searching time for 180 known transcription factors (TFs) regulating 1300 operons. The time is between 160 s (e.g. for SoxS M_w = 12.91 kDa) and 1550 s (e.g. for PepA₆ of M_w = 329.28 kDa). Diffusion coefficients for one-dimensional sliding are between for large proteins up to for small monomers or dimers. Three-dimensional diffusion coefficients in the cytoplasm are 2 orders of magnitude larger than 1D sliding coefficients, nevertheless the sliding enhances the binding rates of TF to specific sites by 1–2 orders of magnitude. The latter effect is due to ubiquitous non-specific binding. We compare the model to experimental data for LacI repressor and find that non-specific binding of the protein to DNA is activation- and not diffusion-limited. We show that the target location rate by LacI repressor is optimized with respect to microscopic rate constant for association to non-specific sites on DNA. We analyse the effect of oligomerization of TFs and DNA looping effects on searching kinetics. We show that optimal searching strategy depends on TF abundance.

Optical Methods to Study Protein-DNA Interactions in Vitro and in Living Cells at the Single-Molecule Level

The maintenance of intact genetic information, as well as the deployment of transcription for specific sets of genes, critically rely on a family of proteins interacting with DNA and recognizing specific sequences or features. The mechanisms by which these proteins search for target DNA are the subject of intense investigations employing a variety of methods in biology. A large interest in these processes stems from the faster-than-diffusion association rates, explained in current models by a combination of 3D and 1D diffusion. Here, we present a review of the single-molecule approaches at the forefront of the study of protein-DNA interaction dynamics and target search in vitro and in vivo. Flow stretch, optical and magnetic manipulation, single fluorophore detection and localization as well as combinations of different methods are described and the results obtained with these techniques are discussed in the framework of the current facilitated diffusion model.

DNA looping in prokaryotes: experimental and theoretical approaches

Transcriptional regulation is at the heart of biological functions such as adaptation to a changing environment or to new carbon sources. One of the mechanisms which has been found to modulate transcription, either positively (activation) or negatively (repression), involves the formation of DNA loops. A DNA loop occurs when a protein or a complex of proteins simultaneously binds to two different sites on DNA with looping out of the intervening DNA. This simple mechanism is central to the regulation of several operons in the genome of the bacterium Escherichia coli, like the lac operon, one of the paradigms of genetic regulation. The aim of this review is to gather and discuss concepts and ideas from experimental biology and theoretical physics concerning DNA looping in genetic regulation. We first describe experimental techniques designed to show the formation of a DNA loop. We then present the benefits that can or could be derived from a mechanism involving DNA looping. Some of these are already experimentally proven, but others are theoretical predictions and merit experimental investigation. Then, we try to identify other genetic systems that could be regulated by a DNA looping mechanism in the genome of Escherichia coli. We found many operons that, according to our set of criteria, have a good chance to be regulated with a DNA loop. Finally, we discuss the proposition recently made by both biologists and physicists that this mechanism could also act at the genomic scale and play a crucial role in the spatial organization of genomes.

Novel insights from hybrid LacI/GalR proteins: family-wide functional attributes and biologically significant variation in transcription repression

LacI/GalR transcription regulators have extensive, non-conserved interfaces between their regulatory domains and the 18 amino acids that serve as ‘linkers’ to their DNA-binding domains. These non-conserved interfaces might contribute to functional differences between paralogs. Previously, two chimeras created by domain recombination displayed novel functional properties. Here, we present a synthetic protein family, which was created by joining the LacI DNA-binding domain/linker to seven additional regulatory domains. Despite ‘mismatched’ interfaces, chimeras maintained allosteric response to their cognate effectors. Therefore, allostery in many LacI/GalR proteins does not require interfaces with precisely matched interactions. Nevertheless, the chimeric interfaces were not silent to mutagenesis, and preliminary comparisons suggest that the chimeras provide an ideal context for systematically exploring functional contributions of non-conserved positions. DNA looping experiments revealed higher order (dimer–dimer) oligomerization in several chimeras, which might be possible for the natural paralogs. Finally, the biological significance of repression differences was determined by measuring bacterial growth rates on lactose minimal media. Unexpectedly, moderate and strong repressors showed an apparent induction phase, even though inducers were not provided; therefore, an unknown mechanism might contribute to regulation of the lac operon. Nevertheless, altered growth correlated with altered repression, which indicates that observed functional modifications are significant.

Helicobacter pylori oriC--the first bipartite origin of chromosome replication in Gram-negative bacteria

Binding of the DnaA protein to oriC leads to DNA melting within the DNA unwinding element (DUE) and initiates replication of the bacterial chromosome. Helicobacter pylori oriC was previously identified as a region localized upstream of dnaA and containing a cluster of DnaA boxes bound by DnaA protein with a high affinity. However, no unwinding within the oriC sequence has been detected. Comprehensive in silico analysis presented in this work allowed us to identify an additional region (oriC2), separated from the original one (oriC1) by the dnaA gene. DnaA specifically binds both regions, but DnaA-dependent DNA unwinding occurs only within oriC2. Surprisingly, oriC2 is bound exclusively as supercoiled DNA, which directly shows the importance of the DNA topology in DnaA-oriC interactions, similarly as previously presented only for initiator-origin interactions in Archaea and some Eukaryota. We conclude that H. pylori oriC exhibits bipartite structure, being the first such origin discovered in a Gram-negative bacterium. The H. pylori mode of initiator-oriC interactions, with the loop formation between the subcomplexes of the discontinuous origin, resembles those discovered in Bacillus subtilis chromosome and in many plasmids, which might suggest a similar way of controlling initiation of replication.