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

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.

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.

Supercoiling and looping promote DNA base accessibility and coordination among distant sites

DNA in cells is supercoiled and constrained into loops and this supercoiling and looping influence every aspect of DNA activity. We show here that negative supercoiling transmits mechanical stress along the DNA backbone to disrupt base pairing at specific distant sites. Cooperativity among distant sites localizes certain sequences to superhelical apices. Base pair disruption allows sharp bending at superhelical apices, which facilitates DNA writhing to relieve torsional strain. The coupling of these processes may help prevent extensive denaturation associated with genomic instability. Our results provide a model for how DNA can form short loops, which are required for many essential processes, and how cells may use DNA loops to position nicks to facilitate repair. Furthermore, our results reveal a complex interplay between site-specific disruptions to base pairing and the 3-D conformation of DNA, which influences how genomes are stored, replicated, transcribed, repaired, and many other aspects of DNA activity.

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.

Biophysical Reviews' "Meet the Editors Series"-a profile of Wilma K. Olson

As one of the five Executive Editors of Biophysical Reviews I have been asked to provide this short biographical sketch for the readers of the journal. I have been a member of the Editorial Board since the inception of the journal in 2008 and an Executive Editor since 2014. I hold a B.S. degree in Chemistry from the University of Delaware and a Ph.D. in Physical Chemistry from Stanford University. Except for a year as a Damon Runyon Postdoctoral Fellow at Columbia University, I have spent my entire professional career at Rutgers, the State University of New Jersey, where I am currently the Mary I. Bunting Professor of Chemistry and Chemical Biology. I served for many years as Founding Director of the Rutgers University Center for Molecular Biophysics and Biophysical Chemistry, and have trained undergraduate, graduate, and postdoctoral students from a variety of academic disciplines.

The mechanics of DNA loops bridged by proteins unveiled by single-molecule experiments

Protein-induced DNA bridging and looping is a common mechanism for various and essential processes in bacterial chromosomes. This mechanism is preserved despite the very different bacterial conditions and their expected influence on the thermodynamic and kinetic characteristics of the bridge formation and stability. Over the last two decades, single-molecule techniques carried out on in vitro DNA systems have yielded valuable results which, in combination with theoretical works, have clarified the effects of different parameters of nucleoprotein complexes on the protein-induced DNA bridging and looping process. In this review, I will outline the features that can be measured for such processes with various single-molecule techniques in use in the field. I will then describe both the experimental results and the theoretical models that illuminate the contribution of the DNA molecule itself as well as that of the bridging proteins in the DNA looping mechanism at play in the nucleoid of E. coli.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

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.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Supercoiling biases the formation of loops involved in gene regulation

The function of DNA as a repository of genetic information is well-known. The post-genomic effort is to understand how this information-containing filament is chaperoned to manage its compaction and topological states. Indeed, the activities of enzymes that transcribe, replicate, or repair DNA are regulated to a large degree by access. Proteins that act at a distance along the filament by binding at one site and contacting another site, perhaps as part of a bigger complex, create loops that constitute topological domains and influence regulation. DNA loops and plectonemes are not necessarily spontaneous, especially large loops under tension for which high energy is required to bring their ends together, or small loops that require accessory proteins to facilitate DNA bending. However, the torsion in stiff filaments such as DNA dramatically modulates the topology, driving it from extended and genetically accessible to more looped and compact, genetically secured forms. Furthermore, there are accessory factors that bias the response of the DNA filament to supercoiling. For example, small molecules like polyamines, which neutralize the negative charge repulsions along the phosphate backbone, enhance flexibility and promote writhe over twist in response to torsion. Such increased flexibility likely pushes the topological equilibrium from twist toward writhe at tensions thought to exist in vivo. A predictable corollary is that stiffening DNA antagonizes looping and bending. Certain sequences are known to be more or less flexible or to exhibit curvature, and this may affect interactions with binding proteins. In vivo all of these factors operate simultaneously on DNA that is generally negatively supercoiled to some degree. Therefore, in order to better understand gene regulation that involves protein-mediated DNA loops, it is critical to understand the thermodynamics and kinetics of looping in DNA that is under tension, negatively supercoiled, and perhaps exposed to molecules that alter elasticity. Recent experiments quantitatively reveal how much negatively supercoiling DNA lowers the free energy of looping, possibly biasing the operation of genetic switches.

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.

The role of structural parameters in DNA cyclization

BACKGROUND

The intrinsic bendability of DNA plays an important role with relevance for myriad of essential cellular mechanisms. The flexibility of a DNA fragment can be experimentally and computationally examined by its propensity for cyclization, quantified by the Jacobson-Stockmayer J factor. In this study, we use a well-established coarse-grained three-dimensional model of DNA and seven distinct sets of experimentally and computationally derived conformational parameters of the double helix to evaluate the role of structural parameters in calculating DNA cyclization.

RESULTS

We calculate the cyclization rates of 86 DNA sequences with previously measured J factors and lengths between 57 and 325 bp as well as of 20,000 randomly generated DNA sequences with lengths between 350 and 4000 bp. Our comparison with experimental data is complemented with analysis of simulated data.

CONCLUSIONS

Our data demonstrate that all sets of parameters yield very similar results for longer DNA fragments, regardless of the nucleotide sequence, which are in agreement with experimental measurements. However, for DNA fragments shorter than 100 bp, all sets of parameters performed poorly yielding results with several orders of magnitude difference from the experimental measurements. Our data show that DNA cyclization rates calculated using conformational parameters based on nucleosome packaging data are most similar to the experimental measurements. Overall, our study provides a comprehensive large-scale assessment of the role of structural parameters in calculating DNA cyclization rates.

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.

Enhanced tethered-particle motion analysis reveals viscous effects

Tethered-particle motion experiments do not require expensive or technically complex hardware, and increasing numbers of researchers are adopting this methodology to investigate the topological effects of agents that act on the tethering polymer or the characteristics of the polymer itself. These investigations depend on accurate measurement and interpretation of changes in the effective length of the tethering polymer (often DNA). However, the bead size, tether length, and buffer affect the confined diffusion of the bead in this experimental system. To evaluate the effects of these factors, improved measurements to calibrate the two-dimensional range of motion (excursion) versus DNA length were carried out. Microspheres of 160 or 240 nm in radius were tethered by DNA molecules ranging from 225 to 3477 basepairs in length in aqueous buffers containing 100 mM potassium glutamate and 8 mM MgCl2 or 10 mM Tris-HCl and 200 mM KCl, with or without 0.5% Tween added to the buffer, and the motion was recorded. Different buffers altered the excursion of beads on identical DNA tethers. Buffer with only 10 mM NaCl and >5 mM magnesium greatly reduced excursion. Glycerol added to increase viscosity slowed confined diffusion of the tethered beads but did not change excursion. The confined-diffusion coefficients for all tethered beads were smaller than those expected for freely diffusing beads and decreased for shorter tethers. Tethered-particle motion is a sensitive framework for diffusion experiments in which small beads on long leashes most closely resemble freely diffusing, untethered beads.

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.

Bacterial promoter repression by DNA looping without protein-protein binding competition

The Escherichia coli lactose operon provides a paradigm for understanding gene control by DNA looping where the lac repressor (LacI) protein competes with RNA polymerase for DNA binding. Not all promoter loops involve direct competition between repressor and RNA polymerase. This raises the possibility that positioning a promoter within a tightly constrained DNA loop is repressive per se, an idea that has previously only been considered in vitro. Here, we engineer living E. coli bacteria to measure repression due to promoter positioning within such a tightly constrained DNA loop in the absence of protein–protein binding competition. We show that promoters held within such DNA loops are repressed ∼100-fold, with up to an additional ∼10-fold repression (∼1000-fold total) dependent on topological positioning of the promoter on the inner or outer face of the DNA loop. Chromatin immunoprecipitation data suggest that repression involves inhibition of both RNA polymerase initiation and elongation. These in vivo results show that gene repression can result from tightly looping promoter DNA even in the absence of direct competition between repressor and RNA polymerase binding.

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.

Screening for protein-DNA interactions by automatable DNA-protein interaction ELISA

DNA-binding proteins (DBPs), such as transcription factors, constitute about 10% of the protein-coding genes in eukaryotic genomes and play pivotal roles in the regulation of chromatin structure and gene expression by binding to short stretches of DNA. Despite their number and importance, only for a minor portion of DBPs the binding sequence had been disclosed. Methods that allow the de novo identification of DNA-binding motifs of known DBPs, such as protein binding microarray technology or SELEX, are not yet suited for high-throughput and automation. To close this gap, we report an automatable DNA-protein-interaction (DPI)-ELISA screen of an optimized double-stranded DNA (dsDNA) probe library that allows the high-throughput identification of hexanucleotide DNA-binding motifs. In contrast to other methods, this DPI-ELISA screen can be performed manually or with standard laboratory automation. Furthermore, output evaluation does not require extensive computational analysis to derive a binding consensus. We could show that the DPI-ELISA screen disclosed the full spectrum of binding preferences for a given DBP. As an example, AtWRKY11 was used to demonstrate that the automated DPI-ELISA screen revealed the entire range of in vitro binding preferences. In addition, protein extracts of AtbZIP63 and the DNA-binding domain of AtWRKY33 were analyzed, which led to a refinement of their known DNA-binding consensi. Finally, we performed a DPI-ELISA screen to disclose the DNA-binding consensus of a yet uncharacterized putative DBP, AtTIFY1. A palindromic TGATCA-consensus was uncovered and we could show that the GATC-core is compulsory for AtTIFY1 binding. This specific interaction between AtTIFY1 and its DNA-binding motif was confirmed by in vivo plant one-hybrid assays in protoplasts. Thus, the value and applicability of the DPI-ELISA screen for de novo binding site identification of DBPs, also under automatized conditions, is a promising approach for a deeper understanding of gene regulation in any organism of choice.

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.