Using Density Functional Theory To Study Neutral and Ionized Stacked Thymine Dimers

Influence of the Number, Nature and Position of Methyl Posttranscriptional Modifications on Nucleobase Stacking in RNA

DFT calculations are employed to quantify the influence of the presence, number, nature, and position of posttranscriptional methylation on stacking strength of RNA bases. We carry out detailed potential energy scans of the variation in stacking energies with characteristic geometrical parameters in three categories of forty stacked dimers - canonical base homodimers (N||N), methylated base homodimers (mN||mN) and heterodimers of canonical bases and methylated counterparts (N||mN). Our analysis reveals that neutral methylation invariably enhances the stacking of bases. Further, N||mN stacking is stronger than mN||mN stacking and charged N||mN exhibit strongest stacking among all dimers. This indicates that methylations greatly enhance stacking when dispersed in RNA sequences containing identical bases. Comparison of stacks involving singly- and doubly-methylated purines reveal that incremental methylation enhances the stacking in neutral dimers. Although methylation at the carbon position of neutral pyrimidine dimers greatly enhances the stacking, methylation on the 5-membered ring imparts better stacking compared to methylation on the 6-membered ring in adenine dimers. However, methylation at the ring nitrogen (N¹ ) provides better stacking than the amino group (N² ) in guanine dimers. Our results thus highlight subtle structural effects of methylation on RNA base stacking and will enhance our understanding of the physicochemical principles of RNA structure and dynamics.