This Article Full Text Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Villescas-Diaz, G. Articles by Zacharias, M. Search for Related Content PubMed PubMed Citation Articles by Villescas-Diaz, G. Articles by Zacharias, M.

Sequence Context Dependence of Tandem Guanine:Adenine Mismatch Conformations in RNA: A Continuum Solvent Analysis

Theoretische Biophysik, Institut für Molekulare Biotechnologie, Jena, Germany; and ^* School of Engineering and Science, International University Bremen, Germany

Correspondence: Address reprint requests to Martin Zacharias, School of Engineering and Science, International University Bremen, Campus Ring 1, D-28759 Bremen, Germany. Tel.: 49-421-200-3541; Fax: 49-421-200-3249; E-mail: m.zacharias{at}iu-bremen.de.

Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.

This article has been cited by other articles:

				K. Reblova, E. Fadrna, J. Sarzynska, T. Kulinski, P. Kulhanek, E. Ennifar, J. Koca, and J. Sponer Conformations of Flanking Bases in HIV-1 RNA DIS Kissing Complexes Studied by Molecular Dynamics Biophys. J., December 1, 2007; 93(11): 3932 - 3949. [Abstract] [Full Text] [PDF]

				K. Seio, T. Sasami, R. Tawarada, and M. Sekine Synthesis of 2'-O-methyl-RNAs incorporating a 3-deazaguanine, and UV melting and computational studies on its hybridization properties Nucleic Acids Res., September 11, 2006; 34(16): 4324 - 4334. [Abstract] [Full Text] [PDF]