Computational model for predicting experimental RNA and DNA nearest-neighbor free energy rankings

Hydrogen-bonding, intrastrand base-stacking, and interstrand base-stacking energies were calculated for RNA and DNA dimers at the MP2(full)/6-311G** level of theory. Standard A-form RNA and B-form DNA geometries from average fiber diffraction data were employed for all base monomer and dimer geometries, and all dimer binding energies were obtained via single-point calculations. The effects of water solvation were considered using the PCM model. The resulting dimer binding energies were used to calculate the 10 unique RNA and 10 unique DNA computational nearest-neighbor energies, and the ranking of these computational nearest neighbor energies are in excellent agreement with the ranking of the experimental nearest-neighbor free energies. These results dispel the notion that average fiber diffraction geometries are insufficient for calculating RNA and DNA stacking energies.