Relationship between curved DNA conformations and slow gel migration

We propose some specific DNA conformations that explain, in terms of molecular conformations, the anomalous gel electrophoretic behavior of the sequences (VA4T4X), and (V2A3T3X2)i where V and X are either G or C. Previously (J. Biomole. Struct. Dyn. 4, 41, 1986) we considered hydrophobic interactions among aliphatic hydrocarbon groups in A/T sequences. In the sequences (T)n.(A)n, the T's are slightly bent to yield structures with tightly stacked methyl groups along one side of the major groove. By folding together the two pairs of stacked methyls on the opposite sides of the major groove. TTAA might yield a relatively sharp bend. On this basis, we show below that the sequences (VT4A4X)i might form a very tightly coiled super-helix whereas the sequences (VA4T4X)i form a broad super-helix of radius approximately 120 A for i = 25. The sequence (V2A3T3X2)i forms a slightly smaller radius super-helix. The time of passage through the gel has been taken to be inversely proportional to the smallest dimension of the molecule. Specifically we are taking the ratio of the apparent molecular weight to the actual molecular weight to be related to the moment of inertia I1 about the smallest principal axis of the molecular conformation. We find a good fit to the experimental gel mobility data of Hagerman (2) if we assume this ratio to be proportional to (I1)1/5.     

Relationship Between Curved DNA Conformations and Slow Gel Migration

Abstract

We propose some specific DNA conformations that explain, in terms of molecular conformations, the anomalous gel electrophoretic behavior of the sequences (VA₄T₄X)₁, and (V₂A₃X₂)₁ where V and X are either G or C. Previously (J. Biomole. Struct. Dyn. 4, 41, 1986) we considered hydrophobic interactions a mong aliphatic hydrocarbon groups in A/T sequences. In the sequences (T)_n · (A)_n, the T's are slightly bent to yield structures with tightly stacked methyl groups along one side of the major groove. By folding together the two pairs of stacked methyls on the opposite sides of the major groove, TTAA might yield a relatively sharp bend. On this basis, we show below that the sequences (VT₄A₄X)₁ might form a very tightly coiled super-helix whereas the sequences (VA₄T₄X)₁ form a broad super-helix of radius ~ 120 A for i = 25. The sequence (V₂A₃T₃X₂)₁ forms a slightly smaller radius super-helix. The time of passage through the gel has been taken to be inversely proportional to the smallesuiimension of the molecule. Specifically we are taking the ratio of the apparent molecular weight to the actual molecular weight to be related to the moment of inertia I₁ about the smallest principal axis of the molecular conformation. We find a good fit to the experimental gel mobility data of Hagerman (2) if we assume this ratio to be proportional to (I₁)^1/5.     

From footprint to toeprint: a close-up of the DnaA box, the binding site for the bacterial initiator protein DnaA.

The Escherichia coli DnaA protein binds as a monomer to the DnaA box, a 9 bp consensus sequence: 5'-TTA/TTNCACA. To assess the contribution of individual bases to protein binding we probed the DnaA-DnaA box complex with the uracil-DNA glycosylase (UDG) footprinting technique. (i) dU at the positions of T2, T4, T7' or T9' completely inhibits DnaA binding to the DnaA box. At these positions the methyl groups of the thymine residues are essential for successful DnaA binding, indicating protein contact with the major groove. Additionally they are positioned exactly on one side of the helix. (ii) dU at the position of T1 or at three T residues adjacent to the 9 bp core sequence of the DnaA box allows DnaA binding. These positions are protected from UDG digestion as revealed by the footprint assay. (iii) dU at the position of T3' on the complementary strand of teh box 5'-TTATCCACA was not protected from UDG digestion in DNA-DnaA complexes. Therefore, DnaA cannot contact the major groove at this position. In addition, a slight bend of the DnaA box towards UDG would help the enzyme to access this site.     

Structure and folding of a rare, natural kink turn in RNA with an A*A pair at the 2b*2n position

The kink turn (k-turn) is a frequently occurring motif, comprising a bulge followed by G•A and A•G pairs that introduces a sharp axial bend in duplex RNA. Natural k-turn sequences exhibit significant departures from the consensus, including the A•G pairs that form critical interactions stabilizing the core of the structure. Kt-23 found in the small ribosomal subunit differs from the consensus in many organisms, particularly in the second A•G pair distal to the bulge (2b•2n). Analysis of many Kt-23 sequences shows that the frequency of occurrence at the 2n position (i.e., on the nonbulged strand, normally G in standard k-turns) is U>C>G>A. Less than 1% of sequences have A at the 2n position, but one such example occurs in Thelohania solenopsae Kt-23. This sequence folds only weakly in the presence of Mg²⁺ ions but is induced to fold normally by the binding of L7Ae protein. Introduction of this sequence into the SAM-I riboswitch resulted in normal binding of SAM ligand, indicating that tertiary RNA contacts have resulted in k-turn folding. X-ray crystallography shows that the T. solenopsae Kt-23 adopts a standard k-turn geometry, making the key, conserved hydrogen bonds in the core and orienting the 1n (of the bulge-proximal A•G pair) and 2b adenine nucleobases in position facing the opposing minor groove. The 2b and 2n adenine nucleobases are not directly hydrogen bonded, but each makes hydrogen bonds to their opposing strands.     

Characterization of divalent cation localization in the minor groove of the A(n)T(n) and T(n)A(n) DNA sequence elements by (1)H NMR spectroscopy and manganese(II)

The localization of Mn(2+) in A-tract DNA has been studied by (1)H NMR spectroscopy using a series of self-complementary dodecamer oligonucleotides that contain the sequence motifs A(n)(n) and T(n)A(n), where n = 2, 3, or 4. Mn(2+) localization in the minor groove is observed for all the sequences that have been studied, with the position and degree of localization being highly sequence-dependent. The site most favored for Mn(2+) localization in the minor groove is near the 5'-most ApA step for both the T(n)A(n) and the A(n)T(n) series. For the T(n)A(n) series, this results in two closely spaced symmetry-related Mn(2+) localization sites near the center of each duplex, while for the A(n)T(n) series, the two symmetry-related sites are separated by as much as one half-helical turn. The degree of Mn(2+) localization in the minor groove of the T(n)A(n) series decreases substantially as the AT sequence element is shortened from T(4)A(4) to T(2)A(2). The A(n)T(n) series also exhibits length-dependent Mn(2+) localization; however, the degree of minor groove occupancy by Mn(2+) is significantly less than that observed for the T(n)A(n) series. For both A(n)T(n) and T(n)A(n) sequences, the 3'-most AH2 resonance is the least broadened of the AH2 resonances. This is consistent with the observation that the minor groove of A-tract DNA narrows in the 5' to 3' direction, apparently becoming too narrow after two base pairs for the entry of a fully hydrated divalent cation. The results that are reported illustrate the delicate interplay that exists between DNA nucleotide sequence, minor groove width, and divalent cation localization. The proposed role of cation localization in helical axis bending by A-tracts is also discussed.     

Dimeric DNA quadruplex containing major groove-aligned A-T-A-T and G-C-G-C tetrads stabilized by inter-subunit Watson-Crick A-T and G-C pairs

We report on an NMR study of unlabeled and uniformly 13C,15N-labeled d(GAGCAGGT) sequence in 1 M NaCl solution, conditions under which it forms a head-to-head dimeric quadruplex containing sequentially stacked G-C-G-C, G-G-G-G and A-T-A-T tetrads. We have identified, for the first time, a slipped A-T-A-T tetrad alignment, involving recognition of Watson-Crick A-T pairs along the major groove edges of opposing adenine residues. Strikingly, both Watson-Crick G-C and A-T pairings within the direct G-C-G-C and slipped A-T-A-T tetrads, respectively, occur between rather than within hairpin subunits of the dimeric d(GAGCAGGT) quadruplex. The hairpin turns in the head-to-head dimeric quadruplex involve single adenine residues and adds to our knowledge of chain reversal involving edgewise loops in DNA quadruplexes. Our structural studies, together with those from other laboratories, definitively establish that DNA quadruplex formation is not restricted to G(n) repeat sequences, with their characteristic stacked uniform G-G-G-G tetrad architectures. Rather, the quadruplex fold is a more versatile and robust architecture, accessible to a range of mixed sequences, with the potential to facilitate G-C-G-C and A-T-A-T tetrad through major and minor groove alignment, in addition to G-G-G-G tetrad formation. The definitive experimental identification of such major groove-aligned mixed A-T-A-T and G-C-G-C tetrads within a quadruplex scaffold, has important implications for the potential alignment of duplex segments during homologous recombination.     

Molecular dynamics of spermine-DNA interactions: sequence specificity and DNA bending for a simple ligand.

We used molecular dynamics to model interactions between the physiologically important polyamine spermine and two B-DNA oligomers, the homopolymer (dG)10-(dC)10 and the heteropolymer (dGdC)5-(dGdC)5. Water and counterions were included in the simulation. Starting coordinates for spermine-DNA complexes were structures obtained by molecular mechanics modeling of spermine with the two oligomers; in these models, spermine binding induced a bend in the heteropolymer but not in the homopolymer. During approximately 40 psec of molecular dynamics simulation, spermine moves away from the floor of the major groove and interacts nospecifically with d(G)10-d(C)10. In contrast, a spermine-induced bend in the helix of (dGdC)5-(dGdC)5 is maintained throughout the simulation and spermine remains closely associated with the major groove. These results provide further evidence that the binding of spermine to nucleic acids can be sequence specific and that bending of alternating purine-pyrimidine sequences may be a physiologically important result of spermine binding.     

Large,sequence-dependent effects on DNA conformation by minor groove binding compounds

To determine what topological changes antiparasitic heterocyclic dications can have on kinetoplast DNA, we have constructed ligation ladders, with phased A5 and ATATA sequences in the same flanking sequence context, as models. Bending by the A5 tract is observed, as expected, while the ATATA sequence bends DNA very little. Complexes of these DNAs with three diamidines containing either furan, thiophene or selenophene groups flanked by phenylamidines were investigated along with netropsin. With the bent A5 ladder the compounds caused either a slight increase or decrease in the bending angle. Surprisingly, however, with ATATA all of the compounds caused significant bending, to values close to or even greater than the A5 bend angle. Results with a mixed cis sequence, which has one A5 and one ATATA, show that the compounds bend ATATA in the same direction as a reference A5 tract, that is, into the minor groove. These results are interpreted in terms of a groove structure for A5 which is largely pre-organized for a fit to the heterocyclic amidines. With ATATA the groove is intrinsically wider and must close to bind the compounds tightly. The conformational change at the binding site then leads to significant bending of the alternating DNA sequence.     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

Crystal structure of the CENP-B protein-DNA complex: the DNA-binding domains of CENP-B induce kinks in the CENP-B box DNA.

The human centromere protein B (CENP-B), one of the centromere components, specifically binds a 17 bp sequence (the CENP-B box), which appears in every other alpha-satellite repeat. In the present study, the crystal structure of the complex of the DNA-binding region (129 residues) of CENP-B and the CENP-B box DNA has been determined at 2.5 A resolution. The DNA-binding region forms two helix-turn-helix domains, which are bound to adjacent major grooves of the DNA. The DNA is kinked at the two recognition helix contact sites, and the DNA region between the kinks is straight. Among the major groove protein-bound DNAs, this 'kink-straight-kink' bend contrasts with ordinary 'round bends' (gradual bending between two protein contact sites). The larger kink (43 degrees ) is induced by a novel mechanism, 'phosphate bridging by an arginine-rich helix': the recognition helix with an arginine cluster is inserted perpendicularly into the major groove and bridges the groove through direct interactions with the phosphate groups. The overall bending angle is 59 degrees, which may be important for the centromere-specific chromatin structure.     

Sequence effects on translesion synthesis of an aminofluorene-DNA adduct: conformational, thermodynamic, and primer extension kinetic studies

The DNA sequence effect is an important structural factor for determining the extent and nature of carcinogen-induced mutational and repair outcomes. In this study, we used two 16-mer template sequences, TG*A [d(5'-CTTCTTG*ACCTCATTC-3')] and CG*A [d(5'-CTTCTCG*ACCTCATTC-3')], to study the impact of the 5'-flanking nucleotide (T vs C) on aminofluorene (AF)-induced stacked (S)/major groove (B)/wedge (W) conformational heterogeneity during a simulated translesion synthesis. In addition, we probed the sequence effect on nucleotide insertion efficiencies catalyzed by the Klenow fragment (exonuclease-deficient) of DNA polymerase I. Our (19)F NMR/ICD/DSC results showed that AF in the CG*A duplex sequence adopts a greater population of S-conformer than the TG*A sequence. We found that the S conformer of CG*A thermodynamically favors insertion of A over C at the lesion site (n). Significant stalling occurred at both the prelesion (n - 1) and lesion (n) sites; however, the effect was more persistent for the S conformer of CG*A than TG*A at the lesion site (n). Kinetics show that relative nucleotide insertion frequencies (f(ins)) were greater for TG*A than the S conformer of CG*A for either dCTP or dATP at the lesion site (n), and the insertion rate was significantly reduced at immediate upstream base pairs (n, n + 1). Taken together, the results provide insight into how the mutagenic AF could exhibit an S/B/W equilibrium in the active site of a polymerase, causing different mutations. This work represents a novel structure-function relationship in which adduct structure is directly linked to nucleotide insertion efficiency in a conformation-specific manner during translesion DNA synthesis.     

Cooperative dimerization of a heterocyclic diamidine determines sequence-specific DNA recognition

In the course of a program aimed at discovering novel DNA-targeted antiparasitic drugs, the phenylfuran-benzimidazole unfused aromatic dication DB293 was identified as the first diamidine capable of forming stacked dimers in the DNA minor groove of GC-containing sequences. Its preferred binding sequence encompasses the tetranucleotide 5'-ATGA.5'-TCAT to which DB293 binds tightly with a strong positive cooperativity. Here we have investigated the influence of the DNA sequence on drug binding using two complementary technical approaches: surface plasmon resonance and DNase I footprinting. The central dinucleotide of the primary ATGA motif was systematically varied to represent all of the eight possible combinations (AXGA and ATYA, where X or Y = A, T, G, or C). Binding affinities for each site were precisely measured by SPR, and the extent of cooperative drug binding was also determined. The sequence recognition process was found to be extremely dependent on the nature of the central dinucleotide pair. Modification of the central TG step decreases binding affinity by a factor varying from 2 to over 500 depending on the base substitution. However, the diminished binding affinity does not affect the unique binding mode. In nearly all cases, the SPR titrations revealed a positive cooperativity in complex formation which reflects the ease of the dication to form stacked dimeric motifs in the DNA minor groove. DNase I footprinting served to identify additional binding sites for DB293 in the context of long DNA sequences offering a large variety of randomly distributed or specifically designed sites. The ATGA motif provided the best receptor for the drug, but lower affinity sequences were also identified. The design of two DNA fragments composed of various targeted tetranucleotide binding sites separated by an "insulator" (nonbinding) sequence allowed us to delineate further the influence of DNA sequence on drug binding and to identify a novel high-affinity site: 5'-ACAA.5'-TTGT. Collectively, the SPR and footprinting results show that the consensus sequence 5'-(A/T)-TG-(A/T) represents the optimal site for cooperative dimerization of the heterocyclic diamidine DB293.     

The extended and eccentric E-DNA structure induced by cytosine methylation or bromination

Cytosine methylation or bromination of the DNA sequence d(GGCGCC)2 is shown here to induce a novel extended and eccentric double helix, which we call E-DNA. Like B-DNA, E-DNA has a long helical rise and bases perpendicular to the helix axis. However, the 3'-endo sugar conformation gives the characteristic deep major groove and shallow minor groove of A-DNA. Also, if allowed to crystallize for a period of time longer than that yielding E-DNA, the methylated sequence forms standard A-DNA, suggesting that E-DNA is a kinetically trapped intermediate in the transition to A-DNA. Thus, the structures presented here chart a crystallographic pathway from B-DNA to A-DNA through the E-DNA intermediate in a single sequence. The E-DNA surface is highly accessible to solvent, with waters in the major groove sitting on exposed faces of the stacked nucleotides. We suggest that the geometry of the waters and the stacked base pairs would promote the spontaneous deamination of 5-methylcytosine in the transition mutation of dm5C-dG to dT-dA base pairs.     

Symmetry elements in DNA structure important for recognition/methylation by DNA [amino]-methyltransferases

The phage T4Dam and EcoDam DNA-[adenine-N6] methyltransferases (MTases) methylate GATC palindromic sequences, while the BamHI DNA-[cytosine-N4] MTase methylates the GGATCC palindrome (which contains GATC) at the internal cytosine residue. We compared the ability of these enzymes to interact productively with defective duplexes in which individual elements were deleted on one chain. A sharp decrease in k_cat was observed for all three enzymes if a particular element of structural symmetry was disrupted. For the BamHI MTase, integrity of the ATCC was critical, while an intact GAT sequence was necessary for the activity of T4Dam, and an intact GA was necessary for EcoDam. Theoretical alignment of the region of best contacts between the protein and DNA showed that in the case of a palindromic interaction site, a zone covering the 5′-symmetric residues is located in the major groove versus a zone of contact covering the 3′-symmetric residues in the minor groove. Our data fit a simple rule of thumb that the most important contacts are aligned around the methylation target base: if the target base is in the 5′ half of the palindrome, the interaction between the enzyme and the DNA occurs mainly in the major groove; if it is in the 3′ half, the interaction occurs mainly in the minor groove.     

Two mutations of basic residues within the N-terminus of HMG-1 B domain with different effects on DNA supercoiling and binding to bent DNA

High mobility group (HMG) 1 protein and its two homologous DNA-binding domains, A and B ("HMG-boxes"), can bend and supercoil DNA in the presence of topoisomerase I, as well as recognize differently bent and distorted DNA structures, including four-way DNA junctions, supercoiled DNA and DNA modified with anticancer drug cisplatin. Here we show that the lysine-rich part of the linker region between A and B domains of HMG-1, the (85)TKKKFKD(91) sequence that is attached to the N-terminus of the B domain within HMG-1, is a prerequisite for a preferential binding of the B domain to supercoiled DNA. The above sequence is also essential for a high-affinity binding of the B domain to DNA containing a site-specific major 1,2-d(GpG) intrastrand DNA adduct of cisplatin. Mutation of Arg(97), but not Lys(90) [Lys(90) forms a specific cross-link with platinum(II) in major groove of cisplatin-modified DNA; Kane, S. A., and Lippard, S. J. (1996) Biochemistry 35, 2180--2188], to alanine significantly (>40-fold) reduces affinity of the B domain to cisplatin-modified DNA, inhibits the ability of the B domain to bend (ligase-mediated circularization) or supercoil DNA, and results in a loss of the preferential binding of the B domain to supercoiled DNA without affecting the structural-specificity of the HMG-box for four-way DNA junctions. Some of the reported activities of the B domain are enhanced when the B domain is covalently linked to the A domain. We propose that binding of the A/B linker region within the major DNA groove helps the two HMG-1 domains to anchor to the minor DNA groove to facilitate their DNA binding and other activities.     

Theoretical exploration of netropsin binding to tRNA(Phe)

Theoretical exploration of the possible interaction of netropsin with tRNAPhe indicates that binding should occur preferentially with the major groove of the T psi C stem of the macromolecule, specifically with the bases G51, U52, G53 and phosphates 52, 53, 61 and 62. This agrees with the recent crystallographic result of Rubin and Sundaralingam. It is demonstrated that the difference with respect to netropsin binding with B-DNA, where it occurs specifically in the minor groove of AT sequences, is due to the differences in the distribution of the electrostatic molecular potential generated by these different types of DNA: this potential is sequence dependent in B-DNA (located in the minor groove of AT sequences and the major groove of GC sequences), while it is sequence independent and always located in the major groove in A-RNA. The result demonstrates the major role of electrostatics in determining the location of the binding site.     

Stacking interaction of guanine with netropsin in the minor groove of d(CGTATATACG)2

We present the structure of the decanucleotide d(CGTATATACG) determined by single crystal X-ray diffraction at 1.58 A resolution. A netropsin drug is found in the minor groove with guanine stacked on a pyrrole ring of the drug, a feature described here for the first time. The stacked guanine is an extra-helical base coming from the end of a neighbour oligonucleotide. This observation may open the way to the development of minor groove binding drugs with a higher sequence selectivity. The oligonucleotide is in the B-conformation, but the terminal base-pairs are disrupted: the cytosine residues are disordered while the guanine residues penetrate into the minor groove of neighbouring duplexes. Four hydrated Ni ions with octahedral co-ordination are found associated with the N7 atoms of each guanine. The high affinity of these ions with guanine suggests that they may be used as probes for specific guanine residues.     

Hx-amides: DNA sequence recognition by the fluorescent Hx (p-anisylbenzimidazole)•pyrrole and Hx•imidazole pairings

Hx-amides are fluorescent hybrids of imidazole (I)- and pyrrole (P)-containing polyamides and Hoechst 33258, and they bind in the minor groove of specific DNA sequences. Synthesis and DNA binding studies of HxII (5) complete our studies on the first set of Hx-amides: Hx–I/P–I/P. HxPP (2), HxIP (3) and HxPI (4) were reported earlier. Results from DNase I footprinting, biosensor-SPR, CD and ΔT_M studies showed that Hx-amides interacted with DNA via the anti-parallel and stacked, side-by-side motif. Hx was found to mimic the DNA recognition properties of two consecutive pyrrole units (PP) in polyamides. Accordingly, the stacked Hx/PP pairing binds preferentially to two consecutive AT base pairs, A/T–A/T; Hx/IP prefers C–A/T; Hx/PI prefers A/T–C; and Hx/II prefers C–C. The results also showed that Hx-amides bound their cognate sequence at a higher affinity than their formamido-triamide counterparts.