Head-to-head bis-hairpin polyamide minor groove binders and their conjugates with triplex-forming oligonucleotides: studies of interaction with target double-stranded DNA

Two hairpin hexa(N-methylpyrrole)carboxamide DNA minor groove binders (MGB) were linked together via their N-termini in head-to-head orientation. Complex formation between these bis-MGB conjugates and target DNA has been studied using DNase I footprinting, circular dichroism, thermal dissociation, and molecular modeling. DNase I footprint revealed binding of these conjugates to all the sites of 492 b.p. DNA fragment containing (A/T)(n)X(m)(A/T)(p) sequences, where n>3, p>3; m=1,2; X = A,T,G, or C. Binding affinity depended on the sequence context of the target. CD experiments and molecular modeling showed that oligo(N-methylpyrrole)carboxamide moieties in the complex form two short antiparallel hairpins rather than a long parallel head-to-head hairpin. Binding of bis-MGB also stabilized a target duplex thermodynamically. Sequence specificity of bis-MGB/DNA binding was validated using bis-conjugates of sequence-specific hairpin (N-methylpyrrole)/(N-methylimidazole) carboxamides. In order to increase the size of recognition sequence, the conjugates of bis-MGB with triplex-forming oligonucleotides (TFO) were synthesized and compared to TFO conjugated with single MGB hairpin unit. Bis-MGB-oligonucleotide conjugates also bind to two blocks of three and more A.T/T.A pairs similarly to bis-MGB alone, independently of the oligonucleotide moiety, but with lower affinity. However, the role of TFO in DNA recognition was demonstrated for mono-MGB-TFO conjugate where the binding was detected mainly in the area of the target sequence consisting of both MGB and TFO recognition sites. Basing on the molecular modeling, three-dimensional models of both target DNA/bis-MGB and target DNA/TFO-bis-MGB complexes were built, where bis-MGB forms two antiparallel hairpins. According to the second model, one MGB hairpin is in the minor groove of 5'-adjacent A/T sequence next to the triplex-forming region, whereas the other one occupies the minor groove of the TFO binding polypurine tract. All these data together give a key information for the construction of MGB-MGB and MGB-oligonucleotide conjugates possessing high specificity and affinity for the target double-stranded DNA.     

Effects of compound structure on carbazole dication-DNA complexes: tests of the minor-groove complex models

Carbazole dications have shown excellent activity against opportunistic infections, but they are quite different in structure from previously studied unfused aromatic cations that function by targeting the DNA minor groove. In a previous report [Tanious, F. A., Ding, D., Patrick, D. A., Tidwell, R. R., and Wilson, W. D. (1997) Biochemistry 36, 15315-15325] we showed that, despite their fused ring structure, the carbazoles also bind in A/T sequences of the DNA minor groove and we proposed models for the carbazole-DNA complexes with the carbazole nitrogen facing out of the groove for 3,6 substituted compounds but into the groove in 2,7 carbazoles. To test and refine the models, carbazole-N-methyl substituted derivatives have been synthesized in both the 3,6 and 2,7 series as well as a new 2,6 substituted NH derivative that is intermediate in structure. Footprinting results indicate a broad AT specificity of carbazole binding and a pattern in agreement with a minor groove complex. Surface plasmon resonance biosensor analysis of carbazole binding to an oligomer with an AATT central sequence indicated that the 2,7 NH compound has the largest binding constant. Both the 3,6 NH and NMe compounds bind with similar equilibrium constants that are less than for the 2,7 NH compound. The 2,7 NMe compound has the lowest binding constant of all the carbazoles. Spectroscopic results are also similar for the two 3,6 derivatives but are quite different for the 2,7 NH and NMe carbazole dications. Structural analysis of carbazole complexes with an AATT sequence by 2D NMR methods also supported a minor groove complex of the carbazoles in orientations in agreement with the previously proposed models. From these results, it is clear that the fused ring carbazoles can bind strongly in the DNA minor groove with a broad A/T specificity and that the 2,7 and 3,6 substituted carbazoles bind to the minor groove in opposite orientations.     

Sequence-specific interaction of the Salmonella Hin recombinase in both major and minor grooves of DNA.

The Hin recombinase of Salmonella catalyzes a site-specific recombination event which leads to flagellar phase variation. Starting with a fully symmetrical recombination site, hixC, a set of 40 recombination sites which vary by pairs of single base substitutions was constructed. This set was incorporated into the Salmonella-specific bacteriophage P22 based challenge phage selection and used to define the DNA sequence determinants for the binding of Hin to DNA in vivo. The critical sequence-specific contacts between a Hin monomer and a 13 bp hix half-site are at two T:A base pairs in the major groove of the DNA which are separated by one base pair, and two consecutive A:T contacts in the minor groove. The base substitutions in the major groove recognition portion which were defective in binding Hin still retained residual binding capability in vivo, while the base pair substitutions affecting the minor groove recognition region lost all in vivo binding. Using in vitro binding assays, Hin was found to bind to hix symmetrical sites with A:T base pairs or I:C base pairs in the minor groove recognition sequences, but not to G:C base pairs. In separate in vitro binding assays, Hin was equally defective in binding to either a G:C or a I:C contact in a major groove recognition sequence. Results from in vitro binding assays to hix sites in which 3-deazaadenine was substituted for adenine are consistent with Hin making a specific contact to either the N3 of adenine or O2 of thymine in the minor groove within the hix recombination site on each symmetric half-site. These results taken with the results of previous studies on the DNA binding domain of Hin suggest a sequence-specific minor groove DNA binding motif.     

Structural effects of DNA sequence on T.A recognition by hydroxypyrrole/pyrrole pairs in the minor groove

Synthetic polyamides composed of three types of aromatic amino acids, N-methylimidazole (Im), N-methylpyrrole (Py) and N-methyl-3-hydroxypyrrole (Hp) bind specific DNA sequences as antiparallel dimers in the minor groove. The side-by-side pairings of aromatic rings in the dimer afford a general recognition code that allows all four base-pairs to be distinguished. To examine the structural consequences of changing the DNA sequence context on T.A recognition by Hp/Py pairs in the minor groove, crystal structures of polyamide dimers (ImPyHpPy)(2) and the pyrrole counterpart (ImPyPyPy)(2) bound to the six base-pair target site 5'-AGATCT-3' in a ten base-pair oligonucleotide have been determined to a resolution of 2.27 and 2.15 A, respectively. The structures demonstrate that the principles of Hp/Py recognition of T.A are consistent between different sequence contexts. However, a general structural explanation for the non-additive reduction in binding affinity due to introduction of the hydroxyl group is less clear. Comparison with other polyamide-DNA cocrystal structures reveals structural themes and differences that may relate to sequence preference.     

Enhanced uranyl photocleavage across the minor groove of all (A/T)4 sequences indicates a similar narrow minor groove conformation

The uranyl(VI)-mediated photocleavage of a Drew–Dickerson sequence oligonucleotide (5′-dGATCACGCGAATTCGCGT) either as the (self-complementary) duplex or cloned into the BamH1 site of pUC19 has been studied. At pH 6.5 in acetate buffer relatively enhanced photocleavage is observed at the 3′-end of the AATT sequence corresponding to maximum cleavage across the minor groove in the A/T tract. Thus maximum cleavage correlates with minimum minor groove width in the crystal structure and also with the largest electronegative potential according to computations. Using plasmid constructs with cloned inserts of the type [CGCG(A/T₄)]_n, we also analysed all possible sequence combinations of the (A/T)₄ tract and in all cases we observed maximum uranyl-mediated photocleavage across the minor groove in the (A/T)₄ tract without any significant differences between the various sequences. From these results we infer that DNA double helices of all (A/T)₄ sequences share the same narrow minor groove helix conformation.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

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.     

Interchange of DNA-binding modes in the deformed and ultrabithorax homeodomains: a structural role for the N-terminal arm

The deformed (Dfd) and ultrabithorax (Ubx) homeoproteins regulate developmental gene expression in Drosophila melanogaster by binding to specific DNA sequences within its genome. DNA binding is largely accomplished via a highly conserved helix-turn-helix DNA-binding domain that is known as a homeodomain (HD). Despite nearly identical DNA recognition helices and similar target DNA sequence preferences, the in vivo functions of the two proteins are quite different. We have previously revealed differences between the two HDs in their interactions with DNA. In an effort to define the individual roles of the HD N-terminal arm and recognition helix in sequence-specific binding, we have characterized the structural details of two Dfd/Ubx chimeric HDs in complex with both the Dfd and Ubx-optimal-binding site sequences. We utilized hydroxyl radical cleavage of DNA to assess the positioning of the proteins on the binding sites. The effects of missing nucleosides and purine methylation on HD binding were also analyzed. Our results show that both the Dfd and Ubx HDs have similar DNA-binding modes when in complex with the Ubx-optimal site. There are subtle but reproducible differences in these modes that are completely interchanged when the Dfd N-terminal arm is replaced with the corresponding region of the Ubx HD. In contrast, we showed previously that the Dfd-optimal site sequence elicits a very different binding mode for the Ubx HD, while the Dfd HD maintains a mode similar to that elicited by the Ubx-optimal site. Our current methylation interference studies suggest that this alternate binding mode involves interaction of the Ubx N-terminal arm with the minor groove on the opposite face of DNA relative to the major groove that is occupied by the recognition helix. As judged by hydroxyl radical footprinting and the missing nucleoside experiment, it appears that interaction of the Ubx recognition helix with the DNA major groove is reduced. Replacing the Dfd N-terminal arm with that of Ubx does not elicit a complete interchange of the DNA-binding mode. Although the position of the chimera relative to DNA, as judged by hydroxyl radical footprinting, is similar to that of the Dfd HD, the missing nucleoside and methylation interference patterns resemble those of the Ubx HD. Repositioning of amino acid side-chains without wholesale structural alteration in the polypeptide appears to occur as a function of N-terminal arm identity and DNA-binding site sequence. Complete interchange of binding modes was achieved only by replacement of the Dfd N-terminal arm and the recognition helix plus 13 carboxyl-terminal residues with the corresponding residues of Ubx. The position of the N-terminal arm in the DNA minor groove appears to differ in a manner that depends on the two base-pair differences between the Dfd and Ubx-optimal-binding sites. Thus, N-terminal arm position dictates the binding mode and the interaction of the recognition helix with nucleosides in the major groove.     

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.     

Syntheses and DNA-binding studies of a series of unsymmetrical cyanine dyes: structural influence on the degree of minor groove binding to natural DNA

Twelve crescent-shaped unsymmetrical dyes have been synthesized and their interactions with DNA have been investigated by spectroscopic methods. A new facile synthetic route to this type of cyanine dyes has been developed, involving the preparation of 6-substituted 2-thiomethyl-benzothiazoles in good yields. The new dyes are analogues to the minor groove binding unsymmetrical cyanine dye, BEBO, recently reported by us. In this dye, the structure of the known intercalating cyanine dye BO was extended with a 6-methylbenzothiazole substituent. Herein we further investigate the role of the extending benzazole heterocycle, as well as of the pyridine or quinoline moiety of the cyanine chromophore, for the binding mode of these crescent-shaped dyes to calf thymus DNA. Flow LD and CD studies of the 12 dyes show that the extent of minor groove binding to mixed sequence DNA varies significantly between the dyes. We find that hydrophobicity and size are the crucial parameters for recognition of the minor groove. The relatively high fluorescence quantum yield of many of these cyanines bound to DNA, combined with their absorption at long wavelengths, may render them useful in biological applications. In particular, two of the benzoxazole containing dyes BOXTO and 2-BOXTO show a high degree of minor groove binding and quantum yields of 0.52 and 0.32, respectively, when bound to DNA.     

Extra-helical guanine interactions in DNA

We review the extra-helical guanine interactions present in many oligonucleotide crystals. Very often terminal guanines interact with other guanines in the minor groove of neighboring oligonucleotides through N2 x N3 hydrogen bonds. In other cases the interaction occurs with the help of Ni2+ ions. Guanine/netropsin stacking in the minor groove has also been found. From these studies we conclude that guanine may have multiple extra-helical interactions. In particular it may be considered a very effective minor groove binder, which could be used in the design of sequence selective binding drugs. Interactions through the major groove are seldom encountered, but might be present when DNA is stretched. Such interactions are also analyzed, since they might be important for homologous chromosome pairing during meiosis.     

DNA curvature in native and modified EcoRI recognition sites and possible influence upon the endonuclease cleavage reaction

The ligation of a decadeoxynucleotide containing the EcoRI recognition site forms a series of multimers which appear to be curved based on observed anomalous gel migration in polyacrylamide gels. The degree of DNA curvature present in the recognition sequence, based upon the observed migration anomaly, can be altered by modifications to the purine functional groups at the 2- and 6-positions. Deletion of the guanine 2-amino group, occurring in the minor groove of the B-DNA helix, is most effective in increasing the observed DNA curvature. Conversely, the displacement of an amino group from the major groove to the minor groove eliminates curvature. DNA curvature is also modulated by the exocyclic group at the purine 6-position with decreasing curvature observed when changing the amino group to a carbonyl or proton substituent. Differences in the kinetic parameters characterizing the cleavage reaction by the endonuclease for many of the modified sequences are the result of modifications of functional groups in the major groove, which are likely to contact the endonuclease during catalysis. However, with two examples, significant decreases in the observed specificity constant (kcat/Km), characterizing the protein-nucleic acid interaction, cannot be easily explained in terms of such functional group contacts. It is more likely in these cases that the functional group modifications affect the efficiency of the endonuclease-DNA interaction by modulation of the structure of the double-stranded DNA helix. With both examples, modifications have been made to minor groove substituents. The extent of DNA curvature is increased significantly for one and decreased for the other, compared with that observed for the native recognition site. The results suggest that curvature of the DNA helix axis is an intrinsic property of the d(GAATTC) sequence which helps to optimize the protein-nucleic acid interactions observed for the EcoRI restriction endonuclease.     

Selective binding to AT sequences in DNA by an acridine-linked peptide containing the SPKK motif.

The sequence selectivity of binding to DNA by an acridine-linked peptide ligand has been investigated by means of footprinting methodologies. The ligand conjugates an anilino-acridine intercalating chromophore with the potentially minor groove binder octapeptide SPKKSPKK. This basic peptide corresponds to a highly conserved DNA recognition motif found in histone H1 and several other nonhistone proteins. Three complementary techniques using DNase I, hydroxyl radicals and osmium tetroxide as sequencing probes have been employed to evaluate both the sequence specificity of binding and the drug-induced conformational changes in DNA. The results converge to demonstrate the AT-selectivity and support a model in which the peptide moiety lies in the minor groove. DNA-binding sites of the conjugate are restricted to a few alternating AT-sequences proximal to GC-rich regions. Binding to homooligomeric runs of A and T is clearly disfavoured by the hybrid whereas such sequences represent preferred binding sites for the unsubstituted basic peptide. These differences reflect the influence of the anilino-acridine chromophore, which evidently contributes to the DNA recognition process allowing the peptide only to contact defined DNA sequences.     

Dissecting direct and indirect readout of cAMP receptor protein DNA binding using an inosine and 2,6-diaminopurine in vitro selection system

The DNA interaction of the Escherichia coli cyclic AMP receptor protein (CRP) represents a typical example of a dual recognition mechanism exhibiting both direct and indirect readout. We have dissected the direct and indirect components of DNA recognition by CRP employing in vitro selection of a random library of DNA-binding sites containing inosine (I) and 2,6-diaminopurine (D) instead of guanine and adenine, respectively. Accordingly, the DNA helix minor groove is structurally altered due to the ‘transfer’ of the 2-amino group of guanine (now I) to adenine (now D), whereas the major groove is functionally intact. The majority of the selected sites contain the natural consensus sequence TGTGAN₆TCACA (i.e. TITIDN₆TCDCD). Thus, direct readout of the consensus sequence is independent of minor groove conformation. Consequently, the indirect readout known to occur in the TG/CA base pair step (primary kink site) in the consensus sequence is not affected by I–D substitutions. In contrast, the flanking regions are selected as I/C rich sequences (mostly I-tracts) instead of A/T rich sequences which are known to strongly increase CRP binding, thereby demonstrating almost exclusive indirect readout of helix structure/flexibility in this region through (anisotropic) flexibility of I-tracts.     

Combinatorial determination of sequence specificity for nanomolar DNA-binding hairpin polyamides

Development of sequence-specific DNA-binding drugs is an important pharmacological goal, given the fact that numerous existing DNA-directed chemotherapeutic drugs rely on the strength and selectivity of their DNA interactions for therapeutic activity. Among the DNA-binding antibiotics, hairpin polyamides represent the only class of small molecules that can practically bind any predetermined DNA sequence. DNA recognition by these ligands depends on their side-by-side amino acid pairings in the DNA minor groove. Extensive studies have revealed that these molecules show extremely high affinity for sequence-directed, minor groove interaction. However, the specificity of such interactions in the presence of a large selection of sequences such as the human genome is not known. We used the combinatorial selection method restriction endonuclease protection, selection, and amplification (REPSA) to determine the DNA binding specificity of two hairpin polyamides, ImPyPyPy-gamma-PyPyPyPy-beta-Dp and ImPyPyPy-gamma-ImPyPyPy-beta-Dp, in the presence of more than 134 million different sequences. These were verified by restriction endonuclease protection assays and DNase I footprinting analysis. Our data showed that both hairpin polyamides preferentially selected DNA sequences having consensus recognition sites as defined by the Dervan pairing rules. These consensus sequences were rather degenerate, as expected, given that the stacked pyrrole-pyrrole amino acid pairs present in both polyamides are unable to discriminate between A.T and T.A base pairs. However, no individual sequence within these degenerate consensus sequences was preferentially selected by REPSA, indicating that these hairpin polyamides are truly consensus-specific DNA-binding ligands. We also discovered a preference for overlapping consensus binding sites among the sequences selected by the hairpin polyamide ImPyPyPy-gamma-PyPyPyPy-beta-Dp, and confirmed by DNase I footprinting that these complex sites provide higher binding affinity. These data suggest that multiple hairpin polyamides can cooperatively bind to their highest-affinity sites.