Understanding the role of base stacking in nucleic acids. MD and QM analysis of tandem GA base pairs in RNA duplexes

Investor logo
Investor logo

Warning

This publication doesn't include Faculty of Sports Studies. It includes Central European Institute of Technology. Official publication website can be found on muni.cz.
Authors

MORGADO Claudio A SVOZIL Daniel TURNER Douglas H ŠPONER Jiří

Year of publication 2012
Type Article in Periodical
Magazine / Source Physical Chemistry Chemical Physics
MU Faculty or unit

Central European Institute of Technology

Citation
Web http://www.ncbi.nlm.nih.gov/pubmed/22722325
Doi http://dx.doi.org/10.1039/c2cp40556c
Field Biophysics
Keywords DENSITY-FUNCTIONAL THEORY; GAUSSIAN-BASIS SETS; DOT-A MISMATCHES; MAGNETIC-RESONANCE-SPECTROSCOPY; RELAXATION MATRIX APPROACH; EMPIRICAL DISPERSION TERM; UNIQUE DINUCLEOTIDE STEPS; MOLECULAR-ORBITAL METHODS; LARGE SUBUNIT 23S-LIKE; PROTEIN DATA-BANK
Attached files
Description Preceding NMR experiments show that the conformation of tandem GA base pairs, an important recurrent non-canonical building block in RNA duplexes, is context dependent. The GA base pairs adopt "sheared" N3(G)-N6(A), N2(G)-N7(A) geometry in the r(CGAG)(2) and r(iGGAiC)(2) contexts while switching to "imino" N1(G)-N1(A), O6(G)-N6(A) geometry in the r(GGAC)(2) and r(iCGAiG)(2) contexts (iC and iG stand for isocytosine and isoguanine, respectively). As base stacking is likely to be one of the key sources of the context dependence of the conformation of GA base pairs, we calculated base stacking energies in duplexes containing such base pairs, to see if this dependence can be predicted by stacking energy calculations. When investigating the context dependence of the GA geometry two different conformations of the same duplex were compared (imino vs. sheared). The geometries were generated via explicit solvent MD simulations of the respective RNA duplexes, while the subsequent QM energy calculations focused on base stacking interactions of the four internal base pairs. Geometrical relaxation of nucleobase atoms prior to the stacking energy computations has a non-negligible effect on the results. The stacking energies were derived at the DFT-D/6-311++G(3df,3pd) level. We show a rather good correspondence between the intrinsic gas-phase stacking energies and the NMR-determined GA geometries. The conformation with more favorable gas-phase stacking is in most cases the one observed in experiments. This correlation is not improved when including solvent effects via the COSMO method. On the other side, the stacking calculations do not predict the relative thermodynamic stability of duplex formation for different sequences.
Related projects:

You are running an old browser version. We recommend updating your browser to its latest version.

More info