Spatial translational motions of base pairs in DNA molecules: application of the extended matrix generator method

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

A-form conformational motifs in ligand-bound DNA structures

Recognition and biochemical processing of DNA requires that proteins and other ligands are able to distinguish their DNA binding sites from other parts of the molecule. In addition to the direct recognition elements embedded in the linear sequence of bases (i.e. hydrogen bonding sites), these molecular agents seemingly sense and/or induce an "indirect" conformational response in the DNA base-pairs that facilitates close intermolecular fitting. As part of an effort to decipher this sequence-dependent structural code, we have analyzed the extent of B-->A conformational conversion at individual base-pair steps in protein and drug-bound DNA crystal complexes. We take advantage of a novel structural parameter, the position of the phosphorus atom in the dimer reference frame, as well as other documented measures of local helical structure, e.g. torsion angles, base-pair step parameters. Our analysis pinpoints ligand-induced conformational changes that are difficult to detect from the global perspective used in other studies of DNA structure. The collective data provide new structural details on the conformational pathway connecting A and B-form DNA and illustrate how both proteins and drugs take advantage of the intrinsic conformational mechanics of the double helix. Significantly, the base-pair steps which exhibit pure A-DNA conformations in the crystal complexes follow the scale of A-forming tendencies exhibited by synthetic oligonucleotides in solution and the known polymorphism of synthetic DNA fibers. Moreover, most crystallographic examples of complete B-to-A deformations occur in complexes of DNA with enzymes that perform cutting or sealing operations at the (O3'-P) phosphodiester linkage. The B-->A transformation selectively exposes sugar-phosphate atoms, such as the 3'-oxygen atom, ordinarily buried within the chain backbone for enzymatic attack. The forced remodeling of DNA to the A-form also provides a mechanism for smoothly bending the double helix, for controlling the widths of the major and minor grooves, and for accessing the minor groove edges of individual base-pairs.

An A-type double helix of DNA having B-type puckering of the deoxyribose rings

DNA usually adopts structure B in aqueous solution, while structure A is preferred in mixtures of trifluoroethanol (TFE) with water. However, the octamer d(CCCCGGGG) and other d(C(n)G(n)) fragments of DNA provide CD spectra that suggest that the base-pairs are stacked in an A-like fashion even in aqueous solution. Yet, d(CCCCGGGG) undergoes a cooperative TFE-induced transition into structure A, indicating that an important part of the aqueous duplex retains structure B. NMR spectroscopy shows that puckering of the deoxyribose rings is of the B-type. Hence, combination of the information provided by CD spectroscopy and NMR spectroscopy suggests an unprecedented double helix of DNA in which A-like base stacking is combined with B-type puckering of the deoxyribose rings. In order to determine whether this combination is possible, we used molecular dynamics to simulate the duplex of d(CCCCGGGG). Remarkably, the simulations, completely unrestrained by the experimental data, provided a very stable double helix of DNA, exhibiting just the intermediate B/A features described above. The double helix contained well-stacked guanine bases but almost unstacked cytosine bases. This generated a hole in the double helix center, which is a property characteristic for A-DNA, but absent from B-DNA. The minor groove was narrow at the double helix ends but wide at the central CG step where the Watson-Crick base-pairs were buckled in opposite directions. The base-pairs stacked tightly at the ends but stacking was loose in the duplex center. The present double helix, in which A-like base stacking is combined with B-type sugar puckering, is relevant to replication and transcription because both of these phenomena involve a local B-to-A transition.

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.

Influence of dynamic fluctuations on DNA curvature

The effect of dynamic fluctuations on physical manifestations of DNA curvature such as electrophoretic retardation, circularization of DNA tracts and nucleosomes positioning is examined. It is shown that in all cases the main features of the processes can be satisfactorily explained by a static curvature model, which appears to be a good representation of time and ensemble averaged superstructures of DNA chains. The dynamic fluctuations around the average curvature appear to influence only the kinetics of these processes. In the case of polyacrylamide gel electrophoretic retardation it is demonstrated that the approximation of the static model holds on the assumption that dynamic fluctuations are independent from intrinsic curvature. The actual validity of the static model we proposed several years ago is satisfactorily demonstrated by the explanation and prediction of different experiments, such as cyclic permutation gel electrophoresis, differential DNAase I cleavage of cyclic versus linear DNA tracts and nucleosome positioning.

Flexing and folding double helical DNA

DNA base sequence, once thought to be interesting only as a carrier of the genetic blueprint, is now recognized as playing a structural role in modulating the biological activity of genes. Primary sequences of nucleic acid bases describe real three-dimensional structures with properties reflecting those structures. Moreover, the structures are base sequence dependent with individual residues adopting characteristic spatial forms. As a consequence, the double helix can fold into tertiary arrangements, although the deformation is much more gradual and spread over a larger molecular scale than in proteins. As part of an effort to understand how local structural irregularities are translated at the macromolecular level in DNA and recognized by proteins, a series of calculations probing the structure and properties of the double helix have been performed. By combining several computational techniques, complementary information as well as a series of built-in checks and balances for assessing the significance of the findings are obtained. The known sequence dependent bending, twisting, and translation of simple dimeric fragments have been incorporated into computer models of long open DNAs of varying length and chemical composition as well as in closed double helical circles and loops. The extent to which the double helix can be forced to bend and twist is monitored with newly parameterized base sequence dependent elastic energy potentials based on the observed configurations of adjacent base pairs in the B-DNA crystallographic literature.

Static contributions to the persistence length of DNA and dynamic contributions to DNA curvature

Long molecules of DNA have the statistical properties of a worm-like coil. Deviations from linearity occur both because of small dynamic bends induced by thermal motion and from a random distribution of static bends. The latter originate in the different conformations of each of the possible base pair sequences. In this paper a statistical theory of the persistence length of DNA is developed which includes both static and dynamic effects for each base pair sequence, as well as the sequence-dependent correlations of bending angles. The result applies to a generic DNA, i.e., the average over an ensemble of all possible sequences. The theory is also applied to the generation of the average properties of curved DNAs by an analytic method that includes dynamic averaging as well as correlated bends. These results provide information which supplements that obtained by others using Monte Carlo methods. The additivity relation 1/P = 1/P(S) + 1/P(d) proposed by Trifonov et al., where P is the persistence length and P(S) and P(d) are the persistence lengths arising from purely static and dynamic effects, respectively, has been verified to be accurate to better than 0.5%. This is true for both a simplified model and one that includes a complete set of static bends at all base pair sequences.