Structure and dynamics of ligand-template interactions of topoisomerase inhibitory analogs of Hoechst 33258: high field 1H-NMR and restrained molecular mechanics studies

Structural studies on ligand-DNA systems: A robust approach in drug design

Molecular docking, molecular mechanics, molecular dynamics and relaxation matrix simulation protocols have been extensively used to generate the structural details of ligand-receptor complexes in order to understand the binding interactions between the two entities. Experimental methods like NMR spectroscopy and X-ray crystallography are known to provide structural information about ligand-receptor complexes. In addition, fluorescence spectroscopy, circular dichroism (CD) spectroscopy and molecular docking have also been utilized to decode the phenomenon of the ligand-DNA interactions, with good correlation between experimental and computational results. The DNA binding affinity was demonstrated by analysing fluorescence spectral data. Structural rigidity of DNA upon ligand binding was identified by CD spectroscopy. Docking is carried out using the DNA-Dock program which results in the binding affinity data along with structural information like interatomic distances and H-bonding, etc. The complete structural analyses of various drug-DNA complexes have afforded results that indicate a specific DNA binding pattern of these ligands. It also exhibited that certain structural features of ligands can make a ligand to be AT- or GC-specific. It was also demonstrated that changing specificity from AT base pairs to GC base pairs further improved the DNA topoisomerase inhibiting activity in certain ligands. Thus, a specific molecular recognition signature encrypted in the structure of ligand can be decoded and can be effectively employed in designing more potent antiviral and antitumour agents.

The Hoechst 33258 covalent dimer covers a total turn of the double-stranded DNA

With the goal to design ligands recognizing extended regions on dsDNA, a covalent dimer of the fluorescent dye Hoechst 33258 [bis-HT(NMe)] composed of two dye molecules linked via the phenol oxygen atoms with a (CH2)3-N+ H(CH3)-(CH2)3 fragment was constructed using computer modeling and then synthesized. Its interactions with the double-stranded DNA (dsDNA) were studied by fluorescent and UV-Vis spectroscopy and circular (CD) and linear dichroism (LD). Based on variations in the affinity to the dsDNA, it was shown that complexes of three types are formed. The first type complexes result from binding of a bis-HT(NMe) monomer in the open conformation; in this case the ligand covers the total dsDNA turn and is located in the minor groove according to the positive value of CD at 370 nm. In addition, the ability to form bis-HT(NMe)-bridges between two dsDNA molecules, i.e., each of the two bis-HT(NMe) ends binds to two different dsDNA molecules, was demonstrated for the first type complexes. Spectral characteristics (maximal absorption at 362 nm, positive sign, and maximal value of CD at 370 nm) of the first type complexes conform to those of the specific Hoechst 33258 complex with poly[d(A-T)] x poly[d(A-T]. The second type complexes correspond to the bis-HT(NMe) sandwich (as an inter- or intramolecular) binding to dsDNA with stoichiometry > or = 5 bp. Thereby, a negative LD at 360 nm and the location of bis-HT(NMe) sandwiches in the minor groove of B form dsDNA seems contradictory. Spectral characteristics (maximal positive CD at 345 nm, a dramatic decrease in fluorescence intensity and the shift of its maximum to 490 nm) of these complexes favor a suggestion that this binding correlates to the formation of nonspecific dimeric Hoechst 33258 complex with dsDNA. The third type complexes are characterized by stoichiometry of one bis-HT(NMe) molecule per approximately 2 bp and the tendency to zero of LD values at 270 and 360 nm. We assume that in these complexes bis-HT(NMe) sandwich dimers are formed on dsDNA. The complexes of this type conform to the aggregation type complex of Hoechst 33258 with dsDNA. The ability of bis-HT(NMe) to cover the whole dsDNA turn or form bridges with two dsDNA upon the formation of the first type complexes essentially distinguishes it from Hoechst 33258, which can only occupy 5 bp and does not form such bridges. This specific property of bis-HT(NMe) may support new biological activities.

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.

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.

Multimode interaction of Hoechst 33258 with eukaryotic DNA; quantitative analysis of the DNA conformational changes

The interaction of the minor groove binding ligand Hoechst 33258 (Hoe) with natural DNA was investigated by high resolution titration rotational viscometry. Analysis of the concomitant DNA conformational changes was performed with two DNA samples of sufficiently different molar mass M, at 4 degrees C, 22 degrees C and 40 degrees C, for Hoe/DNA-P ratios below r = 0.02. In this narrow r range several interaction modes could be resolved. The measured conformational changes were quantified in terms of relative changes of both apparent DNA persistence length, delta a/a, and hydrodynamically operative DNA contour length, deltaL/L. Delta a/a(r) primarily is a measure of ligand-induced DNA helix stiffening, but both, delta a/a(r) and deltaL/L(r), generally depend also on ligand binding induced DNA bending or DNA unbending. The essential difference obviously is that delta a/a(r) is influenced by the randomly distributed helix bends and deltaL/L(r) by phased ones. The measurements performed at different temperatures deliver informations about existence and temperature dependent abolition of intrinsic helix curvature. Both Hoe and netropsin (Nt) prefer binding to AT rich DNA segments, which are candidates for intrinsic DNA helix bends. But our data for Hoe interaction with calf thymus DNA (ctDNA) show characteristic differences to those for Nt-ctDNA interaction. Especially for Hoe, the mode of highest affinity is saturated already at a ligand concentration of roughly 1 nM (r approximately = 0.0015 Hoe/DNA-P). It exhibits an unusually strong temperature dependence of the conformational DNA response. A Hoe-Nt competition experiment shows that Hoe binding to the sites of the very first Hoe mode is almost unaffected by bound Nt. But Hoe binding to the sites of the following Hoe modes does not occur due to the competition with Nt. Thus this mode of strongest Hoe-DNA interaction reflects a unique mechanism, possibly of high relevance for gene regulatory systems.

Crystal structure of the DNA decamer d(CGCAATTGCG) complexed with the minor groove binding drug netropsin

The crystal structure of netropsin bound to the decamer d(CGCAATTGCG) has been determined at 2.4 A resolution. This is the first example of a crystal structure of netropsin bound to decamer DNA. The central eight bases of each DNA single-strand base pair with a self-complementary strand to form an octamer B-DNA duplex. These duplexes lie end to end within the unit cell. The terminal 5'-C and G-3' bases are unpaired and interact with the neighboring duplexes via interactions within both the major and minor groove to form base triplet interactions of the type C(+)-G x C and G*(G x C), respectively. The triplet interaction of the type C(+)-G x C is known to exist within triplex DNA with the C+ base oriented parallel with the Watson-Crick guanine base to which it hydrogen bonds. The netropsin molecule lies within the minor groove of the octamer duplex and assumes a class I type position, with bifurcated hydrogen-bonding interactions from the amide groups of the netropsin to the A x T base pairs of the minor groove. The netropsin molecule fits within a five base pair long minor groove site by bending of the flexible amidinium group at one end of the drug.

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.

Torsional dynamics and orientation of DNA--DAPI complexes

The flexibility of calf thymus DNA and several polynucleotides was measured using the anisotropy decay of DAPI bound to DNA, a minor groove probe. DNA torsional dynamics were analyzed using the Schurr model [Allison, S. A., & Schurr, J. M. (1979) Chem. Phys. 41, 35-44] in the infinite polymer length approximation. Time-resolved fluorescence depolarization was measured using a frequency-double mode-locked dye laser and frequency-domain acquisition methods. At very high P/D ratios, the anisotropy decay is dominated by DNA torsional dynamics. The recovered values of the torsional elastic constant were in good agreement with literature values obtained using other DNA probes. The exact knowledge of the angle between the probe emission dipole transition moment and the helix axis is critical for the determination of the polymer elastic constant. At low P/D ratios, energy transfer between dye molecules strongly contributes to the anisotropy decay. We have developed a statistical model that describes the anisotropy decay, when the correct geometrical factors are included. At low P/D ratios the anisotropy decay is dominated by fluorescence homotransfer. In this regime, it is possible to determine the orientation of the dye molecule with respect to the polymer with accuracy. The values obtained for the distance and orientation of the DAPI molecules in solution using the fluorescence measurements are in excellent agreement with those from the crystal structure of the oligonucleotides-DAPI complex by Dickerson's group [Larsen T.A., Goodsell, D. S., Cascio, D., Grzeskowiak, K., & Dickerson, R. E. (1989) J. Biomol. Struct. Dyn. 7, 477-491].

Molecular mechanics calculations on the complexes between analogues of Hoechst 33258 and d(CGCGAAT-TCGCG)2: influence of bulky group substitution on base pair preference of DNA minor groove binders

Molecular mechanics calculations were performed on the three structures of the complexes formed by the derivatives of Hoechst 33258 and dodecameric DNA duplex d(CGCGAATTCGCG)2. Formation and docking energies of these complexes were compared. It was found that the CG site that is 3′ to the central AATT region can be tolerated by the drugs. This is probably due to the presence of the bulky piperazine ring and, more pronouncedly, by alkylated analogues of the drug that prefer the wider minor groove formed by the GC base pair region of B-DNA. The argument of bulkiness of the piperazine moiety as the origin of enhancement of GC affinity is supported by detailed structural analysis of the intermolecular interface and widening of the DNA minor groove at the binding site. Implications of the results are discussed. Keywords: minor groove binder, docking energy, sequence specificity.

Mode of DNA binding of bis-benzimidazoles and related structures studied by electric linear dichroism

The binding mode of a series of bis-benzimidazole analogues of Hoechst 33258 to a variety of DNAs and polynucleotides has been investigated by electric linear dichroism. Two groups of compounds were examined: (i) benzoxazole and pyridoimidazole derivatives and (ii) pyridoimidazole analogs substituted with an N-alkoxyalkyl group either directed towards the minor groove or directed away from the minor groove. The ELD data indicate that the mode of binding of these drugs varies significantly with the sequence of the target DNA sequence. The DNA binding properties of these drugs are related to their topoisomerase inhibitory properties.

The different binding modes of Hoechst 33258 to DNA studied by electric linear dichroism

The binding mode of the bisbenzimidazole derivative Hoechst 33258 to a series of DNAs and polynucleotides has been investigated by electric linear dichroism. Positive reduced dichroisms were measured for the poly(dA-dT).poly(dA-dT)- and poly(dA).poly(dT)-Hoechst complexes in agreement with a deep penetration of the drug into the minor groove. Similarly, the drug displays positive reduced dichroism in the presence of the DNAs from calf thymus, Clostridium perfringens and Coliphage T4. Conversely, negative reduced dichroisms were obtained when Hoechst 33258 was bound to poly(dG-dC).poly(dG-dC), poly(dA-dC).poly(dG-dT) and poly(dG).poly(dC) as well as with the GC-rich DNA from Micrococcus lysodeikticus indicating that in this case minor groove binding cannot occur. Substitution of guanosines for inosines induces a reversal of the reduced dichroism from negative to positive. Therefore, as anticipated it is the 2-amino group of guanines protruding in this groove which prevents Hoechst 33258 from getting access to the minor groove of GC sequences. The ELD data obtained with the GC-rich biopolymers are consistent with an intercalative binding. Competition experiments performed with the intercalating drug proflavine lend credence to the involvement of an intercalative binding rather than to an external or major groove binding of Hoechst 33258 at GC sequences.