Effect of anisotropy of the bending rigidity on the supercoiling free energy of small circular DNAs

Strong bending of the DNA double helix

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10⁻⁵, it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.

Estimation of the persistence length of DNA from the torsion elastic constant and supercoiling free energy: effect of ethylene glycol

Two different methods are proposed to estimate the persistence length ( P) of DNA from the measured torsion elastic constant (alpha) and the twist energy parameter ( E T ) that governs the supercoiling free energy. The first method involves Monte Carlo simulations and reversible-work calculations of E T for model DNAs that possess the measured alpha and selected trial values of P. Comparison of the computed E T values with the experimental value allows estimation of P (or equivalently the bending elastic constant (kappa beta)) by interpolation. A far simpler, though less accurate, alternative is to solve a previously conjectured analytical relation connecting E T , alpha, kappa beta (or P), and an unknown "constant" ( B). The present simulations are used to ascertain the optimum value of B and to assess the validity and accuracy of that relation. Within the simulation errors, P values obtained from the measured alpha and E T via this analytical expression agree with those determined from the simulations and E T values reckoned from the input alpha and kappa beta by this analytical expression agree with the corresponding simulated values. Although B is found to be insensitive to variation in alpha, it appears to decline slightly with increasing kappa beta. The original analytical expression is modified to take this apparent variation of B with kappa beta into account. By using this modified analytical relation to estimate P (from the measured alpha and E T ) or E T (from the input alpha and kappa beta), much closer agreement is obtained respectively with the values of P or E T obtained from the simulations. As specific examples, these methods are applied to determine P in 0 and 20 w/v % ethylene glycol, which has been shown to induce a structural transition in duplex DNA.

An i-motif-containing DNA device that breaks certain forms of Watson-Crick interactions

An i-motif-containing DNA device is constructed that is able to break certain forms of Watson-Crick interactions under isothermal conditions.

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.

Monte Carlo simulations of locally melted supercoiled DNAs in 20 mM ionic strength

Mesoscopic models of unmelted and locally melted supercoiled DNAs in 20 mM ionic strength are simulated over a range of linking difference from deltal = 0 to -26 turns, or superhelix density from sigma = 0 to -0.062. A domain containing m = 0, 28, or 56 melted basepairs (out of 4349 total) is modeled simply by a region of suitable length with substantially reduced torsion and bending elastic constants. Average structural properties are calculated from the saved configurations, and a reversible work protocol is used to calculate the supercoiling free energy, The cross-writhe between duplex and melted regions (defined herein) is found to be negligibly small. The total writhe, radius of gyration, and ordered elements of the diagonalized inertial tensor are found to be nearly universal functions of the residual linking difference (deltal(r)) associated with the duplex region, independent of m. However, deformability of the tertiary structure, as manifested by the variance of those same properties, is not a universal function of deltal(r)), but depends upon m.delta (SC) varies with deltal(r)) more strongly than deltal(r)) (2)due to the low ionic strength. The twist energy parameter, E (T) obtained from the simulated delta G(SC), deltal(r)), and net twisting strain of the melted region T (D), is found to be independent of m, hence also of the torsion and bending elastic constants of the melted region. However, E(T) increases linearly with -deltalr), which leads to 1). a small overestimation of E (T) for any given deltal(r)) when E(T) is determined from the observed deltal and deltal (r) by the protocol of Bauer and Benham; and 2). a significant enhancement of the apparent slope, -dE(T)/d(T), obtained via the protocol of Bauer and Benham, relative to the actual slope at fixed delta l(r). After taking these two effects into account, the theoretical and experimental values E(T) and -dE(T)/d(T) values agree rather well. For the larger deltal the melted regions are found preferentially in the linker domains between interwound arms, rather than in the apical regions at the ends of interwound arms.

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

DNA cyclization is potentially the most powerful approach for systematic quantitation of sequence-dependent DNA bending and flexibility. We extend the statistical mechanics of the homogeneous DNA circle to a model that considers discrete basepairs, thus allowing for inhomogeneity, and apply the model to analysis of DNA cyclization. The theory starts from an iterative search for the minimum energy configuration of circular DNA. Thermodynamic quantities such as the J factor, which is essentially the ratio of the partition functions of circular and linear forms, are evaluated by integrating the thermal fluctuations around the configuration under harmonic approximation. Accurate analytic expressions are obtained for equilibrium configurations of homogeneous circular DNA with and without bending anisotropy. J factors for both homogeneous and inhomogeneous DNA are evaluated. Effects of curvature, helical repeat, and bending and torsional flexibility in DNA cyclization are analyzed in detail, revealing that DNA cyclization can detect as little as one degree of curvature and a few percent change in flexibility. J factors calculated by our new approach are well consistent with Monte Carlo simulations, whereas the new theory has much greater efficiency in computations. Simulation of experimental results has been demonstrated.

Monte Carlo simulations of supercoiled DNAs confined to a plane

Recent advances in atomic force microscopy (AFM) have enabled researchers to obtain images of supercoiled DNAs deposited on mica surfaces in buffered aqueous milieux. Confining a supercoiled DNA to a plane greatly restricts its configurational freedom, and could conceivably alter certain structural properties, such as its twist and writhe. A program that was originally written to perform Monte Carlo simulations of supercoiled DNAs in solution was modified to include a surface potential. This potential flattens the DNAs to simulate the effect of deposition on a surface. We have simulated transfers of a 3760-basepair supercoiled DNA from solution to a surface in both 161 and 10 mM ionic strength. In both cases, the geometric and thermodynamic properties of the supercoiled DNAs on the surface differ significantly from the corresponding quantities in solution. At 161 mM ionic strength, the writhe/twist ratio is 1.20-1.33 times larger for DNAs on the surface than for DNAs in solution and significant differences in the radii of gyration are also observed. Simulated surface structures in 161 mM ionic strength closely resemble those observed by AFM. Simulated surface structures in 10 mM ionic strength are similar to a minority of the structures observed by AFM, but differ from the majority of such structures for unknown reasons. In 161 mM ionic strength, the internal energy (excluding the surface potential) decreases substantially as the DNA is confined to the surface. Evidently, supercoiled DNAs in solution are typically deformed farther from the minimum energy configuration than are the corresponding surface-confined DNAs. Nevertheless, the work (Delta A(int)) done on the internal coordinates, which include uniform rotations at constant configuration, during the transfer is positive and 2.6-fold larger than the decrease in internal energy. The corresponding entropy change is negative, and its contribution to Delta A(int) is positive and exceeds the decrease in internal energy by 3.6 fold. The work done on the internal coordinates during the solution-to-surface transfer is directed primarily toward reducing their entropy. Evidently, the number of configurations available to the more deformed solution DNA is vastly greater than for the less deformed surface-confined DNA.

The question of long-range allosteric transitions in DNA

The question of long-range allosteric transitions of DNA secondary structure and their possible involvement in transcriptional activation is discussed in the light of new results. A variety of recent evidence strongly supports a fluctuating long-range description of DNA secondary structure. Balanced equilibria between two or more different secondary structures, and the occurrence of very large domain sizes, have been documented in several instances. Long-range allosteric effects stemming from changes in sequence or secondary structure over a small region of the DNA have been observed to extend over distances up to hundreds of base pairs in some cases. The discovery that coherent bending strain beyond a threshold level in small (N < or = 250 base pairs (bp)] circular DNAs significantly alters the DNA secondary structure has important implications, especially for transcriptional activators that either bend the DNA directly or are involved in the formation of DNA loops of sufficiently small size (N < or = 250 bp). Whether the RNA polymerase is activated primarily via protein: protein contacts, as is widely believed, or instead via a bend-induced allosteric transition of the DNA in such a small loop, is now an open question. Binding of the transcriptional activator Sp1 to linear DNA induces a remarkably long-range change in its secondary structure, and catabolite activator protein binding to a supercoiled DNA behaves similarly, though possibly for different reasons. Compelling evidence for a bend-induced long-range structural transmission effect of the transcriptional activator integration host factor on RNA polymerase activity was recently reported. These results may augur a new paradigm in which allosteric transitions of duplex DNA, as well as of the proteins, are involved in the regulation of transcription.

Dynamic bending rigidity of a 200-bp DNA in 4 mM ionic strength: a transient polarization grating study

DNA may exhibit three different kinds of bends: 1) permanent bends; 2) slowly relaxing bends due to fluctuations in a prevailing equilibrium between differently curved secondary conformations; and 3) rapidly relaxing dynamic bends within a single potential-of-mean-force basin. The dynamic bending rigidity (kappa(d)), or equivalently the dynamic persistence length, P(d) = kappa(d)/k(B)T, governs the rapidly relaxing bends, which are responsible for the flexural dynamics of DNA on a short time scale, t < or = 10(-5) s. However, all three kinds of bends contribute to the total equilibrium persistence length, P(tot), according to 1/P(tot) congruent with 1/P(pb) + 1/P(sr) + 1/P(d), where P(pb) is the contribution of the permanent bends and P(sr) is the contribution of the slowly relaxing bends. Both P(d) and P(tot) are determined for the same 200-bp DNA in 4 mM ionic strength by measuring its optical anisotropy, r(t), from 0 to 10 micros. Time-resolved fluorescence polarization anisotropy (FPA) measurements yield r(t) for DNA/ethidium complexes (1 dye/200 bp) from 0 to 120 ns. A new transient polarization grating (TPG) experiment provides r(t) for DNA/methylene blue complexes (1 dye/100 bp) over a much longer time span, from 20 ns to 10 micros. Accurate data in the very tail of the decay enable a model-independent determination of the relaxation time (tau(R)) of the end-over-end tumbling motion, from which P(tot) = 500 A is estimated. The FPA data are used to obtain the best-fit pairs of P(d) and torsion elastic constant (alpha) values that fit those data equally well, and which are used to eliminate alpha as an independent variable. When the relevant theory is fitted to the entire TPG signal (S(t)), the end-over-end rotational diffusion coefficient is fixed at its measured value and alpha is eliminated in favor of P(d). Neither a true minimum in chi-squared nor a satisfactory fit could be obtained for P(d) anywhere in the range 500-5000 A, unless an adjustable amplitude of azimuthal wobble of the methylene blue was admitted. In that case, a well-defined global minimum and a reasonably good fit emerged at P(d) = 2000 A and <deltazeta(2)>(1/2) = 25 degrees. The discrimination against P(d) values <1600 A is very great. By combining the values, P(tot) = 500 A and P(d) = 2000 A with a literature estimate, P(pb) = 1370 A, a value P(sr) = 1300 A is estimated for the contribution of slowly relaxing bends. This value is analyzed in terms of a simple model in which the DNA is divided up into domains containing m bp, each of which experiences an all-or-none equilibrium between a straight and a uniformly curved conformation. With an appropriate estimate of the average bend angle per basepair of the curved conformation, a lower bound estimate, m = 55 bp, is obtained for the domain size of the coherently bent state. Previous measurements suggest that this coherent bend is not directional, or phase-locked, to the azimuthal orientation of the filament.

Comparison of hard-cylinder and screened Coulomb interactions in the modeling of supercoiled DNAs

A 1000 base pair (bp) model supercoiled DNA is simulated using spherical screened Coulomb interactions between subunits on one hand and equivalent hard-cylinder interactions on the other. The amplitudes, or effective charges, of the spherical screened Coulomb electrostatic potentials are chosen so that the electrostatic potential surrounding the middle of a linear array of 2001 subunits (31.8 A diameter) closely matches the solution of the nonlinear Poisson-Boltzmann equation for a cylinder with 12 A radius and the full linear charge density of DNA at all distances beyond the 24 A hard-core diameter. This superposition of spherical screened Coulomb potentials is practically identical to the particular solution of the cylindrical linearized Poisson-Boltzmann equation that matches the solution of the nonlinear Poisson-Boltzmann equation at large distances. The interaction energy between subunits is reckoned from the effective charges according to the standard DLVO expression. The equivalent hard-cylinder diameter is chosen following Stigter's protocol for matching second virial coefficients, but for the full linear charge density of DNA. The electrostatic persistence length of the model with screened Coulomb interactions is extremely sensitive to the (arbitrarily) chosen subunit length at the higher salt concentrations. The persistence length of the hard-cylinder model is adjusted to match that of the screened Coulomb model for each ionic condition. Simulations for a superhelix density sigma = -0.05 using a spherical screened Coulomb interaction plus a 24 A hard-cylinder core (SCPHC) potential indicate that the radius of gyration of this 1000 bp DNA actually undergoes a slight increase as the NaCl concentration is raised from 0.01 to 1.0M. Thus, merely softening the potential from hard-cylinder to screened Coulomb form does not produce a large decrease in radius of gyration with increasing NaCl concentration for DNAs of this size. Radii of gyration, static structure factors, and diffusion coefficients obtained using the equivalent hard-cylinder (EHC) potential agree well with those obtained using the SCPHC potential in 1.0M NaCl, but in 0.1M NaCl the agreement is not as good, and in 0.01M NaCl the agreement is definitely unsatisfactory. These conclusions differ in significant respects from those obtained in previous studies.

Configurational transitions in Fourier series-represented DNA supercoils

A new Fourier series representation of supercoiled DNA is employed in Langevin dynamics simulations to study large-scale configurational motions of intermediate-length chains. The polymer is modeled as an ideal elastic rod subject to long-range van der Waals' interactions. The van der Waals' term prevents the self-contact of distant chain segments and also mimics attractive forces thought to stabilize the association of closely spaced charged rods. The finite Fourier series-derived polymer formulation is an alternative to the piecewise B-spline curves used in past work to describe the motion of smoothly deformed supercoiled DNA in terms of a limited number of independent variables. This study focuses on two large-scale configurational events: the interconversion between circular and figure-8 forms at a relatively low level of supercoiling, and the transformation between branched and interwound structures at a higher superhelical density.

Calculations of nucleic acid conformations

The present computational power and sophistication of theoretical approaches to nucleic acid structural investigation are sufficient for the realization of static and dynamic models that correlate accurately with current crystallographic, NMR and solution-probing structural data, and consequently are able to provide valuable insights and predictions for a variety of nucleic acid conformational families. In molecular dynamics simulations, the year 1995 was marked by the foray of fast Ewald methods, an accomplishment resulting from several years' work in the search for an adequate treatment of the electrostatic long-range forces so primordial in nucleic acid behavior. In very large systems, and particularly in the RNA-folding field, techniques originating from artificial intelligence research, like constraint satisfaction programming or genetic algorithms, have established their utility and potential.

Simulating DNA at low resolution

The past year has witnessed the development of several new mathematical approaches to analyzing the structure of double-helical DNA and to incorporating the sequence-dependent features of the chain in computer simulations of long polymers. Of special interest in this respect are the local and global structural changes induced by the binding of various proteins to DNA, ranging from subtle bending, untwisting and sliding motions at the base-pair level to the apparent organization of supercoiled structure in chains that are thousands residues long. The computational effort has also included both new ways to incorporate the polyelectrolyte character of DNA and other environmental forces in simulations of long chains and new methods to keep track of the multitude of configurations so generated. The collective advances are pointing to ways that will soon connect the sequences of base pairs in large genomes to folded three-dimensional structures based on natural bending, twisting and translational tendencies and in response to deformations produced by the binding of different proteins.