NMR studies of DNA single strands and DNA:RNA hybrids with and without 1-propynylation at C5 of oligopyrimidines

The 1-propynylation at C5 of consecutive pyrimidines in DNA can enhance DNA:RNA hybrid stability at 37 degrees C by over 1 kcal/mol of substitution [Barnes, T. W., III; Turner, D. H. J. Am. Chem. Soc.2001, 123, 4107-4118]. To provide information on the structural consequences of propynylation, two-dimensional NMR spectroscopy was used to study the structures of several oligonucleotides. Intraresidue nuclear Overhauser effect spectroscopy cross peaks were observed at 30 degrees C and a 200 ms mixing time in the H6-H1' region for 5'(dC(P)C(P)U(P)C(P)C(P)U(P)U(P)) (ssPrODN) but not for 5'(dCCUCCUU) (ssODN), suggesting preorganization of the propynylated single strand. NMR structures of the duplexes 5'(dC(P)C(P)U(P)C(P)C(P)U(P)U(P))3':3'(rGAGGAGGAAAU)5' (PrODN:RNA), 5'(dCC(P)U(P)C(P)C(P)U(P)U(P))3':3'(rGAGGAGGAAAU)5' (sPrODN1:RNA), and 5'(dCCUCCUU)3':3'(rGAGGAGGAAAU)5' (ODN:RNA) indicate that their global structures are almost identical. The NMR data, however, suggest that the 5'-end of sPrODN1:RNA is more dynamic than that of PrODN:RNA. In the propynylated duplexes, the propyne group stacks on the aromatic ring of the 5'-base and extends into the major groove. The results suggest that the increased stability of the propynylated duplexes is caused by preorganization of the propynylated single strand and different interactions in the double strand. The propynyl group provides volume exclusion, enhanced stacking, and possibly different solvation.     

Stability and structure of RNA duplexes containing isoguanosine and isocytidine.

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.     

Remarkable stability of hairpins containing 2',5'-linked RNA loops.

We report here the results of a comparative study of hairpin loops that differ in the connectivity of phosphodiester linkages (3',5'- versus 2',5'-linkages). In addition, we have studied the effect of changing the stem composition on the thermodynamic stability of hairpin loops. Specifically, we constructed hairpins containing one of six stem duplex combinations, i.e., DNA:DNA ("DD"), RNA:RNA ("RR"), DNA:RNA ("DR"), 2',5'-RNA:RNA ("RR"), 2',5'-RNA:DNA ("RD"), and 2',5'-RNA:2',5'-RNA ("RR"), and one of three tetraloop compositions, i.e., 2',5'-RNA ("R"), RNA ("R"), and DNA ("D"). All hairpins contained the conserved and well-studied loop sequence 5'-...C(UUCG)G...-3' [Cheong et al. Nature 1990, 346, 680-682]. We show that the 2',5'-linked loop C(UUCG)G, i.e.,...C(3'p5')U(2'p5')U(2'p5')C(2'p5')G(2'p5')G(3'p5')..., like its "normal" RNA counterpart, forms an unusually stable tetraloop structure. We also show that the stability imparted by 2',5'-RNA loops is dependent on base sequence, a property that is shared with the regioisomeric 3',5'-RNA loops. Remarkably, we find that the stability of the UUCG tetraloop is virtually independent of the hairpin stem composition (DD, RR, RR, etc.), whereas the native RNA tetraloop exerts extra stability only when the stem is duplex RNA (R:R). As a result, the relative stabilities of hairpins with a 2',5'-linked tetraloop, e.g. ggac(UUCG)gtcc (T(m) = 61.4 degrees C), are often superior to those with RNA tetraloops, e.g. ggac(UUCG)gtcc (T(m) = 54.6 degrees C). In fact, it has been possible to observe the formation of a 2',5'-RNA:DNA hybrid duplex by linking the hybrid's strands to a (UUCG) loop. These duplexes (RD), which are not stable enough to form in an intermolecular complex [Wasner et al. Biochemistry 1998, 37, 7478-7486], were stable at room temperature (T(m) approximately 50 degrees C). Thus, 2',5'-loops have potentially important implications in the study of nucleic acid complexes where structural data are not yet available. Furthermore, they may be particularly useful as structural motifs for synthetic ribozymes and nucleic acid "aptamers".     

Interconversion of the cis-5R,6S- and trans-5R,6R-thymine glycol lesions in duplex DNA

Thymine glycol (Tg), 5,6-dihydroxy-5,6-dihydrothymine, is formed in DNA by the reaction of thymine with reactive oxygen species. The 5R Tg lesion was incorporated site-specifically into 5'-d(G(1)T(2)G(3)C(4)G(5)Tg(6)G(7)T(8)T(9)T(10)G(11)T(12))-3'; Tg = 5R Tg. The Tg-modified oligodeoxynucleotide was annealed with either 5'-d(A(13)C(14)A(15)A(16)A(17)C(18)A(19)C(20)G(21)C(22)A(23)C(24))-3', forming the Tg(6) x A(19) base pair, corresponding to the oxidative damage of thymine in DNA, or 5'-d(A(13)C(14)A(15)A(16)A(17)C(18)G(19)C(20)G(21)C(22)A(23)C(24))-3', forming the mismatched Tg(6) x G(19) base pair, corresponding to the formation of Tg following oxidative damage and deamination of 5-methylcytosine in DNA. At 30 degrees C, the equilibrium ratio of cis-5R,6S:trans-5R,6R epimers was 7:3 for the duplex containing the Tg(6) x A (19) base pair. In contrast, for the duplex containing the Tg(6) x G(19) base pair, the cis-5R,6S:trans-5R,6R equilibrium favored the cis-5R,6S epimer; the level of the trans-5R,6R epimer remained below the level of detection by NMR. The data suggested that Tg disrupted hydrogen bonding interactions, either when placed opposite to A(19) or G(19). Thermodynamic measurements indicated a 13 degrees C reduction of T(m) regardless of whether Tg was placed opposite dG or dA in the complementary strand. Although both pairings increased the free energy of melting by 3 kcal/mol, the melting of the Tg x G pair was more enthalpically favored than was the melting of the Tg x A pair. The observation that the position of the equilibrium between the cis-5R,6S and trans-5R,6R thymine glycol epimers in duplex DNA was affected by the identity of the complementary base extends upon observations that this equilibrium modulates the base excision repair of Tg [Ocampo-Hafalla, M. T.; Altamirano, A.; Basu, A. K.; Chan, M. K.; Ocampo, J. E.; Cummings, A., Jr.; Boorstein, R. J.; Cunningham, R. P.; Teebor, G. W. DNA Repair (Amst) 2006, 5, 444-454].     

Energetic contributions for the formation of TAT/TAT,TAT/CGC(+), and CGC(+)/CGC(+) base triplet stacks

We used a combination of spectroscopic and calorimetric techniques to determine complete thermodynamic profiles accompanying the folding of a set of triple helices and control duplexes. Specifically, we studied the sequences: d(A(7)C(5)T(7)C(5)T(7)), d(A(6)C(5)T(6)C(5)T(6)), d(A(6)C(5)T(6)), d(AGAGAGAC(5)TCTCTCTC(5)TCTCTCT), d(AGAGAC(5)TCTCTC(5)TCTCT), d(AGAGAC(5)TCTCTC(2)), d(AAGGAC(5)TCCTTC(5)TTCCT), d(AGGAAC(5)TTCCTC(5)TCCTT), and d(GAAAGC(5)CTTTCC(5)CTTTC). Circular dichroism spectroscopy indicated that all triplexes and duplexes are in the "B" conformation. DSC melting experiments revealed that the formation of triplexes is accompanied by a favorable free energy change, which arises from the compensation of a large and favorable enthalpic contribution with an unfavorable entropic contribution. Comparison of the thermodynamic profiles of these triplexes yielded enthalpic contributions of -24 kcal/mol, -23 kcal/mol, and -22 kcal/mol for the formation of TAT/TAT, TAT/CGC(+), and CGC(+)/CGC(+) base triplet stacks, respectively. UV melts as a function of sodium concentration show sodium ions stabilize the triplexes that contain only TAT triplets but destabilize the triplexes that contain CGC(+) triplets. UV melts as a function of pH indicate that the protonation of the third strand and loop cytosines stabilizes the triplexes that contain CGC(+) and TAT triplets, respectively. Our overall results suggest that the triplex to duplex transition of triplexes that contain CGC(+) triplets is accompanied by a release of protons and an uptake of sodium, while their duplex to random coil transition is accompanied by a release of sodium ions. A consequence of this opposite sodium dependence is that their coupled transitions are nearly independent of sodium concentration but are dependent on the experimental pH.     

Effects of fluorine on the structures and energetics of the propynyl and propargyl radicals and their anions

[reaction: see text] The adiabatic electron affinity (EA(ad)) of the CH(3)-C[triple bond]C(*) radical [experiment = 2.718 +/- 0.008 eV] and the gas-phase basicity of the CH(3)-C[triple bond]C:(-) anion [experiment = 373.4 +/- 2 kcal/mol] have been compared with those of their fluorine derivatives. The latter are studied using theoretical methods. It is found that there are large effects on the electron affinities and gas-phase basicities as the H atoms of the alpha-CH(3) group in the propynyl system are substituted by F atoms. The predicted electron affinities are 3.31 eV (FCH(2)-C[triple bond]C(*)), 3.86 eV (F(2)CH-C[triple bond]C(*)), and 4.24 eV (F(3)C-C[triple bond]C(*)), and the predicted gas-phase basicities of the fluorocarbanion derivatives are 366.4 kcal/mol (FCH(2)-C[triple bond]C:(-)), 356.6 kcal/mol (F(2)CH-C[triple bond]C:(-)), and 349.8 kcal/mol (F(3)C-C[triple bond]C:(-)). It is concluded that the electron affinities of fluoropropynyl radicals increase and the gas-phase basicities decrease as F atoms sequentially replace H atoms of the alpha-CH(3) in the propynyl system. The propargyl radicals, lower in energy than the isomeric propynyl radicals, are also examined and their electron affinities are predicted to be 0.98 eV ((*)CH(2)-C[triple bond]CH), 1.18 eV ((*)CFH-C[triple bond]CH), 1.32 eV ((*)CF(2)-C[triple bond] CH), 1.71 eV ((*)CH(2)-C[triple bond]CF), 2.05 eV ((*)CFH-C[triple bond]CF), and 2.23 eV ((*)CF(2)-C[triple bond]CF).     

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ～10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).     

Energetic preference of 8-oxoG eversion pathways in a DNA glycosylase

Base eversion is a fundamental process in the biochemistry of nucleic acids, allowing proteins engaged in DNA repair and epigenetic modifications to access target bases in DNA. Crystal structures reveal end points of these processes, but not the pathways involved in the dynamic process of base recognition. To elucidate the pathway taken by 8-oxoguanine during base excision repair by Fpg, we calculated free energy surfaces during eversion of the damaged base through the major and minor grooves. The minor groove pathway and free energy barrier (6-7 kcal/mol) are consistent with previously reported results (Qi, Y.; Spong, M. C.; Nam, K.; Banerjee, A.; Jiralerspong, S.; Karplus, M.; Verdine, G. L. Nature 2009, 462, 762.) However, eversion of 8-oxoG through the major groove encounters a significantly lower barrier (3-4 kcal/mol) more consistent with experimentally determined rates of enzymatic sliding during lesion search (Blainey, P. C.; van Oijent, A. M.; Banerjee, A.; Verdine, G. L.; Xie, X. S. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 5752.). Major groove eversion has been suggested for other glycosylases, suggesting that in addition to function, dynamics of base eversion may also be conserved.     

NMR solution structure of an oligonucleotide hairpin with a 2'F-ANA/RNA stem: implications for RNase H specificity toward DNA/RNA hybrid duplexes.

The first structure of a 2'-deoxy-2'-fluoro-D-arabinose nucleic acid (2'F-ANA)/RNA duplex is presented. We report the structural characterization by NMR spectroscopy of a small hybrid hairpin, r(GGAC)d(TTCG)2'F-a(GTCC), containing a 2'F-ANA/RNA stem and a four-residue DNA loop. Complete (1)H, (13)C, (19)F, and (31)P resonance assignments, scalar coupling constants, and NOE constraints were obtained from homonuclear and heteronuclear 2D spectra. In the chimeric duplex, the RNA strand adopts a classic A-form structure having C3' endo sugar puckers. The 2'F-ANA strand is neither A-form nor B-form and contains O4' endo sugar puckers. This contrasts strongly with the dynamic sugar conformations previously observed in the DNA strands of DNA/RNA hybrid duplexes. Structural parameters for the duplex, such as minor groove width, x-displacement, and inclination, were intermediate between those of A-form and B-form duplexes and similar to those of DNA/RNA duplexes. These results rationalize the enhanced stability of 2'F-ANA/RNA duplexes and their ability to elicit RNase H activity. The results are relevant for the design of new antisense drugs based on sugar-modified nucleic acids.     

Single-site mutation and secondary structure stability: an isodesmic reaction approach. The case of unnatural amino acid mutagenesis Ala-->Lac

A method is described to evaluate backbone interactions in proteins via computational unnatural amino acid mutagenesis. Several N-acetyl polyalanyl amides (AcA(n)NH(2)) were optimized in the representative helical (3(10)-, 4(13)-, and a "hybrid" kappa-helix, n = 7, 9, 10, 14) and hairpin (two- and three-stranded antiparallel beta-sheets with type I turns betaalphaalphaepsilon, n = 6, 9, 10) conformations, and extended conformers of N-acetyl polyalanyl methylamides (n = 2, 3) were used to derive multistranded beta-sheet fragments. Subsequently, each residue of every model structure was substituted, one at a time, with l-lactic acid. The resulting mutant structures were again optimized, and group-transfer energies DeltaE(GT) were obtained as heats of the isodesmic reactions: AcA(n)NHR + AcOMe --> AcA(x)LacA(y)NHR + AcNHMe (R = H, CH(3)). These group-transfer energies correlate with the degree of charge polarization of the substituted peptide linkages as measured by the difference Deltae in H and O Mulliken populations in HN-C=O and with the H-bond distances in the "wild-type" structures. A good correlation obtains for the HF/3-21G and B3LYP/6-31G* group-transfer energies. The destabilization effects are interpreted in terms of loss of interstrand and intrastrand H-bonds, decrease in Lewis basicity of the C=O group, and O...O repulsion. On the basis of several comparisons of Ala --> Lac DeltaE(GT)'s with heats of the NH --> CH(2) substitutions, the latter contribution is estimated (B3LYP/6-31G*) to range between 1.5 and 2.4 kcal mol(-1), a figure close to the recent experimental DeltaDeltaG(o) value of 2.6 kcal mol(-1) (McComas, C. C.; Crowley, B. M.; Boger, D. L. J. Am.Chem. Soc. 2003, 125, 9314). The partitioning yields the following maximum values of the electronic association energy of H-bonds in the examined sample of model structures (B3LYP/6-31G* estimates): 3(10)-helix D(e) = -1.7 kcal mol(-1), alpha-helix D(e) = -3.8 kcal mol(-1), beta-sheet D(e) = -6.1 kcal mol(-1). The premise of experimental evaluations of the backbone-backbone H-bonding that Ala --> Lac substitution in proteins is isosteric (e.g., Koh, J. T.; Cornish, V. W.; Schultz, P. G. Biochemistry 1997, 36, 11314) is often but not always corroborated. Examination of the integrity of H-bonding pattern and phi(i), psi(i) distribution identified several mutants with significant distortions of the "wild-type" structure resulting inter alia from the transitions between i, i + 3 and i, i + 4 H-bonding in helices, observed previously in the crystallographic studies of depsipeptides (Ohyama, T.; Oku, H.; Hiroki, A.; Maekawa, Y.; Yoshida, M.; Katakai, R. Biopolymers 2000, 54, 375; Karle, I. L.; Das, C.; Balaram, P. Biopolymers 2001, 59, 276). Thus, the isodesmic reaction approach provides a simple way to gauge how conformation of the polypeptide chain and dimensions of the H-bonding network affect the strength of backbone-backbone C=O...HN bonds. The results indicate that the stabilization provided by such interactions increases on going from 3(10)-helix to alpha-helix to beta-sheet.     

Structures and energetics of base flipping of the thymine dimer depend on DNA sequence

The cis,syn-cyclobutane pyrimidine dimer (CPD) is a photoinduced DNA lesion leading to a significant distortion of the DNA structure. Its repair by DNA photolyase requires a flip of the damaged base into an extrahelical position. This base flip is expected to be sequence-dependent, but the structures and energetics as a function of the bases 3' and 5' to the CPD lesion are unknown. Eight-nanosecond MD simulations of four different hexadecamer duplexes with the CPD were performed for the flipped-in and flipped-out structures. Analysis of these results indicates clear sequence-dependent differences. Significant disruptions of the base pairs to the 3' side of the CPD are observed for the flipped-out structures with adjacent A-T pairs, whereas those with G-C pairs adjacent show no such distortions. The conformational spaces occupied by these two duplexes are significantly different. The structural differences correlate well with the free energy differences for base flipping calculated using the previously established 2D potential of mean force (PMF) method. The energy differences for base flipping in duplexes containing A, T, G, and C pairs adjacent to the CPD were found to be 6.25-6.5, 5.25-5.5, 7.25-7.5, and 6.5-6.75 kcal/mol, respectively. These energy differences of up to 2 kcal/mol should be large enough to be detected experimentally using sensitive probes.     

Synthesis and pairing properties of oligodeoxynucleotides containing bicyclo-RNA and bicyclo-ANA modifications

The synthesis of the ribo(bc-rT)- and arabino(bc-araT)-version of bicyclothymidine (bc-dT) has been achieved. A conformational analysis by X-ray and/or (1)H NMR spectroscopy on the corresponding 3',5'-benzyl-protected nucleosides featured a rigid C(2')-endo conformation for the furanose ring, irrespective of the configuration of the OH group at C(2'). The conformation of the carbocyclic ring in these nucleosides was found to be less defined and thus more flexible. Both nucleosides were converted into the corresponding phosphoramidites and incorporated into oligodeoxynucleotides by standard DNA chemistry. T(m)-data of duplexes with cDNA and RNA revealed that a bc-rT unit strongly destabilized duplexes with cDNA and RNA by 6-8 °C/mod, while bc-araT was almost T(m) neutral. A rationale based on a previous structure of a bc-DNA mini duplex suggests that the strong destabilization caused by a bc-rT unit arises from unfavorable steric interactions of the equatorial 2'-OH group with the sugar residue of the 3'-neighboring nucleotide unit.     

A novel RNA motif based on the structure of unusually stable 2',5'-linked r(UUCG) loops

We have recently shown that hairpins containing 2',5'-linked RNA loops exhibit superior thermodynamic stability compared to native hairpins comprised of 3',5'-RNA loops [Hannoush, R. N.; Damha, M. J. J. Am. Chem. Soc. 2001, 123, 12368-12374]. A remarkable feature of the 2',5'-r(UUCG) tetraloop is that, unlike the corresponding 3',5'-linked tetraloop, its stability is virtually independent of the hairpin stem composition. Here, we determine the solution structure of unusually stable hairpins of the sequence 5'-G(1)G(2)A(3)C(4)-(U(5)U(6)C(7)G(8))-G(9)(U/T(10))C(11)C(12)-3' containing a 2',5'-linked RNA (UUCG) loop and either an RNA or a DNA stem. The 2',5'-linked RNA loop adopts a new fold that is completely different from that previously observed for the native 3',5'-linked RNA loop. The 2',5'-RNA loop is stabilized by (a). U5.G8 wobble base pairing, with both nucleotide residues in the anti-conformation, (b). extensive base stacking, and (c). sugar-base and sugar-sugar contacts, all of which contribute to the extra stability of this hairpin structure. The U5:G8 base pair stacks on top of the C4:G9 loop-closing base pair and thus appears as a continuation of the stem. The loop uracil U6 base stacks above U5 base, while the cytosine C7 base protrudes out into the solvent and does not participate in any of the stabilizing interactions. The different sugar pucker and intrinsic bonding interactions within the 2',5'-linked ribonucleotides help explain the unusual stability and conformational properties displayed by 2',5'-RNA tetraloops. These findings are relevant for the design of more effective RNA-based aptamers, ribozymes, and antisense agents and identify the 2',5'-RNA loop as a novel structural motif.     

Synthesis and properties of RNA analogues having amides as interuridine linkages at selected positions

Oligoribonucleotide analogues having amide internucleoside linkages (AM1: 3'-CH(2)CONH-5' and AM2: 3'-CH(2)NHCO-5') at selected positions have been synthesized and the thermal stability of duplexes formed by these analogues with complementary RNA fragments has been evaluated by UV melting experiments. Two series of oligomers with either 2'-OH or 2'-OMe vicinal to the amide linkages were studied. Monomeric synthons (3' and 5'-C amines and carboxylic acids) were synthesized as follows: For synthesis of the AM1 analogue, the known sequence of radical allylation followed by the cleavage of the double bond was adopted. For synthesis of the AM2 analogue, novel routes via addition of nitromethane followed by conversion of the nitro function to either amino or carboxyl groups were developed. Coupling of monomeric amines and carboxylic acids followed by protecting group manipulation and phosphonylation gave dimeric 3'-hydrogenphosphonate building blocks for oligonucleotide synthesis. Monomeric model compounds having 3'-amide and 2'-OH or 2'-OMe groups were also prepared and their conformational equilibrium was determined by (1)H NMR. The AM1 and AM2 models showed equal preferences for the North conformers (at 40 degrees C, 88-89% with 2'-OH, and 92-93% with 2'-OMe). At physiological salt concentration (0.1 M NaCl) the duplexes between AM1 modified oligonucleotides and RNA had stability similar to unmodified RNA-RNA duplexes (Delta t(m)= -0.2 to +0.7 degrees C per modification). However, the AM2 modification resulted in substantial stabilization of duplexes: Delta t(m)= +1 to +2.4 degrees C per modification compared to all RNA. A 2'-O-methyl vicinal to the AM2 linkage further increased the duplex stability. Our results suggest that RNA analogues having amide internucleoside bonds are very promising candidates for medicinal applications.     

Effects of electron attachment on C5'-O5' and C1'-N1 bond cleavages of pyrimidine nucleotides: A theoretical study

Sugar-base C(1')-N(1) and phosphate-sugar C(5')-O(5') bond breakings of 2'-deoxycytidine-5'-monophosphates (dCMP) and 2'-deoxythymidine-5'- monophosphates (dTMP) and their radical anions have been explored theoretically at the B3LYP/DZP++ level of theory. Calculations show that the low-energy electrons attachment to the pyrimidine nucleotides results in remarkable structural and chemical bonding changes. Predicted Gibbs free energies of reaction DeltaG for the C(5')-O(5') bond dissociation process of the radical anions are -14.6 and -11.5 kcal mol(-1), respectively, and such dissociation processes may be intrinsically spontaneous in the gas phase. Furthermore, the C(5')-O(5') bond cleavage processes of the anionic dCMP and dTMP were predicted to have activation energies of 6.9 and 8.0 kcal mol(-1) in the gas phase, respectively, much lower than the barriers for the C(1')-N(1) bond breaking process, showing that the C-O bond dissociation in DNA single strand breaks is a dominant process as observed experimentally.     

Synthesis of 2',4'-propylene-bridged (carba-ENA) thymidine and its analogues: the engineering of electrostatic and steric effects at the bottom of the minor groove for nuclease and thermodynamic stabilities and elicitation of RNase H

2',4'-Propylene-bridged thymidine (carba-ENA-T) and five 8'-Me/NH(2)/OH modified carba-ENA-T analogues have been prepared through intramolecular radical addition to C═N of the tethered oxime-ether. These carba-ENA nucleosides have been subsequently incorporated into 15mer oligodeoxynucleotides (AON), and their affinity toward cDNA and RNA, nuclease resistance, and RNase H recruitment capability have been investigated in comparison with those of the native and ENA counterparts. These carba-ENAs modified AONs are highly RNA-selective since all of them led to slight thermal stabilization effect for the AON:RNA duplex, but quite large destabilization effect for the AON:DNA duplex. It was found that different C8' substituents (at the bottom of the minor groove) on carba-ENA-T only led to rather small variation of thermal stability of the AON:RNA duplexes. We, however, observed that the parent carba-ENA-T modified AONs exhibited higher nucleolytic stability than those of the ENA-T modified counterparts. The nucleolytic stability of carba-ENA-T modified AONs can be further modulated by C8' substituent to variable extents depending on not only the chemical nature but also the stereochemical orientation of the C8' substituents: Thus, (1) 8'S-Me on carba-ENA increases the nucleolytic stability but 8'R-Me leads to a decreased effect; (2) 8'R-OH on carba-ENA had little, if any, effect on nuclease resistance but 8'S-OH resulted in significantly decreased nucleolytic stability; and (3) 8'-NH(2) substituted carba-ENA leads to obvious loss in the nuclease resistance. The RNA strand in all of the carba-ENA derivatives modified AON:RNA hybrid duplexes can be digested by RNase H1 with high efficiency, even at twice the rate of those of the native and ENA modified counterpart.     

Measurement of nucleobase pKa values in model mononucleotides shows RNA-RNA duplexes to be more stable than DNA-DNA duplexes

To understand why the RNA-RNA duplexes in general has a higher thermodynamic stability over the corresponding DNA-DNA duplexes, we have measured the pK(a) values of both nucleoside 3',5'-bis-ethyl phosphates [Etp(d/rN)pEt] and nucleoside 3'-ethyl phosphates [(d/rN)pEt] (N = A, G, C, or T/U), modeling as donors and acceptors of base pairs in duplexes. While the 3',5'-bis-phosphates, Etp(d/rN)pEt, mimic the internucleotidic monomeric units of DNA and RNA, in which the stacking contribution is completely absent, the 3'-ethyl phosphates, (d/rN)pEt, mimic the nucleotide at the 5'-end. The pK(a) values of the nucleobase in each of these model nucleoside phosphates have been determined with low pK(a) error (sigma = +/-0.01 to 0.02) by (1)H NMR (at 500 MHz) with 20-33 different pH measurements for each compound. This study has led us to show the following: (1) All monomeric DNA nucleobases are more basic than the corresponding RNA nucleobases in their respective Etp(d/rN)pEt and (d/rN)pEt. (2) The pK(a) values of the monomeric nucleotide blocks as well as Delta pK(a) values between the donor and acceptor can be used to understand the relative base-pairing strength in the oligomeric duplexes in the RNA and DNA series. (3) The Delta G*(pKa) of the donor and acceptor of the base pair in duplexes enables a qualitative dissection of the relative strength of the base-pairing and stacking in the RNA-RNA over the DNA-DNA duplexes. (4) It is also found that the relative contribution of base-pairing strength and nucleobase stacking in RNA-RNA over DNA-DNA is mutually compensating as the % A-T/U content increases or decreases. This interdependency of stacking and hydrogen bonding can be potentially important in the molecular design of the base-pair mimics to expand the alphabet of the genetic code.     

Short-RNA selective binding of oligonucleotides modified using adenosine and guanosine derivatives that possess cyclohexyl phosphates as substituents

We have developed new artificial oligonucleotides which distinguish short RNA targets from long ones. The modification of the 5' termini of oligonucleotides by using adenosine derivatives that possess a bulky cyclohexyl phosphate moiety at their base moiety and a phosphate group at the position of their 5'-hydroxyl group maximized their short RNA selectivity. The 2'-O-methyl-RNA (5'-XC(m)A(m)A(m)C(m)C(m)U(m)A(m)C(m)U(m)) having these modifications exhibits ca. 10 °C higher T(m) in the duplexes with the complementary short RNA (3'-GUUGGAUGA-5') than with the long RNA (3'-AUUAUAUGUUGGAUGAUGGUUA-5'). The oligodeoxynucleotides having the same modification exhibited similar selectivity. Such short-RNA selective binding of terminally modified oligonucleotides can be employed to distinguish between mature microRNAs and pre-microRNAs.     

The possibility of the decomposition of 2'-deoxyribose moiety of thymidine induced by the low energy electron attachment

To evaluate the possibility of the decomposition of 2-deoxyribose moiety of thymidine induced by low energy electrons (LEE) attachment, the transition states and the energy barriers of the bond breaking processes of the ribose of the nucleoside have been studied theoretically by applying the density functional theory with the double zeta basis sets (DZP++). The energy barriers for the breakage of the C-C bonds (C(1')-C(2'), C(2')-C(3'), C(3')-C(4'), and C(4')-C(5')) of the ribose group of the radical anion of thymidine are found to be high (ca. 42-57 kcal/mol). The total energies of the C-C bond-broken products are significantly higher than that of the radical anion dT(*-). The decomposition of dT(*-) through the C-C bond rupture is unlikely to take place. The rupture of the C(1')-O(4') bond of dT(*-) needs an activation energy as low as 10.4 kcal/mol. However, the reversed reaction (C(1')-O(4') bond formation) needs the activation energy low as 0.3 kcal/mol. Therefore, the intermediate product LM1(C1')-(O4') is unlikely to be stable and the C(1')-O(4') bond-broken is not favored. The activation energy of the C(4')-O(4') bond rupture process amounts to 20.5 kcal/mol. The total energy of the C(4')-O(4') bond broken product is about 6.5 kcal/mol lower than that of the reactant dT(*-). The subsequent N1-glycosidic bond breaking process is found to have a very low energy barrier. Therefore, the LEE-induced base release through the C(4')-O(4') bond rupture might be a possible pathway.     

Architecture of a validated microRNA::target interaction

MicroRNAs are small approximately 22 nucleotide regulators of numerous biological processes and bind target gene messenger RNAs to control gene expression. The C. elegans microRNA let-7 and its target lin-41 were the first microRNA::target interaction to be validated in vivo. let-7 molecules form imperfect duplexes with two required let-7 complementary sites in the lin-41 3' UTR. Here, we show that base pairing at both the 5' and 3' ends of the let-7 binding site, as well as the presence of unpaired RNA residues in the predicted duplexes, are required for lin-41 downregulation. In this study, our model for microRNA::target interactions also demonstrates that the context of a microRNA binding can be critical for function, revealing an unforeseen complexity in microRNA::target interactions.