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

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.     

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.     

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

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.     

Proton NMR and Optical Spectroscopic Studies on the DNA Triplex Formed by d-A-(G-A)7-G and d-C-(T-C)7-T

Abstract

Triplex and duplex formation of two deoxyribohexadecamers d-A-(G-A)₇-G (a) and d-C-(T-C)₇-T (b) have been studied by UV, CD, fluorescence, and proton NMR spectroscopy. Optical studies of a and b at dilute concentrations (μM range) yielded results similar to those seen for polymers of the same sequence, indicating that these hexadecamers have properties similar to the polymers in regard to triplex formation. The CD spectra of concentrated NMR samples (mM range) are similar to those observed at optical concentrations at both low and high pH, making possible a correlation between CD and NMR studies. In NMR spectra, two imido NH-N hydrogen bonded resonance envelopes at 12.6 and 13.7 ppm indicate that only the duplex conformation is present at pH > 7.7. Four new NH-N hydrogen-bonded resonance envelopes at 12.7, 13.5, 14.2. and 14.9 ppm are observed under acidic conditions (pH 5.6) and the two original NH-N resonances gradually disappear as the pH is lowered. Assignment of these four peaks to Watson-Crick G · C, Hoogsteen T · A. Watson-Crick A · T. and Hoogsteen C⁺ · G hydrogen-bonded imidos, respectively, confirm the formation of triple-stranded DNA NMR results also show that triplex is more stable than duplex at the same salt condition and that triplex melts to single strands directly without going through a duplex intermediate. However, in the melting studies, a structural change within the triple-stranded complex is evident at temperatures significantly below the major helix-to-coil transition. These studies demonstrate the feasibility of using NMR spectroscopy and oligonucleotide model compounds a and b for the study of DNA triplex formation.     

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.     

Effect of 5-Methylcytosine on the Structure and Stability of DNA. Formation of Triple-Stranded Concatenamers by Overlapping Oligonucleotides

Abstract

A triple helix can be formed upon binding of a pyrimidine oligonucleotide to the major groove of a homopurine-homopyrimidine (R·Y) double-stranded DNA target site. Here, we report that this reaction can be influenced by base methylation. The pyrimidine strand 5′- TmCTmCTmCTmCTTmCT (mY₁₂), whose cytosine residues are methylated at C5, does not bind the duplex 5′-AGAGAGAGAAGA·3′-TCTCTCTCTTCT (R₁₂·Y₁₂) to yield a 12-triad triplex, as would be expected from these DNA sequences. Rather, a complex of overlapping oligonucleotides, which we define concatenamer, is formed. The concatenamer is clearly evidenced by Polyacrylamide gel electrophoresis (PAGE) since it migrates with a smeared band of very low mobility. The stoichiometry of the concatenamer, determined by both UV mixing curves and electrophoresis, is surprisingly found to be (R₁₂· 2mY₁₂)_n, thus showing that the unmethylated Y₁₂ strand is excluded from the complex. Denaturation experiments performed by ultraviolet absorbance (UV) and differential scanning calorimetry (DSC) show that the concatenamers melt with a single and highly cooperative transition whose Tm strongly depends on pH. Overall, the data point to the conclusion that the concatenamers are in triple helix, where the methylated mY₁₂ strand is engaged in both Watson-Crick and Hoogsteen base pairings, thus displacing the Y₁₂ strand from the R₁₂·Y₁₂ duplex. A possible mechanism of concatenamer formation is proposed. The results presented in this paper show that 5-methylcytosine brings about a strong stabilizing effect on both double and triple DNA helices, and that pyrimidine oligonucleotides containing 5-methylcytosine can displace from R·Y duplexes the analogous non-methylated strand. The advantage of using methylated oligonucleotides in antisense technology is discussed.     

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.     

Structure of the DNA Duplex d(ATTAAT)2 with Hoogsteen Hydrogen Bonds

The traditional Watson-Crick base pairs in DNA may occasionally adopt a Hoogsteen conformation, with a different organization of hydrogen bonds. Previous crystal structures have shown that the Hoogsteen conformation is favored in alternating AT sequences of DNA. Here we present new data for a different sequence, d(ATTAAT)₂, which is also found in the Hoogsteen conformation. Thus we demonstrate that other all-AT sequences of DNA with a different sequence may be found in the Hoogsteen conformation. We conclude that any all-AT sequence might acquire this conformation under appropriate conditions. We also compare the detailed features of DNA in either the Hoogsteen or Watson-Crick conformations.     

Interaction of 9-O-(ω-amino) alkyl ether berberine analogs with poly(dT)·poly(dA)*poly(dT) triplex and poly(dA)·poly(dT) duplex: a comparative study

Isoquinoline alkaloids and their analogs represent an important class of molecules for their broad range of clinical and pharmacological utility. These compounds are of current interest owing to their low toxicity and excellent chemo preventive properties. These alkaloids can play important role in stabilising the nucleic acid triple helices. The present study has focused on the interaction of five 9-O-(ω-amino) alkyl ether berberine analogs with the DNA triplex poly(dT)·poly(dA)*poly(dT) and the parent duplex poly(dA)·poly(dT) studied using various biophysical techniques. Scatchard analysis of the spectral data indicated that the analogs bind both to the duplex and triplex in a non-cooperative manner in contrast to the cooperative binding of berberine to the DNA triplex. Strong intercalative binding to the DNA triplex structure was revealed from ferrocyanide quenching, fluorescence polarization and viscosity results. Thermal melting studies demonstrated higher stabilization of the Hoogsteen base paired third strand of the DNA triplex compared to the Watson–Crick strand. Circular dichroism studies suggested a stronger perturbation of the DNA triplex conformation by the alkaloid analogs compared to the duplex. The binding was entropy-driven in each case and the entropy contribution to free energy increased as the length of the alkyl side chain increased. The analogs exhibited stronger binding affinity to the triple helical structure compared to the parent double helical structure.     

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.     

Novel Hoogsteen-like Bases for Recognition of the C-G 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 that bind 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 modelling of novel unnatural nucleosides that recognize the C-G DNA base pair by Hoogsteen-like major groove interaction. These novel Hoogsteen nucleotides are examined within model A-type and B-type conformation triplex structures since the DNA triplex can be considered to incorporate A-type and/or B-type configurational properties. Using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration, a triplet comprised of a C-G base pair and the novel Hoogsteen nucleotide, Y2, replaces the central T·A-T triplet in the triplex. The presence of any structural or energetic perturbations due to the central triplet in the energy-minimized triplex is assessed with respect to the unmodified energy minimized (T·A-T)₁₁ starting structures. Incorporation of this novel triplet into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets.     

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     

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.     

Strong Binding of Single-stranded DNA by Stem-loop Oligonucleotides

Abstract

We report that oligodeoxynucleotides which form stem-loop hairpin structures and which have pyrimidine-rich loops can form strong complexes with complementary single-stranded DNA sequences. Stem-loop oligonucleotides were constructed with a 25-nt T-rich loop and with variable Watson-Crick stems. The complexes of these oligomers with the sequence dA_gwere studied by thermal denaturation. Evidence is presented that the complexes are one-to-one, bimolecular complexes in which the pyrimidine loop bases comprise the outer strands in a pyr · pur · pyr triplex, in effect chelating the purine strand in the center of the loop. Melting temperatures for the loop complexes are shown to be up to 29 °C higher than Watson- Crick duplex of the same length. It is shown that the presence of a stem increases stability of the triplex relative to an analogous oligomer without a stem. The effect of stem length on the stability of such a complex is examined. Such hairpin oligomers represent a new approach to the sequence-specific binding of single-stranded RNA and DNA. In addition, the finding raises the possibility that such a complex may exist in natural RNA folded sequences.     

Evidence for (PNA)2/DNA triplex structure upon binding of PNA to dsDNA by strand displacement

The binding of PNA (peptide nucleic acid) T₂CT₂CT₄-LysNH₂ to the double-stranded DNA target 5′ -A₂GA₂GA₄ was studied by KMnO₄ and dimethylsulfate (DMS) probing. It is found that upon sequence-specific strand displacement binding of the PNA to the dsDNA target concomitant protection of the N-7 of guanines within the target takes place. It is furthermore shown that the binding of this PNA is more efficient at pH 5.5 that at pH 6.5 and very inefficient at pH 7.5. These results clearly indicate that C⁺G Hoogsteen base pairing is present and important for binding and that the strand displacement complex therefore involves a PNA·DNA-PNA triplex.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

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 usual “cut-off” 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 legends.