Four-stranded DNA structures can be stabilized by two different types of minor groove G:C:G:C tetrads

Four-stranded nucleic acid structures are central to many processes in biology and in supramolecular chemistry. It has been shown recently that four-stranded DNA structures are not only limited to the classical guanine quadruplex but also can be formed by tetrads resulting from the association of Watson-Crick base pairs. Such an association may occur through the minor or the major groove side of the base pairs. Structures stabilized by minor groove tetrads present distinctive features, clearly different from the canonical guanine quadruplex, making these quadruplexes a unique structural motif. Within our efforts to study the sequence requirements for the formation of this unusual DNA motif, we have determined the solution structure of the cyclic oligonucleotide dpCCGTCCGT by two-dimensional NMR spectroscopy and restrained molecular dynamics. This molecule self-associates, forming a symmetric dimer stabilized by two G:C:G:C tetrads with intermolecular G-C base pairs. Interestingly, although the overall three-dimensional structure is similar to that found in other cyclic and linear oligonucleotides of related sequences, the tetrads that stabilize the structure of dpCCGTCCGT are different to other minor groove G:C:G:C tetrads found earlier. Whereas in previous cases the G-C base pairs aligned directly, in this new tetrad the relative position of the two base pairs is slipped along the axis defined by the base pairs. This is the first time that a quadruplex structure entirely stabilized by slipped minor groove G:C:G:C tetrads is observed in solution or in the solid state. However, an analogous arrangement of G-C base pairs occurs between the terminal residues of contiguous duplexes in some DNA crystals. This structural polymorphism between minor groove GC tetrads may be important in stabilization of higher order DNA structures.     

Four-stranded DNA structure stabilized by a novel G:C:A:T tetrad

The solution structure of a cyclic oligonucleotide d has been determined by two-dimensional NMR spectroscopy and restrained molecular dynamics. Under the appropriate experimental conditions, this molecule self-associates, forming a symmetric dimer stabilized by four intermolecular Watson-Crick base pairs. The resulting four-stranded structure consists of two G:C:A:T tetrads, formed by facing the minor groove side of the Watson-Crick base-pairs. Most probably, the association of the base-pairs is stabilized by coordinating a Na(+) cation. This is the first time that this novel G:C:A:T tetrad has been found in an oligonucleotide structure. This observation increases considerably the number of sequences that may adopt a four-stranded architecture. Overall, the three-dimensional structure is similar to those observed previously in other quadruplexes formed by minor groove alignment of Watson-Crick base pairs. This resemblance strongly suggests that we may be observing a general motif for DNA-DNA recognition.     

Ground-state recovery following UV excitation is much slower in G x C-DNA duplexes and hairpins than in mononucleotides

Excited states in double-stranded oligonucleotides containing G.C base pairs were studied by femtosecond transient absorption spectroscopy. Relaxation to the electronic ground state occurs about 10 times more slowly in the duplexes and hairpins studied on average than in the individual mononucleotides of G and C. Detection of long-lived excited states in G.C oligonucleotides complements the earlier observation of slow ground-state recovery in A.T DNA, showing that excited states with picosecond lifetimes are formed in DNAs containing either kind of base pair. The results show further that Watson-Crick G.C base pairs in these base-paired and base-stacked duplexes do not enable subpicosecond relaxation to the electronic ground state. A model is proposed in which fluorescent exciton states decay rapidly and irreversibly to dark exciplex states. This model explains the seemingly contradictory observations of femtosecond fluorescence and slower, picosecond recovery of the ground-state population.     

Spectroscopic characterization of interstrand carbinolamine cross-links formed in the 5'-CpG-3' sequence by the acrolein-derived gamma-OH-1,N2-propano-2'-deoxyguanosine DNA adduct

The interstrand N2,N2-dG DNA cross-linking chemistry of the acrolein-derived gamma-OH-1,N2-propanodeoxyguanosine (gamma-OH-PdG) adduct in the 5'-CpG-3' sequence was monitored within a dodecamer duplex by NMR spectroscopy, in situ, using a series of site-specific 13C- and 15N-edited experiments. At equilibrium 40% of the DNA was cross-linked, with the carbinolamine form of the cross-link predominating. The cross-link existed in equilibrium with the non-crosslinked N2-(3-oxo-propyl)-dG aldehyde and its geminal diol hydrate. The ratio of aldehyde/diol increased at higher temperatures. The 1,N2-dG cyclic adduct was not detected. Molecular modeling suggested that the carbinolamine linkage should be capable of maintaining Watson-Crick hydrogen bonding at both of the tandem C x G base pairs. In contrast, dehydration of the carbinolamine cross-link to an imine (Schiff base) cross-link, or cyclization of the latter to form a pyrimidopurinone cross-link, was predicted to require disruption of Watson-Crick hydrogen bonding at one or both of the tandem cross-linked C x G base pairs. When the gamma-OH-PdG adduct contained within the 5'-CpG-3' sequence was instead annealed into duplex DNA opposite T, a mixture of the 1,N2-dG cyclic adduct, the aldehyde, and the diol, but no cross-link, was observed. With this mismatched duplex, reaction with the tetrapeptide KWKK formed DNA-peptide cross-links efficiently. When annealed opposite dA, gamma-OH-PdG remained as the 1,N2-dG cyclic adduct although transient epimerization was detected by trapping with the peptide KWKK. The results provide a rationale for the stability of interstrand cross-links formed by acrolein and perhaps other alpha,beta-unsaturated aldehydes. These sequence-specific carbinolamine cross-links are anticipated to interfere with DNA replication and contribute to acrolein-mediated genotoxicity.     

Hydrogen-bonded nucleic acid base pairs containing unusual base tautomers: complete basis set calculations at the MP2 and CCSD(T) levels

The total interaction energies of altogether 15 hydrogen-bonded nucleic acid base pairs containing unusual base tautomers were calculated. The geometry properties of all selected adenine-thymine and guanine-cytosine hydrogen-bonded base pairs enable their incorporation into DNA. Unusual base pairing patterns were compared with Watson-Crick H-bonded structures of the adenine-thymine and guanine-cytosine pairs. The complete basis set (CBS) limit of the MP2 interaction energy and the CCSD(T) correction term, determined as the difference between the CCSD(T) and MP2 interaction energies, was evaluated. Extrapolation to the MP2 CBS limit was done using the aug-cc-pVDZ and aug-cc-pVTZ results, and the CCSD(T) correction term was determined with the 6-31G*(0.25) basis set. Final interaction energies were corrected while taking into account both tautomeric penalization determined at the CBS level and solvation/desolvation free energies. The situation for the adenine-thymine pairs is straightforward, and tautomeric pairs are significantly less stable than the Watson-Crick pair consisting of the canonical forms. In the case of the guanine-cytosine pair, the Watson-Crick structure made by canonical forms is again the most stable. The other two structures are, however, energetically rather similar (by 5 and 6 kcal/mol), which provides a very small but non-negligible chance of detecting these structures in the DNA double helix (1:5000). Due to the fact that DNA bases and base pairs incorporated into DNA are solvated less favorably than in isolated systems, this probability represents the very upper limit. The results clearly show how precisely the canonical building blocks of DNA molecules were chosen and how well their stability is maintained.     

How well can new-generation density functional methods describe stacking interactions in biological systems?

We compare the performance of four recently developed DFT methods (MPW1B95, MPWB1K, PW6B95, and PWB6K) and two previous, generally successful DFT methods (B3LYP and B97-1) for the calculation of stacking interactions in six nucleic acid bases complexes and five amino acid pairs and for the calculation of hydrogen bonding interactions in two Watson-Crick type base pairs. We found that the four newly developed DFT methods give reasonable results for the stacking interactions in the DNA base pairs and amino acid pairs, whereas the previous DFT methods fail to describe interactions in these stacked complexes. We conclude that the new generation of DFT methods have greatly improved performance for stacking interaction as compared to previously available methods. We recommend the PWB6K method for investigating large DNA or protein systems where stacking plays an important role.     

Formation of an amino-acid-binding pocket through adaptive zippering-up of a large DNA hairpin loop

Background: In vitro selection has identified DNA aptamers that target cofactors, amino acids, peptides and proteins. Structure determination of such ligand-DNA aptamer complexes should elucidate the details of adaptive DNA structural transitions, binding-pocket architectures and ligand recognition. We have determined the solution structure of the complex of a DNA aptamer containing a guanine-rich 18-residue hairpin loop that binds l-argininamide with ∼ 100μM affinity.Results: The DNA aptamer generates its l-argininamide-binding pocket by adaptive zippering up the 18-residue loop through formation of Watson-Crick pairs, mismatch pairs and base triples, while maximizing stacking interactions. Three of the four base triples involve minor-groove recognition through sheared G·A mismatch formation. The unique fold is also achieved through positioning of an adenine residue deep within the minor groove and through nestling of a smaller loop within the larger loop on complex formation. The accessibility to the unique l-argininamide-binding pocket is restricted by a base pair that bridges across one side of the major-groove-binding site. The guanidinium group of the bound l-argininamide aligns through intermolecular hydrogen-bond formation with the base edges of nonadjacent guanine and cytosine residues while being sandwiched between the planes of nonadjacent guanine residues.Conclusions: The available structures of l-arginine/l-argininamide bound to their DNA and RNA targets define the common principles and patterns associated with molecular recognition, as well as the diversity of intermolecular hydrogen-bonding alignments associated with the distinct binding pockets.     

Oligonucleotide cross-linking with copper ions. Spectral and quantum chemical study

Copper(II) complexes with synthetic oligonucleotides consisting of repeating adenine–thymine and guanine–cytosine complementary base pairs have been studied by UV spectroscopy and simulated by DFT quantum chemical calculations at the B3LYP/6-311G++(d,p) level of theory with inclusion of solvation (hydration) effects. The obtained data suggest selective interaction of copper(II) ions with guanine–cytosine complementary pairs, followed by DNA cross-linking at those sites.     

On the effect of covalently appended quinolones on termini of DNA duplexes

Quinolones are gyrase inhibitors that are widely used as antibiotics in the clinic. When covalently attached to oligonucleotides as 5'-acylamido substituents, quinolones were found to stabilize duplexes of oligonucleotides against thermal denaturation. For short duplexes, such as qu-T*GCGCA, where qu is a quinolone residue and T is a 5'-amino-5'-deoxythymidine residue, an increase in the UV melting point of up to 27.8 degrees C was measured. The stabilizing effect was demonstrated for all quinolones tested, namely nalidixic acid, oxolinic acid, pipemidic acid, cinoxacin, norfloxacin, and ofloxacin. The three-dimensional structure of (oa-T*GCGCA)2, where oa is an oxolinic acid residue, was solved by two-dimensional NMR spectroscopy and restrained molecular dynamics. In this complex, the oxolinic acid residues disrupt the terminal T1:A6 base pairs and stack on the G2:C5 base pairs. The displaced adenosine residues bind in the minor groove of the core duplex, while the thymidine residues pack against the oxolinic acid residues. The "molecular cap" thus formed fits tightly on the G:C base pairs, resulting in increased base-pairing fidelity, as demonstrated in UV melting experiments with the sequence oa-T*GGTTGAC and target strands containing a mismatched nucleobase. The structure of the "molecular cap" with its disrupted terminal base pair may also be helpful for modeling how quinolones block re-ligation of DNA strands in the active site of gyrases.     

MP2 study on the hydrogen bonding interaction between 5-hydroxy-5-methylhydantoin and DNA bases: A,C, G,T

The 5-hydroxy-5-methylhydantoin (5-OH-5-Me-dHyd) is a nucleobase lesion induced by the action of ionizing radiation on thymine residue in DNA. In this study, we present the hydrogen bonding base pairs involving 5-OH-5-Me-dHyd bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys–Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: 5-OH-5-Me-dHyd:G>5-OH-5-Me-dHyd:A>5-OH-5-Me -dHyd:C~5-OH-5-Me-dHyd:T.     

Examination of the Effect of the Annealing Cation on Higher Order Structures Containing Guanine or Isoguanine Repeats

Isoguanine (2‐oxo‐6‐amino‐guanine), a natural but non‐standard base, exhibits unique self‐association properties compared to its isomer, guanine, and results in formation of different higher order DNA structures. In this work, the higher order structures formed by oligonucleotides containing guanine repeats or isoguanine repeats after annealing in solutions containing various cations are evaluated by electrospray ionization mass spectrometry (ESI‐MS) and circular dichroism (CD) spectroscopy. The guanine‐containing strand (G9) consistently formed quadruplexes upon annealing, whereas the isoguanine strand (Ig9) formed both pentaplexes and quadruplexes depending on the annealing cation. Quadruplex formation with G9 showed some dependence on the identity of the cation present during annealing with high relative quadruplex formation detected with six of ten cations. Analogous annealing experiments with Ig9 resulted in complex formation with all ten cations, and the majority of the resulting complexes were pentaplexes. CD results indicated most of the original complexes survived the desalting process necessary for ESI‐MS analysis. In addition, several complexes, especially the pentaplexes, were found to be capable of cation exchange with ammonium ions. Ab initio calculations were conducted for isoguanine tetrads and pentads coordinated with all ten cations to predict the most energetically stable structures of the complexes in the gas phase. The observed preference of forming quadruplexes versus pentaplexes as a function of the coordinated cation can be interpreted by the calculated reaction energies of both the tetrads and pentads in combination with the distortion energies of tetrads.     

MP2 study on the hydrogen-bonding interaction between 5-hydroxymethyl-uracil and DNA bases: A,C, G,T

The 5-hydroxymethyl-uracil (HmU) is a product of oxidative attack on the methyl group of thymine in DNA. In this work, we present the hydrogen bonding complexes formation involving HmU bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys-Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: HmU:A > HmU:G > HmU:C > HmU:T.     

Metal–base pairing in DNA

The use of DNA as a molecular wire in nanoscale electronic architectures would greatly benefit from its capability of sequence-specific self-assembly. Although single electrons and positive charges have been shown to be transmitted by natural DNA over a distance of several base pairs, the high ohmic resistance of unmodified oligonucleotides imposes a serious obstacle. Exchanging some or all of the Watson–Crick base pairs in DNA by metal complexes may solve this problem and evolve DNA-like materials with superior conductivity for future nano-electronic applications. The so-called metal–base pairs are formed from suitable transition metal ions and ligand-like nucleosides which are introduced into both of the two pairing strands by automated DNA synthesis. This review illustrates the basic concepts of metal–base pairing and highlights recent developments in the field.     

Lipophilic oligonucleotides spontaneously insert into lipid membranes, bind complementary DNA strands, and sequester into lipid-disordered domains

For the development of surface functionalized bilayers, we have synthesized lipophilic oligonucleotides to combine the molecular recognition mechanism of nucleic acids and the self-assembly characteristics of lipids in planar membranes. A lipophilic oligonucleotide consisting of 21 thymidine units and two lipophilic nucleotides with an alpha-tocopherol moiety as a lipophilic anchor was synthesized using solid-phase methods with a phosphoramadite strategy. The interaction of the water soluble lipophilic oligonucleotide with vesicular lipid membranes and its capability to bind complementary DNA strands was studied using complementary methods such as NMR, EPR, DSC, fluorescence spectroscopy, and fluorescence microscopy. This oligonucleotide inserted stably into preformed membranes from the aqueous phase. Thereby, no significant perturbation of the lipid bilayer and its stability was observed. However, the non-lipidated end of the oligonucleotide is exposed to the aqueous environment, is relatively mobile, and is free to interact with complementary DNA strands. Binding of the complementary single-stranded DNA molecules is fast and accomplished by the formation of Watson-Crick base pairs, which was confirmed by 1H NMR chemical shift analysis and fluorescence resonance energy transfer. The molecular structure of the membrane bound DNA double helix is very similar to the free double-stranded DNA. Further, the membrane bound DNA double strands also undergo regular melting. Finally, in raft-like membrane mixtures, the lipophilic oligonucleotide was shown to preferentially sequester into liquid-disordered membrane domains.     

Novel benzopyrimidines as widened analogues of DNA bases

We report on the synthesis, stacking, and pairing properties of a new structural class of size-expanded pyrimidine nucleosides, abbreviated dyT and dyC. Their bases are benzo-homologated variants of thymine and cytosine and have a design that is distinct from a previously described class of size-expanded (xDNA) pyrimidines, with a different vector of expansion relative to the sugar. We term this new base geometry "yDNA" (a mnemonic for "wide DNA"). Both C-glycosides were prepared using Pd-mediated coupling of iodinated base derivatives with a deoxyribose precursor. As free deoxynucleosides, both dyT and dyC displayed robust fluorescence, with emission maxima at 375 and 390 nm, respectively. Both widened pyrimidines could be incorporated readily as protected phosphoramidite derivatives into synthetic oligonucleotides. Experiments in "dangling end" DNA contexts revealed that both yT and yC stack more favorably than their natural counterparts. When opposite natural bases in the context of Watson-Crick DNA were paired, the yT nucleotide formed a pair with A that was equally stable as a T-A pair, despite the mismatch in size with the neighboring natural pairs. The yC nucleotide (paired opposite G) was destabilizing by a small amount in the same context. Despite the large size of the pairs, both yT and yC were selective for their Watson-Crick complementary partners A and G, respectively. The pairing properties and fluorescence of yDNA nucleotides may lead to useful applications in the study of steric effects in DNA-protein interactions. In addition, the compounds may serve as building blocks for a large-sized artificial genetic system.     

A Novel DNA Helical Wire Containing HgII‐Mediated T:T and T:G Pairs

Numerous applications of metal‐mediated base pairs (metallo‐base‐pairs) to nucleic acid based nanodevices and genetic code expansion have been extensively studied. Many of these metallo‐base‐pairs are formed in DNA and RNA duplexes containing Watson–Crick base pairs. Recently, a crystal structure of a metal–DNA nanowire with an uninterrupted one‐dimensional silver array was reported. We now report the crystal structure of a novel DNA helical wire containing Hg^II‐mediated T:T and T:G base pairs and water‐mediated C:C base pairs. The Hg‐DNA wire does not contain any Watson–Crick base pairs. Crystals of the Hg‐DNA wire, which is the first DNA wire structure driven by Hg^II ions, were obtained by mixing the short oligonucleotide d(TTTGC) and Hg^II ions. This study demonstrates the potential of metallo‐DNA to form various structural components that can be used for functional nanodevices.     

Structural basis of DNA folding and recognition in an AMP-DNA aptamer complex: distinct architectures but common recognition motifs for DNA and RNA aptamers complexed to AMP

Background: Structural studies by nuclear magnetic resonance (NMR) of RNA and DNA aptamer complexes identified through in vitro selection and amplification have provided a wealth of information on RNA and DNA tertiary structure and molecular recognition in solution. The RNA and DNA aptamers that target ATP (and AMP)' with micromolar affinity exhibit distinct binding site sequences and secondary structures. We report below on the tertiary structure of the AMP-DNA aptamer complex in solution and compare it with the previously reported tertiary structure of the AMP-RNA aptamer complex in solution.Results: The solution structure of the AMP-DNA aptamer complex shows, surprisingly, that two AMP molecules are intercalated at adjacent sites within a rectangular widened minor groove. Complex formation involves adaptive binding where the asymmetric internal bubble of the free DNA aptamer zippers up through formation of a continuous six-base mismatch segment which includes a pair of adjacent three-base platforms. The AMP molecules pair through their Watson-Crick edges with the minor groove edges of guanine residues. These recognition G·A mismatches are flanked by sheared G·A and reversed Hoogsteen G·G mismatch pairs.Conclusions: The AMP-DNA aptamer and AMP-RNA aptamer complexes have distinct tertiary structures and binding stoichiometries. Nevertheless, both complexes have similar structural features and recognition alignments in their binding pockets. Specifically, AMP targets both DNA and RNA aptamers by intercalating between purine bases and through identical G·A mismatch formation. The recognition G·A mismatch stacks with a reversed Hoogsteen G·G mismatch in one direction and with an adenine base in the other direction in both complexes. It is striking that DNA and RNA aptamers selected independently from libraries of 10¹⁴ molecules in each case utilize identical mismatch alignments for molecular recognition with micromolar affinity within binding-site pockets containing common structural elements.     

Extending the language of DNA molecular recognition by polyamides: unexpected influence of imidazole and pyrrole arrangement on binding affinity and specificity

Pyrrole (Py) and imidazole (Im) polyamides can be designed to target specific DNA sequences. The effect that the pyrrole and imidazole arrangement, plus DNA sequence, have on sequence specificity and binding affinity has been investigated using DNA melting (DeltaT(M)), circular dichroism (CD), and surface plasmon resonance (SPR) studies. SPR results obtained from a complete set of triheterocyclic polyamides show a dramatic difference in the affinity of f-ImPyIm for its cognate DNA (K(eq) = 1.9 x 10(8) M(-1)) and f-PyPyIm for its cognate DNA (K(eq) = 5.9 x 10(5) M(-1)), which could not have been anticipated prior to characterization of these compounds. Moreover, f-ImPyIm has a 10-fold greater affinity for CGCG than distamycin A has for its cognate, AATT. To understand this difference, the triamide dimers are divided into two structural groupings: central and terminal pairings. The four possible central pairings show decreasing selectivity and affinity for their respective cognate sequences: -ImPy > -PyPy- > -PyIm- approximately -ImIm-. These results extend the language of current design motifs for polyamide sequence recognition to include the use of "words" for recognizing two adjacent base pairs, rather than "letters" for binding to single base pairs. Thus, polyamides designed to target Watson-Crick base pairs should utilize the strength of -ImPy- and -PyPy- central pairings. The f/Im and f/Py terminal groups yielded no advantage for their respective C/G or T/A base pairs. The exception is with the -ImPy- central pairing, for which f/Im has a 10-fold greater affinity for C/G than f/Py has for T/A.     

DNA duplexes and triplex-forming oligodeoxynucleotides incorporating modified nucleosides forming stable and selective triplexes

We have previously reported DNA triplexes containing the unnatural base triad G-PPI·C3, in which PPI is an indole-fused cytosine derivative incorporated into DNA duplexes and C3 is an abasic site in triplex-forming oligonucleotides (TFOs) introduced by a propylene linker. In this study, we developed a new unnatural base triad A-ψ·C(R1) where ψ and C(R1) are base moieties 2'-deoxypseudouridine and 5-substituted deoxycytidine, respectively. We examined several electron-withdrawing substituents for R1 and found that 5-bromocytosine (C(Br)) could selectively recognize ψ. In addition, we developed a new PPI derivative, PPI(Me), having a methyl group on the indole ring in order to achieve selective triplex formation between DNA duplexes incorporating various Watson-Crick base pairs, such as T-A, C-G, A-ψ, and G-PPI(Me), and TFOs containing T, C, C(Br), and C3. We studied the selective triplex formation between these duplexes and TFOs using UV-melting and gel mobility shift assays.