Drug Recognition of DNA. Proposal for GC Minor Groove Specific Ligands: Vinylexins

Abstract

In a previous publication in this journal (ref. 14) we have proposed an isolexin-like prototype of a GC minor groove specific ligand. The present paper is devoted to refinements of this prototype (increase in specificity and in DNA binding energy). It is shown that only a very limited improvement can be obtained by increasing the proton accepting capabilities of the heteroaromatic ring systems of the prototype, although these rings interact directly with the proton donating NH₂ group of guanine. On the other hand a significant increase both in GC specificity and in DNAbinding energy is obtained by replacing the NH linkers of the isolexin by C=C double bonds (yielding what we term “vinylexins”). Specificity is still largely conserved and the DNA binding energy is significantly increased in monocationic vinylexins, which should thus be efficient GC minor groove specific ligands. The outstanding importance for the GC specificity of the C=C linkers is evidenced by the disappearance of this specificity when these linkers are replaced by peptide bonds (peptilexins). On the other hand vinylexins with proton donating heteroaromatic rings are, as expected, AT specific. The vinylexin family may thus represent universal minor groove binding agents susceptible to bind to any given base pair sequence of DNA following the positioning of their proton donor and proton acceptor rings. This study confirms the insufficiency of purely geometrical and/or hydrogen bonding considerations for the correct estimation of GC versus AT specificity of groove binding ligands. These can only be accounted for by taking into consideration the overall electronic properties of the interacting species and explicitly calculating the energies of complex formation including all the relevant contributions.     

Drug recognition of DNA. Proposal for GC minor groove specific ligands: vinylexins

In a previous publication in this journal we have proposed an isolexin-like prototype of a GC minor groove specific ligand. The present paper is devoted to refinements of this prototype (increase in specificity and in DNA binding energy). It is shown that only a very limited improvement can be obtained by increasing the proton accepting capabilities of the heteroaromatic ring systems of the prototype, although these rings interact directly with the proton donating NH2 group of guanine. On the other hand a significant increase both in GC specificity and in DNA binding energy is obtained by replacing the NH linkers of the isolexin by C = C double bonds (yielding what we term "vinylexins"). Specificity is still largely conserved and the DNA binding energy is significantly increased in monocationic vinylexins, which should thus be efficient GC minor groove specific ligands. The outstanding importance for the GC specificity of the C = C linkers is evidenced by the disappearance of this specificity when these linkers are replaced by peptide bonds (peptilexins). On the other hand vinylexins with proton donating heteroaromatic rings are, as expected, AT specific. The vinylexin family may thus represent universal minor groove binding agents susceptible to bind to any given base pair sequence of DNA, following the positioning of their proton donor and proton acceptor rings. This study confirms the insufficiency of purely geometrical and/or hydrogen bonding considerations for the correct estimation of GC versus AT specificity of groove binding ligands. These can only be accounted for by taking into consideration the overall electronic properties of the interacting species and explicitly calculating the energies of complex formation including all the relevant contributions.     

Groove Width and Depth of B-DNA Structures Depend on Local Variation in Slide

Abstract

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

Groove width and depth of B-DNA structures depend on local variation in slide.

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA. A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

DNA-sequence specific recognition by a thiazole analogue of netropsin: a comparative footprinting study.

Four different footprinting techniques have been used to probe the DNA sequence selectivity of Thia-Net, a bis-cationic analogue of the minor groove binder netropsin in which the N-methylpyrrole moieties are replaced by thiazole groups. In Thia-Net the ring nitrogen atoms are directed into the minor groove where they could accept hydrogen bonds from the exocyclic 2-amino group of guanine. Three nucleases (DNAase I, DNAase II, and micrococcal nuclease) were employed to detect binding sites on the 160bp tyr T fragment obtained from plasmid pKM delta-98, and further experiments were performed with 117mer and 253mer fragments cut out of the plasmid pBS. MPE.Fe(II) was used to footprint binding sites on an EcoRI/HindIII fragment from pBR322. Thia-Net binds to sites in the minor groove containing 4 or 5 base pairs which are predominantly composed of alternating A and T residues, but with significant acceptance of intrusive GC base pairs. Unlike the parent antibiotic netropsin, Thia-Net discriminates against homooligomeric runs of A and T. The evident preference of Thia-Net for AT-rich sites, despite its containing thiazole nitrogens capable of accepting GC sites by hydrogen bonding, supports the view that the biscationic nature of the ligand imposes a bias due to the electrostatic potential differences in the receptor which favour the ligand reading alternating AT sequences.     

Comparative thermodynamics for monomer and dimer sequence-dependent binding of a heterocyclic dication in the DNA minor groove

Phenylamidine cationic groups linked by a furan ring (furamidine) and related symmetric diamidine compounds bind as monomers in the minor groove of AT sequences of DNA. DB293, an unsymmetric derivative with one of the phenyl rings of furamidine replaced with a benzimidazole, can bind to AT sequences as a monomer but binds more strongly to GC-containing minor-groove DNA sites as a stacked dimer. The dimer-binding mode has high affinity, is highly cooperative and sequence selective. In order to develop a better understanding of the correlation between structural and thermodynamic aspects of DNA molecular recognition, DB293 was used as a model to compare the binding of minor-groove agents with AT and mixed sequence DNA sites. Isothermal titration calorimetry and surface plasmon resonance results clearly show that the binding of DB293 and other related compounds into the minor groove of AT sequences is largely entropy-driven while the binding of DB293 as a dimer into the minor groove of GC-containing sequences is largely enthalpy-driven. At 25 degrees C, for example, the AT binding has DeltaG degrees, DeltaH degrees and TDeltaS degrees values of -9.6, -3.6 and 6.0 kcal/mol while the values for dimer binding to a GC-containing site are -9.0, -10.9 and -1.9 kcal/mol (per mol of bound compound), respectively. These results show that the thermodynamic components for binding of compounds of this type to DNA are very dependent on the structure, solvation and sequence of the DNA binding site.     

Specificity in the interaction of non intercalative groove binding ligands with nucleic acids

This paper presents results of theoretical computations on the interaction energies and geometries for the binding to nucleic acids of a number of representative groove binding non intercalating drugs: netropsin, distamycin A, SN 18071, etc. The computations account for the specificity of binding in all cases and demonstrate that the formation of hydrogen bonds is not necessary neither for binding nor for the preference for the minor groove of AT sequences of B-DNA. It appears that if a relatively good steric fit can be obtained in the minor groove, the interaction will be preferentially stabilized there by the favorable electrostatic potential generated in this groove by the AT sequences. The computation of the interaction energies in free space does not reproduce, however, the order of affinities of the ligands studied and yields too great values of the binding energies. The introduction of the solvent effect, through the computation of the hydration and cavitation effects, confirms the specificity, improves the ordering and brings the values of the energies close to the experimental ones. The theoretical account of the “surprising” effect of netrospin binding to the major groove of theTψC stem of tRNA^Phe confirms the decisive significance of the distribution of the molecular electrostatic potential for the selection of the binding site. The inclusion in the computations of the flexibility of DNA enables to predict correctly the main features of the macromolecular deformation upon the binding of the ligand.     

DNA binding of the nonintercalative ligands SN-6132, SN-6131 and SN-6113: minor variations of the ligand structure may cause changes in the base pair preference.

The DNA binding selectivity of three ligands of a series of antitumor agents of bisquaternary ammonium heterocycles has been investigated by means of CD spectroscopy and melting measurements. From the spectroscopic results and binding data it is concluded that the agents SN-6132, SN-6131 and SN-6113 have relatively high affinity to AT base pair sequences whereas the binding to GC pairs is very low. The binding selectivity to AT base pair sequences decreases in the order netropsin > SN-6132 > SN-6113 > SN-6131. Poly(dA).poly(dT) has the highest binding preference for SN-6132 relative to that of SN-6131. The different binding behavior of the ligands is related to their distinct changes in the chemical structure and to the DNA minor groove properties which determines the adaptability of the ligands in the groove.     

A polarized-light spectroscopy study of interactions of a hairpin polyamide with DNA

We here study the interactions of a polyamide with large DNA, and compare to those of minor groove binder distamycin (DST), including high ligand/DNA binding ratios. Specific as well as nonspecific binding is probed using polarized-light spectroscopy combined with singular value decomposition analysis. Circular and linear dichroism data confirm binding geometries consistent with minor groove binding for both of the ligands. Interestingly, at high and intermediate ligand/DNA ratios the polyamide exhibits no significant sequence discrimination between mixed-sequence (calf thymus) and AT DNA as compared to DST. Each ligand is concluded to exhibit two different binding modes depending upon ligand/DNA ratio and nucleo-base sequence. At high binding ratios, distinct differences between the ligands are observed: circular dichroism spectra exciton effects provide evidence of bimolecular interactions of the polyamide when bound to AT-DNA, whereas no effects are seen with DST or mixed-sequence DNA. Also linear dichroism indicates that a change in binding geometry occurs at high polyamide/AT ratios, and that the effect occurs only with polyamide in contrast to DST. Since the effect is insignificant with DST, or with calf thymus DNA, it is concluded that it relates to the sizes of the ligands and the minor grooves, becoming critical in the limit of crowding.     

Design of Sequence-Specific DNA Binding Ligands That Use a Two-Stranded Peptide Motif for DNA Sequence Recognition

Abstract

The design and DNA binding activity of β-structure-forming peptides and netropsin-peptide conjugates are reported. It is found that a pair of peptides - S,S'-bis(Lys-Gly-Val-Cys-Val- NH-NH-Dns) - bridged by an S-S bond binds at least 10 times more strongly to poly(dG)?poly(dC) than to poly(dA)?poly(dT). This peptide can also discriminate between 5′-GpG-3′ and 5′-GpC-3′ steps in the DNA minor groove. Based on these observations, new synthetic ligands, bis-netropsins, were constructed in which two netropsin-like fragments were attached by means of short linkers to a pair of peptides - Gly-Cys-Gly- or Val-Cys-Val - bridged by S-S bonds. These compounds possess a composite binding specificity: the peptide chains recognize 5′-GpG-3′ steps on DNA, whereas the netropsin-like fragments bind preferentially to tuns of 4 AT base pairs. Our data indicate that combining the AT-base-pair specific properties of the netropsin-type structure with the 5′-GpG-3′-specific properties of certain oligopeptides offers a new approach to the synthesis of ligands capable of recognizing mixed sequences of AT- and GC-base pairs in the DNA minor groove. These compounds are potential models for DNA-binding domains in proteins which specifically recognize base pair sequences in the minor groove of DNA.     

Quantitative and sequence-specific analysis of DNA-ligand interaction by means of fluorescent intercalator probes

A novel method of analysis of double-stranded DNA-ligand interaction is presented. The interaction is monitored by the fluorescence of a DNA bis-intercalator oxazole homodimer YoYo-3. The fluorescence intensity or its decay time reflects the modification of the DNA double helix. The DNA sequence is scanned by hybridization with short oligomers having consecutively overlapping complementary sequences to analyse the sequence specificity of binding. In our experiments we used as ligands the minor groove binders netropsin, SN6999 (both with AT-preference), the GC-specific ligand chromomycin A3 as well as the derivative SN6113 (non-specific interaction), which displace the bis-intercalator YoYo-3 or influence the duplex structure in such away that the fluorescence intensity and lifetime decrease in comparison to a ligand-free screening. The changes of fluorescence emission clearly define the binding motif and indicate minor groove interactions with a reduced DNA binding site. Titration of the ligand quantitatively characterizes its binding by determining the dependence of the binding constant on the oligonucleotide sequence.     

Structural analysis of the binding modes of minor groove ligands comprised of disubstituted benzenes

Two-dimensional homonuclear NMR was used to characterize synthetic DNA minor groove-binding ligands in complexes with oligonucleotides containing three different A-T binding sites. The three ligands studied have a C₂ axis of symmetry and have the same general structural motif of a central para-substituted benzene ring flanked by two meta-substituted rings, giving the molecules a crescent shape. As with other ligands of this shape, specificity seems to arise from a tight fit in the narrow minor groove of the preferred A-T-rich sequences. We found that these ligands slide between binding subsites, behavior attributed to the fact that all of the amide protons in the ligand backbone cannot hydrogen bond to the minor groove simultaneously.     

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

Abstract

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5′-YpR-3′ intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Binding of symmetrical cyanine dyes into the DNA minor groove

Optical methods, such as fluorescence, circular dichroism and linear flow dichroism, were used to study the binding to DNA of four symmetrical cyanine dyes, each consisting of two identical quinoline, benzthiazole, indole, or benzoxazole fragments connected by a trimethine bridge. The ligands were shown to form a monomer type complex into the DNA minor groove. The complex of quinoline-containing ligand with calf thymus DNA appeared to be the most resistant to ionic strength, and it did not dissociate completely even in 1 M NaCl. Binding of cyanine dyes to DNA could also be characterized by possibility to form ligand dimers into the DNA minor groove, by slight preference of binding to AT pairs, as well as by possible intercalation between base pairs of poly(dG)-poly(dC). The correlation found between the binding constants to DNA and the extent of cyanine dyes hydrophobicity estimated as the n-octanol/water partition coefficient is indicative of a significant role of hydrophobic interactions for the ligand binding into the DNA minor groove.     

Alternative heterocycles for DNA recognition: a 3-pyrazole/pyrrole pair specifies for G.C base pairs

Synthetic ligands comprising three aromatic amino acids, pyrrole (Py), imidazole (Im), and hydroxypyrrole (Hp), specifically recognize predetermined sequences as side-by-side pairs in the minor groove of DNA. To expand the repertoire of aromatic rings that may be utilized for minor groove recognition, three five-membered heterocyclic rings, 3-pyrazolecarboxylic acid (3-Pz), 4-pyrazolecarboxylic acid (4-Pz), and furan-2-carboxylic acid (Fr), were examined at the N-terminus of eight-ring hairpin polyamide ligands. The DNA binding properties of 3-Pz, 4-Pz, and Fr each paired with Py were studied by quantitative DNase I footprinting titrations on a 283 bp DNA restriction fragment containing four 6-bp binding sites 5'-ATNCCTAA-3' (N = G, C, A, or T; 6-bp polyamide binding site is underlined). The pair 3-Pz/Py has increased binding affinity and sequence specificity for G.C bp compared with Im/Py.     

Modelling basic features of specificity in the binding of a dicationic steroid diamine to double-stranded oligonucleotides.

An investigation of the intrinsically preferred binding modes of a steroid diamine, dipyrandium, to the double-stranded hexanucleotides d(TATATA)2, d(ATATAT)2, and d(CGCGCG)2 is carried out by the energy minimization procedure JUMNA. Several alternative binding modes are compared: groove binding in which the conformation of the oligonucleotide remains close to that of B-DNA, intercalation between base-pairs and interaction with variously kinked structures in which base pairs of dinucleoside steps open towards the groove in which the binding occurs. The favored binding configuration occurs at the d(TpA) step of the AT kinked nucleotides in which the kink opens the base pairs towards the minor groove. Thus, for the d(T1A2T3A4T5A6)2 sequences the preferred complexation involves the kink at the T3A4 step facing the cyclohexane rings A, B, and C of the ligand. For the d(A1T2A3T4A5T6)2 sequence, the kink occurs at the T2A3 step facing the cationic pyrrolidine ring linked to ring A. The binding of dipyrandium to d(CGCGCG)2 is found to be considerably less favourable than for either of the two (AT) sequences.     

DNA sequence preferences of several AT-selective minor groove binding ligands.

We have examined the interaction of distamycin, netropsin, Hoechst 33258 and berenil, which are AT-selective minor groove-binding ligands, with synthetic DNA fragments containing different arrangements of AT base pairs by DNase I footprinting. For fragments which contain multiple blocks of (A/T)4 quantitative DNase I footprinting reveals that AATT and AAAA are much better binding sites than TTAA and TATA. Hoechst 33258 shows that greatest discrimination between these sites with a 50-fold difference in affinity between AATT and TATA. Alone amongst these ligands, Hoechst 33258 binds to AATT better than AAAA. These differences in binding to the various AT-tracts are interpreted in terms of variations in DNA minor groove width and suggest that TpA steps within an AT-tract decrease the affinity of these ligands. The behaviour of each site also depends on the flanking sequences; adjacent pyrimidine-purine steps cause a decrease in affinity. The precise ranking order for the various binding sites is not the same for each ligand.