Tetraplex structure formation in the thrombin-binding DNA aptamer by metal cations measured by vibrational spectroscopy.

Formation of intramolecular tetraplex structures by the thrombin-binding DNA aptamer (TBA) in the presence of K(+), Pb(2+), Ba(2+), Sr(2+) and Mn(2+) has been studied by vibrational spectroscopy. All tetraplex structures contain G-G Hoogsteen type base pairing, both C2'endo/anti and C2'endo/syn deoxyguanosine glycosidic conformations and local B like form DNA phosphate geometries. Addition of Pb(2+) ions modifies the structure by interacting at the level of the guanine carbonyl groups. The very important downshift of the guanine C6=O6 carbonyl vibration mode in the TBA spectrum induced by the addition of one Pb(2+) ion per TBA molecule is in agreement with a localization of the metal ion between both guanine quartets. FTIR melting experiments show an important stabilization of the tetraplex structure upon addition of Pb(2+) ions (DeltaT = 15 degrees C). This strong interaction of lead cations may be correlated with a change in the geometry of the cage formed by the two guanine quartets. A similar but weaker effect is observed for barium and strontium cations.

The three-dimensional structure of the 4:1 mithramycin:d(ACCCGGGT)(2) complex: evidence for an interaction between the E saccharides.

Mithramycin and chromomycin, two antitumor drugs, each having an identical aglycone and nearly identical disaccharide and trisaccharide side chains, have differing binding properties to a small oligonucleotide, d(ACCCGGGT)(2) (M. A. Keniry et al., Journal of Molecular Biology, 1993, Vol. 231, pp. 753-767). In order to understand the forces that induce four mithramycin molecules to bind to d(ACCCGGGT)(2) instead of two drug molecules in the case of chromomycin, the structure of the 4:2:1 mithramycin: Mg(2+):d(ACCCGGGT)(2) complex was investigated by (1)H-nmr and restrained molecular dynamics. The resulting three-dimensional model showed that in order to accommodate the close approach of one neighboring mithramycin dimer, the inwardly directed CDE saccharide chain of the neighboring mithramycin dimer undergoes a conformational change such that the E saccharide no longer spans the minor groove but reorients so that the hydrophilic face of the E saccharides from the two dimers oppose each other. Two hydrogen bonds are formed between the hydroxyl groups of the two opposing E saccharide groups. The results are interpreted in terms of the differences in stereochemistry and functional group substitutions between mithramycin and chromomycin. A mithramycin dimer is able to self-associate on an oligonucleotide template because it has two hydroxyl groups on the same face of its terminal E saccharide. A chromomycin dimer is unable to self-associate because one of these hydroxyl groups is acetylated and the neighboring hydroxyl group has a stereochemistry that cannot permit close contact of the hydroxyl group with a neighbouring chromomycin dimer.

Effect of loop sequence and size on DNA aptamer stability.

The thrombin aptamer is a 15-mer oligodeoxyribonucleotide that folds into a unimolecular quadruplex consisting of a stack of two guanine quartets connected by two external loops and one central loop and possesses a high affinity for thrombin. We have undertaken a systematic examination, in KCl, of the thermodynamic stability of thrombin aptamer analogues containing sequence modifications in one or more of the loops, as well as in the number of quartets. UV melting studies have been carried out to obtain the relevant thermodynamic parameters for these aptamers. van't Hoff analysis of these data, with a two-state model for unimolecular denaturation, gave excellent fits to the experimental observations. Thermodynamic analysis indicates that the central loop sequence in the parent aptamer is optimal for stability. Modifications in this or other loops can effect either DeltaH degrees, DeltaS degrees, or both. Addition of a single G at the 5'-end decreases stability while addition of a G at the 3'-end increases stability. Differential scanning calorimetry experiments on the thrombin aptamer reveal that a heat capacity change, not detected by UV measurements, accompanies the unfolding of the aptamer.

Lead is unusually effective in sequence-specific folding of DNA.

DNA quadruplex structures based on the guanine quartet are typically stabilized by monovalent cations such as K(+), Na(+), or NH(+)(3). Certain divalent cations can also induce quadruplex formation, such as Sr(2+). Here we show that Pb(2+) binds with unusually high affinity to the thrombin binding aptamer, d(GGTTGGTGTGGTTGG), inducing a unimolecular folded structure. At micromolar concentrations the binding is stoichiometric, and a single lead cation suffices to fold the aptamer. The lead-induced changes in UV and CD spectra are characteristic of folded quadruplexes, although the long wavelength CD maximum occurs at 312 nm rather than the typical value of 293 nm. The one-dimensional exchangeable proton NMR spectrum shows resonances expected for imino protons involved in guanine quartet base-pairing. Furthermore, two-dimensional NMR experiments reveal NOE contacts typically seen in folded structures formed by guanine quartets, such as the K(+) form of the thrombin aptamer. Only sequences capable of forming guanine quartets appear to bind Pb(+2) tightly and change conformation. This sequence-specific, tight DNA binding may be relevant to possible genotoxic effects of lead in the environment.

Biological aspects of DNA/RNA quadruplexes.

Among the many unusual conformations of DNA and RNA, quadruplex structures, based on the guanine quartet, possess several unique properties. These properties, along with the general features of guanine quadruplexes, are described in the context of possible roles for these structures in biological systems. A variety of experimental observations supporting the notion that quadruplexes are important in vivo is presented, including proteins known to specifically bind to quadruplex structures, guanine-rich DNA, and RNA sequences endowed with the potential for forming quartet-based structures in telomeres and regulatory regions, such as gene promoters, quadruplexes as DNA aptamer folding motifs arising from in vitro selection experiments, and potential chemotherapeutic, quadruplex-forming oligonucleotides. Taken together, all of these observations argue cogently not only for the presence of quadruplexes in biological systems but also for their significance in terms of their roles in various biological processes.

NMR structure refinement and dynamics of the K+-[d(G3T4G3)]2 quadruplex via particle mesh Ewald molecular dynamics simulations.

The solution structure and dynamical properties of the potassium-stabilized, hairpin dimer quadruplex formed by the oligonucleotide d(G3T4G3) have been elucidated by a combination of high-resolution NMR and molecular dynamics simulations. Refinement calculations were carried out both in vacuo, without internally coordinated K+ cations, and in explicit water, with internally coordinated K+ cations. In the latter case, the electrostatic interactions were calculated using the particle mesh Ewald (PME) method. The NMR restraints indicate that the K+ quadruplex has a folding arrangement similar to that formed by the same oligonucleotide in the presence of sodium, but with significant local differences. Unlike the Na+ quadruplex, the thymine loops found in K+ exhibit considerable flexibility, and appear to interconvert between two preferred conformations. Furthermore, the NMR evidence points toward K+-stabilized guanine quartets of slightly larger diameter relative to the Na+-stabilized structure. The characteristics of the quartet stem are greatly affected by the modeling technique employed: caged cations alter the size and symmetry of the quartets, and explicit water molecules form hydration spines within the grooves. These results provide insight into those factors that determine the overall stability of hairpin dimer quadruplexes and the effects of different cations in modulating the relative stability of the dimeric hairpin and linear, four-stranded, quadruplex forms.

Stability and structure of model DNA triplexes and quadruplexes and their interactions with small ligands.

This review focuses on the structural and thermodynamic characterization of model DNA triplex and quadruplex structures, taking into account effects of stoichiometry and sequence. Methods such as gel electrophoresis, UV melting, and scanning calorimetry, and the results thereof, are described for determination of the thermodynamic stability of such systems. Three classes of triplexes are considered based on the composition of the third strand, while quadruplex systems are limited to those based on the guanine quartet. X-ray crystallography and high resolution NMR studies are also described for these two classes of unusual structures. Ligand binding to triplexes and quadruplexes is also reviewed, with emphasis on specific molecular recognition. The availability of three-dimensional structures for triplex and quadruplex species sets the stage for structure-based development of ligands capable of binding to them specifically. To this end, we consider the application of DOCK, a program for the discovery of small molecules that can recognize macromolecular structures, to the problem of recognizing folded quadruplex structures. Such studies may ultimately lead to pharmaceutically active compounds.

The contribution of thymine-thymine interactions to the stability of folded dimeric quadruplexes.

The loop of four thymines in the sodium form of the dimeric folded quadruplex [d(G3T4G3)]2 assumes a well-defined structure in which hydrogen bonding between the thymine bases appears to contribute to the stability and final conformation of the quadruplex. We have investigated the importance of the loop interactions by systematically replacing each thymine in the loop with a cytosine. The quadruplexes formed by d(G3CT3G3), d(G3TCT2G3), d(G3T2CTG3) and d(G3T3CG3) in the presence of 150 mM Na+ were studied by gel mobility, circular dichroism and 1H NMR spectroscopy. The major species formed by d(G3CT3G3), d(G3TCT2G3) and d(G3T3CG3) at 1 mM strand concentration at neutral pH is a dimeric folded quadruplex. d(G3T2CTG3) has anomalous behaviour and associates into a greater percentage of linear four-stranded quadruplex than the other three oligonucleotides at neutral pH and at the same concentration. The linear four-stranded quadruplex has a greater tendency to oligomerize to larger ill-defined structures, as demonstrated by broad 1H NMR resonances. At pH 4, when the cytosine is protonated, there is a greater tendency for each of the oligonucleotides to form some four-stranded linear quadruplex, except for d(G3T2CTG3), which has the reverse tendency. The experimental results are discussed in terms of hydrogen bonding within the thymine loop.

Structure-based discovery of ligands targeted to the RNA double helix.

Ligands capable of specific recognition of RNA structures are of interest in terms of the principles of molecular recognition as well as potential chemotherapeutic applications. We have approached the problem of identifying small molecules with binding specificity for the RNA double helix through application of the DOCK program [Kuntz, I. D., Meng, E. C., and Shoichet, B. K. (1994) Acc. Chem. Res. 27, 117-123], a structure-based method for drug discovery. A series of lead compounds was generated through a database search for ligands with shape complementarity to the RNA deep major groove. Compounds were then evaluated with regard to their fit into the minor groove of B DNA. Those compounds predicted to have an optimal fit to the RNA groove and strong discrimination against DNA were examined experimentally. Of the 11 compounds tested, 3, all aminoglycosides, exhibited pronounced stabilization of RNA duplexes against thermal denaturation with only marginal effects on DNA duplexes. One compound, lividomycin, was examined further, and shown to facilitate the ethanol-induced B to A transition in calf thymus DNA. Fluorine NMR solvent isotope shift measurements on RNA duplexes containing 5-fluorouracil provided evidence that lividomycin binds in the RNA major groove. Taken together, these results indicate that lividomycin recognizes the general features of the A conformation of nucleic acids through deep groove binding, confirming the predictions of our DOCK analysis. This approach may be of general utility for identifying ligands possessing specificity for additional RNA structures as well as other nucleic acid structural motifs.

Studies on formation and stability of the d[G(AG)5]* d[G(AG)5]. d[C(TC)5] and d[G(TG)5]* d[G(AG)5]. d[C(TC)5] triple helices.

We have targeted the d[G(AG)5]. d[C(TC)5] duplex for triplex formation at neutral pH with either d[G(AG)5] or d[G(TG)5]. Using a combination of gel electrophoresis, uv and CD spectra, mixing and melting curves, along with DNase I digestion studies, we have investigated the stability of the 2:1 pur*pur.pyr triplex, d[G(AG)5]*d[G(AG)5].d[C(TC)5], in the presence of MgCl2. This triplex melts in a monophasic fashion at the same temperature as the underlying duplex. Although the uv spectrum changes little upon binding of the second purine strand, the CD spectrum shows significant changes in the wavelength range 200-230 nm and about a 7 nm shift in the positive band near 270 nm. In contrast, the 1:1:1 pur/pyr*pur.pyr triplex, d[G(TG)5]*d[G(AG)5].d[C(TC)5], is considerably less stable thermally, melting at a much lower temperature than the underlying duplex, and possesses a CD spectrum that is entirely negative from 200 to 300 nm. Ethidium bromide undergoes a strong fluorescence enhancement upon binding to each of these triplexes, and significantly stabilizes the pur/pyr*pur.pyr triplex. The uv melting and differential scanning calorimetry analysis of the alternating sequence duplex and pur*pur.pyr triplex shows that they are lower in thermodynamic stability than the corresponding 10-mer d(G3A4G3). d(C3T4C3) duplex and its pur*pur.pyr triplex under identical solution conditions.

Design of sequence-specific DNA binding ligands that use a two-stranded peptide motif for DNA sequence recognition.

The design and DNA binding activity of beta-structure-forming peptides and netropsin-peptide conjugates are reported. It is found that a pair of peptides-S,S'-bis(Lys-Gly-Val-Cys-Val-NH-NH-Dns)-bridged by an S-S bond binds at least 10 times more strongly to poly(dG).poly(dC) than to poly(dA).poly(dT). This peptide can also discriminate between 5'-GpG-3' and 5'-GpC-3' steps in the DNA minor groove. Based on these observations, new synthetic ligands, bis-netropsins, were constructed in which two netropsin-like fragments were attached by means of short linkers to a pair of peptides-Gly-Cys-Gly- or Val-Cys-Val-bridged by S-S bonds. These compounds possess a composite binding specificity: the peptide chains recognize 5'-GpG-3' steps on DNA, whereas the netropsin-like fragments bind preferentially to runs of 4 AT base pairs. Our data indicate that combining the AT-base-pair specific properties of the netropsin-type structure with the 5'-GpG-3'-specific properties of certain oligopeptides offers a new approach to the synthesis of ligands capable of recognizing mixed sequences of AT- and GC-base pairs in the DNA minor groove. These compounds are potential models for DNA-binding domains in proteins which specifically recognize base pair sequences in the minor groove of DNA.

Calorimetric analysis of triple helices targeted to the d(G3A4G3).d(C3T4C3) duplex.

We present a thermodynamic analysis based on differential scanning calorimetry (DSC) of three short intermolecular DNA triplexes targeted to the same DNA duplex: d(C+3T4C+3)*d- (G3A4G3).d(C3T4C3) (PYR), d(G3A4G3)*d(G3A4G3).d(C3T4C3) (PUR), and d(G3T4G3)*d(G3A4G3).d(C3T4C3) (PUR/PYR). Enthalpies, delta H, and entropies, delta S, are measured by model-free integration of the DSC curves and are compared to the same quantities determined by van't Hoff analysis of the DSC curves and, in the case of the PYR and PUR/PYR triplexes, UV melting curves as well. In the case of the PUR triplex, which exhibits monophasic melting behavior, the calorimetric delta H and the calorimetrically determined van't Hoff delta H are in excellent agreement, indicating an all-or-none transition for this triplex. For the PYR and PUR/PYR triplexes, which melt in a biphasic manner, the calorimetrically determined van't Hoff delta H values are somewhat larger than the model-independent calorimetric delta H values. In those cases, however, good agreement is found between the calorimetric delta H values and the spectrophotometrically determined van't Hoff delta H values. The calorimetrically determined delta H values, expressed per mole of triplet, for the three triplexes are 4.5, 3.8, and 2.4 kcal/mol for the PUR, PYR, and PUR/PYR triplexes, respectively. The same order of stability is observed in terms of delta G and Tm values. The high stability of the PUR triplex at neutral pH indicates that purine oligonucleotides may be the most effective at targeting duplex regions for triple helix formation in vivo.

Spectroscopic recognition of guanine dimeric hairpin quadruplexes by a carbocyanine dye.

Isolated guanine quadruplex structures have been described at high resolution both in solution and in the solid state. The existence of this unusual DNA structure in vivo and its biological significance remain to be determined. We describe the binding of 3,3'-diethyloxadicarbocyanine to dimeric hairpin guanine quadruplexes. This interaction results in a set of unique spectrophotometric signatures, none of which arises from binding to single strands or Watson-Crick duplexes. These unique signatures include a new absorbance peak (lambda max = 534 nm), an induced circular dichroism (lambda = 534-626 nm), a quenching of the dye fluorescence upon excitation with visible light, and strong energy transfer from DNA. This last effect provides the basis for detecting hairpin quadruplex structures in the presence of excess amounts of nonquadruplex DNA structures, such as single strands and Watson-Crick duplexes. The mechanism of quadruplex recognition by this dye is discussed, along with the possibility of using this dye as a probe for hairpin quadruplex structures in vitro and in vivo.

Mono-, di- and trimeric binding of a bis-netropsin to DNA.

An unusual 3:1 stoichiometry for complex formation between an elongated bis-netropsin compound and its binding site on DNA has been observed. Circular dichroism measurements distinguish two types of complexes formed between this bis-netropsin and poly[d(A-T)].poly[d(A-T)]. The first type is characterized by a 1:1 saturating ratio of bound molecules per ten base pairs. Formation of the second type results from the cooperative binding of two additional bis-netropsin molecules to the first type of complex. In contrast to these results observed for binding to the alternating polynucleotide, only the 1:1 type of complex is formed when this ligand binds to the homopolymer poly(dA).poly(dT).

Solution structure of the Na+ form of the dimeric guanine quadruplex [d(G3T4G3)]2.

The solution structure of the DNA quadruplex formed by the association of two strands of the DNA oligonucleotide, d(G3T4G3), in NaCl solution has been determined by 1H two-dimensional NMR techniques, full relaxation matrix calculations and restrained molecular dynamics. The refined structure incorporates the sequences 5'-G1sG2AG3AT4AT5AT6AT7AG8sG9AG10A-3' and 5'-G11sG12AG13AT14AT15AT16AT17AG18sG19sG20A-3' (where S and A denote syn and anti, respectively) in a three-quartet, diagonal-looped structure that we [Strahan, G. D., Shafer, R. H. & Keniry, M. A. (1994) Nucleic Acids Res. 22, 5447-5455] and others [Smith, F. W., Lau, F. W. & Feigon, J. (1994) Proc. Natl. Acad. Sci. USA 91, 10546-10550] have described. The loop structure is compact and incorporates many of the features found in duplex hairpin loops including base stacking, intraloop hydrogen bonding and extensive van der Waals' interactions. The first and third loop thymines stack over the outermost G-quartet and are also associated by hydrogen bonding. The second and the fourth loop thymines fold inwards in order to enhance van der Waals' interactions. The unexpected sequential syn-syn deoxyguanosines in the quadruplex stem appear to be a direct consequence of the way DNA oligonucleotides fold and the subsequent search for the most stable loop structure. The implications of loop sequence and length on the structure of quadruplexes are discussed.

Physicochemical studies of the d(G3T4G3)*d(G3A4G3).d(C3T4C3) triple helix.

We have targeted the d(G3A4G3).d(C3T4C3) duplex for triplex formation with d(G3T4G3) in the presence of MgCl2. The resulting triple helix, d(G3T4G3)*d(G3-A4G3).d(C3T4C3), is considerably weaker than the related triplex, d(G3A4G3)*d(G3A4G3).d(C3T4C3), and melts in a biphasic manner, with the third strand dissociating at temperatures about 20-30 degrees C below that of the remaining duplex. This is in distinct contrast to the d(G3A4G3)*d(G3A4G3).d(C3T4C3) triplex, which melts in essentially a single transition. Gel electrophoresis under non-denaturing conditions shows the presence of the d(G3T4G3)*d(G3A4G3).d(C3T4C3) triplex as a band of low mobility compared to the duplex or the single strand bands. Binding of the d(G3T4G3) third strand and the purine strand of the duplex can be monitored by imino proton NMR spectra. While these spectra are typically very broad for intermolecular triplexes, the line widths can be dramatically narrowed by the addition of two thymines to both ends of the pyrimidine strand. Thermodynamic analysis of UV melting curves shows that this triplex is considerably less stable than related triplexes formed with the same duplex. The orientation of the third strand was addressed by a combination of fluorescence energy transfer and UV melting experiments. Results from these experiments suggest that, in the unlabeled triplex, the preferred orientation of the third strand is parallel to the purine strand of the duplex.

NMR studies of drug-DNA complexes.

The application of high-resolution, multidimensional NMR techniques to the problem of determining the structure of drug-DNA complexes in solution has led to substantial progress in understanding the effect of drugs on DNA at the molecular level. With the development of isotopic labeling methods applied in three- and four-dimensional experiments, we anticipate that more complex drug-DNA systems will become amenable to structural analysis. In addition to implementing these newer techniques, progress will also be made in terms of investigating the structure of drug complexes with more unusual forms of DNA, such as triplexes, quadruplexes, multistranded junctions, and so forth.

Structural properties of the [d(G3T4G3)]2 quadruplex: evidence for sequential syn-syn deoxyguanosines.

Two-dimensional 1H NMR studies on the dimeric hairpin quadruplex formed by d(G3T4G3) in the presence of either NaCl or KCl are presented. In the presence of either salt, the quadruplex structure is characterized by half the guanine nucleosides in the syn conformation about the glycosidic bond, the other half in the anti conformation, as reported for other similar sequences. However, 1H NOESY and 1H-31P heteronuclear correlation experiments demonstrate that the deoxyguanosines do not strictly alternate between syn and anti along individual strands. Thus we find the following sequences with regard to glycosidic bond conformation: 5'-G1SG2SG3AT4AT5A-T6AT7AG8SG9AG10A-3' and 5'-G11SG12AG13AT14AT1 5AT16AT17AG18SG19SG20A-3', where S and A denote syn and anti, respectively. This represents the first experimental evidence for a nucleic acid structure containing two sequential nucleosides in the syn conformation. The stacking interactions of the resulting quadruplex quartets and their component bases have been evaluated using unrestrained molecular dynamics calculations and energy component analysis. These calculations suggest that the sequential syn-syn/anti-anti and syn-anti quartet stacks are almost equal in energy, whereas the anti-syn stack, which is not present in our structure, is energetically less favorable by about 4 kcal/mol. Possible reasons for this energy difference and its implications for the stability of quadruplex structures are discussed.

Solution structure of the octamer motif in immunoglobulin genes via restrained molecular dynamics calculations.

The solution structure of the DNA decamer d(CATTTGCATC)-d(GATGCAAATG), comprising the octamer motif of immunoglobulin genes, is determined by restrained molecular dynamics (rMD) simulations. The restraint data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were previously obtained by a complete relaxation matrix analysis of the two-dimensional nuclear Overhauser enhancement (2D NOE) intensities and by the quantitative simulation of cross-peaks in double-quantum-filtered correlated (2QF-COSY) spectra. The influence of torsion angles and the number of experimental distance restraints on the structural refinement has been systematically examined. Omitting part of the experimental NOE-derived distances results in reduced restraint violations and lower R factors but impairs structural convergence in the rMD refinement. Eight separate restrained molecular dynamics simulations were carried out for 20 ps each, starting from either energy-minimized A- or B-DNA. Mutual atomic root-mean-square (rms) differences among the refined structures are well below 1 A and comparable to the rms fluctuations of the atoms about their average position, indicating convergence to essentially identical structures. The average refined structure was subjected to an additional 100 ps of rMD simulations and analyzed in terms of average torsion angles and helical parameters. The B-type duplex exhibits clear sequence-dependent variations in its geometry with a narrow minor groove at the T3.A3 tract and a large positive roll at the subsequent TG.CA step. This is accompanied by a noticeable bend of the global helix axis into the major groove. There is also evidence of significant flexibility of the sugar-phosphate backbone with rapid interconversion among different conformers.

Differential binding of the enantiomers of chloroquine and quinacrine to polynucleotides: implications for stereoselective metabolism.

Interaction of the antimalarial drugs quinacrine and chloroquine with DNA has been studied extensively in order to understand the origin of their biological activity. These studies have shown that they bind to DNA through an intercalative mode and show little sequence specificity. All previous experiments were carried out using the racemic form of these drugs. We have investigated the binding of the enantiomeric forms of quinacrine and chloroquine to synthetic polynucleotides poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly (dG-dC), and found interesting differences in their binding parameters. Quinacrine enantiomers have a much higher binding affinity for the two polynucleotides compared to those of chloroquine. The negative enantiomers were found to have higher binding affinity than the positive ones. The binding constant for the binding of quinacrine(-) to poly(dG-dC).poly(dG-dC) was found to be about 3 times that of quinacrine(+). The differences in these binding affinities were further confirmed by equilibrium dialysis of the complexes of the polynucleotides with the racemic form of the drugs, which resulted in the enrichment of the dialysate with the positive enantiomer. CD spectra of the enantiomers and their polynucleotide complexes are reported. Changes in the fluorescence properties of quinacrine in the presence of the two polynucleotides are also described. Biological implications of these findings are discussed.