Intrinsic curvature of DNA influences LacR-mediated looping

Environmentally Controlled Curvature of Single Collagen Proteins

The predominant structural protein in vertebrates is collagen, which plays a key role in extracellular matrix and connective tissue mechanics. Despite its prevalence and physical importance in biology, the mechanical properties of molecular collagen are far from established. The flexibility of its triple helix is unresolved, with descriptions from different experimental techniques ranging from flexible to semirigid. Furthermore, it is unknown how collagen type (homo- versus heterotrimeric) and source (tissue derived versus recombinant) influence flexibility. Using SmarTrace, a chain-tracing algorithm we devised, we performed statistical analysis of collagen conformations collected with atomic force microscopy to determine the protein's mechanical properties. Our results show that types I, II, and III collagens-the key fibrillar varieties-exhibit similar molecular flexibilities. However, collagen conformations are strongly modulated by salt, transitioning from compact to extended as KCl concentration increases in both neutral and acidic pH. Although analysis with a standard worm-like chain model suggests that the persistence length of collagen can attain a wide range of values within the literature range, closer inspection reveals that this modulation of collagen's conformational behavior is not due to changes in flexibility but rather arises from the induction of curvature (either intrinsic or induced by interactions with the mica surface). By modifying standard polymer theory to include innate curvature, we show that collagen behaves as an equilibrated curved worm-like chain in two dimensions. Analysis within the curved worm-like chain model shows that collagen's curvature depends strongly on pH and salt, whereas its persistence length does not. Thus, we find that triple-helical collagen is well described as semiflexible irrespective of source, type, pH, and salt environment. These results demonstrate that collagen is more flexible than its conventional description as a rigid rod, which may have implications for its cellular processing and secretion.

Dynamic Model Exposes the Energetics and Dynamics of the Injection Machinery for Bacteriophage T4

Bacteriophage T4 infects the bacterial host (Escherichia coli) using an efficient genomic delivery machine that is driven by elastic energy stored in a contractile tail sheath. Although the atomic structure of T4 is largely known, the dynamics of its fascinating injection machinery is not understood. This article contributes, to our knowledge, the first predictions of the energetics and dynamics of the T4 injection machinery using a novel dynamic model. The model employs an atomistic (molecular dynamics) representation of a fraction of the sheath structure to generate a continuum model of the entire sheath that also couples to a model of the viral capsid and tail tube. The resulting model of the entire injection machine reveals estimates for the energetics, timescale, and pathway of the T4 injection process as well as the force available for cell rupture. It also reveals the large and highly nonlinear conformational changes of the sheath whose elastic energy drives the injection process.

Simulation of DNA Supercoil Relaxation

Several recent single-molecule experiments observe the response of supercoiled DNA to nicking endonucleases and topoisomerases. Typically in these experiments, indirect measurements of supercoil relaxation are obtained by observing the motion of a large micron-sized bead. The bead, which also serves to manipulate DNA, experiences significant drag and thereby obscures supercoil dynamics. Here we employ our discrete wormlike chain model to bypass experimental limitations and simulate the dynamic response of supercoiled DNA to a single strand nick. From our simulations, we make three major observations. First, extension is a poor dynamic measure of supercoil relaxation; in fact, the linking number relaxes so fast that it cannot have much impact on extension. Second, the rate of linking number relaxation depends upon its initial partitioning into twist and writhe as determined by tension. Third, the extensional response strongly depends upon the initial position of plectonemes.

DNA buckling in bacteriophage cavities as a mechanism to aid virus assembly

While relatively simple biologically, bacteriophages are sophisticated biochemical machines that execute a precise sequence of events during virus assembly, DNA packaging, and ejection. These stages of the viral life cycle require intricate coordination of viral components whose structures are being revealed by single molecule experiments and high resolution (cryo-electron microscopy) reconstructions. For example, during packaging, bacteriophages employ some of the strongest known molecular motors to package DNA against increasing pressure within the viral capsid shell. Located upstream of the motor is an elaborate portal system through which DNA is threaded. A high resolution reconstruction of the portal system for bacteriophage ϕ29 reveals that DNA buckles inside a small cavity under large compressive forces. In this study, we demonstrate that DNA can also buckle in other bacteriophages including T7 and P22. Using a computational rod model for DNA, we demonstrate that a DNA buckle can initiate and grow within the small confines of a cavity under biologically-attainable force levels. The forces of DNA-cavity contact and DNA-DNA electrostatic repulsion ultimately limit cavity filling. Despite conforming to very different cavity geometries, the buckled DNA within T7 and P22 exhibits near equal volumetric energy density (∼1kT/nm(3)) and energetic cost of packaging (∼22kT). We hypothesize that a DNA buckle creates large forces on the cavity interior to signal the conformational changes to end packaging. In addition, a DNA buckle may help retain the genome prior to tail assembly through significantly increased contact area with the portal.

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.

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.

Structural ensemble and dynamics of toroidal-like DNA shapes in bacteriophage ϕ29 exit cavity

In the bacteriophage ϕ29, DNA is packed into a preassembled capsid from which it ejects under high pressure. A recent cryo-EM reconstruction of ϕ29 revealed a compact toroidal DNA structure (30-40 basepairs) lodged within the exit cavity formed by the connector-lower collar protein complex. Using multiscale models, we compute a detailed structural ensemble of intriguing DNA toroids of various lengths, all highly compatible with experimental observations. In particular, coarse-grained (elastic rod) and atomistic (molecular dynamics) models predict the formation of DNA toroids under significant compression, a largely unexplored state of DNA. Model predictions confirm that a biologically attainable compressive force of 25 pN sustains the toroid and yields DNA electron density maps highly consistent with the experimental reconstruction. The subsequent simulation of dynamic toroid ejection reveals large reactions on the connector that may signal genome release.

Effect of 8-oxoguanine on DNA structure and deformability

8-Oxoguanine (oxoG) is an abundant product of oxidative DNA damage. It is removed by repair glycosylases, but exactly how the enzymes recognize oxoG in the large surplus of undamaged bases is not fully understood. The lesion may induce changes in the properties of naked DNA that facilitate the recognition. In this work, we assess the effect of oxoG on DNA structure and mechanical deformability. We performed extensive unrestrained, atomic resolution molecular dynamics simulations to parametrize a nonlocal, rigid base mechanical model of DNA. Our data indicate that oxoG induces unwinding of the base pair step at the 5'-side of the lesion. This brings the damaged DNA closer to its conformation in the initial complex with bacterial glycosylase MutM. The untwisting is partially caused by different BII substate populations and is further enhanced by the base-sugar repulsion within oxoG. On the other hand, our analysis shows that damaged and undamaged DNA have very similar harmonic stiffness. These results suggest an indirect readout component of the MutM-DNA initial complex formation. They also help one to understand the effect of oxoG on the formation of nucleosomes and looped gene regulatory complexes.

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.

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

DNA looping mediated by the Lac repressor is an archetypal test case for modeling protein and DNA flexibility. Understanding looping is fundamental to quantitative descriptions of gene expression. Systematic analysis of LacI•DNA looping was carried out using a landscape of DNA constructs with lac operators bracketing an A-tract bend, produced by varying helical phasings between operators and the bend. Fluorophores positioned on either side of both operators allowed direct Förster resonance energy transfer (FRET) detection of parallel (P1) and antiparallel (A1, A2) DNA looping topologies anchored by V-shaped LacI. Combining fluorophore position variant landscapes allows calculation of the P1, A1 and A2 populations from FRET efficiencies and also reveals extended low-FRET loops proposed to form via LacI opening. The addition of isopropyl-β-d-thio-galactoside (IPTG) destabilizes but does not eliminate the loops, and IPTG does not redistribute loops among high-FRET topologies. In some cases, subsequent addition of excess LacI does not reduce FRET further, suggesting that IPTG stabilizes extended or other low-FRET loops. The data align well with rod mechanics models for the energetics of DNA looping topologies. At the peaks of the predicted energy landscape for V-shaped loops, the proposed extended loops are more stable and are observed instead, showing that future models must consider protein flexibility.

DNA modeling reveals an extended lac repressor conformation in classic in vitro binding assays

Protein-mediated DNA looping, such as that induced by the lactose repressor (LacI) of Escherichia coli, is a well-known gene regulation mechanism. Although researchers have given considerable attention to DNA looping by LacI, many unanswered questions about this mechanism, including the role of protein flexibility, remain. Recent single-molecule observations suggest that the two DNA-binding domains of LacI are capable of splaying open about the tetramerization domain into an extended conformation. We hypothesized that if recent experiments were able to reveal the extended conformation, it is possible that such structures occurred in previous studies as well. In this study, we tested our hypothesis by reevaluating two classic in vitro binding assays using a computational rod model of DNA. The experiments and computations evaluate the looping of both linear DNA and supercoiled DNA minicircles over a broad range of DNA interoperator lengths. The computed energetic minima align well with the experimentally observed interoperator length for optimal loop stability. Of equal importance, the model reveals that the most stable loops for linear DNA occur when LacI adopts the extended conformation. In contrast, for DNA minicircles, optimal stability may arise from either the closed or the extended protein conformation depending on the degree of supercoiling and the interoperator length.

A multiscale dynamic model of DNA supercoil relaxation by topoisomerase IB

In this study, we report what we believe to be the first multiscale simulation of the dynamic relaxation of DNA supercoils by human topoisomerase IB (topo IB). We leverage our previous molecular dynamics calculations of the free energy landscape describing the interaction between a short DNA fragment and topo IB. Herein, this landscape is used to prescribe boundary conditions for a computational, elastodynamic continuum rod model of a long length of supercoiled DNA. The rod model, which accounts for the nonlinear bending, twisting, and electrostatic interaction of the (negatively charged) DNA backbone, is extended to include the hydrodynamic drag induced by the surrounding physiological buffer. Simulations for a 200-bp-long DNA supercoil in complex with topo IB reveal a relaxation timescale of ∼0.1-1.0 μs. The relaxation follows a sequence of cascading reductions in the supercoil linking number (Lk), twist (Tw), and writhe (Wr) that follow companion cascading reductions in the supercoil elastic and electrostatic energies. The novel (to our knowledge) multiscale modeling method may enable simulations of the entire experimental setup that measures DNA supercoiling and relaxation via single molecule magnetic trapping.

Roles of DNA looping in enhancer-blocking activity

Enhancer-promoter interactions in eukaryotic genomes are often controlled by sequence elements that block the actions of enhancers. Although the experimental evidence suggests that those sequence elements contribute to forming loops of chromatin, the molecular mechanism of how such looping affects the enhancer-blocking activity is still largely unknown. In this article, the roles of DNA looping in enhancer blocking are investigated by numerically simulating the DNA conformation of a prototypical model system of gene regulation. The simulated results show that the enhancer function is indeed blocked when the enhancer is looped out so that it is separated from the promoter, which explains experimental observations of gene expression in the model system. The local structural distortion of DNA caused by looping is important for blocking, so the ability of looping to block enhancers can be lost when the loop length is much larger than the persistence length of the chain.

Gene repression by minimal lac loops in vivo

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

Wrapped-around models for the lac operon complex

The protein-DNA complex, involved in the lac operon of enteric bacteria, is paradigmatic in understanding the extent of DNA bending and plasticity due to interactions with protein assemblies acting as DNA regulators. For the lac operon, two classes of structures have been proposed: 1), with the protein tetramer lying away from the DNA loop (wrapped-away model); and 2), with the protein tetramer lying inside the DNA loop (wrapped-around model). A recently developed electrostatic analytical model shows that the size and net charge of the Lac protein tetramer allow the bending of DNA, which is consistent with another wrapped-around model from the literature. Coarse-grained models, designed based on this observation, are extensively investigated and show three kinds of wrapped-around arrangements of DNA and a lower propensity for wrapped-away configurations. Molecular dynamics simulations of an all-atom model, built on the basis of the most tightly collapsed coarse-grained model, show that most of the DNA double-helical architecture is maintained in the region between O3 and O1 DNA operators, that the DNA distortion is concentrated in the chain beyond the O1 operator, and that the protein tetramer can adapt the N-terminal domains to the DNA tension.

Design and characterization of molecular tools for a Synthetic Biology approach towards developing cyanobacterial biotechnology

Cyanobacteria are suitable for sustainable, solar-powered biotechnological applications. Synthetic biology connects biology with computational design and an engineering perspective, but requires efficient tools and information about the function of biological parts and systems. To enable the development of cyanobacterial Synthetic Biology, several molecular tools were developed and characterized: (i) a broad-host-range BioBrick shuttle vector, pPMQAK1, was constructed and confirmed to replicate in Escherichia coli and three different cyanobacterial strains. (ii) The fluorescent proteins Cerulean, GFPmut3B and EYFP have been demonstrated to work as reporter proteins in cyanobacteria, in spite of the strong background of photosynthetic pigments. (iii) Several promoters, like P_rnpB and variants of P_rbcL, and a version of the promoter P_trc with two operators for enhanced repression, were developed and characterized in Synechocystis sp. strain PCC6803. (iv) It was shown that a system for targeted protein degradation, which is needed to enable dynamic expression studies, is working in Synechocystis sp. strain PCC6803. The pPMQAK1 shuttle vector allows the use of the growing numbers of BioBrick parts in many prokaryotes, and the other tools herein implemented facilitate the development of new parts and systems in cyanobacteria.

First-principles calculation of DNA looping in tethered particle experiments

We calculate the probability of DNA loop formation mediated by regulatory proteins such as Lac repressor (LacI), using a mathematical model of DNA elasticity. Our model is adapted to calculating quantities directly observable in tethered particle motion (TPM) experiments, and it accounts for all the entropic forces present in such experiments. Our model has no free parameters; it characterizes DNA elasticity using information obtained in other kinds of experiments. It assumes a harmonic elastic energy function (or wormlike chain type elasticity), but our Monte Carlo calculation scheme is flexible enough to accommodate arbitrary elastic energy functions. We show how to compute both the 'looping J factor' (or equivalently, the looping free energy) for various DNA construct geometries and LacI concentrations, as well as the detailed probability density function of bead excursions. We also show how to extract the same quantities from recent experimental data on TPM, and then compare to our model's predictions. In particular, we present a new method to correct observed data for finite camera shutter time and other experimental effects. Although the currently available experimental data give large uncertainties, our first-principles predictions for the looping free energy change are confirmed to within about 1 k(B)T, for loops of length around 300 basepairs. More significantly, our model successfully reproduces the detailed distributions of bead excursion, including their surprising three-peak structure, without any fit parameters and without invoking any alternative conformation of the LacI tetramer. Indeed, the model qualitatively reproduces the observed dependence of these distributions on tether length (e.g., phasing) and on LacI concentration (titration). However, for short DNA loops (around 95 basepairs) the experiments show more looping than is predicted by the harmonic-elasticity model, echoing other recent experimental results. Because the experiments we study are done in vitro, this anomalously high looping cannot be rationalized as resulting from the presence of DNA-bending proteins or other cellular machinery. We also show that it is unlikely to be the result of a hypothetical 'open' conformation of the LacI tetramer.

Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.

Mesoscale modeling of multi-protein-DNA assemblies: the role of the catabolic activator protein in Lac-repressor-mediated looping

DNA looping plays a key role in the regulation of the lac operon in Escherichia coli. The presence of a tightly bent loop (between sequentially distant sites of Lac repressor protein binding) purportedly hinders the binding of RNA polymerase and subsequent transcription of the genetic message. The unexpectedly favorable binding interaction of this protein-DNA assembly with the catabolic activator protein (CAP), a protein that also bends DNA and paradoxically facilitates the binding of RNA polymerase, stimulated extension of our base-pair level theory of DNA elasticity to the treatment of DNA loops formed in the presence of several proteins. Here we describe in detail a procedure to determine the structures and free energies of multi-protein-DNA assemblies and illustrate the predicted effects of CAP binding on the configurations of the wild-type 92-bp Lac repressor-mediated O3-O1 DNA loop. We show that the DNA loop adopts an antiparallel orientation in the most likely structure and that this loop accounts for the published experimental observation that, when CAP is bound to the loop, one of the arms of LacR binds to an alternative site that is displaced from the original site by 5 bp.