Oligonucleotides with conjugated dihydropyrroloindole tripeptides: base composition and backbone effects on hybridization.

The ability of conjugated minor groove binding (MGB) residues to stabilize nucleic acid duplexes was investigated by synthesis of oligonucleotides bearing a tethered dihydropyrroloindole tripeptide (CDPI3). Duplexes bearing one or more of these conjugated MGBs were varied by base composition (AT- or GC-rich oligonucleotides), backbone modifications (phosphodiester DNA, 2'-O-methyl phosphodiester RNA or phosphorothioate DNA) and site of attachment of the MGB moiety (5'- or 3'-end of either duplex strand). Melting temperatures of the duplexes were determined. The conjugated CDPI3 residue enhanced the stability of virtually all duplexes studied. The extent of stabilization was backbone and sequence dependent and reached a maximum value of 40-49 degrees C for d(pT)8. d(pA)8. Duplexes with a phosphorothioate DNA backbone responded similarly on CDPI3 conjugation, although they were less stable than analogous phosphodiesters. Modest stabilization was obtained for duplexes with a 2'-O-methyl RNA backbone. The conjugated CDPI3 residue stabilized GC-rich DNA duplexes, albeit to a lesser extent than for AT-rich duplexes of the same length.     

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.     

Parallel thermodynamic analysis of duplexes on oligodeoxyribonucleotide microchips.

A microchip method has been developed for massive and parallel thermodynamic analyses of DNA duplexes. Fluorescently labeled oligonucleotides were hybridized with oligonucleotides immobilized in the 100 x 100 x 20 mum gel pads of the microchips. The equilibrium melting curves for all microchip duplexes were measured in real time in parallel for all microchip duplexes. Thermodynamic data for perfect and mismatched duplexes that were obtained using the microchip method directly correlated with data obtained in solution. Fluorescent labels or longer linkers between the gel and the oligonucleotides appeared to have no significant effect on duplex stability. Extending the immobilized oligonucleotides with a four-base mixture from the 3'-end or one or two universal bases (5-nitroindole) from the 3'- and/or 5'-end increased the stabilities of their duplexes. These extensions were applied to increase the stabilities of the duplexes formed with short oligonucleotides in microchips, to significantly lessen the differences in melting curves of the AT- and GC-rich duplexes, and to improve discrimination of perfect duplexes from those containing poorly recognized terminal mismatches. This study explored a way to increase the efficiency of sequencing by hybridization on oligonucleotide microchips.     

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.     

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.     

Parallel multiplex thermodynamic analysis of coaxial base stacking in DNA duplexes by oligodeoxyribonucleotide microchips

Parallel thermodynamic analysis of the coaxial stacking effect of two bases localized in one strand of DNA duplexes has been performed. Oligonucleotides were immobilized in an array of three-dimensional polyacrylamide gel pads of microchips (MAGIChips‘). The stacking effect was studied for all combinations of two bases and assessed by measuring the increase in melting temperature and in the free energy of duplexes formed by 5mers stacked to microchip-immobilized 10mers. For any given interface, the effect was studied for perfectly paired bases, as well as terminal mismatches, single base overlaps, single and double gaps, and modified terminal bases. Thermodynamic parameters of contiguous stacking determined by using microchips closely correlated with data obtained in solution. The extension of immobilized oligonucleotides with 5,6-dihydroxyuridine, a urea derivative of deoxyribose, or by phosphate, decreased the stacking effect moderately, while extension with FITC or Texas Red virtually eliminated stacking. The extension of the immobilized oligonucleotides with either acridine or 5-nitroindole increased stacking to mispaired bases and in some GC-rich interfaces. The measurements of stacking parameters were performed in different melting buffers. Although melting temperatures of AT- and GC-rich oligonucleotides in 5 M tetramethylammonium chloride were equalized, the energy of stacking interaction was significantly diminished.     

Recognition of native DNA methylation by the PvuII restriction endonuclease

Recognizing the methylation status of specific DNA sequences is central to the function of many systems in eukaryotes and prokaryotes. Restriction–modification systems have to distinguish between ‘self’ and ‘non-self’ DNA and depend on the inability of restriction endonucleases to cleave their DNA substrates when the DNA is appropriately methylated. These endonucleases thus provide a model system for studying the recognition of DNA methylation by proteins. We have characterized the interaction of R·PvuII with DNA containing the physiologically relevant N4-methylcytosine modification. R·PvuII binds ^N4mC-modified DNA and cleaves it very slowly. Methylated strands in hemimethylated duplexes were cleaved at a higher rate than in fully methylated duplexes, in parallel with a higher binding affinity for hemimethylated DNA. The co-crystal structures of R·PvuII–DNA, together with a mutagenesis study, have implicated specific amino acids in recognition of the methylatable base; one of these is His84. We report that replacing His84 with Ala reduced the rate of cleavage of unmodified DNA but, in contrast, slightly increased the cleavage of ^N4mC-modified DNA.     

Spatial distribution analysis of AT- and GC-rich regions in nuclei using corrected fluorescence resonance energy transfer.

We employed microscopic intensity-based fluorescence resonance energy transfer (FRET) images with correction by donor and acceptor concentrations to obtain unbiased maps of spatial distribution of the AT- and GC-rich DNA regions in nuclei. FRET images of 137 bovine aortic endothelial cells stained by the AT-specific donor Hoechst 33258 and the GC-specific acceptor 7-aminoactinomycin D were acquired and corrected for the donor and acceptor concentrations by the Gordon's method based on the three fluorescence filter sets. The corrected FRET images were quantitatively analyzed by texture analysis to correlate the spatial distribution of the AT- and GC-rich DNA regions with different phases of the cell cycle. Both visual observation and quantitative texture analysis revealed an increased number and size of the low FRET efficiency centers for cells in the G(2)/M-phases, compared to the G(1)-phase cells. We have detected cell cycle-dependent changes of the spatial organization and separation of the AT- and GC-rich DNA regions. Using the corrected FRET (cFRET) technique, we were able to detect early DNA separation stages in late interphase nuclei.     

DNA hybridization kinetics: zippering,internal displacement and sequence dependence

Although the thermodynamics of DNA hybridization is generally well established, the kinetics of this classic transition is less well understood. Providing such understanding has new urgency because DNA nanotechnology often depends critically on binding rates. Here, we explore DNA oligomer hybridization kinetics using a coarse-grained model. Strand association proceeds through a complex set of intermediate states, with successful binding events initiated by a few metastable base-pairing interactions, followed by zippering of the remaining bonds. But despite reasonably strong interstrand interactions, initial contacts frequently dissociate because typical configurations in which they form differ from typical states of similar enthalpy in the double-stranded equilibrium ensemble. Initial contacts must be stabilized by two or three base pairs before full zippering is likely, resulting in negative effective activation enthalpies. Non-Arrhenius behavior arises because the number of base pairs required for nucleation increases with temperature. In addition, we observe two alternative pathways—pseudoknot and inchworm internal displacement—through which misaligned duplexes can rearrange to form duplexes. These pathways accelerate hybridization. Our results explain why experimentally observed association rates of GC-rich oligomers are higher than rates of AT- rich equivalents, and more generally demonstrate how association rates can be modulated by sequence choice.     

Oligodeoxynucleotides containing conformationally constrained abasic sites: a UV and fluorescence spectroscopic investigation on duplex stability and structure

The synthesis and incorporation into oligodeoxy­nucleotides of two novel, conformationally restricted abasic (AB) site analogs are described. The stability of oligonucleotide 18mer duplexes containing one such AB site opposite any of the four natural DNA bases was investigated by UV melting curve analysis and compared to that of duplexes containing a conformationally flexible propanediol unit 1 or a tetrahydrofuran unit 2 as an AB site analog. No major differences in the melting temperatures (ΔT_m 0–3°C) between the different abasic duplexes were observed. All AB duplexes were found to have T_ms that were lower by 9–15°C relative to a fully matched 18mer control duplex, and by 4–10°C relative to the corresponding 19mer duplexes in which the AB site is replaced by a mismatched nucleobase. Thus we conclude that the loss of stability of a duplex that is encountered by removal of a nucleobase from the stack cannot be compensated with conformational restriction of the AB site. From the van’t Hoff transition enthalpies obtained from the melting curves, it appears that melting cooperativity is higher for the duplexes containing the conformationally rigid AB sites. Fluorescence quenching experiments with duplexes containing the fluorescent base 2-amino­purine (2AP) opposite the AB sites showed a weak tendency towards more efficient stacking of this base in duplexes containing the conformationally constrained AB sites. Thus, such AB sites may structurally stabilize the cavity formed by the removal of a base. Potential applications emerging from the properties of such conformationally constrained AB sites in DNA diagnostics are discussed.     

The consequences of base-pair substitution mutations in AT- and GC-rich bacteria

The likely consequences, in terms of premature stop codons, detectable missense mutants, silent missense mutants, and degenerate codon changes, have been determined for all 12 individual base substitution changes. This has been done for the full, 61 sense codon, genetic code and also for the much more limited codon availabilities of AT- or GC-rich DNA. The specificities and outcomes of individual base substitutions are likely to be rather different at AT- or GC-rich extremes, and also from the situation at an intermediate DNA base-ratio where all 61 sense codons are available. In particular, at DNA base-ratio extremes many mutations will be to non-utilized codons, which may well act as nonsense mutants. These in turn will give novel classes of suppressor-containing revertants. Even in bacteria with intermediate DNA base-ratios, particular codons for a given amino acid may be favoured, over alternatives, because their use maximizes, or minimizes, the mutational consequences of one, or more, base substitution changes.     

DNA display III. Solid-phase organic synthesis on unprotected DNA

DNA-directed synthesis represents a powerful new tool for molecular discovery. Its ultimate utility, however, hinges upon the diversity of chemical reactions that can be executed in the presence of unprotected DNA. We present a solid-phase reaction format that makes possible the use of standard organic reaction conditions and common reagents to facilitate chemical transformations on unprotected DNA supports. We demonstrate the feasibility of this strategy by comprehensively adapting solid-phase 9-fluorenylmethyoxycarbonyl–based peptide synthesis to be DNA-compatible, and we describe a set of tools for the adaptation of other chemistries. Efficient peptide coupling to DNA was observed for all 33 amino acids tested, and polypeptides as long as 12 amino acids were synthesized on DNA supports. Beyond the direct implications for synthesis of peptide–DNA conjugates, the methods described offer a general strategy for organic synthesis on unprotected DNA. Their employment can facilitate the generation of chemically diverse DNA-encoded molecular populations amenable to in vitro evolution and genetic manipulation.     

Insights into the DNA stabilizing contributions of a bicyclic cytosine analogue: crystal structures of DNA duplexes containing 7,8-dihydropyrido [2,3-d]pyrimidin-2-one

The incorporation of the bicyclic cytosine analogue 7,8-dihydropyrido[2,3-d]pyrimidin-2-one (X) into DNA duplexes results in a significant enhancement of their stability (3–4 K per modification). To establish the effects of X on the local hydrogen-bonding and base stacking interactions and the overall DNA conformation, and to obtain insights into the correlation between the structure and stability of X-containing DNA duplexes, the crystal structures of [d(CGCGAATT-X-GCG)]₂ and [d(CGCGAAT-X-CGCG)]₂ have been determined at 1.9–2.9 Å resolutions. In all of the structures, the analogue X base pairs with the purine bases on the opposite strands through Watson–Crick and/or wobble type hydrogen bonds. The additional ring of the X base is stacked on the thymine bases at the 5′-side and overall exhibits greatly enhanced stacking interactions suggesting that this is a major contribution to duplex stabilization.     

Protein–DNA footprinting by endcapped duplex oligodeoxyribonucleotides

Oligodeoxyribonucleotides (5′-phosphorylated) of varying lengths were capped using a polyamide linker to form thermodynamically stable, endcapped DNA duplexes containing 8–14 bp. We have employed these endcapped DNA duplexes as tools to determine the DNA footprint of T4 DNA ligase. By high-performance liquid chromatography and PAGE analysis of the ligation mixtures of the endcapped DNA duplexes, we have found that by varying the lengths and the position of the nick, we can determine the minimal DNA-binding site as well as the mode of binding (symmetrical or asymmetrical binding) by the enzyme. The results of the study revealed that a 11 bp endcapped duplex was the shortest duplex effectively ligated. Dependence of ligation efficiency on nick position demonstrates that T4 DNA ligase bound asymmetrically to its DNA substrate. The use of a set of thermodynamically stable endcapped deoxyribonucleoside duplexes as a tool to elucidate the DNA footprint provides an efficient strategy for footprinting, which avoids ambiguities associated with chemical and biochemical footprinting methods.     

Structural study of DNA duplexes containing the (6–4) photoproduct by fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET) experiments have been performed to elucidate the structural features of oligonucleotide duplexes containing the pyrimidine(6–4)pyrimidone photoproduct, which is one of the major DNA lesions formed at dipyrimidine sites by UV light. Synthetic 32mer duplexes with and without the (6–4) photoproduct were prepared and fluorescein and tetramethylrhodamine were attached, as a donor and an acceptor, respectively, to the aminohexyl linker at the C5 position of thymine in each strand. Steady-state and time-resolved analyses revealed that both the FRET efficiency and the fluorescence lifetime of the duplex containing the (6–4) photoproduct were almost identical to those of the undamaged duplex, while marked differences were observed for a cisplatin-modified duplex, as a model of kinked DNA. Lifetime measurements of a series of duplexes containing the (6–4) photoproduct, in which the fluorescein position was changed systematically, revealed a small unwinding at the damage site, but did not suggest a kinked structure. These results indicate that formation of the (6–4) photoproduct induces only a small change in the DNA structure, in contrast to the large kink at the (6–4) photoproduct site reported in an NMR study.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Texture analysis of fluorescence lifetime images of AT- and GC-rich regions in nuclei.

We used intensity and fluorescence lifetime microscopy (FLIM) of 3T3 nuclei to investigate the existence of AT-rich and GC-rich regions of the nuclear DNA. Hoechst 33258 (Ho) and 7-aminoactinomycin D (7-AAD) were used as fluorescence probes specific for AT and GC base pairs, respectively. YOYO-1 (Yo) was used as a dye that displays distinct fluorescence lifetimes when bound to AT or GC base pairs. We combined fluorescence imaging of Ho and 7-AAD with time-resolved measurements of Yo and took advantage of an additional information content of the time-resolved fluorescence. Because a single nucleus could not be stained and measured with all three dyes, we used texture analysis to compare the spatial distribution of AT-rich and GC-rich DNA in 100 nuclei in different phases of the cell cycle. The fluorescence intensity-based analysis of Ho- or 7-AAD-stained images indicates increased number and larger size of the DNA condensation centers in the G2/M-phases compared to G0/1-phases. The lifetime-based study of Yo-stained images suggests spatial separation of the AT- or GC-rich DNA regions in the G2/M-phase. Texture analysis of fluorescence intensity and lifetime images was used to quantitatively study the spatial change of condensation and separation of AT- and GC-rich DNA during the cell cycle.     

DNA sequence analysis by hybridization with oligonucleotide microchips: MALDI mass spectrometry identification of 5mers contiguously stacked to microchip oligonucleotides

Matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) has been applied to increase the informational output from DNA sequence analysis. It has been used to analyze DNA by hybridization with microarrays of gel-immobilized oligonucleotides extended with stacked 5mers. In model experiments, a 28 nt long DNA fragment was hybridized with 10 immobilized, overlapping 8mers. Then, in a second round of hybridization DNA–8mer duplexes were hybridized with a mixture of 10 5mers. The stability of the 5mer complex with DNA was increased to raise the melting temperature of the duplex by 10–15°C as a result of stacking interaction with 8mers. Contiguous 13 bp duplexes containing an internal break were formed. MALDI MS identified one or, in some cases, two 5mers contiguously stacked to each DNA–8mer duplex formed on the microchip. Incorporating a mass label into 5mers optimized MALDI MS monitoring. This procedure enabled us to reconstitute the sequence of a model DNA fragment and identify polymorphic nucleotides. The application of MALDI MS identification of contiguously stacked 5mers to increase the length of DNA for sequence analysis is discussed.     

Synergistic effects between analogs of DNA and RNA improve the potency of siRNA-mediated gene silencing

We report that combining a DNA analog (2′F-ANA) with rigid RNA analogs [2′F-RNA and/or locked nucleic acid (LNA)] in siRNA duplexes can produce gene silencing agents with enhanced potency. The favored conformations of these two analogs are different, and combining them in a 1–1 pattern led to reduced affinity, whereas alternating short continuous regions of individual modifications increased affinity relative to an RNA:RNA duplex. Thus, the binding affinity at key regions of the siRNA duplex could be tuned by changing the pattern of incorporation of DNA-like and RNA-like nucleotides. These heavily or fully modified duplexes are active against a range of mRNA targets. Effective patterns of modification were chosen based on screens using two sequences targeting firefly luciferase. We then applied the most effective duplex designs to the knockdown of the eIF4E binding proteins 4E-BP1 and 4E-BP2. We identified modified duplexes with potency comparable to native siRNA. Modified duplexes showed dramatically enhanced stability to serum nucleases, and were characterized by circular dichroism and thermal denaturation studies. Chemical modification significantly reduced the immunostimulatory properties of these siRNAs in human peripheral blood mononuclear cells.     

Computational studies on full-length Ku70 with DNA duplexes: base interactions and a helical path

The Ku70/80 heterodimer is among the first responding proteins to recognize and bind the DNA double strand breaks (DSBs). Once Ku is loaded at the DSB, it works as a scaffold to recruit other repair factors in non-homologous end joining thereby facilitates the following repair processes. In this work, we characterized the detailed interactions and binding free energies between a Ku70 subunit and several DNA duplexes, by using some well-established computational methods. The results reveal that the structure of the protein may suffer certain contractions without the company of Ku80, and may experience large conformational changes in the presence of different DNA duplexes. Notably, we observe the closest interactions between Ku70 and DNA can be easily strengthened to form H-bonds with the bases in the minor groove, which is unexpected. However, this finding is supported by the presence of a similar bond between Ku80 and DNA in the published crystal structure (PDB code 1JEY). We suggest that these interactions are responsible for the observed pausing sites when Ku translocates along DNA and the subtle difference in binding with AT- and GC-rich DNA ends. Additionally, simulations indicate the inner surface of the ring encircling the DNA is not flat, but contains a delicate clamp like structure, which is ideal to grip the two strands of DNA in the minor groove and confine the movement of the duplex in a unique helical path.