Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved ∼200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44° and 46°, respectively, located within ∼60 bp “curvature modules” containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.     

Cytosine methylation can induce local distortions in the structure of duplex DNA.

Methyl groups at the C5 position of pyrimidines located within oligopurine-oligopyrimidine tracts in DNA have been shown previously to modulate curvature generated by those tracts. However, it was not known whether the influence of such methyl groups is consequent to the altered helical structure within the tracts themselves. In the current study, it is demonstrated that methylation of cytosines up to three base pairs away from a (dA)5.(dT)5 tract (A-tract) can still result in alterations of the net curvature of the A-tract-containing DNA, as measured by alterations in electrophoretic mobility. This latter effect depends strongly on both the sequence of the non-A-tract DNA and the positions of the methylated C residues. The current results lend further support to the notion that the biological consequences of cytosine methylation may be effected through local alterations in DNA structure as well as through direct protein-DNA interactions.     

A-tract polarity dominate the curvature in flanking sequences

It is well-known, but little understood, that the nucleotide sequences between phased A(4-6)-tracts (at 10-11 bp intervals) have only a slight effect on overall curvature. To explore this phenomenon, we have examined the gel-migration properties of sequences containing both A-tracts as well as G-tracts (i.e., sequences of the form G(n)C(m) or C(n)G(m), n + m > 4) in various relative positioning. We show that the composite bend of these sequences depends on their relative arrangement. When G-tracts are placed between two A-tracts, such that both tracts are repeated in phase to themselves (e.g., G(5)A(6)G(5)A(5)), or adjacent to the 3'-side of A-tracts (e.g., A(6)G(5)N(10)), they have minimal influence on the extent of bending of the composite sequence. When G-tracts are placed one helical repeat away from A-tracts (e.g., G(5)N(5)A(6)N(6)), or are adjacent only to the 5'-side of A-tracts (e.g., G(5)A(6)N(10)) their influence on the composite bend is larger. The differential behavior of AG- versus GA-tracts means that A-tracts influence their flanking sequences in a polar manner. Whereas they suppress, or make constant, the intrinsic bending characteristics of any sequence placed immediately 3' to them (and hence by definition any sequence placed between two phased A-tracts), sequences adjoining them on their 5'-side are free to modulate the overall curvature. We interpret these results as evidence for the dominant nature of the unique and nonuniform structure adopted by tracts of four adenines or more. The effects of A-tracts extend at least five base pairs into the adjoining 3' region. This is further evidence for the complexity of DNA structure and the inadequacy of simple nearest-neighbor models to explain all its manifestations.     

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.     

2-Pyrimidinone as a probe for studying the EcoRII DNA methyltransferase-substrate interaction.

EcoRII DNA methyltransferase (M.EcoRII) recognizes the 5' em leader CC*T/AGG em leader 3' DNA sequence and catalyzes the transfer of the methyl group from S-adenosyl-l-methionine to the C5 position of the inner cytosine residue (C*). Here, we study the mechanism of inhibition of M.EcoRII by DNA containing 2-pyrimidinone, a cytosine analogue lacking an NH(2) group at the C4 position of the pyrimidine ring. Also, DNA containing 2-pyrimidinone was used for probing contacts of M.EcoRII with functional groups of pyrimidine bases of the recognition sequence. 2-Pyrimidinone was incorporated into the 5' em leader CCT/AGG em leader 3' sequence replacing the target and nontarget cytosine and central thymine residues. Study of the DNA stability using thermal denaturation of 2-pyrimidinone containing duplexes pointed to the influence of the bases adjacent to 2-pyrimidinone and to a greater destabilizing influence of 2-pyrimidinone substitution for thymine than that for cytosine. Binding of M.EcoRII to 2-pyrimidinone containing DNA and methylation of these DNA demonstrate that the amino group of the outer cytosine in the EcoRII recognition sequence is not involved in the DNA-M.EcoRII interaction. It is probable that there are contacts between the functional groups of the central thymine exposed in the major groove and M.EcoRII. 2-Pyrimidinone replacing the target cytosine in the EcoRII recognition sequence forms covalent adducts with M.EcoRII. In the absence of the cofactor S-adenosyl-l-methionine, proton transfer to the C5 position of 2-pyrimidinone occurs and in the presence of S-adenosyl-l-methionine, methyl transfer to the C5 position of 2-pyrimidinone occurs.     

Evidence from base-pair kinetics for two types of adenine tract structures in solution: their relation to DNA curvature

We have measured the base-pair lifetimes in oligodeoxynucleotides containing tracts of A.T base pairs using imino proton magnetic resonance. When the tract contains more than four consecutive A.T base pairs, possibly including a 5'-AT step but not a 5'-TA step, anomalously long lifetimes are observed. For example, the lifetimes of the central A.T base pairs of the dodecamer 5'-d-CGCAAAAAAGCG are 122 and 91 ms at 15 degrees C whereas, in the same conditions, the lifetime of the central A.T pair of the decamer 5'-d-CGCGATCGCG is only 4 ms, a value similar to those measured in several other B-DNA oligoduplexes [Leroy et al. (1988) J. Mol. Biol. 200, 223-238]. This strongly suggests that, in tracts of four A.T pairs or more, a conformation distinct from standard B-DNA is formed cooperatively. All sequences known to generate curved DNA exhibit anomalously long base-pair lifetimes. This is the first local and physical property shown to correlate with DNA curvature. Our observations suggest that the structure responsible for the long lifetimes is involved in the curvature of DNA.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

Structure of the N6-adenine DNA methyltransferase M.TaqI in complex with DNA and a cofactor analog

The 2.0 A crystal structure of the N6-adenine DNA methyltransferase M.TaqI in complex with specific DNA and a nonreactive cofactor analog reveals a previously unrecognized stabilization of the extrahelical target base. To catalyze the transfer of the methyl group from the cofactor S-adenosyl-l-methionine to the 6-amino group of adenine within the double-stranded DNA sequence 5'-TCGA-3', the target nucleoside is rotated out of the DNA helix. Stabilization of the extrahelical conformation is achieved by DNA compression perpendicular to the DNA helix axis at the target base pair position and relocation of the partner base thymine in an interstrand pi-stacked position, where it would sterically overlap with an innerhelical target adenine. The extrahelical target adenine is specifically recognized in the active site, and the 6-amino group of adenine donates two hydrogen bonds to Asn 105 and Pro 106, which both belong to the conserved catalytic motif IV of N6-adenine DNA methyltransferases. These hydrogen bonds appear to increase the partial negative charge of the N6 atom of adenine and activate it for direct nucleophilic attack on the methyl group of the cofactor.     

Uranyl photoprobing of nonbent A/T- and bent A-tracts. A difference of flexibility?

In this study, we have systematically compared the uranyl photocleavage of a range of bent A-tracts and nonbent TA-tracts as well as interrupted A-tracts. We demonstrate that uranyl photocleavage of A-tracts and TA-tracts is almost identical, indicating a very similar minor groove conformation. Furthermore, a 10 base pair A-tract is divided into two independent tracts by an intervening TA or GC step. Uranyl probing also clearly distinguishes the bent A4T4 and the nonbent T4A4 sequences as adopting different structures, and our interpretation of the data is consistent with a structure for the bent A4T4 sequence that resembles a continuous A-tract, whereas the nonbent T4A4 sequences are closer to two independent and opposite A-tracts that cancel each other in terms of macroscopic bending. Finally, we also note that even single TA and TAT steps are highly sensitive to uranyl photocleavage and propose that in addition to average minor groove width, uranyl also senses DNA helix flexibility/deformability. Thus, the structural difference of TA-tracts and A-tracts may to a large extent reflect a difference in flexibility, and DNA curvature may consequently require a rigid narrow minor groove conformation that creates distinct A-tract-B-DNA junctions as the predominant cause of the bending.     

M.HhaI binds tightly to substrates containing mismatches at the target base.

The (cytosine-5) DNA methyltransferase M.HhaI causes its target cytosine base to be flipped completely out of the DNA helix upon binding. We have investigated the effects of replacing the target cytosine by other, mismatched bases, including adenine, guanine, thymine and uracil. We find that M.HhaI binds more tightly to such mismatched substrates and can even transfer a methyl group to uracil if a G:U mismatch is present. Other mismatched substrates in which the orphan guanine is changed exhibit similar behavior. Overall, the affinity of DNA binding correlates inversely with the stability of the target base pair, while the nature of the target base appears irrelevant for complex formation. The presence of a cofactor analog. S-adenosyl-L-homocysteine, greatly enhances the selectivity of the methyltransferase for cytosine at the target site. We propose that the DNA methyltransferases have evolved from mismatch binding proteins and that base flipping was, and still is, a key element in many DNA-enzyme interactions.     

Human DNA polymerase kappa bypasses and extends beyond thymine glycols during translesion synthesis in vitro,preferentially incorporating correct nucleotides

Human polymerase kappa (polkappa), the product of the human POLK (DINB1) gene, is a member of the Y superfamily of DNA polymerases that support replicative bypass of chemically modified DNA bases (Ohmori, H., Friedberg, E. C., Fuchs, R. P., Goodman, M. F., Hanaoka, F., Hinkle, D., Kunkel, T. A., Lawrence, C. W., Livneh, Z., Nohmi, T., Prakash, L., Prakash, S., Todo, T., Walker, G. C., Wang, Z., and Woodgate, R. (2001) Mol. Cell 8, 7-8; Gerlach, V. L., Aravind, L., Gotway, G., Schultz, R. A., Koonin, E. V., and Friedberg, E. C. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 11922-11927). Polkappa is shown here to bypass 5,6-dihydro-5,6-dihydroxythymine (thymine glycol) generated in two different DNA substrate preparations. Polkappa inserts the correct base adenine opposite thymine glycol in preference to the other three bases. Additionally, the enzyme correctly extends beyond the site of the thymine glycol lesion when presented with adenine opposite thymine glycol at the primer terminus. However, steady state kinetic analysis of nucleotides incorporated opposite thymine glycol demonstrates different misincorporation rates for guanine with each of the two DNA substrates. The two substrates differ only in the relative proportions of thymine glycol stereoisomers, suggesting that polkappa distinguishes among stereoisomers and exhibits reduced discrimination between purines when incorporating a base opposite a 5R thymine glycol stereoisomer. When extending beyond the site of the lesion, the misincorporation rate of polkappa for each of the three incorrect nucleotides (adenine, guanine, and thymine) is dramatically increased. Our findings suggest a role for polkappa in both nonmutagenic and mutagenic bypass of oxidative damage.     

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.     

Improved quantitation of DNA curvature using ligation ladders.

It is often desirable to estimate accurately the local shape of DNA molecules. Such measurements are useful in understanding the intrinsic contribution of DNA sequence to curvature, as well as in assessing the effects of chemical modifications. We have been investigating the effects of asymmetric phosphate neutralization on DNA shape using the well-characterized ligation ladder approach developed by Crothers and co-workers [D.M.Crothers and J.Drak (1992) Meth. Enzymol.,212, 46-71]. This technique is remarkably sensitive to differences in DNA shape. We now report a general quantitative assay of DNA curvature that we have validated using a set of phased A(5)tract standards. This approach allows simultaneous estimation of helix axis deflection magnitude and direction when a test sequence is monitored in at least three phasings relative to a reference A(5-6)tract in short DNA duplexes. Analysis using this improved approach confirms our published data on DNA curvature due to electrostatic effects.     

Detection of an unusual distortion in A-tract DNA using KMnO4: effect of temperature and distamycin on the altered conformation.

The chemical probes potassium permanganate (KMnO4) and diethylpyrocarbonate (DEPC) can be used to study the conformational flexibility of short tracts of adenine (A-tracts) present in DNA. With these probes, we demonstrate that a novel distortion is induced in a 5 base pair A-tract at low temperature. Formation of this distorted A-tract structure, which occurs in a DNA fragment from the promoter region of the plasmid pBR322, is distinguished by a dramatic increase in the KMnO4 reactivity of the central thymines in this tract at 12 degrees C. This alteration occurs in the absence of any detectable rearrangement in the conformation of the adenines in the complementary strand. Induction of this low temperature A-tract structure is blocked by the minor groove binding drug distamycin. Hydroxyl radical footprinting of distamycin binding to the fragment containing the d(A)5 tract at 12 degrees C suggests that this drug has two different modes of binding to DNA in agreement with recent NMR data. These experiments show that short A-tracts are capable of forming more than one structural variant of B DNA in solution. The possible relationship between the intrinsic bending of DNA containing short phased A-tracts and the low temperature A-tract conformation is discussed.     

Structural constraints regulating triple helix formation by A-tracts

The study concerns the propensity of triple helix formation by different DNA oligonucleotides containing long A-tracts with and without flanking GxC base pairs in order to probe the role of length of the A-tract and the flanking sequences. From nuclear magnetic resonance (NMR) studies of imino proton spectra and circular dichroism (CD) spectroscopy of samples composed of potential triplex forming strand sequences in correct stoichiometries, we have concluded that 8-mer A-tracts flanked by GxC base pairs exert significant steric hindrance to triple helix formation. When as much as 50 mM Mg2+ was added, no triple helix formation was observed in these samples. In contrast, open-ended 8-mer A-tracts formed triplex with the corresponding two T8 strands under relatively mild ionic conditions (100 mM Na+). Moreover, the shorter the length of the A-tract, the less is the hindrance to form a triple helix.     

DNA duplexes containing photoactive derivatives of 2'-deoxyuridine as photocrosslinking probes for EcoRII DNA methyltransferase-substrate interaction

EcoRII DNA methyltransferase (M.EcoRII) recognizes the DNA sequence 5'.CC*T/AGG.3' and catalyzes the transfer of the methyl group from S-adenosyl-L-methionine to the C5 position of the inner cytosine residue (C*). We obtained several DNA duplexes containing photoactive 5-iodo-2'-deoxyuridine (i(5)dU) or 5-[4-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenyl]-2'-deoxyuridine (Tfmdp-dU) to characterize regions of M.EcoRII involved in DNA binding and to investigate the DNA double helix conformational changes that take place during methylation. The efficiencies of methylation, DNA binding affinities and M.EcoRII-DNA photocrosslinking yields strongly depend on the type of modification and its location within the EcoRII recognition site. The data obtained agree with the flipping of the target cytosine out of the DNA double helix for catalysis. To probe regions of M.EcoRII involved in DNA binding, covalent conjugates M.EcoRII-DNA were cleaved by cyanogen bromide followed by analysis of the oligonucleotide-peptides obtained. DNA duplexes containing i(5)dU or Tfmdp-dU at the central position of the recognition site, or instead of the target cytosine were crosslinked to the Gly(268)-Met(391) region of the EcoRII methylase. Amino acid residues from this region may take part both in substrate recognition and stabilization of the extrahelical target cytosine residue.     

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.     

Increased temperature and 2-methyl-2,4-pentanediol change the DNA structure of both curved and uncurved adenine/thymine-rich sequences

DNA curvature is affected by elevated temperature and dehydrating agents such as 2-methyl-2,4-pentanediol (MPD) (used in crystallization). This effect of MPD has been ascribed to a specific distortion of the structure of adenine tracts (A-tracts), probably through a deformation of the characteristic narrow minor groove. Uranyl photoprobing indicates that a narrowed minor groove is present in all A/T regions containing four or more A/T base pairs. Consequently, this technique may be employed to study conformational changes in other A/T-rich sequences than pure A-tracts. In this study we use uranyl photoprobing to demonstrate that the effect of elevated temperature and MPD is analogous on both "normal" and curve-inducing A/T-rich sequences. The results therefore indicate that under these conditions the minor groove is widened in all A/T sequences and not only in pure A-tracts as previously suggested. Thus, the rather subtle structural difference of AT regions and A-tracts in nonbent DNA versus A-tracts in bent DNA may be quantitative rather than qualitative; i.e., the structure is more persistent and/or rigid in bent DNA.     

A structural analysis of the bent kinetoplast DNA from Crithidia fasciculata by high resolution chemical probing.

The chemical probes potassium permanganate (KMnO4) and diethylpyrocarbonate (DEPC) have been used to study the conformation of bent kinetoplast DNA from Crithidia fasciculata at different temperatures. Chemical reactivity data shows that the numerous short A-tracts of this bent DNA adopt a similar structure at 43 degrees C. This conformation appears to be very similar to the conformation of A-tracts in DNA exhibiting normal gel mobility. The A-tract structure detected by chemical probing is characterized by a high degree of base stacking on the thymine strand, and by an abrupt conformational change at the 3' end of the adenine strand. In general, no major alteration of this A-tract specific structure was detected between 4-53 degrees C. However, probing with KMnO4 revealed two unusual features of the C. fasciculata sequence that may contribute to the highly aberrant gel mobility of this DNA: 1) the B DNA/A-tract junction 5' dC/A3-6 3'. 5' dT3-6/G 3' is disproportionately represented and is conformationally distinct from other 5' end junctions, and 2) low temperature favors a novel strand-specific conformational distortion over a 20 base pair region of the bent kinetoplast DNA. Presence of the minor groove binding drug distamycin had little detectable effect on the A-tract conformation. However, distamycin did inhibit formation of the novel KMnO4 sensitive low temperature structure and partially eliminated the anomalous gel mobility of the kinetoplast DNA. Finally, we describe a simple and reproducible procedure for the production of an adenine-specific chemical DNA sequence ladder.