The base sequence dependent flexibility of linear single-stranded oligoribonucleotides correlates with the reactivity of the phosphodiester bond

The effect of base sequence on the structure and flexibility of linear single-stranded RNA molecules and the influence of the base sequence on phosphodiester bond reactivity have been studied. Molecular dynamics simulations of 2.1 ns were carried out for nine chimeric oligonucleotides containing only one unsubstituted ribo unit, all the rest of sugars being 2'-O-methylated. The base sequence has recently been reported to make a big contribution to the reactivity of these compounds. A detailed examination of the interaction energies between the base moieties shows that base stacking is strongly context-dependent and cooperative. The strength of stacking at the site susceptible to chain cleavage by intramolecular transesterification was observed to be dependent on both the flanking bases of the cleavage site and those further apart in the molecule. The interaction energies between the bases in the vicinity of the scissile linkage were found to correlate well with the experimental phosphodiester bond cleavage rates: the stronger the bases close to the cleavage site are stacked, the slower the cleavage rate is.