Interactions of 4'', 6-diamidine-2-phenylindole with synthetic polynucleotides.

4', 6-Diamidine-2-phenylindole forms fluorescent complexes with synthetic DNA duplexes containing AT, AU and IC base pairs; no fluorescent complexes were observed with duplexes containing GC base pairs or with duplexes containing a single AT base pair sandwiched between GC pairs. The binding site size is one molecule of dye per 3 base pairs. The intrinsic binding constants are higher for alternating sequence duplexes than for the corresponding homopolymer pairs. With the exception of the four-stranded helical poly rI which exhibits considerable fluorescence enhancement upon binding of the ligand, none of the single- or multi- stranded polyribonucleotides and ribo-deoxyribonucleotide hybrid structures form fluorescent complexes with the dye. Poly rI is the only RNA which forms a DNA B-like structure (Arnott et al. (1974) Biochem. J. 141, 537). The B conformation of the helix and the absence of guanine appear to be the major determinants of the specificity of the fluorescent binding mode of the dye. Nonfluorescent interactions of the dye with polynucleotides are nonspecific; UV absorption and circular dichroic spectra demonstrate binding to synthetic single- and double-stranded DNA and RNA analogs, including those containing GC base pairs.     

Solution structure of a steroid-DNA complex with cholic acid residues sealing the termini of a Watson-Crick duplex

The three-dimensional structure of a covalent hybrid between cholic acid and the self-complementary DNA hexamer 5'-TGCGCA-3' was solved via two-dimensional NMR and restrained torsion angle molecular dynamics. In the complex, refined to a pairwise rmsd of 0.64 A, the steroid binds to the terminal T:A base pairs via extensive van der Waals contacts but without any hydrogen bonds or detectable dipole-dipole interactions. The contacts involve the methyl groups as well as one edge of the streoid's sterane skeleton and both nucleobases and the deoxyriboses of the terminal base pair of the DNA. The surprising shape complementarity between steroid and the undisturbed DNA termini explains the increase in fidelity and affinity observed for hybridization probes bearing bile acid residues. Since the hydroxyl groups of the steroid do not contribute to the binding of the DNA, they may be derivatized, potentially giving access to a new class of specific binders for blunt ends of Watson-Crick duplexes.     

The stability and number of nucleating interactions determine DNA hybridization rates in the absence of secondary structure

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson–Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.     

End-to-end attraction of duplex DNA

Recent experiments [Nakata, M. et al., End-to-end stacking and liquid crystal condensation of 6 to 20 basepair DNA duplexes. Science 2007; 318:1276-1279] have demonstrated spontaneous end-to-end association of short duplex DNA fragments into long rod-like structures. By means of extensive all-atom molecular dynamic simulations, we characterized end-to-end interactions of duplex DNA, quantitatively describing the forces, free energy and kinetics of the end-to-end association process. We found short DNA duplexes to spontaneously aggregate end-to-end when axially aligned in a small volume of monovalent electrolyte. It was observed that electrostatic repulsion of 5'-phosphoryl groups promoted the formation of aggregates in a conformation similar to the B-form DNA double helix. Application of an external force revealed that rupture of the end-to-end assembly occurs by the shearing of the terminal base pairs. The standard binding free energy and the kinetic rates of end-to-end association and dissociation processes were estimated using two complementary methods: umbrella sampling simulations of two DNA fragments and direct observation of the aggregation process in a system containing 458 DNA fragments. We found the end-to-end force to be short range, attractive, hydrophobic and only weakly dependent on the ion concentration. The relation between the stacking free energy and end-to-end attraction is discussed as well as possible roles of the end-to-end interaction in biological and nanotechnological systems.     

Hybridization kinetics of out-of-equilibrium mixtures of short RNA oligonucleotides

Hybridization and strand displacement kinetics determine the evolution of the base paired configurations of mixtures of oligonucleotides over time. Although much attention has been focused on the thermodynamics of DNA and RNA base pairing in the scientific literature, much less work has been done on the time dependence of interactions involving multiple strands, especially in RNA. Here we provide a study of oligoribonucleotide interaction kinetics and show that it is possible to calculate the association, dissociation and strand displacement rates displayed by short oligonucleotides (5nt–12nt) that exhibit no expected secondary structure as simple functions of oligonucleotide length, CG content, ΔG of hybridization and ΔG of toehold binding. We then show that the resultant calculated kinetic parameters are consistent with the experimentally observed time dependent changes in concentrations of the different species present in mixtures of multiple competing RNA strands. We show that by changing the mixture composition, it is possible to create and tune kinetic traps that extend by orders of magnitude the typical sub-second hybridization timescale of two complementary oligonucleotides. We suggest that the slow equilibration of complex oligonucleotide mixtures may have facilitated the nonenzymatic replication of RNA during the origin of life.     

An intercalation-locked parallel-stranded DNA tetraplex

DNA has proved to be an excellent material for nanoscale construction because complementary DNA duplexes are programmable and structurally predictable. However, in the absence of Watson–Crick pairings, DNA can be structurally more diverse. Here, we describe the crystal structures of d(ACTCGGATGAT) and the brominated derivative, d(AC^BrUCGGA^BrUGAT). These oligonucleotides form parallel-stranded duplexes with a crystallographically equivalent strand, resulting in the first examples of DNA crystal structures that contains four different symmetric homo base pairs. Two of the parallel-stranded duplexes are coaxially stacked in opposite directions and locked together to form a tetraplex through intercalation of the 5′-most A–A base pairs between adjacent G–G pairs in the partner duplex. The intercalation region is a new type of DNA tertiary structural motif with similarities to the i-motif. ¹H–¹H nuclear magnetic resonance and native gel electrophoresis confirmed the formation of a parallel-stranded duplex in solution. Finally, we modified specific nucleotide positions and added d(GAY) motifs to oligonucleotides and were readily able to obtain similar crystals. This suggests that this parallel-stranded DNA structure may be useful in the rational design of DNA crystals and nanostructures.     

Pyrene is highly emissive when attached to the RNA duplex but not to the DNA duplex: the structural basis of this difference

Through binding and fluorescence studies of oligonucleotides covalently attached to a pyrene group via one carbon linker at the sugar residue, we previously found that pyrene-modified RNA oligonucleotides do not emit well in the single-stranded form, yet the attached pyrene emits with a significantly high quantum yield upon binding to a complementary RNA strand. In sharp contrast, similarly modified pyrene–DNA probes exhibit very weak fluorescence both in the double-stranded and single-stranded forms. The pyrene-modified RNA oligonucleotides therefore provide a useful tool for monitoring RNA hybridization. The purpose of this paper is to present the structural basis for the different fluorescence properties of pyrene-modified RNA/RNA and pyrene-modified DNA/DNA duplexes. The results of absorption, fluorescence anisotropy and circular dichroism studies all consistently indicated that the pyrene attached to the RNA duplex is located outside of the duplex, whereas the pyrene incorporated into the DNA duplex intercalates into the double helix. ¹H NMR measurements unambiguously confirmed that the pyrene attached to the DNA duplex indeed intercalates between the base pairs of the duplex. Molecular dynamics simulations support these differences in the local structural elements around the pyrene between the pyrene–RNA/RNA and the pyrene–DNA/DNA duplexes.     

Hybridization of 2'-ribose modified mixed-sequence oligonucleotides: thermodynamic and kinetic studies

In this study, we characterize the thermodynamics of hybridization, binding kinetics and conformations of four ribose-modified (2′-fluoro, 2′-O-propyl, 2′-O-methoxyethyl and 2′-O-aminopropyl) decameric mixed-sequence oligonucleotides. Hybridization to the complementary non-modified DNA or RNA decamer was probed by fluorescence and circular-dichroism spectroscopy and compared to the same duplex formed between two non-modified strands. The thermal melting points of DNA–DNA duplexes were increased by 1.8, 2.2, 0.3 and 1.3°C for each propyl, methoxyethyl, aminopropyl and fluoro modification, respectively. In the case of DNA–RNA duplexes, the melting points were increased by 3.1, 4.1 and 1.0°C for each propyl, methoxyethyl and aminopropyl modification, respectively. The high stability of the duplexes formed with propyl-, methoxyethyl- and fluoro-modified oligonucleotides correlated with high preorganization in these single-strands. Despite higher thermodynamic duplex stability, hybridization kinetics to complementary DNA or RNA was slower for propyl- and methoxyethyl-modified oligonucleotides than for the non-modified control. In contrast, the positively-charged aminopropyl-modified oligonucleotide showed rapid binding to the complementary DNA or RNA.     

Thermodynamics of sequence-specific binding of PNA to DNA

For further characterization of the hybridization properties of peptide nucleic acids (PNAs), the thermodynamics of hybridization of mixed sequence PNA-DNA duplexes have been studied. We have characterized the binding of PNA to DNA in terms of binding affinity (perfectly matched duplexes) and sequence specificity of binding (singly mismatched duplexes) using mainly absorption hypochromicity melting curves and isothermal titration calorimetry. For perfectly sequence-matched duplexes of varying lengths (6-20 bp), the average free energy of binding (DeltaG degrees ) was determined to be -6.5+/-0.3 kJ mol(-1) bp(-1), corresponding to a microscopic binding constant of about 14 M(-1) bp(-1). A variety of single mismatches were introduced in 9- and 12-mer PNA-DNA duplexes. Melting temperatures (T(m)) of 9- and 12-mer PNA-DNA duplexes with a single mismatch dropped typically 15-20 degrees C relative to that of the perfectly matched sequence with a corresponding free energy penalty of about 15 kJ mol(-1) bp(-1). The average cost of a single mismatch is therefore estimated to be on the order of or larger than the gain of two matched base pairs, resulting in an apparent binding constant of only 0.02 M(-1) per mismatch. The impact of a mismatch was found to be dependent on the neighboring base pairs. To a first approximation, increasing the stability of the surrounding region, i.e., the distribution of A.T and G.C base pairs, decreases the effect of the introduced mismatch.     

Massive parallel analysis of DNA-Hoechst 33258 binding specificity with a generic oligodeoxyribonucleotide microchip.

A generic oligodeoxyribonucleotide microchip was used to determine the sequence specificity of Hoechst 33258 binding to double-stranded DNA. The generic microchip contained 4096 oxctadeoxynucleo-tides in which all possible 4(6)= 4096 hexadeoxy-nucleotide sequences are flanked on both the 3'- and 5'-ends with equimolar mixtures of four bases. The microchip was manufactured by chemical immobilization of presynthesized 8mers within polyacrylamide gel pads. A selected set of immobilized 8mers was converted to double-stranded form by hybridization with a mixture of fluorescently labeled complementary 8mers. Massive parallel measurements of melting curves were carried out for the majority of 2080 6mer duplexes, in both the absence and presence of the Hoechst dye. The sequence-specific affinity for Hoechst 33258 was calculated as the increase in melting temperature caused by ligand binding. The dye exhibited specificity for A:T but not G:C base pairs. The affinity is low for two A:T base pairs, increases significantly for three, and reaches a plateau for four A:T base pairs. The relative ligand affinity for all trinucleotide and tetranucleotide sequences (A/T)(3)and (A/T)(4)was estimated. The free energy of dye binding to several duplexes was calculated from the equilibrium melting curves of the duplexes formed on the oligonucleotide microchips. This method can be used as a general approach for massive screening of the sequence specificity of DNA-binding compounds.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.     

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2′-O-methyl (2′-OMe) oligonucleotides GCGAAUUCGC, (UA)₆ and (CG)₃ using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Δn_W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Δn_W was 3–4 for RNA duplexes and 2–3 for DNA and 2′-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2′-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2′-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.     

Massive parallel analysis of the binding specificity of histone-like protein HU to single- and double-stranded DNA with generic oligodeoxyribonucleotide microchips

A generic hexadeoxyribonucleotide microchip has been applied to test the DNA-binding properties of HU histone-like bacterial protein, which is known to have a low sequence specificity. All 4096 hexamers flanked within 8mers by degenerate bases at both the 3′- and 5′-ends were immobilized within the 100 × 100 × 20 mm polyacrylamide gel pads of the microchip. Single-stranded immobilized oligonucleotides were converted in some experiments to the double-stranded form by hybridization with a specified mixture of 8mers. The DNA interaction with HU was characterized by three type of measurements: (i) binding of FITC-labeled HU to microchip oligonucleotides; (ii) melting curves of complexes of labeled HU with single-stranded microchip oligonucleotides; (iii) the effect of HU binding on melting curves of microchip double-stranded DNA labeled with another fluorescent dye, Texas Red. Large numbers of measurements of these parameters were carried out in parallel for all or many generic microchip elements in real time with a multi-wavelength fluorescence microscope. Statistical analysis of these data suggests some preference for HU binding to G/C-rich single-stranded oligonucleotides. HU complexes with double-stranded microchip 8mers can be divided into two groups in which HU binding either increased the melting temperature (T_m) of duplexes or decreased it. The stabilized duplexes showed some preference for presence of the sequence motifs AAG, AGA and AAGA. In the second type of complex, enriched with A/T base pairs, the destabilization effect was higher for longer stretches of A/T duplexes. Binding of HU to labeled duplexes in the second type of complex caused some decrease in fluorescence. This decrease also correlates with the higher A/T content and lower T_m. The results demonstrate that generic microchips could be an efficient approach in analysis of sequence specificity of proteins.     

Incorporation of thio-pseudoisocytosine into triplex-forming peptide nucleic acids for enhanced recognition of RNA duplexes

Peptide nucleic acids (PNAs) have been developed for applications in biotechnology and therapeutics. There is great potential in the development of chemically modified PNAs or other triplex-forming ligands that selectively bind to RNA duplexes, but not single-stranded regions, at near-physiological conditions. Here, we report on a convenient synthesis route to a modified PNA monomer, thio-pseudoisocytosine (L), and binding studies of PNAs incorporating the monomer L. Thermal melting and gel electrophoresis studies reveal that L-incorporated 8-mer PNAs have superior affinity and specificity in recognizing the duplex region of a model RNA hairpin to form a pyrimidine motif major-groove RNA₂–PNA triplex, without appreciable binding to single-stranded regions to form an RNA–PNA duplex or, via strand invasion, forming an RNA–PNA₂ triplex at near-physiological buffer condition. In addition, an L-incorporated 8-mer PNA shows essentially no binding to single-stranded or double-stranded DNA. Furthermore, an L-modified 6-mer PNA, but not pseudoisocytosine (J) modified or unmodified PNA, binds to the HIV-1 programmed −1 ribosomal frameshift stimulatory RNA hairpin at near-physiological buffer conditions. The stabilization of an RNA₂–PNA triplex by L modification is facilitated by enhanced van der Waals contacts, base stacking, hydrogen bonding and reduced dehydration energy. The destabilization of RNA–PNA and DNA–PNA duplexes by L modification is due to the steric clash and loss of two hydrogen bonds in a Watson–Crick-like G–L pair. An RNA₂–PNA triplex is significantly more stable than a DNA₂–PNA triplex, probably because the RNA duplex major groove provides geometry compatibility and favorable backbone–backbone interactions with PNA. Thus, L-modified triplex-forming PNAs may be utilized for sequence-specifically targeting duplex regions in RNAs for biological and therapeutic applications.     

Mismatching base-pair dependence of the kinetics of DNA-DNA hybridization studied by surface plasmon fluorescence spectroscopy

Two single-stranded DNAs consisting of complementary base pairs except for one mismatching base pair (MM1) can form double-stranded DNA by molecular recognition. This type of duplex is not as stable as that formed by MM0. In order to add to a better understanding of the physical mechanism of the hybridization and dissociation processes at sensor (chip) surfaces, we studied the kinetics of the MM1 hybridization by surface plasmon fluorescence spectroscopy. Target DNA strands labelled with a fluorescent molecule Cy5 at the 5′ end and hybridizing with the surface-attached probe DNA can be excited by the strong optical field of a surface plasmon resonance mode. The emitted fluorescence can be detected with high sensitivity. The affinity of a duplex was found to depend on the chemical nature, i.e. G–G, G–T etc., and on the position of the mismatching base pair along the 15mer duplex.     

Simultaneous DNA binding, bending, and base flipping: evidence for a novel M.EcoRI methyltransferase-DNA complex

We measured the kinetics of DNA bending by M.EcoRI using DNA labeled at both 5'-ends and observed changes in fluorescence resonance energy transfer. Although known to bend its cognate DNA site, energy transfer is decreased upon enzyme binding. This unanticipated effect is shown to be robust because we observe the identical decrease with different dye pairs, when the dye pairs are placed on the respective 3'-ends, the effect is cofactor- and protein-dependent, and the effect is observed with duplexes ranging from 14 through 17 base pairs. The same labeled DNA shows the anticipated increased energy transfer with EcoRV endonuclease, which also bends this sequence, and no change in energy transfer with EcoRI endonuclease, which leaves this sequence unbent. We interpret these results as evidence for an increased end-to-end distance resulting from M.EcoRI binding, mediated by a mechanism novel for DNA methyltransferases, combining DNA bending and an overall expansion of the DNA duplex. The M.EcoRI protein sequence is poorly accommodated into well defined classes of DNA methyltransferases, both at the level of individual motifs and overall alignment. Interestingly, M.EcoRI has an intercalation motif observed in the FPG DNA glycosylase family of repair enzymes. Enzyme-dependent changes in anisotropy and fluorescence resonance energy transfer have similar rate constants, which are similar to the previously determined rate constant for base flipping; thus, the three processes are nearly coincidental. Similar fluorescence resonance energy transfer experiments following AdoMet-dependent catalysis show that the unbending transition determines the steady state product release kinetics.     

Design of LNA probes that improve mismatch discrimination

Locked nucleic acids (LNA) show remarkable affinity and specificity against native DNA targets. Effects of LNA modifications on mismatch discrimination were studied as a function of sequence context and identity of the mismatch using ultraviolet (UV) melting experiments. A triplet of LNA residues centered on the mismatch was generally found to have the largest discriminatory power. An exception was observed for G–T mismatches, where discrimination decreased when the guanine nucleotide at the mismatch site or even the flanking nucleotides were modified. Fluorescence experiments using 2-aminopurine suggest that LNA modifications enhance base stacking of perfectly matched base pairs and decrease stabilizing stacking interactions of mismatched base pairs. LNAs do not change the amount of counterions (Na⁺) that are released when duplexes denature. New guidelines are suggested for design of LNA probes, which significantly improve mismatch discrimination in comparison with unmodified DNA probes.     

Use of hybridization melting kinetics for detecting Phytophthora species using three-dimensional microarrays: demonstration of a novel concept for the differentiation of detection targets

Microarray-based detection is limited by variable and inconsistent hybridization intensities across the diversity of probes used in each array. In this paper, we introduce a novel concept for the differentiation of detection targets using duplex melting kinetics. A microarray assay was developed on a PamChip microarray enabling the differentiation of target Phytophthora species using the melting kinetics of probe-target duplexes. In the majority of cases the hybridization kinetics of target and non-target duplexes differed significantly. Analysis of the melting kinetics of duplexes formed by probes with target and non-target DNA was found to be an effective method for determining specific hybridization and was independent of fluctuations in hybridization signal intensity. This form of analysis was more robust than the traditional approach based on hybridization intensity, and enabled the detection of individual Phytophthora species and mixtures thereof.