Supercoiled DNA energetics and dynamics by computer simulation.

A new formulation is presented for investigating supercoiled DNA configurations by deterministic techniques. Thus far, the computational difficulties involved in applying deterministic methods to supercoiled DNA studies have generally limited computer simulations to stochastic approaches. While stochastic methods, such as simulated annealing and Metropolis-Monte Carlo sampling, are successful at generating a large number of configurations and estimating thermodynamic properties of topoisomer ensembles, deterministic methods offer an accurate characterization of the minima and a systematic following of their dynamics. To make this feasible, we model circular duplex DNA compactly by a B-spline ribbon-like model in terms of a small number of control vertices. We associate an elastic deformation energy composed of bending and twisting integrals and represent intrachain contact by a 6-12 Lennard Jones potential. The latter is parameterized to yield an energy minimum at the observed DNA-helix diameter inclusive of a hydration shell. A penalty term to ensure fixed contour length is also included. First and second partial derivatives of the energy function have been derived by using various mathematical simplifications. First derivatives are essential for Newton-type minimization as well as molecular dynamics, and partial second-derivative information can significantly accelerate minimization convergence through preconditioning. Here we apply a new large-scale truncated-Newton algorithm for minimization and a Langevin/implicit-Euler scheme for molecular dynamics. Our truncated-Newton method exploits the separability of potential energy functions into terms of differing complexity. It relies on a preconditioned conjugate gradient method that is efficient for large-scale problems to solve approximately for the search direction at every step. Our dynamics algorithm is numerically stable over large time steps. It also introduces a frequency-discriminating mechanism so that vibrational modes with frequencies greater than a chosen cutoff frequency are essentially frozen by the method. With these tools, we rapidly identify corresponding circular and interwound energy minima for small DNA rings for a series of imposed linking-number differences. These structures are consistent with available electron microscopy data. The energetic exchange of stability between the circle and the figure-8, in very good agreement with analytical results, is also detailed. Molecular dynamics trajectories at 100 femtosecond time steps then reveal the rapid folding of the unstable circular state into supercoiled forms. Significant bending and twisting motions of the interwound structures are also observed. Such information may be useful for understanding transition states along the folding pathway and the role of enzymes that regulate supercoiling.(ABSTRACT TRUNCATED AT 400 WORDS)     

The effect of the superhelicity on the double helix twist angle in DNA.

An attempt to estimate the relative contributions of twisting and bending to the free energy of superhelix formation from the relaxed DNA is undertaken. The extent of teritiary ordering (number of DNA axis turns tau) and that of secondary ordering (duplex twist angle beta) have been taken as thermodynamical parameters, which characterize the state of the supercoild DNA at the fixed linking number (Lk) value. Such a thermodynamical approach implies the phenomenological parameters of rigidities of twisting and supercoiling (Gbeta, Gtau). Gtau/Gbeta ratio is estimated from the presented experimental data on the winding of the double helix upon increasing the ionic strength when twist alterations are followed by circular dichroism method. The adequacy of such interpretation of CD spectra changes are discussed. The values of Gtau and Gbeta are estimated to be of the same order of magnitude.     

DNA bending and twisting properties of integration host factor determined by DNA cyclization

The binding of many proteins to DNA is profoundly affected by DNA bending, twisting, and supercoiling. When protein binding alters DNA conformation, interaction between inherent and induced DNA conformation can affect protein binding affinity and specificity. Integration host factor (IHF), a sequence-specific, DNA-binding protein of Escherichia coli, strongly bends the DNA upon binding. To assess the influence of inherent DNA bending on IHF binding, we took advantage of the high degree of natural static curvature associated with an IHF site on a 163-bp minicircle and measured the binding affinity of IHF for its recognition site contained on this DNA in both circular and linear form. IHF showed a higher affinity for the circular form of the DNA when compared to the linear form. In addition, the presence of IHF during DNA cyclization changed the topology of cyclization products and their ability to bind IHF, consistent with IHF untwisting DNA. These results show that inherent DNA conformation anisotropy is an important determinant of IHF binding affinity and suggests a mechanism for modulation of IHF activity by local DNA conformation.     

DNA,Flexibly Flexible

Investigators have constructed dsDNA molecules with several different base modifications and have characterized their bending and twisting flexibilities using atomic force microscopy, DNA ring closure, and single-molecule force spectroscopy with optical tweezers. The three methods provide persistence length measurements that agree semiquantitatively, and they show that the persistence length is surprisingly similar for all of the modified DNAs. The circular dichroism spectra of modified DNAs differ substantially. Simple explanations based on base stacking strength, polymer charge, or groove occupancy by functional groups cannot explain the results, which will guide further high-resolution theory and experiments.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric A_nT_n and asymmetric A_2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.     

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.     

Origin of Twist-Bend Coupling in Actin Filaments

Actin filaments are semiflexible polymers that display large-scale conformational twisting and bending motions. Modulation of filament bending and twisting dynamics has been linked to regulatory actin-binding protein function, filament assembly and fragmentation, and overall cell motility. The relationship between actin filament bending and twisting dynamics has not been evaluated. The numerical and analytical experiments presented here reveal that actin filaments have a strong intrinsic twist-bend coupling that obligates the reciprocal interconversion of bending energy and twisting stress. We developed a mesoscopic model of actin filaments that captures key documented features, including the subunit dimensions, interaction energies, helicity, and geometrical constraints coming from the double-stranded structure. The filament bending and torsional rigidities predicted by the model are comparable to experimental values, demonstrating the capacity of the model to assess the mechanical properties of actin filaments, including the coupling between twisting and bending motions. The predicted actin filament twist-bend coupling is strong, with a persistence length of 0.15-0.4 μm depending on the actin-bound nucleotide. Twist-bend coupling is an emergent property that introduces local asymmetry to actin filaments and contributes to their overall elasticity. Up to 60% of the filament subunit elastic free energy originates from twist-bend coupling, with the largest contributions resulting under relatively small deformations. A comparison of filaments with different architectures indicates that twist-bend coupling in actin filaments originates from their double protofilament and helical structure.     

Ring closure probabilities for DNA fragments by Monte Carlo simulation

The rate of ligation of DNA molecules into circular forms depends on the ring closure probability, commonly called the j-factor, which is a sensitive measure of the extent to which thermal fluctuations contribute to bending and twisting of DNA molecules in solution. We present a theoretical treatment of the cyclization equilibria of DNA that employs a special Monte Carlo method for generating large ensembles of model DNA chains. Using this method, the chain length dependence of the j-factor was calculated for molecules. in the size range 250 to 2000 base-pairs. The Monte Carlo results are compared with recent analytical theory and experimental data. We show that a value of 475 A for the persistence length of DNA, close to values measured by a number of other methods, is in excellent agreement with the cyclization results. Preliminary applications of the Monte Carlo method to the problem of systematically bent DNA molecules are presented. The calculated j-factor is shown to be very sensitive to the amount of bending in these fragments. This fact suggests that ligase closure measurements of systematically bent DNA molecules should be a useful method for studying sequence-directed bending in DNA.     

DNA bending by charged peptides: electrophoretic and spectroscopic analyses

We are testing the idea that placement of fixed charges near one face of the DNA double helix can induce DNA bending by a purely electrostatic mechanism. If stretching forces between DNA phosphates are significant, fixed charges should induce DNA bending by asymmetrically modulating these forces. We have previously tested this hypothesis by adding charged residues to small bZIP DNA binding peptides and monitoring DNA bending using electrophoretic phasing assays. Our results were consistent with an electrostatic model of DNA bending in predicted directions. We now confirm these observations with fluorescence resonance energy transfer (FRET). Using a "U"-shaped DNA probe, we report that DNA bending by charged bZIP peptides is readily detected by FRET. We further show that charged bZIP peptides cause DNA bending rather than DNA twisting.     

Structure of plectonemically supercoiled DNA

Using electron microscopy and topological methods, we have deduced an average structure for negatively supercoiled circular DNA in solution. Our data suggest that DNA has a branched plectonemic (interwound) form over the range of supercoiling tested. The length of the superhelix axis is constant at 41% of the DNA length, whereas the superhelix radius decreases essentially hyperbolically as supercoiling increases. The number of supercoils is 89% of the linking deficit. Both writhe and twist change with supercoiling, but the ratio of the change in writhe to the change in twist is fixed at 2.6:1. The extent of branching of the superhelix axis is proportional to the length of the plasmid, but is insensitive to superhelix density. The relationship between DNA flexibility constants for twisting and bending calculated using our structural data is similar to that deduced from previous studies. The extended thin form of plectonemically supercoiled DNA offers little compaction for cellular packaging, but promotes interaction between cis-acting sequence elements that may be distant in primary structure. We discuss additional biological implications of our structural data.     

Toroidal elastic supercoiling of DNA

Covalently closed circular DNA can exist in different configurations known as circular, toroidal, and interwound. Changes among these forms can be made in several ways, including the insertion of dye molecules between adjacent base pairs, which tends to untwist the double-helical structure. The aim of this paper is to discuss these configurations, and the changes among them, in the context of classical elastomechanics. The concepts of twisting, linkage and writhing are explained. Simple experiments on a twisted linear-elastic rod are described, and it is shown that although the circular and interwound forms may be modeled in this way, the toroidal form does not occur, being mechanically unstable. Theoretical energy calculations by Levitt on bent and twisted DNA show that DNA exhibits a particular kind of nonlinear elasticity in which there is an unusual coupling between bending and twisting. The aim of the paper is to show qualitatively that this special kind of elasticity can stabilize the toroidal form of closed circular DNA.     

Correction to: Three-dimensional quantification of twisting in the Arabidopsis petiole

Organisms have a variety of three-dimensional (3D) structures that change over time. These changes include twisting, which is 3D deformation that cannot happen in two dimensions. Twisting is linked to important adaptive functions of organs, such as adjusting the orientation of leaves and flowers in plants to align with environmental stimuli (e.g. light, gravity). Despite its importance, the underlying mechanism for twisting remains to be determined, partly because there is no rigorous method for quantifying the twisting of plant organs. Conventional studies have relied on approximate measurements of the twisting angle in 2D, with arbitrary choices of observation angle. Here, we present the first rigorous quantification of the 3D twisting angles of Arabidopsis petioles based on light sheet microscopy. Mathematical separation of bending and twisting with strict definition of petiole cross-sections were implemented; differences in the spatial distribution of bending and twisting were detected via the quantification of angles along the petiole. Based on the measured values, we discuss that minute degrees of differential growth can result in pronounced twisting in petioles.

     

Mechanical properties of DNA-like polymers

The molecular structure of the DNA double helix has been known for 60 years, but we remain surprisingly ignorant of the balance of forces that determine its mechanical properties. The DNA double helix is among the stiffest of all biopolymers, but neither theory nor experiment has provided a coherent understanding of the relative roles of attractive base stacking forces and repulsive electrostatic forces creating this stiffness. To gain insight, we have created a family of double-helical DNA-like polymers where one of the four normal bases is replaced with various cationic, anionic or neutral analogs. We apply DNA ligase-catalyzed cyclization kinetics experiments to measure the bending and twisting flexibilities of these polymers under low salt conditions. Interestingly, we show that these modifications alter DNA bending stiffness by only 20%, but have much stronger (5-fold) effects on twist flexibility. We suggest that rather than modifying DNA stiffness through a mechanism easily interpretable as electrostatic, the more dominant effect of neutral and charged base modifications is their ability to drive transitions to helical conformations different from canonical B-form DNA.     

On the toroidal condensed state of closed circular DNA

The influence of double helix torsional elasticity on the compaction and structure of circular DNA compact form is studied theoretically in the case when the compact (globular) form has torus shape. For closed circular DNA the topological invariant, the linking number, yields a strict connection between conformation of the double helix considered as unifilar homopolymer and elastic energy of torsional twisting. The contribution of torsional elasticity to the free energy of the toruslike globule is calculated. This contribution is shown to be proportional to the square of superhelical density. Allowance of the torsional elasticity decreases the equilibrium radius of the toruslike globule formed by circular DNA. Closure of linear DNA into a ring widens the stability range of the relatively short DNA compact form and tightens it for long DNA.     

In vitro evaluation of combined graft deformation in anterior cruciate ligament reconstructions

In this paper, we present a simple method of calculating deformation of an anterior cruciate ligament graft in combined elongation, bending and total twisting. We also report our results on these deformations for three types of ligament reconstruction in cadaver knees: two isometric reconstructions, using either a biological bone-patellar tendon-bone autograft or an artificial Trevira prosthesis, and modified over-the-top reconstruction with the prosthesis. The data show that the modified OTT technique produced elongation of the graft equivalent to that of the synthetic isometric technique, but significantly less than isometric reconstruction with the biological graft (p<0.05, ANOVA). Moreover, the grafts were subjected to bending and twisting.     

Structure of a B-DNA decamer with a central T-A step: C-G-A-T-T-A-A-T-C-G.

The X-ray crystal structure analysis of the decamer C-G-A-T-T-A-A-T-C-G has been carried out to a resolution of 1.5 A. The crystals are space group P2(1)2(1)2(1), cell dimensions a = 38.60 A, b = 39.10 A, c = 33.07 A. The structure was solved by molecular replacement and refined with X-PLOR and NUCLSQ. The final R factor for a model with 404 DNA atoms, 108 water molecules and one magnesium hexahydrate cation is 15.7%. The double helix is essentially isostructural with C-G-A-T-C-G-A-T-C-G, with closely similar local helix parameters. The structure of the T-T-A-A center differs from that found in C-G-C-G-T-T-A-A-C-G-C-G in that the minor groove in our decamer is wide at the central T-A step rather than narrow, and the twist angle of the T-A step is small (31.1 degrees) rather than large. Whereas the tetrad model provides a convenient framework for discussing local DNA helix structure, it cannot be the entire story. The articulated helix model of DNA structure proposes that certain sequence regions of DNA show preferential twisting or bending properties, whereas other regions are less capable of deformation, in a manner that may be useful in sequence recognition by drugs and protein. Further crystal structure analyses should help to delineate the precise nature of sequence-dependent articulation in the DNA double helix.     

Toroidal globular state of circular DNA

The influence of torsional elasticity of the double helix on compactization and structure of circular DNA in a compact form is studied in the case when the compact (globular) particle has a torus shape. For closed circular DNA the topological invariant, linking number of two strains, yields strict connection between conformation of double helix, considered as a unifilar homopolymer, and elastic energy of torsional twisting. The contribution of torsional elasticity to free energy of the toruslike globule is calculated. This contribution is shown to be proportional to the square of superturn's density. Torsional elasticity decreases the equilibrium radius of the toruslike globule formed by circular DNA in comparison with the case of linear DNA. Closure of linear DNA into a ring widens the stability range of the relatively short DNA compact form and tightens it for long DNA.