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

Systematic Evaluation of Benchmark G4 Probes and G4 Clinical Drugs using three Biophysical Methods: A Guideline to Evaluate Rapidly G4-Binding Affinity

G-quadruplex DNA structures (G4) are proven to interfere with most genetic and epigenetic processes. Small molecules binding these structures (G4 ligands) are invaluable tools to probe G4-biology and address G4-druggability in various diseases (cancer, viral infections). However, the large number of reported G4 ligands (>1000) could lead to confusion while selecting one for a given application. Herein we conducted a systematic affinity ranking of 11 popular G4 ligands vs 5 classical G4 sequences using FRET-melting, G4-FID assays and SPR. Interestingly SPR data globally align with the rankings obtained from the two semi-quantitative assays despite discrepancies due to limits and characteristics of each assay. In the whole, PhenDC3 emerges as the most potent binder irrespective of the G4 sequence. Immediately below PDS, PDC-360A, BRACO19, TMPyP4 and RHPS4 feature strong to medium binding again with poor G4 topology discrimination. More strikingly, the G4 drugs Quarfloxin, CX5461 and c-PDS exhibit weak affinity with all G4s studied. Finally, NMM and Cu-ttpy showed heterogeneous behaviors due, in part, to their physicochemical particularities poorly compatible with screening conditions. The remarkable properties of PhenDC3 led us to propose its use for benchmarking FRET-melting and G4-FID assays for rapid G4-affinity evaluation of newly developed ligands.

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.

Contrasting Effect of Salts on the Binding of Antimalarial Drug Hydroxychloroquine with Different Sequences of Duplex DNA

Hydroxychloroquine (HCQ) is an important antimalarial drug which functions plausibly by targeting the DNA of parasites. Salts play a crucial role in the functionality of various biological processes. Hence, the effect of salts (NaCl and MgCl₂) on the binding of HCQ with AT- and CG-DNAs as well as the binding-induced stability of both sequences of DNAs have been investigated using the spectroscopic and molecular dynamics (MD) simulation methods. It has been found that the effect of salts on the binding of HCQ is highly sensitive to the nature of ions as well as DNA sequences. The effect of ions is opposite for the binding of AT- and CG-DNAs as the presence of Mg²⁺ ions enhances the binding of HCQ with AT-DNA, whereas the binding of HCQ with CG-DNA gets decreased on the addition of both ions. Similarly, the presence of Mg²⁺ enhances the stabilization of HCQ-bound AT-DNA, whereas the effect is opposite for the CG-DNA in the presence of both the ions. The MD simulation study suggests that the hydration states of both ions are different and they interact differently in the minor and major grooves of both the sequences of DNA which may be one of the reasons for the different binding of HCQ with these two sequences of DNA in the presence of salts. The information about the effect of salts on the binding of HCQ with DNAs in a sequence-specific manner may be useful in understanding the mechanism of the action and toxicity effect of HCQ against malaria.

Thermodynamic Factors That Drive Sequence-Specific DNA Binding of Designed, Synthetic Minor Groove Binding Agents

Ken Breslauer began studies on the thermodynamics of small cationic molecules binding in the DNA minor groove over 30 years ago, and the studies reported here are an extension of those ground-breaking reports. The goals of this report are to develop a detailed understanding of the binding thermodynamics of pyridine-based sequence-specific minor groove binders that have different terminal cationic groups. We apply biosensor-surface plasmon resonance and ITC methods to extend the understanding of minor groove binders in two directions: (i) by using designed, heterocyclic dicationic minor groove binders that can incorporate a G•C base pair (bp), with flanking AT base pairs, into their DNA recognition site, and bind to DNA sequences specifically; and (ii) by using a range of flanking AT sequences to better define molecular recognition of the minor groove. A G•C bp in the DNA recognition site causes a generally more negative binding enthalpy than with most previously used pure AT binding sites. The binding is enthalpy-driven at 25 °C and above. The flanking AT sequences also have a large effect on the binding energetics with the -AAAGTTT- site having the strongest affinity. As a result of these studies, we now have a much better understanding of the effects of the DNA sequence and compound structure on the molecular recognition and thermodynamics of minor groove complexes.

A fluorescent reporter on electrostatic DNA-ligand interactions

Among the various types of interactions between biomolecules, electrostatic interactions dominate as these are long-range interactions and are often a generic first step in the recruitment of specific ligands. DNA, being a highly charged molecule, attracts a plethora of molecules. Interactions between DNA and proteins or small molecules shape the overall function of the cell. Various processes such as DNA replication, DNA repair, synthesis of mRNA, and packaging of DNA are mediated by interactions between protein molecules and DNA that are predominantly electrostatic. Here, we present a fluorescence resonance energy transfer (FRET)-based probe which can report on the electrostatic interactions between the negatively-charged DNA and positively-charged metal ions, oligopeptides, as well as DNA groove-binding drug molecules. The simplicity, sensitivity, and versatility of the DNA-based probe makes it suited for applications where specific protein-DNA interactions can be probed, and DNA-binding drugs can be discovered in high-throughput screens of compound libraries. This is particularly relevant given that some of the most potent antitumor and antimicrobial drugs associate with DNA electrostatically.

Engineered modular heterocyclic-diamidines for sequence-specific recognition of mixed AT/GC base pairs at the DNA minor groove

This report describes a breakthrough in a project to design minor groove binders to recognize any sequence of DNA. A key goal is to invent synthetic chemistry for compound preparation to recognize an adjacent GG sequence that has been difficult to target. After trying several unsuccessful compound designs, an N-alkyl-benzodiimidazole structure was selected to provide two H-bond acceptors for the adjacent GG-NH groups. Flanking thiophenes provide a preorganized structure with strong affinity, DB2831, and the structure is terminated by phenyl-amidines. The binding experimental results for DB2831 with a target AAAGGTTT sequence were successful and include a high ΔT m, biosensor SPR with a K D of 4 nM, a similar K D from fluorescence titrations and supporting competition mass spectrometry. MD analysis of DB2831 bound to an AAAGGTTT site reveals that the two unprotonated N of the benzodiimidazole group form strong H-bonds (based on distance) with the two central G-NH while the central -CH of the benzodiimidazole is close to the -C[double bond, length as m-dash]O of a C base. These three interactions account for the strong preference of DB2831 for a -GG- sequence. Surprisingly, a complex with one dynamic, interfacial water is favored with 75% occupancy.

Impact of Self-Association on the Architectural Properties of Bacterial Nucleoid Proteins

The chromosomal DNA of bacteria is folded into a compact body called the nucleoid, which is composed essentially of DNA (∼80%), RNA (∼10%), and a number of different proteins (∼10%). These nucleoid proteins act as regulators of gene expression and influence the organization of the nucleoid by bridging, bending, or wrapping the DNA. These so-called architectural properties of nucleoid proteins are still poorly understood. For example, the reason why certain proteins compact the DNA coil in certain environments but make the DNA more rigid instead in other environments is the subject of ongoing debates. Here, we address the question of the impact of the self-association of nucleoid proteins on their architectural properties and try to determine whether differences in self-association are sufficient to induce large changes in the organization of the DNA coil. More specifically, we developed two coarse-grained models of proteins, which interact identically with the DNA but self-associate differently by forming either clusters or filaments in the absence of the DNA. We showed through Brownian dynamics simulations that self-association of the proteins dramatically increases their ability to shape the DNA coil. Moreover, we observed that cluster-forming proteins significantly compact the DNA coil (similar to the DNA-bridging mode of H-NS proteins), whereas filament-forming proteins significantly increase the stiffness of the DNA chain instead (similar to the DNA-stiffening mode of H-NS proteins). This work consequently suggests that the knowledge of the DNA-binding properties of the proteins is in itself not sufficient to understand their architectural properties. Rather, their self-association properties must also be investigated in detail because they might actually drive the formation of different DNA-protein complexes.

Dissecting Dynamic and Hydration Contributions to Sequence-Dependent DNA Minor Groove Recognition

Sequence selectivity is a critical attribute of DNA-binding ligands and underlines the need for detailed molecular descriptions of binding in representative sequence contexts. We investigated the binding and volumetric properties of DB1976, a model bis(benzimidazole)-selenophene diamidine compound with emerging therapeutic potential in acute myeloid leukemia, debilitating fibroses, and obesity-related liver dysfunction. To sample the scope of cognate DB1976 target sites, we evaluated three dodecameric duplexes spanning >10³-fold in binding affinity. The attendant changes in partial molar volumes varied substantially, but not in step with binding affinity, suggesting distinct modes of interactions in these complexes. Specifically, whereas optimal binding was associated with loss of hydration water, low-affinity binding released more hydration water. Explicit-atom molecular dynamics simulations showed that minor groove binding perturbed the conformational dynamics and hydration at the termini and interior of the DNA in a sequence-dependent manner. The impact of these distinct local dynamics on hydration was experimentally validated by domain-specific interrogation of hydration with salt, which probed the charged axial surfaces of oligomeric DNA preferentially over the uncharged termini. Minor groove recognition by DB1976, therefore, generates dynamically distinct domains that can make favorable contributions to hydration release in both high- and low-affinity binding. Because ligand binding at internal sites of DNA oligomers modulates dynamics at the termini, the results suggest both short- and long-range dynamic effects along the DNA target that can influence their effectiveness as low-MW competitors of protein binding.

Thermodynamics and site stoichiometry of DNA binding by a large antiviral hairpin polyamide

PA1 (dIm-PyPyβPyPyPy-γ-PyPyβPyPyPyPyβ-Ta) is a large (14-ring) hairpin polyamide that was designed to recognize the DNA sequence 5'-W₂GW₇-3', where W is either A or T. As is common among the smaller 6-8-ring hairpin polyamides (PAs), it binds its target recognition sequence with low nM affinity. However, in addition to its large size, it is distinct from these more extensively characterized PAs in its high tolerance for mismatches and antiviral properties. In ongoing attempts to understand the basis for these distinctions, we conducted thermodynamics studies of PA1-DNA interactions. The temperature dependence of binding affinity was measured using TAMRA-labeled hairpin DNAs containing a single target sequence. PA1 binding to either an ATAT/TATA or an AAAA/TTTT pattern is consistently entropically driven. This is in contrast to the A/T pattern-dependent driving forces for DNA binding by netropsin, distamycin, and smaller hairpin polyamides. Analysis of the salt dependence of PA1-DNA binding reveals that within experimental error, there is no dependence on ionic strength, indicating that the polyelectrolyte effect does not contribute to PA1-DNA binding energetics. This is similar to that observed for smaller PAs. PA1-DNA recognition sequence binding stoichiometries were determined at both nM (fluorescence) and μM (circular dichroism) concentrations. With all sequences and under both conditions, multiple PA1 molecules bind the small DNA hairpin that contains only a single recognition sequence. Implications for these observations are discussed.

Role of Salt Valency in the Switch of H-NS Proteins between DNA-Bridging and DNA-Stiffening Modes

This work investigates the interactions of H-NS proteins and bacterial genomic DNA through computer simulations performed with a coarse-grained model. The model was developed specifically to study the switch of H-NS proteins from the DNA-stiffening to the DNA-bridging mode, which has been observed repeatedly upon addition of multivalent cations to the buffer but is still not understood. Unraveling the corresponding mechanism is all the more crucial, as the regulation properties of H-NS proteins, as well as other nucleoid proteins, are linked to their DNA-binding properties. The simulations reported here support a mechanism, according to which the primary role of multivalent cations consists in decreasing the strength of H-NS/DNA interactions compared to H-NS/H-NS interactions, with the latter ones becoming energetically favored with respect to the former ones above a certain threshold of the effective valency of the cations of the buffer. Below the threshold, H-NS dimers form filaments, which stretch along the DNA molecule but are quite inefficient in bridging genomically distant DNA sites (DNA-stiffening mode). In contrast, just above the threshold, H-NS dimers form three-dimensional clusters, which are able to connect DNA sites that are distant from the genomic point of view (DNA-bridging mode). The model provides clear rationales for the experimental observations that the switch between the two modes is a threshold effect and that the ability of H-NS dimers to form higher order oligomers is crucial for their bridging capabilities.

Impact of Linker Length and Composition on Fragment Binding and Cell Permeation: Story of a Bisbenzimidazole Dye Fragment

Small molecules that modulate biological functions are targets of modern day drug discovery efforts. In a common platform fragment-based drug discovery, two fragments that bind to adjacent sites on a target are identified and are then linked together using different linkers to identify the linkage for optimum activity. What are not known from these studies are the effects these linkers, which typically contain C, H, and O atoms, have on the properties of the individual fragment. Herein, we investigate such effects in a bisbenzimidazole fragment whose derivatives have a wide range of therapeutic applications in nucleic acid recognition, sensing, and photodynamic therapy and as cellular probes. We report a dramatic effect of linker length and composition of alkynyl (clickable) Hoechst 33258 derivatives in target binding and cell uptake. We show that the binding of Hoechst 33258-modeled bisbenzimidazoles (1-9) that contain linkers of varying lengths (3-21 atoms) display length- and composition-dependent variation in B-DNA stabilization using a variety of spectroscopic methods. For a dodecamer DNA duplex, the thermal stabilization varied from 0.3 to 9.0 °C as the linker length increased from 3 to 21 atoms, respectively. Compounds with linker lengths of ≤11 atoms (such as compounds 1 and 5) are localized in the nucleus, while compounds with long linkers (such as compounds 8 and 9) are distributed in the extranuclear space, as well, with possible interactions with extranuclear targets. These findings provide insights into future drug design by revealing how linkers can influence the biophysical and cellular properties of individual drug fragments.

Determination of Ion Atmosphere Effects on the Nucleic Acid Electrostatic Potential and Ligand Association Using AH+·C Wobble Formation in Double-Stranded DNA

The high charge density of nucleic acids and resulting ion atmosphere profoundly influence the conformational landscape of RNA and DNA and their association with small molecules and proteins. Electrostatic theories have been applied to quantitatively model the electrostatic potential surrounding nucleic acids and the effects of the surrounding ion atmosphere, but experimental measures of the potential and tests of these models have often been complicated by conformational changes and multisite binding equilibria, among other factors. We sought a simple system to further test the basic predictions from electrostatics theory and to measure the energetic consequences of the nucleic acid electrostatic field. We turned to a DNA system developed by Bevilacqua and co-workers that involves a proton as a ligand whose binding is accompanied by formation of an internal AH⁺·C wobble pair [Siegfried, N. A., et al. Biochemistry, 2010, 49, 3225]. Consistent with predictions from polyelectrolyte models, we observed logarithmic dependences of proton affinity versus salt concentration of −0.96 ± 0.03 and −0.52 ± 0.01 with monovalent and divalent cations, respectively, and these results help clarify prior results that appeared to conflict with these fundamental models. Strikingly, quantitation of the ion atmosphere content indicates that divalent cations are preferentially lost over monovalent cations upon A·C protonation, providing experimental indication of the preferential localization of more highly charged cations to the inner shell of the ion atmosphere. The internal AH⁺·C wobble system further allowed us to parse energetic contributions and extract estimates for the electrostatic potential at the position of protonation. The results give a potential near the DNA surface at 20 mM Mg²⁺ that is much less substantial than at 20 mM K⁺ (−120 mV vs −210 mV). These values and difference are similar to predictions from theory, and the potential is substantially reduced at higher salt, also as predicted; however, even at 1 M K⁺ the potential remains substantial, counter to common assumptions. The A·C protonation module allows extraction of new properties of the ion atmosphere and provides an electrostatic meter that will allow local electrostatic potential and energetics to be measured within nucleic acids and their complexes with proteins.

Investigation of the electrostatic and hydration properties of DNA minor groove-binding by a heterocyclic diamidine by osmotic pressure

Previous investigations of sequence-specific DNA binding by model minor groove-binding compounds showed that the ligand/DNA complex was destabilized in the presence of compatible co-solutes. Inhibition was interpreted in terms of osmotic stress theory as the uptake of significant numbers of excess water molecules from bulk solvent upon complex formation. Here, we interrogated the AT-specific DNA complex formed with the symmetric heterocyclic diamidine DB1976 as a model for minor groove DNA recognition using both ionic (NaCl) and non-ionic cosolutes (ethylene glycol, glycine betaine, maltose, nicotinamide, urea). While the non-ionic cosolutes all destabilized the ligand/DNA complex, their quantitative effects were heterogeneous in a cosolute- and salt-dependent manner. Perturbation with NaCl in the absence of non-ionic cosolute showed that preferential hydration water was released upon formation of the DB1976/DNA complex. As salt probes counter-ion release from charged groups such as the DNA backbone, we propose that the preferential hydration uptake in DB1976/DNA binding observed in the presence of osmolytes reflects the exchange of preferentially bound cosolute with hydration water in the environs of the bound DNA, rather than a net uptake of hydration waters by the complex.

Modulation of DNA-polyamide interaction by β-alanine substitutions: a study of positional effects on binding affinity, kinetics and thermodynamics

Hairpin polyamides (PAs) are an important class of sequence-specific DNA minor groove binders, and frequently employ a flexible motif, β-alanine (β), to reduce the molecular rigidity to maintain the DNA recognition register. To better understand the diverse effects that β can have on DNA-PA binding affinity, selectivity, and especially kinetics, which have rarely been reported, we have initiated a detailed study for an eight-heterocyclic hairpin PA and its β derivatives with their cognate and mutant sequences. With these derivatives, all internal pyrroles of the parent PA are systematically substituted with single or double βs. A set of complementary experiments have been conducted to evaluate the molecular interactions in detail: UV-melting, biosensor-surface plasmon resonance, circular dichroism and isothermal titration calorimetry. The β substitutions generally weaken the binding affinities of these PAs with cognate DNA, and have large and diverse influences on PA binding kinetics in a position- and number-dependent manner. The DNA base mutations have also shown positional effects on the binding of a single PA. Besides the β substitutions, the monocationic Dp group [3-(dimethylamino)propylamine] in parent PA has been modified into a dicationic Ta group (3,3'-diamino-N-methyldipropylamine) to minimize the frequently observed PA aggregation with ITC experiments. The results clearly show that the Ta modification not only maintains the DNA binding mode and affinity of PA, but also significantly reduces PA aggregation and allows the complete thermodynamic signature of eight-ring hairpin PA to be determined for the first time. This combined set of results significantly extends our understanding of the energetic basis of specific DNA recognition by PAs.

Understanding mixed sequence DNA recognition by novel designed compounds: the kinetic and thermodynamic behavior of azabenzimidazole diamidines

Sequence-specific recognition of DNA by small organic molecules offers a potentially effective approach for the external regulation of gene expression and is an important goal in cell biochemistry. Rational design of compounds from established modules can potentially yield compounds that bind strongly and selectively with specific DNA sequences. An initial approach is to start with common A·T bp recognition molecules and build in G·C recognition units. Here we report on the DNA interaction of a synthetic compound that specifically binds to a G·C bp in the minor groove of DNA by using an azabenzimidazole moiety. The detailed interactions were evaluated with biosensor-surface plasmon resonance (SPR), isothermal calorimetric (ITC), and mass spectrometry (ESI-MS) methods. The compound, DB2277, binds with single G·C bp containing sequences with sub-nanomolar potency and displays slow dissociation kinetics and high selectivity. A detailed thermodynamic and kinetic study at different experimental salt concentrations and temperatures shows that the binding free energy is salt concentration dependent but essentially temperature independent under our experimental conditions, and binding enthalpy is temperature dependent but salt concentration independent. The results show that in the proper compound structural context novel heterocyclic cations can be designed to strongly recognize complex DNA sequences.

The unusual monomer recognition of guanine-containing mixed sequence DNA by a dithiophene heterocyclic diamidine

DB1255 is a symmetrical diamidinophenyl-dithiophene that exhibits cellular activity by binding to DNA and inhibiting binding of ERG, an ETS family transcription factor that is commonly overexpressed or translocated in leukemia and prostate cancer [Nhili, R., Peixoto, P., Depauw, S., Flajollet, S., Dezitter, X., Munde, M. M., Ismail, M. A., Kumar, A., Farahat, A. A., Stephens, C. E., Duterque-Coquillaud, M., Wilson, W. D., Boykin, D. W., and David-Cordonnier, M. H. (2013) Nucleic Acids Res. 41, 125–138]. Because transcription factor inhibition is complex but is an attractive area for anticancer and antiparasitic drug development, we have evaluated the DNA interactions of additional derivatives of DB1255 to gain an improved understanding of the biophysical chemistry of complex function and inhibition. DNase I footprinting, biosensor surface plasmon resonance, and circular dichroism experiments show that DB1255 has an unusual and strong monomer binding mode in minor groove sites that contain a single GC base pair flanked by AT base pairs, for example, 5′-ATGAT-3′. Closely related derivatives, such as compounds with the thiophene replaced with furan or selenophane, bind very weakly to GC-containing sequences and do not have biological activity. DB1255 is selective for the ATGAT site; however, a similar sequence, 5′-ATGAC-3′, binds DB1255 more weakly and does not produce a footprint. Molecular docking studies show that the two thiophene sulfur atoms form strong, bifurcated hydrogen bond-type interactions with the G-N-H sequence that extends into the minor groove while the amidines form hydrogen bonds to the flanking AT base pairs. The central dithiophene unit of DB1255 thus forms an excellent, but unexpected, single-GC base pair recognition module in a monomer minor groove complex.

DNA Binding Polyamides and the Importance of DNA Recognition in their use as Gene-Specific and Antiviral Agents

There is a long history for the bioorganic and biomedical use of N-methyl-pyrrole-derived polyamides (PAs) that are higher homologs of natural products such as distamycin A and netropsin. This work has been pursued by many groups, with the Dervan and Sugiyama groups responsible for many breakthroughs. We have studied PAs since about 1999, partly in industry and partly in academia. Early in this program, we reported methods to control cellular uptake of polyamides in cancer cell lines and other cells likely to have multidrug resistance efflux pumps induced. We went on to discover antiviral polyamides active against HPV31, where SAR showed that a minimum binding size of about 10 bp of DNA was necessary for activity. Subsequently we discovered polyamides active against two additional high-risk HPVs, HPV16 and 18, a subset of which showed broad spectrum activity against HPV16, 18 and 31. Aspects of our results presented here are incompatible with reported DNA recognition rules. For example, molecules with the same cognate DNA recognition properties varied from active to inactive against HPVs. We have since pursued the mechanism of action of antiviral polyamides, and polyamides in general, with collaborators at NanoVir, the University of Missouri-St. Louis, and Georgia State University. We describe dramatic consequences of β-alanine positioning even in relatively small, 8-ring polyamides; these results contrast sharply with prior reports. This paper was originally presented by JKB as a Keynote Lecture in the 2^nd International Conference on Medicinal Chemistry and Computer Aided Drug Design Conference in Las Vegas, NV, October 2013.