Computational evaluation of the stability of 2′-O-methyl-RNA/RNA duplexes incorporating 3-deazaguanine derivatives by ab initio calculations and a molecular dynamics simulation

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.

Synthesis and hybridization properties of 2'-O-methylated oligoribonucleotides incorporating 2'-O-naphthyluridines

2'-O-(1-Naphthyl)uridine and 2'-O-(2-naphthyl)uridine were synthesized by a microwave-mediated reaction of 2,2'-anhydrouridine with naphthols. Using the 3'-phosphoramidite building blocks, these 2'-O-aryluridine derivatives were incorporated into 2'-O-methylated oligoribonucleotides. Incorporation of five 2'-O-(2-naphthyl)uridines into a 2'-O-methylated RNA sense strand significantly increased the thermostability of the duplex with a 2'-O-methylated RNA antisense strand. Circular dichroism spectroscopy and molecular dynamic simulation of the duplexes formed between the modified RNAs and 2'-O-methyl RNAs suggested that there are π-π interactions between two neighboring naphthyl groups in a sequence of the five consecutively modified nucleosides.