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.     

Solution structure of a five-adenine bulge loop within a DNA duplex.

The three-dimensional solution structure of a DNA molecule of the sequence 5'-d(GCATCGAAAAAGCTACG)-3' paired with 5'-d(CGTAGCCGATGC)-3' containing a five-adenine bulge loop (dA(5)-bulge) between two double helical stems was determined by 2D (1)H and (31)P NMR, infrared, and Raman spectroscopy. The DNA in both stems adopt a classical B-form double helical structure with Watson-Crick base pairing and C2'-endo sugar conformation. In addition, the two dG/dC base pairs framing the dA(5)-bulge loop are formed and are stable at least up to 30 degrees C. The five adenine bases of the bulge loop are localized at intrahelical positions within the double helical stems. Stacking on the double helical stem is continued for the first four 5'-adenines in the bulge loop. The total rise (the height) of these four stacked adenines roughly equals the diameter of the double helical stem. The stacking interactions are broken between the last of these four 5'-adenines and the fifth loop adenine at the 3'-end. This 3'-adenine partially stacks on the other stem. The angle between the base planes of the two nonstacking adenines (A10 and A11) in the bulge loop reflects the kinking angle of the global DNA structure. The neighboring cytosines opposite the dA(5)-bulge (being parts of the bulge flanking base pairs) do not stack on one another. This disruption of stacking is characterized by a partial shearing of these bases, such that certain sequential NOEs for this base step are preserved. In the base step opposite the loop, an extraordinary hydrogen bond is observed between the phosphate backbone of the 5'-dC and the amino proton of the 3'-dC in about two-thirds of the conformers. This hydrogen bond probably contributes to stabilizing the global DNA structure. The dA(5)-bulge induces a local kink into the DNA molecule of about 73 degrees (+/-11 degrees ). This kinking angle and the mutual orientation of the two double helical stems agree well with results from fluorescence resonance energy transfer measurements of single- and double-bulge DNA molecules.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Base pair opening within B-DNA: free energy pathways for GC and AT pairs from umbrella sampling simulations

The conformational pathways and the free energy variations for base opening into the major and minor grooves of a B-DNA duplex are studied using umbrella sampling molecular dynamics simulations. We compare both GC and AT base pair opening within a double-stranded d(GAGAGAGAGAGAG)· d(CTCTCTCTCTCTC) oligomer, and we are also able to study the impact of opening on the conformational and dynamic properties of DNA and on the surrounding solvent. The results indicate a two-stage opening process with an initial coupling of the movements of the bases within the perturbed base pair. Major and minor groove pathways are energetically comparable in the case of the pyrimidine bases, but the major groove pathway is favored for the larger purine bases. Base opening is coupled to changes in specific backbone dihedrals and certain helical distortions, including untwisting and bending, although all these effects are dependent on the particular base involved. Partial opening also leads to well defined water bridging sites, which may play a role in stabilizing the perturbed base pairs.     

Interaction between hydroxystilbamidine and DNA. I. Binding isotherms and thermodynamics of the association.

Isotherms describing the binding of hydroxystilbamidine to DNA and polydeoxyribonucleotides were obtained by means of sedimentation or dialysis experiments and fluorescence measurements, over a large range of ionic strengths, temperatures and base compositions. Two different sets of binding sites are necessary to explain the shapes of the isotherms. The first one is characterized by a higher binding constant, a topological specificity for the A-T pair, exclusion of four base pairs per bound dye molecule, the involvement of two ion-pairs, an almost purely entropic free energy of binding and a large enhancement of the blue fluorescence (450 nm) when the site corresponds to three adjacent A-T pairs. The latter does not present any specificity nor enhancement of fluorescence and only one ion-pair is formed. From the geometry of the dye and its selective binding to a double stranded structure, the hydroxystilbamidine molecule in the first set of sites is likely to be situated in the small groove astride the two complementary strands and slightly distorting the helical structure. The angle of the dye axis with the helix axis has a value close to 47 degrees. No definite explanation could be given for the specific binding of hydroxystilbamidine but the phenolic hydroxyl group is likely to play a major role. The hydroxystilbamidine molecule can be considered as a useful tool for checking the accessibility of the small groove.     

Polyelectrolyte theory and chromatin-DNA quaternary structure: Role of ionic strength and Hl histone

Under the assumptions outlined in this paper and those of Manning's theory of polyelectrolyte screening the electrostatic interaction of histone 1 and DNA is examined and shown to lead to spontaneous bending of DNA. Critical variables for this bending are the ionic strength, the length of DNA interacting with each histone 1 molecule, and the fraction of DNA phosphate charges neutralized by the histone 1 molecule. This interaction is postulated as accounting for the observed folding of polynucleosomes into the “solenoid” regular structure. Wherever possible comparison is made between theoretical predictions and available experimental results. Finally, biological implications are briefly discussed.     

Molecular structure of an A-DNA decamer d(ACCGGCCGGT)

The molecular structure of the DNA decamer d(ACCGGCCGGT) has been solved and refined by single-crystal X-ray-diffraction analysis at 0.20 nm to a final R-factor of 18.0%. The decamer crystallizes as an A-DNA double helical fragment with unit-cell dimensions of a = b = 3.923 nm and c = 7.80 nm in the space group P6(1)22. The overall conformation of this A-DNA decamer is very similar to that of the fiber model for A-DNA which has a large average base-pair tilt and hence a wide and shallow minor groove. This structure is in contrast to that of several A-DNA octamers in which the molecules all have low base-pair-tilt angles (8-12 degrees) resulting in an appearance intermediate between B-DNA and A-DNA. The average helical parameters of this decamer are typical of A-DNA with 10.9 base pairs/turn of helix, an average helical twist angle of 33.1 degrees, and a base-pair-tilt angle of 18.2 degrees. However, the CpG step in this molecule has a low local-twist angle of 24.5 degrees, similar to that seen in other A-DNA oligomers, and therefore appears to be an intrinsic stacking pattern for this step. The molecules pack in the crystal using a recurring binding motif, namely, the terminal base pair of one helix abuts the surface of the shallow minor groove of another helix. In addition, the GC base pairs have large propeller-twist angles, unlike those found most other A-DNA structures.     

NMR analysis of helix I from the 5S RNA of Escherichia coli.

The structure of helix I of the 5S rRNA from Escherichia coli has been determined using a nucleolytic digest fragment of the intact molecule. The fragment analyzed, which corresponds to bases (-1)-11 and 108-120 of intact 5S rRNA, contains a G-U pair and has unpaired bases at its termini. Its proton resonances were assigned by two-dimensional NMR methods, and both NOE distance and coupling constant information have been used to calculate structural models for it using the full relaxation matrix algorithm of the molecular dynamics program XPLOR. Helix I has A-type helical geometry, as expected. Its most striking departure from regular helical geometry occurs at its G-U, which stacks on the base pair to the 5' side of its G but not on the base pair to its 3' side. This stacking pattern maximizes interstrand guanine-guanine interactions and explains why the G-U in question fails to give imino proton NOE's to the base pair to 5' side of its G. These results are consistent with the crystal structures that have been obtained for wobble base pairs in tRNAPhe [Mizuno, H., & Sundaralingam, M. (1978) Nucleic Acids Res. 5, 4451-4461] and A-form DNA [Rabbinovich, D., Haran, T., Eisenstein, M., & Shakked, Z. (1988) J. Mol. Biol. 200, 151-161]. The conformations of the terminal residues of helix I, which corresponds to bases (-1)-11 and 108-120 of native 5S RNA, are less well-determined, and their sugar puckers are intermediate between C2' and C3'-endo, on average.     

Anisotropic flexibility of DNA and the nucleosomal structure.

Potential energy calculations of the DNA duplex dimeric subunit show that the double helix may be bent in the direction of minor and major grooves much more easily than in other directions. It is found that the total winding angle of DNA decreases upon such bending. A new model for DNA folding in the nucleosome is proposed on the basis of these findings according to which the DNA molecule is kinked each fifth base pair to the side of the minor and major grooves alternatively. The model explains the known contradiction between a C-like circular dichroism for the nucleosomal DNA and the nuclease digestion data, which testify to the B-form of DNA.     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be -0.5 degrees per base pair (bp). This is equivalent to a change of -0.15 +/- 0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55 +/- 0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38 +/- 0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained .     

The response of DNA length and twist to changes in ionic strength

We examine twist‐stretch coupling of unconstrained DNA using polyelectrolyte theory as applied to a line‐charge model along with published data on the ionic‐strength dependence of the twist angle. We conclude that twist‐stretch coupling is negative: environmental changes that stretch free DNA, unconstrained by externally applied pulling or twisting forces, are accompanied by unwinding of the double helix. We also analyze a helical model and conclude that the observed unwinding of the DNA helix when ionic strength is decreased is driven by radial swelling of the helix. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 223–226, 2015.     

Studies on the structure of chemically methylated DNA

CD and melting temperature measurements on the nature of DNA with chemically methylated guanine-rich sites indicate that the stable secondary structure of DNA depicted by Ramstein et al- involves considerable distortions resulting from decreased base-base stacking interaction. Besides that quantum chemical data gained from PPP calculations are in favor of a weaker hydrogen bonding interaction in the methylated guanine-cytosine base pair. CD measurements demonstrate that methylated DNA-regions differ from the nonmethylated helical structure, since formation of a condensed conformation as occurs in the transition from B to the C-uke structure is prevented by positively charged methylated guanine residues. An increase in helix winding angle, however, can not be excluded. Binding ability of the dyes acridine orange, phenosafranine, and the antibiotic actinomycin C is lowered for methylated DNA, while binding of proflavine is, in accordance with the results of Ramstein and Leng, slightly enhanced. The reason for the opposite behavior of proflavine is at present not fully understood. In particular changes in the binding ability with dyes could not be correlated with base specificity of complex formation. It is discussed that structural changes in DNA towards a loose conformation decrease the binding tendency for acridine orange, phenosafranine, and actinomycin C.     

Dynamic opening of DNA during the enzymatic search for a damaged base

Uracil DNA glycosylase (UDG) removes uracil from U.A or U.G base pairs in genomic DNA by extruding the aberrant uracil from the DNA base stack. A question in enzymatic DNA repair is whether UDG and related glycosylases also use an extrahelical recognition mechanism to inspect the integrity of undamaged base pairs. Using NMR imino proton exchange measurements we find that UDG substantially increases the equilibrium constant for opening of T-A base pairs by almost two orders of magnitude relative to free B-DNA. This increase is brought about by enzymatic stabilization of an open state of the base pair without increasing the rate constant for spontaneous base pair opening. These findings indicate a passive search mechanism in which UDG uses the spontaneous opening dynamics of DNA to inspect normal base pairs in a rapid genome-wide search for uracil in DNA.     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be −0.5° per base pair (bp). This is equivalent to a change of −0.15±0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55±0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38±0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained.     

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.     

Origin of DNA helical structure and its sequence dependence

Conformational analysis of DNA shows that the origin of the B-form double helix can be attributed in large part to the atomic charge pattern in the base pairs. The charge patterns favor specific helical stacking of the base pairs. Base pairs alone--without backbones--have a strong tendency to form helix, indicating that the backbones play a rather passive role in determining the basic helical structure of DNA. It is mainly the electrostatic interactions determined by the charge pattern on base pairs that stabilize a particular helical conformation. The charge pattern in the base pairs appears to be responsible for much of the sequence dependence of DNA conformation, rather than steric clashes.