Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

DNA sequence context conceals α-anomeric lesions

DNA sequence context has long been known to modulate detection and repair of DNA damage. Recent studies using experimental and computational approaches have sought to provide a basis for this observation. We have previously shown that an α-anomeric adenosine (αA) flanked by cytosines (5'CαAC-3') resulted in a kinked DNA duplex with an enlarged minor groove. Comparison of different flanking sequences revealed that a DNA duplex containing a 5'CαAG-3' motif exhibits unique substrate properties. However, this substrate was not distinguished by unusual thermodynamic properties. To understand the structural basis of the altered recognition, we have determined the solution structure of a DNA duplex with a 5'CαAG-3' core, using an extensive set of restraints including dipolar couplings and backbone torsion angles. The NMR structure exhibits an excellent agreement with the data (total R(X) <5.3%). The αA base is intrahelical, in a reverse Watson-Crick orientation, and forms a weak base pair with a thymine of the opposite strand. In comparison to the DNA duplex with a 5'CαAC-3' core, we observe a significant reduction of the local perturbation (backbone, stacking, tilt, roll, and twist), resulting in a straighter DNA with narrower minor groove. Overall, these features result in a less perturbed DNA helix and obscure the presence of the lesion compared to the 5'CαAC-3' sequence. The improved stacking of the 5'CαAG-3' core also affects the energetics of the DNA deformation that is required to form a catalytically competent complex. These traits provide a rationale for the modulation of the recognition by endonuclease IV.     

DNA fine structure and dynamics in crystals and in solution: the impact of BI/BII backbone conformations

Sugar-phosphate backbone conformations are an important structural element for a complete understanding of specific recognition in nucleic acid-protein interactions. They can be involved both in early stages of target discrimination and in structural adaptation upon binding. In the first part of this study, we have analyzed high-resolution structures of double-stranded B-DNA either isolated or bound to proteins, and explored the impact of both the standard BI and the unusual BII phosphate backbone conformations on neighboring sugar puckers and on selected helical parameters. Correlations are found to be similar for free and bound DNA, and in both categories, the possible facing backbone conformations (BI.BI, BI.BII, and BII.BII) define well-characterized substates in the B-DNA conformational space. Notably, BII.BII steps are characterized by specific, and sequence-independent, structural effects involving reduced standard deviations for almost all conformational parameters. In the second part of this work, we analyze four 10 ns molecular dynamics simulations in explicit solvent on the DNA targets of NF-kappaB and bovine papillomavirus E2 proteins, highlighting the multiplicity of backbone dynamical behavior. These results show sequence effects on the percentages of BI and BII conformers, the preferential state of facing backbones, the occurrence of coupled transitions. The backbone states can consequently be seen as a mechanism for transmitting information from the bases to the phosphate groups and thus for modulating the overall structural properties of the target DNA.     

Index to Subjects,Vol. 27

Abstract

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (χ) and sugar torsions (both the ring puckering and the exocyclic C5′-C4′ (ψ) torsion) are treated as independent variables and the resulting O3′…O5′ distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles- φ,ω,ω,ω and φ about the C3′-O3′-P-O5′-C5′ bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of thelocal sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Energy hyperspace for stacking interaction in AU/AU dinucleotide step: Dispersion‐corrected density functional theory study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence‐directed features, such as small roll values for the purine–pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near‐zero slide to be most favorable for the purine–pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar–phosphate backbone with C3′‐endo sugars and this demands C1′–C1′ distance of about 5.4 Å along the chains. Consideration of an energy penalty term for deviation of C1′–C1′ distance from the mean value, to the recent DFT‐D functionals, specifically ωB97X‐D appears to predict reliable energy contour for AU/AU step. Such distance‐based penalty improves energy contours for the other purine–pyrimidine sequences also. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 107–120, 2014.     

Base interaction and conformation manifestations of repetitive nucleotide sequences

To understand why different nucleotide sequences prefer different double helical conformations and to predict conformational behaviour of definite sequences the base-base interaction energy in regular helices consisting of A:U, A:T, G:C and I:C (hypoxanthine-cytosine) base pairs was calculated. Interaction energy was assumed to be a function of eight conformational parameters: H, the distance between adjacent pairs along helix axes; tau, turn angle of one pair relative to the neighbouring one; angles between base planes in a pair (TW, propeller twist and BL, buckle) and position of pairs with respect to helix axes (D and SL, displacements in the plane normal to helix axes, and TL and RL, inclinations to this plane, tilt and roll, respectively). For H and tau characteristic of A- and B-families of nucleic acid conformations (2.5 A less than H less than or equal to 3.5 A, 30 degrees less than or equal to tau less than or equal to 45 degrees) the ranges of conformational parameters corresponding to energy values close to minimal ones (valleys) and correlations between conformational parameters were revealed. Valleys for different sequences largely coincide but have distinctive characteristics for each sequence. Reasons for base pair planarity distortion in double-stranded helices were considered. The calculations permit to account for A-phility of G:C sequences and B-phility for A:T sequences. The valley for I:C sequence branches. This corresponds to A:T-like behaviour in some cases and G:C-like in the others.     

Determining local conformational variations in DNA. Nuclear magnetic resonance structures of the DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC) generated using back-calculation of the nuclear Overhauser effect spectra, a distance geometry algorithm and constrained molecular dynamics

Two-dimensional nuclear magnetic resonance (n.m.r.) spectroscopy and a variety of computational techniques have been used to generate three-dimensional structures of the two DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC). The central six base-pairs in these two decamers contain all ten dinucleotide pairs in DNA and thus, represent a model system for investigating how the local structure of DNA varies with base sequence. Resonance assignments were made for the non-exchangeable base protons and most of the C-1'-C-4' sugar protons in both decamers. Three-dimensional structures were generated using a distance geometry algorithm and these initial structures were refined by optimizing the fit of back-calculated spectra against the experimental two-dimensional nuclear Overhauser effect (NOE) spectra. This back-calculation procedure consists of calculating NOE cross relaxation rates for a given structure by solution of the Bloch equations, and directly accounts for spin diffusion effects. Use of this refinement procedure eliminates some assumptions that have been invoked when generating structures of DNA oligomers from n.m.r. data. Constrained energy minimization and constrained quenched molecular dynamics calculation were also performed on both decamers to help generate energetically favorable structures consistent with the experimental data. Analysis of the local conformational parameters of helical twist, helical rise, propeller twist, displacement and the alpha, beta, gamma, epison and zeta backbone torsion angles in these structures shows that these parameters span a large range of values relative to the X-ray data of nucleic acids. However, the glycosidic and pseudorotation angles are quite well defined in these structures. The implications that these results have for determination of local structural variations of DNA in solution, such as those predicted by Callidine's rules, are discussed. Our results differ significantly from some previous studies on determining local conformations of nucleic acids and comparisons with these studies are made.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

The crystal structure of an ‘All Locked’ nucleic acid duplex

‘Locked nucleic acids’ (LNAs) are known to introduce enhanced bio- and thermostability into natural nucleic acids rendering them powerful tools for diagnostic and therapeutic applications. We present the 1.9 Å X-ray structure of an ‘all LNA’ duplex containing exclusively modified β-d-2′-O-4′C-methylene ribofuranose nucleotides. The helix illustrates a new type of nucleic acid geometry that contributes to the understanding of the enhanced thermostability of LNA duplexes. A notable decrease of several local and overall helical parameters like twist, roll and propeller twist influence the structure of the LNA helix and result in a widening of the major groove, a decrease in helical winding and an enlarged helical pitch. A detailed structural comparison to the previously solved RNA crystal structure with the corresponding base pair sequence underlines the differences in conformation. The surrounding water network of the RNA and the LNA helix shows a similar hydration pattern.     

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.     

Nucleic Acid Structure Analysis: A Users Guide to a Collection of New Analysis Programs

Abstract

Common nomenclature describing the geometry of nucleic acid structures was established at a 1988 EMBO Workshop on DNA Curvature and Bending (Diekmann, S. (1988) J. Mol. Biol. 208, 787–791; Diekmann, S. (1989) The EMBO Journal 8, 1–4; Sarma, RH. (1988) J. Biomol. Structure & Dynamics 6, 391–395; Dickerson, R.E. (1989) J. Biomol. Structured Dynamics 6, 627–634; Dickerson, RE. et al. (1989)Nuc. Acids Res. 17, 1979–1803). We have subsequently developed and incorporated sophisticated mathematics in a computer program to calculate the parameters described by the guidelines. The program calculates all the local parameters relating complementary bases and neighboring base and base pairs in both Cartesian and helical coordinate frames. In addition, the main mathematical property requested by the EMBO guidelines—that the magnitude of the parameters be independent of strand or direction of measurement—is accomplished without the use of a midway coordinate frame for the rotational parameters. The mathematics preserve the physical intuition used in defining the parameters; in particular, the rotational parameters are true rotations based on a simple physical model (rotation at constant angular velocity for a unit amount of time), not Euler angles or angles between vectors and planes as is the case with other approaches. As a result, the mathematical equations are symmetric with the property that a 5° tilt is the same as a 5° roll or a 5° twist, except that the rotations take place about different axes. In other approaches, a 5° tilt can mean a different amount of net rotation from a 5° roll or a 5° twist. In addition, a great deal of flexibility is built into the program so that the user has control over the analysis, including the input format, the coordinate frame used for the base pairing relationship, the point about which the rotations are performed, and which geometric relationships are analyzed. While there is a great deal of flexibility, the program is easy to use. Interactive queries and user accessible files make the options in the program very convenient to tailor to individual needs. In addition, there is also a program that calculates bond lengths, valence angles, and torsion angles along the nucleic acid backbone, and within the sugar and base rings. Another program ‘learns’ the identities of the bond lengths, valence angles, and torsion angles that the user would like to determine. This last program is especially useful for calculating the hydrogen bonds between atoms in complementary strands as well as the unusual hydrogen bonds found in recently determined nucleic acid NMR structures or within protein/nucleic acid complexes.     

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.     

Influence of Drug Binding on DNA Flexibility: A Normal Mode Analysis

Abstract

DNA-drug complexes are important because of their pharmacological interest but, in addition, they provide a useful model to study the essential aspects of DNA recognition processes. In order to investigate the influence of ligand binding on the dynamic properties of DNA we have carried out normal mode analysis for complexes with drugs of two types: a typical intercalator, 9-aminoacridine, and a typical groove binder, netropsin. Normal modes are analysed in terms of helicoidal parameter variations with special attention being paid to global deformations of the double helix. The results show that the influence of these two drugs is very different. Intercalation of 9-aminoacridine leads to an increase in the flexibility of the intercalated dinucleotide step, with notably larger vibrational amplitudes for both roll and twist parameters compared to free DNA. In contrast, the groove binding of netropsin induces a stiffening of the DNA segment which is in contact with the drug reflected by decreased vibrational amplitudes for backbone angles and inter base pair helicoidal parameters and an increase in vibrations for adjacent base pairs in terms of buckle and propeller twist.     

Alpha/gamma transitions in the B-DNA backbone

In the crystal structures of protein complexes with B-DNA, α and γ DNA backbone torsion angles often exhibit non-canonical values. It is not known if these alternative backbone conformations are easily accessible in solution and can contribute to the specific recognition of DNA by proteins. We have analysed the coupled transition of the α and γ torsion angles within the central GpC step of a B-DNA dodecamer by computer simulations. Five stable or metastable non-canonical α/γ sub-states are found. The most favourable pathway from the canonical α/γ structure to any unusual form involves a counter-rotation of α and γ, via the trans conformation. However, the corresponding free energy indicates that spontaneous flipping of the torsions is improbable in free B-DNA. This is supported by an analysis of the available high resolution crystallographic structures showing that unusual α/γ states are only encountered in B-DNA complexed to proteins. An analysis of the structural consequences of α/γ transitions shows that the non-canonical backbone geometry influences essentially the roll and twist values and reduces the equilibrium dispersion of structural parameters. Our results support the hypothesis that unusual α/γ backbones arise during protein–DNA complexation, assisting the fine structural adjustments between the two partners and playing a role in the overall complexation free energy.     

Solution structure of [d(GCGTATACGC)]2.

The solution structure of the alternating pyrimidine-purine DNA duplex [d(GCGTATACGC)]2 has been determined using two-dimensional nuclear magnetic resonance techniques and distance geometry methods. Backbone distance constraints derived from experimental nuclear Overhauser enhancement and J-coupling torsion angle constraints were required to adequately define the conformation of the inter-residue backbone linkages and to avoid underwinding of the duplex. The distance geometry structures were further refined by back-calculation of the two-dimensional nuclear Overhauser enhancement spectra to correct spin-diffusion distance errors. Fifteen final structures for [d(GCGTATACGC)]2 were generated from the refined experimental distance bounds. These structures all exhibit fully wound B-form geometry with small penalty values (< 1.5 A) against the distance bounds and small pair-wise root-mean-square deviation values (typically 0.6 A to 1.5 A). The final structures exhibit positive base-pair inclination with respect to the helix axis, a marked alternation in rise and twist, and are shorter and wider than classical fiber B-form DNA. The purines were found to adopt a sugar pucker close to the C-2'-endo conformation while pyrimidine sugars exhibited significantly lower pseudorotation phase angles in the C-1'-exo to C-2'-endo range. The minor groove cross-strand steric clashes at pyrimidine-purine steps that would exist in pure B-DNA are attenuated by an increased rise at these steps (and an increased roll angle at TpA steps). Concomitantly the backbone torsion angles of the pyrimidine moieties have larger gamma values, larger epsilon values, and smaller zeta values than the purines. The structures generated by distance geometry methods were also compared with those obtained from restrained molecular dynamics with empirical force-field potentials. The results indicate that the nuclear magnetic resonance/distance geometry approach alone is capable of elucidating most of the salient structural features of double-stranded helical nucleic acids in solution without resorting to empirical energy potentials and without using any structural assumptions from crystallographic data.     

Structural features of B-DNA dodecamer crystal structures: influence of crystal packing versus base sequence

We have analyzed the set of inter and intra base pair parameters for each dinucleotide step in single crystal structures of dodecamers, solved at high and medium resolution and all crystallized in P2(1)2(1)2(1) space group. The objective was to identify whether all the structures which have either the Drew-Dickerson (DD) sequence d[CGCGAATTCGCG] with some base modification or related sequence (non-DD), would display the same sequence dependent structural variability about its palindromic sequence, despite the molecule being bent at one end because of similar crystal lattice packing effect. Most of the local doublet parameters for base pairs steps G2-C3 and G10-C11 positions, symmetrically situated about the lateral two-fold, were significantly correlated between themselves. In non-DD sequences, significant correlations between these positional parameters were absent. The different range of local step parameter values at each sequence position contributed to the gross feature of smooth helix axis bending in all structures. The base pair parameters in some of the positions, for medium resolution DD sequence, were quite unlike the high-resolution set and encompassed a higher range of values. Twist and slide are the two main parameters that show wider conformational range for the middle region of non-DD sequence structures in comparison to DD sequence structures. On the contrary, the minor and major groove features bear good resemblance between DD and non-DD sequence crystal structure datasets. The sugar-phosphate backbone torsion angles are similar in all structures, in sharp contrast to base pair parameter variation for high and low resolution DD and non-DD sequence structures, consisting of unusual (epsilon = g-, xi = t) BII conformation at the 10th position of the dodecamer sequence. Thus examining DD and non-DD sequence structures packed in the same crystal lattice arrangement, we infer that inter and intra base pair parameters are as symmetrically equivalent in its value as the symmetry related step for the palindromic DD sequence about lateral two-fold axis. This feature would lead us to agree with the conclusion that DNA conformation is not substantially affected by end-to-end or lateral inter-molecular interaction due to crystal lattice packing effect. Non-DD sequence structures acquire step parameter values which reflect the altered sequence at each of the dodecamer sequence position in the orthorhombic lattice while showing similar gross features of DD sequence structures.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

DNA bending: the prevalence of kinkiness and the virtues of normality.

DNA bending in 86 complexes with sequence-specific proteins has been examined using normal vector plots, matrices of normal vector angles between all base pairs in the helix, and one-digit roll/slide/twist tables. FREEHELIX, a new program especially designed to analyze severely bent and kinked duplexes, generates the foregoing quantities plus local roll, tilt, twist, slide, shift and rise parameters that are completely free of any assumptions about an overall helix axis. In nearly every case, bending results from positive roll at pyrimidine-purine base pair steps: C-A (= T-G), T-A, or less frequently C-G, in a direction that compresses the major groove. Normal vector plots reveal three well-defined types of bending among the 86 examples: (i) localized kinks produced by positive roll at one or two discrete base pairs steps, (ii) three-dimensional writhe resulting from positive roll at a series of adjacent base pairs steps, or (iii) continuous curvature produced by alternations of positive and negative roll every 5 bp, with side-to-side zig-zag roll at intermediate position. In no case is tilt a significant component of the bending process. In sequences with two localized kinks, such as CAP and IHF, the dihedral angle formed by the three helix segments is a linear function of the number of base pair steps between kinks: dihedral angle = 36 degrees x kink separation. Twenty-eight of the 86 examples can be described as major bends, and significant elements in the recognition of a given base sequence by protein. But even the minor bends play a role in fine-tuning protein/DNA interactions. Sequence-dependent helix deformability is an important component of protein/DNA recognition, alongside the more generally recognized patterns of hydrogen bonding. The combination of FREEHELIX, normal vector plots, full vector angle matrices, and one-digit roll/slide/twist tables affords a rapid and convenient method for assessing bending in DNA.