Structural Properties of Hybrid Triplex of Polycation Deoxyribonucleic S-methylthiourea (DNmt) Strands with a Complementary DNA Strand,Probed by Nanosecond Molecular Dynamics

Abstract

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)₈ with a complementary DNA oligomer strand d(Ap)₈ have been carried out with explicit water solvent and Na counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt·DNA·DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermody- namic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20° for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to ?30°, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3′-endo domain and C2′-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt·DNA·DNmt hybrid triplex is compared to the DNG·DNA·DNG hybrid triplex (In DNG the -O-(PO₂-)-O- linkers in DNA is replaced by -NH-C(=N₂)-NH-).     

Base-pair dynamics in an antiparallel DNA triplex measured by catalyzed imino proton exchange monitored via 1H NMR spectroscopy

Using (1)H NMR spectroscopy, the base-pair opening dynamics of an antiparallel foldback DNA triplex and the corresponding duplex has been characterized via catalyzed imino proton exchange. The triplex system was found to be in an equilibrium between a duplex and a triplex form. The exchange rate between the two forms (i.e., the on/off-rate of the third strand) was measured to be 5 s(-1) at 1 degrees C, and the base-pair dynamics of both forms were investigated separately. Both Watson-Crick and reverse Hoogsteen base pairs were found to have base-pair lifetimes in the order of milliseconds. The stability of the Watson-Crick base pairs was, however, substantially increased in the presence of the third strand. In the DNA triplex, the opening dynamics of the reverse Hoogsteen base pairs was significantly faster than the dynamics of the Watson-Crick pairs. We were able to conclude that, for both Watson-Crick and reverse Hoogsteen base pairs, spontaneous and individual opening from within the closed base triplet is the dominating opening pathway.     

Proton nuclear magnetic resonance assignments and structural characterization of an intramolecular DNA triplex.

Two-dimensional 1H n.m.r. spectroscopy has been used to study the 31-base DNA oligonucleotide 5'-dAGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3', which folds to form a stable intramolecular triplex in solution at acidic pH. This structure is considerably more difficult to assign than short B-DNA duplexes and requires new assignment methods. The assignment strategy and assignments of almost all of the exchangeable and nonexchangeable resonances are presented. Seven base triplets and one Watson-Crick base-pair form the core of the structure and are connected by a four C and four T loop at either end. The second pyrimidine "strand" (bases 24 to 31) in this intramolecular pyrimidine-purine-pyrimidine triplex binds via Hoogsteen base-pairs in the major groove and is parallel to the purine "strand" (bases 1 to 8). Analysis of the sugar puckers reveals that, contrary to widely accepted belief, the triplex sugars are not predominantly in the N-type (close to C3'-endo) conformation. Except for some of the C nucleotides, all sugars are predominantly S-type (close to C2'-endo). Thus, the duplex DNA does not assume N-type sugar conformations to accommodate a third strand in the major groove. A preliminary model of the triplex structure is presented.     

Spectroscopic studies on the interaction of aristololactam-beta-D-glucoside with DNA and RNA double and triple helices: A comparative study.

The interaction of aristololactam-beta-D-glucoside (ADG), a DNA intercalating alkaloid, with the DNA triplexes, poly(dT). poly(dA)xpoly(dT) and poly(dC).poly(dG)xpoly(dC+), and the RNA triplex poly(rU).poly(rA)xpoly(rU) was investigated by circular dichroic, UV melting profile, spectrophotometric, and spectrofluorimetric techniques. Comparative interaction with the corresponding Watson-Crick duplexes has also been examined under identical experimental conditions. Triplex formation has been confirmed from biphasic thermal melting profiles and analysis of temperature-dependent circular dichroic measurements. The binding of ADG to triplexes and duplexes is characterized by the typical hypochromic and bathochromic effects in the absorption spectrum, quenching of steady-state fluorescence intensity, a decrease in fluorescence quantum yield, an increase or decrease of thermal melting temperatures, and perturbation in the circular dichroic spectrum. Scatchard analysis indicates that ADG binds both to the triplexes and the duplexes in a noncooperative manner. Binding parameters obtained from spectrophotometric measurements are best fit by the neighbor exclusion model. The binding affinity of ADG to the DNA triplexes is substantially stronger than to the RNA triplex. Thermal melting study further indicates that ADG stabilizes the Hoogsteen base-paired third strand of the DNA triplexes whereas it destabilizes the same strand of RNA triplex but stabilizes its Watson-Crick strands. Comparative data reveal that ADG exhibits a stronger binding to the triple helical structures than to the respective double helical structures.     

Proton NMR studies of 5'-d-(TC)(3) (CT)(3) (AG)(3)-3'--a paperclip triplex: the structural relevance of turns

In this study, we present the results of structural analysis of an 18-mer DNA 5'-T(1)C(2)T(3)C(4)T(5)C(6)C(7)T(8)C(9)T(10)C(11)T(12)A(13)G(14)A(15)G(16)A(17)G(18)-3' by proton nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The NMR data are consistent with characteristics for triple helical structures of DNA: downfield shifting of resonance signals, typical for the H3(+) resonances of Hoogsteen-paired cytosines; pH dependence of these H3(+) resonance; and observed nuclear Overhauser effects consistent with Hoogsteen and Watson-Crick basepairing. A three-dimensional model for the triplex is developed based on data obtained from two-dimensional NMR studies and molecular modeling. We find that this DNA forms an intramolecular "paperclip" pyrimidine-purine-pyrimidine triple helix. The central triads resemble typical Hoogsteen and Watson-Crick basepairing. The triads at each end region can be viewed as hairpin turns stabilized by a third base. One of these turns is comprised of a hairpin turn in the Watson-Crick basepairing portion of the 18-mer with the third base coming from the Hoogsteen pairing strand. The other turn is comprised of two bases from the continuous pyrimidine portion of the 18-mer, stabilized by a hydrogen-bond from a purine. This "triad" has well defined structure as indicated by the number of nuclear Overhauser effects and is shown to play a critical role in stabilizing triplex formation of the internal triads.     

NMR studies of triple-strand formation from the homopurine-homopyrimidine deoxyribonucleotides d(GA)4 and d(TC)4

The complexes formed by the homopurine and homopyrimidine deoxyribonucleotides d(GA)4 and d(TC)4 have been investigated by one- and two-dimensional 1H NMR. Under appropriate conditions [low pH, excess d(TC)4 strand] the oligonucleotides form a triplex containing one d(GA)4 and two d(TC)4 strands. The homopurine and one of the homopyrimidine strands are Watson-Crick base paired, and the second homopyrimidine strand is Hoogsteen base paired in the major groove to the d(GA)4 strand. Hoogsteen base pairing in GC base pairs requires hemiprotonation of C; we report direct observation of the C+ imino proton in these base pairs. Both homopyrimidine strands have C3'-endo sugar conformations, but the purine strand does not. The major triplex formed appears to have four TAT and three CGC+ triplets formed by binding of the second d(TC)4 strand parallel to the d(GA)4 strand with a 3' dangling end. In addition to the triplexes formed, at least one other heterocomplex is observed under some conditions.     

Single-Stranded DNA and RNA Targeted Triplex-Formation: UV,CD and Molecular Modeling Studies of Foldback Triplexes Containing Different RNA, 2′-OMe-RNA and DNA Strand Combinations

Abstract

We studied the influence of different 2′-OMe-RNA and DNA strand combinations on single strand targeted foldback triplex formation in the Py.Pu:Py motif using ultraviolet (UV) and circular dichroism (CD) spectroscopy, and molecular modeling. The study of eight combinations of triplexes (D D:D, R* D:D, D D:R*, R* D:R*, D R:D, R* R:D, DR:R*, and R*-R:R*; where the first, middle, and last letters stand for the Hoogsteen Pyrimidine, Watson-Crick [WC] purine and WC pyrimidine strands, respectively, and D, R and R* stand for DNA, RNA and 2′-OMe-RNA strands, respectively) indicate more stable foldback triplex formation with a DNA purine strand than with an RNA purine strand. Of the four possible WC duplexes with RNA/DNA combinations, the duplex with a DNA purine strand and a 2′-O-Me-RNA pyrimidine strand forms the most thermally stable triplex, although its thermal stability is the lowest of all four duplexes. Irrespective of the duplex combination, a 2′-OMe-RNA Hoogsteen pyrimidine strand forms a stable foldback triplex over a DNA Hoogsteen pyrimidine strand confirming the earlier reports with conventional and circular triplexes. The CD studies suggest a B-type conformation for an all DNA homo-foldback triplex (D.D.D), while hetero-foldback triplex spectra suggest intermediate conformation to both Atype and B-type structures. A novel molecular modeling study has been carried out to understand the stereochemical feasibility of all the combinations of foldback triplexes using a geometric approach. The new approach allows use of different combinations of chain geometries depending on the nature of the chain (RNA vs. DNA).     

Evidence from CD spectra and melting temperatures for stable Hoogsteen-paired oligomer duplexes derived from DNA and hybrid triplexes.

The pyr*pur.pyr type of nucleic acid triplex has a purine strand that is Hoogsteen-paired with a parallel pyrimidine strand (pyr*pur pair) and that is Watson-Crick-paired with an antiparallel pyrimidine strand (pur.pyr pair). In most cases, the Watson-Crick pair is more stable than the Hoogsteen pair, although stable formation of DNA Hoogsteen-paired duplexes has been reported. Using oligomer triplexes of repeating d(AG)12 and d(CT)12 or r(CU)12 sequences that were 24 nt long, we found that hybrid RNA*DNA as well as DNA*DNA Hoogsteen-paired strands of triplexes can be more stable than the Watson-Crick-paired strands at low pH. The structures and relative stabilities of these duplexes and triplexes were evaluated by circular dichroism (CD) spectroscopy and UV absorption melting studies of triplexes as a function of pH. The CD contributions of Hoogsteen-paired RNA*DNA and DNA*DNA duplexes were found to dominate the CD spectra of the corresponding pyr*pur.pyr triplexes.     

The paperclip triplex: understanding the role of apex residues in tight turns

In this study, we investigate the role of the apex nucleotides of the two turns found in the intramolecular "paperclip" type triplex DNA formed by 5'-TCTCTCCTCTCTAGAGAG-3'. Our previously published structure calculations show that residues C7-A18 form a hairpin turn via Watson-Crick basepairing and residues T1-C6 bind into the major groove of the hairpin via Hoogsteen basepairing resulting in a broad turn of the T1-T12 5'-pyrimidine section of the DNA. We find that only the C6C7/G18 apex triad (and not the T12A13/T1 apex triad) is required for intramolecular triplex formation, is base independent, and occurs whether the purine section is located at the 5' or 3' end of the sequence. NMR spectroscopy and molecular dynamics simulations are used to investigate a bimolecular complex (which retains only the C6C7/G18 apex) in which a pyrimidine strand 5'- TCTCTCCTCTCT-3' makes a broad fold stabilized by the purine strand 5'-AGAGAG-3' via Watson Crick pairing to the T8-T12 and Hoogsteen basepairing to T1-T5 of the pyrimidine strand. Interestingly, this investigation shows that this 5'-AGAGAG-3' oligo acts as a new kind of triplex forming oligonucleotide, and adds to the growing number of triplex forming oligonucleotides that may prove useful as therapeutic agents.     

Effect of a 5'-phosphate on the stability of triple helix.

An effect of 5'-phosphorylation on the stability of triple helical DNA containing pyrimidine:purine:pyrimidine strands has been demonstrated by both gel electrophoresis and UV melting. A 5'-phosphate on the purine-rich middle strand of a triple helix lowers the stability of triple helix formation by approximately 1 kcal/mol at 25 degrees C. The middle strand is involved in both Watson-Crick and Hoogsteen base pairing. In contrast, a 5'-phosphate on the pyrimidine-rich strands, which are involved in either Watson-Crick or Hoogsteen base pairing, has a smaller effect on the stability of triple helix. The order of stability is: no phosphate on either strand > phosphate on both pyrimidine strands > phosphate on purine strand > phosphate on all three strands. Differential stability of triple helix species is postulated to stem from an increase in rigidity due to steric hindrance from the 5'-phosphate. This result indicates that labelling with 32P affect equilibrium in triplex formation.     

pH and cation effects on the properties of parallel pyrimidine motif DNA triplexes

The effects of cytosine protonation and various cations on the properties of parallel pyrimidine motif DNA triplexes were intensively investigated and characterized by several different techniques, such as circular dichroism (CD) conformation, ultraviolet (UV) melting, differential scanning calorimetry (DSC) thermal denaturation, and surface plasmon resonance (SPR) real-time dynamics. The comparative CD spectra of the triplex and the corresponding homoduplexes showed that the negative peak at approximately 218 nm would be the eigenpeak of the Hoogsteen paired strand, and moreover, the formation pathway of a triplex was significantly pH-dependent and fell into three groups: under acidic conditions, the triplex is formed by a one-step docking, under near physiological conditions, the Watson-Crick duplex is first structured and then accepts the Hoogsteen third strand into its major groove, and under basic conditions, the triplex is not formed. The pH-dependent thermodynamics of the global triplex, the Watson-Crick antiparallel duplex, and the Crick-Hoogsteen parallel duplex were comparatively discussed for the first time. These data revealed that the thermodynamic stabilities of the Watson-Crick-Hoogsteen triplex and the Crick-Hoogsteen duplex would be strongly dependent on cytosine protonation, but a low-pH environment somewhat destabilized the Watson-Crick duplex. The binding energy of triplex formation would be different from the unfolding energy of triplex melting under acidic conditions due to the disparity in the pathway between the formation and unfolding of a triplex. Real-time dynamic measurements showed that the association and dissociation rate constants of a duplex-to-triplex formation are (1.98 +/- 0.24) x 10(3) M(-1) s(-1) and (4.09 +/- 0.96) x 10(-4) s(-1) at 20 degrees C and pH 6.0, respectively. The formation energy of the duplex-to-triplex transition derived from SPR measurements was in agreement with the unfolding energy of the free Hoogsteen paired duplex derived from UV measurements. The calorimetric enthalpies of the triplex-to-duplex-to-single transition were 39.3 and 75.3 kcal/mol under near physiological conditions (pH 7.0), respectively, which were underestimated relative to the van't Hoff enthalpies. In addition, the effects of various cations, ionic strength, mixed-valent cations, and the position of the C(+)xG.C triplets on the thermodynamics of the triplexes were addressed under near physiological conditions. The interaction of metal ions with the triplexes clearly depended on the type and ionic strength of the cations, and the efficiency with which the cations stabilized the global triplex was in the order Mg(2+) > Mn(2+) > Ca(2+) > Ba(2+) > Na(+). These observations would be useful for the design of triplex-forming oligonucleotides for antigene drugs and therapeutic purposes.     

Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

Influence of 5-bromodeoxycytosine substitution on triplex DNA stability and conformation

Three triple-helical hairpin DNAs with substitution of 5-bromocytosine for cytosine in different strands have been investigated by molecular mechanics and Raman spectroscopy. The stability of the three substituted triplexes were compared with the corresponding unsubstituted triplex DNA by the molecular mechanics method. Base stacking interactions and strand--strand interactions of each triplex were analyzed in detail. Sugar conformations in these triplexes have been determined by both vibrational spectroscopy and molecular dynamics simulation. The hairpin triplexes with substitution occurring in strand I or both in strands I and III have the main sugar conformation of C3'-endo, while the triplex with substitution occurring in strand III is the combination of C3'-endo and C2'-endo sugar conformation. Theoretical results are basically in agreement with experiments.     

Homopurine and homopyrimidine strands complementary in parallel orientation form an antiparallel duplex at neutral pH with A-C,G-T,and T-C mismatched base pairs

DNA sequences d-TGAGGAAAGAAGGT (a 14-mer) and d-CTCCTTTCTTCC (a 12-mer) are complementary in parallel orientation forming either Donahue (reverse Watson-Crick) base pairing at neutral pH or Hoogsteen base pairing at slightly acidic pH. The structure of the complex formed by dissolving the two strands in equimolar ratio in water has been investigated by nmr. At neutral pH, the system forms an ordered antiparallel duplex with five A : T and four G : C Watson-Crick base pairs and three mismatches, namely G-T, A-C, and T-C. The nuclear Overhauser effect cross-peak pattern suggests an overall B-DNA conformation with major structural perturbations near the mismatches. The duplex has a low melting point and dissociates directly into single strands with a broad melting profile. The hydrogen-bonding schemes in the mismatched base pairs have been investigated. It has been shown earlier that in acidic pH, the system prefers a triple-stranded structure with two pyrimidine strands and one purine strand. One of the pyrimidine strands has protonated cytosines, forms Hoogsteen base pairing, and is aligned parallel to the purine strand; the other has nonprotonated cytosines and has base-pairing scheme similar to the one discussed in this paper. The parallel duplex is therefore less stable than either the antiparallel duplex or the triplex, in spite of its perfect complementarity. © 1997 John Wiley & Sons, Inc. Biopoly 41: 773–784, 1997     

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations.

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.     

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.     

Structural characterization of separated H DNA conformers

Polypyrimidine/polypurine DNA sequences in plasmids can adopt protonated triplex-containing structures (H DNA) in response to negative superhelical stress and low pH. A d(TC)17-d(GA)17 insert adopts two isomeric protonated structures, which differ in degree of helical unwinding. The variant forms of individual topoisomers were separated by agarose gel electrophoresis and their reactivities to permanganate and acid-induced depurination were compared. Depurination patterns of the individual conformers indicate that in the more mobile form (H-y5) the 5'-half of the d(GA)n strand participates in a triplex while in the other (H-y3) the 3'-half forms the triplex. The H-y5 form is more stable than the H-y3 form at low negative superhelix densities. Because of the difference in helical unwinding, the H-y5 form becomes relatively less stable as the superhelix density increases. Topological models of the two forms show that providing there is no linkage at the tips of the triple helical segments one more positive twist is localized in the H-y5 form than in the H-y3 form. The foldback in the pyrimidine strand of the H-y5 form is less accessible to solvent than that of the H-y3 form as assessed by its lower reactivity to permanganate. Consideration of a pyrimidine loop model (Harvey, S. C., Luo, J., & Lavery, R. (1988) Nucleic Acids Res. 16, 11795-11809) suggests that the unique stability of the H-y5 form results from Watson-Crick base pairs between residues of the d(TC)n loop and the d(GA)n strand as it exits the triplex.     

Comparison of positively charged DNG with DNA duplexes: a computational approach

Molecular dynamics is used to investigate the structural properties of the cationic DNA analogue deoxynucleic guanidine (DNG), in which a guanidinium group replaces the phosphate moiety of DNA. This study examines the DNG duplex dodecamers d(Ag)(12).d(Tg)(12) and d(Gg)(12).d(Cg)(12), as well as their DNA counterparts. Watson-Crick base-pairing is maintained in the solvated DNG duplex models during the 5ns simulations. The idealized DNG dodecamers assume many parameters characteristic of the corresponding native DNA, assuming B-DNA conformations. Several helical parameters are rather unique to DNG, including buckle, slide, inclination, propeller, and X-displacement. Fewer transitions in backbone torsions occur in the DNG duplexes compared to those of the DNA, which may result from the greater rigidity of the sp(2) hybridized guanidinium group verses the flexible sp(3) phosphate group. The DNG helices have exceptionally shallow major grooves and very deep minor grooves. The major and minor groove widths of DNG are narrower than those of the respective DNA counterparts.     

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.