Effect of supercoiling on formation of protein-mediated DNA loops

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.

Allosteric interactions in a birod model of DNA

Allosteric interactions between molecules bound to DNA at distant locations have been known for a long time. The phenomenon has been studied via experiments and numerical simulations, but a comprehensive understanding grounded in a theory of DNA elasticity remains a challenge. Here, we quantify allosteric interactions between two entities bound to DNA by using the theory of birods. We recognize that molecules bound to DNA cause local deformations that can be captured in a birod model which consists of two elastic strands interacting via an elastic web representing the basepairs. We show that the displacement field caused by bound entities decays exponentially with distance from the binding site. We compute the interaction energy between two proteins on DNA as a function of distance between them and find that it decays exponentially while oscillating with the periodicity of the double helix, in excellent agreement with experiments. The decay length of the interaction energy can be determined in terms of the mechanical properties of the strands and the webbing in our birod model, and it varies with the GC content of the DNA. Our model provides a framework for viewing allosteric interactions in DNA within the ambit of configurational forces of continuum elasticity.

Analytical expressions for the closure probability of a stiff wormlike chain for finite capture radius

Estimating the probability that two monomers of the same polymer chain are close together is a key ingredient to characterize intramolecular reactions and polymer looping. In the case of stiff wormlike polymers (rigid fluctuating elastic rods), for which end-to-end encounters are rare events, we derive an explicit analytical formula for the probability η(r_{c}) that the distance between the chain extremities is smaller than some capture radius r_{c}. The formula is asymptotically exact in the limit of stiff chains, and it leads to the identification of two distinct scaling regimes for the closure factor, originating from a strong variation of the fluctuations of the chain orientation at closure. Our theory is compatible with existing analytical results from the literature that cover the cases of a vanishing capture radius and of nearly fully extended chains.

Structural transitions in torsionally constrained DNA and their dependence on solution electrostatics

Experimental studies on single molecules of DNA have reported a rich variety of structural transitions, including coexistence of three phases, in a torsionally constrained molecule. A comprehensive knowledge of these structural transitions is useful for unraveling the in vivo and in vitro behavior of DNA. Our objective is to understand the structural transitions in a torsionally constrained DNA molecule when it is pulled using optical or magnetic tweezers. We use foundational concepts from the Zimm-Bragg helix-coil transition theory and merge them with ideas from the theory of fluctuating elastic rods to model the mechanics of DNA. We also account for the electrostatic interactions between the ions and the negatively charged phosphate backbone of DNA. Using our model, we calculate the force and torque corresponding to the overstretching transition characterized by a 70% jump in the contour length of the molecule and examine the effect of salt concentration on this transition. We also deduce conditions under which the co-existence of B-, S- and P-DNA is possible. We examine how the cooperativity parameter for each transition affects the force-extension curve or torque-rotation curve. We attempt to rationalize the non-monotonic dependence of external work done on the ion concentration by connecting it to the electrostatic dependence of the interfacial energy between two phases of DNA. Our theoretical results are in agreement with multiple experiments documented in the literature and they generate falsifiable predictions that can be tested in new experiments.

STATEMENT OF SIGNIFICANCE

The overarching objective of this paper is to explore the implications of variation in ion concentration on the structural transitions driven by external forces in a torsionally constrained DNA molecule. A comprehensive understanding of the phase behavior of torsionally constrained DNA is useful because DNA in cells is tightly packaged and is acted upon by molecular machines in different ionic environments. We examine the mechanics of the overstretching transition, characterized by a 70% jump in contour length, wherein a mixture of B- and S-DNA converts into a mixture of S- and P-DNA through a triple point in the phase diagram. Our results are corroborated by experimental data at every step and we make predictions that are experimentally verifiable.

Functional relationship between high mobility group A1 (HMGA1) protein and insulin-like growth factor-binding protein 3 (IGFBP-3) in human chondrocytes

Introduction

Insulin-like growth factor I (IGF-I) regulates articular cartilage homeostasis. During osteoarthritis (OA), the anabolic responses of chondrocytes to IGF-I are likely to be prevented by the enhanced production of IGF-binding proteins (IGFBPs), especially IGFBP-3. The aim of this study is to evaluate whether the architectural transcription factor high mobility group A1 (HMGA1) influences IGFBP-3 overexpression in vitro, in cultured chondrocytic cell lines, and ex vivo, in human osteoarthritic cartilage compared to healthy human cartilage controls.

Methods

Quantitative real-time reverse transcription-PCR (qRT-PCR) was performed to assess the relative transcript levels of HMGA1 and IGFBP-3 in vitro, in the human chondrocytic cell lines T/C-28a4 and C-28/I2. An electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) and transient transfection assays were performed to investigate the HMGA1-IGFBP-3 gene interaction. Samples of articular cartilage were harvested from osteoarthritic patients and controls and analyzed by qRT-PCR for HMGA1 and IGFBP-3 mRNA levels.

Results

A parallelism between HMGA1 protein levels and IGFBP-3 gene expression has been observed in T/C-28a4 and C-28/I2 cells. The interaction of HMGA1 with the IGFBP-3 gene promoter has been demonstrated by EMSA and ChIP. In transient transfections, IGFBP-3 promoter activity increased in cells overexpressing HMGA1 and decreased in cells pretreated with siRNA detected against HMGA1. IGFBP-3 mRNA expression was higher in cartilage from patients with OA, where the increased expression of IGFBP-3 closely paralleled the increased expression of HMGA1 mRNA.

Conclusions

Our observations indicate that increased HMGA1 expression in human chondrocytes is associated with increased expression of IGFBP-3. It is tempting to speculate that, through the regulation of IGFBP3 expression, HMGA1 may act as a pathogenetic factor for OA.

DNA looping by FokI: the impact of twisting and bending rigidity on protein-induced looping dynamics

Protein-induced DNA looping is crucial for many genetic processes such as transcription, gene regulation and DNA replication. Here, we use tethered-particle motion to examine the impact of DNA bending and twisting rigidity on loop capture and release, using the restriction endonuclease FokI as a test system. To cleave DNA efficiently, FokI bridges two copies of an asymmetric sequence, invariably aligning the sites in parallel. On account of the fixed alignment, the topology of the DNA loop is set by the orientation of the sites along the DNA. We show that both the separation of the FokI sites and their orientation, altering, respectively, the twisting and the bending of the DNA needed to juxtapose the sites, have profound effects on the dynamics of the looping interaction. Surprisingly, the presence of a nick within the loop does not affect the observed rigidity of the DNA. In contrast, the introduction of a 4-nt gap fully relaxes all of the torque present in the system but does not necessarily enhance loop stability. FokI therefore employs torque to stabilise its DNA-looping interaction by acting as a ‘torsional’ catch bond.

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.

Looping probabilities of elastic chains: a path integral approach

We consider an elastic chain at thermodynamic equilibrium with a heat bath, and derive an approximation to the probability density function, or pdf, governing the relative location and orientation of the two ends of the chain. Our motivation is to exploit continuum mechanics models for the computation of DNA looping probabilities, but here we focus on explaining the novel analytical aspects in the derivation of our approximation formula. Accordingly, and for simplicity, the current presentation is limited to the illustrative case of planar configurations. A path integral formalism is adopted, and, in the standard way, the first approximation to the looping pdf is obtained from a minimal energy configuration satisfying prescribed end conditions. Then we compute an additional factor in the pdf which encompasses the contributions of quadratic fluctuations about the minimum energy configuration along with a simultaneous evaluation of the partition function. The original aspects of our analysis are twofold. First, the quadratic Lagrangian describing the fluctuations has cross-terms that are linear in first derivatives. This, seemingly small, deviation from the structure of standard path integral examples complicates the necessary analysis significantly. Nevertheless, after a nonlinear change of variable of Riccati type, we show that the correction factor to the pdf can still be evaluated in terms of the solution to an initial value problem for the linear system of Jacobi ordinary differential equations associated with the second variation. The second novel aspect of our analysis is that we show that the Hamiltonian form of these linear Jacobi equations still provides the appropriate correction term in the inextensible, unshearable limit that is commonly adopted in polymer physics models of, e.g. DNA. Prior analyses of the inextensible case have had to introduce nonlinear and nonlocal integral constraints to express conditions on the relative displacement of the end points. Our approximation formula for the looping pdf is of quite general applicability as, in contrast to most prior approaches, no assumption is made of either uniformity of the elastic chain, nor of a straight intrinsic shape. If the chain is uniform the Jacobi system evaluated at certain minimum energy configurations has constant coefficients. In such cases our approximate pdf can be evaluated in an entirely explicit, closed form. We illustrate our analysis with a planar example of this type and compute an approximate probability of cyclization, i.e., of forming a closed loop, from a uniform elastic chain whose intrinsic shape is an open circular arc.

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.

Direct demonstration and quantification of long-range DNA looping by the lambda bacteriophage repressor

Recently, it was proposed that DNA looping by the λ repressor (CI protein) strengthens repression of lytic genes during lysogeny and simultaneously ensures efficient switching to lysis. To investigate this hypothesis, tethered particle motion experiments were performed and dynamic CI-mediated looping of single DNA molecules containing the λ repressor binding sites separated by 2317 bp (the wild-type distance) was quantitatively analyzed. DNA containing all three intact operators or with mutated o3 operators were compared. Modeling the thermodynamic data established the free energy of CI octamer-mediated loop formation as 1.7 kcal/mol, which decreased to –0.7 kcal/mol when supplemented by a tetramer (octamer+tetramer-mediated loop). These results support the idea that loops secured by an octamer of CI bound at oL1, oL2, oR1 and oR2 operators must be augmented by a tetramer of CI bound at the oL3 and oR3 to be spontaneous and stable. Thus the o3 sites are critical for loops secured by the CI protein that attenuate cI expression.

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.

Single-molecule manipulation reveals supercoiling-dependent modulation of lac repressor-mediated DNA looping

Gene expression regulation is a fundamental biological process which deploys specific sets of genomic information depending on physiological or environmental conditions. Several transcription factors (including lac repressor, LacI) are present in the cell at very low copy number and increase their local concentration by binding to multiple sites on DNA and looping the intervening sequence. In this work, we employ single-molecule manipulation to experimentally address the role of DNA supercoiling in the dynamics and stability of LacI-mediated DNA looping. We performed measurements over a range of degrees of supercoiling between -0.026 and +0.026, in the absence of axial stretching forces. A supercoiling-dependent modulation of the lifetimes of both the looped and unlooped states was observed. Our experiments also provide evidence for multiple structural conformations of the LacI-DNA complex, depending on torsional constraints. The supercoiling-dependent modulation demonstrated here adds an important element to the model of the lac operon. In fact, the complex network of proteins acting on the DNA in a living cell constantly modifies its topological and mechanical properties: our observations demonstrate the possibility of establishing a signaling pathway from factors affecting DNA supercoiling to transcription factors responsible for the regulation of specific sets of genes.

Geometry of mediating protein affects the probability of loop formation in DNA

Recent single molecule experiments have determined the probability of loop formation in DNA as a function of the DNA contour length for different types of looping proteins. The optimal contour length for loop formation as well as the probability density functions have been found to be strongly dependent on the type of looping protein used. We show, using Monte Carlo simulations and analytical calculations, that these observations can be replicated using the wormlike-chain model for double-stranded DNA if we account for the nonzero size of the looping protein. The simulations have been performed in two dimensions so that bending is the only mode of deformation available to the DNA while the geometry of the looping protein enters through a single variable which is representative of its size. We observe two important effects that seem to directly depend on the size of the enzyme: 1), the overall propensity of loop formation at any given value of the DNA contour length increases with the size of the enzyme; and 2), the contour length corresponding to the first peak as well as the first well in the probability density functions increases with the size of the enzyme. Additionally, the eigenmodes of the fluctuating shape of the looped DNA calculated from simulations and theory are in excellent agreement, and reveal that most of the fluctuations in the DNA occur in regions of low curvature.

Biological consequences of tightly bent DNA: the other life of a macromolecular celebrity

The mechanical properties of DNA play a critical role in many biological functions. For example, DNA packing in viruses involves confining the viral genome in a volume (the viral capsid) with dimensions that are comparable to the DNA persistence length. Similarly, eukaryotic DNA is packed in DNA-protein complexes (nucleosomes), in which DNA is tightly bent around protein spools. DNA is also tightly bent by many proteins that regulate transcription, resulting in a variation in gene expression that is amenable to quantitative analysis. In these cases, DNA loops are formed with lengths that are comparable to or smaller than the DNA persistence length. The aim of this review is to describe the physical forces associated with tightly bent DNA in all of these settings and to explore the biological consequences of such bending, as increasingly accessible by single-molecule techniques.