Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2)

Hydrophobic Amino Acids as Universal Elements of Protein-Induced DNA Structure Deformation

Interaction with the DNA minor groove is a significant contributor to specific sequence recognition in selected families of DNA-binding proteins. Based on a statistical analysis of 3D structures of protein–DNA complexes, we propose that distortion of the DNA minor groove resulting from interactions with hydrophobic amino acid residues is a universal element of protein–DNA recognition. We provide evidence to support this by associating each DNA minor groove-binding amino acid residue with the local dimensions of the DNA double helix using a novel algorithm. The widened DNA minor grooves are associated with high GC content. However, some AT-rich sequences contacted by hydrophobic amino acids (e.g., phenylalanine) display extreme values of minor groove width as well. For a number of hydrophobic amino acids, distinct secondary structure preferences could be identified for residues interacting with the widened DNA minor groove. These results hold even after discarding the most populous families of minor groove-binding proteins.

Cooperativity effect of the ππ interaction between drug and DNA on intercalative binding induced by H-bonds: a QM/QTAIM investigation of the curcuminadenineH2O model system

In order to reveal the nature of intercalative binding of drug to DNA, the cooperativity effect of the ππ interaction was investigated in the curcuminadenineH2O model system by applying a combined QM and QTAIM computational approach. The H-bonds between the electron-donating group of curcumin and adenine induce the formation of the ππ stacking. The introduction of H2O weakens the H-bonding and ππ interactions, leading to an anti-cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that the anti-cooperative effect is the main driving force for the intercalative binding of drug to DNA bases, which is in agreement with many experimental phenomena. Therefore, the designed DNA-targeted intercalating drugs should possess not only hydrophobic moieties, but also strong electron-donating groups bound to the DNA bases with H-bonds, which can slow the variation rates of the strengths of the H-bonding and ππ interactions between drug and DNA bases in the anti-cooperative process, leading to the intercalation formation. The enthalpy change is the major factor driving the positive thermodynamic cooperativity.

Cooperativity effect involving drug-DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H2O system

In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug-DNA/RNA intermolecular interaction in ketoprofen⋯cytosine⋯H₂O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N-H⋯O, O-H⋯O, O-H⋯N, C-H⋯N, and C-H⋯O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket⋯cyt⋯H₂O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H₂O, the anti-cooperativity effect plays a dominant role in the drug-DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H₂O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.

Minor groove binding compounds that jump a gc base pair and bind to adjacent AT base pair sites

Most A/T specific heterocyclic diamidine derivatives need at least four A/T base pairs for tight binding to the DNA minor groove. Addition of a GC base pair to A/T sequences typically causes a large decrease in binding constant. The ability to target biologically important sequences of DNA could be significantly increased if compounds that could recognize A/T sites with an intervening GC base pair could be designed. The kinetoplast DNA sequence of parasitic microorganisms, for example, contains numerous three A/T binding sites that are separated by a single G. A series of compounds were prepared to target the AAAGTTT sequence as a model system for discovery of "G-jumpers". The new synthetic compounds have two aromatic-amidine groups for A/T recognition, and these are connected through an oxy-methylene linker to cross the GC. CD experiments indicated a minor groove binding mode, as expected, for these compounds. T(max), surface plasmon resonance, and isothermal titration calorimetry experiments revealed 1:1 binding to the AAAGTTT sequence with an affinity that depends on compound structure. Benzimidazole derivatives gave the strongest binding and had generally good solution properties. The binding affinities to the classical AATT sequence were similar to that for AAAGTTT for these extended compounds, but binding was weaker to the AAAGCTTT sequence with two intervening GC base pairs. Binding to both AAAGTTT and AATT was enthalpy driven for strong binding benzimidazole derivatives.

Hoechst 33258--poly(dG-dC).poly(dG-dC) complexes of three types

It was found recently that Hoechst 33258, a dsDNA fluorescent dye used in cytological studies, is an efficient inhibitor of the interaction of TATA-box-binding protein with DNA, DNA topoisomerase I, and DNA helicases. In addition it proved to be a radioprotector. Biological activity of Hoechst 33258 may be associated with dsDNA complexes of not only monomeric, but also dimeric type. In this work, the Hoechst 33258 interaction with poly(dG-dC).poly(dG-dC) was studied using UV-vis and fluorescent spectroscopy, circular and flow-type linear dichroism. It was found that Hoechst 33258 formed with poly(dG-dC).poly(dG-dC) complexes of three types, namely, monomeric, dimeric, and, apparently, tetrameric, and their spectral properties were studied. Complexes of monomeric and dimeric types competed with distamycin A, a minor groove ligand, for binding to poly(dG-dC).poly(dG-dC). We proposed that Hoechst 33258 both monomers and dimers form complexes of the external type with poly(dG-dC).poly(dG-dC) from the side of the minor groove.

Influence of DNA structure on adjacent site cooperative binding

Previous NMR studies of Hoechst 33258 with the d(CTTTTGCAAAAG)2 sequence have shown very strong (K2 >> K1) cooperativity between two adjacent binding sites (Searle, M. S.; Embrey, K. J. Nucleic Acids Res. 1990, 18 (13), 3753- 3762). In contrast, surface plasmon resonance (SPR) results with the hairpin analog of the same sequence show significantly reduced cooperativity. In an effort to explain the difference, two-dimensional (2-D) NMR experiments were done on both duplex and hairpin. Hoechst 33258 and an amidine analog, DB183, show very strong cooperativity with the duplex DNA but much weaker cooperativity with the hairpin. The significantly lower thermal melting temperature (Tm) of the duplex (34.8 degrees C) in comparison to its hairpin analog (62.3 degrees C) supports the idea of a dynamic difference between the two DNA structures that can influence cooperativity in binding. These results confirm the role of conformational entropy in positive cooperativity in some DNA interactions.

Calorimetric and spectroscopic studies of Hoechst 33258: self-association and binding to non-cognate DNA

Sequence and structure-specific molecular recognition of DNA by small molecules is an important goal in biophysical chemistry and drug discovery. Many candidate ligands possess flat aromatic surfaces and other molecular features that allow them to self-associate. In addition, non-specific binding to the target is a complicating feature of these interactions. Therefore, multiple equilibria are present and need to be accounted for in data analysis in order to obtain meaningful thermodynamic parameters. In order to address these issues we have systematically examined the bis-benzimidazole dye Hoechst 33258 (H33258) in terms of self-aggregation and binding to DNA oligonucleotides lacking any cognate minor groove A.T sites. This model system has been interrogated using isothermal titration calorimetry (ITC), circular dichroism (CD), fluorescence spectroscopy and pulsed gradient spin echo NMR. Three distinct binding events and ligand self-aggregation have been identified and, where possible, quantified. H33258 self-aggregation involves a step-wise aggregation mechanism, driven by stacking interactions. The DNA binding process includes two specific binding modes and non-specific DNA-templated H33258 stacking. We have written novel ITC data-fitting software (IC-ITC; freely available to the biophysics community), which simultaneously fits ligand aggregation and ligand-DNA binding. Here, this numerical analysis, which uses simulated annealing of complex calorimetric data representing multiple coupled equilibria, is described.

Expanding the specificity of DNA targeting by harnessing cooperative assembly

The precise control of developmental and regulatory processes in the cell requires accurate recognition of specific DNA sites. For genomes as large as that of humans, single-molecule-DNA binders have difficulties accurately and specifically recognizing the intended targets. Natural transcription factors overcome these difficulties by forming non-covalent complexes on the DNA with other transcription factors. These cooperative complexes overcome the difficulties of single-molecule transcription factors, allowing specific, combinatorial control of a range of transcriptional targets. Artificial transcription factors have been designed to take advantage of this technique of cooperative assembly, facilitating future studies in whole genome targeting. In contrast to a simple model of component independence, cooperative complexes as a whole often display slightly altered DNA specificity from what would be expected from the analysis of their separate components. The true sequence specificity of cooperative complexes, and thus their presumed in vivo targets, have to be experimentally probed. A number of techniques, such as the cognate site identity array, now allow for rapid, high-throughput determination of the specificity of cooperative complexes.

Multiple Binding Modes for Dicationic Hoechst 33258 to DNA

The binding of dicationic Hoechst 33258 (ligand) to DNA was characterized by means of the fluorescence spectra, fluorescence intensity titration, time-resolved fluorescence decay, light scattering, circular dichroism, and fluorescence thermal denaturation measurements, and two binding modes were distinguished by the experimental results. Type 1 binding has the stoichiometry of one ligand to more than 12 base pairs, and it is defined as quasi-minor groove binding which has the typical prolonged fluorescence lifetime of about 4.4 ns. In type 1 binding, planar conformation of the ligand is favorable. Type 2 binding with phosphate to ligand ratio (P/L) < 2.5 has the stoichiometry of one ligand to two phosphates. It is defined as a highly dense and orderly stacked binding with DNA backbone as the template. Electrostatic interactions between doubly protonated ligands and negatively charged DNA backbone play a predominant role in the type 2 binding mode. The characteristics of this type of binding result in a twisted conformation of the ligand that has a fluorescence lifetime of less than 1 ns. The results also indicate that the binding is in a cooperative manner primarily by stacking of the aromatic rings of the neighboring ligands. Type 1 binding is only observed for double-stranded DNA (dsDNA) with affinity constant of 1.83 x 10(7) M-1. In the type 2 binding mode, the binding affinity constants are 4.9 x 10(6) and 4.3 x 10(6) M-1 for dsDNA and single-stranded DNA (ssDNA), respectively. The type 2 binding is base pair independent while the type 1 binding is base pair related. The experiments described in this paper revealed that the dication bindings are different from the monocation bindings reported by previous study. The dication binding leads to stronger aggregation at low ligand concentration and results in orderly arrangements of the ligands along DNA chains. Furthermore the dication binding is demonstrated to be beneficial for enhancing the DNA's stability.

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.

Impedance sensing of DNA binding drugs using gold substrates modified with gold nanoparticles

Interfacial interactions between immobilized DNA probes and DNA-specific sequence binding drugs were investigated using impedance spectroscopy toward the development of a novel biosensing scheme. The impedance measurements are based on the charge-transfer kinetics of the [Fe(CN)6]3-/4- redox couple. Compared to bare gold surfaces, the immobilization of DNA and then the DNA-drug interaction on electrode surfaces altered the capacitance and the interfacial electron resistance and thus diminished the charge-transfer kinetics by reducing the active area of the electrode or by preventing the redox species from approaching the electrode. Electrochemical deposition of gold nanoparticles on a gold electrode surface showed significant improvement in sensitivity. DNA-capped gold nanoparticles on electrodes act as selective sensing interfaces with tunable sensitivity due to higher amounts of DNA probes and the concentric orientation of the DNA self-assembled monolayer. The specificity of the interactions of two classical minor groove binders, mythramycin, a G-C specific-DNA binding anticancer drug, netropsin, an A-T specific-DNA binding drug and an intercalator, nogalamycin on AT-rich DNA-modified substrate and GC-rich DNA-modified substrate are compared. Using gold nanoparticle-deposited substrates, impedance spectroscopy resulted in a 20-40-fold increase in the detection limit. Arrays of deposited gold nanoparticles on gold electrodes offered a convenient tool to subtly control probe immobilization to ensure suitably adsorbed DNA orientation and accessibility of other binding molecules.

Configurational entropy change of netropsin and distamycin upon DNA minor-groove binding

Binding of a small molecule to a macromolecular target reduces its conformational freedom, resulting in a negative entropy change that opposes the binding. The goal of this study is to estimate the configurational entropy change of two minor-groove-binding ligands, netropsin and distamycin, upon binding to the DNA duplex d(CGCGAAAAACGCG).d(CGCGTTTTTCGCG). Configurational entropy upper bounds based on 10-ns molecular dynamics simulations of netropsin and distamycin in solution and in complex with DNA in solution were estimated using the covariance matrix of atom-positional fluctuations. The results suggest that netropsin and distamycin lose a significant amount of configurational entropy upon binding to the DNA minor groove. The estimated changes in configurational entropy for netropsin and distamycin are -127 J K(-1) mol(-1) and -104 J K(-1) mol(-1), respectively. Estimates of the configurational entropy contributions of parts of the ligands are presented, showing that the loss of configurational entropy is comparatively more pronounced for the flexible tails than for the relatively rigid central body.

Exploring the potential of complex formation between a mutant DNA and the wild type protein counterpart: a MM and MD simulation approach

We have demonstrated that the methods of molecular modeling and molecular dynamics simulation might be used to assess whether a specific mutation in the DNA would destabilize a known DNA-protein complex. The approach is based on probing the changes in the interaction that would be induced into the complex if within the already formed wild type complex the mutation could be introduced. We have used Hoxc8-DNA complex as a test system where it is known that the Hoxc8 binding affinity of the DNA is completely lost upon mutation of the DNA by replacing TAAT stretch to GCCG. Mutation was obtained by changing the relevant base pairs into the DNA of the model of the corresponding wild type complex developed by homology modeling and MD simulation in water for 2.0 ns. Comparison of the structure, dynamics and interactions between the hypothetical mutant model with those of the similarly refined wild type model shows that the loss of affinity of the mutant DNA to Hoxc8 has two different origins: (i) loss of several strong H-bonds as the direct consequences of mutation and (ii) reduced H-bonds in the common parts due to a net loss or inferior H-bonding geometry induced by the mutation as indirect effects. The net change in the interaction energy between the DNA and the protein in the best possible configuration indicated the experimentally observed destabilization effects. No significant change in the groove width was observed and no correlation was found between the water-bridges and the loss of affinity.

Principal component analysis of DNA oligonucleotide structural data

The microstructure of a DNA helix is characterized by several base pair and base step parameters such as twist, rise, roll, propeller twist, etc., in addition to conformational parameters such as the backbone and the glycosidic torsion angles. Among these only a few, which are independent of all others and of each other, may be used to precisely characterize the helix. The problem however is to identify these independent parameters. We have used principal component analysis to identify a relatively small set of independent parameters, with which to characterize each DNA helix. We show that these principal components clearly discriminate between A and B DNA helical types. The calculations further suggest that the microstructure of a DNA helix is better characterized using dinucleotides.

Acid-Base and Electronic Structure-Dependent Properties of Hoechst 33342

An ab initio quantum mechanical study of the bis-benzimidazole derivative Hoechst 33342 is presented. Specifically, we report on the molecular geometry, proton affinity, gas-phase basicity, hydration, pKa, and molecular interaction potential map of Hoechst 33342 and the forms obtained by its protonation and deprotonation. Results have been compared with those obtained for Hoechst 33258 to analyze the influence of the chemical structure on the properties of this bis-benzimidazole derivative. Finally, the connection of some theoretical predictions with experimental evidences has been examined.

DNA Sequence-Specific Ligands: XII. Synthesis and Cytological Studies of Dimeric Hoechst 33258 Molecules

We synthesized dimeric Hoechst dye molecules composed of two moieties of the Hoechst 33258 fluorescent dye phenolic hydroxy groups of which were tethered via pentamethylene, heptamethylene, or triethylene oxide linkers. A characteristic pattern of differential staining of chromosome preparations from human premonocytic leukemia HL60 cells was observed for all the three fluorescent dyes. The most contrast pattern was obtained for the bis-Hoechst analogue with the heptamethylene linker; its quality was comparable with the picture obtained in the case of chromosome staining with 4',6-diamidino-2-phenylindole. The ability to penetrate into the live human fibroblasts was studied for the three bis-Hoechst compounds. The fluorescence intensity of nuclei of live and fixed cells stained with the penta- and heptamethylene-linked bis-Hoechst analogues was found to differ only slightly, whereas the fluorescence of the nuclei of live cells stained with triethylene oxide-linked bis-Hoechst was considerably weaker than that of the fixed cells. The bis-Hoechst molecules are new promising fluorescent dyes that can both differentially stain chromosome preparations and penetrate through cell and nuclear membranes and effectively stain cell nuclei.

Compensating bends in a 16-base-pair DNA oligomer containing a T(3)A(3) segment: A NMR study of global DNA curvature

In-phase ligated DNA containing T(n)A(n) segments fail to exhibit the retarded polyacrylamide gel electrophoresis (PAGE) migration observed for in-phase ligated A(n)T(n) segments, a behavior thought to be correlated with macroscopic DNA curvature. The lack of macroscopic curvature in ligated T(n)A(n) segments is thought to be due to cancellation of bending in regions flanking the TpA steps. To address this issue, solution-state NMR, including residual dipolar coupling (RDC) restraints, was used to determine a high-resolution structure of [d(CGAGGTTTAAACCTCG)2], a DNA oligomer containing a T3A3 tract. The overall magnitude and direction of bending, including the regions flanking the central TpA step, was measured using a radius of curvature, Rc, analysis. The Rc for the overall molecule indicated a small magnitude of global bending (Rc = 138 +/- 23 nm) towards the major groove, whereas the Rc for the two halves (72 +/- 33 nm and 69 +/- 14 nm) indicated greater localized bending into the minor groove. The direction of bending in the regions flanking the TpA step is in partial opposition (109 degrees), contributing to cancellation of bending. The cancellation of bending did not correlate with a pattern of roll values at the TpA step, or at the 5' and 3' junctions, of the T3A3 segment, suggesting a simple junction/roll model is insufficient to predict cancellation of DNA bending in all T(n)A(n) junction sequence contexts. Importantly, Rc analysis of structures refined without RDC restraints lacked the precision and accuracy needed to reliably measure bending.

A host-guest approach for determining drug-DNA interactions: an example using netropsin

Netropsin is a well-characterized DNA minor groove binding compound that serves as a model for the study of drug-DNA interactions. Our laboratory has developed a novel host-guest approach to study drug-DNA interactions in which the host, the N-terminal fragment of Moloney murine leukemia virus reverse transcriptase (MMLV RT) is co-crystallized with a DNA oligonucleotide guest in the presence and absence of drug. We have co-crystallized netropsin with the RT fragment bound to the symmetric 16mer d(CTTAATTCGAATTAAG)2 and determined the structure of the complex at 1.85 A. In contrast to previously reported netropsin-DNA structures, our oligonucleotide contains two AATT sites that bind netropsin with flanking 5' and 3' sequences that are not symmetric. The asymmetric unit of the RT fragment-DNA-netropsin crystals contains one protein molecule and one-half of the 16mer with one netropsin molecule bound. The guanidinium moiety of netropsin binds in a narrow part of the minor groove, while the amidinium is bound in the widest region within the site. We compare this structure to other Class I netropsin-DNA structures and find that the asymmetry of minor groove widths in the AATT site contributes to the orientation of netropsin within the groove while hydrogen bonding patterns vary in the different structures.

Large-scale molecular dynamics simulation of DNA: implementation and validation of the AMBER98 force field in LAMMPS

Molecular modelling played a central role in the discovery of the structure of DNA by Watson and Crick. Today, such modelling is done on computers: the more powerful these computers are, the more detailed and extensive can be the study of the dynamics of such biological macromolecules. To fully harness the power of modern massively parallel computers, however, we need to develop and deploy algorithms which can exploit the structure of such hardware. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) is a scalable molecular dynamics code including long-range Coulomb interactions, which has been specifically designed to function efficiently on parallel platforms. Here we describe the implementation of the AMBER98 force field in LAMMPS and its validation for molecular dynamics investigations of DNA structure and flexibility against the benchmark of results obtained with the long-established code AMBER6 (Assisted Model Building with Energy Refinement, version 6). Extended molecular dynamics simulations on the hydrated DNA dodecamer d(CTTTTGCAAAAG)(2), which has previously been the subject of extensive dynamical analysis using AMBER6, show that it is possible to obtain excellent agreement in terms of static, dynamic and thermodynamic parameters between AMBER6 and LAMMPS. In comparison with AMBER6, LAMMPS shows greatly improved scalability in massively parallel environments, opening up the possibility of efficient simulations of order-of-magnitude larger systems and/or for order-of-magnitude greater simulation times.

Determination of DNA minor groove width in distamycin-DNA complexes by solid-state NMR

We have performed solid-state 31P-19F REDOR nuclear magnetic resonance (NMR) experiments to monitor changes in minor groove width of the oligonucleotide d(CGCAAA2'FUTGGC)*d(GCCAAT(pS)TT GCG) (A3T2) upon binding of the drug distamycin A at different stoichiometries. In the hydrated solid-state sample, the minor groove width for the unbound DNA, measured as the 2'FU7-pS19 inter-label distance, was 9.4 +/- 0.7 A, comparable to that found for similar A:T-rich DNAs. Binding of a single drug molecule is observed to cause a 2.4 A decrease in groove width. Subsequent addition of a second drug molecule results in a larger conformational change, expanding this minor groove width to 13.6 A, consistent with the results of a previous solution NMR study of the 2:1 complex. These 31P-19F REDOR results demonstrate the ability of solid-state NMR to measure distances of 7-14 A in DNA-drug complexes and provide the first example of a direct spectroscopic measurement of minor groove width in nucleic acids.

Influence of phenyl ring disubstitution on bisbenzimidazole and terbenzimidazole cytotoxicity: synthesis and biological evaluation as radioprotectors

DNA minor groove binders, Hoechst 33258 and Hoechst 33342, have been reported to protect against radiation-induced DNA-strand breakage, but their mutagenicity and cytotoxicity limit their use as protectors of normal tissue during radiotherapy and as biological radioprotectors during accidental radiation exposure. On the basis of these observations, two new nontoxic disubstituted benzimidazoles were synthesized, one having two methoxy groups (5-(4-methylpiperazin-1-yl)-2-[2'-(3,4-dimethoxyphenyl)-5'-benzimidazolyl]benzimidazole, 5) and another having a methoxy and a hydroxyl group (5-(4-methylpiperazin-1-yl)-2-[2'[2''-(4-hydroxy-3-methoxyphenyl)-5' '-benzimidazolyl]-5'-benzimidazolyl]benzimidazole, 6) ortho to each other on the phenyl ring. The radiomodifying effects of these nontoxic ligands were investigated with a human glioma cell line exposed to low linear energy transfer radiation by determining cell survival and cell proliferation compared with effects of the parent compound, Hoechst 33342. Cytotoxicity assayed by analyzing clonogenicity, cell growth, and metabolic viability showed that both 5 and 6 were nontoxic at 100 microM after 72 h of exposure, whereas Hoechst 33342 resulted in lysis of 77% of these cells in 24 h. Macrocolony assay (clonogenicity) showed that 73%, 92%, and 10% of the cells survived when treated with 100 microM 5, 6, and Hoechst 33342, respectively. Both 5 and 6 did not affect the growth of BMG-1 cells. At 10 microM, 5 and 6 showed 82% and 37% protection against radiation-induced cell death (macrocolony assay) while 100% protection was observed against growth inhibition. Disubstitution of the phenyl ring has not only reduced cytotoxicity but also enhanced DNA-ligand stability, conferring high degree of radioprotection.

Hydration of hydroxypyrrole influences binding of ImHpPyPy-beta-Dp polyamide to DNA

Ligands which are able to recognize DNA sequence specifically are of fundamental interest as transcription controlling drugs. Recently a polyamide ligand was developed (ImHpPyPy-beta-Dp) which differentiates in a dimeric arrangement between all four possible base pair steps in the minor groove. This is a landmark for the design of DNA binding drugs because it was believed that such a recognition could only be possible in the major groove of DNA. Although the OH groups of the hydroxypyrrole (Hp) moieties of the ligands are responsible for this sequence discrimination, experiments showed that this OH group also reduces the absolute binding constant. We performed a free energy calculation by means of thermodynamic integration in order to find out the influence of this single hydroxyl on DNA binding. In our simulation, we found that the hydroxyl group reduces binding by about 1.3 kcal/mol, which is in excellent agreement with the experimentally determined value of 1.2 kcal/mol. In further MD simulations, the structural reasons for this reduction was estimated. The results of these simulations qualitatively agree with the X-ray structures, but in contrast, in the simulations both (ImHpPyPy-beta-Dp and ImPyPyPy-beta-Dp) ligand-DNA (d(CCAGTACTGG)(2)) complexes exhibit only slight structural differences. This is consistent with a recently published second pair of similar polyamide DNA crystal structures. Thus, we believe that the explanations resulting from the X-ray structures must be modified. We attribute the large structural differences between the two polyamide DNA complexes to a buffer molecule which binds only in the case of the ImHpPyPy-beta-Dp-DNA complex at the region of interest. We propose that the differential hydration of both ligands in the unbound state is responsible for the reduction of the binding constant. Additionally, we suggest an indirect readout of DNA, because of a lengthening of the Watson-Crick base pairs, which possibly contributes to the differentiation between T.A, A.T from G.C, C.G base pairs.

NMR characterization of the DNA binding properties of a novel Hoechst 33258 analogue peptide building block

A novel aryl-bis-benzimidazole amino acid analogue of the DNA-binding compound Hoechst 33258 has recently been designed for incorporation in peptide combinatorial libraries by replacing the N-methylpiperazine group with a carboxyl group and the hydroxy group with an amino-methyl group. The DNA-binding properties of the aryl-bis-benzimidazole monomer with the C-terminus derivatized with 3-(dimethylamino)-propylamine has been investigated in this paper by (1)H NMR studies of two different complexes with two different DNA sequences: A(5) d(5'-GCCA(5)CG-3'):d(5'-CGT(5)GGC-3') and A(3)T(3) d(5'-CGA(3)T(3)CG-3')(2). Chemical shift footprinting shows that the ligand binds at the center of the A(3)T(3) sequence but at the 3'-end of A(5). A large number of NOEs show a well-defined complex with the ligand situated at the center of the palindromic A(3)T(3) but with the asymmetric A(5) the ligand binds with an orientational preference with the bis-benzimidazole moiety displaced toward the 3'-end from the center of the duplex. Two families of models of the complexes with A(5) and A(3)T(3) were derived with restrained molecular dynamics based on a large set of 70 and 61, respectively, intermolecular ligand NOEs. Both models give a picture of a tightly fitting ligand with close van der Waals contacts with the walls of the minor groove and with the two benzimidazole and the amide hydrogens involved in bifurcated cross-strand hydrogen bonds to adenine N3 and thymine O2. The minor groove width of the models correlate well with the binding site of the ligand, and the orientational preference is argued to be a consequence of the minor groove width and hydrogen bonding.

Determination of affinity, stoichiometry and sequence selectivity of minor groove binder complexes with double-stranded oligodeoxynucleotides by electrospray ionization mass spectrometry

Electrospray mass spectrometry was evaluated regarding the reliability of the determination of the stoichiometries and equilibrium association constants from single spectra. Complexes between minor groove binders (Hoechst 33258, Hoechst 33342, DAPI, netropsin and berenil) and 12mer oligonucleotide duplexes with a central sequence (A/T)4 flanked by G/C base pairs were chosen as model systems. To validate the electrospray ionization mass spectrometry (ESI-MS) method, comparisons were made with circular dichroism and fluorescence spectroscopy measurements. ESI-MS allowed the detection of minor (2 drug + DNA) species for Hoechst 33258, Hoechst 33342, DAPI and berenil with duplex d(GGGG(A/T)4GGGG). d(CCCC(A/T)4CCCC), which were undetectable with the other techniques. Assuming that the duplexes and the complexes have the same electrospray response factors, the equilbrium association constants of the 1:1 and 2:1 complexes were determined by ESI-MS, and the values show a good quantitative agreement with fluorescence determined constants for Hoechst 33258 and Hoechst 33342. It is also shown that ESI-MS can quickly give reliable information on the A/T sequence selectivity of a drug: the signal of a complex is directly related to the affinity of the drug for that particular duplex. The potential of ESI-MS as a qualitative and quantitative affinity screening method is emphasized.

Thermodynamics of drug-DNA interactions

Many anticancer, antibiotic, and antiviral drugs exert their primary biological effects by reversibly interacting with nucleic acids. Therefore, these biomolecules represent a major target in drug development strategies designed to produce next generation therapeutics for diseases such as cancer. In order to improve the clinical efficacy of existing drugs and also to design new ones it is necessary to understand the molecular basis of drug-DNA interactions in structural, thermodynamic, and kinetic detail. The past decade has witnessed an increase in the number of rigorous biophysical studies of drug-DNA systems and considerable knowledge has been gained in the energetics of these binding reactions. This is, in part, due to the increased availability of high-sensitivity calorimetric techniques, which have allowed the thermodynamics of drug-DNA interactions to be probed directly and accurately. The focus of this article is to review thermodynamic approaches to examining drug-DNA recognition. Specifically, an overview of a recently developed method of analysis that dissects the binding free energy of these reactions into five component terms is presented. The results of applying this analysis to the DNA binding interactions of both minor groove drugs and intercalators are discussed. The solvent water plays a key role in nucleic acid structure and consequently in the binding of ligands to these biomolecules. Any rational approach to DNA-targeted drug design requires an understanding of how water participates in recognition and binding events. Recent studies examining hydration changes that accompany DNA binding by intercalators will be reviewed. Finally some aspects of cooperativity in drug-DNA interactions are described and the importance of considering cooperative effects when examining these reactions is highlighted.

Circular dichroism and thermal melting differentiation of Hoechst 33258 binding to the curved (A(4)T(4)) and straight (T(4)A(4)) DNA sequences

The ability of the B-DNA minor groove ligand Hoechst 33258 to discriminate between prototype curved and straight duplex DNA sequences was investigated by circular dichroism (CD) titrations at the wavelengths of absorbance of the ligand. The sequences were studied either within the framework of the ligated decamers (CA(4)T(4)G)(n) and (CT(4)A(4)G)(n), or within that of the single dodecamers GCA(4)T(4)GC and GCT(4)A(4)GC, to confirm and extend our earlier results based on fluorescence titrations of ligated decamers. A unique, strong binding site is invariantly present in both sequence units. The binding affinity of the drug for the site in the curved A(4)T(4) sequence was found 3- to 4-fold higher compared to the straight sequence. All these features hold true irrespective of the sequence framework, thus confirming that they reflect specific properties of the binding to the two sequences. Ligand binding increases the thermal stability of straight and curved duplex dodecamers to the same extent, thus maintaining the melting temperature differential between the two sequences. However, the different melting patterns and the difference between [total ligand]:[site] ratios needed for site saturation in the two duplexes are in agreement with the difference between binding constants derived from CD measurements.

Tris-benzimidazole derivatives: design, synthesis and DNA sequence recognition

Two tris-benzimidazole derivatives have been designed and synthesized based on the known structures of the bis-benzimidazole stain Hoechst 33258 complexed to short oligonucleotide duplexes derived from single crystal X-ray studies and from NMR. In both derivatives the phenol group has been replaced by a methoxy-phenyl substituent. Whereas one tris-benzimidazole carries a N-methyl-piperazine at the 6-position, the other one has this group replaced by a 2-amino-pyrrolidine ring. This latter substituent results in stronger DNA binding. The optimized synthesis of the drugs is described. The two tris-benzimidazoles exhibit high AT-base pair (bp) selectivity evident in footprinting experiments which show that five to six base pairs are protected by the tris-benzimidazoles as compared to four to five protected by the bis-benzimidazoles. The tris-benzimidazoles bind well to sequences like 5'-TAAAC, 5'-TTTAC and 5'-TTTAT, but it is also evident that they can bind weakly to sequences such as 5'-TATGTT-3' where the continuity of an AT stretch is interrupted by a single G*C base pair.

A simple, high-resolution method for establishing DNA binding affinity and sequence selectivity

Full details of the development of a simple, nondestructive, and high-throughput method for establishing DNA binding affinity and sequence selectivity are described. The method is based on the loss of fluorescence derived from the displacement of ethidium bromide or thiazole orange from the DNA of interest or, in selected instances, the change in intrinsic fluorescence of a DNA binding agent itself and is applicable for assessing relative or absolute DNA binding affinities. Enlisting a library of hairpin deoxyoligonucleotides containing all five base pair (512 hairpins) or four base pair (136 hairpins) sequences displayed in a 96-well format, a compound's rank order binding to all possible sequences is generated, resulting in a high-resolution definition of its sequence selectivity using this fluorescent intercalator displacement (FID) assay. As such, the technique complements the use of footprinting or affinity cleavage for the establishment of DNA binding selectivity and provides the information at a higher resolution. The merged bar graphs generated by this rank order binding provide a qualitative way to compare, or profile, DNA binding affinity and selectivity. The 96-well format assay (512 hairpins) can be conducted at a minimal cost (presently ca. $100 for hairpin deoxyoligonucleotides/assay with ethiduim bromide or less with thiazole orange), with a rapid readout using a fluorescent plate reader (15 min), and is adaptable to automation (Tecan Genesis Workstation 100 robotic system). Its use in generating a profile of DNA binding selectivity for several agents including distamycin A, netropsin, DAPI, Hoechst 33258, and berenil is described. Techniques for establishing binding constants from quantitative titrations are compared, and recommendations are made for use of a Scatchard or curve fitting analysis of the titration binding curves as a reliable means to quantitate the binding affinity.

Complex of B-DNA with polyamides freezes DNA backbone flexibility

The development of sequence-specific minor groove binding ligands is a modern and rapidly growing field of research because of their extraordinary importance as transcription-controlling drugs. We performed three molecular dynamics simulations in order to clarify the influence of minor groove binding of two ImHpPyPy-beta-Dp polyamides to the d(CCAGTACTGG)(2) decamer in the B-form. This decamer contains the recognition sequence for the trp repressor (5'-GTACT-3'), and it was investigated recently by X-ray crystallography. On one hand we are able to reproduce X-ray-determined DNA--drug contacts, and on the other hand we provide new contact information which is important for the development of potential ligands. The new insights show how the beta-tail of the polyamide ligands contributes to binding. Our simulations also indicate that complexation freezes the DNA backbone in a specific B(I) or B(II) substate conformation and thus optimizes nonbonded contacts. The existence of this distinct B(I)/B(II) substate pattern also allows the formation of water-mediated contacts. Thus, we suggest the B(I) <==> B(II) substate behavior to be an important part of the indirect readout of DNA.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)(2) and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar-phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving 'lock and key' recognition.

Structural constraints regulating triple helix formation by A-tracts

The study concerns the propensity of triple helix formation by different DNA oligonucleotides containing long A-tracts with and without flanking GxC base pairs in order to probe the role of length of the A-tract and the flanking sequences. From nuclear magnetic resonance (NMR) studies of imino proton spectra and circular dichroism (CD) spectroscopy of samples composed of potential triplex forming strand sequences in correct stoichiometries, we have concluded that 8-mer A-tracts flanked by GxC base pairs exert significant steric hindrance to triple helix formation. When as much as 50 mM Mg2+ was added, no triple helix formation was observed in these samples. In contrast, open-ended 8-mer A-tracts formed triplex with the corresponding two T8 strands under relatively mild ionic conditions (100 mM Na+). Moreover, the shorter the length of the A-tract, the less is the hindrance to form a triple helix.