Models of triple-stranded polynucleotides with optimised stereochemistry

Folding Double-Stranded DNA into Designed Shapes with Triplex-Forming Oligonucleotides

The compaction and organization of genomic DNA is a central mechanism in eukaryotic cells, but engineered architectural control over double-stranded DNA (dsDNA) is notably challenging. Here, long dsDNA templates are folded into designed shapes via triplex-mediated self-assembly. Triplex-forming oligonucleotides (TFOs) bind purines in dsDNA via normal or reverse Hoogsteen interactions. In the triplex origami methodology, these non-canonical interactions are programmed to compact dsDNA (linear or plasmid) into well-defined objects, which demonstrate a variety of structural features: hollow and raster-filled, single- and multi-layered, with custom curvatures and geometries, and featuring lattice-free, square-, or honeycomb-pleated internal arrangements. Surprisingly, the length of integrated and free-standing dsDNA loops can be modulated with near-perfect efficiency; from hundreds down to only six bp (2 nm). The inherent rigidity of dsDNA promotes structural robustness and non-periodic structures of almost 25.000 nt are therefore formed with fewer unique starting materials, compared to other DNA-based self-assembly methods. Densely triplexed structures also resist degradation by DNase I. Triplex-mediated dsDNA folding is methodologically straightforward and orthogonal to Watson-Crick-based methods. Moreover, it enables unprecedented spatial control over dsDNA templates.

A single natural RNA modification can destabilize a U•A-T-rich RNA•DNA-DNA triple helix

Recent studies suggest noncoding RNAs interact with genomic DNA, forming RNA•DNA-DNA triple helices, as a mechanism to regulate transcription. One way cells could regulate the formation of these triple helices is through RNA modifications. With over 140 naturally occurring RNA modifications, we hypothesize that some modifications stabilize RNA•DNA-DNA triple helices while others destabilize them. Here, we focus on a pyrimidine-motif triple helix composed of canonical U•A-T and C•G-C base triples. We employed electrophoretic mobility shift assays and microscale thermophoresis to examine how 11 different RNA modifications at a single position in an RNA•DNA-DNA triple helix affect stability: 5-methylcytidine (m⁵C), 5-methyluridine (m⁵U or rT), 3-methyluridine (m³U), pseudouridine (Ψ), 4-thiouridine (s⁴U), N ⁶-methyladenosine (m⁶A), inosine (I), and each nucleobase with 2'-O-methylation (Nm). Compared to the unmodified U•A-T base triple, some modifications have no significant change in stability (Um•A-T), some have ∼2.5-fold decreases in stability (m⁵U•A-T, Ψ•A-T, and s⁴U•A-T), and some completely disrupt triple helix formation (m³U•A-T). To identify potential biological examples of RNA•DNA-DNA triple helices controlled by an RNA modification, we searched RMVar, a database for RNA modifications mapped at single-nucleotide resolution, for lncRNAs containing an RNA modification within a pyrimidine-rich sequence. Using electrophoretic mobility shift assays, the binding of DNA-DNA to a 22-mer segment of human lncRNA Al157886.1 was destabilized by ∼1.7-fold with the substitution of m⁵C at known m⁵C sites. Therefore, the formation and stability of cellular RNA•DNA-DNA triple helices could be influenced by RNA modifications.

Insights into the structural stability of major groove RNA triplexes by WAXS-guided MD simulations

RNA triple helices are commonly observed tertiary motifs that are associated with critical biological functions, including signal transduction. Because the recognition of their biological importance is relatively recent, their full range of structural properties has not yet been elucidated. The integration of solution wide-angle X-ray scattering (WAXS) with molecular dynamics (MD) simulations, described here, provides a new way to capture the structures of major-groove RNA triplexes that evade crystallographic characterization. This method yields excellent agreement between measured and computed WAXS profiles and allows for an atomically detailed visualization of these motifs. Using correlation maps, the relationship between well-defined features in the scattering profiles and real space characteristics of RNA molecules is defined, including the subtle conformational variations in the double-stranded RNA upon the incorporation of a third strand by base triples. This readily applicable approach has the potential to provide insight into interactions that stabilize RNA tertiary structure that enables function.

Molecular structure of a U•A-U-rich RNA triple helix with 11 consecutive base triples

Three-dimensional structures have been solved for several naturally occurring RNA triple helices, although all are limited to six or fewer consecutive base triples, hindering accurate estimation of global and local structural parameters. We present an X-ray crystal structure of a right-handed, U•A-U-rich RNA triple helix with 11 continuous base triples. Due to helical unwinding, the RNA triple helix spans an average of 12 base triples per turn. The double helix portion of the RNA triple helix is more similar to both the helical and base step structural parameters of A′-RNA rather than A-RNA. Its most striking features are its wide and deep major groove, a smaller inclination angle and all three strands favoring a C3′-endo sugar pucker. Despite the presence of a third strand, the diameter of an RNA triple helix remains nearly identical to those of DNA and RNA double helices. Contrary to our previous modeling predictions, this structure demonstrates that an RNA triple helix is not limited in length to six consecutive base triples and that longer RNA triple helices may exist in nature. Our structure provides a starting point to establish structural parameters of the so-called ‘ideal’ RNA triple helix, analogous to A-RNA and B-DNA double helices.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

3D-NuS: A Web Server for Automated Modeling and Visualization of Non-Canonical 3-Dimensional Nucleic Acid Structures

The inherent conformational flexibility of nucleic acids facilitates the formation of a range of conformations such as duplex, triplex, quadruplex, etc. that play crucial roles in biological processes. Elucidation of the influence of non-canonical base pair mismatches on DNA/RNA structures at different sequence contexts to understand the mismatch repair, misregulation of alternative splicing mechanisms and the sequence-dependent effect of RNA-DNA hybrid in relevance to antisense strategy demand their three-dimensional structural information. Furthermore, structural insights about nucleic acid triplexes, which are generally not tractable to structure determination by X-ray crystallography or NMR techniques, are essential to establish their biological function(s). A web server, namely 3D-NuS (http://iith.ac.in/3dnus/), has been developed to generate energy-minimized models of 80 different types of triplexes, 64 types of G-quadruplexes, left-handed Z-DNA/RNA duplexes, and RNA-DNA hybrid duplex along with inter- and intramolecular DNA or RNA duplexes comprising a variety of mismatches and their chimeric forms for any user-defined sequence and length. It also generates an ensemble of conformations corresponding to the modeled structure. These structures may serve as good starting models for docking proteins and small molecules with nucleic acids, NMR structure determination, cryo-electron microscope modeling, DNA/RNA nanotechnology applications and molecular dynamics simulation studies.

Detection of an en masse and reversible B- to A-DNA conformational transition in prokaryotes in response to desiccation

The role that DNA conformation plays in the biochemistry of cells has been the subject of intensive research since DNA polymorphism was discovered. B-DNA has long been considered the native form of DNA in cells although alternative conformations of DNA are thought to occur transiently and along short tracts. Here, we report the first direct observation of a fully reversible en masse conformational transition between B- and A-DNA within live bacterial cells using Fourier transform infrared (FTIR) spectroscopy. This biospectroscopic technique allows for non-invasive and reagent-free examination of the holistic biochemistry of samples. For this reason, we have been able to observe the previously unknown conformational transition in all four species of bacteria investigated. Detection of this transition is evidence of a previously unexplored biological significance for A-DNA and highlights the need for new research into the role that A-DNA plays as a cellular defence mechanism and in stabilizing the DNA conformation. Such studies are pivotal in understanding the role of A-DNA in the evolutionary pathway of nucleic acids. Furthermore, this discovery demonstrates the exquisite capabilities of FTIR spectroscopy and opens the door for further investigations of cell biochemistry with this under-used technique.

The multiple personalities of Watson and Crick strands

Background

In genetics it is customary to refer to double-stranded DNA as containing a "Watson strand" and a "Crick strand." However, there seems to be no consensus in the literature on the exact meaning of these two terms, and the many usages contradict one another as well as the original definition. Here, we review the history of the terminology and suggest retaining a single sense that is currently the most useful and consistent.

Proposal

The Saccharomyces Genome Database defines the Watson strand as the strand which has its 5'-end at the short-arm telomere and the Crick strand as its complement. The Watson strand is always used as the reference strand in their database. Using this as the basis of our standard, we recommend that Watson and Crick strand terminology only be used in the context of genomics. When possible, the centromere or other genomic feature should be used as a reference point, dividing the chromosome into two arms of unequal lengths. Under our proposal, the Watson strand is standardized as the strand whose 5'-end is on the short arm of the chromosome, and the Crick strand as the one whose 5'-end is on the long arm. Furthermore, the Watson strand should be retained as the reference (plus) strand in a genomic database. This usage not only makes the determination of Watson and Crick unambiguous, but also allows unambiguous selection of reference stands for genomics.

Reviewers

This article was reviewed by John M. Logsdon, Igor B. Rogozin (nominated by Andrey Rzhetsky), and William Martin.

Calorimetric and spectroscopic studies of aminoglycoside binding to AT-rich DNA triple helices

Calorimetric and fluorescence techniques were used to characterize the binding of aminoglycosides-neomycin, paromomycin, and ribostamycin, with 5'-dA(12)-x-dT(12)-x-dT(12)-3' intramolecular DNA triplex (x = hexaethylene glycol) and poly(dA).2poly(dT) triplex. Our results demonstrate the following features: (1) UV thermal analysis reveals that the T(m) for triplex decreases with increasing pH value in the presence of neomycin, while the T(m) for the duplex remains unchanged. (2) The binding affinity of neomycin decreases with increased pH, although there is an increase in observed binding enthalpy. (3) ITC studies conducted in two buffers (sodium cacodylate and MOPS) yield the number of protonated drug amino groups (Deltan) as 0.29 and 0.40 for neomycin and paromomycin interaction with 5'-dA(12)-x-dT(12)-x-dT(12)-3', respectively. (4) The specific heat capacity change (DeltaC(p)) determined by ITC studies is negative, with more negative values at lower salt concentrations. From 100 mM to 250 mM KCl, the DeltaC(p) ranges from -402 to -60 cal/(mol K) for neomycin. At pH 5.5, a more positive DeltaC(p) is observed, with a value of -98 cal/(mol K) at 100 mM KCl. DeltaC(p) is not significantly affected by ionic strength. (5) Salt dependence studies reveal that there are at least three amino groups of neomycin participating in the electrostatic interactions with the triplex. (6) FID studies using thiazole orange were used to derive the AC(50) (aminoglycoside concentration needed to displace 50% of the dye from the triplex) values. Neomycin shows a seven fold higher affinity than paromomycin and eleven fold higher affinity than ribostamycin at pH 6.8. (7) Modeling studies, consistent with UV and ITC results, show the importance of an additional positive charge in triplex recognition by neomycin. The modeling and thermodynamic studies indicate that neomycin binding to the DNA triplex depends upon significant contributions from charge as well as shape complementarity of the drug to the DNA triplex Watson-Hoogsteen groove.

Design of sequence-specific bifunctional nucleic acid ligands

Homopyrimidine oligodeoxynucleotides have been covalently linked to intercalating agents. These bifunctional nucleic acid ligands bind to the major groove of DNA at homopurine.homopyrimidine sequences, where they form triple helices. The homopyrimidine oligonucleotide binds parallel to the purine strand of the double helix. Two hydrogen bonds are formed between bases of the oligonucleotide and the purines engaged in Watson-Crick base pairs. The intercalating agent inserts its aromatic ring at the triplex-duplex junction, resulting in a strong stabilization of the triple helical structure. Bifunctional oligonucleotide-intercalator conjugates provide new tools for a selective control of gene expression. In addition, irreversible reactions can be targeted to the oligonucleotide recognition sequence. Cleavage reactions can be induced by a copper-phenanthroline chelate or an ellipticine derivative covalently linked to the triple helix-forming oligonucleotide.

Molecular aspects on the interaction of protoberberine, benzophenanthridine, and aristolochia group of alkaloids with nucleic acid structures and biological perspectives

Alkaloids occupy an important position in chemistry and pharmacology. Among the various alkaloids, berberine and coralyne of the protoberberine group, sanguinarine of the benzophenanthridine group, and aristololactam-beta-d-glucoside of the aristolochia group have potential to form molecular complexes with nucleic acid structures and have attracted recent attention for their prospective clinical and pharmacological utility. This review highlights (i) the physicochemical properties of these alkaloids under various environmental conditions, (ii) the structure and functional aspects of various forms of deoxyribonucleic acid (DNA) (B-form, Z-form, H(L)-form, protonated form, and triple helical form) and ribonucleic acid (RNA) (A-form, protonated form, and triple helical form), and (iii) the interaction of these alkaloids with various polymorphic DNA and RNA structures reported by several research groups employing various analytical techniques like absorbance, fluorescence, circular dichroism, and NMR spectroscopy; electrospray ionization mass spectrometry, thermal melting, viscosity, and DNase footprinting as well as molecular modeling and thermodynamic studies to provide detailed binding mechanism at the molecular level for structure-activity relationship. Nucleic acids binding properties of these alkaloids are interpreted in relation to their biological activity.

Structure, Recognition Properties, and Flexibility of the DNA·RNA Hybrid

Molecular dynamics is used to investigate the properties of the DNA.RNA hybrid in aqueous solution at room temperature. The structure of the hybrid is intermediate between A and B forms but, in general, closer to the canonical A-type helix. All the riboses exhibit North puckerings, while 2'-deoxyriboses exist in North, East, and South puckerings, the latter being the most populated one. The molecular recognition pattern of the DNA.RNA hybrid is a unique combination of those of normal DNA and RNA duplexes. Finally, the results obtained from essential dynamics and stiffness analysis demonstrate the large and very asymmetric flexibility of the hybrid and the strong predilection that each strand (DNA or RNA) has on the nature of their intrinsic motions in the corresponding homoduplexes. The implications of the unique structural and dynamic properties of the DNA.RNA hybrid on the mechanism of cleavage by RNase H are discussed.

Relative Flexibility of DNA and RNA: a Molecular Dynamics Study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.

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.

The structure and dynamics of DNA in the gas phase

The impact of the transfer from water to the vacuum in the duplex DNA has been explored by using long molecular dynamic simulations. In opposition to chemical intuition, it is found that vaporization of DNA, even at high temperatures, does not lead to a total disruption of the double helix. Rather, the DNA duplex preserves gross structural, energetic, and dynamic features of the conformation of the double helix in aqueous solution. Thus, the two strands remain bound, the global structure has a slight helicity, and the total number of DNA-DNA interactions is not dramatically different from that found in solution. The results provide detailed structural and dynamic information useful to complement current mass spectroscopy studies of DNA structures.

Aminoglycoside complexation with a DNA.RNA hybrid duplex: the thermodynamics of recognition and inhibition of RNA processing enzymes

Spectroscopic and calorimetric techniques were employed to characterize and contrast the binding of the aminoglycoside paromomycin to three octamer nucleic acid duplexes of identical sequence but different strand composition (a DNA.RNA hybrid duplex and the corresponding DNA.DNA and RNA.RNA duplexes). In addition, the impact of paromomycin binding on both RNase H- and RNase A-mediated cleavage of the RNA strand in the DNA.RNA duplex was also determined. Our results reveal the following significant features: (i) Paromomycin binding enhances the thermal stabilities of the RNA.RNA and DNA.RNA duplexes to similar extents, with this thermal enhancement being substantially greater in magnitude than that of the DNA.DNA duplex. (ii) Paromomycin binding to the DNA.RNA hybrid duplex induces CD changes consistent with a shift from an A-like to a more canonical A-conformation. (iii) Paromomycin binding to all three octamer duplexes is linked to the uptake of a similar number of protons, with the magnitude of this number being dependent on pH. (iv) The affinity of paromomycin for the three host duplexes follows the hierarchy, RNA.RNA > DNA.RNA >> DNA.DNA. (v) The observed affinity of paromomycin for the RNA.RNA and DNA.RNA duplexes decreases with increasing pH. (vi) The binding of paromomycin to the DNA.RNA hybrid duplex inhibits both RNase H- and RNase A-mediated cleavage of the RNA strand. We discuss the implications of our combined results with regard to the specific targeting of DNA.RNA hybrid duplex domains and potential antiretroviral applications.

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.

Direct photocleavage of HIV-DNA by quinacridine derivatives triggered by triplex formation

Amino-p-quinacridine compounds (PQs) have been shown to stabilize strongly and specifically triple-helical DNA. Moreover, these derivatives display photoactive properties that make them efficient DNA cleavage agents. We exploited these two properties (triplex-specific binding and photoactivity) to selectively cleave a double-stranded (ds)DNA sequence present in the HIV-1 genome. Cleavage was first carried out on a linearized plasmid (3300 bp) containing the HIV polypurine tract (PPT) that allowed targeting by a triplex-forming oligonucleotide (TFO). PQ(3)(), the most active compound of the series, efficiently cleaved double-stranded DNA in the vicinity of the PPT when this sequence had formed a triplex with a 16-mer TFO. Investigation of the cleavage at the molecular level was addressed on a short DNA fragment (56 bp); the photoinduced cleavage by PQ(3)() occurred only in the presence of the triple helix. Nevertheless, unusual cleavage patterns were observed: damage was observed at guanines located 6-9 bp away from the end of the triple helical site. This cleavage is very efficient (up to 60%), does not require alkaline treatment, and is observed on both strands. A quinacridine-TFO conjugate produced the same cleavage pattern. This observation, along with others, excludes the hypothesis of a triplex-induced allosteric binding site of PQ(3 )()adjacent to the damaged sequence and indicates that PQ(3 )()preferentially binds in the vicinity of the 5'-triplex junction. Irradiation in the presence of TFO-conjugates with acridine (an intercalative agent) and with the tripeptide lys-tryp-lys led to a complete inhibition of the photocleavage reaction. These results are interpreted in terms of competitive binding and of electron-transfer quenching. Together with the findings of simple mechanistic investigations, they led to the conclusion that the photoinduced damage proceeds through a direct electron transfer between the quinacridine and the guanines. This study addresses the chemical mechanism leading to strand breakage and characterizes the particular photosensitivity of the HIV-DNA target sequence which could be an oxidative hot spot for addressed photoinduced strand scission by photosensitizers.

Sedimentation analysis of novel DNA structures formed by homo-oligonucleotides

Sedimentation velocity analysis has been used to examine the base-specific structural conformations and unusual hydrogen bonding patterns of model oligonucleotides. Homo-oligonucleotides composed of 8-28 residues of dA, dT, or dC nucleotides in 100 mM sodium phosphate, pH 7.4, at 20 degrees C behave as extended monomers. Comparison of experimentally determined sedimentation coefficients with theoretical values calculated for assumed helical structures show that dT and dC oligonucleotides are more compact than dA oligonucleotides. For dA oligonucleotides, the average width (1.7 nm), assuming a cylindrical model, is smaller than for control duplex DNA whereas the average rise per base (0.34 nm) is similar to that of B-DNA. For dC and dT oligonucleotides, there is an increase in the average widths (1.8 nm and 2.1 nm, respectively) whereas the average rise per base is smaller (0.28 nm and 0.23 nm, respectively). A significant shape change is observed for oligo dC(28) at lower temperatures (10 degrees C), corresponding to a fourfold decrease in axial ratio. Optical density, circular dichroism, and differential scanning calorimetry data confirm this shape change, attributable from nuclear magnetic resonance analysis to i-motif formation. Sedimentation equilibrium studies of oligo dG(8) and dG(16) reveal extensive self-association and the formation of G-quadruplexes. Continuous distribution analysis of sedimentation velocity data for oligo dG(16) identifies the presence of discrete dimers, tetramers, and dodecamers. These studies distinguish the conformational and colligative properties of the individual bases in DNA and their inherent capacity to promote specific folding pathways.

RNA folding: beyond Watson-Crick pairs

Several crystal structures of RNA fragments, alone or in complex with a specific protein, have been recently solved. In addition, the structures of an artificial ribozyme, the leadzyme, and the cleavage product of a human pathogen ribozyme, have extended the structural diversity of ribozyme architectures. The attained set of folding rules and motifs expand the repertoire seen previously in tRNA structures.

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.

Multistranded DNA structures

DNA oligonucleotides can form multistranded helices through either the folding of a single strand or the association of two, three or four strands of DNA. Structures of several new DNA triplexes, G-quartet DNA quadruplexes and I-motif DNA quadruplexes have been reported recently. These structures provide new insights into helix stability and folding, loop conformations and cation interactions.

Can G-C Hoogsteen-wobble pairs contribute to the stability of d(G. C-C) triplexes?

Quantum mechanics, molecular dynamics and statistical mechanics methods are used to analyze the importance of neutral Hoogsteen-wobble G.C pairing in the stabilization of triple helices based on the poly-(G.C-C) trio at neutral pH and low ionic strength. In spite of the existence of a single hydrogen bond, the Hoogsteen-wobble G.C pair is found to be quite stable both in gas phase and solvated DNA. Molecular dynamics simulations of different triplexes based on the d(G.C-C) trio leads to stable structures if the neutral d(G.C-C) steps stabilized by Hoogsteen-wobble pairs are mixed with d(G.C-C+) steps. Finally, high level ab initio calculations and thermodynamic integration techniques are used to determine the relative stability of G.C wobble and G.C imino pairings. It is found that triplexes containing the imino pairing are slightly more stable structures than those with the wobble one, due mainly to a better stacking.

The high-resolution structure of the triplex formed by the GAA/TTC triplet repeat associated with Friedreich's ataxia

Expansions of the triplet repeat, GAA/TTC, inside the first intron of the frataxin gene causes Friedreich's ataxia (FRDA). It was of interest to us to examine whether the FRDA repeat forms an unusual DNA structure, since formation of such structure during replication may cause its expansion. Here, we show that the FRDA repeat forms a triplex in which the TTC strand folds on either side of the same GAA strand. We have determined the high-resolution NMR structures of two intramolecularly folded FRDA triplexes, (GAA)2T4(TTC)2T4(CTT)2 and (GAA)2T4(TTC)2T2CT2(CTT)2 with T.A.T and C+.G.C triads. T4 represents a synthetic loop sequence, whereas T2CT2 is the natural loop-folding sequence of the TTC strand. We have also made use of site-specific 15N-labeling of the cytosine residues to investigate their protonation status and their interaction with other protons. We show that the cytosine residues of the Hoogsteen C+.G pairs in this triplex are protonated close to physiological pH. Therefore, it appears that the triplex formation offers a plausible explanation for the expansion of the GAA/TTC repeats in FRDA.

Nucleic acid vibrational circular dichroism, absorption, and linear dichroism spectra. II. A DeVoe theory approach

The DeVoe polarizability theory is used to calculate vibrational circular dichroism (VCD) and infrared (IR) absorption spectra of four polyribonucleotides: poly(rA) x poly(rU), poly(rU) x poly(rA) x poly(rU), poly(rG) x poly(rC), and poly(rC+) x poly(rI) x poly(rC). This is the first report on the use of the DeVoe theory to calculate VCD, oriented VCD, IR absorption, and IR linear dichroism (LD) spectra of double- and triple-stranded polyribonucleotides. Results are reported for DeVoe theory calculations--within the base-stretching 1750-1550 cm(-1) spectral region--on several proposed multistranded polyribonucleotide geometries. The calculated spectra obtained from these proposed geometries are compared with previously reported measured and calculated VCD and IR spectral results. Base-base hydrogen-bonding effects on the frequencies and magnitudes of the base carbonyl stretching modes are explicitly considered. The good agreements found between calculated and measured spectra are proposed to be further evidence of the usefulness of the DeVoe theory in drawing three-dimensional structural conclusions from measured polyribonucleotide VCD and IR spectra.

Interactions of the antiviral quinoxaline derivative 9-OH-B220 [2, 3-dimethyl-6-(dimethylaminoethyl)- 9-hydroxy-6H-indolo-[2, 3-b]quinoxaline] with duplex and triplex forms of synthetic DNA and RNA

The binding of an antiviral quinoxaline derivative, 2,3-dimethyl- 6 - (dimethylaminoethyl) - 9 - hydroxy - 6H - indolo - [2,3 - b]quinoxaline (9-OH-B220), to synthetic double and triple helical DNA (poly(dA).poly(dT) and poly(dA).2poly(dT)) and RNA (poly(rA). poly(rU) and poly (rA).2poly(rU)) has been characterized using flow linear dichroism (LD), circular dichroism (CD), fluorescence spectroscopy, and thermal denaturation. When either of the DNA structures or the RNA duplex serve as host polymers a strongly negative LD is displayed, consistent with intercalation of the chromophoric ring system between the base-pairs/triplets of the nucleic acid structures. Evidence for this geometry also includes weak induced CD signals and strong increments of the fluorescence emission intensities upon binding of the drug to each of these polymer structures. In agreement with intercalative binding, 9-OH-B220 is found to effectively enhance the thermal stability of both the double and triple helical states of DNA as well as the RNA duplex. In the case of poly(dA).2poly(dT), the drug provides an unusually large stabilization of its triple helical state; upon binding of 9-OH-B220 the triplex-to-duplex equilibrium is shifted towards higher temperature by 52.5 deg. C in a 10 mM sodium cacodylate buffer (pH 7.0) containing 100 mM NaCl and 1 mM EDTA. When triplex RNA serves as host structure, LD indicates that the average orientation angle between the drug chromophore plane and the helix axis of the triple helical RNA is only about 60 to 65 degrees. Moreover, the thermal stabilizing capability, as well as the fluorescence increment, CD inducing power and perturbations of the absorption envelope, of 9-OH-B220 in complex with the RNA triplex are all less pronounced than those observed for the complexes with DNA and duplex RNA. These features indicate binding of 9-OH-B220 in the wide and shallow minor groove of poly(rA).2poly(rU). Based on the present results, some implications for the applications of this low-toxic, antiviral and easily administered drug in an antigene strategy, as well as its potential use as an antiretroviral agent, are discussed.

Effect of selective substitution of 5-bromocytosine on conformation of DNA triple helices

Three triplex DNAs containing 5-bromocytosine[BrC] were studied by vibrational spectroscopy and molecular modelling. Firstly, three oligodeoxypyrimidines of 5'-(TC)3-T4-(BrCT)3 [CBrC], 5'-(TBrC)3-T4-(CT)3 [BrCC] and 5'-(TBrC)3-T4-(BrCT)3 [BrCBrC] were synthesized and then reacted with an oligodeoxypurine of 5'-(AG)3 at pH=4.5 in phosphate buffer respectively to form three comparative hairpin triplex named CY,YC and YY. The results of FT-Raman and IR revealed that YY is almost in A-like form, CY and YC are combinations of A-like form and B-like form, but A-form dominates in CY while B-form is equivalent as A-form in YC. The result is consistent with the theoretical analysis.

Rational design of a triple helix-specific intercalating ligand

DNA triple helices offer new perspectives toward oligonucleotide-directed gene regulation. However, the poor stability of some of these structures might limit their use under physiological conditions. Specific ligands can intercalate into DNA triple helices and stabilize them. Molecular modeling and thermal denaturation experiments suggest that benzo[f]pyrido[3, 4-b]quinoxaline derivatives intercalate into triple helices by stacking preferentially with the Hoogsteen-paired bases. Based on this model, it was predicted that a benzo[f]quino[3,4-b]quinoxaline derivative, which possesses an additional aromatic ring, could engage additional stacking interactions with the pyrimidine strand of the Watson-Crick double helix upon binding of this pentacyclic ligand to a triplex structure. This compound was synthesized. Thermal denaturation experiments and inhibition of restriction enzyme cleavage show that this new compound can indeed stabilize triple helices with great efficiency and specificity and/or induce triple helix formation under physiological conditions.

Stability and structure of model DNA triplexes and quadruplexes and their interactions with small ligands

This review focuses on the structural and thermodynamic characterization of model DNA triplex and quadruplex structures, taking into account effects of stoichiometry and sequence. Methods such as gel electrophoresis, UV melting, and scanning calorimetry, and the results thereof, are described for determination of the thermodynamic stability of such systems. Three classes of triplexes are considered based on the composition of the third strand, while quadruplex systems are limited to those based on the guanine quartet. X-ray crystallography and high resolution NMR studies are also described for these two classes of unusual structures. Ligand binding to triplexes and quadruplexes is also reviewed, with emphasis on specific molecular recognition. The availability of three-dimensional structures for triplex and quadruplex species sets the stage for structure-based development of ligands capable of binding to them specifically. To this end, we consider the application of DOCK, a program for the discovery of small molecules that can recognize macromolecular structures, to the problem of recognizing folded quadruplex structures. Such studies may ultimately lead to pharmaceutically active compounds.

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

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

Modulation of nucleic acid structure by ligand binding: induction of a DNA.RNA.DNA hybrid triplex by DAPI intercalation

The aromatic diamidine, DAPI (4',6-diamidino-2-phenylindole), is used as an important biological and cytological tool since it forms highly fluorescent complexes with nucleic acid duplexes via minor groove-directed/intercalative modes of interaction. In this study, we find that DAPI binding can induce the formation of an RNA-DNA hybrid triplex that would not otherwise form. More specifically, through application of a broad range of spectroscopic, viscometric, and molecular modeling techniques, we demonstrate that DAPI intercalation induces the formation of the poly(dT).poly(rA).poly(dT) hybrid triple helix, a structure which does not form in the absence of the ligand. Using UV mixing studies, we demonstrate that, in the presence of DAPI, the poly(rA).poly(dT) duplex and the poly(dT) single strand form a 1:1 complex (a triplex) that does not form in the absence of DAPI. Through temperature-dependent absorbance measurements, we show that the poly(dT).poly(rA).poly(dT) triplex melts via two distinct transitions: initial conversion of the triplex to the duplex state, with the DAPI remaining bound, followed by denaturation of the duplex-DAPI complex to its component single strands and free DAPI. Using optical melting profiles, we show that DAPI binding enhances the thermal stability of the poly(dT).poly(rA).poly(dT) triplex, an observation consistent with the preferential binding of the ligand to the triplex versus the duplex and single-stranded states. Our differential scanning calorimetric measurements reveal melting of the DAPI-saturated poly(dT).poly(rA).poly(dT) triplex to be associated with a lower enthalpy but greater cooperativity than melting of the corresponding DAPI-saturated poly(rA).poly(dT) duplex. Our flow linear dichroism and viscometric data are consistent with an intercalative mode of binding when DAPI interacts with both the poly(dT).poly(rA).poly(dT) triplex and the poly(rA).poly(dT) duplex. Finally, computer modeling studies suggest that a combination of both stacking and electrostatic interactions between the intercalated ligand and the host nucleic acid play important roles in the DAPI-induced stabilization of the poly(dT).poly(rA).poly(dT) triplex. In the aggregate, our results demonstrate that ligand binding can be used to induce the formation of triplex structures that do not form in the absence of the ligand. This triplex-inducing capacity has potentially important implications in the design of novel antisense, antigene, antiviral, and diagnostic strategies.

Molecular modelling study of the netropsin complexation with a nucleic acid triple helix

A detailed molecular mechanical study has been made on the complexes of netropsin with the double stranded oligonucleotide (dA)12.(dT)12 and with the triple helix (dA)12.(dT)12.(dT)12. The complexes were built using computer graphics and energy refined using JUMNA program. In agreement with circular dichroism experiments we have shown that 3 netropsins can bind the minor grooves of the triple helix and of the double helix. The groove geometry in the duplex and in the triplex is very similar. However a detailed analysis of the energetic terms shows, in agreement with thermal denaturation studies, that the affinity of netropsin toward the double helices is larger than towards triple helices.

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)11 proposed starting structures. Incorporation of these novel triplets into both A-type and B-type natural tiplex 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.