DNA polymorphism and local variation in base-pair orientation: a theoretical rationale.

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B----A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.     

Groove width and depth of B-DNA structures depend on local variation in slide.

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA. A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.     

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

The solution conformation of a carbocyclic analog of the Dickerson-Drew dodecamer: comparison with its own X-ray structure and that of the NMR structure of the native counterpart

The NMR conformation of a carbocyclic analog of the Dickerson-Drew dodecamer [d(CGC-GAAT*T*CGCG)]2 containing 6'-alpha-Me carbocyclic thymidines (T*) has been determined and compared with that of its X-ray structure. The solution structure of the 6'-alpha-Me carbocyclic thymidine modified duplex has also been compared with the solution structure of the corresponding unmodified Dickerson-Drew duplex solved by us under the same experimental conditions. The NMR structures have been based on 24 experimental distance and torsion constraints per residue for [d(CGCGAAT*T*CGCG)]2 (1) and on 21 constraints per residue for the natural counterpart. In general, both final NMR structures are more close to the B-type DNA. The cyclopentane moieties of the carbocyclic thymidine residues adopt C1'-exo B-DNA type puckers (the phase angles P = 136-139 degrees and the puckering amplitudes psi = 36-37 degrees) that are close to their previously published crystal C1'-exo or C2'-endo puckers. The main differences between the two NMR structures are for beta(T*8) and epsilon, xi(T*7) backbone torsions (27-50 degrees ), for basepair twist for the 7-8 and 8-9 basepair steps (5-6 degrees), tilt for the 8-9 step (7 degrees), roll for the 7-8 step (7 degrees), shift for the 7-8 step (0.9A) and slide for the 9-10 step (0.6A). The relatively small deviations of helical structure parameters lead to structural isomorphism of these duplexes in aqueous solutions (atomic RMSD = 1.0A). The difference of the minor groove widths (less than 0.7A) in the core part of the modified duplex in comparison with the native one is much smaller than the difference between the X-ray structures of these duplexes. A detailed comparison of NMR and X-ray structure parameters showed significant monotonic differences (0.9-2.5A) for all basepair slides in both duplexes. Deviations between NMR and X-ray structure parameters for the modified duplex were also found for basepair tilt of the 4-5 step (13 degrees), rolls for the 8-9 and 10-11 steps (16 and 19 degrees), twist of the 3-4 step (8 degrees) and shift of the 9-10 step (0.9A).     

Groove Width and Depth of B-DNA Structures Depend on Local Variation in Slide

Abstract

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

Sequence-dependent DNA structure: a database of octamer structural parameters

We have constructed the potential energy surfaces for all unique tetramers, hexamers and octamers in double helical DNA, as a function of the two principal degrees of freedom, slide and shift at the central step. From these potential energy maps, we have calculated a database of structural and flexibility properties for each of these sequences. These properties include: the values of each of the six step parameters (twist roll, tilt, rise, slide and shift), for each step of the sequence; flexibility measures for both decrease and increase in each property value from the minimum energy conformation for the central step; and the deviation from the path of a hypothetical straight octamer. In an analysis of structural change as a function of sequence length, we observe that almost all DNA tends to B-DNA and becomes less flexible. A more detailed analysis of octamer properties has allowed us to determine the structural preferences of particular sequence elements. GGC and GCC sequences tend to confer bistability, low stability and a predisposition to A-form DNA, whereas AA steps strongly prefer B-DNA and inhibit A-structures. There is no correlation between flexibility and intrinsic curvature, but bent DNA is less stable than straight. The most difficult deformation is undertwisting. The TA step stands out as the most flexible sequence element with respect to decreasing twist and increasing roll. However, as with the structural properties, this behavior is highly context-dependent and some TA steps are very straight.     

Overview of the structure of all-AT oligonucleotides: organization in helices and packing interactions

We present the crystalline organization of 33 all-AT deoxyoligonucleotide duplexes, studied by x-ray diffraction. Most of them have very similar structures, with Watson-Crick basepairs and a standard average twist close to 36 degrees. The molecules are organized as parallel columns of stacked duplexes in a helical arrangement. Such organization of duplexes is very regular and repetitive: all sequences show the same pattern. It is mainly determined by the stacking of the terminal basepairs, so that the twist in the virtual TA base step between neighbor duplexes is always negative, approximately -22 degrees. The distance between the axes of parallel columns is practically identical in all cases, approximately 26 A. Interestingly, it coincides with that found in DNA viruses and fibers in their hexagonal phase. It appears to be a characteristic distance for ordered parallel DNA molecules. This feature is due to the absence of short range intermolecular forces, which are usually due to the presence of CG basepairs at the end of the oligonucleotide sequence. The duplexes apparently interact only through their diffuse ionic atmospheres. The results obtained can thus be considered as intermediate between liquid crystals, fibers, and standard crystal structures. They provide new information on medium range DNA-DNA interactions.     

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.     

Molecular mechanics studies on poly(purine).poly(pyrimidine) sequences in DNA: polymorphism and local variability

Energy minimization has been carried out on three poly(purine).poly(pyrimidine) sequences--d(G)10.d(C)10, d(A)10.d(T)10, and d(AG)5.d(CT)5--using the molecular mechanics program AMBER (Assisted Model Building and Energy Refinement). In order to extensively scan the conformational space available, five different helical models were studied, three of them being right-handed helices while the other two were left helical. For all three sequences the right-handed A- and B-type helices are energetically slightly preferred over the left helices, but the energy difference between the various right-handed helices is only marginal. A detailed analysis has been carried out to characterize the local structural variability in the refined structures, both in terms of torsion angles as well as other parameters such as base-pair tilt, wedge roll, and wedge tilt, etc. All three sequences exhibit similar structural features for a particular form, but both the forms A and B show significant deviations from fiber models. In particular, the A-form structures have higher unit rise (2.7 A), and lower unit twist (31 degrees) and base-pair tilt (12 degrees), compared to the fiber model, which has corresponding values of 2.56 A, 32.7 degrees, and 20 degrees, respectively. All these changes indicate that the refined models are closer to the A-form structure observed in crystals of oligonucleotides. In the refined B-for models, the helical parameters are close to the fiber B-form, although the torsion angles show considerable variations. None of the three sequences examined, including the d(A)n.d(T)n sequence, show any pronounced curvature for the B-form structure.     

The solution structure of a locked nucleic acid (LNA) hybridized to DNA

LNA (Locked Nucleic Acids) is a novel oligonucleotide analogue containing a conformationally restricted nucleotide with a 2'-O, 4'-C-methylene bridge that induces unprecedented thermal affinities when mixed with complementary single stranded DNA and RNA. We have used two-dimensional 1H NMR spectroscopy obtained at 750 and 500 MHz to determine a high resolution solution structure of an LNA oligonucleotide hybridized to the complementary DNA strand. The determination of the structure was based on a complete relaxation matrix analysis of the NOESY cross peaks followed by restrained molecular dynamics calculations. Forty final structures were generated for the duplex from A-type and B-type dsDNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the forty structures of the complex was 0.32A. The structures were analysed by use of calculated helix parameters. This showed that the values for rise and buckle in the LNA duplex is markedly different from canonical B-DNA at the modification site. A value of twist similar to A-DNA is also observed at the modification site. The overall length of the helix which is 27.3 A. The average twist over the sequence are 35.9 degrees +/- 0.3 degrees. Consequently, the modification does not cause the helix to unwind. The bis-intercalation of the thiazole orange dye TOTO to the LNA duplex was also investigated by 1H NMR spectroscopy to sense the structural change from the unmodified oligonucleotide. We observed that the bis-intercalation of TOTO is much less favourable in the 5'-CT(L)AG-3' site than in the unmodified 5'-CTAG-3' site. This was related to the change in the base stacking of the LNA duplex compared to the unmodified duplex.     

Sensitivity of NMR internucleotide distances to B-DNA conformation: underlying mechanics.

Nuclear magnetic resonance (NMR) spectroscopy, combining correlated spectroscopy (COSY) coupling constant measurements with nuclear Overhauser effect spectroscopy (NOESY) interatomic distances, should make it possible to determine an averaged solution structure for DNA oligomers. However, even if such data could be obtained with high accuracy, it is not clear which structural parameters of DNA would be determined. Here, the relationships between measurable internucleotide distances and helical parameters are systematically studied through molecular modelling. Investigations are carried out using four representative sequences, (ACGT)n, (TCGA)n, (AGCT)n and (TGCA)n, composed of repeated tetranucleotides belonging to oligomers previously studied by NMR. Correlations between interatomic distances become evident and strong connections between distances and inter-base helical parameters are observed. Results imply that twist, roll, shift and slide values can be accurately determined from NMR data. Sequence independent mechanical coupling which link backbone and sugar conformations to helical twist are also described.     

Local sequential minimization of double stranded B-DNA using Monte Carlo annealing

A software algorithm has been developed to investigate the folding process in B-DNA structures in vacuum under a simple and accurate force field. This algorithm models linear double stranded B-DNA sequences based on a local, sequential minimization procedure. The original B-DNA structures were modeled using initial nucleotide structures taken from the Brookhaven database. The models contain information at the atomic level allowing one to investigate as accurately as possible the structure and characteristics of the resulting DNA structures. A variety of DNA sequences and sizes were investigated containing coding and non-coding, random and real, homogeneous or heterogeneous sequences in the range of 2 to 40 base pairs. The force field contains terms such as angle bend, Lennard-Jones, electrostatic interactions and hydrogen bonding which are set up using the Dreiding II force field and defined to account for the helical parameters such as twist, tilt and rise. A close comparison was made between this local minimization algorithm and a global one (previously published) in order to find out advantages and disadvantages of the different methods. From the comparison, this algorithm gives better and faster results than the previous method, allowing one to minimize larger DNA segments. DNA segments with a length of 40 bases need approximately 4 h, while 2.5 weeks are needed with the previous method. After each minimization the angles between phosphate–oxygen-carbon A₁, the oxygen–phosphate–oxygen A₂ and the average helical twists were calculated. From the generated fragments it was found that the bond angles are A₁=150°±2°and A₂=130°±10°, while the helical twist is 36.6°±2° in the A strand and A₁=150°±6° and A₂=130±6° with helical twist 39.6°±2° in the B strand for the DNA segment with the same sequence as the Dickerson dodecamer.Figure The final minimized DNA segment of the Dickerson dodecamer sequence represented by ball drawings and viewed (left) perpendicular and (right) down the helical axis     

Effect of neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.     

Effect of sequence on the conformation of DNA holliday junctions

Structures of the DNA sequences d(CCGGCGCCGG) and d(CCAGTACbr(5)UGG) are presented here as four-way Holliday junctions in their compact stacked-X forms, with antiparallel alignment of the DNA strands. Thus, the ACC-trinucleotide motif, previously identified as important for stabilizing the junction, is now extended to PuCPy, where Pu is either an adenine or guanine, and Py is either a cytosine, 5-methylcytosine, or 5-bromouracil but not thymine nucleotide. We see that both sequence and base substituents affect the geometry of the junction in terms of the interduplex angle as well as the previously defined conformational variables, J(roll) (the rotation of the stacked duplexes about their respective helical axis) and J(slide) (the translational displacement of the stacked duplexes along their respective helical axis). The structures of the GCC and parent ACC containing junctions fall into a distinct conformational class that is relatively undistorted in terms of J(slide) and J(roll), with interduplex angles of 40-43 degrees. The substituted ACbr(5)U structure, however, is more akin to that of the distorted methylated ACm(5)C containing junction, with J(slide) (>or=2.3 A) and a similar J(roll) (164 degrees) opening the major groove-side of the junction, but shows a reduced interduplex angle. In contrast, the analogous d(CCAGTACTGG) sequence has to date been crystallized only as resolved B-DNA duplexes. This suggests that there is an electronic effect of substituents at the pyrimidine Py position on the stability of four-stranded junctions. The single-crystal structures presented here, therefore, show how sequence affects the detailed geometry, and subsequently, the associated stability and conformational dynamics of the Holliday junction.     

B-DNA under stress: over- and untwisting of DNA during molecular dynamics simulations

The twist flexibility of DNA is central to its many biological functions. Explicit solvent molecular dynamics simulations in combination with an umbrella sampling restraining potential have been employed to study induced twist deformations in DNA. Simulations allowed us to extract free energy profiles for twist deformations and were performed on six DNA dodecamer duplexes to cover all 10 possible DNA basepair steps. The shape of the free energy curves was similar for all duplexes. The calculated twist deformability was in good agreement with experiment and showed only modest variation for the complete duplexes. However, the response of the various basepair steps on twist stress was highly nonuniform. In particular, pyrimidine/purine steps were much more flexible than purine/purine steps followed by purine/pyrimidine steps. It was also possible to extract correlations of twist changes and other helical as well as global parameters of the DNA molecules. Twist deformations were found to significantly alter the local as well as global shape of the DNA modulating the accessibility for proteins and other ligands. Severe untwisting of DNA below an average of 25 degrees per basepair step resulted in the onset of a global structural transition with a significantly smaller twist at one end of the DNA compared to the other.     

Stacking geometry for non‐canonical G:U wobble base pair containing dinucleotide sequences in RNA: dispersion‐corrected DFT‐D study

Emergence of thousands of crystal structures of noncoding RNA molecules indicates its structural and functional diversity. RNA function is based upon a large variety of structural elements which are specifically assembled in the folded molecules. Along with the canonical Watson‐Crick base pairs, different orientations of the bases to form hydrogen‐bonded non‐canonical base pairs have also been observed in the available RNA structures. Frequencies of occurrences of different non‐canonical base pairs in RNA indicate their important role to maintain overall structure and functions of RNA. There are several reports on geometry and energetic stabilities of these non‐canonical base pairs. However, their stacking geometry and stacking stability with the neighboring base pairs are not well studied. Among the different non‐canonical base pairs, the G:U wobble base pair (G:U W:WC) is most frequently observed in the RNA double helices. Using quantum chemical method and available experimental data set we have studied the stacking geometry of G:U W:WC base pair containing dinucleotide sequences in roll‐slide parameters hyperspace for different values of twist. This study indicates that the G:U W:WC base pair can stack well with the canonical base pairs giving rise to large interaction energy. The overall preferred stacking geometry in terms of roll, twist and slide for the eleven possible dinucleotide sequences is seen to be quite dependent on their sequences. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 328–338, 2015.