Statistical mechanics of sequence-dependent circular DNA and its application for DNA cyclization

Stability and kinetics of G-quadruplex structures

In this review, we give an overview of recent literature on the structure and stability of unimolecular G-rich quadruplex structures that are relevant to drug design and for in vivo function. The unifying theme in this review is energetics. The thermodynamic stability of quadruplexes has not been studied in the same detail as DNA and RNA duplexes, and there are important differences in the balance of forces between these classes of folded oligonucleotides. We provide an overview of the principles of stability and where available the experimental data that report on these principles. Significant gaps in the literature have been identified, that should be filled by a systematic study of well-defined quadruplexes not only to provide the basic understanding of stability both for design purposes, but also as it relates to in vivo occurrence of quadruplexes. Techniques that are commonly applied to the determination of the structure, stability and folding are discussed in terms of information content and limitations. Quadruplex structures fold and unfold comparatively slowly, and DNA unwinding events associated with transcription and replication may be operating far from equilibrium. The kinetics of formation and resolution of quadruplexes, and methodologies are discussed in the context of stability and their possible biological occurrence.

Effects of kink and flexible hinge defects on mechanical responses of short double-stranded DNA molecules

We predict various detectable mechanical responses to the presence of local DNA defects which are defined as short DNA segments exhibiting mechanical properties obviously different from the 50 nm persistence length based semiflexible polymer model. The defects discussed are kinks and flexible hinges either permanently fixed on DNA or thermally excited. Their effects on extension shift, the effective persistence length, the end-to-end distance distribution, and the cyclization probability are computed using a transfer-matrix method. Our predictions will be useful in future experimental designs to study DNA nicks or mismatch base pairs, mechanics of specific DNA sequences, and specific DNA-protein interaction using magnetic tweezer, fluorescence resonance energy transfer, plasmon resonance techniques, and the traditional biochemistry cyclization probability measurements.

Unique Physical Signature of DNA Curvature and Its Implications for Structure and Dynamics

A particularly sensitive birefringence technique is used to analyze a curved DNA fragment with 118 bp and a standard DNA with 119 bp. At salt concentrations from 0.5 to 10 mM, both fragments show the usual negative stationary birefringence and monotonic transients - differences are relatively small. At 100 mM salt the curved DNA shows a positive stationary birefringence and non-monotonic transients with processes having amplitudes of opposite sign, whereas signals of the standard DNA remain as usual. Transients induced by reversal of the field vector indicate the existence of a permanent dipole for the curved DNA. 2-MHz-ac pulses induce a negative stationary birefringence in both DNAs. These results are consistent with calculations on models for curved DNA predicting a quasi-permanent dipole and a positive dichroism/birefringence. The quasi-permanent dipole results from the loss of symmetry in the charge distribution of the curved polyelectrolyte. The appearance of the unique signature of curvature at high salt is mainly due to a strong decrease of the polarizability by about 2 orders of magnitude. The special mode of orientation resulting from the quasi-permanent dipole is expected to contribute to the gel migration anomaly. The time constants of birefringence decay for the curved fragment are shorter than those of the 119 bp fragment by a factor of approximately 1.10 at 0.6 mM salt, whereas this factor is approximately 1.20 at 100 mM Na+. If both fragments were normal DNA with 3.4 A rise per base pair, the factor would be approximately 1.02. At high salt and high electric field strengths the factor increases up to 1.37. The implications for the bending dynamics and the potential to distinguish static from dynamic persistence by field reversal experiments are discussed. The dependence of the curvature on the salt concentration indicated by the time constants is consistent with a clear decrease of the electrophoretic anomaly at decreasing salt concentration.

Intrinsic curvature of DNA influences LacR-mediated looping

Protein-mediated DNA looping is a common mechanism for regulating gene expression. Loops occur when a protein binds to two operators on the same DNA molecule. The probability of looping is controlled, in part, by the basepair sequence of inter-operator DNA, which influences its structural properties. One structural property is the intrinsic or stress-free curvature. In this article, we explore the influence of sequence-dependent intrinsic curvature by exercising a computational rod model for the inter-operator DNA as applied to looping of the LacR-DNA complex. Starting with known sequences for the inter-operator DNA, we first compute the intrinsic curvature of the helical axis as input to the rod model. The crystal structure of the LacR (with bound operators) then defines the requisite boundary conditions needed for the dynamic rod model that predicts the energetics and topology of the intervening DNA loop. A major contribution of this model is its ability to predict a broad range of published experimental data for highly bent (designed) sequences. The model successfully predicts the loop topologies known from fluorescence resonance energy transfer measurements, the linking number distribution known from cyclization assays with the LacR-DNA complex, the relative loop stability known from competition assays, and the relative loop size known from gel mobility assays. In addition, the computations reveal that highly curved sequences tend to lower the energetic cost of loop formation, widen the energy distribution among stable and meta-stable looped states, and substantially alter loop topology. The inclusion of sequence-dependent intrinsic curvature also leads to nonuniform twist and necessitates consideration of eight distinct binding topologies from the known crystal structure of the LacR-DNA complex.

Nuclear protein‐induced bending and flexing of the hypoxic response element of the rat vascular endothelial growth factor promoter

Bending and flexing of DNA may contribute to transcriptional regulation. Because hypoxia and other physiological signals induce formation of an abasic site at a key base within the hypoxic response element (HRE) of the vascular endothelial growth factor (VEGF) gene (FASEB J. 19, 387-394, 2005) and because abasic sites can introduce flexibility in model DNA sequences, in the present study we used a fluorescence resonance energy transfer-based reporter system to assess topological changes in a wild-type (WT) sequence of the HRE of the rat VEGF gene and in a sequence harboring a single abasic site mimicking the effect of hypoxia. Binding of the hypoxia-inducible transcriptional complex present in hypoxic pulmonary artery endothelial cell nuclear extract to the WT sequence failed to alter sequence topology whereas nuclear protein binding to the modified HRE engendered considerable sequence flexibility. Topological effects of nuclear proteins on the modified VEGF HRE were dependent on the transcription factor hypoxia-inducible factor-1 and on formation of a single-strand break at the abasic site mediated by the coactivator, Ref-1/Ape1. These observations suggest that oxidative base modifications in the VEGF HRE evoked by physiological signals could be a precursor to single-strand break formation that has an impact on gene expression by modulating sequence flexibility.

T4 DNA ligase is more than an effective trap of cyclized dsDNA

T4 DNA ligase is used in standard cyclization assays to trap double-stranded DNA (dsDNA) in low-probability, cyclic or highly bent conformations. The cyclization probability, deduced from the relative yield of cyclized product, can be used in conjunction with statistical mechanical models to extract the bending stiffness of dsDNA. By inserting the base analog 2-aminopurine (2-AP) at designated positions in 89 bp and 94 bp dsDNA fragments, we find that T4 DNA ligase can have a previously unknown effect. Specifically, we observe that addition of T4 ligase to dsDNA in proportions comparable to what is used in the cyclization assay leads to a significant increase in fluorescence from 2-AP. This effect is believed to originate from stabilization of local base-pair opening by formation of transient DNA-ligase complexes. Non-specific binding of T4 ligase to dsDNA is also confirmed using fluorescence correlation spectroscopy (FCS) experiments, which reveal a systematic reduction of dsDNA diffusivity in the presence of ligase. ATP competes with regular DNA for non-covalent binding to the T4 ligase and is found to significantly reduce DNA-ligase complexation. For short dsDNA fragments, however, the population of DNA-ligase complexes at typical ATP concentrations used in DNA cyclization studies is determined to be large enough to dominate the cyclization reaction.

Kinking occurs during molecular dynamics simulations of small DNA minicircles

Recent experiments on minicircle formation suggest that a conformational mechanism other than smooth deformation may be playing a role in enhancing DNA flexibility. Both local base unpairing and kink formation have been suggested as possible explanations. Although kinks within isolated DNA were proposed 30 years ago, they have, until now, only been observed within DNA complexed with proteins. In order to test how DNA behaves in the strong bending regime, we have carried out molecular dynamics simulations of a 94 base pair minicircle in explicit solvent with two different linking numbers, corresponding to a torsionally relaxed state and a positively supercoiled state. The simulations suggest that sharp kinks can indeed arise in small minicircles. The relaxed minicircle is generally associated with a single kink, while two kinks occur with the supercoiled state. No evidence is seen of base unpaired regions.

Analysis of in-vivo LacR-mediated gene repression based on the mechanics of DNA looping

Interactions of E. coli lac repressor (LacR) with a pair of operator sites on the same DNA molecule can lead to the formation of looped nucleoprotein complexes both in vitro and in vivo. As a major paradigm for loop-mediated gene regulation, parameters such as operator affinity and spacing, repressor concentration, and DNA bending induced by specific or non-specific DNA-binding proteins (e.g., HU), have been examined extensively. However, a complete and rigorous model that integrates all of these aspects in a systematic and quantitative treatment of experimental data has not been available. Applying our recent statistical-mechanical theory for DNA looping, we calculated repression as a function of operator spacing (58-156 bp) from first principles and obtained excellent agreement with independent sets of in-vivo data. The results suggest that a linear extended, as opposed to a closed v-shaped, LacR conformation is the dominant form of the tetramer in vivo. Moreover, loop-mediated repression in wild-type E. coli strains is facilitated by decreased DNA rigidity and high levels of flexibility in the LacR tetramer. In contrast, repression data for strains lacking HU gave a near-normal value of the DNA persistence length. These findings underscore the importance of both protein conformation and elasticity in the formation of small DNA loops widely observed in vivo, and demonstrate the utility of quantitatively analyzing gene regulation based on the mechanics of nucleoprotein complexes.

Torsional rigidities of weakly strained DNAs

Measurements on unstrained linear and weakly strained large (> or =340 bp) circular DNAs yield torsional rigidities in the range C = 170-230 fJ fm. However, larger values, in the range C = 270-420 fJ fm, are typically obtained from measurements on sufficiently small (< or =247 bp) circular DNAs, and values in the range C = 300-450 fJ fm are obtained from experiments on linear DNAs under tension. A new method is proposed to estimate torsional rigidities of weakly supercoiled circular DNAs. Monte Carlo simulations of the supercoiling free energies of solution DNAs, and also of the structures of surface-confined supercoiled plasmids, were performed using different trial values of C. The results are compared with experimental measurements of the twist energy parameter, E(T), that governs the supercoiling free energy, and also with atomic force microscopy images of surface-confined plasmids. The results clearly demonstrate that C-values in the range 170-230 fJ fm are compatible with experimental observations, whereas values in the range C > or = 269 fJ fm, are incompatible with those same measurements. These results strongly suggest that the secondary structure of DNA is altered by either sufficient coherent bending strain or sufficient tension so as to enhance its torsional rigidity.

Effect of supercoiling on formation of protein-mediated DNA loops

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short (50-100 basepair), protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay.

Spontaneous sharp bending of DNA: role of melting bubbles

The role of centrally located and distributed base pair mismatches (‘melting bubbles’) on localized bending and stiffness of short dsDNA fragments is evaluated using time-dependent fluorescence lifetime measurements. Distributed melting bubbles are found to induce larger bending angles and decreased levels of stiffness in DNA than centrally located ones of comparable overall size. Our results indicate that spontaneous local opening-up of the DNA duplex could facilitate sharp bending of short DNA strands even in the absence of DNA binding proteins. We also find that the occurrence of two closely spaced melting bubbles will generally be favored when a large energetic barrier must be overcome in forming the desired bent DNA structure.

A computational study of nucleosomal DNA flexibility

Molecular dynamics simulations of the nucleosome core particle and its isolated DNA free in solution are reported. The simulations are based on the implicit solvent methodology and provide insights into the nature of large-scale structural fluctuations and flexibility of the nucleosomal DNA. In addition to the kinked regions previously identified in the x-ray structure of the nucleosome, the simulations support the existence of a biochemically identified distorted region of the DNA. Comparison of computed relative free energies shows that formation of the kinks is associated with little, if any, energy cost relative to a smooth, ideal conformation of the DNA superhelix. Isolated nucleosomal DNA is found to be considerably more flexible than expected for a 147 bp stretch of DNA based on its canonical persistence length of 500 A. Notably, the significant bending of the DNA observed in our simulations occurs without breaking of Watson-Crick bonds. The computed relative stability of bent conformations is sensitive to the ionic strength of the solution in the physiological range; the sensitivity suggests possible experiments that might provide further insights into the structural origins of the unusual flexibility of the DNA.

Modeling the Lac repressor-operator assembly: the influence of DNA looping on Lac repressor conformation

Repression of transcription of the Escherichia coli Lac operon by the Lac repressor (LacR) is accompanied by the simultaneous binding of LacR to two operators and the formation of a DNA loop. A recently developed theory of sequence-dependent DNA elasticity enables one to relate the fine structure of the LacR-DNA complex to a wide range of heretofore-unconnected experimental observations. Here, that theory is used to calculate the configuration and free energy of the DNA loop as a function of its length and base-pair sequence, its linking number, and the end conditions imposed by the LacR tetramer. The tetramer can assume two types of conformations. Whereas a rigid V-shaped structure is observed in the crystal, EM images show extended forms in which two dimer subunits are flexibly joined. Upon comparing our computed loop configurations with published experimental observations of permanganate sensitivities, DNase I cutting patterns, and loop stabilities, we conclude that linear DNA segments of short-to-medium chain length (50-180 bp) give rise to loops with the extended form of LacR and that loops formed within negatively supercoiled plasmids induce the V-shaped structure.

DNA looping by two-site restriction endonucleases: heterogeneous probability distributions for loop size and unbinding force

Proteins interacting at multiple sites on DNA via looping play an important role in many fundamental biochemical processes. Restriction endonucleases that must bind at two recognition sites for efficient activity are a useful model system for studying such interactions. Here we used single DNA manipulation to study sixteen known or suspected two-site endonucleases. In eleven cases (BpmI, BsgI, BspMI, Cfr10I, Eco57I, EcoRII, FokI, HpaII, NarI, Sau3AI and SgrAI) we found that substitution of Ca2+ for Mg2+ blocked cleavage and enabled us to observe stable DNA looping. Forced disruption of these loops allowed us to measure the frequency of looping and probability distributions for loop size and unbinding force for each enzyme. In four cases we observed bimodal unbinding force distributions, indicating conformational heterogeneity and/or complex binding energy landscapes. Measured unlooping events ranged in size from 7 to 7500 bp and the most probable size ranged from less than 75 bp to nearly 500 bp, depending on the enzyme. In most cases the size distributions were in much closer agreement with theoretical models that postulate sharp DNA kinking than with classical models of DNA elasticity. Our findings indicate that DNA looping is highly variable depending on the specific protein and does not depend solely on the mechanical properties of DNA.

Do femtonewton forces affect genetic function? A review

Protein-Mediated DNA looping is intricately related to gene expression. Therefore any mechanical constraint that disrupts loop formation can play a significant role in gene regulation. Polymer physics models predict that less than a piconewton of force may be sufficient to prevent the formation of DNA loops. Thus, it appears that tension can act as a molecular switch that controls the much larger forces associated with the processive motion of RNA polymerase. Since RNAP can exert forces over 20 pN before it stalls, a 'substrate tension switch' could offer a force advantage of two orders of magnitude. Evidence for such a mechanism is seen in recent in vitro micromanipulation experiments. In this article we provide new perspective on existing theory and experimental data on DNA looping in vitro and in vivo. We elaborate on the connection between tension and a variety of other intracellular mechanical constraints including sequence specific curvature and supercoiling. In the process, we emphasize that the richness and versatility of DNA mechanics opens up a whole new paradigm of gene regulation to explore.

Formation and positioning of nucleosomes: effect of sequence-dependent long-range correlated structural disorder

The understanding of the long-range correlations (LRC) observed in DNA sequences is still an open and very challenging problem. In this paper, we start reviewing recent results obtained when exploring the scaling properties of eucaryotic, eubacterial and archaeal genomic sequences using the space-scale decomposition provided by the wavelet transform (WT). These results suggest that the existence of LRC up to distances approximately 20-30 kbp is the signature of the nucleosomal structure and dynamics of the chromatin fiber. Actually the LRC are mainly observed in the DNA bending profiles obtained when using some structural coding of the DNA sequences that accounts for the fluctuations of the local double-helix curvature within the nucleosome complex. Because of the approximate planarity of nucleosomal DNA loops, we then study the influence of the LRC structural disorder on the thermodynamical properties of 2D elastic chains submitted locally to mechanical/topological constraint as loops. The equilibrium properties of the one-loop system are derived numerically and analytically in the quite realistic weak-disorder limit. The LRC are shown to favor the spontaneous formation of small loops, the larger the LRC, the smaller the size of the loop. We further investigate the dynamical behavior of such a loop using the mean first passage time (MFPT) formalism. We show that the typical short-time loop dynamics is superdiffusive in the presence of LRC. For displacements larger than the loop size, we use large-deviation theory to derive a LRC-dependent anomalous-diffusion rule that accounts for the lack of disorder self-averaging. Potential biological implications on DNA loops involved in nucleosome positioning and dynamics in eucaryotic chromatin are discussed.

Modeling DNA loops using the theory of elasticity

An elastic rod model of a protein-bound DNA loop is adapted for application in multi-scale simulations of protein-DNA complexes. The classical Kirchhoff system of equations which describes the equilibrium structure of the elastic loop is modified to account for the intrinsic twist and curvature, anisotropic bending properties, and electrostatic charge of DNA. The effects of bending anisotropy and electrostatics are studied for the DNA loop clamped by the lac repressor protein. For two possible lengths of the loop, several topologically different conformations are predicted and extensively analyzed over the broad range of model parameters describing DNA bending and electrostatic properties. The scope and applications of the model in already accomplished and in future multi-scale studies of protein-DNA complexes are discussed.

Statistical-mechanical theory of DNA looping

The lack of a rigorous analytical theory for DNA looping has caused many DNA-loop-mediated phenomena to be interpreted using theories describing the related process of DNA cyclization. However, distinctions in the mechanics of DNA looping versus cyclization can have profound quantitative effects on the thermodynamics of loop closure. We have extended a statistical mechanical theory recently developed for DNA cyclization to model DNA looping, taking into account protein flexibility. Notwithstanding the underlying theoretical similarity, we find that the topological constraint of loop closure leads to the coexistence of multiple classes of loops mediated by the same protein structure. These loop topologies are characterized by dramatic differences in twist and writhe; because of the strong coupling of twist and writhe within a loop, DNA looping can exhibit a complex overall helical dependence in terms of amplitude, phase, and deviations from uniform helical periodicity. Moreover, the DNA-length dependence of optimal looping efficiency depends on protein elasticity, protein geometry, and the presence of intrinsic DNA bends. We derive a rigorous theory of loop formation that connects global mechanical and geometric properties of both DNA and protein and demonstrates the importance of protein flexibility in loop-mediated protein-DNA interactions.

End-to-end distance vector distribution with fixed end orientations for the wormlike chain model

We find exact expressions for the end-to-end distance vector distribution function with fixed end orientations for the wormlike chain model. This function in Fourier-Laplace space adopts the form of infinite continued fractions, which emerges upon exploiting the hierarchical structure of the moment-based expansion. Our results are used to calculate the root-mean-square end displacement in a given direction for a chain with both end orientations fixed. We find that the crossover from rigid to flexible chains is marked by the root-mean-square end displacement slowly losing its angular dependence as the coupling between chain conformation and end orientation wanes. However, the coupling remains strong even for relatively flexible chains, suggesting that the end orientation strongly influences chain conformation for chains that are several persistence lengths long. We then show the behavior of the distribution function by a density plot of the probability as a function of the end-to-end distance vector for a wormlike chain in two dimensions with one end pointed in a fixed direction and the other end free (in its orientation). As we progress from high to low rigidity, the distribution function shifts from being peaked at a location near the full contour length of the chain in the forward direction, corresponding to a straight configuration, to being peaked near zero end separation, as in the Gaussian limit. The function exhibits double peaks in the crossover between these limiting behaviors.

DNA dynamics: bubble 'n' flip for DNA cyclisation?

A recent demonstration of the facile in vitro formation of DNA microcircles of fewer than 100 base pairs throws new light on the basis of DNA flexibility.

Thermodynamics of DNA loops with long-range correlated structural disorder

We study the influence of a structural disorder on the thermodynamical properties of 2D-elastic chains submitted to mechanical/topological constraint as loops. The disorder is introduced via a spontaneous curvature whose distribution along the chain presents either no correlation or long-range correlations (LRC). The equilibrium properties of the one-loop system are derived numerically and analytically for weak disorder. LRC are shown to favor the formation of small loop, larger the LRC, smaller the loop size. We use the mean first passage time formalism to show that the typical short time loop dynamics is superdiffusive in the presence of LRC. Potential biological implications on nucleosome positioning and dynamics in eukaryotic chromatin are discussed.

Statistics of loop formation along double helix DNAs

We compute relative position distributions of distant sites along discretized semiflexible polymers, focusing on encounter statistics for pairs of sites along a double-stranded DNA molecule (dsDNA), using a transfer-matrix approach. We generalize the usual semiflexible polymer, considering nonlinear elasticity effects arising from inhomogeneities which either appear at any position via thermal fluctuation, or which occur at specific "quenched" locations. We apply our theory to two problems associated with dsDNA looping. First, we discuss how local flexible defects in double-helix structure facilitate cyclization of short dsDNA molecules. Flexible defects greatly enhance cyclization rate, and strongly modify its dependence on the closure orientational boundary condition. This effect is relevant to free-solution cyclization experiments, and to loop formation in vivo. Second, we present calculations of force dependence of the probability of formation of loops along single dsDNAs which show how the probability of loop formation is suppressed by tension.

Cyclization of short DNA fragments and bending fluctuations of the double helix

Cloutier and Widom [Cloutier, T. E. & Widom, J. (2004) Mol. Cell 14, 355-362] recently reported that the cyclization efficiency of short DNA fragments, about 100 bp in length, exceeds theoretical expectations by three orders of magnitude. In an effort to resolve this discrepancy, we tried modifying the theory. We investigated how the distribution of the angles between adjacent base pairs of the double helix affects the cyclization efficiency. We found that only the incorporation of sharp kinks in the angle distribution provides the desired increase of the cyclization efficiency. We did not find a model, however, that fits all cyclization data for DNA fragments of different lengths. Therefore, we carefully reinvestigated the cyclization of 100-bp DNA fragments experimentally and found their cyclization efficiency to be in remarkable agreement with the traditional model of DNA bending. We also found an explanation for the discrepancy between our results and those of Cloutier and Widom.

DNA twisting flexibility and the formation of sharply looped protein-DNA complexes

Gene-regulatory complexes often require that pairs of DNA-bound proteins interact by looping-out short (often approximately 100-bp) stretches of DNA. The loops can vary in detailed length and sequence and, thus, in total helical twist, which radically alters their geometry. How this variability is accommodated structurally is not known. Here we show that the inherent twistability of 89- to 105-bp DNA circles exceeds theoretical expectation by up to 400-fold. These results can be explained only by greatly enhanced DNA flexibility, not by permanent bends. They invalidate the use of classic theories of flexibility for understanding sharp DNA looping but support predictions of two recent theories. Our findings imply an active role for DNA flexibility in loop formation and suggest that variability in the detailed helical twist of regulatory loops is accommodated naturally by the inherent twistability of the DNA.

Exact theory of kinkable elastic polymers

The importance of nonlinearities in material constitutive relations has long been appreciated in the continuum mechanics of macroscopic rods. Although the moment (torque) response to bending is almost universally linear for small deflection angles, many rod systems exhibit a high-curvature softening. The signature behavior of these rod systems is a kinking transition in which the bending is localized. Recent DNA cyclization experiments by Cloutier and Widom have offered evidence that the linear-elastic bending theory fails to describe the high-curvature mechanics of DNA. Motivated by this recent experimental work, we develop a simple and exact theory of the statistical mechanics of linear-elastic polymer chains that can undergo a kinking transition. We characterize the kinking behavior with a single parameter and show that the resulting theory reproduces both the low-curvature linear-elastic behavior which is already well described by the worm-like chain model, as well as the high-curvature softening observed in recent cyclization experiments.

Disruption of protein-mediated DNA looping by tension in the substrate DNA

Protein-mediated DNA looping is important in a variety of biological processes, including gene regulation and genetic transformation. Although the biochemistry of loop formation is well established, the mechanics of loop closure in a constrained cellular environment has received less attention. Recent single molecule measurements show that mechanical constraints have a significant impact on DNA looping and motivate the need for a more comprehensive characterization of the effects of tension. By modeling DNA as a wormlike chain, we calculate how continuous stretching of the substrate DNA affects the loop formation probability. We find that when the loop size is >100 bp, a tension of 500 fN can increase the time required for loop closure by two orders of magnitude. This force is small compared to the piconewton forces that are associated with RNA polymerases and other molecular motors, indicating that intracellular mechanical forces might affect transcriptional regulation. In contrast to existing theory, we find that for loops <200 bp, the effect of tension is partly dependent on the relative orientation of the DNA-binding domains in the linker protein. Our results provide perspective on recent DNA looping experiments and suggestions for future micromechanical studies.

DNA Codes for Nanoscience

The nanometer scale is a special place where all sciences meet and develop a particularly strong interdisciplinarity. While biology is a source of inspiration for nanoscientists, chemistry has a central role in turning inspirations and methods from biological systems to nanotechnological use. DNA is the biological molecule by which nanoscience and nanotechnology is mostly fascinated. Nature uses DNA not only as a repository of the genetic information, but also as a controller of the expression of the genes it contains. Thus, there are codes embedded in the DNA sequence that serve to control recognition processes on the atomic scale, such as the base pairing, and others that control processes taking place on the nanoscale. From the chemical point of view, DNA is the supramolecular building block with the highest informational content. Nanoscience has therefore the opportunity of using DNA molecules to increase the level of complexity and efficiency in self-assembling and self-directing processes.

Localized single-stranded bubble mechanism for cyclization of short double helix DNA

Recent experiments indicate that double-stranded DNA molecules of approximately 100 base pairs in length have a probability of cyclization which is up to 10(5) times larger than that expected based on the known bending modulus of the double-helix. We argue that for short molecules, the formation of a few base pairs of single-stranded DNA can provide a "flexible hinge" that facilitates loop formation. A detailed calculation shows that this mechanism explains the experimental data.

Spontaneous sharp bending of double-stranded DNA

Sharply bent DNA is essential for gene regulation in prokaryotes and is a major feature of eukaryotic nucleosomes and viruses. The explanation normally given for these phenomena is that specific proteins sharply bend DNA by application of large forces, while the DNA follows despite its intrinsic inflexibility. Here we show that DNAs that are 94 bp in length-comparable to sharply looped DNAs in vivo-spontaneously bend into circles. Proteins can enhance the stability of such loops, but the loops occur spontaneously even in naked DNA. Random DNA sequences cyclize 10(2)-10(4) times more easily than predicted from current theories of DNA bending, while DNA sequences that position nucleosomes cyclize up to 10(5) times more easily. These unexpected results establish DNA as an active participant in the formation of looped regulatory complexes in vivo, and they point to a need for new theories of DNA bending.

Predicting indirect readout effects in protein-DNA interactions

Recognition of DNA by proteins relies on direct interactions with specific DNA-functional groups, along with indirect effects that reflect variable energetics in the response of DNA sequences to twisting and bending distortions induced by proteins. Predicting indirect readout requires knowledge of the variations in DNA curvature and flexibility in the affected region, which we have determined for a series of DNA-binding sites for the E2 regulatory protein by using the cyclization kinetics method. We examined 16 sites containing different noncontacted spacer sequences, which vary by more than three orders of magnitude in binding affinity. For 15 of these sites, the variation in affinity was predicted within a factor of 3, by using experimental curvature and flexibility values and a statistical mechanical theory. The sole exception was traced to differential magnesium ion binding.

DNA sequence-dependent deformability?insights from computer simulations

The article reviews some recent developments in studying DNA sequence-dependent deformability, with emphasis on computer modeling. After a brief outline of available experimental techniques, we proceed to computational methods and focus on atomic-resolution molecular dynamics (MD) simulations. A sequence-dependent local (base-pair step) force field inferred from MD is compared with force fields obtained by other techniques. Various methods for establishing global (flexible-rod) DNA elastic constants are reviewed, including an approach based on atomic resolution MD. The problem of defining the global deformation variables, as well as the question of anisotropy and nonlocal effects, are discussed. As an example, both local and global deformability calculations from atomic-resolution MD of EcoRI dodecamer are presented.

Link, twist, energy, and the stability of DNA minicircles

We describe how the stability properties of DNA minicircles can be directly read from plots of various biologically intuitive quantities along families of equilibrium configurations. Our conclusions follow from extensions of the mathematical theory of distinguished bifurcation diagrams that are applied within the specific context of an elastic rod model of minicircles. Families of equilibria arise as a twisting angle alpha is varied. This angle is intimately related to the continuously varying linking number Lk for nicked DNA configurations that is defined as the sum of Twist and Writhe. We present several examples of such distinguished bifurcation diagrams involving plots of the energy E, linking number Lk, and a twist moment m3, along families of cyclized equilibria of both intrinsically straight and intrinsically curved DNA fragments.

Comparison of analyses of DNA curvature

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.

High-throughput approach for detection of DNA bending and flexibility based on cyclization

We have developed a high-throughput approach to the labor-intensive problems of DNA cyclization, which we use to characterize DNA curvature and mechanical properties. The method includes a combinatorial approach to make the DNA constructs needed and automated real-time measurement of the kinetics using fluorescence. We validated the approach and investigated the flexibility of two kinds of nicked DNA and AT dinucleotide repeats. We found that, although the nicks hardly alter the bending flexibility, they significantly increase the torsional flexibility, and that the AT repeat has 28% (+/-12%) lower bending rigidity than a generic DNA sequence.