Hydration sites in purine.purine.pyrimidine and pyrimidine.purine.pyrimidine DNA triplexes in aqueous solution

Systematic study of hybrid triplex topology and stability suggests a general triplex-mediated regulatory mechanism

By combining in silico, biophysical, and in vitro experiments, we decipher the topology, physical, and potential biological properties of hybrid-parallel nucleic acids triplexes, an elusive structure at the basis of life. We found that hybrid triplex topology follows a stability order: r(Py)-d(Pu)·r(Py) > r(Py)-d(Pu)·d(Py) > d(Py)-d(Pu)·d(Py) > d(Py)-d(Pu)·r(Py). The r(Py)-d(Pu)·d(Py) triplex is expected to be preferred in the cell as it avoids the need to open the duplex reducing the torsional stress required for triplex formation in the r(Py)-d(Pu)·r(Py) topology. Upon a massive collection of melting data, we have created the first predictor for hybrid triplex stability. Leveraging this predictor, we conducted a comprehensive scan to assess the likelihood of the human genome and transcriptome to engage in triplex formation. Our findings unveil a remarkable inclination-of both the human genome and transcriptome-to generate hybrid triplex formation, particularly within untranslated (UTRs) and regulatory regions, thereby corroborating the existence of a triplex-mediated regulatory mechanism. Furthermore, we found a correlation between nucleosome linkers and Triplex-forming sequence (TFS) which agree with a putative role of triplexes in arranging chromatin structure.

New insights into DNA triplexes: residual twist and radial difference as measures of base triplet non-isomorphism and their implication to sequence-dependent non-uniform DNA triplex

DNA triplexes are formed by both isomorphic (structurally alike) and non-isomorphic (structurally dissimilar) base triplets. It is espoused here that (i) the base triplet non-isomorphism may be articulated in structural terms by a residual twist (Delta(t) degrees), the angle formed by line joining the C1'...C1' atoms of the adjacent Hoogsteen or reverse Hoogsteen (RH) base pairs and the difference in base triplet radius (Delta(r) A), and (ii) their influence on DNA triplex is largely mechanistic, leading to the prediction of a high (t + Delta(t))degrees and low (t - Deltat)degrees twist at the successive steps of Hoogsteen or RH duplex of a parallel or antiparallel triplex. Efficacy of this concept is corroborated by molecular dynamics (MD) simulation of an antiparallel DNA triplex comprising alternating non-isomorphic G*GC and T*AT triplets. Conformational changes necessitated by base triplet non-isomorphism are found to induce an alternating (i) high anti and anti glycosyl and (ii) BII and an unusual BIII conformation resulting in a zigzag backbone for the RH strand. Thus, base triplet non-isomorphism causes DNA triplexes into exhibiting sequence-dependent non-uniform conformation. Such structural variations may be relevant in deciphering the specificity of interaction with DNA triplex binding proteins. Seemingly then, residual twist (Delta(t) degrees) and radial difference (Deltar A) suffice as indices to define and monitor the effect of base triplet non-isomorphism in nucleic acid triplexes.

Antiparallel triple helices. Structural characteristics and stabilization by 8-amino derivatives

The structural, dynamical, and recognition properties of antiparallel DNA triplexes formed by the antiparallel d(G#G.C), d(A#A.T), and d(T#A.T) motifs (the pound sign and dot mean reverse-Hoogsteen and Watson-Crick hydrogen bonds, respectively) are studied by means of "state of the art" molecular dynamics simulations. Once the characteristics of the helix are defined, molecular dynamics and thermodynamic integration calculations are used to determine the expected stabilization of the antiparallel triplex caused by the introduction of 8-aminopurines. Finally, oligonucleotides containing 8-aminopurine derivatives are synthesized and tested experimentally using several approaches in a variety of systems. A very large stabilization of the triplex is found experimentally, as predicted by simulations. These results open the possibility for the use of oligonucleotides carrying 8-aminopurines to bind single-stranded nucleic acids by formation of antiparallel triplexes.

Triplex hydration: nanosecond molecular dynamics simulation of the solvated triplex formed by mixed sequences

A theoretical model for the hydration pattern and motion of ions around the triple helical DNA with mixed sequences d(GACTGGTGAC)d(GTCACCAGTC)*d(GACTGGTGAC) in solution, during MD simulation, using the particle mesh Ewald sum method, is elaborated here. The AMBER 5.0 force field has been used during the simulation in solvent. The simulation studies support a dynamically stable atmosphere around the DNA triplex in solution over the entire length of the trajectory. The results have been compared with Hoogsteen triplexes and examined in the context of the observed behaviour of hydration in crystallographic data of duplexes. The dynamical organization of counterions and water molecules around the triplex formed by mixed sequences is described here. It has been observed that cations prefer to bind between two adjoining purines of the second and the third strands. The idea of localized complexes (mobile counterions in unspecific electronegative pockets around the DNA triplex with water molecules) may have important implications for understanding the specificity of the interactions of nucleic acids with proteins and other ligands.

Investigation of Overhauser effects between pseudouridine and water protons in RNA helices

The inherent chemical properties of RNA molecules are expanded by posttranscriptional modification of specific nucleotides. Pseudouridine (psi), the most abundant of the modified bases, features an additional imino group, NH1, as compared with uridine. When psi forms a Watson-Crick base pair with adenine in an RNA helix, NH1 is positioned within the major groove. The presence of psi often increases thermal stability of the helix or loop in which it is found [Hall, K. B. & McLaughlin, L. (1992) Nucleic Acids Res. 20, 1883-1889]. X-ray crystal structures of transfer RNAs [e.g., Arnez, J. & Steitz, T. (1994) Biochemistry 33, 7560-7567] have depicted water molecules bridging psiNH1 groups and nearby phosphate oxygen atoms, but direct evidence for this interaction in solution has not been acquired. Toward this end, we have used a rotating-frame Overhauser effect spectroscopy-type NMR pulse sequence with a CLEAN chemical-exchange spectroscopy spin-lock pulse train [Hwang, T.-L., Mori, S., Shaka, A. J. & van Zijl, P. C. M. (1997) J. Am. Chem. Soc. 119, 6203-6204] to test for psiNH1-water cross-relaxation effects within two RNA helices: (i) a complementary duplex, in which psi is not associated with structural change, and (ii) an RNA duplex representing the eukaryotic pre-mRNA branch-site helix from Saccharomyces cerevisiae, in which a conserved psi extrudes the branch-site adenosine from the helix. Our data implicate a water-psiNH1 hydrogen bond both in stabilizing the complementary helix and in favoring formation of the unique structure of the branch-site helix.

Internal dynamics in a DNA triple helix probed by (1)H-(15)N-NMR spectroscopy

The amino group of adenine plays a key role in maintaining DNA triple helical structures, being the only functional group in DNA that is involved in both Watson-Crick and Hoogsteen hydrogen bonds. In the present work we have probed the internal dynamics of the adenine amino group in the intramolecular YRY triple helix formed by the 31-mer DNA oligonucleotide d(AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT). The DNA triple helix was specifically labeled with (15)N at the amino group of the adenine in the fifth position. The rotation rate of the labeled amino group was measured as a function of temperature using (1)H-(15)N heteronuclear NMR spectroscopy. The results indicate that, in the DNA triple helix, the rotation of the adenine amino group is greatly slowed relative to that in a DNA double helix. The temperature dependence of the rotation rate suggests a large entropic contribution to this effect, which may originate from different hydration patterns of the adenine amino group in the two structures.

Hoogsteen-based parallel-stranded duplexes of DNA. Effect of 8-amino-purine derivatives

The structure of parallel-stranded duplexes of DNA-containing a mixture of guanines (G) and adenines (A) is studied by means of molecular dynamics (MD) simulation, as well as NMR and circular dichroism (CD) spectroscopy. Results demonstrate that the structure is based on the Hoogsteen motif rather than on the reverse Watson-Crick one. Molecular dynamics coupled to thermodynamic integration (MD/TI) calculations and melting experiments allowed us to determine the effect of 8-amino derivatives of A and G and of 8-amino-2'-deoxyinosine on the stability of parallel-stranded duplexes. The large stabilization of the parallel-stranded helix upon 8-amino substitution agrees with a Hoogsteen pairing, confirming MD, NMR, and CD data, and suggests new methods to obtain DNA triplexes for antigene and antisense purposes.

Temperature dependence of the Raman spectrum of DNA. II. Raman signatures of premelting and melting transitions of poly(dA).poly(dT) and comparison with poly(dA-dT).poly(dA-dT)

The temperature dependence of the Raman spectrum of poly(dA).poly(dT) (dA: deoxyadenosine; dT: thymidine), a model for DNA containing consecutive adenine.thymine (A.T) pairs, has been analyzed using a spectrometer of high spectral precision and sensitivity. Three temperature intervals are distinguished: (a) premelting (10 < t < 70 degrees C), in which the native double helix is structurally altered but not dissociated into single strands; (b) melting (70 < t < 80 degrees C), in which the duplex is dissociated into single strands; and (c) postmelting (80 < t degrees C), in which no significant structural change can be detected. The distinctive Raman difference signatures observed between 10 and 70 degrees C and between 70 and 80 degrees C are interpreted in terms of the structural changes specific to premelting and melting transitions, respectively. Premelting alters the low-temperature conformation of the deoxyribose-phosphate backbone and eliminates base hydrogen bonding that is distinct from canonical Watson-Crick hydrogen bonding; these premelting perturbations occur without disruption of base stacking. Conversely, melting eliminates canonical Watson-Crick pairing and base stacking. The results are compared with those reported previously on poly(dA-dT).poly(dA-dT), the DNA structure consisting of alternating A.T and T.A pairs (L. Movileanu, J. M. Benevides, and G. J. Thomas, Jr. Journal of Raman Spectroscopy, 1999, Vol. 30, pp. 637-649). Poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) exhibit strikingly dissimilar temperature-dependent Raman profiles prior to the onset of melting. However, the two duplexes exhibit very similar melting transitions, including the same Raman indicators of ruptured Watson-Crick pairing, base unstacking and collapse of backbone order. A detailed analysis of the data provides a comprehensive Raman assignment scheme for adenosine and thymidine residues of B-DNA, delineates Raman markers diagnostic of consecutive A.T and alternating A.T/T.A tracts of DNA, and identifies the distinct Raman difference signatures for premelting and melting transitions in the two types of sequences.

The effect of amino groups on the stability of DNA duplexes and triplexes based on purines derived from inosine

The effect of amino groups attached at positions 2 and 8 of the hypoxanthine moiety in the structure, reactivity and stability of DNA duplexes and triplexes is studied by means of quantum mechanical calculations, as well as extended molecular dynamics (MD) and thermodynamic integration (MD/TI) simulations. Theoretical estimates of the change in stability related to 2'-deoxyguanosine (G) --> 2'-deoxyinosine (I) --> 8-amino-2'-deoxyinosine (8AI) mutations have been experimentally verified, after synthesis of the corresponding compounds. An amino group placed at position 2 stabilizes the duplex, as expected, and surprisingly also the triplex. The presence of an amino group at position 8 of the hypoxanthine moiety stabilizes the triplex but, surprisingly, destabilizes the duplex. The subtle electronic redistribution occurring upon the introduction of an amino group on the purine seems to be responsible for this surprising behavior. Interesting 'universal base' properties are found for 8AI.

Stabilisation of TG- and AG-containing antiparallel DNA triplexes by triplex-binding ligands

We have used DNase I footprinting to examine the interaction of several triplex-binding ligands with antiparallel TG- and AG-containing triplexes. We find that although a 17mer TG-containing oligonucleotide on its own fails to produce a footprint at concentrations as high as 30 microM, this interaction can be stabilised by several ligands. Within a series of disubstituted amidoanthraquinones we find that the 2,7- regioisomer affords the best stabilisation of this TG triplex, though the 1,8- isomer also stabilises this interaction to some extent. By contrast the 1,5- and 2,6- regioisomers show no interaction with TG triplexes. Similar studies with a 13mer AG-containing oligonucleotide show the opposite pattern of stabilisation: the 2,6- and 1,5- isomers stabilise this triplex, but the 2,7- and 1,8-compounds do not. The polycyclic compound BePI strongly stabilises TG- but not AG-containing triplexes, while a substituted naphthylquinoline interacts with both antiparallel triplex motifs.

Unexpected formation of parallel duplex in GAA and TTC trinucleotide repeats of Friedreich's ataxia

The onset and progress of Friedreich's ataxia (FRDA) is associated with the genetic instability of the (GAA).(TTC) trinucleotide repeats located within the frataxin gene. The instability of these repeats may involve the formation of an alternative DNA structure. Poly-purine (R)/poly-pyrimidine (Y) sequences typically form triplex DNA structures which may contribute to genetic instability. Conventional wisdom suggested that triplex structures formed by these poly-purine (R)/poly-pyrimidine (Y) sequences may contribute to their genetic instability. Here, we report the characterization of the single-stranded GAA and TTC sequences and their mixtures using NMR, UV-melting, and gel electrophoresis, as well as chemical and enzymatic probing methods. We show that the FRDA GAA/TTC, repeats are capable of forming various alternative structures. The most intriguing is the observation of a parallel (GAA).(TTC) duplex in equilibrium with the antiparallel Watson-Crick (GAA).(TTC) duplex. We also show that the GAA strands form self-assembled structures, whereas the TTC strands are essentially unstructured. Finally, we demonstrate that the FRDA repeats form only the YRY triplex (but not the RRY triplex) at neutral pH and the complete formation of the YRY triplex requires the ratio of GAA to TTC strand larger than 1:2. The structural features presented here and in other studies distinguish the FRDA (GAA)¿(TTC) repeats from the fragile X (CGG).CCG), myotonic dystrophy (CTG).(CAG) and the Huntington (CAG).(CTG) repeats.

The effect of structural water molecules on the normal mode spectrum of dTn . dAn x dTn DNA

In this paper we present a theoretical treatment of triplex B type DNA hydration using normal mode calculation techniques. Discrete solvent is added as spines of hydration in the Watson-Crick and Crick-Hoogsteen grooves as well as water bridges between the Phosphate groups. The effect of binding the discrete structural waters on the normal mode of vibration of the system was studied by introducing a parameter, Xw, that is proportional to the degree of water binding and inversely proportional to the relative humidity (RH) of the system. We examined the variation of the dipole moments of characteristic modes with Xw. The results show that there is a direct relationship between the degree of binding of the water molecules to the atoms in the triple helix, the relative humidity of the system and the conformation and stability of the triple helix. At high RH and Xw = 0:0 the triple helix has mostly B type conformation characteristics, with C'2 -endo sugars. The emergence of normal modes of vibration characteristic to the A type conformation (C'3 - endo sugars) at Xw = 0:4 and 60% RH indicates a conformational shift towards A-type for some of the sugars between Xw = 0.2 (80% RH) and Xw = 0.4 (60% RH). These results are in agreement with the "economy of hydration hypothesis" of Saenger (Saenger et al., 1986) which maintains that the main difference in the hydration of A- and B- forms of DNA is the presence of water bridges between adjacent Phosphate groups in the low-hydration A-form but not in the B- form. Free energy calculations for the triplex DNA with structural waters show that there is a minimum of the free energy at Xw = 0.2 and the free energy increases with Xw and becomes larger than the free energy of the B conformation without structural waters for Xw equal to and larger than 0.4. This result indicates that the B conformation is more stable with bound structural water molecules (for degrees of water binding that are not over 20% higher than the degree of binding between bulk water molecules). The structural water molecules are bound much tighter in the A conformation than in the B conformation. The model predicts that the B to A transition occurs at higher relative humidities in D2O than in H2O. Part of these results (Dadarlat, 1997) have been subsequently confirmed by the experimental work and MD simulations of Ouali (Ouali et al., 1997). The experimental results showed that the N-type sugars corresponding to the A conformation are clearly detected below 75% RH.

Spectroscopic investigation of an intramolecular DNA triplex containing both G.G:C and T.A:T triads and its complex with netropsin

The triple helix formation by the oligonucleotide 5'd(G4T4G4-[T4]-G4A4G4-[T4]-C4T4C4) ([T4] represents a stretch of 4 thymine residues) has been investigated by UV absorption spectroscopy and circular dichroism. In a 10 mM sodium cacodylate, 0.2 mM disodium EDTA (pH 7) buffer, we show the following significant results: i) In the absence of MgCl2, the oligonucleotide adopts a hairpin duplex structure with the dangling tail 5'd(G4T4G4-[T4]). This 5' extremity, which contains separated runs of four guanine residues, does not assume the expected tetraplex conformation observed when this sequence is free. ii) In the presence of MgCl2, the oligonucleotide folds back on itself twice to give a triple helix via a double hairpin formation, with [T4] single-strand loops. iii) The addition of high concentration of KCl to the preformed triplex does not disrupt the structure. Nevertheless, if the oligonucleotide is allowed to fold back in the presence of K+, triplex formation is inhibited. Circular dichroism studies demonstrate that the oligonucleotide adopts a dimeric conformation, resulting from the association of two hairpin duplexes, via the formation of an antiparallel G-quadruplex by the telomeric 5'd(G4T4G4-[T4]) extremities. iv) Under the experimental conditions used in this report, the triplex melts in a monophasic manner. v) Netropsin, a DNA minor groove ligand, binds to the central site A4/T4 of the duplex and to that of the triplex in an equimolar stoichiometry. In contrast with previous studies concerning pyr.pur:pyr triplexes, thermal denaturation experiments demonstrate that the netropsin binding stabilizes the intramolecular triplex.

Different effects of nonintercalative antitumor drugs on DNA triple helix stability: SN-18071 promotes triple helix formation

The interaction of the nonintercalating bisquaternary ammonium heterocyclic drugs SN-18071 and SN-6999 with a DNA triple helix has been studied using thermal denaturation and CD spectroscopy. Our data show, that both minor groove binders can bind to the triple helix of poly(dA).2poly(dT) under comparable ionic conditions, but they influence the stability of the triplex relative to the duplex structure of poly(dA).poly(dT) in a different manner. SN-18071, a ligand devoid of forming hydrogen bonds, can promote triplex formation and thermally stabilizes it up to 500 mM Na+ concentration. SN-6999 destabilizes the triplex to duplex equibilirium whereas it stabilizes the duplex. The binding constant of SN-18071 is found to be greater than that to the duplex. The stabilizing effect of SN-18071 is explained by electrostatic interactions of three ligand molecules with the three grooves of the triple stranded structure. From the experiments it is concluded that SN-6999 binds to the triplex minor groove thereby destabilizing the triplex similar as previously reported for netropsin.

Hydration of the RNA duplex r(CGCAAAUUUGCG)2 determined by NMR

The so-called spine of hydration in the minor groove of AnTn tracts in DNA is thought to stabilise the structure, and kinetically bound water detected in the minor groove of such DNA species by NMR has been attributed to a narrow minor groove [Liepinsh, E., Leupin, W. and Otting, G. (1994) Nucleic Acids Res. 22, 2249-2254]. We report here an NMR study of hydration of an RNA dodecamer which has a wide, shallow minor groove. Complete assignments of exchangeable protons, and a large number of non-exchangeable protons in r(CGCAAAUUUGCG)2 have been obtained. In addition, ribose C2'-OH resonances have been detected, which are probably involved in hydrogen bonds. Hydration at different sites in the dodecamer has been measured using ROESY and NOESY experiments at 11.75 and 14.1 T. Base protons in both the major and minor grooves are in contact with water, with effective correlation times for the interaction of approximately 0.5 ns, indicating weak hydration, in contrast to the hydration of adenine C2H in the homologous DNA sequence. NOEs to H1' in the minor groove are consistent with hydration water present that is not observed in the analogous DNA sequence. Hydration kinetics in nucleic acids may be determined by chemical factors such as hydrogen-bonding more than by simple conformational factors such as groove width.

Hydration and solution structure of nucleic acids

Of particular interest among the revelations from recent new DNA structures is the finding that both strands of the repeated DNA sequences found in telomeres and centromeres may adopt alternative conformations. High-definition NMR studies yielded information on the residence time of the water molecules interacting with nucleic acids. A better knowledge of the residence time of the hydration water may be useful in assessing its contribution to nucleic acid structure.