Structure of the alpha-homo-DNA:RNA duplex and the function of twist and slide to catalogue nucleic acid duplexes

High-resolution NMR studies of an alpha-homo-DNA:RNA duplex reveal the formation of a right-handed parallel-oriented helix. It differs significantly from a standard A- or B-type helix by a small twist value (26.2 degrees ), which leads to a helical pitch of 13.7 base pairs per helical turn, a negative inclination (-1.78 Angstrom) and a large x displacement (5.90 Angstrom). The rise (3.4 Angstrom) is similar to that found in B-DNA. The solution of this new helix structure has stimulated us to develop a mathematical and geometrical model based on slide and twist parameters to describe nucleic acid duplexes. All existing duplexes can be positioned within this landscape, which can be used to understand the helicalization process.     

Structural studies of LNA:RNA duplexes by NMR: conformations and implications for RNase H activity

We have used NMR and CD spectroscopy to study the conformations of modified oligonucleotides (locked nucleic acid, LNA) containing a conformationally restricted nucleotide (T(L)) with a 2'-O,4'-C-methylene bridge. We have investigated two LNA:RNA duplexes, d(CTGAT(L)ATGC):r(GCAUAUCAG) and d(CT(L)GAT(L)AT(L)GC):r(GCAUAUCAG), along with the unmodified DNA:RNA reference duplex. Increases in the melting temperatures of +9.6 degrees C and +8.1 degrees C per modification relative to the unmodified duplex were observed for these two LNA:RNA sequences. The three duplexes all adopt right-handed helix conformations and form normal Watson-Crick base pairs with all the bases in the anti conformation. Sugar conformations were determined from measurements of scalar coupling constants in the sugar rings and distance information derived from 1H-1H NOE measurements; all the sugars in the RNA strands of the three duplexes adopt an N-type conformation (A-type structure), whereas the sugars in the DNA strands change from an equilibrium between S- and N-type conformations in the unmodified duplex towards more of the N-type conformation when modified nucleotides are introduced. The presence of three modified T(L) nucleotides induces drastic conformational shifts of the remaining unmodified nucleotides of the DNA strand, changing all the sugar conformations except those of the terminal sugars to the N type. The CD spectra of the three duplexes confirm the structural changes described above. On the basis of the results reported herein, we suggest that the observed conformational changes can be used to tune LNA:RNA duplexes into substrates for RNase H: Partly modified LNA:RNA duplexes may adopt a duplex structure between the standard A and B types, thereby making the RNA strand amenable to RNase H-mediated degradation.     

Bicyclo-DNA: a Hoogsteen-selective pairing system

Background: The natural nucleic acids (DNA and RNA) can adopt a variety of structures besides the antiparallel double helix described by Watson and Crick, depending on base sequence and solvent conditions. Specifically base-paired DNA structures with regular backbone units include left-handed and parallel duplexes and triple and quadruple helical arrangements. Given the base-pairing pattern of the natural bases, preferences for how single strands associate are determined by the structure and flexibility of the sugar-phosphate backbone. We set out to determine the role of the backbone in complex formation by designing DNA analogs with well defined modifications in backbone structure.Results: We recently developed a DNA analog (bicyclo-DNA) in which one (γ) of the six torsion angles (a-ζ) describing the DNA-backbone conformation is fixed in an orientation that deviates from that observed in B-DNA duplexes by about +100°, a shift from the synclinal to the antiperiplanar range. Upon duplex formation between homopurine and homopyrimidine sequences, this analog preferentially selects the Hoogsteen and reversed Hoogsteen mode, forming A-T and G-C⁺ base pairs. Base-pair formation is highly selective, but degeneracy is observed with respect to strand orientation in the duplex.Conclusions: The flexibility and orientation of the DNA backbone can influence the preferences of the natural bases for base-pairing modes, and can alter the relative stability of duplexes and triplexes.     

2-Aminopurine: a probe of structural dynamics and charge transfer in DNA and DNA:RNA hybrids

Spectroscopic techniques are employed to probe relationships between structural dynamics and charge transfer (CT) efficiency in DNA duplexes and DNA:RNA hybrids containing photoexcited 2-aminopurine (Ap). To better understand the variety of interactions and reactions, including CT, between Ap and DNA, the fluorescence behavior of Ap is investigated in a full series of redox-inactive as well as redox-active assemblies. Thus, Ap is developed as a dual reporter of structural dynamics and base-base CT reactions in nucleic acid duplexes. CD, NMR, and thermal denaturation profiles are consistent with the family of DNA duplexes adopting a distinct conformation versus the DNA:RNA hybrids. Fluorescence measurements establish that the d(A)-r(U) tract of the DNA:RNA hybrid exhibits enhanced structural flexibility relative to that of the d(A)-d(T) tract of the DNA duplexes. The yield of CT from either G or 7-deazaguanine (Z) to Ap in the assemblies was determined by comparing Ap emission in redox-active G- or Z-containing duplexes to otherwise identical duplexes in which the G or Z is replaced by inosine (I), the redox-inactive nucleoside analogue. Investigations of CT not only demonstrate efficient intrastrand base-base CT in the DNA:RNA hybrids but also reveal a distance dependence of CT yield that is more shallow through the d(A)-r(U) bridge of the A-form DNA:RNA hybrids than through the d(A)-d(T) bridge of the B-form DNA duplexes. The shallow distance dependence of intrastrand CT in DNA:RNA hybrids correlates with the increased conformational flexibility of bases within the hybrid duplexes. Measurements of interstrand base-base CT provide another means to distinguish between the A- and B-form helices. Significantly, in the A-form DNA:RNA hybrids, a similar distance dependence is obtained for inter- and intrastrand reactions, while, in B-DNA, a more shallow distance dependence is evident with interstrand CT reactions. These observations are consistent with evaluations of intra- and interstrand base overlap in A- versus B-form duplexes. Overall, these data underscore the sensitivity of CT chemistry to nucleic acid structure and structural dynamics.     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ～10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).     

Consecutive GC base pairs determine the energy barrier of DNA duplex formation under molecularly crowded conditions

The kinetics of DNA duplex formation was affected by the addition of PEGs with different masses (MW = 200-8000) to an aqueous solution; for each condition, two duplexes (5'-TAGGTTATAA-3'/5'-TTATAACCTA-3' and 5'-CAGGTCACAG-3'/5'-CTGTGACCTG-3') with different stabilities were formed after overcoming the same association activation energy barrier, suggesting that the formation of consecutive GC base pairs in the helices rather than the helix terminus is the initiation nucleus for DNA duplex formation not only in the absence, but also in the presence of PEGs.     

Chemical ligation as a method for the assembly of double-stranded nucleic acids: Modifications and local structure studies

A new procedure was developed as an alternative to the enzymatic assembly of natural and modified double-stranded DNAs using chemical reagent (chemical ligation). BrCN was suggested as an efficient coupling reagent, which induces superfast reactions in DNA duplexes. The physicochemical properties and the structure of new types of DNA duplexes, which are the substrates for chemical ligation, with breaks in phosphodiester chains, including concatemers, were studied. Chemical ligation was applied to prepare biologically active 17–200 base-pair double-stranded DNAs and DNA-RNA block-copolymers, to incorporate various modifications into DNA duplexes including pyrophosphate and phosphoramidate unnatural internucleotide bonds. The unique possibilities of this approach were demonstrated in the development of methods for circularization of oligodeoxy ribonucleotides and assembly of branched DNAs. The structural-kinetic concept of chemicalligation was created and the relationship between the reactivity of interacting groups and sequence-dependent local conformation of the ligation site in B-DNA was established. The lesser efficiency of chemical ligation of RNA fragments in comparison to that of DNA analogs was demonstrated and rationalized. This approach was used as a sensitive monitor of a stable double helix formation and third-strand binding to a DNA duplex.Translated from Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1889–1911, August, 1996.     

Stabilizing or destabilizing oligodeoxynucleotide duplexes containing single 2'-deoxyuridine residues with 5-alkynyl substituents

The 5-position of pyrimidines in DNA duplexes offers a site for introducing alkynyl substituents that protrude into the major groove and thus do not sterically interfere with helix formation. Substituents introduced at the 5-position of the deoxyuridine residue of dU:dA base pairs may stabilize duplexes and reinforce helices weakened by a low G/C content, which would otherwise lead to false negative results in DNA chip experiments. Here we report on a method for preparing oligonucleotides with a 5-alkynyl substituent at a 2'-deoxyuridine residue by on-support Sonogashira coupling involving the fully assembled oligonucleotide. A total of 25 oligonucleotides with 5-alkynyl substituents were prepared. The substituents either decrease the UV melting point of the duplex with the complementary strand or increase it by up to 7.1 degrees C, compared with that of the unmodified control duplex. The most duplex-stabilizing substituent, a pyrenylbutyramidopropyne moiety, is likely to intercalate but does not prevent sequence-specific base pairing of the modified deoxyuridine residue or the neighboring nucleotides. It also increases the signal for a target strand when employed on a small oligonucleotide microarray. The ability to tune the melting point of a DNA dodecamer duplex with a single side chain over a temperature range of >11 degrees C may prove useful when developing DNA sequences for biomedical applications.     

Solution structure of a DNA duplex containing a biphenyl pair

Hydrogen-bonding and stacking interactions between nucleobases are considered to be the major noncovalent interactions that stabilize the DNA and RNA double helices. In recent work we found that one or multiple biphenyl pairs, devoid of any potential for hydrogen bond formation, can be introduced into a DNA double helix without loss of duplex stability. We hypothesized that interstrand stacking interactions of the biphenyl residues maintain duplex stability. Here we present an NMR structure of the decamer duplex d(GTGACXGCAG) d(CTGCYGTCAC) that contains one such X/Y biaryl pair. X represents a 3',5'-dinitrobiphenyl- and Y a 3',4'-dimethoxybiphenyl C-nucleoside unit. The experimentally determined solution structure shows a B-DNA duplex with a slight kink at the site of modification. The biphenyl groups are intercalated side by side as a pair between the natural base pairs and are stacked head to tail in van der Waals contact with each other. The first phenyl rings of the biphenyl units each show tight intrastrand stacking to their natural base neighbors on the 3'-side, thus strongly favoring one of two possible interstrand intercalation structures. In order to accommodate the biphenyl units in the duplex the helical pitch is widened while the helical twist at the site of modification is reduced. Interestingly, the biphenyl rings are not static in the duplex but are in dynamic motion even at 294 K.     

Synthesis of [Ru(phen)(2)dppz](2+)-tethered oligo-DNA and studies on the metallointercalation mode into the DNA duplex

To explore the binding properties of [Ru(phen)(2)dppz](2+) complex (phen = 1,10-phenanthroline, dppz = dipyrido[3,2-a:2',3'-c]phenazine) in a sequence-specific manner in DNA duplex, it was tethered through the dppz ligand to a central position as well as both at the 3'- and 5'-ends of oligodeoxyribonucleotide (ODN). The middle [Ru(phen)(2)dppz](2+)-ODN tethered was resolved and isolated as four pure diastereomers, while the 3'- or 5'-[Ru(phen)(2)dppz](2+)-ODNs were inseparable on RP-HPLC. Thermal stability of the (Ru(2+)-ODN).DNA duplexes is found to increase considerably (DeltaT(m) = 12.8-23.4 degrees C), depending upon the site of the covalent attachment of the tethered [Ru(phen)(2)dppz](2+) complex, or the chirality of the [Ru(phen)(2)dppz](2+)-linker tethered at the middle of the ODN, compared to the unlabeled counterpart. Gross differences in CD between the [Ru(phen)(2)dppz](2+)-tethered and the native DNA duplexes showed that the global duplex conformation of the former has considerably altered from the B-type, but is still recognized by DNase I. The thermal melting studies, CD measurements, as well as DNase I digestion data, are interpreted as a result of intercalation of the dppz moiety, which is realized by threading of the Ru(phen)(2) complex part through the DNA duplex core. DNase I footprinting with four diastereomerically pure middle ([Ru(phen)(2)dppz](2+)-ODN).DNA duplexes furthermore showed that the tethered [Ru(phen)(2)dppz](2+)-linker chirality dictates the stereochemical accessibility of various phosphodiester moieties (around the intercalation site) toward the cleavage reaction by the enzyme. The diastereomerically pure ruthenium-modified duplexes, with the well-defined pi-stack, will be useful to explore stereochemistry-dependent energy- and electron-transfer chemistry to understand oxidative damage to the DNA double helix as well as the long-range energy- and electron-transfer processes with DNA as a reactant.     

Stabilization by extra-helical thymines of a DNA duplex with Hoogsteen base pairs

We present the crystal structure of the DNA duplex formed by d(ATATATCT). The crystals contain seven stacked antiparallel duplexes in the asymmetric unit with A.T Hoogsteen base pairs. The terminal CT sequences bend over so that the thymines enter the minor groove and form a hydrogen bond with thymine 2 of the complementary strand in the Hoogsteen duplex. Cytosines occupy extra-helical positions; they contribute to the crystal lattice through various kinds of interactions, including a unique CAA triplet. The presence of thymine in the minor groove apparently contributes to the stability of the DNA duplex in the Hoogsteen conformation. These observations open the way toward finding under what conditions the Hoogsteen duplex may be stabilized in vivo. The present crystal structure also confirms the tendency of A.T-rich oligonucleotides to crystallize as long helical stacks of duplexes.     

DNA double helices recognize mutual sequence homology in a protein free environment

The structure and biological function of the DNA double helix are based on interactions recognizing sequence complementarity between two single strands of DNA. A single DNA strand can also recognize the double helix sequence by binding in its groove and forming a triplex. We now find that sequence recognition occurs between intact DNA duplexes without any single-stranded elements as well. We have imaged a mixture of two fluorescently tagged, double helical DNA molecules that have identical nucleotide composition and length (50% GC; 294 base pairs) but different sequences. In electrolytic solution at minor osmotic stress, these DNAs form discrete liquid-crystalline aggregates (spherulites). We have observed spontaneous segregation of the two kinds of DNA within each spherulite, which reveals that nucleotide sequence recognition occurs between double helices separated by water in the absence of proteins, consistent with our earlier theoretical hypothesis. We thus report experimental evidence and discuss possible mechanisms for the recognition of homologous DNAs from a distance.     

Synthesis and properties of nucleic acid analogues consisting of a benzene-phosphate backbone

[Chemical reaction: See text] The synthesis and properties of a nucleic acid analogue consisting of a benzene-phosphate backbone are described. The building blocks of the nucleic acid analogue are composed of bis(hydroxymethyl)benzene residues connected to nucleobases via the biaryl-like axis. Stabilities of the duplexes were studied by thermal denaturation. It was found that the thermal stabilities of the duplexes composed of the benzene-phosphate backbone are highly dependent on their sequences. The duplexes with the benzene-phosphate backbone comprised of the mixed sequences were thermally less stable than the natural DNA duplexes, whereas that composed of the homopyrimidine and homopurine sequences was thermally and thermodynamically more stable than the corresponding natural DNA duplex. It was suggested that the analogues more efficiently stabilize the duplexes in a B-form duplex rather than in an A-form duplex. Thus, the duplexes consisting of the benzene-phosphate backbone, especially composed of the homopyrimidine and homopurine sequences, may offer a novel structural motif useful for developing novel materials applicable in the fields of bio- and nanotechnologies.     

Synthesis and biophysical characterization of oligonucleotides containing a 4'-selenonucleotide

The first synthesis of oligonucleotides containing 4'-selenium-modified ribonucleotides (4'-Se-rN) is described. Four sequences containing 4'-Se-rT were successfully synthesized and compared with DNA and RNA oligonucleotides containing a dT, rT, or LNA insert in place of the 4'-Se-rT. The 4'-Se-rT behaved more like rT than dT in its effects on binding affinity, despite the DNA-like structure previously observed for the nucleoside, suggesting that a conformational switch occurs upon incorporation into an oligonucleotide. Incorporation of 4'-Se-rT into A-RNA and hybrid duplexes led to increased binding affinity, while incorporation into B-DNA destabilized the duplex to the same extent as an rT nucleotide.     

NMR and UV studies of 3'-S-phosphorothiolate modified DNA in a DNA : RNA hybrid dodecamer duplex; implications for antisense drug design

High-resolution NMR spectroscopy has been used to establish the conformational consequences of the introduction of a single 3[prime or minute]-S-phosphorothiolate link in the DNA strand of a DNA : RNA hybrid. These systems are of interest as potential antisense therapeutic agents. Previous studies on similarly modified dinucleotides have shown that the conformation of the sugar to which the sulfur is attached shifts to the north (C(3[prime or minute])-endo/C(2[prime or minute])-exo). Comparisons made between NOESY cross-peak intensities, and coupling constants from PE-COSY spectra, for both non-modified and modified duplexes confirm that this conformational shift is also present in the double helical oligonucleotide system. In addition it is noted that in both the dinucleotides and the modified duplex, the conformation of the sugar ring 3[prime or minute] to the site of modification is also shifted to the north. That this pattern is observed in the small monomeric system as well as the larger double helix is suggestive of some pre-ordering of the sequences. The conclusion is supported by consideration of the (1)H chemical shifts of the heterocyclic bases near the site of the modification. The enhanced stability that these conformational changes should bring was confirmed by UV thermal melting studies. Subsequently a series of singly and doubly 3[prime or minute]-S-phosphorothiolate-modified duplexes were investigated by UV. The results are indicative of an additive effect of the modification with thermodynamic benefit being derived from alternate spacing of two modified linkers.     

The structure of LNA:DNA hybrids from molecular dynamics simulations: the effect of locked nucleotides

Locked nucleic acids (LNAs) exhibit a modified sugar fragment that is restrained to the C3'-endo conformation. LNA-containing duplexes are rather stable and have a more rigid structure than DNA duplexes, with a propensity for A-conformation of the double helix. To gain detailed insight into the local structure of LNA-modified DNA oligomers (as a foundation for subsequent exploration of the electron-transfer capabilities of such modified duplexes), we carried out molecular dynamics simulations on a set of LNA:DNA 9-mer duplexes and analyzed the resulting structures in terms of base step parameters and the conformations of the sugar residues. The perturbation introduced by a single locked nucleotide was found to be fairly localized, extending mostly to the first neighboring base pairs; such duplexes featured a B-type helix. With increasing degree of LNA modification the structure gradually changed; the duplex with one complete LNA strand assumed a typical A-DNA structure. The relative populations of the sugar conformations agreed very well with NMR data, lending credibility to the validity of the computational protocol.     

DNA sequence context and multiplex hybridization reactions: melting studies of heteromorphic duplex DNA complexes

Heteromorphic hybrid duplex DNA complexes are duplex states, other than perfectly matched duplexes, that can form when single strands comprising several different perfectly matched duplexes are simultaneously present in solution. Such cross-hybridization "side reactions" are of particular nuisance in multiplex reaction schemes, where many strands are designed to hybridize in parallel fashion with only their corresponding perfect complement strand. Relative to the perfect match duplexes, the sequence dependent features of these heteromorphic duplex states and their thermodynamic stability are an important consideration for multiplex hybridization reaction design. We have measured absorbance versus temperature melting curves and performed differential scanning calorimetry measurements on various mixtures of eight different 24 base single strands. When perfect complementary pairs of strands are mixed in single reactions, four perfectly matched duplexes form. When mixtures of strands that are not perfectly matched are prepared and analyzed, melting transitions for cross-hybridization are observed along with significant hyperchromicity changes. This is indicative of a melting hybrid, heteromorphic duplex states formed from two nonperfectly matched strands. In addition, when both the perfectly matched and noncomplementary strands are mixed together (in multiplex hybridization reactions) at molar ratios of 1:1, 3:1, and 1:3, evidence of perfect duplex and heteromorphic duplex complexes is found in all cases. A new analytical tool for considering heterogeneous, duplex complexes in multiplex hybridization mixtures is presented and employed to interpret the acquired melting data.     

How do size-expanded DNA nucleobases enhance duplex stability? Computational analysis of the hydrogen-bonding and stacking ability of xDNA bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.     

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.