Parallel and antiparallel (G.GC)2 triple helix fragments in a crystal structure

Structure and stability of different triplets involving artificial nucleobases: clues for the formation of semisynthetic triple helical DNA

A triple helical DNA can control gene expression, help in homologous recombination, induce mutations to facilitate DNA repair mechanisms, suppress oncogene formations, etc. However, the structure and function of semisynthetic triple helical DNA are not known. To understand this, various triplets formed between eight artificial nucleobases (P, Z, J, V, B, S, X, and K) and four natural DNA bases (G, C, A, and T) are studied herein by employing a reliable density functional theoretic (DFT) method. Initially, the triple helix-forming artificial nucleobases interacted with the duplex DNA containing GC and AT base pairs, and subsequently, triple helix-forming natural bases (G and C) interacted with artificial duplex DNA containing PZ, JV, BS, and XK base pairs. Among the different triplets formed in the first category, the C-JV triplet is found to be the most stable with a binding energy of about - 31 kcal/mol. Similarly, among the second category of triplets, the Z-GC and V-GC triplets are the most stable. Interestingly, Z-GC and V-GC are found to be isoenergetic with a binding energy of about - 30 kcal/mol. The C-JV, and Z-GC or V-GC triplets are about 12-14 kcal/mol more stable than the JV and GC base pairs respectively. Microsolvation of these triplets in 5 explicit water molecules further enhanced their stability by 16-21 kcal/mol. These results along with the consecutive stacking of the C-JV triplet (C-JV/C-JV) data indicate that the synthetic nucleobases can form stable semisynthetic triple helical DNA. However, consideration of a full-length DNA containing one or more semisynthetic bases or base pairs is necessary to understand the formation of semisynthetic DNA in living cells.

The intrinsic ability of double-stranded DNA to carry out D-loop and R-loop formation

Double-stranded (ds)DNA, not dsRNA, has an ability to form a homologous complex with single-stranded (ss)DNA or ssRNA of homologous sequence. D-loops and homologous triplexes are homologous complexes formed with ssDNA by RecA/Rad51-family homologous-pairing proteins, and are a key intermediate of homologous (genetic/DNA) recombination. R-loop formation independent of transcription (R-loop formation in trans) was recently found to play roles in gene regulation and development of mammals and plants. In addition, the crRNA-Cas effector complex in CRISPR-Cas systems also relies on R-loop formation to recognize specific target. In homologous complex formation, ssDNA/ssRNA finds a homologous sequence in dsDNA by Watson-Crick base-pairing. crRNA-Cas effector complexes appear to actively melt dsDNA to make its bases available for annealing to crRNA. On the other hand, in D-loop formation and homologous-triplex formation, it is likely that dsDNA recognizes the homologous sequence before the melting of its double helix by using its intrinsic molecular function depending on CH2 at the 2'-position of the deoxyribose, and that the major role of RecA is the extension of ssDNA and the holding dsDNA at a position suitable for homology search. This intrinsic dsDNA function would also play a role in R-loop formation. The dependency of homologous-complex formation on 2'-CH2 of the deoxyribose would explain the absence of homologous complex formation by dsRNA, and dsDNA as sole genome molecule in all cellular organisms.

Parallel triplex structure formed between stretched single-stranded DNA and homologous duplex DNA

The interaction between the single-stranded DNA and the homologous duplex DNA is essential for DNA homologous repair. Here, we report that parallel triplex structure can form spontaneously between a mechanically extended ssDNA and a homologous dsDNA in protein-free condition. The triplex has a contour length close to that of a B-form DNA duplex and remains stable after force is released. The binding energy between the ssDNA and the homologous dsDNA in the triplex is estimated to be comparable to the basepairing energy in a B-form dsDNA. As ssDNA is in a similar extended conformation within recombinase-coated nucleoprotein filaments, we propose that the parallel triplex may form and serve as an intermediate during recombinase-catalyzed homologous joint formation.

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Cytosine, the double helix and DNA self-assembly

DNA self-assembly has crucial implications in reading out the genetic information in the cell and in nanotechnological applications. In a recent paper, self-assembled DNA crystals displaying spectacular triangular motifs have been described (Zheng et al., 2009). The authors claimed that their data demonstrate the possibility to rationally design well-ordered macromolecular 3D DNA lattice with precise spatial control using sticky ends. However, the authors did not recognize the fundamental features that control DNA self-assembly in the lateral direction. By analysing available crystallographic data and simulating a DNA triangle, we show that the double helix geometry, sequence-specific cytosine–phosphate interactions and divalent cations are in fact responsible for the precise spatial assembly of DNA.

Designed nucleotide binding motifs

Gaps in the central strand of oligonucleotide triplexes bind nucleoside phosphates tightly. Watson-Crick and Hoogsteen base pairing as design principle yield motifs with high affinity for nucleoside phosphates with A or G as nucleobase, including ATP. The second messenger 3',5'-cAMP is bound with nanomolar affinity. A designed DNA motif accommodates seven nucleotides at a time. The design was implemented for both DNA and RNA.

A non-canonical DNA structure enables homologous recombination in various genetic systems

Homologous recombination, which is critical to genetic diversity, depends on homologous pairing (HP). HP is the switch from parental to recombinant base pairs, which requires expansion of inter-base pair spaces. This expansion unavoidably causes untwisting of the parental double-stranded DNA. RecA/Rad51-catalyzed ATP-dependent HP is extensively stimulated in vitro by negative supercoils, which compensates for untwisting. However, in vivo, double-stranded DNA is relaxed by bound proteins and thus is an unfavorable substrate for RecA/Rad51. In contrast, Mhr1, an ATP-independent HP protein required for yeast mitochondrial homologous recombination, catalyzes HP without the net untwisting of double-stranded DNA. Therefore, we questioned whether Mhr1 uses a novel strategy to promote HP. Here, we found that, like RecA, Mhr1 induced the extension of bound single-stranded DNA. In addition, this structure was induced by all evolutionarily and structurally distinct HP proteins so far tested, including bacterial RecO, viral RecT, and human Rad51. Thus, HP includes the common non-canonical DNA structure and uses a common core mechanism, independent of the species of HP proteins. We discuss the significance of multiple types of HP proteins.

Are there three polynucleotide strands in the catalytic centre of DNA polymerases?

Mitochondrial DNA may undergo large-scale rearrangements, thus leading to diseases. The mechanisms of these rearrangements are still the matter of debates. Several lines of evidence indicate that breakpoints are characterized by direct repeats (DR), one of them being eliminated from the normal genome. Analysis of DR showed their skewed nucleotide content compatible with the formation of known triple helices. Here, I propose a novel mechanism involving the formation of triplex structures that result from the dissociation of the [synthesized repeat-DNA polymerase] complex. Upon binding to the homologous sequence, replication is initiated from the primer bound in a triple helix manner. This feature implies the initiation of replication on the double-stranded DNA from the triple helix primer. Hereby, I review evidences supporting this model. Indeed, all short d(G)-rich primers 10 nucleotides long can be elongated on double-stranded DNA by phage, bacterial, reverse transcriptases and eukaryotic DNA polymerases. Mismatches may be tolerated between the primer and its double-stranded binding site. In contrast to previous studies, evidences for the parallel binding of the triple helix to its homologous strand are provided. This suggest the displacement of the non-template strand by the triple helix primer upon binding within the DNA polymerase catalytic centre. Computer modelling indicates that the triple helix primer lies within the major groove of the double helix, with its 3' hydroxyl end nearby the catalytic amino acids. Taken together, I bring new concepts on DNA rearrangements, and novel features of triple helices and DNA polymerases that can bind three polynucleotide strands similar to RNA polymerases.

Heteroduplex joint formation free of net topological change by Mhr1, a mitochondrial recombinase

Homologous pairing, an essential process for homologous recombination, is the formation of a heteroduplex joint by an invading single-stranded DNA tail and a complementary sequence within double-stranded DNA (dsDNA). The base rotation of the parental dsDNA, to switch from parental base pairs to heteroduplex ones with the invading single-stranded DNA, sterically requires vertical extension between adjacent base pairs, which inevitably induces untwisting of the dsDNA. RecA is a prototype of the RecA/Rad51/Dmc1 family proteins, which promote ATP-dependent homologous pairing in homologous DNA recombination in vivo, except in mitochondria. As predicted by the requirement for the untwisting, dsDNA bound to RecA is extended and untwisted, and homologous pairing by RecA in vitro is extensively stimulated by the negative supercoils of dsDNA substrates. D-loop formation in negatively supercoiled dsDNA, which serves as an assay for homologous pairing, is also catalyzed in an ATP-independent manner by proteins structurally unrelated to RecA, such as Mhr1. Mhr1 is required for yeast mitochondrial DNA recombination instead of RecA family proteins. Inconsistent with the topological requirements, tests for the effects of negative supercoils revealed that Mhr1 catalyzes homologous pairing with relaxed closed circular dsDNA much more efficiently than with negatively supercoiled dsDNA. Topological analyses indicated that neither the process nor the products of homologous pairing by Mhr1 involve a net topological change of closed circular dsDNA. This would be favorable for homologous recombination in mitochondria, where dsDNA is unlikely to be under topological stress toward unwinding. We propose a novel topological mechanism wherein Mhr1 induces untwisting without net topological change.

Interguanine hydrogen-bonding patterns in adducts with water and Zn–purine complexes (purine is 9-methyladenine and 9-methylguanine). Unexpected preference of Zn(II) for adenine-N7 over guanine-N7

Guanine-guanine hydrogen bonding involving the Watson-Crick edge [N(1)H, N(2)H2] of one base and the Hoogsteen edge (N7, O6) of the other is the dominant association pattern in the solid-state structures of two hydrates of 9-ethylguanine (9-EtGH), and in adducts of 9-methylguanine (9-MeGH) with the Zn compounds [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) as well as [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) (9-MeA is 9-methyladenine). The structures of 9-EtGH*2H2O and 9-EtGH*3.5H2O are dominated by polymeric tape structures of the guanine and extended water clusters. In [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) the metalated guanine is involved in hydrogen bonding (GG3 motif) with a free 9-MeGH, which in turn is centrosymmetrically related to itself via hydrogen bonds involving N2H2 and N3 (GG4 motif). In [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) the metalated adenine base interacts via its Watson-Crick edge [N1, N(6)H2] with the sugar edge [N(2)H2, N3] of one of the guanine nucleobases of the GG pair. Crystallization of [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) from an aqueous solution containing 9-MeGH, 9-MeA, and ZnCl2 is fully unexpected in that the anticipated preference of Zn(II) for guanine-N7 is not realized and instead coordination to adenine-N7 is observed. The relevance of [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) and [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) for metal-containing nucleic acid triplex structures is discussed.

Exploration of triple-helical fragments: crystallization and preliminary X-ray diffraction of d(TGGCCTTAAGG)

The nonamer d(GCGAATTCG) and decamer d(GGCCAATTGG), containing one and two overhanging guanines, respectively, form G x GC triplets in their crystal packing. In order to introduce a third subsequent T x AT triplet, the decamer was further extended by one overhanging thymine residue. Two different crystal morphologies of the sequence d(TGGCCTTAAGG) were obtained by hanging-drop vapour diffusion and diffracted to 2.5 and 2.3 A resolution, respectively. However, both crystals belong to the orthorhombic space group P2(1)2(1)2(1), with similar unit-cell parameters. Therefore, the two data sets could be merged to a resolution of 2.4 A with unit-cell parameters a = 26.97, b = 41.12, c = 52.72 A.

A binding site for Pur alpha and Pur beta is structurally unstable and is required for replication in vivo from the rat aldolase B origin

The rat aldolase B promoter acts as a replication origin in vivo, as well as an autonomously replicating sequence (ARS). Here, we examined roles of a polypurine stretch (site PPu) in this origin, which is indispensable to the ARS activity. Purification of site PPu-binding protein revealed that site PPu binds Puralpha and Purbeta, i.e., single-stranded DNA-binding proteins whose roles in replication have been implicated, but less clear. Biochemical analyses showed that site PPu even in a longer DNA fragment is unstable in terms of double-helix, implying that Puralpha/beta may stabilize single-stranded state. Deletion of site PPu from the origin DNA, which was ectopically positioned in the mouse chromosome, significantly reduced replicator activity. Chromatin immunoprecipitation experiments showed that deletion of site PPu abolishes binding of the Puralpha/beta proteins to the origin. These observations suggest functional roles of site PPu and Puralpha/beta proteins in replication initiation.

Synthesis of novel poly(dG)-poly(dG)-poly(dC) triplex structure by Klenow exo- fragment of DNA polymerase I

The extension of the G-strand of long (700 bp) poly(dG)-poly(dC) by the Klenow exo(-) fragment of DNA polymerase I yields a complete triplex structure of the H-DNA type. High-performance liquid chromatography analysis demonstrates that the length of the G-strand is doubled during the polymerase synthesis. Fluorescence resonance energy transfer analysis shows that the 5' ends of the G- and the C-strands, labeled with fluorescein and TAMRA, respectively, are positioned close to each other in the product of the synthesis. Atomic force microscopy morphology imaging shows that the synthesized structures lack single-stranded fragments and have approximately the same length as the parent 700 bp poly(dG)-poly(dC). CD spectrum of the polymer has a large negative peak at 278 nm, which is characteristic of the poly(dG)-poly(dG)-poly(dC) triplex. The polymer is resistant to DNase and interacts much more weakly with ethidium bromide as compared with the double-stranded DNA.

Netropsin interactions in the minor groove of d(GGCCAATTGG) studied by a combination of resolution enhancement and ab initio calculations

The structure of the complex between the minor groove binder netropsin and d(GGCCAATTGG) was determined via single-crystal X-ray techniques. The structure was refined to completion using refmac5.1.24, resulting in a residual R-factor of 20.0% (including 68 water molecules). Using crystal engineering and cryocooling techniques, the resolution could be enhanced to 1.75 A, resulting in an unambiguous determination of the drug conformation and orientation. As previously noticed, bifurcated hydrogen bonds are formed between the amide nitrogen atoms of the drug and the N3 and O2 atoms of A and T base pairs, respectively, clearly cataloging the structure to class I. As the bulky NH2 group on guanine was believed to prevent binding of the drug in the minor groove, the detailed nature of several of the amidinium and guanidinium end contacts were further investigated by ab initio quantum chemical methods.

Structural precursor of the hemideprotonated guanine pair

A pairing scheme (GG7) between two neutral guanine nucleobases taking place via the Watson-Crick faces and involving two bridging water molecules is presented, which upon removal of a H5O2+ entity could convert into the hemideprotonated guanine pair (GG5).

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.

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.

Base composition at mtDNA boundaries suggests a DNA triple helix model for human mitochondrial DNA large-scale rearrangements

Different mechanisms have been proposed to account for mitochondrial DNA (mtDNA) instability based on the presence of short homologous sequences (direct repeats, DR) at the potential boundaries of mtDNA rearrangements. Among them, slippage-mispairing of the replication complex during the asymmetric replication cycle of the mammalian mitochondrial DNA has been proposed to account for the preferential localization of deletions. This mechanism involves a transfer of the replication complex from the first neo-synthesized heavy (H) strand of the DR1, to the DR2, thus bypassing the intervening sequence and producing a deleted molecule. Nevertheless, the nature of the bonds between the DNA strands remains unknown as the forward sequence of DR2, beyond the replication complex, stays double-stranded. Here, we have analyzed the base composition of the DR at the boundaries of mtDNA deletions and duplications and found a skewed pyrimidine content of about 75% in the light-strand DNA template. This suggests the possible building of a DNA triple helix between the G-rich neo-synthesized DR1 and the base-paired homologous G.C-rich DR2. In vitro experiments with the purified human DNA polymerase gamma subunits enabled us to show that the third DNA strand may be used as a primer for DNA replication, using a template with the direct repeat forming a hairpin, with which the primer could initiate DNA replication. These data suggest a novel molecular basis for mitochondrial DNA rearrangements through the distributive nature of the DNA polymerase gamma, at the level of the direct repeats. A general model accounting for large-scale mitochondrial DNA deletion and duplication is proposed. These experiments extend to a DNA polymerase from an eucaryote source the use of a DNA triple helix strand as a primer, like other DNA polymerases from phage and bacterial origins.

Two 1 : 1 binding modes for distamycin in the minor groove of d(GGCCAATTGG)

Single-crystal X-ray structure determinations of the complex between the minor-groove binder distamycin and d(GGCCAATTGG) reveal two 1 : 1 binding modes which differ in the orientation of the drug molecule in the minor groove. The two crystals were grown from different crystallization conditions and found to diffract to 2.38 and 1.85 A, respectively. The structures were refined to completion using SHELXL-93, resulting in a residual R factor of 20.30% for the 2.38-A resolution structure (including 46 water molecules) and 19.74% for the 1.85-A resolution structure (including 74 water molecules). In both orientations, bifurcated hydrogen bonds are formed between the amide nitrogen atoms of the drug and AT base pairs. With a binding site of at least five base pairs, close contacts between the terminal distamycin atoms and guanine amino groups are inevitable. The detailed nature of several of these interactions was further investigated by ab initio quantum chemical methods.

Parallel and antiparallel G*G.C base triplets in pur*pur.pyr triple helices formed with (GA) third strands

Triple helices with G*G.C and A*A.T base triplets with third GA strands either parallel or antiparallel with respect to the homologous duplex strand have been formed in presence of Na (+) or Mg(2+) counterions. Antiparallel triplexes are more stable and can be obtained even in presence of only monovalent Na(+) counterions. A biphasic melting has been observed, reflecting third strand separation around 20 degrees C followed by the duplex -> coil transition around 63 degrees C. Parallel triplexes are far less stable than the antiparallel ones. Their formation requires divalent ions and is observed at low temperature and in high concentration conditions. Different FTIR signatures of G*G.C triplets in parallel and antiparallel triple helices with GA rich third strands have been obtained allowing the identification of such base triplets in triplexes formed by nucleic acids with heterogeneous compositions. Only S-type sugars are found in the antiparallel triplex while some N-type sugar conformation is detected in the parallel triplex.

Initiation of DNA replication by DNA polymerases from primers forming a triple helix

Despite extensive studies on oligonucleotide-forming triple helices, which were discovered in 1957, their possible relevance in the initiation of DNA replication remains unknown. Using sequences forming triple helices, we have developed a DNA polymerisation assay by using hairpin DNA templates with a 3' dideoxynucleotide end and an unpaired 5'-end extension to be replicated. The T7 DNA polymerase successfully elongated nucleotides to the expected size of the template from the primers forming triple helices composed of 9-14 deoxyguanosine-rich residues. The triple helix-forming primer required for this reaction has to be oriented parallel to the homologous sequence of the hairpin DNA template. Substitution of the deoxyguanosine residues by N7 deazadeoxyguanosines in the hairpin of the template prevented primer elongation, suggesting that the formation of a triple helix is a prerequisite for primer elongation. Furthermore, DNA sequencing could be achieved with the hairpin template through partial elongation of the third DNA strand forming primer. The T4 DNA polymerase and the Klenow fragment of DNA polymerase I provided similar DNA elongation to the T7 polymerase-thioredoxin complex. On the basis of published crystallographic data, we show that the third DNA strand primer fits within the catalytic centre of the T7 DNA polymerase, thus underlying this new property of several DNA polymerases which may be relevant to genome rearrangements and to the evolution of the genetic apparatus, namely the DNA structure and replication processes.

Parallel intramolecular DNA triple helix with G and T bases in the third strand stabilized by Zn(2+) ions

We present evidence of formation of an intramolecular parallel triple helix with T*A.T and G*G.C base triplets (where * represents the hydrogen bonding interaction between the third strand and the duplex while. represents the Watson-Crick interactions which stabilize the duplex). The third GT strand, containing seven GpT/TpG steps, targets the polypurine sequence 5'-AGG-AGG-GAG-GAG-3'. The triple helix is obtained by the folding back twice of a 36mer, formed by three dodecamers tethered by hydroxyalkyl linkers (-L-). Due to the design of the oligonucleotide, the third strand orientation is parallel with respect to the polypurine strand. Triple helical formation has been studied in concentration conditions in which native gel electrophoresis experiments showed the absence of intermolecular structures. Circular dichroism (CD) and UV spectroscopy have been used to evidence the triplex structure. A CD spectrum characteristic of triple helical formation as well as biphasic UV and CD melting curves have been obtained in high ionic strength NaCl solutions in the presence of Zn(2+) ions. Specific interactions with Zn(2+) ions in low water activity conditions are necessary to stabilize the parallel triplex.

Extra-helical guanine interactions in DNA

We review the extra-helical guanine interactions present in many oligonucleotide crystals. Very often terminal guanines interact with other guanines in the minor groove of neighboring oligonucleotides through N2 x N3 hydrogen bonds. In other cases the interaction occurs with the help of Ni2+ ions. Guanine/netropsin stacking in the minor groove has also been found. From these studies we conclude that guanine may have multiple extra-helical interactions. In particular it may be considered a very effective minor groove binder, which could be used in the design of sequence selective binding drugs. Interactions through the major groove are seldom encountered, but might be present when DNA is stretched. Such interactions are also analyzed, since they might be important for homologous chromosome pairing during meiosis.

Energetics of the lattice: packing elements in crystals of four-stranded intercalated cytosine-rich DNA molecules

Condensation of single molecules from solution into crystals represents a transition between distinct energetic states. In solution, the atomic interactions within the molecule dominate. In the crystalline state, however, a set of additional interactions are formed between molecules in close contact in the lattice--these are the packing interactions. The crystal structures of d(CCCT), d(TAACCC), d(CCCAAT), and d(AACCCC) have in common a four-stranded intercalated cytosine segment, built by stacked layers of cytosine.cytosine+ (C.C+) base pairs coming from two parallel duplexes that intercalate into each other with opposite polarity. The intercalated cytosine segments in these structures are similar in their geometry, even though the sequences crystallized in different space groups. In the crystals, adenine and thymine residues of the sequences are used to build the three-dimensional crystal lattice by elaborately interacting with symmetry-related molecules. The packing elements observed provide novel insight about the copious ways in which nucleic acid molecules can interact with each other--for example, when folded in more complicated higher order structures, such as mRNA and chromatin.

Radioprobing of a RecA-three-stranded DNA complex with iodine 125: evidence for recognition of homology in the major groove of the target duplex

A fundamental problem in homologous recombination is how homology between DNAs is recognized. In all current models, a recombination protein loads onto a single strand of DNA and scans another duplex for homology. When homology is found, a synaptic complex is formed, leading to strand exchange and a heteroduplex. A novel technique based on strand cleavage by the Auger radiodecay of iodine 125, allows us to determine the distances between (125)I on the incoming strand and the target sugars of the duplex DNA strands in an Escherichia coli RecA protein-mediated synaptic complex. Analysis of these distances shows that the complex represents a post-strand exchange intermediate in which the heteroduplex is located in the center, while the outgoing strand forms a relatively wide helix intertwined with the heteroduplex and located in its minor groove. The structure implies that homology is recognized in the major groove of the duplex.

Heavy metal mutagenicity: insights from bioinorganic model chemistry

The mutagenicity of metal species may be the result of a direct interaction with the target molecule DNA. Possible scenarios leading to nucleobase mispairing are discussed, and selected examples are presented. They include changes in nucleobase selectivity as a consequence of alterations in acid-base properties of nucleobase atoms and groups involved in complementary H bond formation, guanine deprotonation, and stabilization of rare nucleobase tautomers by metal ions. Oxidative nucleobase damage brought about by metal species will not be considered.

Dimerization of double-stranded RNA possessing a 3'-overhanging single-stranded end

A simple system derived from the acceptor stem of tRNA(Ala) is presented which undergoes a pH-dependent dimerization. This is brought about by formation of C+-G-C base triples of the pyrimidine motif type between protonated cytidines of the 3' single-stranded end and regularly paired G-C pairs in the double-helical stem. In addition, an unusual interaction between a protonated adenine and a regular G-C pair is suggested. The equilibrium between monomer and dimer forms can be monitored via NMR spectroscopy and UV melting curve analysis. A dimerization enthalpy of 159 kJ mol(-1) was found at pH 5.0. The system could serve as a model for inter- and intra-molecular association, respectively, of single-stranded and double-helical regions to enable optimal packing of large RNA molecules.

Correction of chromosomal point mutations in human cells with bifunctional oligonucleotides

A sequence-specific genomic delivery system for the correction of chromosomal mutations was designed by incorporating two different binding domains into a single-stranded oligonucleotide. A repair domain (RD) contained the native sequence of the target region. A third strand-forming domain (TFD) was designed to form a triplex by Hoogsteen interactions. The design was based upon the premise that the RD will rapidly form a heteroduplex that is anchored synergistically by the TFD. Deoxyoligonucleotides were designed to form triplexes in the human adenosine deaminase (ADA) and p53 genes adjacent to known point mutations. Transfection of ADA-deficient human lymphocytes corrected the mutant sequence in 1-2% of cells. Neither the RD or TFD individually corrected the mutation. Transfection of p53 mutant human glioblastoma cells corrected the mutation and induced apoptosis in 7.5% of cells.

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions. Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.

Crystal structure of a double-stranded DNA containing a cisplatin interstrand cross-link at 1.63 A resolution: hydration at the platinated site

cis-diamminedichloroplatinum (II) (cisplatin) is a powerful anti-tumor drug whose target is cellular DNA. In the reaction between DNA and cisplatin, covalent intrastrand and interstrand cross-links (ICL) are formed. Two solution structures of the ICL have been published recently. In both models the double-helix is bent and unwound but with significantly different angle values. We solved the crystal structure at 100K of a double-stranded DNA decamer containing a single cisplatin ICL, using the anomalous scattering (MAD) of platinum as a unique source of phase information. We found 47 degrees for double-helix bending and 70 degrees for unwinding in agreement with previous electrophoretic assays. The crystals are stabilized by intermolecular contacts involving two cytosines extruded from the double-helix, one of which makes a triplet with a terminal G.C pair. The platinum coordination is nearly square and the platinum residue is embedded into a cage of nine water molecules linked to the cross-linked guanines, to the two amine groups, and to the phosphodiester backbone through other water molecules. This water molecule organization is discussed in relation with the chemical stability of the ICL.

Parallel and antiparallel triple helices with G, A-containing third strands

A triple helix, formed by a 13 nucleotide (nt) all-purine oligonucleotide, containing six contiguous guanines, oriented parallel to a homopurine strand present in the polypurine tract of Friend leukemia virus, was obtained in 0.1 M LiCl. Its dissociation constant at 25 degrees C, given by electrophoretic titration, of the order of 50 nM, is at least ten times lower than that of the corresponding antiparallel triplex formed on the same target. At 4 degrees C, the parallel orientation of the homopurine strands is favored to the point that the guanine block of 6 nt, present in the 'antiparallel' oligonucleotide, attaches in a parallel fashion to the corresponding block in the target strand, to generate a partial, parallel triplex, that coexists with the antiparallel one.

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.

Dissociation kinetics of RecA protein-three-stranded DNA complexes reveals a low fidelity of RecA-assisted recognition of homology

We determined that the incorporation of one mismatch into RecA mediated synaptic complexes between oligonucleotide single-stranded DNAs and target duplex DNAs destabilizes the complex by 0.8 to 1.9 kcal/mol. This finding supports our previous result, that RecA binding per se can significantly decrease the loss in free energy associated with mismatch incorporation even in the absence of ATP hydrolysis. We show that the specificity is mostly driven by the dissociation process. We found that the relative destabilization induced by different mismatches depends on their position. Thus, while there is a good correlation between the ranking order of mismatches at the 5' end of synaptic complexes and mismatches in heteroduplexes (D-loops), there is no correlation between the ranking order for mismatches at the 3' end and mismatches in various DNA structures. This difference between the 5' and 3' ends of synaptic complexes agrees well with the established 5' to 3' polarity of the strand exchange promoted by RecA protein. The lack of a correlation between mismatches at the 3' end of synaptic complexes and mismatches in D-loops suggests the intermediate is probably not a canonical protein-free D-loop.

Synthesis and triplex forming properties of an acyclic N7-glycosylated guanine nucleoside

A chiral acyclic nucleoside, one in which the ribose carbohydrate has been replaced with a glycerol-based linker, is prepared by glycosylating guanine at the N7-nitrogen. The stereochemically pure derivative is converted to a DMT-protected phosphoramidite for incorporation into DNA sequences. Sequence containing the acyclic N7-dG nucleoside are capable of forming DNA triplexes in which it is likely that the N1-H and N2-amino groups of the N7-dG are involved in recognition of the guanine base in G-C base pairs.

Stabilizing and destabilizing effects of intercalators on DNA triplexes

Oligonucleotide-directed triplex formation attracts much attention due to its potential usefulness in diagnostic and biotechnological applications. Among other aspects, the research embraces numerous studies probing the influence of intercalating ligands on triplex stability. The effect of the intercalator on triplex formation and stability is known to depend on nucleotide sequence, type of intercalator and solution conditions. The present work is aimed at determining the average number of intercalated ethidium bromide (EtBr) and acridine orange (AO) molecules leading to the most effective stabilization of triplexes. First, fluorescing complexes of intramolecular parallel (recombinant) triplex 5'-d(CATGCTAACT)-L-d(AGTTAGCATG)-L-d(CATGCTAACT)-3' (parARB) and classical antiparallel 5'-(dA)10-L-(dT)10-L-(dT)10-3'(antiATT) (L = -pO(CH2CH2O)3p-) with EtBr and AO were characterized, binding constants were obtained and compared to those for homologous DNA duplexes. Then the total EtBr and AO concentrations corresponding to an average of one, two or three intercalated molecules per oligonucleotide were estimated. Thermal denaturation of parARB and antiATT complexes with an average of one, two or three bound molecules was carried out, thermodynamic parameters of the triplex-to-duplex and duplex-to-open-strand transitions were evaluated using a three-state model. The ability of EtBr and AO to stabilize or destabilize both parallel (recombinant) and classical antiparallel triplexes was found to depend strongly on the concentration of bound intercalator. The triplexes were shown to be stabilized by intercalation of the first and second EtBr or AO molecules, while binding of the third intercalator molecule to 10 nucleotide long triplex resulted in significant triplex destabilization.

A model for the [C+-GxC]n triple helix derived from observation of the C+-GxC base triplet in a crystal structure

A molecular modelling study on the [C+-GxC]n triple helix is reported. We have observed the C+-GxC base triplet in the crystal structure of an oligonucleotide-drug complex, between the minor-groove drug netropsin and the decanucleotide d(CGCAATTGCG)2. The complex was crystallised at pH 7.0, but the crystal structure, at a resolution of 2.4 A, shows that a terminal cytosine has become protonated and participates in a parallel C+-GxC base triplet. The structure of this triplet and its associated sugar-phosphate backbones have been energy-refined and then used to generate a triple helix. This has characteristics of the B-type family of DNA structures for two strands, with the third, the C+ strand, having backbone conformations closer to the A family.

Correlation of hydrophobicity and packing in A-DNA oligonucleotides

A-DNA oligomers pack in a slanted fashion with the terminal base pairs abutting into the minor groove of neighboring molecules unlike the other forms of DNA which pack by vertically stacking one over the other into helical columns. To explain the differences in packing we have advanced a hypothesis that the orientation of the sugar-phosphate backbone is different in A-DNA from that in the other forms of DNA, mainly due to the differences in the sugar puckering.

Crystal structure of d(GCGCGCG) with 5'-overhang G residues

The crystal structure of the DNA heptamer d(GCGCGCG) has been solved at 1.65 A resolution by the molecular replacement method and refined to an R-value of 0.184 for 3598 reflections. The heptamer forms a Z-DNA d(CGCGCG)2 with 5'-overhang G residues instead of an A-DNA d(GCGCGC)2 with 3'-overhang G residues. The overhang G residues from parallel strands of two adjacent duplexes form a trans reverse Hoogsteen G x G basepair that stacks on the six Z-DNA basepairs to produce a pseudocontinuous helix. The reverse Hoogsteen G x G basepair is unusual in that the displacement of one G base relative to the other allows them to participate in a bifurcated (G1)N2 . . . N7(G8) and an enhanced (G8)C8-H . . . O6(G1) hydrogen bond, in addition to the two usual hydrogen bonds. The 5'-overhang G residues are anti and C2'-endo while the 3'-terminal G residues are syn and C2'-endo. The conformations of both G residues are different from the syn/C3'-endo for the guanosine in a standard Z-DNA. The two cobalt hexammine ions bind to the phosphate groups in both GpC and CpG steps in Z(I) and Z(II) conformations. The water structure motif is similar to the other Z-DNA structures.

Progress in developments of triplex-based strategies

Recognition of B-DNA by oligonucleotides that form triple helices is a unique method to specifically recognize sequences of double-stranded DNA. Recently, some significant limitations of the triple-based applications have been overcome. Stable intermolecular triplexes can be formed under physiologic conditions. Binding affinities of modified oligonucleotides to their target sequence due to Hoogsteen or reverse Hoogsteen hydrogen bonding interactions are now in the range of those obtained for duplex formation via Watson-Crick hydrogen bonding interactions even if the kinetics may be quite different. Progress has been made toward developing general procedures to determine the molecular mechanisms of action of triplex-forming oligonucleotides (TFO) administered to cultured cells to provide a rational proof-of-concept for antigene strategies. The antigene strategy has reached a point where TFOs can be used to interfere with several biologic progresses (replication, transcription, recombination, repair) in relevant systems both in vitro and ex vivo.

DNA damage-dependent transcriptional arrest and termination of RNA polymerase II elongation complexes in DNA template containing HIV-1 promoter

We have developed a new biochemical method to isolate a homogeneous population of RNA polymerase II (RNA pol II) elongation complexes arrested at a DNA damage site. The method involves triple-helix formation at a predetermined site in DNA template with a third strand labeled with psoralen at its 5'-end and a biotin at the 3'-end. After triplex formation and near-ultraviolet irradiation (360 nm), DNA templates modified with psoralen were immobilized on streptavidin-coated magnetic beads and used for in vitro transcription reactions with HeLa nuclear extracts. Separation of magnetic beads from solution results in isolation of arrested elongation complexes on the immobilized DNA templates. We have applied the method to arrest RNA pol II elongation complexes on a DNA template containing HIV-1 promoter. Our results indicate that psoralen crosslink in the template strand efficiently arrests elongation complexes, and psoralen monoadducts terminate transcription. Our results also demonstrate that a triple-helical structure stabilized by an intercalator, acridine, attached to the third strand of the helix inhibits transcription by a termination pathway. Isolation of stable RNA pol II elongation complexes arrested at DNA damage sites is a remarkable finding. This result demonstrates that arrested elongation complexes are impervious to DNA damage repair machinery and other regulatory proteins present in HeLa nuclear extracts. The method of delivering site-specific psoralen damage by a triplex structure and isolation of arrested RNA pol II elongation complexes should be generalizable to any promoter and DNA template sequences. This strategy provides a new approach to study the mechanism of transcription elongation and transcription-coupled DNA damage repair.

Unusual DNA conformations

DNA is on the move across conformational space. Duplexes diversity and, joined by triplexes, quadruplexes, loops, bulges and multiarmed junctions, open the route to a bewildering array of increasingly complex conformations. In addition to this structural growth, DNA has come under increasing scrutiny thanks to the development of chemical and physical techniques for deforming its conformation and probing its properties. These investigations help us to learn more about the mechanics and the activity of this remarkably versatile macromolecule.