FRET studies of a landscape of Lac repressor-mediated DNA loops

Designed architectural proteins that tune DNA looping in bacteria

Architectural proteins alter the shape of DNA. Some distort the double helix by introducing sharp kinks. This can serve to relieve strain in tightly-bent DNA structures. Here, we design and test artificial architectural proteins based on a sequence-specific Transcription Activator-like Effector (TALE) protein, either alone or fused to a eukaryotic high mobility group B (HMGB) DNA-bending domain. We hypothesized that TALE protein binding would stiffen DNA to bending and twisting, acting as an architectural protein that antagonizes the formation of small DNA loops. In contrast, fusion to an HMGB domain was hypothesized to generate a targeted DNA-bending architectural protein that facilitates DNA looping. We provide evidence from Escherichia coli Lac repressor gene regulatory loops supporting these hypotheses in living bacteria. Both data fitting to a thermodynamic DNA looping model and sophisticated molecular modeling support the interpretation of these results. We find that TALE protein binding inhibits looping by stiffening DNA to bending and twisting, while the Nhp6A domain enhances looping by bending DNA without introducing twisting flexibility. Our work illustrates artificial approaches to sculpt DNA geometry with functional consequences. Similar approaches may be applicable to tune the stability of small DNA loops in eukaryotes.

SerraNA: a program to determine nucleic acids elasticity from simulation data

AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly.

Base-Pair Mismatch Can Destabilize Small DNA Loops through Cooperative Kinking

Base-pair mismatch can relieve mechanical stress in highly strained DNA molecules, but how it affects their kinetic stability is not known. Using single-molecule fluorescence resonance energy transfer, we measured the lifetimes of tightly bent DNA loops with and without base-pair mismatch. Surprisingly, for loops captured by stackable sticky ends which leave single-stranded DNA breaks (or nicks) upon annealing, the mismatch decreased the loop lifetime despite reducing the overall bending stress, and the decrease was largest when the mismatch was placed at the DNA midpoint. These findings suggest that base-pair mismatch increases bending stress at the opposite side of the loop through an allosteric mechanism known as cooperative kinking. Based on this mechanism, we present a three-state model that explains the apparent dichotomy between thermodynamic and kinetic stability.

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.

Optimized expression and purification of biophysical quantities of Lac repressor and Lac repressor regulatory domain

The recombinant production of Lac repressor (LacI) in Escherichia coli is complicated by its ubiquitous use as a regulatory element in commercially-available expression vectors and host strains. While LacI-regulated expression systems are often used to produce recombinant LacI, the product can be heterogeneous and unsuitable for some studies. Alternative approaches include using unregulated vectors which typically suffer from low yield or vectors with promoters induced by metabolically active sugars which can dilute isotope labels necessary for certain biophysical studies. Here, an optimized expression system and isolation protocol for producing various constructs of LacI is introduced which eliminates these complications. The expression vector is an adaptation of the pASK backbone wherein expression of the lacI gene is regulated by an anhydrotetracyline inducible tetA promoter and the host strain lacks the lacI gene. Typical yields in highly deuterated minimal medium are nearly 2-fold greater than those previously reported. Notably, the new expression system is also able to produce the isolated regulatory domain of LacI without co-expression of the full-length protein and without any defects in cell viability, eliminating the inconvenient requirement for strict monitoring of cell densities during pre-culturing. Typical yields in highly deuterated minimal medium are significantly greater than those previously reported. Characterization by solution NMR shows that LacI constructs produced using this expression system are highly homogenous and functionally active.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

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.

Tethered particle motion reveals that LacI·DNA loops coexist with a competitor-resistant but apparently unlooped conformation

The lac repressor protein (LacI) efficiently represses transcription of the lac operon in Escherichia coli by binding to two distant operator sites on the bacterial DNA and causing the intervening DNA to form a loop. We employed single-molecule tethered particle motion to observe LacI-mediated loop formation and breakdown in DNA constructs that incorporate optimized operator binding sites and intrinsic curvature favorable to loop formation. Previous bulk competition assays indirectly measured the loop lifetimes in these optimized DNA constructs as being on the order of days; however, we measured these same lifetimes to be on the order of minutes for both looped and unlooped states. In a range of single-molecule DNA competition experiments, we found that the resistance of the LacI-DNA complex to competitive binding is a function of both the operator strength and the interoperator sequence. To explain these findings, we present what we believe to be a new kinetic model of loop formation and DNA competition. In this proposed new model, we hypothesize a new unlooped state in which the unbound DNA-binding domain of the LacI protein interacts nonspecifically with nonoperator DNA adjacent to the operator site at which the second LacI DNA-binding domain is bound.

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.

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.

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.

Rationally designed coiled-coil DNA looping peptides control DNA topology

Artificial DNA looping peptides were engineered to study the roles of protein and DNA flexibility in controlling the geometry and stability of protein-mediated DNA loops. These LZD (leucine zipper dual-binding) peptides were derived by fusing a second, C-terminal, DNA-binding region onto the GCN4 bZip peptide. Two variants with different coiled-coil lengths were designed to control the relative orientations of DNA bound at each end. Electrophoretic mobility shift assays verified formation of a sandwich complex containing two DNAs and one peptide. Ring closure experiments demonstrated that looping requires a DNA-binding site separation of 310 bp, much longer than the length needed for natural loops. Systematic variation of binding site separation over a series of 10 constructs that cyclize to form 862-bp minicircles yielded positive and negative topoisomers because of two possible writhed geometries. Periodic variation in topoisomer abundance could be modeled using canonical DNA persistence length and torsional modulus values. The results confirm that the LZD peptides are stiffer than natural DNA looping proteins, and they suggest that formation of short DNA loops requires protein flexibility, not unusual DNA bendability. Small, stable, tunable looping peptides may be useful as synthetic transcriptional regulators or components of protein–DNA nanostructures.

Structural insights into the role of architectural proteins in DNA looping deduced from computer simulations

Bacterial gene expression is regulated by DNA elements that often lie far apart along the genomic sequence, but come close together during genetic processing. The intervening residues form loops, which are organized by the binding of various proteins. For example, the Escherichia coli Lac repressor protein binds DNA operators, separated by 92 or 401 bp, and suppresses the formation of gene products involved in the metabolism of lactose. The system also includes several highly abundant architectural proteins, such as the histone-like (heat-unstable) HU protein, which severely deform the double helix upon binding. In order to gain a better understanding of how the naturally stiff DNA double helix forms the short loops detected in vivo, we have developed new computational methods to study the effects of various non-specific binding proteins on the three-dimensional configurational properties of DNA sequences. The present article surveys the approach that we use to generate ensembles of spatially constrained protein-decorated DNA structures (minicircles and Lac repressor-mediated loops) and presents some of the insights gained from the correspondence between computation and experiment about the potential contributions of architectural and regulatory proteins to DNA looping and gene expression.

Interplay of protein and DNA structure revealed in simulations of the lac operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.

LacI-DNA-IPTG loops: equilibria among conformations by single-molecule FRET

The E. coli Lac repressor (LacI) tetramer binds simultaneously to a promoter-proximal DNA binding site (operator) and an auxiliary operator, resulting in a DNA loop, which increases repression efficiency. Induction of the lac operon by allolactose reduces the affinity of LacI for DNA, but induction does not completely prevent looping in vivo. Our previous work on the conformations of LacI loops used a hyperstable model DNA construct, 9C14, that contains a sequence directed bend flanked by operators. Single-molecule fluorescence resonance energy transfer (SM-FRET) on a dual fluorophore-labeled LacI-9C14 loop showed that it adopts a single, stable, high-FRET V-shaped LacI conformation. Ligand-induced changes in loop geometry can affect loop stability, and the current work assesses loop population distributions for LacI-9C14 complexes containing the synthetic inducer IPTG. SM-FRET confirms that the high-FRET LacI-9C14 loop is only partially destabilized by saturating IPTG. LacI titration experiments and FRET fluctuation analysis suggest that the addition of IPTG induces loop conformational dynamics and re-equilibration between loop population distributions that include a mixture of looped states that do not exhibit high-efficiency FRET. The results show that repression by looping even at saturating IPTG should be considered in models for regulation of the operon. We propose that persistent DNA loops near the operator function biologically to accelerate rerepression upon exhaustion of inducer.

Mechanism of promoter repression by Lac repressor-DNA loops

The Escherichia coli lactose (lac) operon encodes the first genetic switch to be discovered, and lac remains a paradigm for studying negative and positive control of gene expression. Negative control is believed to involve competition of RNA polymerase and Lac repressor for overlapping binding sites. Contributions to the local Lac repressor concentration come from free repressor and repressor delivered to the operator from remote auxiliary operators by DNA looping. Long-standing questions persist concerning the actual role of DNA looping in the mechanism of promoter repression. Here, we use experiments in living bacteria to resolve four of these questions. We show that the distance dependence of repression enhancement is comparable for upstream and downstream auxiliary operators, confirming the hypothesis that repressor concentration increase is the principal mechanism of repression loops. We find that as few as four turns of DNA can be constrained in a stable loop by Lac repressor. We show that RNA polymerase is not trapped at repressed promoters. Finally, we show that constraining a promoter in a tight DNA loop is sufficient for repression even when promoter and operator do not overlap.