Sequence dependence of transcription factor-mediated DNA looping

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.

Coarse-grained modeling reveals the impact of supercoiling and loop length in DNA looping kinetics

Measurements of protein-mediated DNA looping reveal that in vivo conditions favor the formation of loops shorter than those that occur in vitro, yet the precise physical mechanisms underlying this shift remain unclear. To understand the extent to which in vivo supercoiling may explain these shifts, we develop a theoretical model based on coarse-grained molecular simulation and analytical transition state theory, enabling us to map out looping energetics and kinetics as a function of two key biophysical parameters: superhelical density and loop length. We show that loops on the scale of a persistence length respond to supercoiling over a much wider range of superhelical densities and to a larger extent than longer loops. This effect arises from a tendency for loops to be centered on the plectonemic end region, which bends progressively more tightly with superhelical density. This trend reveals a mechanism by which supercoiling favors shorter loop lengths. In addition, our model predicts a complex kinetic response to supercoiling for a given loop length, governed by a competition between an enhanced rate of looping due to torsional buckling and a reduction in looping rate due to chain straightening as the plectoneme tightens at higher superhelical densities. Together, these effects lead to a flattening of the kinetic response to supercoiling within the physiological range for all but the shortest loops. Using experimental estimates for in vivo superhelical densities, we discuss our model's ability to explain available looping data, highlighting both the importance of supercoiling as a regulatory force in genetics and the additional complexities of looping phenomena in vivo.

OUP accepted manuscript

DNA mechanical properties play a critical role in every aspect of DNA-dependent biological processes. Recently a high throughput assay named loop-seq has been developed to quantify the intrinsic bendability of a massive number of DNA fragments simultaneously. Using the loop-seq data, we develop a software tool, DNAcycP, based on a deep-learning approach for intrinsic DNA cyclizability prediction. We demonstrate DNAcycP predicts intrinsic DNA cyclizability with high fidelity compared to the experimental data. Using an independent dataset from in vitro selection for enrichment of loopable sequences, we further verified the predicted cyclizability score, termed C-score, can well distinguish DNA fragments with different loopability. We applied DNAcycP to multiple species and compared the C-scores with available high-resolution chemical nucleosome maps. Our analyses showed that both yeast and mouse genomes share a conserved feature of high DNA bendability spanning nucleosome dyads. Additionally, we extended our analysis to transcription factor binding sites and surprisingly found that the cyclizability is substantially elevated at CTCF binding sites in the mouse genome. We further demonstrate this distinct mechanical property is conserved across mammalian species and is inherent to CTCF binding DNA motif.

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles' local geometry consists of straight segments connected by sharp bends compressing the DNA's inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles' inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

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.

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.

Characterization of Gene Repression by Designed Transcription Activator-like Effector Dimer Proteins

Gene regulation by control of transcription initiation is a fundamental property of living cells. Much of our understanding of gene repression originated from studies of the Escherichia coli lac operon switch, in which DNA looping plays an essential role. To validate and generalize principles from lac for practical applications, we previously described artificial DNA looping driven by designed transcription activator-like effector dimer (TALED) proteins. Because TALE monomers bind the idealized symmetrical lac operator sequence in two orientations, our prior studies detected repression due to multiple DNA loops. We now quantitatively characterize gene repression in living E. coli by a collection of individual TALED loops with systematic loop length variation. Fitting of a thermodynamic model allows unequivocal demonstration of looping and comparison of the engineered TALED repression system with the natural lac repressor system.

Deciphering the regulatory genome of Escherichia coli, one hundred promoters at a time

Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most well-studied of organisms, the bacterium Escherichia coli, for ≈65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, Reg-Seq, that links massively parallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than a E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

Sequence-dependent dynamics of synthetic and endogenous RSSs in V(D)J recombination

Developing lymphocytes of jawed vertebrates cleave and combine distinct gene segments to assemble antigen-receptor genes. This process called V(D)J recombination that involves the RAG recombinase binding and cutting recombination signal sequences (RSSs) composed of conserved heptamer and nonamer sequences flanking less well-conserved 12- or 23-bp spacers. Little quantitative information is known about the contributions of individual RSS positions over the course of the RAG-RSS interaction. We employ a single-molecule method known as tethered particle motion to track the formation, lifetime and cleavage of individual RAG-12RSS-23RSS paired complexes (PCs) for numerous synthetic and endogenous 12RSSs. We reveal that single-bp changes, including in the 12RSS spacer, can significantly and selectively alter PC formation or the probability of RAG-mediated cleavage in the PC. We find that some rarely used endogenous gene segments can be mapped directly to poor RAG binding on their adjacent 12RSSs. Finally, we find that while abrogating RSS nicking with Ca2+ leads to substantially shorter PC lifetimes, analysis of the complete lifetime distributions of any 12RSS even on this reduced system reveals that the process of exiting the PC involves unidentified molecular details whose involvement in RAG-RSS dynamics are crucial to quantitatively capture kinetics in V(D)J recombination.

Protein-mediated looping of DNA under tension requires supercoiling

Protein-mediated DNA looping is ubiquitous in chromatin organization and gene regulation, but to what extent supercoiling or nucleoid associated proteins promote looping is poorly understood. Using the lac repressor (LacI), a paradigmatic loop-mediating protein, we measured LacI-induced looping as a function of either supercoiling or the concentration of the HU protein, an abundant nucleoid protein in Escherichia coli. Negative supercoiling to physiological levels with magnetic tweezers easily drove the looping probability from 0 to 100% in single DNA molecules under slight tension that likely exists in vivo. In contrast, even saturating (micromolar) concentrations of HU could not raise the looping probability above 30% in similarly stretched DNA or 80% in DNA without tension. Negative supercoiling is required to induce significant looping of DNA under any appreciable tension.

Tethered multifluorophore motion reveals equilibrium transition kinetics of single DNA double helices

Understanding cellular functions and dysfunctions often begins with quantifying the interactions between the binding partners involved in the processes. Learning about the kinetics of the interactions is of particular importance to understand the dynamics of cellular processes. We created a tethered multifluorophore motion assay using DNA origami that enables over 1-hour-long recordings of the statistical binding and unbinding of single pairs of biomolecules directly in equilibrium. The experimental concept is simple and the data interpretation is very direct, which makes the system easy to use for a wide variety of researchers. Due to the modularity and addressability of the DNA origami-based assay, our system may be readily adapted to study various other molecular interactions.

We describe a tethered multifluorophore motion assay based on DNA origami for revealing bimolecular reaction kinetics on the single-molecule level. Molecular binding partners may be placed at user-defined positions and in user-defined stoichiometry; and binding states are read out by tracking the motion of quickly diffusing fluorescent reporter units. Multiple dyes per reporter unit enable singe-particle observation for more than 1 hour. We applied the system to study in equilibrium reversible hybridization and dissociation of complementary DNA single strands as a function of tether length, cation concentration, and sequence. We observed up to hundreds of hybridization and dissociation events per single reactant pair and could produce cumulative statistics with tens of thousands of binding and unbinding events. Because the binding partners per particle do not exchange, we could also detect subtle heterogeneity from molecule to molecule, which enabled separating data reflecting the actual target strand pair binding kinetics from falsifying influences stemming from chemically truncated oligonucleotides. Our data reflected that mainly DNA strand hybridization, but not strand dissociation, is affected by cation concentration, in agreement with previous results from different assays. We studied 8-bp-long DNA duplexes with virtually identical thermodynamic stability, but different sequences, and observed strongly differing hybridization kinetics. Complementary full-atom molecular-dynamics simulations indicated two opposing sequence-dependent phenomena: helical templating in purine-rich single strands and secondary structures. These two effects can increase or decrease, respectively, the fraction of strand collisions leading to successful nucleation events for duplex formation.

Tethered Particle Motion: An Easy Technique for Probing DNA Topology and Interactions with Transcription Factors

Tethered Particle Motion (TPM) is a versatile in vitro technique for monitoring the conformations a linear macromolecule, such as DNA, can exhibit. The technique involves monitoring the diffusive motion of a particle anchored to a fixed point via the macromolecule of interest, which acts as a tether. In this chapter, we provide an overview of TPM, review the fundamental principles that determine the accuracy with which effective tether lengths can be used to distinguish different tether conformations, present software tools that assist in capturing and analyzing TPM data, and provide a protocol which uses TPM to characterize lac repressor-induced DNA looping. Critical to any TPM assay is the understanding of the timescale over which the diffusive motion of the particle must be observed to accurately distinguish tether conformations. Approximating the tether as a Hookean spring, we show how to estimate the diffusion timescale and discuss how it relates to the confidence with which tether conformations can be distinguished. Applying those estimates to a lac repressor titration assay, we describe how to perform a TPM experiment. We also provide graphically driven software which can be used to speed up data collection and analysis. Lastly, we detail how TPM data from the titration assay can be used to calculate relevant molecular descriptors such as the J factor for DNA looping and lac repressor-operator dissociation constants. While the included protocol is geared toward studying DNA looping, the technique, fundamental principles, and analytical methods are more general and can be adapted to a wide variety of molecular systems.

Programmable DNA looping using engineered bivalent dCas9 complexes

DNA looping is a ubiquitous and critical feature of gene regulation. Although DNA looping can be efficiently detected, tools to readily manipulate DNA looping are limited. Here we develop CRISPR-based DNA looping reagents for creation of programmable DNA loops. Cleavage-defective Cas9 proteins of different specificity are linked by heterodimerization or translational fusion to create bivalent complexes able to link two separate DNA regions. After model-directed optimization, the reagents are validated using a quantitative DNA looping assay in E. coli. Looping efficiency is ~15% for a 4.7 kb loop, but is significantly improved by loop multiplexing with additional guides. Bivalent dCas9 complexes are also used to activate endogenous norVW genes by rewiring chromosomal DNA to bring distal enhancer elements to the gene promoters. Such reagents should allow manipulation of DNA looping in a variety of cell types, aiding understanding of endogenous loops and enabling creation of new regulatory connections.

DNA loops are a ubiquitious feature of gene regulation across the kingdoms of life. Here the authors design a Cas9-based dimerization system for inducing DNA loops in E. coli, allowing activation and rewiring of gene expression.

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.

Insights into Genome Architecture Deduced from the Properties of Short Lac Repressor-mediated DNA Loops

Genomic DNA is vastly longer than the space allotted to it in a cell. The molecule must fold with a level of organization that satisfies the imposed spatial constraints as well as allows for the processing of genetic information. Key players in this organization include the negative supercoiling of DNA, which facilitates the unwinding of the double-helical molecule, and the associations of DNA with proteins, which partition the DNA into isolated loops, or domains. In order to gain insight into the principles of genome organization and to visualize the folding of spatially constrained DNA, we have developed new computational methods to identify the preferred three-dimensional pathways of protein-mediated DNA loops and to characterize the topological properties of these structures. Here we focus on the levels of supercoiling and the spatial arrangements of DNA in model nucleoprotein systems with two topological domains. We construct these systems by anchoring DNA loops in opposing orientations on a common protein-DNA assembly, namely the Lac repressor protein with two bound DNA operators. The linked pieces of DNA form a covalently closed circle such that the protein attaches to two widely spaced sites along the DNA. We examine the effects of operator spacing, loop orientation, and long-range contacts on overall chain configuration and topology and discuss our findings in the context of classic experiments on the effects of supercoiling and operator spacing on Lac repressor-mediated looping and recent work on the role of proteins as barriers that divide genomes into independent topological domains.

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.

Single-molecule fluorescence studies on DNA looping

Structure and dynamics of DNA impact how the genetic code is processed and maintained. In addition to its biological importance, DNA has been utilized as building blocks of various nanomachines and nanostructures. Thus, understanding the physical properties of DNA is of fundamental importance to basic sciences and engineering applications. DNA can undergo various physical changes. Among them, DNA looping is unique in that it can bring two distal sites together, and thus can be used to mediate interactions over long distances. In this paper, we introduce a FRET-based experimental tool to study DNA looping at the single molecule level. We explain the connection between experimental measurables and a theoretical concept known as the J factor with the intent of raising awareness of subtle theoretical details that should be considered when drawing conclusions. We also explore DNA looping-assisted protein diffusion mechanism called intersegmental transfer using protein induced fluorescence enhancement (PIFE). We present some preliminary results and future outlooks.

Interplay of Protein Binding Interactions, DNA Mechanics, and Entropy in DNA Looping Kinetics

DNA looping plays a key role in many fundamental biological processes, including gene regulation, recombination, and chromosomal organization. The looping of DNA is often mediated by proteins whose structural features and physical interactions can alter the length scale at which the looping occurs. Looping and unlooping processes are controlled by thermodynamic contributions associated with mechanical deformation of the DNA strand and entropy arising from thermal fluctuations of the conformation. To determine how these confounding effects influence DNA looping and unlooping kinetics, we present a theoretical model that incorporates the role of the protein interactions, DNA mechanics, and conformational entropy. We show that for shorter DNA strands the interaction distance affects the transition state, resulting in a complex relationship between the looped and unlooped state lifetimes and the physical properties of the looped DNA. We explore the range of behaviors that arise with varying interaction distance and DNA length. These results demonstrate how DNA deformation and entropy dictate the scaling of the looping and unlooping kinetics versus the J-factor, establishing the connection between kinetic and equilibrium behaviors. Our results show how the twist-and-bend elasticity of the DNA chain modulates the kinetics and how the influence of the interaction distance fades away at intermediate to longer chain lengths, in agreement with previous scaling predictions.

Distinct structural features of TFAM drive mitochondrial DNA packaging versus transcriptional activation

TFAM (transcription factor A, mitochondrial) is a DNA-binding protein that activates transcription at the two major promoters of mitochondrial DNA (mtDNA)--the light strand promoter (LSP) and the heavy strand promoter 1 (HSP1). Equally important, it coats and packages the mitochondrial genome. TFAM has been shown to impose a U-turn on LSP DNA; however, whether this distortion is relevant at other sites is unknown. Here we present crystal structures of TFAM bound to HSP1 and to nonspecific DNA. In both, TFAM similarly distorts the DNA into a U-turn. Yet, TFAM binds to HSP1 in the opposite orientation from LSP explaining why transcription from LSP requires DNA bending, whereas transcription at HSP1 does not. Moreover, the crystal structures reveal dimerization of DNA-bound TFAM. This dimerization is dispensable for DNA bending and transcriptional activation but is important in DNA compaction. We propose that TFAM dimerization enhances mitochondrial DNA compaction by promoting looping of the DNA.

DNA, flexibly flexible

Investigators have constructed dsDNA molecules with several different base modifications and have characterized their bending and twisting flexibilities using atomic force microscopy, DNA ring closure, and single-molecule force spectroscopy with optical tweezers. The three methods provide persistence length measurements that agree semiquantitatively, and they show that the persistence length is surprisingly similar for all of the modified DNAs. The circular dichroism spectra of modified DNAs differ substantially. Simple explanations based on base stacking strength, polymer charge, or groove occupancy by functional groups cannot explain the results, which will guide further high-resolution theory and experiments.

Single-molecule analysis of RAG-mediated V(D)J DNA cleavage

The recombination-activating gene products, RAG1 and RAG2, initiate V(D)J recombination during lymphocyte development by cleaving DNA adjacent to conserved recombination signal sequences (RSSs). The reaction involves DNA binding, synapsis, and cleavage at two RSSs located on the same DNA molecule and results in the assembly of antigen receptor genes. We have developed single-molecule assays to examine RSS binding by RAG1/2 and their cofactor high-mobility group-box protein 1 (HMGB1) as they proceed through the steps of this reaction. These assays allowed us to observe in real time the individual molecular events of RAG-mediated cleavage. As a result, we are able to measure the binding statistics (dwell times) and binding energies of the initial RAG binding events and characterize synapse formation at the single-molecule level, yielding insights into the distribution of dwell times in the paired complex and the propensity for cleavage on forming the synapse. Interestingly, we find that the synaptic complex has a mean lifetime of roughly 400 s and that its formation is readily reversible, with only ∼40% of observed synapses resulting in cleavage at consensus RSS binding sites.

Modulation of DNA loop lifetimes by the free energy of loop formation

Storage and retrieval of the genetic information in cells is a dynamic process that requires the DNA to undergo dramatic structural rearrangements. DNA looping is a prominent example of such a structural rearrangement that is essential for transcriptional regulation in both prokaryotes and eukaryotes, and the speed of such regulations affects the fitness of individuals. Here, we examine the in vitro looping dynamics of the classic Lac repressor gene-regulatory motif. We show that both loop association and loop dissociation at the DNA-repressor junctions depend on the elastic deformation of the DNA and protein, and that both looping and unlooping rates approximately scale with the looping J factor, which reflects the system's deformation free energy. We explain this observation by transition state theory and model the DNA-protein complex as an effective worm-like chain with twist. We introduce a finite protein-DNA binding interaction length, in competition with the characteristic DNA deformation length scale, as the physical origin of the previously unidentified loop dissociation dynamics observed here, and discuss the robustness of this behavior to perturbations in several polymer parameters.

DNA topology confers sequence specificity to nonspecific architectural proteins

Topological constraints placed on short fragments of DNA change the disorder found in chain molecules randomly decorated by nonspecific, architectural proteins into tightly organized 3D structures. The bacterial heat-unstable (HU) protein builds up, counter to expectations, in greater quantities and at particular sites along simulated DNA minicircles and loops. Moreover, the placement of HU along loops with the "wild-type" spacing found in the Escherichia coli lactose (lac) and galactose (gal) operons precludes access to key recognition elements on DNA. The HU protein introduces a unique spatial pathway in the DNA upon closure. The many ways in which the protein induces nearly the same closed circular configuration point to the statistical advantage of its nonspecificity. The rotational settings imposed on DNA by the repressor proteins, by contrast, introduce sequential specificity in HU placement, with the nonspecific protein accumulating at particular loci on the constrained duplex. Thus, an architectural protein with no discernible DNA sequence-recognizing features becomes site-specific and potentially assumes a functional role upon loop formation. The locations of HU on the closed DNA reflect long-range mechanical correlations. The protein responds to DNA shape and deformability—the stiff, naturally straight double-helical structure—rather than to the unique features of the constituent base pairs. The structures of the simulated loops suggest that HU architecture, like nucleosomal architecture, which modulates the ability of regulatory proteins to recognize their binding sites in the context of chromatin, may influence repressor-operator interactions in the context of the bacterial nucleoid.

What controls DNA looping?

The looping of DNA provides a means of communication between sequentially distant genomic sites that operate in tandem to express, copy, and repair the information encoded in the DNA base sequence. The short loops implicated in the expression of bacterial genes suggest that molecular factors other than the naturally stiff double helix are involved in bringing the interacting sites into close spatial proximity. New computational techniques that take direct account of the three-dimensional structures and fluctuations of protein and DNA allow us to examine the likely means of enhancing such communication. Here, we describe the application of these approaches to the looping of a 92 base-pair DNA segment between the headpieces of the tetrameric Escherichia coli Lac repressor protein. The distortions of the double helix induced by a second protein--the nonspecific nucleoid protein HU--increase the computed likelihood of looping by several orders of magnitude over that of DNA alone. Large-scale deformations of the repressor, sequence-dependent features in the DNA loop, and deformability of the DNA operators also enhance looping, although to lesser degrees. The correspondence between the predicted looping propensities and the ease of looping derived from gene-expression and single-molecule measurements lends credence to the derived structural picture.

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.

Efficient search, mapping, and optimization of multi-protein genetic systems in diverse bacteria

Developing predictive models of multi-protein genetic systems to understand and optimize their behavior remains a combinatorial challenge, particularly when measurement throughput is limited. We developed a computational approach to build predictive models and identify optimal sequences and expression levels, while circumventing combinatorial explosion. Maximally informative genetic system variants were first designed by the RBS Library Calculator, an algorithm to design sequences for efficiently searching a multi-protein expression space across a > 10,000-fold range with tailored search parameters and well-predicted translation rates. We validated the algorithm's predictions by characterizing 646 genetic system variants, encoded in plasmids and genomes, expressed in six gram-positive and gram-negative bacterial hosts. We then combined the search algorithm with system-level kinetic modeling, requiring the construction and characterization of 73 variants to build a sequence-expression-activity map (SEAMAP) for a biosynthesis pathway. Using model predictions, we designed and characterized 47 additional pathway variants to navigate its activity space, find optimal expression regions with desired activity response curves, and relieve rate-limiting steps in metabolism. Creating sequence-expression-activity maps accelerates the optimization of many protein systems and allows previous measurements to quantitatively inform future designs.

Principles of genetic circuit design

Cells navigate environments, communicate and build complex patterns by initiating gene expression in response to specific signals. Engineers seek to harness this capability to program cells to perform tasks or create chemicals and materials that match the complexity seen in nature. This Review describes new tools that aid the construction of genetic circuits. Circuit dynamics can be influenced by the choice of regulators and changed with expression 'tuning knobs'. We collate the failure modes encountered when assembling circuits, quantify their impact on performance and review mitigation efforts. Finally, we discuss the constraints that arise from circuits having to operate within a living cell. Collectively, better tools, well-characterized parts and a comprehensive understanding of how to compose circuits are leading to a breakthrough in the ability to program living cells for advanced applications, from living therapeutics to the atomic manufacturing of functional materials.

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

Comparison of the theoretical and real-world evolutionary potential of a genetic circuit

With the development of next-generation sequencing technologies, many large scale experimental efforts aim to map genotypic variability among individuals. This natural variability in populations fuels many fundamental biological processes, ranging from evolutionary adaptation and speciation to the spread of genetic diseases and drug resistance. An interesting and important component of this variability is present within the regulatory regions of genes. As these regions evolve, accumulated mutations lead to modulation of gene expression, which may have consequences for the phenotype. A simple model system where the link between genetic variability, gene regulation and function can be studied in detail is missing. In this article we develop a model to explore how the sequence of the wild-type lac promoter dictates the fold-change in gene expression. The model combines single-base pair resolution maps of transcription factor and RNA polymerase binding energies with a comprehensive thermodynamic model of gene regulation. The model was validated by predicting and then measuring the variability of lac operon regulation in a collection of natural isolates. We then implement the model to analyze the sensitivity of the promoter sequence to the regulatory output, and predict the potential for regulation to evolve due to point mutations in the promoter region.

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.

Weak operator binding enhances simulated Lac repressor-mediated DNA looping

The 50th anniversary of Biopolymers coincides closely with the like celebration of the discovery of the Escherichia coli (lac) lactose operon, a classic genetic system long used to illustrate the influence of biomolecular structure on function. The looping of DNA induced by the binding of the Lac repressor protein to sequentially distant operator sites on DNA continues to serve as a paradigm for understanding long-range genomic communication. Advances in analyses of DNA structures and in incorporation of proteins in computer simulations of DNA looping allow us to address long-standing questions about the role of protein-mediated DNA loop formation in transcriptional control. Here we report insights gained from studies of the sequence-dependent contributions of the natural lac operators to Lac repressor-mediated DNA looping. Novel superposition of the ensembles of protein-bound operator structures derived from NMR measurements reveals variations in DNA folding missed in conventional structural alignments. The changes in folding affect the predicted ease with which the repressor induces loop formation and the ways that DNA closes between the protein headpieces. The peeling of the auxiliary operators away from the repressor enhances the formation of loops with the 92-bp wildtype spacing and hints of a structural reason behind their weak binding.

DNA sequence-dependent mechanics and protein-assisted bending in repressor-mediated loop formation

As the chief informational molecule of life, DNA is subject to extensive physical manipulations. The energy required to deform double-helical DNA depends on sequence, and this mechanical code of DNA influences gene regulation, such as through nucleosome positioning. Here we examine the sequence-dependent flexibility of DNA in bacterial transcription factor-mediated looping, a context for which the role of sequence remains poorly understood. Using a suite of synthetic constructs repressed by the Lac repressor and two well-known sequences that show large flexibility differences in vitro, we make precise statistical mechanical predictions as to how DNA sequence influences loop formation and test these predictions using in vivo transcription and in vitro single-molecule assays. Surprisingly, sequence-dependent flexibility does not affect in vivo gene regulation. By theoretically and experimentally quantifying the relative contributions of sequence and the DNA-bending protein HU to DNA mechanical properties, we reveal that bending by HU dominates DNA mechanics and masks intrinsic sequence-dependent flexibility. Such a quantitative understanding of how mechanical regulatory information is encoded in the genome will be a key step towards a predictive understanding of gene regulation at single-base pair resolution.

Poly(dA:dT)-rich DNAs are highly flexible in the context of DNA looping

Large-scale DNA deformation is ubiquitous in transcriptional regulation in prokaryotes and eukaryotes alike. Though much is known about how transcription factors and constellations of binding sites dictate where and how gene regulation will occur, less is known about the role played by the intervening DNA. In this work we explore the effect of sequence flexibility on transcription factor-mediated DNA looping, by drawing on sequences identified in nucleosome formation and ligase-mediated cyclization assays as being especially favorable for or resistant to large deformations. We examine a poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong to a class of highly inflexible DNAs; two strong nucleosome positioning sequences that share a set of particular sequence features common to nucleosome-preferring DNAs; and a CG-rich sequence representative of high G+C-content genomic regions that correlate with high nucleosome occupancy in vivo. To measure the flexibility of these sequences in the context of DNA looping, we combine the in vitro single-molecule tethered particle motion assay, a canonical looping protein, and a statistical mechanical model that allows us to quantitatively relate the looping probability to the looping free energy. We show that, in contrast to the case of nucleosome occupancy, G+C content does not positively correlate with looping probability, and that despite sharing sequence features that are thought to determine nucleosome affinity, the two strong nucleosome positioning sequences behave markedly dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich DNA that is often characterized as highly inflexible in fact exhibits one of the highest propensities for looping that we have measured. These results argue for a need to revisit our understanding of the mechanical properties of DNA in a way that will provide a basis for understanding DNA deformation over the entire range of biologically relevant scenarios that are impacted by DNA deformability.

Strong bending of the DNA double helix

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10⁻⁵, it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.

Multiple interferon regulatory factor and NF-κB sites cooperate in mediating cell-type- and maturation-specific activation of the human CD83 promoter in dendritic cells

CD83 is one of the best-known surface markers for fully mature dendritic cells (mature DCs), and its cell-type- and maturation-specific regulation makes the CD83 promoter an interesting tool for the genetic modulation of DCs. To determine the mechanisms regulating this DC- and maturation-specific CD83 expression, chromatin immunoprecipitation (ChIP)-on-chip microarray, biocomputational, reporter, electrophoretic mobility shift assay (EMSA), and ChIP analyses were performed. These studies led to the identification of a ternary transcriptional activation complex composed of an upstream regulatory element, a minimal promoter, and an enhancer, which have not been reported in this arrangement for any other gene so far. Notably, these DNA regions contain a complex framework of interferon regulatory factor (IRF)- and NF-κB transcription factor-binding sites mediating their arrangement. Mutation of any of the IRF-binding sites resulted in a significant loss of promoter activity, whereas overexpression of NF-κB transcription factors clearly enhanced transcription. We identified IRF-1, IRF-2, IRF-5, p50, p65, and cRel to be involved in regulating maturation-specific CD83 expression in DCs. Therefore, the characterization of this promoter complex not only contributes to the knowledge of DC-specific gene regulation but also suggests the involvement of a transcriptional module with binding sites separated into distinct regions in transcriptional activation as well as cell-type- and maturation-specific transcriptional targeting of DCs.

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.

Theoretical and experimental dissection of DNA loop-mediated repression

Transcriptional networks across all domains of life feature a wide range of regulatory architectures. Theoretical models now make clear predictions about how key parameters describing those architectures modulate gene expression, and the ability to construct genetic circuits with tunable parameters enables precise tests of such models. We dissect gene regulation through DNA looping by tuning network parameters such as repressor copy number, DNA binding strengths, and loop length in both thermodynamic models and experiments. Our results help clarify the short-length mechanical properties of DNA.