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.     

A critical role for linker DNA in higher-order folding of chromatin fibers

Nucleosome-nucleosome interactions drive the folding of nucleosomal arrays into dense chromatin fibers. A better physical account of the folding of chromatin fibers is necessary to understand the role of chromatin in regulating DNA transactions. Here, we studied the unfolding pathway of regular chromatin fibers as a function of single base pair increments in linker length, using both rigid base-pair Monte Carlo simulations and single-molecule force spectroscopy. Both computational and experimental results reveal a periodic variation of the folding energies due to the limited flexibility of the linker DNA. We show that twist is more restrictive for nucleosome stacking than bend, and find the most stable stacking interactions for linker lengths of multiples of 10 bp. We analyzed nucleosomes stacking in both 1- and 2-start topologies and show that stacking preferences are determined by the length of the linker DNA. Moreover, we present evidence that the sequence of the linker DNA also modulates nucleosome stacking and that the effect of the deletion of the H4 tail depends on the linker length. Importantly, these results imply that nucleosome positioning in vivo not only affects the phasing of nucleosomes relative to DNA but also directs the higher-order structure of chromatin.     

Studies of DNA dumbbells. I. Melting curves of 17 DNA dumbbells with different duplex stem sequences linked by T4 endloops: evaluation of the nearest-neighbor stacking interactions in DNA.

Seventeen DNA dumbbells were constructed that have duplex sequences ranging in length from 14 to 18 base pairs linked on the ends by T4 single-strand loops. Fifteen of the molecules have the core duplexes with the sequences 5'G-T-A-T-C-C-(W-X-Y-Z)-G-G-A-T-A-C3', where (W-X-Y-Z) represents a unique combination of A.T, T.A, G.C, and C.G base pairs. The remaining two molecules have the central sequence (W-X-Y-Z) = A-C and A-C-A-C-A-C. These duplex sequences were designed such that the central sequences include different combinations of the 10 possible nearest-neighbor (n-n) stacks in DNA. In this sense the set of molecules is complete and serves as a model system for evaluating sequence-dependent local stability of DNA. Optical melting curves of the samples were collected in 25, 55, 85, and 115 mM [Na+], and showed, regardless of solvent ionic strength, that the transition temperatures of the dumbbells vary by as much as 14 degrees for different molecules of the set. Results of melting experiments analyzed in terms of a n-n sequence-dependent model allowed evaluation of nine independent linear combinations of the n-n stacking interactions in DNA as a function of solvent ionic strength. Although there are in principle 10 possible different n-n interactions in DNA, these 10 are not linearly independent and therefore can not be uniquely determined. For molecules with ends, there are 9 linearly independent combinations, as opposed to circular or semiinfinite repeating copolymers where only 8 linear combinations of the 10 possible n-n interactions are linearly independent. The n-n interactions are presented as combinations of the deviations from average stacking for the 5'-3' base-pair doublets, delta Gi, and reveal several interesting features: (1) Titratable changes in the values of delta Gi with changing salt environment are observed. In all salts the most stable unique combination is delta G4 = (delta GGpC+delta GCpG)/2, and the least stable is the GpG/CpC stack, delta G2 = delta GGpG/CpC. (2) The chi 2 values of the fits of the evaluated delta Gi's to experimental data increased with decreasing [Na+], suggesting that significant interactions beyond nearest neighbors become more pronounced, particularly at 25 nM Na+.(ABSTRACT TRUNCATED AT 400 WORDS)     

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.     

DNA binding properties of the yeast Msh2-Msh6 and Mlh1-Pms1 heterodimers

We describe here our recent studies of the DNA binding properties of Msh2-Msh6 and Mlh1-Pms1, two protein complexes required to repair mismatches generated during DNA replication. Mismatched DNA binding by Msh2-Msh6 was probed by mutagenesis based on the crystal structure of the homologous bacterial MutS homodimer bound to DNA. The results suggest that several amino acid side chains inferred to interact with the DNA backbone near the mismatch are critical for repair activity. These contacts, which are different in Msh2 and Msh6, likely facilitate stacking and hydrogen bonding interactions between side chains in Msh6 and the mismatched base, thus stabilizing a kinked DNA conformation that permits subsequent repair steps coordinated by the Mlh1-Pms1 heterodimer. Mlh1-Pms1 also binds to DNA, but independently of a mismatch. Mlh1-Pms1 binds short DNA substrates with low affinity and with a slight preference for single-stranded DNA. It also binds longer duplex DNA molecules, but with a higher affinity indicative of cooperative binding. Indeed, imaging by atomic force microscopy reveals cooperative DNA binding and simultaneous interaction with two DNA duplexes. The novel DNA binding properties of Mlh1-Pms1 may be relevant to signal transduction during DNA mismatch repair and to recombination, meiosis and cellular responses to DNA damage.     

Properties of the nucleic-acid bases in free and Watson-Crick hydrogen-bonded states: computational insights into the sequence-dependent features of double-helical DNA

The nucleic-acid bases carry structural and energetic signatures that contribute to the unique features of genetic sequences. Here we review the connection between the chemical structure of the constituent nucleotides and the polymeric properties of DNA. The sequence-dependent accumulation of charge on the major- and minor-groove edges of the Watson-Crick base pairs, obtained from ab initio calculations, presents unique motifs for direct sequence recognition. The optimization of base interactions generates a propellering of base-pair planes of the same handedness as that found in high-resolution double-helical structures. The optimized base pairs also deform along conformational pathways, i.e., normal modes, of the same type induced by the binding of proteins. Empirical energy computations that incorporate the properties of the base pairs account satisfactorily for general features of the next level of double-helical structure, but miss key sequence-dependent differences in dimeric structure and deformability. The latter discrepancies appear to reflect factors other than intrinsic base-pair structure.     

Three-dimensional structures of bulge-containing DNA fragments.

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.     

Structure and stability of DNA containing an aristolactam II-dA lesion: implications for the NER recognition of bulky adducts

Aristolochic acids I and II are prevalent plant toxicants found in the Aristolochiaceae plant family. Metabolic activation of the aristolochic acids leads to the formation of a cyclic N-hydroxylactam product that can react with the peripheral amino group of purine bases generating bulky DNA adducts. These lesions are mutagenic and established human carcinogens. Interestingly, although AL-dG adducts progressively disappear from the DNA of laboratory animals, AL-dA lesions has lasting persistence in the genome. We describe here NMR structural studies of an undecameric duplex damaged at its center by the presence of an ALII-dA adduct. Our data establish a locally perturbed double helical structure that accommodates the bulky adduct by displacing the counter residue into the major groove and stacking the ALII moiety between flanking bases. The presence of the ALII-dA perturbs the conformation of the 5'-side flanking base pair, but all other pairs of the duplex adopt standard conformations. Thermodynamic studies reveal that the lesion slightly decreases the energy of duplex formation in a sequence-dependent manner. We discuss our results in terms of its implications for the repair of ALII-dA adducts in mammalian cells.     

A computational and experimental study of the bending induced at a double-triple helix junction.

We have studied the conformation of a 17 base-pair homopyrimidine.homopurine triple helix formed on a fragment of duplex DNA derived from Simian Virus SV40. Gel retardation assays indicate that an 80 base-pair fragment has an altered conformation when the triple helix is formed, which is most likely to result from an induced bend in the DNA. Investigation of the detailed conformation of the double helix-triple helix junctions has been performed by means of molecular modelling. Bending on the 5' and 3' sides of the third strand oligonucleotide are not located at equivalent positions with respect to the junctions, which is explained in terms of base stacking. The junction effects on DNA structure, induced by the requirement for cytosine protonation in the Hoogsteen-bonded strand to form CGC+ base triplets, are also discussed.     

Stacking and stability of RNA duplexes containing fluorobenzene and fluorobenzimidazole nucleosides

Six different fluorobenzene or fluorobenzimidazole ribonucleosides and one abasic site were incorporated in oligoribonucleotides. Individual contributions of base stacking and solvation of the modified nucleosides could be determined. In fluorobenzene.fluorobenzimidazole-modified base pairs a duplex stabilizing force was found that points to a weak F...H hydrogen bond. The lipophilicity of the unprotected nucleosides were investigated by determination of 1-octanol water partition coefficients.     

The dynamic structural basis of differential enhancement of conformational stability by 5'- and 3'-dangling ends in RNA

Unpaired bases at the end of an RNA duplex (dangling ends) can stabilize the core duplex in a sequence-dependent manner and are important determinants of RNA folding, recognition, and functions. Using 2-aminopurine as a dangling end purine base, we have employed femtosecond time-resolved fluorescence spectroscopy, combined with UV optical melting, to quantitatively investigate the physical and structural nature of the stacking interactions between the dangling end bases and the terminal base pairs. A 3'-dangling purine base has a large subpopulation that stacks on the guanine base of the terminal GC or UG pair, either intrastrand or cross-strand depending on the orientation of the pair, thus providing stabilization of different magnitudes. On the contrary, a 5'-dangling purine base only has a marginal subpopulation that stacks on the purine of the same strand (intrastrand) but has little cross-strand stacking. Thus a 5'-dangling purine does not provide significant stabilization. These stacking structures are not static, and a dangling end base samples a range of stacked and unstacked conformations with respect to the terminal base pair. Femtosecond time-resolved anisotropy decay reveals certain hindered base conformational dynamics that occur on the picosecond to nanosecond time scales, which allow the dangling base to sample these substates. When the dangling purine is opposite to a U and is able to form a potential base pair at the end of the duplex, there is an interplay of base stacking and hydrogen-bonding interactions that depends on the orientation of the base pair relative to the adjacent GC pair. By resolving these populations that are dynamically exchanging on fast time scales, we elucidated the correlation between dynamic conformational distributions and thermodynamic stability.     

STACKING AND STABILITY OF RNA DUPLEXES CONTAINING FLUOROBENZENE AND FLUOROBENZIMIDAZOLE NUCLEOSIDES

Six different fluorobenzene or fluorobenzimidazole ribonucleosides and one abasic site were incorporated in oligoribonucleotides. Individual contributions of base stacking and solvation of the modified nucleosides could be determined. In fluorobenzene · fluorobenzimidazole-modified base pairs a duplex stabilizing force was found that points to a weak F…H hydrogen bond. The lipophilicity of the unprotected nucleosides were investigated by determination of 1-octanol water partition coefficients.     

Fluorescent alpha-anomeric 1,N(6)etheno-deoxyadenosine in DNA duplexes. The alpha-epsilondA / dG pair

Structural properties of the fluorescent alpha-anomeric 1,N(6)ethenodeoxyadenosine residue placed in opposition to all four canonical deoxynucleotide units within 11-mer DNA duplexes have been studied. The duplex with alpha-epsilondA / dG pairing is most thermodynamically stable while the alpha-epsilondA / dC one is the least stable. Fluorescence measurements confirm the thermodynamic data and indicate base-pair dependent stacking properties of alpha-epsilondA within duplex structures. Results of molecular dynamics (MD) simulations in aqueous solution for the most stable duplex point to the presence of different conformational states of the alpha-1,N(6)etheno-deoxyadenosine residue, including formation of a hydrogen bonded pair with the dG and possible occurrence of severe kinking in the duplex.     

Nucleosomes are positioned on mouse satellite DNA in multiple highly specific frames that are correlated with a diverged subrepeat of nine base-pairs

Nucleosome phasing on highly repetitive DNA was investigated using a novel strategy. Nucleosome cores were prepared from mouse liver nuclei with micrococcal nuclease, exonuclease III and nuclease S1. The core DNA population that contains satellite sequences was then purified from total core DNA by denaturation of the DNA, reassociation to a low Cot value and hydroxyapatite chromatography to separate the renatured satellite fraction. After end-labeling, the termini of the satellite core DNA fragments were mapped with an accuracy of +/- 1 base-pair relative to known restriction sites on the satellite DNA. Sixteen dominant nucleosome positions were detected. There is a striking correlation between these nucleosome frames and an internal highly diverged 9 base-pair subrepeat of the satellite DNA. The results are consistent with a sequence-dependent association of histone octamers with the satellite DNA. Our finding that histone octamers can interact with a given DNA in a number of different defined frames has important implications for the possible biological significance of nucleosome phasing.     

UV resonance Raman and circular dichroism studies of a DNA duplex containing an A(3)T(3) tract: evidence for a premelting transition and three-centered H-bonds.

The presence of A(n) and A(n)T(n) tracts in double-helical sequences perturbs the structural properties of DNA molecules, resulting in the formation of an alternate conformation to standard B-DNA known as B'-DNA. Evidence for a transition occurring prior to duplex melting in molecules containing A(n) tracts was previously detected by circular dichroism (CD) and calorimetric studies. This premelting transition was attributed to a conformational change from B'- to B-DNA. Structural features of A(n) and A(n)T(n) tracts revealed by X-ray crystallography include a large degree of propeller twisting of adenine bases, narrowed minor grooves, and the formation of three-centered H-bonds between dA and dT bases. We report UV resonance Raman (UVRR) and CD spectroscopic studies of two related DNA dodecamer duplexes, d(CGCAAATTTGCG)(2) (A(3)T(3)) and d(CGCATATATGCG)(2) [(AT)(3)]. These studies address the presence of three-centered H-bonds in the B' conformation and gauge the impact of these putative H-bonds on the structural and thermodynamic properties of the A(3)T(3) duplex. UVRR and CD spectra reveal that the premelting transition is only observed for the A(3)T(3) duplex, is primarily localized to the dA and dT bases, and is associated with base stacking interactions. Spectroscopic changes associated with the premelting transition are not readily detectable for the sugar-phosphate backbone or the cytosine and guanosine bases. The temperature-dependent concerted frequency shifts of dA exocyclic NH(2) and dT C4=O vibrational modes suggest that the A(3)T(3) duplex forms three-centered hydrogen bonds at low temperatures, while the (AT)(3) duplex does not. The enthalpy of this H-bond, estimated from the thermally induced frequency shift of the dT C4=O vibrational mode, is approximately 1.9 kJ/mol or 0.46 kcal/mol.     

Solution structure of an 11-mer duplex containing the 3, N(4)-ethenocytosine adduct opposite 2'-deoxycytidine: implications for the recognition of exocyclic lesions by DNA glycosylases

Lipid peroxidation products, as well as the metabolic products of vinyl chloride, react with cellular DNA producing the mutagenic adduct 3,N(4)-etheno-2'-deoxycytidine (epsilondC), along with several other exocyclic derivatives. High-resolution NMR spectroscopy and restrained molecular dynamics simulations were used to establish the solution structure of an 11-mer duplex containing an epsilondC.dC base-pair at its center. The NMR data suggested a regular right-handed helical structure having all residues in the anti orientation around the glycosydic torsion angle and Watson-Crick alignments for all canonical base-pairs of the duplex. Restrained molecular dynamics generated a three-dimensional model in excellent agreement with the spectroscopic data. The (epsilondC. dC)-duplex structure is a regular right-handed helix with a slight bend at the lesion site and no severe distortions of the sugar-phosphate backbone. The epsilondC adduct and its partner dC were displaced towards opposite grooves of the helix, resulting in a lesion-containing base-pair that was highly sheared but stabilized to some degree by the formation of a single hydrogen bond. Such a sheared base-pair alignment at the lesion site was previously observed for epsilondC.dG and epsilondC.T duplexes, and was also present in the crystal structures of duplexes containing dG.T and dG. U mismatches. These observations suggest the existence of a substrate structural motif that may be recognized by specific DNA glycosylases during the process of base excision repair.     

DNA curvature does not require bifurcated hydrogen bonds or pyrimidine methyl groups.

Short tracts of the homopolymer dA.dT confer intrinsic curvature on the axis of the DNA double helix. This phenomenon is assumed to be a consequence of such tracts adopting a stable B'-DNA conformation that is distinct from B-form structure normally assumed by other DNA sequences. The more stable B' structure of dA.dT tracts has been attributed to several possible stabilizing factors: (1) optimal base stacking interactions consequent upon the high propeller twist, (2) bifurcated hydrogen bonds between adjacent dA.dT base-pairs, (3) stacking interactions involving the dT methyl groups, and finally (4) a putative spine of ordered water molecules in the minor groove. DNA oligodeoxynucleotides have been synthesized that enable these hypotheses to be tested; of particular interest is the combination of effects due to bifurcation (2) and methylation of the pyrimidines nucleotides (3). The data indicate that neither bifurcated hydrogen bonds nor pyrimidine methyl groups nor both are essential for DNA curvature. The data further suggest that the influence of the minor groove spine of hydration on the B'-formation is small. The experiments favor the hypothesis that base stacking interactions are the dominant force in stabilizing the B'-form structure.