A thermodynamic and structural analysis of DNA minor-groove complex formation

Copper(II) complexes containing hydrazone and bipyridine/phenanthroline ligands for anticancer application against breast cancer cells

In this work, mixed ligand Cu(II) complexes containing hydrazone and bipyridine ligands (CB1-CB5), or hydrazone and phenanthroline ligands (CP1-CP5) have been synthesized and characterized by spectroscopic and analytical techniques. Single crystal X-ray structure of complex CB1 revealed that two nitrogen atoms from bipyridine, one carbonyl oxygen, one azomethine nitrogen and one hydroxyl oxygen from the hydrazone ligand coordinated to Cu(II) ion, adopting a distorted square pyramidal geometry. Interaction of these complexes with calf thymus (CT) DNA and bovine serum albumin (BSA) was analyzed by absorption and emission studies. Further, the in vitro anticancer activity of the complexes was investigated exclusively against the breast cancer cells namely MCF7, T47D and MDA MB 231, and a normal breast MCF 10a cell line. The phenanthroline bearing complexes (CP1-CP5) displayed better activity than their bipyridine counterparts as seen from the IC50 values. In addition, the most active complex CP1 having an IC50 value of 5.8 ± 0.3 μM against T47D cancer cells was investigated for its mode of cell death through acridine orange/ethidium bromide(AO/EB), 4',6-diamidino-2-phenylindole (DAPI) and Annexin-V fluorescein isothiocyanate (FITC) staining assays which revealed apoptosis. Lastly, the cell cycle analysis revealed that complex CP1 induced cell death in T47D cancer cells at the G0/G1 phase.

Novel Amidine Derivative K1586 Sensitizes Colorectal Cancer Cells to Ionizing Radiation by Inducing Chk1 Instability

Checkpoint kinase 1 (Chk1) is a key mediator of the DNA damage response that regulates cell cycle progression, DNA damage repair, and DNA replication. Small-molecule Chk1 inhibitors sensitize cancer cells to genotoxic agents and have shown preclinical activity as single agents in cancers characterized by high levels of replication stress. However, the underlying genetic determinants of Chk1-inhibitor sensitivity remain unclear. Although treatment options for advanced colorectal cancer are limited, radiotherapy is effective. Here, we report that exposure to a novel amidine derivative, K1586, leads to an initial reduction in the proliferative potential of colorectal cancer cells. Cell cycle analysis revealed that the length of the G2/M phase increased with K1586 exposure as a result of Chk1 instability. Exposure to K1586 enhanced the degradation of Chk1 in a time- and dose-dependent manner, increasing replication stress and sensitizing colorectal cancer cells to radiation. Taken together, the results suggest that a novel amidine derivative may have potential as a radiotherapy-sensitization agent that targets Chk1.

Drug design and DNA structural research inspired by the Neidle laboratory: DNA minor groove binding and transcription factor inhibition by thiophene diamidines

The understanding of sequence-specific DNA minor groove interactions has recently made major steps forward and as a result, the goal of development of compounds that target the minor groove is an active research area. In an effort to develop biologically active minor groove agents, we are preparing and exploring the DNA interactions of diverse diamidine derivatives with a 5'-GAATTC-3' binding site using a powerful array of methods including, biosensor-SPR methods, and X-ray crystallography. The benzimidazole-thiophene module provides an excellent minor groove recognition component. A central thiophene in a benzimidazole-thiophene-phenyl aromatic system provides essentially optimum curvature for matching the shape of the minor groove. Comparison of that structure to one with the benzimidazole replaced with an indole shows that the two structures are very similar, but have some interesting and important differences in electrostatic potential maps, the DNA minor groove binding structure based on x-ray crystallographic analysis, and inhibition of the major groove binding PU.1 transcription factor complex. The binding KD for both compounds is under 10 nM and both form amidine H-bonds to DNA bases. They both have bifurcated H-bonds from the benzimidazole or indole groups to bases at the center of the -AATT- binding site. Analysis of the comparative results provides an excellent understanding of how thiophene compounds recognize the minor groove and can act as transcription factor inhibitors.

Biological Activity of Newly Synthesized Benzimidazole and Benzothizole 2,5-Disubstituted Furane Derivatives

Newly designed and synthesized cyano, amidino and acrylonitrile 2,5-disubstituted furane derivatives with either benzimidazole/benzothiazole nuclei have been evaluated for antitumor and antimicrobial activity. For potential antitumor activity, the compounds were tested in 2D and 3D cell culture methods on three human lung cancer cell lines, A549, HCC827 and NCI-H358, with MTS cytotoxicity and BrdU proliferation assays in vitro. Compounds 5, 6, 8, 9 and 15 have been proven to be compounds with potential antitumor activity with high potential to stop the proliferation of cells. In general, benzothiazole derivatives were more active in comparison to benzimidazole derivatives. Antimicrobial activity was evaluated with Broth microdilution testing (according to CLSI (Clinical Laboratory Standards Institute) guidelines) on Gram-negative Escherichia coli and Gram-positive Staphylococcus aureus. Additionally, Saccharomyces cerevisiae was included in testing as a eukaryotic model organism. Compounds 5, 6, 8, 9 and 15 showed the most promising antibacterial activity. In general, the compounds showed antitumor activity, higher in 2D assays in comparison with 3D assays, on all three cell lines in both assays. In natural conditions, compounds with such an activity profile (less toxic but still effective against tumor growth) could be promising new antitumor drugs. Some of the tested compounds showed antimicrobial activity. In contrast to ctDNA, the presence of nitro group or chlorine in selected furane-benzothiazole structures did not influence the binding mode with AT-DNA. All compounds dominantly bound inside the minor groove of AT-DNA either in form of monomers or dimer and higher-order aggregates.

Biosensor-Surface Plasmon Resonance: Label-Free Method for Investigation of Small Molecule-Quadruplex Nucleic Acid Interactions

Biosensor-surface plasmon resonance (SPR) technology is now well established as a quantitative approach for the study of nucleic acid interactions in real time, without the need for labeling any components of the interaction. The method provides real-time equilibrium and kinetic characterization for quadruplex DNA interactions and requires small amounts of materials and no external probe. A detailed protocol for quadruplex-DNA interaction analyses with a variety of binding molecules using biosensor-SPR methods is presented. Explanations of the SPR method with basic fundamentals for use and analysis of results are described with recommendations on the preparation of the SPR instrument, sensor chips, and samples. Details of experimental design, quantitative and qualitative data analyses, and presentation are described. Some specific examples of small molecule-DNA quadruplex interactions are presented with results evaluated by both kinetic and steady-state SPR methods.

Small molecule binders recognize DNA microstructural variations via an induced fit mechanism

Regulation of gene-expression by specific targeting of protein-nucleic acid interactions has been a long-standing goal in medicinal chemistry. Transcription factors are considered "undruggable" because they lack binding sites well suited for binding small-molecules. In order to overcome this obstacle, we are interested in designing small molecules that bind to the corresponding promoter sequences and either prevent or modulate transcription factor association via an allosteric mechanism. To achieve this, we must design small molecules that are both sequence-specific and able to target G/C base pair sites. A thorough understanding of the relationship between binding affinity and the structural aspects of the small molecule-DNA complex would greatly aid in rational design of such compounds. Here we present a comprehensive analysis of sequence-specific DNA association of a synthetic minor groove binder using long timescale molecular dynamics. We show how binding selectivity arises from a combination of structural factors. Our results provide a framework for the rational design and optimization of synthetic small molecules in order to improve site-specific targeting of DNA for therapeutic uses in the design of selective DNA binders targeting transcription regulation.

Intermolecular interactions between DNA and methamphetamine, amphetamine, ecstasy and their major metabolites

In this work, we carried out a theoretical investigation regarding amphetamine-type stimulants, which can cause central nervous system degeneration, interacting with human DNA. These include amphetamine, methamphetamine, 3,4-Methylenedioxymethamphetamine (also known as ecstasy), as well as their main metabolites. The studies were performed through molecular docking and molecular dynamics simulations, where molecular interactions of the receptor-ligand systems, along with their physical-chemical energies, were reported. Our results show that 3,4-Methylenedioxymethamphetamine and 3,4-Dihydroxymethamphetamine (ecstasy) present considerable reactivity with the receptor (DNA), suggesting that these molecules may cause damage due to human-DNA. These results were indicated by free Gibbs change of bind (ΔG_bind) values referring to intermolecular interactions between the drugs and the minor grooves of DNA, which were predominant for all simulations. In addition, it was observed that 3,4-Dihydroxymethamphetamine (ΔG_bind = -13.15 kcal/mol) presented greater spontaneity in establishing interactions with DNA in comparison to 3,4-Methylenedioxymethamphetamine (ΔG_bind = -8.61 kcal/mol). Thus, according with the calculations performed our results suggest that the 3,4-Methylenedioxymethamphetamine and 3,4-Dihydroxymethamphetamine have greater probability to provide damage to human DNA fragments.

Shape readout of AT-rich DNA by carbohydrates

Gene expression can be altered by small molecules that target DNA; sequence as well as shape selectivities are both extremely important for DNA recognition by intercalating and groove-binding ligands. We have characterized a carbohydrate scaffold (1) exhibiting DNA "shape readout" properties. Thermodynamic studies with 1 and model duplex DNAs demonstrate the molecule's high affinity and selectivity towards B* form (continuous AT-rich) DNA. Isothermal titration calorimetry (ITC), circular dichroism (CD) titration, ultraviolet (UV) thermal denaturation, and Differential Scanning Calorimetry were used to characterize the binding of 1 with a B* form AT-rich DNA duplex d[5'-G2 A6 T6 C2 -3']. The binding constant was determined using ITC at various temperatures, salt concentrations, and pH. ITC titrations were fit using a two-binding site model. The first binding event was shown to have a 1:1 binding stoichiometry and was predominantly entropy-driven with a binding constant of approximately 10(8) M(-1) . ITC-derived binding enthalpies were used to obtain the binding-induced change in heat capacity (ΔCp ) of -225 ± 19 cal/mol·K. The ionic strength dependence of the binding constant indicated a significant electrolytic contribution in ligand:DNA binding, with approximately four to five ion pairs involved in binding. Ligand 1 displayed a significantly higher affinity towards AT-tract DNA over sequences containing GC inserts, and binding experiments revealed the order of binding affinity for 1 with DNA duplexes: contiguous B* form AT-rich DNA (d[5'-G2 A6 T6 C2 -3']) >B form alternate AT-rich DNA (d[5'-G2 (AT)6 C2- 3']) > A form GC-rich DNA (d[5'-A2 G6 C6 T2 -3']), demonstrating the preference of ligand 1 for B* form DNA.

Validation of a computational docking methodology to identify the non-covalent binding site of ligands to DNA

Despite the biomedical consequences of carcinogen-DNA interactions and the potential of DNA as a drug target in medicinal chemistry, only a small number of studies have validated or used docking methods for the prediction of the physical binding of small molecules to DNA. Knowledge of the DNA-physically-bound ligand geometry can lead to the elucidation of the molecular-level mechanism of drugs as well as predicting the subsequent chemical interactions that lead to DNA damage from carcinogens. We sought to validate AutoDock 4.2, a docking method that includes a physics-based free energy function and a Lamarckian Genetic Algorithm, for the prediction of ligand geometries upon physical binding to DNA. We performed simulations by systematically changing the length of the search process for a comprehensive set of 32 ligand-DNA molecular systems with different physico-chemical properties, and we used a free-energy-based convergence criterion to terminate our simulations. For 11 out of 28 molecular systems for which convergence was achieved, the lowest binding free energy geometries were within 2 Å of the experimentally determined geometry. Considering all predicted sites with free energy changes within 20% of the lowest binding free energy site, we found a site within 2 Å of the experimentally determined geometry for 24 out of the 28 systems. However, the predicted hydrogen bonding interactions were different for most molecular systems compared to the same interactions in the experimentally determined geometry. We discuss reasons for the successes and failures, implications, and the importance of ensuring an adequate search in docking calculations. Overall, we concluded that AutoDock 4.2 can be used to predict the non-covalent binding geometry of a small molecule to DNA with some limitations.

Bisimidazoline arylamides binding to the DNA minor groove: N1-hydroxylation enhances binding affinity and selectivity to AATT sites

Selective binding of N-hydroxy bisimidazolines to dsDNA GCAATTGC is derived from a tighter fit to this narrower minor groove.

Sensitivities to parameterization in the size-modified Poisson-Boltzmann equation

Experimental results have demonstrated that the numbers of counterions surrounding nucleic acids differ from those predicted by the nonlinear Poisson-Boltzmann equation, NLPBE. Some studies have fit these data against the ion size in the size-modified Poisson-Boltzmann equation, SMPBE, but the present study demonstrates that other parameters, such as the Stern layer thickness and the molecular surface definition, can change the number of bound ions by amounts comparable to varying the ion size. These parameters will therefore have to be fit simultaneously against experimental data. In addition, the data presented here demonstrate that the derivative, SK, of the electrostatic binding free energy, ΔGel, with respect to the logarithm of the salt concentration is sensitive to these parameters, and experimental measurements of SK could be used to parameterize the model. However, although better values for the Stern layer thickness and ion size and better molecular surface definitions could improve the model's predictions of the numbers of ions around biomolecules and SK, ΔGel itself is more sensitive to parameters, such as the interior dielectric constant, which in turn do not significantly affect the distributions of ions around biomolecules. Therefore, improved estimates of the ion size and Stern layer thickness to use in the SMPBE will not necessarily improve the model's predictions of ΔGel.

Anticancer, antioxidant activities, and DNA affinity of novel monocationic bithiophenes and analogues

A series of 15 monocationic bithiophenes and isosteres were prepared and subjected to in vitro antiproliferative screening using the full National Cancer Institute (NCI)-60 cell line panel, representing nine types of cancer. Among the nine types of cancer involved in a five-dose screen, non-small cell lung and breast cancer cell lines were the most responsive to the antiproliferative effect of the tested compounds, especially cell lines A549/ATCC, NCI-H322M, and NCI-H460, whereas compounds 1a, 1c, 1d, and 7 exhibited potent activity, with GI50 values (drug concentration that causes 50% inhibition of cell growth) from less than 10 nM to 102 nM. In addition, compounds 1c and 1d gave GI50 values of 73 nM and 79 nM, respectively, against the MDA-MB-468 breast cancer cell line. Structure-activity relationship findings indicated that the mononitriles were far less active than their corresponding monoamidines and, within the amidines series, the bioisosteric replacement of a thiophene ring by a furan led to a reduction in antiproliferative activity. Also, molecular manipulations, involving substitution on the phenyl ring, or its replacement by a pyridyl, or alteration of the position of the amidine group, led to significant alteration in antiproliferative activity. On the other hand, DNA studies demonstrated that these monoamidine bichalcophenes have promising ability to cleave the genomic DNA. These monoamidines show a wide range of DNA affinities, as judged from their DNA cleavage effect, which are remarkably sensitive to all kinds of structural modifications. Finally, the novel bichalcophenes were tested for their antioxidant property by the ABTS (2,2'-azino- bis(3-ethylbenzthiazoline-6-sulfonic acid) diammonium salt) assay, as well as lipid and nitric oxide scavenging techniques, and were found to exhibit good-to-potent antioxidant abilities.

Different thermodynamic signatures for DNA minor groove binding with changes in salt concentration and temperature

The effects of salt concentration and temperature on the thermodynamics of DNA minor groove binding have quite different signatures: binding enthalpy is salt concentration independent but temperature dependent. Conversely, binding free energy is salt dependent but essentially temperature independent through enthalpy-entropy compensation.

Design, synthesis, and DNA binding characteristics of a group of orthogonally positioned diamino, N-formamido, pyrrole- and imidazole-containing polyamides

Orthogonally positioned diamino/dicationic polyamides (PAs) have good water solubility and enhanced binding affinity, whilst retaining DNA minor groove and sequence specificity compared to their monoamino/monocationic counterparts. The synthesis and DNA binding properties of the following diamino PAs: f-IPI (3a), f-IPP (4), f-PIP (5), and f-PPP (6) are described. P denotes the site where a 1-propylamino group is attached to the N1-position of the heterocycle. Binding of the diamino PAs to DNA was assessed by DNase I footprinting, thermal denaturation, circular dichroism titration, biosensor surface plasmon resonance (SPR), and isothermal titration calorimetry (ITC) studies. According to SPR studies, f-IPI (3a) bound more strongly (K(eq)=2.4×10(8) M(-1)) and with comparable sequence selectivity to its cognate sequence 5'-ACGCGT-3' when compared to its monoamino analog f-IPI (1). The binding of f-IPI (3a) to 5'-ACGCGT-3' via the stacked dimer motif was balanced between enthalpy and entropy, and that was quite different from the enthalpy-driven binding of its monoamino parent f-IPI (1). f-IPP (4) also bound more strongly to its cognate sequence 5'-ATGCAT-3' (K(eq)=7.4×10(6) M(-1)) via the side-by-side stacked motif than its monoamino analog f-IPP (2a). Although f-PPP (6) bound via a 1:1 motif, it bound strongly to its cognate sequence 5'-AAATTT-3' (K(eq)=4.8×10(7) M(-1)), 15-times higher than the binding of its monoamino analog f-PPP (2c), albeit f-PPP bound via the stacked motif. Finally, f-PIP (5) bound to its target sequence 5'-ATCGAT-3' as a stacked dimer and it has the lowest affinity among the diamino PAs tested (Keq <1×10(5) M(-1)). This was about two times lower in affinity than the binding of its monoamino analog f-PIP (2b). The results further demonstrated that the 'core rules' of DNA recognition by monoamino PAs also apply to their diamino analogs. Specifically, PAs that contain a stacked IP core structure bind most strongly (highest binding constants) to their cognate GC doublet, followed by the binding of PAs with a stacked PP structure to two degenerate AT base pairs, and finally the binding of PAs with a PI core to their cognate CG doublet.

Effect of G-quadruplex polymorphism on the recognition of telomeric DNA by a metal complex

The physiological role(s) played by G-quadruplexes renders these 'non-canonical' DNA secondary structures interesting new targets for therapeutic intervention. In particular, the search for ligands for selective recognition and stabilization of G-quadruplex arrangements has led to a number of novel targeted agents. An interesting approach is represented by the use of metal-complexes, their binding to DNA being modulated by ligand and metal ion nature, and by complex stoichiometry. In this work we characterized thermodynamically and stereochemically the interactions of a Ni(II) bis-phenanthroline derivative with telomeric G-quadruplex sequences using calorimetric, chiroptical and NMR techniques. We employed three strictly related sequences based on the human telomeric repeat, namely Tel22, Tel26 and wtTel26, which assume distinct conformations in potassium containing solutions. We were able to monitor specific enthalpy/entropy changes according to the structural features of the target telomeric sequence and to dissect the binding process into distinct events. Interestingly, temperature effects turned out to be prominent both in terms of binding stoichiometry and ΔH/ΔS contributions, while the final G-quadruplex-metal complex architecture tended to merge for the examined sequences. These results underline the critical choice of experimental conditions and DNA sequence for practical use of thermodynamic data in the rational development of effective G-quadruplex binders.

Correlation of local effects of DNA sequence and position of β-alanine inserts with polyamide-DNA complex binding affinities and kinetics

To improve our understanding of the effects of β-alanine (β) substitution and the number of heterocycles on DNA binding affinity and selectivity, we investigated the interactions of an eight-ring hairpin polyamide (PA) and two β derivatives as well as a six-heterocycle analogue with their cognate DNA sequence, 5'-TGGCTT-3'. Binding selectivity and the effects of β have been investigated with the cognate and five mutant DNAs. A set of powerful and complementary methods have been employed for both energetic and structural evaluations: UV melting, biosensor surface plasmon resonance, isothermal titration calorimetry, circular dichroism, and a DNA ligation ladder global structure assay. The reduced number of heterocycles in the six-ring PA weakens the binding affinity; however, the smaller PA aggregates significantly less than the larger PAs and allows us to obtain the binding thermodynamics. The PA-DNA binding enthalpy is large and negative with a large negative ΔC(p) and is the primary driving component of the Gibbs free energy. The complete SPR binding results clearly show that β substitutions can substantially weaken the binding affinity of hairpin PAs in a position-dependent manner. More importantly, the changes in the binding of PA to the mutant DNAs further confirm the position-dependent effects on the PA-DNA interaction affinity. Comparison of mutant DNA sequences also shows a different effect in recognition of T·A versus A·T base pairs. The effects of DNA mutations on binding of a single PA as well as the effects of the position of β substitution on binding tell a clear and very important story about sequence-dependent binding of PAs to DNA.

Multi-spectroscopic methods combined with molecular modeling dissect the interaction mechanisms of ractopamine and calf thymus DNA

The toxic interaction of ractopamine (RAC) with calf thymus DNA (ct DNA) was studied in vitro using multi-spectroscopic methods and molecular modeling methods. The hypochromic effect without a noticeable shift in UV-vis absorption indicated that the minor groove binding mode existed in the interaction between RAC and DNA. The fluorescence quenching of RAC was observed with the increasing addition of DNA and was proved to be the static quenching. The binding constant and the binding site sizes were 4.13 × 10(3) and 0.97, respectively. The thermodynamic calculation demonstrated that the hydrogen bond and van der Waals were main acting forces. This result further confirmed the existence of groove binding mode. Afterwards, we found another interaction mode, electrostatic binding mode through the fluorescence polarization, ionic effects and denatured DNA experiments. Circular dichroism spectroscopy (CD) was then employed to monitor the conformation changes of DNA. Molecular modeling studies illustrated the visual display of the binding mode and the detailed information of the H-bond.

The efficacy of novel arylimidamides against Trypanosoma cruzi in vitro

The present study aimed to determine the in vitro biological efficacy and selectivity of 7 novel AIAs upon bloodstream trypomastigotes and intracellular amastigotes of Trypanosoma cruzi. The biological activity of these aromatic compounds was assayed for 48 and 24 h against intracellular parasites and bloodstream forms of T. cruzi (Y strain), respectively. Additional assays were also performed to determine their potential use in blood banks by treating the bloodstream parasites with the compounds diluted in mouse blood for 24 h at 4°C. Toxicity against mammalian cells was evaluated using primary cultures of cardiac cells incubated for 24 and 48 h with the AIAs and then cellular death rates were determined by MTT colorimetric assays. Our data demonstrated the outstanding trypanocidal effect of AIAs against T. cruzi, especially DB1853, DB1862, DB1867 and DB1868, giving IC50 values ranging between 16 and 70 nanomolar against both parasite forms. All AIAs presented superior efficacy to benznidazole and some, such as DB1868, also demonstrated promising activity as a candidate agent for blood prophylaxis. The excellent anti-trypanosomal efficacy of these novel AIAs against T. cruzi stimulates further in vivo studies and justifies the screening of new analogues with the goal of establishing a useful alternative therapy for Chagas disease.

Specificity of RSG-1.2 peptide binding to RRE-IIB RNA element of HIV-1 over Rev peptide is mainly enthalpic in origin

Rev is an essential HIV-1 regulatory protein which binds to the Rev responsive element (RRE) present within the env gene of HIV-1 RNA genome. This binding facilitates the transport of the RNA to the cytoplasm, which in turn triggers the switch between viral latency and active viral replication. Essential components of this complex have been localized to a minimal arginine rich Rev peptide and stem IIB region of RRE. A synthetic peptide known as RSG-1.2 binds with high binding affinity and specificity to the RRE-IIB than the Rev peptide, however the thermodynamic basis of this specificity has not yet been addressed. The present study aims to probe the thermodynamic origin of this specificity of RSG-1.2 over Rev Peptide for RRE-IIB. The temperature dependent melting studies show that RSG-1.2 binding stabilizes the RRE structure significantly (ΔT_m = 4.3°C), in contrast to Rev binding. Interestingly the thermodynamic signatures of the binding have also been found to be different for both the peptides. At pH 7.5, RSG-1.2 binds RRE-IIB with a K_a = 16.2±0.6×10⁷ M⁻¹ where enthalpic change ΔH = −13.9±0.1 kcal/mol is the main driving force with limited unfavorable contribution from entropic change TΔS = −2.8±0.1 kcal/mol. A large part of ΔH may be due to specific stacking between U72 and Arg15. In contrast binding of Rev (K_a = 3.1±0.4×10⁷ M⁻¹) is driven mainly by entropy (ΔH = 0 kcal/mol and TΔS = 10.2±0.2 kcal/mol) which arises from major conformational changes in the RNA upon binding.

Water-mediated binding of agents that target the DNA minor groove

Small molecule complexes with DNA that incorporate linking water molecules are rare, and the DB921-DNA complex has provided a unique and well-defined system for analysis of water-mediated binding in the context of a DNA complex. DB921 has a benzimidazole-biphenyl system with terminal amidines that results in a linear conformation that does not possess the appropriate radius of curvature to match the minor groove shape and represents a new paradigm that does not fit the classical model of minor groove interactions. To better understand the role of the bound water molecule observed in the X-ray crystal structure of the DB921 complex, synthetic modifications have been made in the DB921 structure, and the interactions of the new compounds with DNA AT sites have been evaluated with an array of methods, including DNase I footprinting, biosensor-surface plasmon resonance, isothermal titration microcalorimetry, and circular dichroism. The interaction of a key compound, which has the amidine at the phenyl shifted from the para position in DB921 to the meta position, has also been examined by X-ray crystallography. The detailed structural, thermodynamic, and kinetic results provide valuable new information for incorporation of water molecules in the design of new lead scaffolds for targeting DNA in chemical biology and therapeutic applications.

Interaction of minor groove ligands with G-quadruplexes: thermodynamic contributions of the number of quartets, T-U substitutions, and conformation

In the presence of specific metal ions, DNA oligonucleotides containing guanine repeat sequences can adopt G-quadruplex structures. In this work, we used a combination of spectroscopic and calorimetric techniques to investigate the conformation and unfolding thermodynamics of the K(+)-form of five G-quadruplexes with sequences: d(G(2)T(2)G(2)TGTG(2)T(2)G(2)), G2, d(G(3)T(2)G(3)TGTG(3)T(2)G(3)), G3, their analogs where T is replaced with U, G2-U and G3-U, and r(G(2)U(2)G(2)UGUG(2)U(2)G(2)), rG2. These G-quadruplexes show CD spectra characteristic of the "chair" conformation (G2 and G2-U), or "basket" conformation (rG2); or a mixture of these two conformers (G3 and G3-U). Thermodynamic profiles show that the favorable folding of each G-quadruplex results from the typical compensation of a favorable enthalpy and unfavorable entropy contributions. G-quadruplex stability increase in the following order (in ΔG°(20)): rG2 (-1.3 kcal/mol) < G2 < G2-U <G3-U (chair) < G3 (chair) <G3-U (basket) < G3 (basket) (-8.6 kcal/mol), due to favorable enthalpy contribution from the stacking of G-quartets. We used ITC to determine thermodynamic binding profiles for the interaction of the minor groove ligands, netropsin and distamycin, with each G-quadruplex. Both ligands bind with high exothermic enthalpies (~ -10.8 kcal/mol), 1:1 stoichiometries, and weak affinities (~8 × 10(4) M(-1)). The similarity of the binding thermodynamic profiles, together with the absence of induced Cotton effects, indicates a surface or outside binding mode. We speculate that the top and bottom surfaces of the G-quadruplex comprise the potential MGBL binding sites, where the ligand lies on the surface forming van der Waals interactions with the guanines of the G-quartets and loop nucleotides.

Revisiting the association of cationic groove-binding drugs to DNA using a Poisson-Boltzmann approach

Proper modeling of nonspecific salt-mediated electrostatic interactions is essential to understanding the binding of charged ligands to nucleic acids. Because the linear Poisson-Boltzmann equation (PBE) and the more approximate generalized Born approach are applied routinely to nucleic acids and their interactions with charged ligands, the reliability of these methods is examined vis-à-vis an efficient nonlinear PBE method. For moderate salt concentrations, the negative derivative, SK(pred), of the electrostatic binding free energy, DeltaG(el), with respect to the logarithm of the 1:1 salt concentration, [M(+)], for 33 cationic minor groove drugs binding to AT-rich DNA sequences is shown to be consistently negative and virtually constant over the salt range considered (0.1-0.4 M NaCl). The magnitude of SK(pred) is approximately equal to the charge on the drug, as predicted by counterion condensation theory (CCT) and observed in thermodynamic binding studies. The linear PBE is shown to overestimate the magnitude of SK(pred), whereas the nonlinear PBE closely matches the experimental results. The PBE predictions of SK(pred) were not correlated with DeltaG(el) in the presence of a dielectric discontinuity, as would be expected from the CCT. Because this correlation does not hold, parameterizing the PBE predictions of DeltaG(el) against the reported experimental data is not possible. Moreover, the common practice of extracting the electrostatic and nonelectrostatic contributions to the binding of charged ligands to biopolyelectrolytes based on the simple relation between experimental SK values and the electrostatic binding free energy that is based on CCT is called into question by the results presented here. Although the rigid-docking nonlinear PB calculations provide reliable predictions of SK(pred), at least for the charged ligand-nucleic acid complexes studied here, accurate estimates of DeltaG(el) will require further development in theoretical and experimental approaches.

DNA minor groove induced dimerization of heterocyclic cations: compound structure, binding affinity, and specificity for a TTAA site

With the increasing number and variations of genome sequences available, control of gene expression with synthetic, cell-permeable molecules is within reach. The variety of sequence-specific binding agents is, however, still quite limited. Many minor groove binding agents selectivity recognize AT over GC sequences but have less ability to distinguish among different AT sequences. The goal with this article is to develop compounds that can bind selectively to different AT sequences. A number of studies indicate that AATT and TTAA sequences have significantly different physical and interaction properties and different requirements for minor groove recognition. Although it has been difficult to get minor groove binding at TTAA, DB293, a phenyl-furan-benzimidazole diamidine, was found to bind as a strong, cooperative dimer at TTAA but with no selectivity over AATT. In order to improve selectivity, we made modifications to each unit of DB293. Binding affinities and stoichiometries obtained from biosensor-surface plasmon resonance experiments show that DB1003, a furan-furan-benzimidazole diamidine, binds strongly to TTAA as a dimer and has selectivity (K(TTAA)/K(AATT)=6). CD and DNase I footprinting studies confirmed the preference of this compound for TTAA. In summary, (i) a favorable stacking surface provided by the pi system, (ii) H-bond donors to interact with TA base pairs at the floor of the groove provided by a benzimidazole (or indole) -NH and amidines, and (iii) appropriate curvature of the dimer complex to match the curvature of the minor groove play important roles in differentiating the TTAA and AATT minor grooves.

Calorimetric and spectroscopic studies of aminoglycoside binding to AT-rich DNA triple helices

Calorimetric and fluorescence techniques were used to characterize the binding of aminoglycosides-neomycin, paromomycin, and ribostamycin, with 5'-dA(12)-x-dT(12)-x-dT(12)-3' intramolecular DNA triplex (x = hexaethylene glycol) and poly(dA).2poly(dT) triplex. Our results demonstrate the following features: (1) UV thermal analysis reveals that the T(m) for triplex decreases with increasing pH value in the presence of neomycin, while the T(m) for the duplex remains unchanged. (2) The binding affinity of neomycin decreases with increased pH, although there is an increase in observed binding enthalpy. (3) ITC studies conducted in two buffers (sodium cacodylate and MOPS) yield the number of protonated drug amino groups (Deltan) as 0.29 and 0.40 for neomycin and paromomycin interaction with 5'-dA(12)-x-dT(12)-x-dT(12)-3', respectively. (4) The specific heat capacity change (DeltaC(p)) determined by ITC studies is negative, with more negative values at lower salt concentrations. From 100 mM to 250 mM KCl, the DeltaC(p) ranges from -402 to -60 cal/(mol K) for neomycin. At pH 5.5, a more positive DeltaC(p) is observed, with a value of -98 cal/(mol K) at 100 mM KCl. DeltaC(p) is not significantly affected by ionic strength. (5) Salt dependence studies reveal that there are at least three amino groups of neomycin participating in the electrostatic interactions with the triplex. (6) FID studies using thiazole orange were used to derive the AC(50) (aminoglycoside concentration needed to displace 50% of the dye from the triplex) values. Neomycin shows a seven fold higher affinity than paromomycin and eleven fold higher affinity than ribostamycin at pH 6.8. (7) Modeling studies, consistent with UV and ITC results, show the importance of an additional positive charge in triplex recognition by neomycin. The modeling and thermodynamic studies indicate that neomycin binding to the DNA triplex depends upon significant contributions from charge as well as shape complementarity of the drug to the DNA triplex Watson-Hoogsteen groove.

DNA intercalation of methylene blue and quinacrine: new insights into base and sequence specificity from structural and thermodynamic studies with polynucleotides

The binding of the known DNA intercalators methylene blue and quinacrine with four sequence specific polynucleotides, viz. poly(dG-dC).poly(dG-dC), poly(dG).poly(dC), poly(dA-dT).poly(dA-dT) and poly(dA).poly(dT), have been compared using absorbance, fluorescence, competition dialysis and thermal melting and the thermodynamic aspects of the interaction studied. In all the cases, non-cooperative binding phenomena obeying neighbor exclusion principle was observed though the affinity was remarkably higher for quinacrine and the nature of the binding was characterized to be true intercalation. The data on the salt dependence of binding derived from the plot of log Kvs. log[Na(+)] revealed a slope of around 1.0, consistent with the values predicted by the theories for the binding of monovalent cations, and contained contributions from polyelectrolytic and non-polyelectrolytic forces. The bindings were characterized by strong stabilization of the polynucleotides against thermal strand separation in both optical melting as well as differential scanning calorimetry studies. The data analyzed from the thermal melting and isothermal titration calorimetry studies were in close proximity to those obtained from absorption spectral titration data. Isothermal titration calorimetry results revealed the bindings to poly(dG-dC).poly(dG-dC), poly(dG).poly(dC) and poly(dA-dT).poly(dA-dT) to be exothermic and favoured by both negative enthalpy and large favourable positive entropy changes, while that to poly(dA).poly(dT) was endothermic and entropy driven. The heat capacity changes obtained from temperature dependence of enthalpy gave negative values to all polynucleotides. New insights on the molecular aspects of interaction of these molecules to DNA have emerged from these studies.

Large, sequence-dependent effects on DNA conformation by minor groove binding compounds

To determine what topological changes antiparasitic heterocyclic dications can have on kinetoplast DNA, we have constructed ligation ladders, with phased A5 and ATATA sequences in the same flanking sequence context, as models. Bending by the A5 tract is observed, as expected, while the ATATA sequence bends DNA very little. Complexes of these DNAs with three diamidines containing either furan, thiophene or selenophene groups flanked by phenylamidines were investigated along with netropsin. With the bent A5 ladder the compounds caused either a slight increase or decrease in the bending angle. Surprisingly, however, with ATATA all of the compounds caused significant bending, to values close to or even greater than the A5 bend angle. Results with a mixed cis sequence, which has one A5 and one ATATA, show that the compounds bend ATATA in the same direction as a reference A5 tract, that is, into the minor groove. These results are interpreted in terms of a groove structure for A5 which is largely pre-organized for a fit to the heterocyclic amidines. With ATATA the groove is intrinsically wider and must close to bind the compounds tightly. The conformational change at the binding site then leads to significant bending of the alternating DNA sequence.

The diamidine DB75 targets the nucleus of Plasmodium falciparum

Background

DB289, [2,5-bis(4-amidinophenyl)furan bis-O-methylamidoxime], is a broad spectrum anti-parasitic compound which has been shown to be effective against malaria in recent clinical trials. DB75, [2,5-bis(4-amidinophenyl)furan], is the active metabolite of this drug. The objective of this study was to determine the mechanism of action of DB75 in Plasmodium falciparum.

Methods

Live parasites were observed by confocal microscopy after treatment with organelle specific dyes and DB75, an inherently fluorescent compound. Parasites were exposed to DB75 and assessed for growth and morphological changes over time using blood smears and light microscopy. Also, to determine if DB75 affects gene transcription, real time PCR was used to monitor transcript levels over time for six developmentally expressed genes, including trophozoite antigen R45-like (PFD1175w), lactate dehydrogenase (PF13_0141), DNA primase (PFI0530c), isocitrate dehydrogenase (PF13_0242), merozoite surface protein-1 (PFI1475w), and merozoite surface protein-7 (PF13_0197).

Results

The results show that DB75 localizes in the parasite nucleus but not in other organelles. Once rings are exposed, parasites mature to the trophozoite stage and stall. No stage-dependent or gene-specific inhibition of transcription was seen. However, DB75 delayed peak transcription of trophozoite-stage genes.

Conclusion

Taken together, DB75 appears to concentrate in the nucleus and delay parasite maturation.

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.

Effect of 3,4-ethylenedioxy-extension of thiophene core on the DNA/RNA binding properties and biological activity of bisbenzimidazole amidines

Novel bisbenzimidazoles (4-6), characterized by 3,4-ethylenedioxy-extension of thiophene core, revealed pronounced affinity and strong thermal stabilization effect toward ds-DNA. They interact within ds-DNA grooves as dimmers or even oligomers and agglomerate along ds-RNA. Compounds 4-6 have shown moderate to strong antiproliferative effect toward panel of eight carcinoma cell lines. Compound 5 displayed the best inhibitory potential and in equitoxic concentration (IC(50)=1 x 10(-6)M) induced accumulation of cells in G2/M phase after 48 h of incubation. Fluorescence microscopy showed that 5 entered into live HeLa cells within 30 min, but did not accumulate in nuclei even after 2.5h. Compound 5 inhibited the growth of Trypanosome cruzi epimastigotes (IC(50)=4.3 x 10(-6)M).

A role for water molecules in DNA-ligand minor groove recognition

Targeting the minor groove of DNA through binding to a small molecule has long been considered an important molecular-recognition strategy in biology. A wide range of synthetic heterocyclic molecules bind noncovalently in the minor groove of the double helix and are also effective against a number of human and animal diseases. A classic structural concept, the isohelicity principle, has guided much of this work: such heterocyclic molecules require a shape that complements the convex surface of the minor groove. Researchers have used this principle to design molecules that can read DNA sequences. This principle also predicts that molecules that lack the complementary shape requirement would only bind weakly to DNA. Recently, however, researchers have unexpectedly found that some essentially linear compounds, which do not have this feature, can have high DNA affinity. In this Account, we discuss an alternative recognition concept based on these new findings. We demonstrate that highly structured water molecules can play a key role in mediating between the ligand and DNA minor groove without loss of binding affinity. Combined structural and thermodynamic approaches to understanding the behavior of these molecules have shown that there are different categories of bound water in their DNA complexes. For example, application of this water-bridging concept to the phenylamidine platform has resulted in the discovery of molecules with high levels of biological activity and low nonspecific toxicity. Some of these molecules are now in advanced clinical trials.

Influence of substituent modifications on the binding of 2-amino-1,8-naphthyridines to cytosine opposite an AP site in DNA duplexes: thermodynamic characterization

Here, we report on a significant effect of substitutions on the binding affinity of a series of 2-amino-1,8-naphthyridines, i.e., 2-amino-1,8-naphthyridine (AND), 2-amino-7-methyl-1,8-naphthyridine (AMND), 2-amino-5,7-dimethyl-1,8-naphthyridine (ADMND) and 2-amino-5,6,7-trimethyl-1,8-naphthyridine (ATMND), all of which can bind to cytosine opposite an AP site in DNA duplexes. Fluorescence titration experiments show that the binding affinity for cytosine is effectively enhanced by the introduction of methyl groups to the naphthyridine ring, and the 1:1 binding constant (10⁶ M⁻¹) follows in the order of AND (0.30) < AMND (2.7) < ADMND (6.1) < ATMND (19) in solutions containing 110 mM Na⁺ (pH 7.0, at 20°C). The thermodynamic parameters obtained by isothermal titration calorimetry experiments indicate that the introduction of methyl groups effectively reduces the loss of binding entropy, which is indeed responsible for the increase in the binding affinity. The heat capacity change (ΔC_p), as determined from temperature dependence of the binding enthalpy, is found to be significantly different between AND (−161 cal/mol K) and ATMND (−217 cal/mol K). The hydrophobic contribution appears to be a key force to explain the observed effect of substitutions on the binding affinity when the observed binding free energy (ΔG_obs) is dissected into its component terms.

Novel linear triaryl guanidines, N-substituted guanidines and potential prodrugs as antiprotozoal agents

A series of triaryl guanidines and N-substituted guanidines designed to target the minor groove of DNA were synthesized and evaluated as antiprotozoal agents. Selected carbamate prodrugs of these guanidines were assayed for their oral efficacy. The linear triaryl bis-guanidines 6a,b were prepared from their corresponding diamines 4a,b through the intermediate BOC protected bis-guanidines 5a,b followed by acid catalyzed deprotection. The N-substituted guanidino analogues 9c-f were obtained in three steps starting by reacting the diamines 4a,b with ethyl isothiocyanatoformate to give the carbamoyl thioureas 7a,b. Subsequent condensation of 7a,b with various amines in the presence of EDCI provided the carbamoyl N-substituted guanidine intermediates 8a-f which can also be regarded as potential prodrugs for the guanidino derivatives. Compounds 9c-f were obtained via the base catalyzed decarbamoylation of 8a-f. The DNA binding affinities for the target dicationic bis-guanidines were assessed by DeltaT(m) values. In vitro antiprotozoal screening of the new compounds showed that derivatives 6a, 9c and 9e possess high to moderate activity against Trypanosoma brucei rhodesiense (T.b.r.) and Plasmodium falciparum (P.f.). While the prodrugs did not yield cures upon oral administration in the antitrypanosomal STIB900 mouse model, compounds 8a and 8c prolonged the survival of the treated mice.

Heterocyclic diamidine interactions at AT base pairs in the DNA minor groove: effects of heterocycle differences, DNA AT sequence and length

Given the increasing significance of diamidines as DNA-targeted therapeutics and biotechnology reagents, it is important to establish the variations in thermodynamic quantities that characterize the interactions of closely related compounds to different sequence AT binding sites. In this study, an array of methods including biosensor-surface plasmon resonance (SPR), isothermal titration microcalorimetry (ITC), circular dichroism (CD), thermal melting (Tm) and molecular modeling have been used to characterize the binding of dicationic diamidines related to DB75 (amidine-phenyl-furan-phenyl-amidine) with alternating and nonalternating AT sequences. Conversion of the central furan of DB75 to other similar groups, such as thiophene or selenophene, can yield compounds with increased affinity and sequence binding selectivity for the minor groove. Calorimetric measurements revealed that the thermodynamic parameters (Delta G, Delta H, Delta S) that drive diamidine binding to alternating and nonalternating oligomers can be quite different and depend on both DNA sequence and length. Small changes in a compound can have major effects on DNA interactions. By choosing an appropriate central group it is possible to "tune" the shape of the molecule to match DNA for enhanced affinity and sequence recognition.