Interaction of minor-groove-binding diamidine ligands with an asymmetric DNA duplex. NMR and molecular modelling studies

First Structure of a Designed Minor Groove Binding Heterocyclic Cation that Specifically Recognizes Mixed DNA Base Pair Sequences

The high-resolution NMR structure of the first heterocyclic, non-amide, organic cation that strongly and selectively recognizes mixed AT/GC bp (bp=base pair) sequences of DNA in a 1:1 complex is described. Compound designs of this type provide essential methods for control of functional, non-genomic DNA sequences and have broad cell uptake capability, based on studies from animals to humans. The high-resolution structural studies described in this report are essential for understanding the molecular basis for the sequence-specific binding as well as for new ideas for additional compound designs for sequence-specific recognition. The molecular features, in this report, explain the mechanism of recognition of both A⋅T and G⋅C bps and are an interesting molecular recognition story. Examination of the experimental structure and the NMR restrained molecular dynamics model suggests that recognition of the G⋅C base pair involves two specific H-bonds. The structure illustrates a wealth of information on different DNA interactions and illustrates an interfacial water molecule that is a key component of the complex.

Imino proton NMR guides the reprogramming of A•T specific minor groove binders for mixed base pair recognition

Sequence-specific binding to DNA is crucial for targeting transcription factor-DNA complexes to modulate gene expression. The heterocyclic diamidine, DB2277, specifically recognizes a single G•C base pair in the minor groove of mixed base pair sequences of the type AAAGTTT. NMR spectroscopy reveals the presence of major and minor species of the bound compound. To understand the principles that determine the binding affinity and orientation in mixed sequences of DNA, over thirty DNA hairpin substrates were examined by NMR and thermal melting. The NMR exchange dynamics between major and minor species shows that the exchange is much faster than compound dissociation determined from biosensor–surface plasmon resonance. Extensive modifications of DNA sequences resulted in a unique DNA sequence with binding site AAGATA that binds DB2277 in a single orientation. A molecular docking result agrees with the model representing rapid flipping of DB2277 between major and minor species. Imino spectral analysis of a ¹⁵N-labeled central G clearly shows the crucial role of the exocyclic amino group of G in sequence-specific recognition. Our results suggest that this approach can be expanded to additional modules for recognition of more sequence-specific DNA complexes. This approach provides substantial information about the sequence-specific, highly efficient, dynamic nature of minor groove binding agents.

Metal-catalyzed uncaging of DNA-binding agents in living cells†Electronic supplementary information (ESI) available: Synthesis and characterization of the studied molecules and required precursors. NMR, UV, and fluorescence spectra, titrations, control experiments, and detailed procedures for cell uptake and co-staining experiments. See DOI: 10.1039/c3sc53317dClick here for additional data file

Ruthenium-catalyzed activation of DNA-binding compounds in aqueous buffers and in cellular environments.

Attachment of alloc protecting groups to the amidine units of fluorogenic DNA-binding bisbenzamidines or to the amino groups of ethidium bromide leads to a significant reduction of their DNA affinity. More importantly, the active DNA-binding species can be readily regenerated by treatment with ruthenium catalysts in aqueous conditions, even in cell cultures. The catalytic chemical uncaging can be easily monitored by fluorescence microscopy, because the protected products display both different emission properties and cell distribution to the parent compounds.

NMR and UV studies of 3'-S-phosphorothiolate modified DNA in a DNA : RNA hybrid dodecamer duplex; implications for antisense drug design

High-resolution NMR spectroscopy has been used to establish the conformational consequences of the introduction of a single 3[prime or minute]-S-phosphorothiolate link in the DNA strand of a DNA : RNA hybrid. These systems are of interest as potential antisense therapeutic agents. Previous studies on similarly modified dinucleotides have shown that the conformation of the sugar to which the sulfur is attached shifts to the north (C(3[prime or minute])-endo/C(2[prime or minute])-exo). Comparisons made between NOESY cross-peak intensities, and coupling constants from PE-COSY spectra, for both non-modified and modified duplexes confirm that this conformational shift is also present in the double helical oligonucleotide system. In addition it is noted that in both the dinucleotides and the modified duplex, the conformation of the sugar ring 3[prime or minute] to the site of modification is also shifted to the north. That this pattern is observed in the small monomeric system as well as the larger double helix is suggestive of some pre-ordering of the sequences. The conclusion is supported by consideration of the (1)H chemical shifts of the heterocyclic bases near the site of the modification. The enhanced stability that these conformational changes should bring was confirmed by UV thermal melting studies. Subsequently a series of singly and doubly 3[prime or minute]-S-phosphorothiolate-modified duplexes were investigated by UV. The results are indicative of an additive effect of the modification with thermodynamic benefit being derived from alternate spacing of two modified linkers.

Pentamidine-induced alteration in restriction endonuclease cleavage of plasmid DNA

We have used restriction enzymes and DNaseI as probes to determine the specificity of pentamidine binding to plasmid DNA. Cleavage of plasmid pAZ130 by EcoRI, EcoRV and ApaI is inhibited by pentamidine, cleavage by XbaI, NotI and AvaI is unaffected, while cleavage by XhoI, which recognizes the same sequence as AvaI, is stimulated. DNaseI footprinting of DNA containing these restriction sites revealed that pentamidine protection is not strictly limited to AT-rich regions. We suggest that perturbation of the DNA micro- environment by pentamidine binding is responsible for its effect on nucleases.

Heterogeneous DNA binding modes of berenil

Isothermal titration calorimetry (ITC) profiles of berenil bound to different DNAs show that, despite the strong preference of berenil for AT-rich regions in DNA, it can bind to other DNA sequences significantly. The ITC results were used to quantify the binding of berenil, and the thermodynamic profiles were obtained using natural DNAs as well as synthetic polynucleotides. ITC binding isotherms cannot be simply described when a single set of identical binding sites is considered, except for poly[d(A-T)2]. Ultraviolet melting of DNA and differential scanning calorimetry were also used to quantify several aspects of the binding of berenil to salmon testes DNA. We present evidence for secondary binding sites for berenil in DNA, corresponding to G+C rich sites. Berenil binding to poly[d(G-C)2] is also observed. Circular dichroism experiments showed that binding to GC-rich sites involves drug intercalation. Using a molecular modeling approach we demonstrate that intercalation of berenil into CpG steps is sterically feasible.

Synthetic DNA minor groove-binding drugs

In this review, both cationic and neutral synthetic ligands that bind in the minor groove of DNA are discussed. Certain bis-distamycins and related lexitropsins show activities against human immunodeficiency virus (HIV)-1 and HIV-2 at low nanomolar concentrations. DAPI binds strongly to AT-containing polymers and is located in the minor groove of DNA. DAPI intercalates in DNA sequences that do not contain at least three consecutive AT bp. Berenil can also exhibit intercalative, as well as minor groove binding, properties depending on sequence. Furan-containing analogues of berenil play an important role in their activities against Pneumocystis carinii and Cryptosporidium parvuam infections in vivo. Pt(II)-berenil conjugates show a good activity profile against HL60 and U-937 human leukemic cells. Pt-pentamidine shows higher antiproliferative activity against small cell lung, non-small cell lung, and melanoma cancer cell lines compared with many other tumor cell lines. trans-Butenamidine shows good anti-P. carinii activity in rats. Pentamidine is used against P. carinii pneumonia in individuals infected with HIV who are at high risk from this infection. A comparison of the cytotoxic potencies of adozelesin, bizelesin, carzelesin, cisplatin, and doxorubicin indicates that adozelesin is a potent analog of CC-1065. Naturally occurring pyrrolo[2,1-c][l,4]benzodiazepines such as anthramycin have a 2- to 3-bp sequence specificity, but a synthetic PBD dimer spans 6 bp, actively recognizing a central 5'-GATC sequence. The crosslinking efficiency of PBD dimers is much greater than that of other major groove crosslinkers, such as cisplatin, melphalan, etc. Neothramycin is used clinically for the treatment of superficial carcinoma of the bladder.

Effects of liposome-encapsulated drugs on macrophages: comparative activity of the diamidine 4',6-diamidino-2-phenylindole and the phenanthridinium salts ethidium bromide and propidium iodide

Liposomes can be used for the intracellular delivery of drugs into macrophages. Previously, we developed a liposome-mediated macrophage 'suicide' technique based on the intraphagocytic accumulation of the liposomally delivered bisphosphonate clodronate. Later we found that the diamidine propamidine is even more effective in this approach. In the present study it is shown that liposome-encapsulated 4', 6-diamidino-2-phenylindole (L-DAPI), another well known DNA-binding diamidine, is the most effective drug in killing liver macrophages (Kupffer cells), when intravenously administered in rat. Compared to liposome-encapsulated propamidine (L-propamidine) it showed about 10-fold more activity on a molar basis. Furthermore, L-DAPI was found to induce cell death by inducing apoptosis. The structurally strongly related phenanthridinium salts ethidium bromide (EB) and propidium iodide (PI) exert marked differences in their efficacy. Whereas liposome-encapsulated PI (L-PI) was about 5 times more active in killing macrophages than L-propamidine, liposome-encapsulated EB (L-EB) showed a strongly reduced activity (10 times less than L-PI). As is shown here, PI remains mainly encapsulated in liposomes, while substantial amounts of EB leak out of liposomes. This may very well explain the differences in in vivo activity between L-EB and L-PI.

1-[(omega-aminoalkyl)amino]-4-[N-(omega-aminoalkyl)carbamoyl]-9-oxo-9, 10-dihydroacridines as intercalating cytotoxic agents: synthesis, DNA binding, and biological evaluation

A series of DNA-intercalating potential antitumor agents, 1-[(omega-aminoalkyl)amino]-4-[N-(omega-aminoalkyl)carbamoyl]-9-oxo-9, 10-dihydroacridines, has been prepared by aminolysis of the corresponding 4-[N-(omega-aminoalkyl)carbamoyl]-1-chloro derivative with a suitable omega-aminoalkylamine. The noncovalent DNA-binding properties of these bis-functionalized compounds have been examined using a combination of fluorometric and thermal denaturation techniques and are compared with the behaviors for established DNA intercalants and cationic minor groove ligands. The results indicate that (i) the agents are considerably more DNA-affinic than less functionalized acridinones, with 'apparent' binding constants of (0.1-2.1) x 10(7) and (0.3-7.5) x 10(7) M-1 at pH 5 and 7, respectively, (ii) overall affinity is sensitive to both the length of the flexible side chain and the complexity of the attached amine substituents, and (iii) the pendant side chains effect a switch to moderate AT-preferential binding. In vitro cytotoxic potencies toward six tumor cell lines broadly parallel the observed DNA affinities, although poor correlation is evident for certain compounds. The octanol/water partition coefficients have been also calculated, but there is no correlation with cytotoxicity values. Two highly DNA-affinic analogs, 10 and 13, have been identified with a useful broad spectrum of cytotoxic activity.

Solution structure of the ATF-2 recognition site and its interaction with the ATF-2 peptide

The effect of leucine zipper proteins binding to the DNA recognition site is controversial. Results from crystallography, gel and solution methods have led to opposite conclusions about the conformation of the DNA in the complex. The role of the DNA binding site in the recognition process and in the gene induction mediated by transcription factors needs to be investigated further. In this article the self-complementary 16 bp oligodeoxynucleotide (CATGTGACGTCACATG)2, which contains the cAMP response element recognised by numerous transcription factors of the leucine zipper family, has been examined free from proteins and in its interaction with the mammalian activating transcription factor 2. The recognition process has been investigated by circular dichroism analysis, which has revealed conformational changes in both DNA and protein upon binding. The solution structure of the 16mer, important in order to define the effects induced by binding of leucine zipper proteins and the intrisic bending properties of DNA, has been determined from NMR data using direct refinement against NOE intensities, analysis of scalar coupling constants and restrained molecular dynamics calculations. Final structures starting from the A and B forms of DNA agreed to a pairwise root mean square deviation (r.m.s.d.) of 1.04 +/- 0.3 A (0.7 +/- 0.2 A to the average) for all atoms. The terminal base pairs were less well determined, and the pairwise deviation of the 12 core bp was 0.83 +/- 0.27 A (0.55 +/- 0.19 A to the average). The final structures are within the B-family with an average helical twist of 36+/-2 degrees. No significant intrinsic DNA bend is shown in the activating transcription factor regulatory site. However, there are substantial deviations from the canonical B-DNA (r.m.s.d. = 3.6 A) in the core of the molecule, associated with relatively large base inclinations.

Sequence-selective binding to DNA of bis(amidinophenoxy)alkanes related to propamidine and pentamidine

The DNA sequences targeted by a complete homologous series of aromatic diamidines have been determined at single-nucleotide resolution via protection from cutting by the endonucleases DNase I, DNase II and micrococcal nuclease. Propamidine, pentamidine and to a lesser extent hexamidine bind selectively to nucleotide sequences composed of at least four consecutive A-T base pairs. In contrast, the binding to DNA of butamidine, heptamidine, octamidine and nonamidine is poorly sequence-selective. Sequences composed of only three consecutive A-T base pairs do not afford a potential binding site for propamidine or the longer homologues, and none of the drugs tolerate the presence of a G-C base pair within the binding site. Experiments with DNA molecules containing inosine in place of guanosine and 2,6-diaminopurine in place of adenine reveal that the lack of binding of propamidine to GC-containing sites is attributable to an obstructive effect of the exocyclic 2-amino group of guanosine. The present data support the view that the local conformation of the double helix (in particular the width of the minor groove) plays a dominant role in the binding reaction and that the capacity of diamidines to recognize AT-rich sequences selectively varies considerably depending on the length of the alkyl chain. The evidence indicates that binding to AT-tracts in DNA must play a role in the biological activity of these diamidines, but there is no simple correlation between binding and pharmacological efficacy.

Hydration and solution structure of d(CGCAAATTTGCG)2 and its complex with propamidine from NMR and molecular modelling

The hydration of the d(CGCAAATTTGCG)2 duplex and its complex with a propamidine reporter ligand has been examined in aqueous solution by two-dimensional NMR at two spectrometer frequencies and three temperatures. Quantitative analysis of ROESY and NOESY cross-peaks showed effective correlation times of approximately 0.5 ns at 283 K for DNA-water interactions in the major groove. In some cases the sign of the NOE inverts on changing either the temperature or spectrometer frequency. Larger effective correlation times of approximately 1 ns were observed for water interactions with A5(H2) and A6(H2) atoms located in the minor groove. Interproton NOEs and changes in chemical shifts showed that propamidine binds in the minor groove 5'-AATTT region of the host duplex, but does not displace waters adjacent to either A5(H2) or A6(H2). In the complex, the effective correlation times of these waters increase more than two-fold, possibly as a result of stabilisation due to H-bonded interaction with the amidine groups of the ligand. Hydration of the bound molecule was also found, suggesting that water may contribute to the DNA binding process for bis(amidine) drugs. Structure refinement by a NOE-restrained dynamic annealing procedure revealed that ligand binding is non-centrosymmetric with respect to the duplex, in accordance with the energetically favoured 5'-ATT (= 5'-AAT) sites predicted by analytical molecular modelling. In particular, the bound propamidine spans 3-4 base pairs in the A6-T7-T8 tract and makes close H-bonded contacts with A(N3/O4) acceptors positioned close to the minor groove floor. The refined NMR structure for the DNA-propamidine complex is compared with that determined recently using X-ray crystallographic methods.

AT selectivity and DNA minor groove binding: modelling, NMR and structural studies of the interactions of propamidine and pentamidine with d(CGCGAATTCGCG)2

A molecular modelling strategy has been developed to identify potential binding sites for bis(amidine) ligands in the minor groove of duplex DNA. Calculations of interaction energy for propamidine and pentamidine with d(CGCGAAT TCGCG)2 show that this duplex contains two symmetrically equivalent binding sites of identical affinity, each displaced by 0.3-0.4 bp from the centre of the AT segment. The ligands occupy groove sites spanning approximately 4 and 4-5 bp, respectively with asymmetric binding to the 5'-AATT sequence. The DNA-bis(amidine) interactions have been examined by high-resolution 1H-NMR. The patterns of induced changes in DNA proton chemical shift and the DNA-ligand NOEs confirm that both agents bind in the AT minor groove in a non-centrosymmetric fashion. Detailed structures were determined for each complex using a NOE-restrained simulated annealing procedure, showing that the B-type DNA conformation is not significantly altered upon complexation with either ligand. The free DNA duplex has previously been shown to be extensively hydrated in the minor groove [Kubinec, M.G. and Wemmer, D.E. (1992) J. Am, Chem. Soc. 114, 8739-8740 Liepinsh, E. Otting, G. and Wüthrich, K. (1992) Nucleic Acids Res. 20. 6549-6553]. We detect hydration water close to the A(H2) protons in the presence of propamidine, which may stabilise certain waters against exchange. This conclusion supports recent crystallographic analyses, suggesting that such ligands may use water molecules to bridge between amidinium protons and host DNA bases Details of the ligand interactions with AT-tract DNA duplexes can now be compared for the subsequences 5'-AAT, 5'-AATT and 5'-AAATTT.

The dependence of DNase I activity on the conformation of oligodeoxynucleotides

We have developed a sensitive continuous assay for nucleases using proton release. The assay has been applied to the determination of the kinetics of DNase I acting on short, defined deoxyoligonucleotides. The dependence of Kcat/K(m) on sequence and structure of short oligonucleotide substrates has been measured: increasing lengths of AnTn sequences decrease the rate of cleavage. G.A mismatches in which the bases pair using imino protons are cleaved quite effectively by DNase I. In contrast, tandem G.A mismatches which use amino pairing and have BII phosphodiesters, are refractory to DNase I. Also, the DNA strands of DNA.RNA hybrid duplexes are not cleaved by DNase I. These results show that the global conformation of a duplex and the details of its minor groove affect the cleavage efficiency by DNase I. The assay has also been used to measure the inhibition constant of the minor-groove-binding ligand propamidine. A value of 3 microM has been determined for binding to the sequence d(CGCGAATTCGCG)2, showing that dissociation constants can be determined even when there are no convenient optical signals for titrations.

Simultaneous and different binding mechanisms of 4',6-diamidino-2-phenylindole to DNA hexamer (d(CGATCG))2. A 1H NMR study

The solution structure of the complex between 4', 6-diamidino-2-phenylindole (DAPI) and DNA oligomer (d(CGATCG))2 at a 2:1 drug/duplex ratio has been characterized by combined use of proton one- and two-dimensional NMR spectroscopy, molecular mechanics, and molecular dynamics computations. Intermolecular nuclear Overhauser effects (NOEs), DNA structure perturbations, and resonance shifts induced by binding provide evidence that DAPI interacts with DNA hexamer by two different binding mechanisms, in fast exchange on the NMR time scale, without any significant distortion of the B-type conformation of DNA hexamer. The results indicate that the ligand binds into the minor groove of the central 5'-ATC-3' region of the hexamer and on the outside of the oligomer by a pi,pi-stacking interaction with the terminal C1:G6 base pairs. A model for both binding mechanisms that accounts for all experimental data was generated by molecular mechanics and dynamics calculations based on experimental NOEs. In the minor groove binding, N2 amino group of G2 precludes a deep insertion of phenyl ring of DAPI into the groove. Position and orientation of the drug in the external stacking interaction resemble those suggested for intercalation of DAPI between C:G base pairs.

Apoptosis of macrophages induced by liposome-mediated intracellular delivery of clodronate and propamidine

Liposomes can be used as vehicles for intracellular delivery of drugs into phagocytic cells. Clodronate and propamidine, delivered into macrophages in this way, will kill these cells as a result of intracellular accumulation and irreversible metabolic damage. The so-called liposome-mediated macrophage 'suicide' approach, which is based on this principle, is now frequently applied in studies aimed at unravelling macrophage function. In the present study, the mechanism of phagocytic cell death induced by liposome encapsulated drugs was investigated 'in vitro'. Peritoneal macrophages and macrophages of the RAW 264 cell line were cultured in the presence of the liposome encapsulated drugs clodronate, propamidine and several forms of ethylenediaminetetraacetic acid (EDTA). The results obtained suggest that apoptotic death is induced in phagocytic cells both by liposomally delivered clodronate and by liposomally delivered propamidine. Although intracellular EDTA did induce apoptosis in a minority of the experiments, the results support earlier findings that EDTA does not deplete macrophages as effectively as clodronate and propamidine.

Binding of Net-Fla, a netropsin-flavin hybrid molecule, to DNA: molecular mechanics and dynamics studies in vacuo and in water solution

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A. T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the drug-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.