Jonathan B. Chaires

Criteria for the Mode of Binding of DNA Binding Agents

A complete characterization of DNA binding agents requires that their mode of binding to DNA be established. In the absence of high resolution structural data, the mode of binding is, of necessity, usually inferred indirectly from various solution studies. The purpose of this study is to show that only certain methods can be used reliably to infer the DNA binding mode. Comparative fluorescence and hydrodynamic studies using the proven intercalator ethidium and the groove binder Hoechst 33258 are described. The results of our studies show that while fluorescence intensity, polarization, and quenching measurements can detect a binding interaction of the ligand with DNA, none are sensitive indicators of the binding mode. Fluorescence contact energy transfer studies can reliably indicate intercalation, as can viscosity measurements. Our results illustrate reliable criteria that may be used to distinguish intercalation from groove binding in the absence of high resolution structural data.

ENERGETICS OF DRUG-DNA INTERACTIONS

Understanding the thermodynamics of drug binding to DNA is of both practical and fundamental interest. The practical interest lies in the contribution that thermodynamics can make to the rational design process for the development of new DNA targeted drugs. Thermodynamics offer key insights into the molecular forces that drive complex formation that cannot be obtained by structural or computational studies alone. The fundamental interest in these interactions lies in what they can reveal about the general problems of parsing and predicting ligand binding free energies. For these problems, drug-DNA interactions offer several distinct advantages, among them being that the structures of many drug-DNA complexes are known at high resolution and that such structures reveal that in many cases the drug acts as a rigid body, with little conformational change upon binding. Complete thermodynamic profiles (delta G, delta H, delta S, delta Cp) for numerous drug-DNA interactions have been obtained, with the help of high-sensitivity microcalorimetry. The purpose of this article is to offer a perspective on the interpretation of these thermodynamics parameters, and in particular how they might be correlated with known structural features. Obligatory conformational changes in the DNA to accommodate intercalators and the loss of translational and rotational freedom upon complex formation both present unfavorable free energy barriers for binding. Such barriers must be overcome by favorable free energy contributions from the hydrophobic transfer of ligand from solution into the binding site, polyelectrolyte contributions from coupled ion release, and molecular interactions (hydrogen and ionic bonds, van der Waals interactions) that form within the binding site. Theoretical and semiempirical tools that allow estimates of these contributions to be made will be discussed, and their use in dissecting experimental data illustrated. This process, even at the current level of approximation, can shed considerable light on the drug-DNA binding process.

Calorimetry and Thermodynamics in Drug Design

Modern instrumentation for calorimetry permits direct determination of enthalpy values for binding reactions and conformational transitions in biomolecules. Complete thermodynamic profiles consisting of free energy, enthalpy, and entropy may be obtained for reactions of interest in a relatively straightforward manner. Such profiles are of enormous value in drug design because they provide information about the balance of driving forces that cannot be obtained from structural or computational methods alone. This perspective shows several examples of the insight provided by thermodynamic data in drug design.

Drug—DNA interactions

Significant progress has been made over the past few years in studies of drug-DNA interactions. Structure-based design strategies have yielded new DNA-binding agents with clinical promise. The hairpin polyamides represent the result of a design strategy with outstanding potential. One specific molecule of this class has now been proven to inhibit the expression of a specific gene in vivo. A new bisintercalating anthracycline antibiotic binds with high affinity to DNA, and appears to overcome a specific form of multidrug resistance. Progress in fundamental studies of drug binding to DNA continues, with detailed thermodynamic studies providing new insights into the forces that drive complex formation. New tools have been developed in order to characterize both the binding mode and the sequence specificity of drug binding to DNA, tools that will enable the fundamental aspects of these biologically important reactions to be understood in more detail.

Thermal difference spectra: a specific signature for nucleic acid structures

We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (Tm). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.

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.

Kinetics and mechanism of K+- and Na+-induced folding of models of human telomeric DNA into G-quadruplex structures

Cation-induced folding into quadruplex structures for three model human telomeric oligonucleotides, d[AGGG(TTAGGG)3], d[TTGGG(TTAGGG)3A] and d[TTGGG(TTAGGG)3], was characterized by equilibrium titrations with KCl and NaCl and by multiwavelength stopped flow kinetics. Cation binding was cooperative with Hill coefficients of 1.5–2.2 in K+ and 2.4–2.9 in Na+ with half-saturation concentrations of 0.5–1 mM for K+ and 4–13 mM for Na+ depending on the oligonucleotide sequence. Oligonucleotide folding in 50 mM KCl at 25°C consisted of single exponential processes with relaxation times τ of 20–60 ms depending on the sequence. In contrast, folding in100 mM NaCl consisted of three exponentials with τ-values of 40–85 ms, 250–950 ms and 1.5–10.5 s. The folding rate constants approached limiting values with increasing cation concentration; in addition, the rates of folding decreased with increasing temperature over the range 15–45°C. Taken together, these results suggest that folding of G-rich oligonucleotides into quadruplex structures proceeds via kinetically significant intermediates. These intermediates may consist of antiparallel hairpins in rapid equilibrium with less ordered structures. The hairpins may subsequently form nascent G-quartets stabilized by H-bonding and cation binding followed by relatively slow strand rearrangements to form the final completely folded topologies. Fewer kinetic intermediates were evident with K+ than Na+, suggesting a simpler folding pathway in K+ solutions.

Hydration is a major determinant of the G-quadruplex stability and conformation of the human telomere 3' sequence of d(AG3(TTAG3)3).

The factors that determine the conformation and stability of G-quadruplex forming sequences remain poorly understood. Here we demonstrate the influence of cosolvents on the conformation and stability of the human telomeric sequence d(A(GGGTTA)3GGG)) in both K(+) and Na(+) containing solutions using a combination of circular dichroism, NMR, and thermodynamics. Molecular crowding arguments have previously been used to suggest that the parallel quadruplex form may be biologically relevant. However, the small cosolvents previously used, PEG 200 and 400, are actually dehydrating agents. We have used acetonitrile as a non-hydrogen-bonding dehydrating agent; similar conformational transitions were observed in K(+) solution. Moreover, NMR analysis shows that the resulting structure contains non-anti guanine glycosyl torsion angles suggesting that the conformation present in acetonitrile is not identical to the all-parallel crystal structure, despite the supposed parallel type CD spectrum.

Molecular Docking of Intercalators and Groove-Binders to Nucleic Acids Using Autodock and Surflex

The molecular docking tools Autodock and Surflex accurately reproduce the crystallographic structures of a collection of small molecule ligands that have been shown to bind nucleic acids. Docking studies were performed with the intercalators daunorubicin and ellipticine and the minor groove binders distamycin and pentamidine. Autodock and Surflex dock daunorubicin and distamycin to their nucleic acid targets within a resolution of approximately 2 A, which is similar to the limit of the crystal structure resolution. However, for the top ranked poses, Autodock and Surflex both dock ellipticine into the correct site but in a different orientation compared to the crystal structure. This appears not only to be partly related to the symmetry of the target nucleic acid, as ellipticine is able to dock from either side of the intercalation site, but also due to the shape of the ligand and docking accuracy. Surflex docks pentamidine in a symmetrically equivalent orientation relative to the crystal structure, while Autodock was able to dock this molecule in the original orientation. In the case of the Surflex docking of pentamidine, the initial rmsd is misleading, given the symmetrical structure of pentamidine. Importantly, the ranking functions of both of these programs are able to return a top pose within approximately 2 A rmsd for daunorubicin, distamycin, and pentamidine and approximately 3 A rmsd for ellipticine compared to their respective crystal structures. Some docking challenges and potential pitfalls are explored, such as the importance of hydrogen treatment on ligands as well as the scoring functions of Autodock and Surflex. Overall for this set of complexes, Surflex is preferred over Autodock for virtual screening, as although the results are comparable, Surflex has significantly faster performance and ease of use under the optimal software conditions tested. These experiments show that molecular docking techniques can be successfully extended to include nucleic acid targets, a finding which has important implications for virtual screening applications and in the design of new small molecules to target therapeutically relevant morphologies of nucleic acids.

The Tail of the Telomere

The structure of a higher-order G-quadruplex structure for human telomeric DNA is presented. The structure was determined by a novel integrated approach in which molecular dynamics simulations were used to produce a stable structure, from which specific experimentally accessible properties were predicted. These properties were tested by sedimentation velocity and steady-state fluorescence measurements. The structure that emerges is a dimeric structure with two quadruplex units, each with a different structure. The interface between the quadruplex units is stabilized by specific stacking interactions of loop nucleotides. The interface is a unique structure and a unique target for drug design.