Structure of the deoxytetranucleotide d-pApTpApT and a sequence-dependent model for poly(dA-dT)

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.

DNA coiled coils

We report the coiled-coil structure of DNA, which is generated by the dodecanucleotide d(ATATATATATAT). The structure has been determined by single-crystal x-ray crystallography. The molecules form duplexes with single-stranded overhangs, which associate with neighbor molecules and give rise to infinite double helices in a coiled-coil conformation, with staggered nicks in both strands. The coiled coils have six dodecamer duplexes per turn. Despite the presence of nicks, the structure is very rigid. Statistical disorder is present, which gives rise to continuous scattering along the layer lines. This observation can be interpreted by the classical theory of coiled coils developed for proteins, which we apply here to a DNA structure. A clear splitting of the original layer lines of the DNA double helix is detected. The structure we have found adds a previously unrecognized element to the architectures that can be built from DNA oligonucleotides. Any duplex with complementary single-strand overhangs should be expected to give rise to regular coiled-coil structures.

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.

Thymine-methyl/pi interaction implicated in the sequence-dependent deformability of DNA

The crystal structures of deoxy-oligonucleotides were retrieved from the Nucleic Acid Database and analyzed with the use of our program CHPI. The structure of 5'-ApTpApT-3' has been shown to be stabilized by the 5-methyl group in the thymidine moiety that favorably interacts with the adenine pi-ring preceding it. H2' of the deoxyribose in adenine also interacts with the thymine ring next to it. Since a 5'-ApT-3' sequence is accompanied by another 5'-ApT-3' in the complementary strand, the interaction is duplicated, thus forming a 'twin A/T-Me interaction'. Coordinates of oligonucleotides with A-T rich sequences were retrieved and analyzed. In almost every case, the thymidine 5-methyl group favorably interacts with an adenine ring in the same strand. The structure of duplexes incorporating A-tracts was also analyzed. The 5-methyl group in the thymidine moiety has been found to interact favorably with the base pi-ring before it. Since an A-tract is lined with an oligo-T sequence in the complementary strand, a successive N/T-Me stacking may contribute in making the A-tracts robust and straight. The possible involvement of the N/T-Me and the twin A/T-Me motif in the deformability of DNA has been suggested. The role of methyl groups in modified DNA has been discussed on a similar basis.

Crystal structure of an antiparallel DNA fragment with Hoogsteen base pairing

We report here an alternative double-helical structure of the DNA molecule. It has been found in the d(ATA(Br)UAT) and d(ATATAT) sequences by single-crystal x-ray crystallography. This sequence is found not only in TATA boxes, but also in other regulatory regions of DNA. Bases of the two antiparallel strands form Hoogsteen pairs, with adenines in the syn conformation. The structure is related neither to those found in triple helices nor to parallel DNA duplexes. Its conformational parameters are very similar to those of duplex DNA in the B form. Both forms may coexist under physiological conditions, although the Hoogsteen pairing greatly influences the recognition sites on DNA. Our results demonstrate that an alternative to the classical B-DNA double helix is possible.

Analysis of nucleic acids by on-line liquid chromatography-mass spectrometry

The numerous problems posed by modern biochemistry, biology, and medicine, as well as the growing significance of genetic engineering require the application of fast and reliable methods of utmost sensitivity and selectivity for the analysis of nucleic acids. High-performance liquid chromatography (HPLC) and mass spectrometry (MS) represent established analytical techniques for the characterization and structural elucidation of single- and double-stranded nucleic acids, ranging in size from a few nucleotides to several thousand base pairs. Although both techniques are independently applicable for nucleic acid analysis, the on-line hyphenation significantly enhances their potential for the robust and fully automable routine analysis of minute amounts of biological samples. Among the various chromatographic and mass spectrometric modes available in principle, ion-pair reversed-phase HPLC and electrospray ionization mass spectrometry (ESI-MS) have been shown to be the most suitable for the direct interfacing of liquid chromatography (LC) and MS. Instrumental setup, as well as chromatographic and mass spectrometric experimental conditions, need to be carefully selected in order to maximize the performance of the hyphenated analytical system. Applications of HPLC-ESI-MS include the characterization of oligodeoxynucleotides synthesized by solid-phase synthesis, the analysis of antisense oligodeoxynucleotides, oligonucleotide metabolites, and DNA adducts, the analysis of genomic segments specifically amplified by the polymerase chain reaction (PCR), the characterization of ribonucleic acids, the sizing of double-stranded DNA restriction fragments, the genotyping of short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs), the detection of mutations in nucleic acid sequences, and the sequencing of nucleic acids.

On the metastability of left-handed DNA motifs

Alternating purine-pyrimidine DNA sequences such as poly [d(C-G)] or poly[d(m5C-G)] undergo a cooperative, salt-induced structural transition from a right-handed B conformation, which prevails at relatively low ionic strength, into a left-handed Z form, generally believed to be stabilized by high salt concentrations. We report here that upon a monotonous increase of the ionic strength, the well-established B to Z transition is followed by a second, hitherto unobserved conformational change leading from Z-DNA back into a right-handed B-like form. This observation indicates that, in contrast with the current convention, the Z motif represents an unstable configuration relative to the B form at both low and high salt concentrations and that the occurrence of a left-handed DNA structure, presently depicted as a step function of the ionic strength, should rather be treated in terms of a pulse. The reported transition underscores the inherent metastability of the Z configuration, and indicates, consequently, that this motif is ideally suited to act as a structural regulatory element, as such an element should be endowed with a large susceptibility to cellular parameters.

A structural model for sequence-specific proflavin-DNA interactions during in vitro frameshift mutagenesis

Molecular models describing intermediates that may lead to proflavin-induced 1 bp deletions during in vitro polymerization by E. coli DNA polymerase I Klenow fragment are proposed. The models provide structural explanations for the fact that the induced frameshifts always occur opposite template bases that are adjacent to 5' pyrimidines and are based on the underlying hypothesis that the deletions arise because the polymerase passes by a template base without copying it. Because the most frequent mutations are opposite Pu in the template sequence 5' Py Pu 3', a single-strand loop-out model was constructed for this sequence and proflavin was added, using structures found in crystalline oligonucleotides and their complexes with proflavin. The model seeks to rationalize the roles of the 5' pyrimidine and proflavin in facilitating the bypass. Four potential roles for proflavin in mutagenesis are described: 1) stacking on the looped-out base; 2) stacking on the base pair immediately preceding the site of mutation; 3) hydrogen bonding with the 5' pyrimidine; 4) hydrogen bonding with the phosphate backbone. These models point to the possibility that a number of proflavin-DNA interactions may be involved. In contrast, modeling does not suggest a role for classically intercalated proflavin in frameshift mutagenesis arising during in vitro DNA polymerization.

Hydration of DNA bases: analysis of crystallographic data

We present a systematic analysis of water structure around nucleic acid bases. We have examined 28 crystal structures of oligonucleotides, and have studied the patterns of water around the four bases, guanine, cytosine, adenine, and thymine. The geometries of water positions were calculated up to 4.00 A from base atoms. We have found conformation-dependent differences in both the geometry and extent of hydration of the bases.

Investigations on the dynamic structures of adenine- and thymine-containing DNA

The structures of poly(dA-dT), poly(dA-dBr5U) and of poly(dA).poly(dT) have been investigated in solution and in fibers, by Raman spectroscopy. Both the alternating poly(dA-dT), poly(dA-dBr5U) and non-alternating poly(dA).poly(dT) exhibit, in the region of sugar phosphate backbone vibrations, two bands of almost equal intensity at about 841 cm-1 and 817 cm-1. The analysis of the characteristic bands of thymine residues that are sensitive to sugar puckers gives indication of a significant displacement from the C(2')-endo conformer suggesting the adoption of alternative conformers such as O(4')-endo. In contrast, the diagnostic Raman bands for the sugar pucker of adenine residues suggest, instead, predominant adoption of C(2')-endo conformations. These Raman results are compatible with rapid dynamic changes of sugar puckers between C(2')-endo and O(4')-endo for the thymidine (and uridine) residues, whereas in adenine residues the sugar puckers fluctuate around the C(2')-endo pucker in all synthetic DNA molecules studied. Molecular dynamics simulations, performed on six different starting models using two distance-dependent dielectric functions epsilon(r) = 4 r and a sigmoidal dependence), all gave similar dynamic behavior in agreement with these Raman data and their interpretation. The mean calculated pseudorotation phases of the adenine residues are systematically higher (around C(2')-endo) than those of the thymine residues (close to O(4')-endo-C(1')-exo). Besides, the mean lifetimes of the thymine residues are 1.5 to 2.0-fold higher in the O(4')-endo than in the C(2')-endo domain, while those of the adenine residues are two to threefold higher in the C(2')-endo than in the O(4')-endo domain. In the Raman spectra of the alternating poly(dA-dBr5U), the splitting of a band into two components arising from the two contributions of ApBr5U and Br5UpA provides strong evidence for a repeating dinucleotide structure in solution. The calculated twist values averaged over the simulation runs are also systematically higher in the 5'T-A3' step (39 degrees) than in the 5'A-T3' step (33 degrees). Simultaneously, the calculated roll values are positive in the 5'T-A3' step (6 degrees) and negative in the 5'A-T3' step (-9 degrees), while the propeller twist values are about the same (-11 degrees to -16 degrees). On the other hand, in the homopolymer, the average twist value is close to 36 degrees with the roll angle close to 0 degrees and large propeller twist values (-20 degrees).

Time-resolved fluorescence emission and excitation spectroscopy of d(TA) and d(AT) using synchrotron radiation

The photophysics of the sequence isomers d(TA) and d(AT) has been investigated at room temperature in 5 x 10(-5) M neutral aqueous solution using pulsed ultraviolet excitation from the ACO synchrotron and detection by time correlation or gated single-photon counting. Decay profiles of the emissions at 350, 400 and 460 have been analyzed both independently and globally by reiterative non-linear least-squares fitting to models of two and three independently emitting species. No evidence has been observed for excited-state reaction. Time-windowed spectra, both emission and excitation, have been collected for three time windows and have been deconvoluted to give time-resolved spectra using the lifetimes resulting from the decay analyses. Spectra are separated into two classes, with picosecond and nanosecond lifetimes, respectively. The picosecond spectra have the emission and excitation spectral characteristics of mixed monomer (A and T) fluorescences and are assigned as originating from the unstacked fractions of d(TA) and d(AT). The nanosecond emission spectra from d(TA) and d(AT) are both two-component, with lambda max approximately 350 and approximately 425 nm and lifetimes of 2.3 and 6.1 ns, respectively. The time-resolved excitation spectra for the nanosecond emissions are quite different from the isotropic absorption spectra of d(TA) and d(AT) but correlate with the anisotropic absorption for out-of-plane transitions between stacked bases of co-crystals of 9-methyladenine and 1-methylthymine reported by Stewart and Davidson. The nanosecond spectra thus represent the direct excitation and emission of stacked pairs of bases. These results provide no evidence for energy transfer and are probably related to sequence-specific photo-adduct formation.

Local supercoil-stabilized DNA structures

The DNA double helix exhibits local sequence-dependent polymorphism at the level of the single base pair and dinucleotide step. Curvature of the DNA molecule occurs in DNA regions with a specific type of nucleotide sequence periodicities. Negative supercoiling induces in vitro local nucleotide sequence-dependent DNA structures such as cruciforms, left-handed DNA, multistranded structures, etc. Techniques based on chemical probes have been proposed that make it possible to study DNA local structures in cells. Recent results suggest that the local DNA structures observed in vitro exist in the cell, but their occurrence and structural details are dependent on the DNA superhelical density in the cell and can be related to some cellular processes.

Molecular modeling and energy refinement of supercoiled DNA

A method is presented for constructing the complete atomic structure of supercoiled DNA starting from a linear description of the double helical pathway. The folding pathway is defined by piecewise B-spline curves and the atoms are initially positioned with respect to the local Frenet trihedra determined by the equations of the curves. The resulting chemical structure is corrected and refined with an energy minimization procedure based on standard potential expressions. The refined molecular structure is then used to study the effects of supercoiling on the local secondary structure of DNA. The minimized structure is found to differ from an isotropic elastic rod model of the double helix, with the base pairs bending in an asymmetric fashion along the supercoiled trajectory. The starting trajectory is chosen so that the refined supercoiled structure is either underwound (10.37 base pairs per turn) or overwound (9.65 base pairs per turn) compared to the standard tenfold B-DNA fiber diffraction model. The underwound supercoil is also lower in energy than the overwound duplex. The variation of base pair sequence in poly(dA).poly(dT).poly(dAT).poly(dTA) and poly(dA5T5).poly(dT5A5) is additionally found to influence the secondary structural features along a given supercoiled pathway. Finally, the detailed features of the refined structures are found to be in agreement with known X-ray crystallographic structures of DNA oligomers.

Oligonucleotide structure: a decade of results from single crystal X-ray diffraction studies

Writing a review gives authors a splendid opportunity to view developments in a particular area of science from a very personal angle. They are at liberty to select material, emphasize aspects of direct interest to their own work and air speculations which, wisely or not, referees have caused to have removed from their publications. Such personal accounts often make good reading, but may be somewhat misleading especially for readers seeking an introduction to the field. One remedy is to aim at a comprehensive review with equal weight given to all publications but boredom, if not bias, is then likely to creep in.

Structure of an alternating-B DNA helix and its relationship to A-tract DNA

The crystal structure of the synthetic DNA dodecamer CGCATATATGCG has been solved at 2.2-A resolution. Its central 6 base pairs adopt the alternating-B-DNA helix structure proposed nearly a decade ago. This alternating poly(AT) structure contrasts with the four known examples of what can be termed a poly(A) subfamily of B-DNA structures: CGCAAAAAAGCG, CGCAAATTTGCG, CGCGAATTCGCG, and CGCGAATTbrCGCG, their defining characteristic being a succession of two or more adenines along one strand, in a region of 4 or more A.T base pairs. All five helices show a characteristically narrow minor groove in their AT centers, but the mean propeller twist at A.T base pairs is lower in the alternating poly(AT) helix than in the poly(A) subfamily of helices. Three general principles emerge from x-ray analyses of B-DNA oligonucleotides: (i) GC and mixed-sequence B-DNA have a wide minor groove, whereas the minor groove is narrow in heteropolymer or homopolymer AT sequences. (ii) G.C base pairs have low propeller twist; A.T pairs can adopt a high propeller twist but need not do so. A high propeller twist can be stabilized by cross-strand hydrogen bonds in the major or minor groove, examples being the minor groove bonds seen in CCAAGATTGG and the major groove bonds that can accompany AA sequences in the poly(A) family. (iii) Homopolymer poly(A) tracts may be stiffer than are alternating AT or general-sequence DNA because of these cross-strand major groove hydrogen bonds. Poly(A) tracts appear internally unbent, but bends may occur at junctions with mixed-sequence DNA because of differences in propeller twist, base pair inclination, and base stacking on the two sides of the junction. Bending occurs most easily via base roll, favoring compression of the broad major groove.

Relative stabilities and transitions of DNA conformations in 1:1 electrolytes: a theoretical study

We use a recently developed formalism (1) to calculate the salt dependent part of the free energy determining DNA conformational stability in 1:1 electrolytes. The conformations studied are the A, B, C and alternating-B right-handed forms and the ZI, ZII left-handed forms of DNA. In the case of the B-ZI transition of d(G-C).d(G-C) helices in NaCl solution, the free energy contribution considered suffices to describe the transition in a quantitative manner. The theory also predicts the occurrence of salt-induced B-A transitions which have been recently observed with poly[d(n2 A-T)] and poly[d(G-C)]. In other cases, additional terms in the free energy balance, particularly due to hydration effects, must be at least as important as salt effects in determining conformational stability and structural transitions in solution. If diffuse ionic cloud electrostatic effects alone would dominate in all cases, the relative helical stabilities at 0.2 M monovalent salt would decrease in the order C greater than B greater than A greater than ZII greater than ZI greater than alternating-B. At high salt concentrations (2.0 M-5.0 M), the order would be alternating-B greater than ZI greater than A greater than ZII greater than B greater than C.

Local and overall conformations of DNA double helices with the A - T base pairs

Raman spectra have been observed of two different poly[d(A-T)].poly[d(A-T)] fibers, whose X-ray diffractions were confirmed to be purely of A and B forms. A number of spectral differences were found between the A and B forms of this DNA duplex, and they were ascribed to local conformational differences in the adenosine, thymidine and phosphodiester portions. The ascription was made on the basis of a separate series of Raman examinations on six crystals involving adenosine or thymidine, and fifteen other nucleotide crystals, whose structures are all known by previous crystallographic works. By taking these structure-spectrum correlations thus obtained into account, a Raman spectroscopic investigation was made of a few double-helical DNAs in aqueous solutions. It has been concluded that both poly[d(A-T)].poly[d(A-T)] and poly(dA).poly(dT) have a C2'endo-anti adenosine, C2'endo-anti thymidine, a b-type mainchain (beta = 160 +/- 15 degrees, gamma = 45 +/- 15 degrees, delta = 140 +/- 10 degrees) and an a2-type mainchain (beta = 210 +/- 10 degrees, gamma = 45 +/- 15 degrees, delta = 140 +/- 10 degrees) not only in low-salt medium but also in 6.6 M CsF solutions, where beta, gamma and delta are the torsion angles around O5'-C5', C5'-C4' and C4'-C3' axes, respectively. Poly(rA).poly(dT), on the other hand, was considered to have a heteronomous duplex structure, in which the poly(rA) strand has a C3'endo-anti adenosine and a1-type mainchain (beta = 175 +/- 25 degrees, gamma = 45 +/- 15 degrees, delta = 80 +/- 10 degrees) whereas the poly(dT) strand has a C2'endo-anti thymidine and b-type mainchain.

Conformational studies of nucleic acids: III. Empirical multiple correlation functions for nucleic acid torsion angles

There are seven significantly variable torsion angles in each monomer unit of a polynucleotide. Because of this, it is computationally infeasible to consider the energetics of all conformations available to a nucleic acid without the use of simplifications. In this paper, we develop functions suggested by and regression fit to crystallographic data which allow three of these torsion angles, alpha (O3'-P-O5'-C5'), delta (C5'-C4'-C3'-O3') and epsilon (C4'-C3'-O3'-P), to be calculated as dependent variables of those remaining. Using these functions, the seven independent torsions are reduced to four, a reduction in complexity sufficient to allow an examination of the global conformational energetics of a nucleic acid for the remaining independent torsion angles. These functions are the first to quantitatively relate a dependent nucleic acid torsion angle to several different independent angles. In all three cases the data are fit reasonably well, and in one case, alpha, the fit is exceptionally good, lending support for the suitability of the functions in conformational searches. In addition, an examination of the most significant terms in each of the correlation functions allows insight into the physical basis for the correlations.

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.

Nuclease recognition of an alternating structure in a d(AT)14 plasmid insert

The nuclease reactivity and specificity of a cloned tract of poly X (dA-dT) X poly(dA-dT) has been explored. Digestion with DNAse I, Mung Bean nuclease, S1 nuclease, DNAse II, and copper (1,10-phenanthroline)2 on a 256 base pair restriction fragment containing d(AT)14A revealed a dinucleotide repeat structure for the alternating sequence. Furthermore, conditions which wind or unwind the linear DNA had little effect on the reactivity of the AT insert. These preferred cleavages offer insights to structural alterations within the DNA helix which differ from A, B, or Z-DNA. Nucleation into flanking sequences by this structural alteration was not observed.

Ordered water structure in an A-DNA octamer at 1.7 A resolution

The crystal structure of the deoxyoctamer d(G-G-Br U-A-BrU-A-C-C) was refined to a resolution of 1.7 A using combined diffractometer and synchrotron data. The analysis was carried out independently in two laboratories using different procedures. Although the final results are identical the comparison of the two approaches highlights potential problems in the refinement of oligonucleotides when only limited data are available. As part of the analysis the positions of 84 solvent molecules in the asymmetric unit were established. The DNA molecule is highly solvated, particularly the phosphate-sugar back-bone and the functional groups of the bases. The major groove contains, in the central BrU-A-BrU-A region, a ribbon of water molecules forming closed pentagons with shared edges. These water molecules are linked to the base O and N atoms and to the solvent chains connecting the O-1 phosphate oxygen atoms on each strand. The minor groove is also extensively hydrated with a continuous network in the central region and other networks at each end. The pattern of hydration is briefly compared with that observed in the structure of a B-dodecamer.

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 4. The torsion angle distribution about the C3'-O3' bond in DNA constituents

Carbon-13 and proton NMR spectra of a series of oligodeoxynucleotides (d(CT), d(CC), d(TA), d(AT), d(CG), d(GC), d(AG), d(AAA), d(TATA) and d(GGTAAT] were measured at various temperatures. The three coupling constants that are related to the magnitude of backbone angle epsilon (J(C4'-P), J(C2'-P) and J(H3'-P] are analyzed in terms of a three-state equilibrium about this bond. Two epsilon (trans) angles occur, which differ in magnitude depending on the conformation (N or S) of the adjoining deoxyribose ring. The S-type deoxyribose ring is associated with a smaller epsilon (trans) angle: epsilon (t,S) = 192 degrees. The N-type deoxyribose ring is associated with a larger epsilon (trans) angle epsilon (t,N) = 212 degrees. The third rotamer participating in the conformational equilibrium, is a gauche(-) (epsilon (-] conformer and occurs exclusively in combination with the S-type sugar ring (epsilon (-,S) = 266 degrees). Within the limits of experimental error, the magnitude of these three angles appears to be independent of the particular base sequence, except in the case of d(CG) where a slightly larger epsilon (t,S) angle (197 degrees) is indicated. A simple equation is proposed which may be used to calculate the population of epsilon (t,S) conformer in cases where only J(H3'-P) is known.

The flexibility of alternating dA-dT sequences

The flexibility of alternating poly (dA-dT) has been investigated by the technique of transient electric dichroism. Rotational relaxation times, which are very sensitive to changes in the end-to-end length of flexible polymers, are determined from the field free dichroism decay curves of four, well defined fragments of poly (dA-dT) ranging in size from 136 to 270 base pairs. Persistence lengths, calculated from the results of Hagerman and Zimm (Biopolymers (1981) 29, 1481-1502), are in the range 200-250 A. This makes alternating dA-dT sequences about twice as flexible as naturally occurring, "random" sequence DNA. Considering a bend around a nucleosome, for example, this difference in persistence length translates to an energy difference between poly (dA-dT) and random sequence DNA of 0.17 kT/base pair or 1 kcal per 10 base pair stretch. This energy difference is sufficiently large to suggest that dA-dT sequences could serve as markers in DNA packaging, for example, at sites where DNA must tightly bend to accommodate structures.

Sequence-dependent structural variation in single-helical DNA. Proton NMR studies of d(T-A-T-A) and d(A-T-A-T) in aqueous solution

The two deoxyribotetranucleoside triphosphates d(T-A-T-A) and d(A-T-A-T) were investigated in aqueous solution by one- and two-dimensional proton NMR at 300 and 500 MHz. It is demonstrated that both compounds occur predominantly in the single-helical form. Accurate coupling constants are obtained by computer simulation of several 500-MHz spectra. The data are interpreted in terms of N and S pseudorotational ranges. The geometry of the major S-type conformers displays a clear sequence dependence, as expressed by variation of the endocyclic backbone angle delta (C5'-C4'-C3'-O3'). A simple sum rule is proposed to predict delta variation in single-helical DNA fragments. Comparisons are made with other sequence-dependent geometries as observed in a double-helical B-DNA fragment in the crystalline state. Furthermore, one- and two-dimensional nuclear Overhauser effect (NOE) spectroscopy was carried out on d(T-A-T-A). An inventory is made of the observed intra- and inter-residue NOEs. The NOE data confirm the presence of a highly stacked single-helical conformation of d(T-A-T-A) in solution. No indications are found for the formation of a bulge-out structure as observed for analogous alternating purine-pyrimidine oligoribonucleotides.

Base-stacking interactions in double-helical DNA structures: experiment versus theory

Atom-atom potential energy calculations have been undertaken for deriving stacking energies in double-helical structures. A comparison between the energy patterns of A- and B-type double-helical fragments determined by single-crystal X-ray diffraction methods versus idealized uniform models based on fibre diffraction data shows that the van der Waals stacking energy is largely sensitive to local changes in the relative orientation of adjacent base pairs. The sequence-dependent conformational variability observed in the high-resolution structures appears to be a consequence of the equipartitioning of the stacking energy along the double helix. The large energy variations expected for a uniform structure are dampened considerably in the observed structures by means of local changes in conformational features such as helix rotation and roll angles between base pairs.

Deoxydodecanucleotide heteroduplex d(TTTTATAATAAA). d(TTTATTATAAAA) containing the promoter Pribnow sequence TATAAT. I. Double-helix stability by UV spectrophotometry and calorimetry

Thermally induced structural transition in the d(TTTTATAATAAA) d(TTTATTATAAAA) heteroduplex is characterized by UV-spectroscopy and differential scanning calorimetry. At low salt (less than 0.1 M) the occurrence of a cooperative transition in the lower temperature range, followed by a broad transition connected with small increase in absorbance is observed. At high salt (greater than or equal to 0.2 M) a single, monophasic transition appears. Linear dependence of the latter on log of salt concentration (dTm:dlogM = 14.2 degrees C) and of 1/Tm on log of oligomer concentration [derived therefrom delta H (v.H.) = 77.1 kcal/mole (duplex)] allows relating it to the melting of the heteroduplex helix. The non-cooperative transition, independent of oligomer concentration and similar to that of the single chain, was attributed to melting of short hairpin helices upon heteroduplex dissociation. Calorimetric enthalpy: 75.6 kcal/mole (duplex) proved significantly lower than predicted from known calorimetric data for poly[d(AT)] and poly d(A) X poly d(T).

Local mobility of nucleic acids as determined from crystallographic data. I. RNA and B form DNA

The local mobility of DNA and RNA can be described well by a segmented rigid body model in which the phosphates, riboses and bases are treated as independent groups. We have developed a computer program for extracting information about the translational and rotational mobility of these groups from X-ray diffraction data of single crystals of DNA and RNA fragments. We plan to extend our studies from the B DNA helix studied here to A and Z form helices for which diffraction data are available. The mobilities described here may be important in allowing the flexibility necessary for interaction of nucleic acids with proteins, ligands or other nucleic acids, and are crucial to our understanding of the factors governing both structural stability and motion of nucleic acids on a larger scale.

NMR studies of conformational states and dynamics of DNA

The application of high resolution NMR techniques to the investigation of DNA double helices in solution is currently in a rapid state of change as a result of advances in three different fields. First, new methods (cloning, enzymatic degradation, sonication, and chemical synthesis) have been developed for producing large quantities of short DNA suitable for NMR studies. Second, there have been major advances in the field of NMR in terms of the introduction of new pulse techniques and improvements in instrumentation. Finally, as a result of recent X-ray diffraction studies on short DNA helices and the discovery of left-handed Z-DNA there is heightened interest in the study of DNA structures in solution and the effect of sequence on structure. In the present review, we discuss the way in which NMR techniques have been used to probe various aspects of the DNA properties, including base pairing structure, dynamics of breathing, effect of sequence on DNA structure, internal molecular motions, the effect of environment on the DNA, and the interaction of DNA with small ligands.