DNA minor-groove recognition by small molecules

Sequence specific and high affinity recognition of 5'-ACGCGT-3' by rationally designed pyrrole-imidazole H-pin polyamides: thermodynamic and structural studies

Imidazole (Im) and Pyrrole (Py)-containing polyamides that can form stacked dimers can be programmed to target specific sequences in the minor groove of DNA and control gene expression. Even though various designs of polyamides have been thoroughly investigated for DNA sequence recognition, the use of H-pin polyamides (covalently cross-linked polyamides) has not received as much attention. Therefore, experiments were designed to systematically investigate the DNA recognition properties of two symmetrical H-pin polyamides composed of PyImPyIm (5) or f-ImPyIm (3e, f=formamido) tethered with an ethylene glycol linker. These compounds were created to recognize the cognate 5'-ACGCGT-3' through an overlapped and staggered binding motif, respectively. Results from DNaseI footprinting, thermal denaturation, circular dichroism, surface plasmon resonance and isothermal titration microcalorimetry studies demonstrated that both H-pin polyamides bound with higher affinity than their respective monomers. The binding affinity of formamido-containing H-pin 3e was more than a hundred times greater than that for the tetraamide H-pin 5, demonstrating the importance of having a formamido group and the staggered motif in enhancing affinity. However, compared to H-pin 3e, tetraamide H-pin 5 demonstrated superior binding preference for the cognate sequence over its non-cognates, ACCGGT and AAATTT. Data from SPR experiments yielded binding constants of 1.6x10(8)M(-1) and 2.0x10(10)M(-1) for PyImPyIm H-pin 5 and f-ImPyIm H-pin 3e, respectively. Both H-pins bound with significantly higher affinity (ca. 100-fold) than their corresponding unlinked PyImPyIm 4 and f-ImPyIm 2 counterparts. ITC analyses revealed modest enthalpies of reactions at 298 K (DeltaH of -3.3 and -1.0 kcal mol(-1) for 5 and 3e, respectively), indicating these were entropic-driven interactions. The heat capacities (DeltaC(p)) were determined to be -116 and -499 cal mol(-1)K(-1), respectively. These results are in general agreement with DeltaC(p) values determined from changes in the solvent accessible surface areas using complexes of the H-pins bound to (5'-CCACGCGTGG)(2). According to the models, the H-pins fit snugly in the minor groove and the linker comfortably holds both polyamide portions in place, with the oxygen atoms pointing into the solvent. In summary, the H-pin polyamide provides an important molecular design motif for the discovery of future generations of programmable small molecules capable of binding to target DNA sequences with high affinity and selectivity.

Antiparasitic compounds that target DNA

Designed, synthetic heterocyclic diamidines have excellent activity against eukaryotic parasites that cause diseases such as sleeping sickness and leishmania and adversely affect millions of people each year. The most active compounds bind specifically and strongly in the DNA minor groove at AT sequences. The compounds enter parasite cells rapidly and appear first in the kinetoplast that contains the mitochondrial DNA of the parasite. With time the compounds are also generally seen in the cell nucleus but are not significantly observed in the cytoplasm. The kinetoplast decays over time and disappears from the mitochondria of treated cells. At this point the compounds begin to be observed in other regions of the cell, such as the acidocalcisomes. The cells typically die in 24-48h after treatment. Active compounds appear to selectively target extended AT sequences and induce changes in kinetoplast DNA minicircles that cause a synergistic destruction of the catenated kinetoplast DNA network and cell death.

The usefulness of cyclic diamidines with different core-substituents as antitumor agents

A series of related polycationic compounds has been screened for potential antitumor activity by the NCI's in vitro testing (one dose primary anticancer assay and the NCI-60 full panel screening). The GI50 values of triazines 3 and 4 are on average 1.9 microM and 2.4 microM, respectively. Furan 8 deserves mention too (1.9 microM). The biological test results showed that carbazole 10 possessed cytotoxic activity in the nanomolar range, much better than the other compounds tested, only against several cancer cell lines: CCRF-CEM, HL-60(TB), MOLT-4, NCI-H522, COLO 205, SF-268, but the average GI50 value was higher (15 microM). The activity appears closely dependent on the core-shape and length of the bisimidazoline molecules (important for both high cytotoxicity and DNA binding). The mechanism of DNA minor-groove binding of diamidines 1-12, based on the anticancer parameters, is highly probable.

High mobility group A1 protein mediates human nitric oxide synthase 2 gene expression

Nitric oxide synthase (NOS)2, an inducible enzyme that produces NO during inflammation, is transcriptionally regulated. Our goal was to determine whether high mobility group (HMG)A1 contributes to human (h)NOS2 gene regulation. Using a small molecule inhibitor of HMGA1 binding to DNA, or a dominant-negative form of HMGA1, we blunted the induction of hNOS2 by pro-inflammatory stimuli. Binding of HMGA1 in the region -3506 to -3375 of the hNOS2 promoter, a region not previously known to be involved in hNOS2 regulation, contributed to the induction of hNOS2 promoter in conjunction with upstream enhancer regions. We demonstrate a previously unknown role for HMGA1 in the regulation of hNOS2.

Inhibition of high-mobility-group A2 protein binding to DNA by netropsin: a biosensor-surface plasmon resonance assay

The design of small synthetic molecules that can be used to affect gene expression is an area of active interest for development of agents in therapeutic and biotechnology applications. Many compounds that target the minor groove in AT sequences in DNA are well characterized and are promising reagents for use as modulators of protein-DNA complexes. The mammalian high-mobility-group transcriptional factor HMGA2 also targets the DNA minor groove and plays critical roles in disease processes from cancer to obesity. Biosensor-surface plasmon resonance methods were used to monitor HMGA2 binding to target sites on immobilized DNA, and a competition assay for inhibition of the HMGA2-DNA complex was designed. HMGA2 binds strongly to the DNA through AT hook domains with KD values of 20-40 nM depending on the DNA sequence. The well-characterized minor groove binder netropsin was used to develop and test the assay. The compound has two binding sites in the protein-DNA interaction sequence, and this provides an advantage for inhibition. An equation for analysis of results when the inhibitor has two binding sites in the biopolymer recognition surface is presented with the results. The assay provides a platform for discovery of HMGA2 inhibitors.

Circular dichroism to determine binding mode and affinity of ligand-DNA interactions

Circular dichroism (CD) is a useful technique for an assessment of DNA-binding mode, being a more accessible, low-resolution complement to NMR and X-ray diffraction methods. Ligand-DNA interactions can be studied by virtue of the interpretation of induced ligand CD signals resulting from the coupling of electric transition moments of the ligand and DNA bases within the asymmetric DNA environment. This protocol outlines methods to determine the binding mode and affinity of ligand-DNA interactions and takes approximately 7.5 h.

Microarray analysis of the in vivo sequence preferences of a minor groove binding drug

Background

Minor groove binding drugs (MGBDs) interact with DNA in a sequence-specific manner and can cause changes in gene expression at the level of transcription. They serve as valuable models for protein interactions with DNA and form an important class of antitumor, antiviral, antitrypanosomal and antibacterial drugs. There is a need to extend knowledge of the sequence requirements for MGBDs from in vitro DNA binding studies to living cells.

Results

Here we describe the use of microarray analysis to discover yeast genes that are affected by treatment with the MGBD berenil, thereby allowing the investigation of its sequence requirements for binding in vivo. A novel approach to sequence analysis allowed us to address hypotheses about genes that were directly or indirectly affected by drug binding. The results show that the sequence features of A/T richness and heteropolymeric character discovered by in vitro berenil binding studies are found upstream of genes hypothesized to be directly affected by berenil but not upstream of those hypothesized to be indirectly affected or those shown to be unaffected.

Conclusion

The data support the conclusion that effects of berenil on gene expression in yeast cells can be explained by sequence patterns discovered by in vitro binding experiments. The results shed light on the sequence and structural rules by which berenil binds to DNA and affects the transcriptional regulation of genes and contribute generally to the development of MGBDs as tools for basic and applied research.

Determination of chemically reduced pyrrolobenzodiazepine SJG-136 in human plasma by HPLC-MS/MS: application to an anticancer phase I dose escalation study

SJG-136 1,1'-[[(propane-1,3-diyl)dioxy]bis[(11aS)-7-methoxy-2-methylidene-1,2,3,11a-tetrahydro-5H-pyr- rolo[2,1-c][1,4]benzodiazepin-5-one]] (NSC 694501), is a bifunctional pyrrolobenzodiazepine (PBD) dimer that forms selective, irreversible, interstrand DNA cross-links via exocyclic N2 atoms of two guanine bases, with a preference for 5'PuGATCPy binding sites. SJG-136 is highly cytotoxic in human tumor cells in vitro and in human tumor xenograft models in vivo at subnanomolar concentrations and is currently in anticancer phase I clinical trials in the United Kingdom and United States. To support correlative pharmacokinetics studies, a highly sensitive HPLC-MS/MS assay was developed and validated for the reliable quantitation of SJG-136 in human plasma, using the structurally similar PBD dimer DSB-120 as an internal standard. Chemical reduction of SJG-136 to its corresponding amine (SJG-136-H(4), [M + H](+)m/z 561) improved HPLC peak resolution and sensitivity by minimizing complications that arose from the reactivity of the labile imine moieties. Plasma samples were processed by protein precipitation and centrifugal membrane dialysis; components were separated by HPLC using an Agilent Rapid Resolution HT 1.8 mm (2.1 mm x 50 mm) analytical column. The total analysis time from injection to injection was 11 min. Electrospray MS/MS detection of SJG-136-H(4) was based on the selected reaction monitoring (SRM) transition [M + H](+)m/z 561 --> 301. The analytical response ratio was linearly proportional to the plasma concentration of SJG-136 over the nominal concentration range of 25 pg/ml to 250 ng/ml, with a coefficient of determination of r > or = 0.999. The intrarun absolute %RE was < or =19.6, 14.2, and 14.0% at 0.056, 2.83, and 56.3 ng/ml, respectively. The corresponding %RSD was < or =14.9%, 9.01, and 4.59%. The interday %RSD was < or =2.72, 3.46, and 5.20%. The lower and upper limits of quantitation were 0.056 and 56 ng/ml, respectively; recovery of SJG-136 from plasma was > or = 62% across the validated concentration range. The sensitivity of the validated assay was sufficient to detect SJG-136 in human subjects for up to 6 h after intravenous administration of 6 microg/m(2), the starting dose of an NCI-sponsored dose escalation study.

Sequence and length dependent thermodynamic differences in heterocyclic diamidine interactions at AT base pairs in the DNA minor groove

With the goal of developing a better understanding of the antiparasitic biological action of DB75, we have evaluated its interaction with duplex alternating and nonalternating sequence AT polymers and oligomers. These DNAs provide an important pair of sequences in a detailed thermodynamic analysis of variations in interaction of DB75 with AT sites. The results for DB75 binding to the alternating and nonalternating AT sequences are quite different at the fundamental thermodynamic level. Although the Gibbs energies are similar, the enthalpies for DB75 binding with poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) are +3.1 and -4.5 kcal/mol, respectively, while the binding entropies are 41.7 and 15.2 cal/mol.K, respectively. The underlying thermodynamics of binding to AT sites in the minor groove plays a key role in the recognition process. It was also observed that DB75 binding with poly(dA).poly(dT) can induce T.A.T triplet formation and the compound binds strongly to the dT.dA.dT triplex.

Design of DNA minor groove binding diamidines that recognize GC base pair sequences: a dimeric-hinge interaction motif

The classical model of DNA minor groove binding compounds is that they should have a crescent shape that closely fits the helical twist of the groove. Several compounds with relatively linear shape and large dihedral twist, however, have been found recently to bind strongly to the minor groove. These observations raise the question of how far the curvature requirement could be relaxed. As an initial step in experimental analysis of this question, a linear triphenyl diamidine, DB1111, and a series of nitrogen tricyclic analogues were prepared. The goal with the heterocycles is to design GC binding selectivity into heterocyclic compounds that can get into cells and exert biological effects. The compounds have a zero radius of curvature from amidine carbon to amidine carbon but a significant dihedral twist across the tricyclic and amidine-ring junctions. They would not be expected to bind well to the DNA minor groove by shape-matching criteria. Detailed DNase I footprinting studies of the sequence specificity of this set of diamidines indicated that a pyrimidine heterocyclic derivative, DB1242, binds specifically to a GC-rich sequence, -GCTCG-. It binds to the GC sequence more strongly than to the usual AT recognition sequences for curved minor groove agents. Other similar derivatives did not exhibit the GC specificity. Biosensor-surface plasmon resonance and isothermal titration calorimetry experiments indicate that DB1242 binds to the GC sequence as a highly cooperative stacked dimer. Circular dichroism results indicate that the compound binds in the minor groove. Molecular modeling studies support a minor groove complex and provide an inter-compound and compound-DNA hydrogen-bonding rational for the unusual GC binding specificity and the requirement for a pyrimidine heterocycle. This compound represents a new direction in the development of DNA sequence-specific agents, and it is the first non-polyamide, synthetic compound to specifically recognize a DNA sequence with a majority of GC base pairs.

Metal-Based Netropsin Mimics Showing AT-Selective DNA Binding and DNA Cleavage Activity at Red Light

Copper(II) bis-arginate [Cu(l-arg)2](NO3)2 (1) and [Cu(l-arg)(phen)Cl]Cl (2) as mimics of the minor-groove-binding natural antibiotic netropsin show preferential binding to the AT-rich region of double-stranded DNA. The complexes with a d-d band near 600 nm display oxidative DNA cleavage activity on photoirradiation at UV-A light of 365 nm and at red light of 647.1 nm (Ar-Kr laser) in a metal-assisted photoexcitation process forming singlet oxygen (1O2) species in a type-2 pathway.

Spectroscopic studies of DNA binding modes of cation-substituted anthrapyrazoles derived from emodin

The DNA binding properties of three cation-substituted anthrapyrazole derivatives of emodin with calf thymus DNA were characterized by spectroscopic methods and the specific binding modes were elucidated. At low drug and high DNA concentrations, compound 1 with a mono-cationic amino side chain exhibited an intercalative binding mode, 2 with a much longer and more flexible di-cationic side chain exhibited an external binding mode, and 3 with a rigid di-cationic side chain exhibited both intercalative and external binding modes. The DNA binding mode of compounds was altered after structural modification. The molecular structure-DNA binding relationships found from this study may be useful for the design of anthrapyrazole derivatives with desired binding characteristics.

Strong and weak hydrogen bonds in drug-DNA complexes: A statistical analysis

A statistical analysis of strong and weak hydrogen bonds in the minor groove of DNA was carried out for a set of 70 drug-DNA complexes. The terms 'strong' and 'weak' pertain to the inherent strengths and weakness of the donor and acceptor fragments rather than to any energy considerations. The dataset was extracted from the protein data bank (PDB). The analysis was performed with an in-house software, hydrogen bond analysis tool (HBAT). In addition to strong hydrogen bonds such as O-H...O and N-H...O, the ubiquitous presence of weak hydrogen bonds such as C-H...O is implicated in molecular recognition. On an average, there are 1.4 weak hydrogen bonds for every strong hydrogen bond. For both categories of interaction, the N3 of purine and the O2 of pyrimidine are favoured acceptors. Donor multifurcation is common with the donors generally present in the drug molecules, and shared by hydrogen bond acceptors in the minor groove. Bifurcation and trifurcation are most commonly observed. The metrics for strong hydrogen bonds are consistent with established trends. The geometries are variable for weak hydrogen bonds. A database of recognition geometries for 26 literature amidinium-based inhibitors of Human African Trypanosomes (HAT) was generated with a docking study using seven inhibitors which occur in published crystal structures included in the list of 70 complexes mentioned above, and 19 inhibitors for which the drug-DNA complex crystal structures are unknown. The virtual geometries so generated correlate well with published activities for these 26 inhibitors, justifying our assumption that strong and weak hydrogen bonds are optimized in the active site.

Molecular recognition of DNA by small molecules: AT base pair specific intercalative binding of cytotoxic plant alkaloid palmatine

The base dependent binding of the cytotoxic alkaloid palmatine to four synthetic polynucleotides, poly(dA).poly(dT), poly(dA-dT).poly(dA-dT), poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) was examined by competition dialysis, spectrophotometric, spectrofluorimetric, thermal melting, circular dichroic, viscometric and isothermal titration calorimetric (ITC) studies. Binding of the alkaloid to various polynucleotides was dependent upon sequences of base pairs. Binding data obtained from absorbance measurements according to neighbour exclusion model indicated that the intrinsic binding constants decreased in the order poly(dA).poly(dT)>poly(dA-dT).poly(dA-dT)>poly(dG-dC).poly(dG-dC)>poly(dG).poly(dC). This affinity was also revealed by the competition dialysis, increase of steady state fluorescence intensity, increase in fluorescence quantum yield, stabilization against thermal denaturation and perturbations in circular dichroic spectrum. Among the polynucleotides, poly(dA).poly(dT) showed positive cooperativity at binding values lower than r=0.05. Viscosity studies revealed that in the strong binding region, the increase of contour length of DNA depended strongly on the sequence of base pairs being higher for AT polymers and induction of unwinding-rewinding process of covalently closed superhelical DNA. Isothermal titration calorimetric data showed a single entropy driven binding event in the AT homo polymer while that with the hetero polymer involved two binding modes, an entropy driven strong binding followed by an enthalpy driven weak binding. These results unequivocally established that the alkaloid palmatine binds strongly to AT homo and hetero polymers by mechanism of intercalation.

Unusually strong binding to the DNA minor groove by a highly twisted benzimidazole diphenylether: induced fit and bound water

RT29 is a dicationic diamidine derivative that does not obey the classical "rules" for shape and functional group placement that are expected to result in strong binding and specific recognition of the DNA minor groove. The compound contains a benzimidazole diphenyl ether core that is flanked by the amidine cations. The diphenyl ether is highly twisted and gives the entire compound too much curvature to fit well to the shape of the minor groove. DNase I footprinting, fluorescence intercalator displacement studies, and circular dichroism spectra, however, indicate that the compound is an AT specific minor groove binding agent. Even more surprisingly, quantitative biosensor-surface plasmon resonance and isothermal titration calorimetric results indicate that the compound binds with exceptional strength to certain AT sequences in DNA with a large negative enthalpy of binding. Crystallographic results for the DNA complex of RT29 compared to calculated results for the free compound show that the compound undergoes significant conformational changes to enhance its minor groove interactions. In addition, a water molecule is incorporated directly into the complex to complete the compound-DNA interface, and it forms an essential link between the compound and base pair edges at the floor of the minor groove. The calculated DeltaCp value for complex formation is substantially less than the experimentally observed value, which supports the idea of water being an intrinsic part of the complex with a major contribution to the DeltaCp value. Both the induced fit conformational changes of the compound and the bound water are essential for strong binding to DNA by RT29.

Capillary electrophoresis for studying drug–DNA interactions

The development of new drugs to treat disease by binding directly to DNA offers much promise but is reliant on methods to determine the relative affinity of the putative drug for different DNA sequences. Such methods should ideally be rapid and inexpensive as well as reliable. Use of capillary electrophoresis in simple silica columns offers such a method. The development of systems in which the solvent carries a soluble polymer allows the reliable separation of DNA oligomers, of 12-20 bp in length, which can then be titrated with the ligand in competition experiments. The results obtained are comparable with those obtained by footprinting and give direct graphical output, easily analysed for relative binding affinity.

Induced fit conformational changes of a "reversed amidine" heterocycle: optimized interactions in a DNA minor groove complex

To better understand the molecular basis for recognition of the DNA minor groove by heterocyclic cations, a series of "reversed amidine" substituted heterocycles has been prepared. Amidine derivatives for targeting the minor groove have the amidine carbon linked to a central heterocyclic system, whereas in the reverse orientation, an amidine nitrogen provides the link. The reverse system has a larger dihedral angle as well as a modified spatial relationship with the groove relative to amidines. Because of the large dihedral, the reversed amidines should have reduced binding to DNA relative to similar amidines. Such a reduction is observed in footprinting, circular dichroism (CD), biosensor-surface plasmon resonance (SPR), and isothermal titration calorimetric (ITC) experiments with DB613, which has a central phenyl-furan-phenyl heterocyclic system. The reduction is not seen when a pyrrole (DB884) is substituted for the furan. Analysis of a number of derivatives defines the pyrrole and a terminal phenyl substituent on the reversed amidine groups as critical components in the strong binding of DB884. ITC and SPR comparisons showed that the better binding of DB884 was due to a more favorable binding enthalpy and that it had exceptionally slow dissociation from DNA. Crystallographic analysis of DB884 bound to an AATT site shows that the compound was bound in the minor groove in a 1:1 complex as suggested by CD solution studies. Surprisingly, unlike the amidine derivative, the pyrrole -NH of DB884 formed an H-bond with a central T of the AATT site and this accounts for the enthalpy-driven strong binding. The structural results and molecular modeling studies provide an explanation for the differences in binding affinities for related amidine and reversed amidine analogues.

Amidine Analogues of Melphalan: Synthesis, Cytotoxic Activity, and DNA Binding Properties

Design, synthesis, and cytotoxic activity of amidine derivatives of melphalan are described and structure-activity relationships are discussed. Evaluation of the cytotoxicity of these compounds employing a MTT assay and inhibition of [(3)H]thymidine incorporation into DNA in both MDA-MB-231 and MCF-7 human breast cancer cells demonstrated that these compounds were more active than melphalan. Data from the ethidium displacement assay showed that these compounds were able to bind in the minor groove-binding mode in AT sequences of DNA. The cytotoxic properties of the amidine analogues of melphalan towards cultured human breast cancer cells correlate with topoisomerase II inhibitory properties but not with DNA-binding properties.

Hydration pattern of A4T4 and T4A4 DNA: A molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.

Pentamidine is an antiparasitic and apoptotic drug that selectively modifies ubiquitin

We have determined the cytotoxic properties of pentamidine isethionate (2) towards the promastigotes of the protozoan parasite Leishmania infantum. The leishmanicidal activity of 2 was 60 times higher after 72 h of incubation than that of cisplatin (4). The pentamidine salt 2 induced a higher amount of programmed cell death (PCD) than cisplatin, which is associated with inhibition of DNA synthesis and cell-cycle arrest in the G2/M phase. Circular dichroism (CD) data indicate that binding of 2 to calf-thymus DNA (CT-DNA) induces conformational changes in the DNA double helix, consistent with a B-->A transition. Moreover, the interaction of 2 with ubiquitin led to a 6% increase in the beta-sheet content of the protein as observed by CD spectroscopy. Fluorescence-spectroscopy studies agreed with the CD data, showing that the pentamidine portion of 2 induces a significant decrease in the fluorescence of the Ub residues Phe4 and Phe45 located on the beta-cluster of the molecule, but not of Tyr59 on the alpha-cluster. These data indicate that pentamidine specifically modifies the beta-cluster, i.e., the 'basic face' of ubiquitin. Our results suggest that the biochemical mechanism of action of pentamidine may be a consequence of its dual binding to DNA and proteins.

Break in the heat capacity change at 303 K for complex binding of netropsin to AATT containing hairpin DNA constructs

Studies performed in our laboratory demonstrated the formation of two thermodynamically distinct complexes on binding of netropsin to a number of hairpin-forming DNA sequences containing AATT-binding regions. These two complexes were proposed to differ only by a bridging water molecule between the drug and the DNA in the lower affinity complex. A temperature-dependent isothermal titration calorimetry (ITC)-binding study was performed using one of these constructs (a 20-mer hairpin of sequence 5'-CGAATTCGTCTCCGAATTCG) and netropsin. This study demonstrated a break in the heat capacity change for the formation of the complex containing the bridging water molecule at approximately 303 K. In the plot of the binding enthalpy change versus temperature, the slope (DeltaCp) was -0.67 kcal mol-1 K-1 steeper after the break at 303 K. Because of the relatively low melting temperature of the 20-mer hairpin (341 K (68 degrees C)), the enthalpy change for complex formation might have included some energy of refolding of the partially denatured hairpin, giving the suggestion of a larger DeltaCp. Studies done on the binding of netropsin to similar constructs, a 24-mer and a 28-mer, with added GC basepairs in the hairpin stem to increase thermal stability, exhibit the same nonlinearity in DeltaCp over the temperature range of from 275 to 333 K. The slopes (DeltaCp) were -0.69 and -0.64 kcal mol-1 K-1 steeper after 303 K for the 24-mer and 28-mer, respectively. This observation strengthens the argument regarding the presence of a bridging water molecule in the lower affinity netropsin/DNA complex. The DeltaCp data seem to infer that because the break in the heat capacity change function for the lower affinity binding occurs at the isoequilibrium temperature for water, water may be included or trapped in the complex. The fact that this break does not occur in the heat capacity change function for formation of the higher affinity complex can similarly be taken as evidence that water is not included in the higher affinity complex.

Hydration changes accompanying the binding of minor groove ligands with DNA

4',6-diamidino-2-phenylindole (DAPI), netropsin, and pentamidine are minor groove binders that have terminal -C(NH2)2+ groups. The hydration changes that accompany their binding to the minor groove of the (AATT)2 sequence have been studied using the osmotic stress technique with fluorescence spectroscopy. The affinity of DAPI for the binding site decreases with the increasing osmolality of the solution, resulting in acquisition of 35+/-1 waters upon binding. A competition fluorescence assay was utilized to measure the binding constants and hydration changes of the other two ligands, using the DNA-DAPI complex as the fluorescence reporter. Upon their association to the (AATT)2 binding site, netropsin and pentamidine acquire 26+/-3 and 34+/-2 additional waters of hydration, respectively. The hydration changes are discussed in the context of the terminal functional groups of the ligands and conformational changes in the DNA.

DNA sequence recognition in the minor groove by hairpin microgonotropens

Two novel microgonotropens (MGTs) comprised of hairpin N-propylaminepyrrole polyamides linked to a Hoechst 33258 (Ht) analogue (3 and 4) were synthesized on solid phase by adopting an Fmoc technique using a series of HOBt mediated coupling reactions. The dsDNA-binding properties of MGTs 3 and 4 were determined by thermal denaturation experiments. Both MGTs were found to be selective for their nine-bp match dsDNA sequence 9 and were less tolerant of G/C bp substitutions in the binding region than linear progenitor MGT 1. MGT 3 was intolerant of a G/C substitution located in the middle of the binding region and did not bind to sequences 13 and 14. MGT 4 also did not bind to sequence 13, and its linker-bound Ht moiety was found to be more sensitive to a G/C substitution in the Ht-binding target, as demonstrated by the lack of binding to sequence 16.

Rational selection of small molecules that increase transcription through the GAA repeats found in Friedreich's ataxia

Friedreich's ataxia (FRDA) is an autosomal recessive trinucleotide repeat disease with no effective therapy. Expanded GAA repeats in the first intron of the FRDA gene are thought to form unusual non-B DNA conformations that decrease transcription and subsequently reduce levels of the encoded protein, frataxin. Frataxin plays a crucial role in iron metabolism and detoxification. To discover small molecules that increase transcription through the GAA repeat region in FRDA, we have made stable cell lines containing a portion of expanded intron 1 fused to a GFP reporter. Small molecules identified using the competition dialysis method were found to increase FRDA-intron 1-reporter gene expression. One of these compounds, pentamidine, increases frataxin levels in patient cells. Thus our approach can be used to detect small molecules of potential therapeutic value in FRDA.

Contributions of Hydroxyethyl Groups to the DNA Binding Affinities of Anthracene Probes

Contributions of hydroxyethyl functions to the DNA binding affinities of substituted anthracenes are evaluated by calorimetry and spectroscopy. Isothermal titration calorimetry indicated that binding of the ligands to calf thymus DNA (5 mM Tris buffer, 50 mM NaCl, pH 7.2, 25 degrees C) is exothermic. The binding constants increased from 1.5 x 10(4) to 1.7 x 10(6) M(-1) as a function of increase in the number of hydroxyethyl functions (0-4). DNA binding was accompanied by red-shifted absorption (approximately 630 cm(-1)), strong hypochromism (>65%), positive induced-circular dichroism bands, and negative linear dichroism signals. DNA binding, in general, increased the helix stabilities to a significant extent (DeltaT(m) approximately 7 degrees C, DeltaDeltaH approximately 3 kcal/mol, DeltaDeltaS approximately 6-20 cal/K.mol). The binding constants showed a strong correlation with the number of hydroxyethyl groups present on the anthracene ring system. Analysis of the binding data using the hydrophobicity parameter (Log P) showed a poor correlation between the binding affinity and hydrophobicity. This observation was also supported by a comparison of the affinities of probes carrying N-ethyl (Kb = 0.8 x 10(5) M(-1)) versus N-hydroxyethyl side chains (Kb = 5.5 x 10(5) M(-1)). These are the very first examples of a strong quantitative correlation between the DNA binding affinity of a probe and the number of hydroxyethyl groups present on the probe. These quantitative findings are useful in the rational design of new ligands for high-affinity binding to DNA.

Virtual screening of DNA minor groove binders

The utility of the two docking programs DOCK and AutoDock in studying the binding of small molecules to the minor groove of B-DNA is examined. The AutoDock program is found to be more effective in both pose prediction and ranking of known binders over random compounds, and this superior performance is shown to be because of the scoring functions rather than the sampling algorithms.

A thermodynamic signature for drug–DNA binding mode

A number of small molecules bind directly and selectively to DNA, acting as chemotherapeutic agents by inhibiting replication, transcription or topoisomerase activity. Two common binding modes for these small molecules are intercalation or groove-binding. Intercalation results from insertion of a planar aromatic substituent between DNA base pairs, with concomitant unwinding and lengthening of the DNA helix. Groove binding, in contrast, does not perturb the duplex structure to any great extent. Groove-binders are typically crescent-shaped, and fit snugly into the minor groove with little distortion of the DNA structure. Recent calorimetric studies have determined the enthalpic and entropic contributions to the DNA binding of representative DNA binding compounds. Analysis of such thermodynamic data culled from the literature reveals distinctive thermodynamic signatures for groove-binding and intercalating compounds. Plots of the binding enthalpy (DeltaH) against binding entropy (-TDeltaS) for 26 drug-DNA interactions reveal that groove-binding interactions are clustered in a region of the graph with favorable entropy contributions to the free energy, while intercalators are clustered in a region with unfavorable entropy but favorable enthalpy contributions. Groove-binding is predominantly entropically driven, while intercalation in enthalpically driven. The molecular basis of the contrasting thermodynamic signatures for the two binding modes is by no means clear, but the pattern should be of use in categorizing new DNA binding agents.

DNA binding affinity of bisguanidine and bis(2-aminoimidazoline) derivatives with in vivo antitrypanosomal activity

A new antitrypanosomal hit compound that cures an acute (STIB 900) mouse model of Trypanosoma brucei rhodesiense trypanosomiasis is described. This bis(2-aminoimidazolinium) dicationic compound proved to be an excellent DNA minor groove binder, suggesting a possible mechanism for its trypanocidal activity. From these studies, the 4,4'-diaminodiphenylamine skeleton emerged as a good scaffold for antitrypanosomal drugs.

DNA sequence recognition by Hoechst 33258 conjugates of hairpin pyrrole/imidazole polyamides

A series of hairpin pyrrole/imidazole polyamides linked to a Hoechst 33258 (Ht) analogue (5-7) were synthesized on solid-phase by adopting an Fmoc technique using a series of PyBOP/HOBt mediated coupling reactions. The dsDNA binding properties of Ht-polyamides 5-7 were determined by thermal denaturation experiments. Hairpin Ht-polyamides 5-7 bound to dsDNA sequences 16 and 18 show DeltaTm values that are 14-18 degrees higher than linear Ht-polyamides bound to the same sequences. All three Ht-polyamides were found to be selective for their 9-bp match dsDNA sequences, supporting a relative stronger interaction of an Im/Py anti-parallel dimer with an appropriately positioned G/Cbp rather than sequences containing only A/Tbps. In addition, Ht-polyamides 5 and 7 showed a 20-fold preference for a properly placed G/Cbp over a C/Gbp, while 6 showed a 10-fold preference.

DNA sequence recognition by an imidazole-containing isopropyl-substituted thiazole polyamide (thiazotropsin B)

We have used DNA footprinting and fluorescence melting experiments to study the sequence specific binding of an imidazole-containing isopropyl-substituted thiazole polyamide (thiazotropsin B) to DNA. While the parent compound (thiazotropsin A) binds to the hexanucleotide sequence ACTAGT, changing one of the N-methylpyrrole groups to N-methylimidazole changes the preferred binding sequence to (A/T)CGCG(T/A). Experiments with DNA fragments that contain variants of this sequence suggest that the ligand can also bind, with lower affinity, to sequences which differ from this by 1bp in any position.

Binding of netropsin to several DNA constructs: evidence for at least two different 1:1 complexes formed from an -AATT-containing ds-DNA construct and a single minor groove binding ligand

Isothermal titration calorimetry, ITC, has been used to determine the thermodynamics (DeltaG, DeltaH, and -TDeltaS) for binding netropsin to a number of DNA constructs. The DNA constructs included: six different 20-22mer hairpin forming sequences and an 8-mer DNA forming a duplex dimer. All DNA constructs had a single -AT-rich netropsin binding with one of the following sequences, (A(2)T(2))(2), (ATAT)(2), or (AAAA/TTTT). Binding energetics are less dependent on site sequence than on changes in the neighboring single stranded DNA (hairpin loop size and tail length). All of the 1:1 complexes exhibit an enthalpy change that is dependent on the fractional saturation of the binding site. Later binding ligands interact with a significantly more favorable enthalpy change (partial differential DeltaH(1-2) from 2 to 6 kcal/mol) and a significantly less favorable entropy change (partial differential (-TDeltaS(1-2))) from -4 to -9 kcal/mol). The ITC data could only be fit within expected experimental error by use of a thermodynamic model that includes two independent binding processes with a combined stoichiometry of 1 mol of ligand per 1 mol of oligonucleotide. Based on the biophysical evidence reported here, including theoretical calculations for the energetics of "trapping" or structuring of a single water molecule and molecular docking computations, it is proposed that there are two modes by which flexible ligands can bind in the minor groove of duplex DNA. The higher affinity binding mode is for netropsin to lay along the floor of the minor groove in a bent conformation and exclude all water from the groove. The slightly weaker binding mode is for the netropsin molecule to have a slightly more linear conformation and for the required curvature to be the result of a water molecule that bridges between the floor of the minor groove and two of the amidino nitrogens located at one end of the bound netropsin molecule.

A high-throughput, high-resolution strategy for the study of site-selective DNA binding agents: analysis of a "highly twisted" benzimidazole-diamidine

A general strategy for the rapid structural analysis of DNA binding ligands is described as it was applied to the study of RT29, a benzimidazole-diamidine compound containing a highly twisted diphenyl ether linkage. By combining the existing high-throughput fluorescent intercalator displacement (HT-FID) assay developed by Boger et al. and a high-resolution (HR) host-guest crystallographic technique, a system was produced that was capable of determining detailed structural information pertaining to RT29-DNA interactions within approximately 3 days. Our application of the HT/HR strategy immediately revealed that RT29 has a preference for 4-base pair (bp), A.T-rich sites (AATT) and a similar tolerance and affinity for three A-T-bp sites (such as ATTC) containing a G.C bp. On the basis of these selectivities, oligonucleotides were designed and the host-guest crystallographic method was used to generate diffraction quality crystals. Analysis of the resulting crystal structures revealed that the diphenyl ether moiety of RT29 undergoes conformational changes that allow it to adopt a crescent shape that now complements the minor groove structure. The presence of a G.C bp in the RT29 binding site of ATTC did not overly perturb its interaction with DNA-the compound adjusted to the nucleobases that were available through water-mediated interactions. Our analyses suggest that the HT/HR strategy may be used to expedite the screening of novel minor groove binding compounds leading to a direct, HR structural determination.

Structure-specific recognition of quadruplex DNA by organic cations: influence of shape, substituents and charge

Combining structure-specific recognition of nucleic acids with limited sequence reading is a promising method to reduce the size of the recognition unit required to achieve the necessary selectivity and binding affinity to control function. It has been demonstrated recently that G-quadruplex DNA structures can be targeted by organic cations in a structure-specific manner. Structural targets of quadruplexes include the planar end surfaces of the G-tetrad stacked columns and four grooves. These provide different geometries and functional groups relative to duplex DNA. We have used surface plasmon resonance and isothermal titration calorimetry to show that binding affinity and selectivity of a series of quadruplex end-stacking molecules to human telomeric DNA are sensitive to compound shape as well as substituent type and position. ITC results indicate that binding is largely enthalpy driven. Circular dichroism was also used to identify a group of structurally related compounds that selectively target quadruplex grooves.

Spectroscopic identification of binding modes of anthracene probes and DNA sequence recognition

The binding properties of two anthracene derivatives with calf thymus DNA (CT DNA), poly(dA-dT), and poly(dG) x poly(dC) are reported. One contained bulky, cyclic cationic substituents at the 9 and 10 positions, and the other carried acylic, branched, cationic substituents. Binding of the probes to the DNA was examined by calorimetry, spectroscopy and helix melting studies. The cyclic derivative indicated exothermic binding, strong hypochromism, bathochromism, positive induced circular dichroism (CD, 300-400 nm), significant unwinding of the helix, large increases in the helix melting temperature, strong but negative linear dichroism (LD, 300-400 nm) and considerable stabilization of the helix. In contrast, the acyclic analog indicated thermoneutral binding, smaller hypochromism, no bathochromism, very weak induced CD, and no change in the helix melting temperature with any of the DNA polymers. A sharp distinction between the binding properties of the two probes is indicated, and both have intrinsic binding constants of approximately 10(6) M(-1) for the three polymers. However, when the ionic strength of the medium was lowered (10 mM NaCl), the absorption as well as CD spectral changes associated with the binding of the acyclic derivative corresponded with those of the cyclic derivative. The acyclic derivative showed large preference (10-fold) for poly(dG) x poly(dC) over poly(dA-dT), whereas the cyclic analog showed no preference. The characteristic spectroscopic signatures of the two distinct binding modes of these probes will be helpful in deciphering the interaction of other anthracene derivatives with DNA.

Novel amidine analogue of melphalan as a specific multifunctional inhibitor of growth and metabolism of human breast cancer cells

A novel amidine analogue of melphalan (AB4) was compared to its parent drug, melphalan in respect to cytotoxicity, DNA and collagen biosynthesis in MDA-MB-231 and MCF-7 human breast cancer cells. It was found that AB4 was more active inhibitor of DNA and collagen synthesis as well more cytotoxic agent than melphalan. The topoisomerase I/II inhibition assay indicated that AB4 is a potent catalytic inhibitor of topoisomerase II. Data from the ethidium displacement assay showed that AB4 intercalated into the minor-groove at AT sequences of DNA. The greater potency of AB4 to suppress collagen synthesis was found to be accompanied by a stronger inhibition of prolidase activity and expression compared to melphalan. The phenomenon was related to the inhibition of beta(1)-integrin and IGF-I receptor mediated signaling caused by AB4. The expression of beta(1)-integrin receptor, as well as Sos-1 and phosphorylated MAPK, ERK(1) and ERK(2) but not FAK, Shc, and Grb-2 was significantly decreased in cells incubated for 24h with 20 microM AB4 compared to the control, not treated cells, whereas in the same conditions melphalan did not evoke any changes in expression of all these signaling proteins, as shown by Western immunoblot analysis. These results indicate the amidine analogue of melphalan, AB4 represent multifunctional inhibitor of breast cancer cells growth and metabolism.

3D QSAR on a library of heterocyclic diamidine derivatives with antiparasitic activity

African trypanosomes, Trypanosoma brucei rhodesiense (TBR) and Trypanosoma brucei gambiense (TBG), affect hundreds of thousands of lives in tropical regions of the world. The toxicity of the diamidine pentamidine, an effective drug against TBG, necessitates the design of better drugs. An orally effective prodrug of the diamidine, furamidine (DB75), presently scheduled for phase III clinical trials, has excellent activity against TBG with toxicity lower than that of pentamidine. As part of an effort to develop additional and improved diamidines against African trypanosomes, CoMFA and CoMSIA 3D QSAR analyses have been conducted with furamidine and a set of 25 other structurally related compounds. Two different alignment strategies, based on a putative kinetoplast DNA minor groove target, were used. Due to conserved electrostatic properties across the compounds, models that used only steric and electronic properties did not perform well in predicting biological results. An extended CoMSIA model with additional descriptors for hydrophobic, donor, and acceptor properties had good predictive ability with a q2=0.699, r2=0.974, SEE, standard error of estimate=0.1, and F=120.04. The results have been used as a guide to design compounds that, potentially, have better activity against African trypanosomes.

Diastereoselective DNA cleavage recognition by Ni(II) x Gly-Gly-His-derived metallopeptides

Site-selective DNA cleavage by diastereoisomers of Ni(II) x Gly-Gly-His-derived metallopeptides was investigated through high-resolution gel analyses and molecular dynamics simulations. Ni(II) x L-Arg-Gly-His and Ni(II) x D-Arg-Gly-His (and their respective Lys analogues) targeted A/T-rich regions; however, the L-isomers consistently modified a subset of available nucleotides within a given minor groove site, while the D-isomers differed in both their sites of preference and their ability to target individual nucleotides within some sites. In comparison, Ni(II) x L-Pro-Gly-His and Ni(II) x D-Pro-Gly-His were unable to exhibit a similar diastereoselectivity. Simulations of the above systems, along with Ni(II) x Gly-Gly-His, indicated that the stereochemistry of the amino-terminal amino acid produces either an isohelical metallopeptide that associates stably at individual DNA sites (L-Arg or L-Lys) or, with D-Arg and D-Lys, a noncomplementary metallopeptide structure that cannot fully employ its side chain nor amino-terminal amine as positional stabilizing moieties. In contrast, amino-terminal Pro-containing metallopeptides of either stereochemistry, lacking an extended side chain directed toward the minor groove, did not exhibit a similar diastereoselectivity. While the identity and stereochemistry of amino acids located in the amino-terminal peptide position influenced DNA cleavage, metallopeptide diastereoisomers containing L- and D-Arg (or Lys) within the second peptide position did not exhibit diastereoselective DNA cleavage patterns; simulations indicated that a positively charged amino acid in this location alters the interaction of the metallopeptide equatorial plane and the minor groove leading to an interaction similar to Ni(II) x Gly-Gly-His.