The force-displacement response of a single duplex DNA molecule was measured. The force saturates at a plateau around 70 piconewtons, which ends when the DNA has been stretched about 1.7 times its contour length. This behavior reveals a highly cooperative transition to a state here termed S-DNA. Addition of an intercalator suppresses this transition. Molecular modeling of the process also yields a force plateau and suggests a structure for the extended form. These results may shed light on biological processes involving DNA extension and open the route for mechanical studies on individual molecules in a previously unexplored range.

An algorithm is presented which solves the problem of obtaining a rigorous helicoidal description of an irregular nucleic acid segment. Central to this approach is the definition of a function describing simultaneously the curvature of the nucleic acid segment in question and the corresponding stepwise variation of helicoidal parameters along the segment. Minimisation of this function leads to an optimal distribution of the conformational irregularity of the segment between these two components. Further, it is shown that this approach can be applied equally easily to single or double stranded nucleic acids. The results of this analysis yield both the absolute helicoidal parameters of individual bases/base pairs and the relative helicoidal parameters between successive bases/base pairs as well as the overall locus of the helical axis. The possibilities of this mathematical approach are demonstrated with the help of a computer program termed “Curves” which is applied to the study of a number of different nucleic acid structures.

We describe Curves+, a new nucleic acid conformational analysis program which is applicable to a wide range of nucleic acid structures, including those with up to four strands and with either canonical or modified bases and backbones. The program is algorithmically simpler and computationally much faster than the earlier Curves approach, although it still provides both helical and backbone parameters, including a curvilinear axis and parameters relating the position of the bases to this axis. It additionally provides a full analysis of groove widths and depths. Curves+ can also be used to analyse molecular dynamics trajectories. With the help of the accompanying program Canal, it is possible to produce a variety of graphical output including parameter variations along a given structure and time series or histograms of parameter variations during dynamics.

The algorithm “Curves”, that we have recently presented in this journal (J. Biolmol. Str. Dynam. 6, 63–91 (1988)), is updated to take into account the conventions developed at the Cambridge meeting on DNA curvature (September 1988) and extended to the calculation of local parameters. In addition, the principles which govern the choices made in establishing the Curves algorithm are compared with the approaches adopted by other authors.

We investigate structural transitions within a single stretched and supercoiled DNA molecule. With negative supercoiling, for a stretching force >0.3 pN, we observe the coexistence of B-DNA and denatured DNA from σ ≈ −0.015 down to σ = −1. Surprisingly, for positively supercoiled DNA ( σ > +0.037) stretched by 3 pN, we observe a similar coexistence of B-DNA and a new, highly twisted structure. Experimental data and molecular modeling suggest that this structure has ≈2.62 bases per turn and an extension 75% larger than B-DNA. This structure has tightly interwound phosphate backbones and exposed bases in common with Pauling’s early DNA structure [Pauling, L. & Corey, R. B. (1953), Proc. Natl. Acad. Sci. USA 39, 84–97] and an unusual structure proposed for the Pf1 bacteriophage [Liu, D. J. & Day, L. A. (1994) Science 265, 671–674].

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.