The influence of intercalator binding on DNA triplex stability: correlation with effects on A-tract duplex structure

The triplex form of DNA is of interest because of a possible biological role as well as the potential therapeutic use of this structure. In this paper the stabilizing effects of two intercalating drugs, ethidium and the quinoxaline derivative 9-OH-B220, on DNA triplexes have been studied by thermal denaturation measurements. The corresponding duplex structures of the DNA triplex systems investigated are either A-tract or normal B-DNA. The largest increases in the triplex melting temperatures caused by the intercalators were found for sequences having A-tract duplex structures. Inserting a single base pair with an N2-amino group in the minor groove, e.g. a G-C pair, breaks up the A-tract duplex structure and also reduces the stabilizing effect of the drugs on the triplex melting temperatures. The large drug-induced increase in triplex melting temperature for complexes having an original duplex A-tract structure is correlated with a low initial melting point of the triplex, not with the triplex being unusually stable in the presence of the drug. Hence, we conclude that the large thermal stabilizing effect exhibited by ethidium and 9-OH-B220 on dTn.dAn-dTn triplexes is partly caused by the intercalators breaking up the intrinsic A-tract structure of the underlying duplex.     

Stabilization of DNA triple helices by a series of mono- and disubstituted amidoanthraquinones.

We have used quantitative DNase I footprinting to measure the relative affinities of four disubstituted and two monosubstituted amidoanthraquinone compounds for intermolecular DNA triplexes, and have examined how the position of the attached base-functionalized substituents affects their ability to stabilize DNA triplexes. All four isomeric disubstituted derivatives examined stabilize DNA triplexes at micromolar or lower concentrations. Of the compounds studied the 2,7-disubstituted amidoanthraquinone displayed the greatest triplex affinity. The order of triplex affinity for the other disubstituted ligands decreases in the order 2,7 > 1,8 = 1,5 > 2,6, with the equivalent monosubstituted compounds being at least an order of magnitude less efficient. The 1,5-disubstituted derivative also shows some interaction with duplex DNA. These results have been confirmed by molecular modelling studies, which provide a rational basis for the structure-activity relationships. These suggest that, although all of the compounds bind through an intercalative mode, the 2,6, 2,7 and 1,5 disubstituted isomers bind with their two side groups occupying adjacent triplex grooves, in contrast with the 1,8 isomer which is positioned with both side groups in the same triplex groove.     

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-).     

Structural properties of hybrid triplex of polycation deoxyribonucleic S-methylthiourea (DNmt) strands with a complementary DNA strand, probed by nanosecond molecular dynamics

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)8 with a complementary DNA oligomer strand d(Ap)8 have been carried out with explicit water solvent and Na+Cl- 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 x DNA x 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 thermodynamic 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 degrees for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to -30 degrees, 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 x DNA x DNmt hybrid triplex is compared to the DNG x DNA x DNG hybrid triplex (In DNG the -O-(PO2-)-O- linkers in DNA is replaced by -NH-C(=N+H2)-NH-).     

Quantitative analysis of the ion-dependent folding stability of DNA triplexes

A DNA triplex is formed through binding of a third strand to the major groove of a duplex. Due to the high charge density of a DNA triplex, metal ions are critical for its stability. We recently developed the tightly bound ion (TBI) model for ion-nucleic acids interactions. The model accounts for the potential correlation and fluctuations of the ion distribution. We now apply the TBI model to analyze the ion dependence of the thermodynamic stability for DNA triplexes. We focus on two experimentally studied systems: a 24-base DNA triplex and a pair of interacting 14-base triplexes. Our theoretical calculations for the number of bound ions indicate that the TBI model provides improved predictions for the number of bound ions than the classical Poisson-Boltzmann (PB) equation. The improvement is more significant for a triplex, which has a higher charge density than a duplex. This is possibly due to the higher ion concentration around the triplex and hence a stronger ion correlation effect for a triplex. In addition, our analysis for the free energy landscape for a pair of 14-mer triplexes immersed in an ionic solution shows that divalent ions could induce an attractive force between the triplexes. Furthermore, we investigate how the protonated cytosines in the triplexes affect the stability of the triplex helices.     

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

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

DNA intramolecular triplexes containing dT --> dU substitutions: unfolding energetics and ligand binding

We used a combination of optical and calorimetric techniques to investigate the incorporation of deoxythymidine --> deoxyuridine (dT --> dU) substitutions in the duplex and third strand of the parallel intramolecular triplex d(A(7)C(5)T(7)C(5)T(7)) (ATT). UV and differential scanning calorimetry melting experiments show that the incorporation of two substitutions yielded triplexes with lower thermal stability and lower unfolding enthalpies. The enthalpies decrease with an increase in salt concentration, indirectly yielding a heat capacity effect, and the magnitude of this effect was lower for the substituted triplexes. The combined results indicate that the destabilizing effect is due to a decrease in the level of stacking interactions. Furthermore, the minor groove ligand netropsin binds to the minor groove and to the hydrophobic groove, created by the double chain of thymine methyl groups in the major groove of these triplexes. Binding of netropsin to the minor groove yielded thermodynamic profiles similar to that of a DNA duplex with a similar sequence. However, and relative to ATT, binding of netropsin to the hydrophobic groove has a decreased binding affinity and lower binding enthalpy. This shows that the presence of uridine bases disrupts the hydrophobic groove and lowers its cooperativity toward ligand binding. The overall results suggest that the stabilizing effect of methyl groups may arise from the combination of both hydrophobic and electronic effects.     

Structural basis for the unusual properties of 2',5' nucleic acids and their complexes with RNA and DNA

To provide insights into the unusual properties of 2',5' nucleic acids (iso nucleic acids), that includes their rejection by Nature as information molecules, modeling studies have been carried out to examine if they indeed possess the stereochemical ability to form helical duplexes and triplexes, just as their 3',5' linked constitutional isomers. The results show that the formation of helical duplexes with 2',5' linkages demands a mandatory displacement of the Watson and Crick base pairs from the helical axis, as a direct consequence of the lateral shift of the sugar-phosphate backbone from the periphery towards the interior of the helix. Thus, both duplexes and triplexes formed with a 2',5'-sugar-phosphate backbone possess this intrinsic trait, manifested normally only in A type duplexes of DNA and RNA. It was found that only a 10-fold symmetric parallel triplex with isomorphous T.AT triplets is stereochemically favorable for isoDNA with 'extended' nucleotide repeats, unlike the 12-fold symmetric triplex favored by DNA. The wider nature of a 12-fold triplex, concomitant with mandatory slide requirement for helix formation in isoDNA, demands even larger displacement, especially with 'extended' nucleotide structural repeats, thereby violating symmetry. However, a symmetric triplex possessing higher twist, can be naturally formed for isoDNA with a 'compact' nucleotide repeat. Two nanosecond molecular dynamics simulation of a 2',5'-B DNA duplex, formed with an intrinsic base pair displacement of -3.3 A, does not seem to favor a total transition to a typical A type duplex, although enhanced slide, X-displacement, decrease in helical rise and narrowing of the major groove during simulation seem to indicate a trend. Modeling of the interaction between the chimeric isoDNA.RNA duplex and E. coli RNase H has provided a structural basis for the inhibitory action of the enzyme. Interaction of residues Gln 80, Trp 81, Asn 16 and Lys 99, of E. coli RNase H with DNA of the DNA.RNA hybrid, are lost when the DNA backbone is replaced by isoDNA. Based on modeling and experimental observations, it is argued that 2',5' nucleic acids possess restricted conformational flexibility for helical polymorphism. The inability of isoDNA to favor the biologically relevant B form duplex and the associated topological inadequacies related to nucleic acid compaction and interactions with regulatory proteins may be some of the factors that might have led to the rejection of 2',5' links.     

Heteropolymeric triplex-based genomic assay to detect pathogens or single-nucleotide polymorphisms in human genomic samples

Human genomic samples are complex and are considered difficult to assay directly without denaturation or PCR amplification. We report the use of a base-specific heteropolymeric triplex, formed by native duplex genomic target and an oligonucleotide third strand probe, to assay for low copy pathogen genomes present in a sample also containing human genomic duplex DNA, or to assay human genomic duplex DNA for Single Nucleotide Polymorphisms (SNP), without PCR amplification. Wild-type and mutant probes are used to identify triplexes containing FVL G1691A, MTHFR C677T and CFTR mutations. The specific triplex structure forms rapidly at room temperature in solution and may be detected without a separation step. YOYO-1, a fluorescent bis-intercalator, promotes and signals the formation of the specific triplex. Genomic duplexes may be assayed homogeneously with single base pair resolution. The specific triple-stranded structures of the assay may approximate homologous recombination intermediates, which various models suggest may form in either the major or minor groove of the duplex. The bases of the stable duplex target are rendered specifically reactive to the bases of the probe because of the activity of intercalated YOYO-1, which is known to decondense duplex locally 1.3 fold. This may approximate the local decondensation effected by recombination proteins such as RecA in vivo. Our assay, while involving triplex formation, is sui generis, as it is not homopurine sequence-dependent, as are "canonical triplexes". Rather, the base pair-specific heteropolymeric triplex of the assay is conformation-dependent. The highly sensitive diagnostic assay we present allows for the direct detection of base sequence in genomic duplex samples, including those containing human genomic duplex DNA, thereby bypassing the inherent problems and cost associated with conventional PCR based diagnostic assays.     

Enhanced stabilization of the triplexes d(C(+)-T)(6):d(A-G)(6);d(C-T)(6), d(T)(21):d(A)(21);d(T)(21) and poly r(U:AU) by water structure-making solutes

A variety of organic cations, cationic lipids, low molecular weight alcohols, sodium dodecylsulfate, trehalose, glycerol, low molecular weight polyethylene glycols, and DMSO were tested for their ability to modulate the stability of the triplexes d(C(+)-T)(6):d(A-G)(6);d(C-T)(6), d(T)(21):d(A)(21);d(T)(21), poly r(U:A U) and their respective core duplexes, d(A-G)(6);d(C-T)(6), d(A)(21);d(T)(21), poly r(A-U). Very substantial enhancement of triplex stability over that in a physiological salt buffer at pH 7 is obtained with different combinations of triplex and high concentrations of these additives, e.g. trimethylammonium chloride and d(C(+)-T)(6):d(A-G)(6);d(C-T)(6); 2-propanol and d(T)(21):d(A)(21);d(T)(21); ethanol and poly r(U:A;U). Triplex formation is even observed with a 1:1 strand mixture of d(A-G)(6) and d(C-T)(6) in the presence of dimethylammonium, tetramethylammonium, and tetraethylammonium-chloride, as well as methanol, ethanol, and 2-propanol. Triplex stability follows the water structure-making ability (and in some cases the duplex unwinding ability) of the organic cations, the low molecular weight alcohols and other neutral organic compounds, whereas water structure-breaking additives decrease triplex stability. These findings are consistent with those reported in the accompanying paper that triplex formation occurs with a net uptake of water. Since the findings suggest that third strand-binding is facilitated by unwinding of the target duplex, it is inferred that triplex formation may be enhanced by nucleic acid binding proteins operating similarly.     

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.     

Three-stranded DNA structure; is this the secret of DNA homologous recognition?

A novel type of triple-stranded DNA structure was proposed by several groups to play a crucial role in homologous recognition between single- and double-stranded DNA molecules. In this still putative structure a duplex DNA was proposed to co-ordinate a homologous single strand in its major groove side. In contrast to the well-characterized pyrimidine-purine-pyrimidine triplexes in which the two like strands are antiparallel and which are restricted to poly-pyrimidine-containing stretches, the homology-specific triplexes would have like strands in parallel orientation and would not be restricted to any particular sequence provided that there is a homology between interacting DNA molecules. For many years the stereo-chemical possibility of forming homology-dependent three- or four-stranded DNA structures during the pairing stage of recombination reactions was seriously considered in published papers. However, only recently has there been a marked increase in the number of papers that have directly tested the formation of triple-stranded DNA structures during the actual pairing stage of the recombination reaction. Unfortunately the results of these tests are not totally clear cut; while some laboratories presented experimental evidence consistent with the formation of triplexes, others studying the same or very similar systems offered alternative explanations. The aim of this review is to present the current state of the central question in the mechanism of homologous recombination, namely, what kind of DNA structure is responsible for DNA homologous recognition. Is it a novel triplex structure or just a classical duplex?     

Structural diversity of target-specific homopyrimidine peptide nucleic acid-dsDNA complexes

Sequence-selective recognition of double-stranded (ds) DNA by homopyrimidine peptide nucleic acid (PNA) oligomers can occur by major groove triplex binding or by helix invasion via triplex P-loop formation. We have compared the binding of a decamer, a dodecamer and a pentadecamer thymine–cytosine homopyrimidine PNA oligomer to a sequence complementary homopurine target in duplex DNA using gel-shift and chemical probing analyses. We find that all three PNAs form stable triplex invasion complexes, and also conventional triplexes with the dsDNA target. Triplexes form with much faster kinetics than invasion complexes and prevail at lower PNA concentrations and at shorter incubation times. Furthermore, increasing the ionic strength strongly favour triplex formation over invasion as the latter is severely inhibited by cations. Whereas a single triplex invasion complex is formed with the decameric PNA, two structurally different target-specific invasion complexes were characterized for the dodecameric PNA and more than five for the pentadecameric PNA. Finally, it is shown that isolated triplex complexes can be converted to specific invasion complexes without dissociation of the Hoogsteen base-paired triplex PNA. These result demonstrate a clear example of a ‘triplex first’ mechanism for PNA helix invasion.     

Novel Hoogsteen-like bases for configurational recognition of the T-A base pair by DNA triplex formation

Effective sequence-specific recognition of duplex DNA is possible by triplex formation with natural oligonucleotides via Hoogsteen H-bonding. However, triplex formation is in practice limited to pyrimidine oligonucleotides binding duplex A-T or G-C base-pair DNA sequences specifically at homopurine sites in the major groove as T·A-T and C⁺·G-C triplets. Here we report the successful modeling of novel unnatural nucleosides that recognize the T-A DNA base pair by Hoogsteen interaction. Since the DNA triplex can be considered to assume an A-type or B-type conformation, these novel Hoogsteen nucleotides are tested within model A-type and B-type conformation triplex structures. A triplet consisting of the T-A base pair and one of the novel Hoogsteen nucleotides replaces the central T·A-T triplet in the triplex using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration. The entire triplex is energy minimized and the presence of any structural or energetic perturbations due to the central triplet is assessed with respect to the unmodified energy-minimized (T·A-T)₁₁ proposed starting structures. Incorporation of these novel triplets into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets. The plan is to produce a series of Hoogsteen-like bases that preferentially bind the T-A major groove in either an A-type or B-type conformation. Selective recognition of the T-A major groove with respect to the G-C major groove, which presents similar keto and amine placement, is also assessed with configurational preference. Evaluation of the triplex solution structure by using these unnatural bases as binding conformational probes is a prerequisite to the further design of triplet forming bases. © 1996 John Wiley & Sons, Inc.     

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.     

Protonated pyrimidine-purine-purine triplex.

We have studied a protonated pyrimidine-purine-purine (Py-Pu-Pu) triplex, which is formed between the d(C)nd(G)n duplex and the d(AG)m oligonucleotide as the third strand and carries the CG*A+ protonated base-triads. We have observed such an intermolecular complex between a plasmid carrying the d(C)18 d(G)18 insert and the d(AG)5 oligonucleotide without bivalent cations in 200 mM of Na+ at pH4.0. Bivalent cations additionally stabilize the complex. We propose the structures for nearly isomorphous base-triads TA*A, CG*G and CG*A+. To identify the H-DNA-like structure, which includes the triplex between d(C)n d(G)n duplex and the AG-strand, we have cloned in a superhelical plasmid the insert: G10TTAA(AG)5. The data on photofootprinting and chemical modification with diethyl pyrocarbonate, potassium permanganate and dimethyl sulfate demonstrate that the H-like structure with triplex carrying CG*G and CG*A+ base triads is actually formed under acid conditions. In the course of this study we have come across unexpected results on probing of Py-Pu-Pu triplexes by dimethyl sulfate (DMS): the protection effect is observed not only for guanines entering the duplex but also for guanines in the third strand lying in the major groove. We have demonstrated this effect not only for the case the novel protonated Py-Pu-Pu triplex but also for the traditional non-protonated Py-Pu-Pu intramolecular triplex (H*-DNA) formed by the d(C)37 d(G)37 insert in supercoiled plasmid in the presence of Mg2+ ions.     

Molecular dynamics simulation study of DNA triplex formed by mixed sequences in solution

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC).d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usualcut-of of the long range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other ligands.     

Molecular dynamics simulation of the DNA triplex d(TC)5.d(GA)5.d(C+T)5.

A molecular dynamics simulation of the DNA triple helix d(TC)5.d(GA)5.d(C+T)5 is described (C+ represents a protonated cytosine residue). The simulation has been performed using the program AMBER 3.1 and includes counterions and explicit solvent under periodic boundary conditions. Both the dynamic and time-averaged behaviour of the system has been analysed. Considerable deviations from the fibre-diffraction model for DNA triple helix structure are observed, including the repuckering of the purine strand sugars that has been identified in some nuclear magnetic resonance (n.m.r.) studies. The simulation suggests that this conformational change may be driven by the possibility of improved interactions between the phosphate groups of this strand and both the solvent and counterions. Several examples of a particular conformational transition are observed, involving correlated changes in the backbone angles alpha and gamma. These transitions provide a possible explanation for some unusual n.m.r. data that have been reported. The structure of the triple helix major groove also suggests an explanation for the observed stabilization of DNA triplexes by polyvalent cations, and their ability to interact with drugs that bind in the minor groove of DNA duplexes.