Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure

Detailed examination of the structure of the B-DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G, obtained by single-crystal X-ray analysis (Drew et al., 1981), reveals that the local helix parameters, twist, tilt and roll, are much more strongly influenced by base sequence than by crystal packing or any other external forces. The central EcoRI restriction endonuclease recognition site, G-A-A-T-T-C, is a B helix with an average of 9.8 base-pairs per turn. It is flanked on either side by single-base-pair steps having aspects of an A-like helix character. The dodecamer structure suggests several general principles, whose validity must be tested by other B-DNA analyses. (1) When an external bending moment is applied to a B-DNA double helix, it bends smoothly, without kinks or breaks, and with relatively little effect on local helix parameters. (2) Purine-3′,5′-pyrimidine steps open their base planes towards the major groove, pyrimidine-purine steps open toward the minor groove, and homopolymer (Pur-Pur, Pyr-Pyr) steps resist rolling in either direction. This behavior is related to the preference of pyrimidines for more negative glycosyl torsion angles. (3) CpG steps have smaller helical twist angles than do GpC, as though in compensation for their smaller intrinsic base overlap. Data on A-T steps are insufficient for generalization. (4) G.C base-pairs have smaller propellor twist than A · T, and this arises mainly from interstrand base overlap rather than the presence of the third hydrogen bond. (5) DNAase I cuts preferentially at positions of high helical twist, perhaps because of increased exposure of the backbone to attack. The correlation of the digestion patterns in solution and helical twist in the crystal argues for the essential identity of the helix structure in the two environments. (6) In the two places where the sequence TpCpG occurs, the C slips from under T in order to stack more efficiently over G. At the paired bases of this CpG step, the G and C are tilted so the angle between base planes is splayed out to the outside of the helix. This TpC is the most favored cutting site for DNAase I by a factor of 4.5 (Lomonossoff et al., 1981). (7) The EcoRI restriction endonuclease and methylase both appear to prefer a cutting site of the type purine-purine-A-T-T-pyrimidine, involving two adjacent homopolymer triplets, and this may be a consequence of the relative stiffness of homopolymer base-stacking observed in the dodecamer.     

Sequence-dependent variation in the conformation of DNA

The specificity of action of the enzyme DNAase I on double-stranded DNA polymers of defined sequence has been investigated. The results obtained with the alternating copolymers poly[d(A-T)] · poly[d(A-T)] and poly[d(G-C)] · poly[d(G-C)] support the suggestion of Klug et al. (1979) that regions of double-stranded DNA containing alternating purine-pyrimidine sequences may exist as structural variants of the classical B-form under physiological salt conditions. Digestion of defined oligomers containing alternating dG-dC sequences indicate that these too exist in some “alternating-B” structure in solution under similar conditions. The results obtained with the oligomers also provide a number of insights into the mode of action of DNAase I.In the case of the B-DNA dodecamer d(C-G-C-G-A-A-T-T-C-G-C-G), for which the crystal structure has been solved (Dickerson &; Drew, 1981), there is a very good correlation between the sites of rapid DNAase I cutting and positions of high local helical twist.     

A sequence analysis system encompassing rules for DNA helical distortion.

Recently, it has been shown by Calladine (1982) and Dickerson (1983) that DNA distortions due to steric clashes between opposing purines and pyrimidines can be quantitated based upon four sum functions. The distortions involve helical twist, roll, torsion angle variations and propeller twist. It is the contention of the authors that these perturbations in structure act as information carriers for various external DNA interactions. This paper describes a system that incorporates these four rules and various other functions that permit the systematic interactive exploration for significant patterns as a consequence of these steric clashes.     

Locally oscillatory motion of RNA helix derived from linear relationships of backbone torsion angles

A linear relationship in each of the torsion angle pairs, α-β, β-?, ?-ζ, and α-γ, has been found by applying a statistical method based on the concept of circular variates to backbone torsion angle data of helical in yeast tTNA^Phe. A series of helical dimer models generated with these relationships have been found to be stereochemically acceptable, and the models also indicate that the backbone unit in the RNA helix is geometrically capable of an oscillatory motion with the distance of about 3.4 Å between adjacent bases. The motion of the backbone unit is analogous to that of a helical spring. The adjacent bases, because of being attached to the backbone, oscillate in a manner similar to the oscillatory dimer model proposed by Davis and Tinoco [Davis, R. C. & Tinoco, I., Jr. (1968) Biopolymers 6 , 223–242]. Here, the oscillation of the backbone unit in the RNA helix is discussed in terms of two geometrical quantities: the torsion (τ) and curvature (κ) of the helix. On these lines, a stereochemical model of RNA strand separation is proposed.     

Refinement of the solution structure of the B DNA hexamer 5'd(C-G-T-A-C-G)2 on the basis of inter-proton distance data

A restrained least-squares refinement of the solution structure of the self-complementary B DNA hexamer 5'd(C-G-T-A-C-G)2 is presented. The structure is refined on the basis of 190 inter-proton distances determined by pre-steady-state nuclear Overhauser enhancement measurements. Two refinements were carried out starting from two initial B DNA structures differing by an overall root-mean-square (r.m.s.) difference of 0.32 A. In both cases, the final r.m.s. difference between the experimental and calculated inter-proton distances was 0.12 A compared to 0.61 A and 0.58 A for the two initial structures. The difference between the two refined structures is small, with an overall r.m.s. difference of 0.16 A, and represents the error in the refined co-ordinates. The refined structures have a B-type conformation with local structural variations in backbone and glycosidic bond torsion angles, and base-pair propellor twist, base roll, base tilt and local helical twist angles.     

Two dimensional NMR studies on the solution structure of d-CTCGAGCTCGAG

Two dimensional (2D) FT-NMR investigations have been carried out on the self-complementary dodecanucleotide d-CTCGAGCTCGAG, which has cleavage sites for the restriction enzyme Xho I (between C and T). The central TCG portion is also known to show a preference for DNAase activity. Complete resonance assignments have been obtained for the non-exchangeable sugar and base protons of the oligonucleotide. Information regarding sugar geometries, glycosidic torsion angles and other structural parameters has been obtained from the relative intensities of the cross peaks in the COSY and NOESY spectra. The results indicate that deoxyribose rings of C1 and C7 adopt a conformation different from the remaining sugars in the double helical oligonucleotide. The central TCG portion also exhibits variations in the backbone structure. The base stacking in the double helix shows interesting sequence dependent effects suggesting that the sequence effects are not localised to nearest neighbours but extended over longer stretches.     

Structural specificities of five commonly used DNA nucleases

Five commonly used nucleases were surveyed for their ability to distinguish among several different DNA backbone configurations. The digestion data suggest that: (1) DNAase I binds across the minor groove; whereas (2) nuclease S1 and (3) micrococcal nuclease bind to an exposed single strand; (4) copper/phenanthroline seeks a base-pair step; and (5) DNAase II requires just a stacked single strand of limited exposure. Only micrococcal nuclease is demonstrably base-specific, with a strong preference for T, A over C, G in any structural context.     

Base sequence and helix structure variation in B and A DNA

The observed propeller twist in base-pairs of crystalline double-helical DNA oligomers improves the stacking overlap along each individual helix strand. But, as proposed by Calladine, it also leads to clash or steric hindrance between purines at adjacent base-pairs on opposite strands of the helix. This clash can be relieved by: (1) decreasing the local helix twist angle between base-pairs; (2) opening up the roll angle between base-pairs on the side on which the clash occurs; (3) separating purines by sliding base-pairs along their long axes so that the purines are partially pulled out of the stack (leading to equal but opposite alterations in main-chain torsion angle delta at the two ends of the base-pair); and (4) flattening the propeller twist of the offending base-pairs. Simple sum functions, sigma 1 through sigma 4, are defined, by which the expected local variation in helix twist, base roll angle, torsion angle delta and propeller twist may be calculated from base sequence. All four functions are quite successful in predicting the behavior of B DNA. Only the helix twist and base roll functions are applicable to A DNA, and the helix twist function begins to fail for an A helical RNA/DNA hybrid. Within these limits, the sequence-derived sum functions match the observed helix parameter variation quite closely, with correlation coefficients greater than 0.900 in nearly all cases. Implications of this sequence-derived helix parameter variation for repressor-operator interactions are considered.     

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.     

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

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.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

DNAase I,DNAase II and staphylococcal nuclease cut at different,yet symmetrically located,sites in the nucleosome core

We have determined the relative location of pancreatic DNAase (DNAase I), spleen acid DNAase (DNAase II) and staphylococcal nuclease cleavage sites in the nucleosome core. Each of these three enzymes cleaves the DNA of chromatin at 10. n nucleotide intervals (n integer); this specificity presumably reflects the internal structure of the nucleosome. We have already reported that DNAase I cleaves nucleosomal DNA so that nearest adjacent cuts on opposite strands are staggered by 2 nucleotides, 3′ end extending (Sollner-Webb and Felsenfeld, 1977). Here we show that the nearest cuts made by DNAase II in nucleosomal DNA are staggered by 4 nucleotides, 3′ end extending, while cuts made by staphylococcal nuclease have a stagger of 2 nucleotides, 5′ end extending. The cutting sites of the three enzymes thus do not coincide. Each pair of staggered cuts, however, is symmetrically located about a common axis-that is, the “dyad axes” that bisect nearest pairs of cutting sites coincide for all three enzymes. This result is consistent with the presence of a true dyad axis in the nucleosome core.Our results support the conclusion that a structural feature of the nucleosome, having a 10 nucleotide periodicity, is the common recognition site for all three nucleases. The position of the cut is determined, however, by the individual characteristics of each enzyme. Sites potentially available to nuclease cleavage span a region of 4 nucleotides out of this 10 nucleotide repeat, and a large fraction of these sites are actually cut. Thus much of the nucleosomal DNA must in some sense be accessible to the environment.     

The conformation of constitutive DNA interaction sites for eukaryotic DNA topoisomerase I on intrinsically curved DNAs.

The analysis of the sites which are cleaved constitutively and preferentially by eukaryotic DNA topoisomerase I on two intrinsically curved DNAs reveals the conformational features that provoke the cleavage reaction on the curve-inducing sequence elements in the absence of supercoiling. This analysis is based on the observation (Caserta et al. (1989) Nucleic Acids Res. 17, 8521-8532 and (1990) Biochemistry 29, 8152-8157) that the reaction of eukaryotic DNA topoisomerase I occurs on two types of DNA sites: sites S (Supercoiled induced) and sites C (Constitutive, whose presence is topology-independent). We report that sites C are abundant on the intrinsically curved DNAs analyzed. The DNAs studied were two intrinsically curved segments of different origin: the Crithidia fasciculata kinetoplast DNA and the bent-containing domain B of the Saccharomyces cerevisiae ARS1. On these DNA segments DNA topoisomerase I cleaves at the junctions between the poly(A) tracts and mixed-sequence DNA. Analysis of the conformation of the double helix around the cleavage sites has revealed that the reaction occurs in correspondence of a defined DNA conformational motif. This motif is described by the set of Eulerian angular values that define the axial path of DNA (helical twist, deflection angle, direction) and of the orthogonal components of wedge (roll and tilt).     

Lethal short limb dwarfism with dysmorphic face, omphalocele and severe ossification defect: Piepkorn syndrome or severe "boomerang dysplasia"?

The authors report a case of lethal neonatal dwarfism characterized by striking micromelia, fused rudimentary and supernumerary digits, large, soft head, pronounced hypertelorism, protruding eyes set laterally, enormous omphalocele and severe deficiency of tubular bone and spine ossification. Histologic examination showed lack of ossification of the cartilaginous anlage of many tubular bones. The cartilage had irregularly distributed chondrocytes. The matrix contained hypocellular and degenerated areas with scattered large chondrocytes. In a few bones a very disorganized growth cartilage was present. The case is similar to that described by Piepkorn et al. (1977) and may represent a severe form of "boomerang dysplasia" (Kozlowski et al., 1981; Tenconi et al., 1983; Kozlowski et al., 1985; Winship et al., 1990).     

Sequence-dependent conformation of an A-DNA double helix. The crystal structure of the octamer d(G-G-T-A-T-A-C-C)

The crystal structures of the synthetic self-complementary octamer d(G-G-T-A-T-A-C-C) and its 5-bromouracil-containing analogue have been refined to R values of 20% and 14% at resolutions of 1.8 and 2.25 A, respectively. The molecules adopt and A-DNA type double-helical conformation, which is minimally affected by crystal forces. A detailed analysis of the structure shows a considerable influence of the nucleotide sequence on the base-pair stacking patterns. In particular, the electrostatic stacking interactions between adjacent guanine and thymine bases produce symmetric bending of the double helix and a major-groove widening. The sugar-phosphate backbone appears to be only slightly affected by the base sequence. The local variations in the base-pair orientation are brought about by correlated adjustments in the backbone torsion angles and the glycosidic orientation. Sequence-dependent conformational variations of the type observed here may contribute to the specificity of certain protein-DNA interactions.