Insights into Stability and Folding of GNRA and UNCG Tetra loops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics

HALDAR, Susanta, Petra KUHROVÁ, Pavel BANÁŠ, Vojtěch SPIWOK, Jiří ŠPONER, Pavel HOBZA a Michal OTYEPKA. Insights into Stability and Folding of GNRA and UNCG Tetra loops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics. Journal of Chemical Theory and Computation. WASHINGTON: AMER CHEMICAL SOC, roč. 11, č. 8, s. 3866-3877. ISSN 1549-9618. doi:10.1021/acs.jctc.5b00010. 2015.

Základní údaje

• Originální název: Insights into Stability and Folding of GNRA and UNCG Tetra loops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics
• Autoři: HALDAR, Susanta, Petra KUHROVÁ, Pavel BANÁŠ, Vojtěch SPIWOK, Jiří ŠPONER, Pavel HOBZA a Michal OTYEPKA.
• Vydání: Journal of Chemical Theory and Computation, WASHINGTON, AMER CHEMICAL SOC, 2015, 1549-9618.

Další údaje

• Originální jazyk: angličtina
• Typ výsledku: Článek v odborném periodiku
• Obor: 10403 Physical chemistry
• Stát vydavatele: Spojené státy
• Utajení: není předmětem státního či obchodního tajemství
• WWW: URL
• Impakt faktor: Impact factor: 5.301
• Organizační jednotka: Středoevropský technologický institut
• Doi: http://dx.doi.org/10.1021/acs.jctc.5b00010
• UT WoS: 000359500000034
• Klíčová slova anglicky: FREE-ENERGY LANDSCAPE; DI-GMP RIBOSWITCH; RNA HAIRPINS; DISTINCTIVE FEATURES; ANGSTROM RESOLUTION; HAMMERHEAD RIBOZYME; BACTERIAL RIBOSOME; REPLICA-EXCHANGE; SIMULATIONS; KINETICS
• Štítky: neMU
• Příznaky: Mezinárodní význam, Recenzováno

Anotace

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15,us-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the G(L4) nudeobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that G(L4)-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.