help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
This Article
Right arrow Full Text
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sarzynska, J.
Right arrow Articles by Kulinski, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sarzynska, J.
Right arrow Articles by Kulinski, T.
Biophysical Journal 85:3445-3459 (2003)
© 2003 The Biophysical Society

Effects of Base Substitutions in an RNA Hairpin from Molecular Dynamics and Free Energy Simulations

Joanna Sarzynska *, Lennart Nilsson {dagger} and Tadeusz Kulinski *

* Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan, Poland; and {dagger} Center for Structural Biochemistry, Department of Biosciences at NOVUM, Karolinska Institutet, Huddinge, Sweden

Correspondence: Address reprint requests to Tadeusz Kulinski, E-mail: tadkul{at}ibch.poznan.pl.

Contributions of individual interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using base substitutions. The G in the first tetraloop position was replaced by inosine (I) or adenosine (A), and the G in the C-G basepair closing the tetraloop was replaced by I. These substitutions eliminate particular hydrogen bonds proposed in the nuclear magnetic resonance model of the GCAA tetraloop. Molecular dynamics simulations of the GCAA tetraloop in aqueous solvent displayed a well-defined hydrogen pattern between the first and last loop nucleotides (G and A) stabilized by a bridging water molecule. Substitution of G->I in the basepair closing the tetraloop did not significantly influence the loop structure and dynamics. The ICAA loop maintained the overall structure, but displayed variation in the hydrogen-bond network within the tetraloop itself. Molecular dynamics simulations of the ACAA loop led to conformational heterogeneity of the resulting structures. Changes of hairpin formation free energy associated with substitutions of individual bases were calculated by the free energy perturbation method. The calculated decrease of the hairpin stability upon G->I substitution in the C-G basepair closing the tetraloop was in good agreement with experimental thermodynamic data. Our theoretical estimates for G->I and G->A mutations located in the tetraloop suggest larger loop destabilization than corresponding experimental results. The extent of conformational sampling of the structures resulting from base substitutions and its impact on the calculated free energy was discussed.




This article has been cited by other articles:


Home page
 J. Gen. Virol.Home page
J. R. Thompson, E. Buratti, M. de Wispelaere, and M. Tepfer
Structural and functional characterization of the 5' region of subgenomic RNA5 of cucumber mosaic virus
J. Gen. Virol., July 1, 2008; 89(7): 1729 - 1738.
[Abstract] [Full Text] [PDF]

Home page
 Biophys. JHome page
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]

Home page
 Nucleic Acids ResHome page
N. Spackova and J. Sponer
Molecular dynamics simulations of sarcin-ricin rRNA motif
Nucleic Acids Res., February 2, 2006; 34(2): 697 - 708.
[Abstract] [Full Text] [PDF]


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2003 by the Biophysical Society.