1
|
Ormazábal A, Palma J, Pierdominici-Sottile G. Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods. Methods Mol Biol 2024; 2741:207-238. [PMID: 38217656 DOI: 10.1007/978-1-0716-3565-0_12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2024]
Abstract
Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.
Collapse
Affiliation(s)
- Agustín Ormazábal
- Departmento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Bernal, Buenos Aires, Argentina
- Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, Godoy Cruz, CABA, Argentina
| | - Juliana Palma
- Departmento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Bernal, Buenos Aires, Argentina
- Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, Godoy Cruz, CABA, Argentina
| | - Gustavo Pierdominici-Sottile
- Departmento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Bernal, Buenos Aires, Argentina.
- Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, Godoy Cruz, CABA, Argentina.
| |
Collapse
|
2
|
Yamamoto S, Kuwahara R, Tanaka K. Dynamic behaviour of water molecules in heterogeneous free space formed in an epoxy resin. SOFT MATTER 2021; 17:6073-6080. [PMID: 34132297 DOI: 10.1039/d1sm00529d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Although an epoxy resin is a stable material, it absorbs moisture over a long period of time, causing deterioration of its material properties. We here applied a full-atomistic molecular dynamics (MD) simulation to study where water molecules exist in an epoxy resin and how they dynamically behave. First, the curing reaction was simulated to obtain a network structure so that the time course of the density, and thereby the free space, in the resin were obtained. The results made it possible to discuss the formation and size distribution of the free spaces which were not connected to each other. Then, a few percent of water were inserted into the free space of the cured epoxy resin to examine the location and dynamics of their molecules. We found that several water molecules were clustered at a preferred site, where hydrogen bonds can be formed with hydroxy, ether and amino groups of the network, in the free space, and they heterogeneously moved from there to other sites.
Collapse
Affiliation(s)
- Satoru Yamamoto
- Centre for Polymer Interface and Molecular Adhesion Science, Kyushu University, Fukuoka 819-0395, Japan.
| | | | | |
Collapse
|
3
|
Hurst T, Chen SJ. Deciphering nucleotide modification-induced structure and stability changes. RNA Biol 2021; 18:1920-1930. [PMID: 33586616 DOI: 10.1080/15476286.2021.1882179] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Nucleotide modification in RNA controls a bevy of biological processes, including RNA degradation, gene expression, and gene editing. In turn, misregulation of modified nucleotides is associated with a host of chronic diseases and disorders. However, the molecular mechanisms driving these processes remain poorly understood. To partially address this knowledge gap, we used alchemical and temperature replica exchange molecular dynamics (TREMD) simulations on an RNA duplex and an analogous hairpin to probe the structural effects of modified and/or mutant nucleotides. The simulations successfully predict the modification/mutation-induced relative free energy change for complementary duplex formation, and structural analyses highlight mechanisms driving stability changes. Furthermore, TREMD simulations for a hairpin-forming RNA with and without modification provide reliable estimations of the energy landscape. Illuminating the impact of methylated and/or mutated nucleotides on the structure-function relationship and the folding energy landscape, the simulations provide insights into modification-induced alterations to the folding mechanics of the hairpin. The results here may be biologically significant as hairpins are widespread structure motifs that play critical roles in gene expression and regulation. Specifically, the tetraloop of the probed hairpin is phylogenetically abundant, and the stem mirrors a miRNA seed region whose modification has been implicated in epilepsy pathogenesis.
Collapse
Affiliation(s)
- Travis Hurst
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, MO, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, MO, USA
| |
Collapse
|
4
|
Leonarski F, Jasiński M, Trylska J. Thermodynamics of the fourU RNA thermal switch derived from molecular dynamics simulations and spectroscopic techniques. Biochimie 2018; 156:22-32. [PMID: 30244136 DOI: 10.1016/j.biochi.2018.09.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2018] [Accepted: 09/10/2018] [Indexed: 10/28/2022]
Abstract
Bacterial 5' untranslated regions of mRNA, termed thermal switches or thermometers, change their structure in response to temperature change. This structural change provides for the regulation of gene expression. One of such thermal switches, called fourU, is present in the Salmonella species. Mutations of fourU were found to abrogate its regulatory properties. We investigated the thermodynamics of the fourU fragment responsible for its structural changes. All-atom molecular dynamics simulations at various temperatures and spectroscopic experiments in solution were performed for the wild-type fourU and its mutants. We found that the U11C and A8C mutations stabilize the fourU structure in comparison to the wild-type fourU, and the double-point G14A/C25U mutant has the most destabilizing effect on the fourU hairpin 2 responsible for temperature sensing. The G14A/C25U mutant is also the easiest to strand-invade by a complementary oligonucleotide as indicated by fluorescence spectroscopy experiments.
Collapse
Affiliation(s)
- Filip Leonarski
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland
| | - Maciej Jasiński
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland
| | - Joanna Trylska
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097, Warsaw, Poland.
| |
Collapse
|
5
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 327] [Impact Index Per Article: 54.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
6
|
RNA structure refinement using NMR solvent accessibility data. Sci Rep 2017; 7:5393. [PMID: 28710477 PMCID: PMC5511288 DOI: 10.1038/s41598-017-05821-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2017] [Accepted: 06/02/2017] [Indexed: 12/21/2022] Open
Abstract
NMR spectroscopy is a powerful technique to study ribonucleic acids (RNAs) which are key players in a plethora of cellular processes. Although the NMR toolbox for structural studies of RNAs expanded during the last decades, they often remain challenging. Here, we show that solvent paramagnetic relaxation enhancements (sPRE) induced by the soluble, paramagnetic compound Gd(DTPA-BMA) provide a quantitative measure for RNA solvent accessibility and encode distance-to-surface information that correlates well with RNA structure and improves accuracy and convergence of RNA structure determination. Moreover, we show that sPRE data can be easily obtained for RNAs with any isotope labeling scheme and is advantageous regarding sample preparation, stability and recovery. sPRE data show a large dynamic range and reflect the global fold of the RNA suggesting that they are well suited to identify interaction surfaces, to score structural models and as restraints in RNA structure determination.
Collapse
|
7
|
Pathak AK, Bandyopadhyay T. Water isotope effect on the thermostability of a polio viral RNA hairpin: A metadynamics study. J Chem Phys 2017; 146:165104. [DOI: 10.1063/1.4982049] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Arup K. Pathak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai 400 085, India
| | - Tusar Bandyopadhyay
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai 400 085, India
| |
Collapse
|
8
|
Borkar AN, Vallurupalli P, Camilloni C, Kay LE, Vendruscolo M. Simultaneous NMR characterisation of multiple minima in the free energy landscape of an RNA UUCG tetraloop. Phys Chem Chem Phys 2017; 19:2797-2804. [PMID: 28067358 PMCID: PMC6529357 DOI: 10.1039/c6cp08313g] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
RNA molecules in solution tend to undergo structural fluctuations of relatively large amplitude and to populate a range of different conformations some of which with low populations. It is still very challenging, however, to characterise the structures of these low populated states and to understand their functional roles. In the present study, we address this problem by using NMR residual dipolar couplings (RDCs) as structural restraints in replica-averaged metadynamics (RAM) simulations. By applying this approach to a 14-mer RNA hairpin containing the prototypical UUCG tetraloop motif, we show that it is possible to construct the free energy landscape of this RNA molecule. This free energy landscapes reveals the surprisingly rich dynamics of the UUCG tetraloop and identifies the multiple substates that exist in equilibrium owing to thermal fluctuations. The approach that we present is general and can be applied to the study of the free energy landscapes of other RNA or RNA-protein systems.
Collapse
Affiliation(s)
- Aditi N Borkar
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, UK.
| | - Pramodh Vallurupalli
- Departments of Molecular Genetics, Biochemistry, and Chemistry, University of Toronto, Toronto, Canada M5S 1A8
| | - Carlo Camilloni
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, UK.
| | - Lewis E Kay
- Departments of Molecular Genetics, Biochemistry, and Chemistry, University of Toronto, Toronto, Canada M5S 1A8
| | | |
Collapse
|
9
|
Šponer J, Krepl M, Banáš P, Kührová P, Zgarbová M, Jurečka P, Havrila M, Otyepka M. How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes? WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27863061 DOI: 10.1002/wrna.1405] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 01/01/2023]
Abstract
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| |
Collapse
|
10
|
Prostova MA, Gmyl AP, Bakhmutov DV, Shishova AA, Khitrina EV, Kolesnikova MS, Serebryakova MV, Isaeva OV, Agol VI. Mutational robustness and resilience of a replicative cis-element of RNA virus: Promiscuity, limitations, relevance. RNA Biol 2016; 12:1338-54. [PMID: 26488412 DOI: 10.1080/15476286.2015.1100794] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Since replication of RNA-viruses is generally a low-fidelity process, it would be advantageous, if specific interactions of their genomic cis-elements with dedicated ligands are relatively tolerant to mutations. The specificity/promiscuity trade-off of such interactions was addressed here by investigating structural requirements of the oriL (also known as the clover leaf-like element), of poliovirus RNA, a replicative cis-element containing a conserved essential tetraloop functionally interacting with the viral protein 3CD. The sequence of this tetraloop and 2 adjacent base-pairs was randomized in the viral genome, and viable viruses were selected in susceptible cells. Strikingly, each position of this octanucleotide in 62 investigated viable viruses could be occupied by any nucleotide (with the exception of one position, which lacked U), though with certain sequence preferences, confirmed by engineering mutant viral genomes whose phenotypic properties were found to correlate with the strength of the cis-element/ligand interaction. The results were compatible with a hypothesis that functional recognition by 3CD requires that this tetraloop should stably or temporarily adopt a YNMG-like (Y=U/C, N=any nucleotide, M=A/C) fold. The fitness of "weak" viruses could be increased by compensatory mutations "improving" the tetraloops. Otherwise, the recognition of "bad" tetraloops might be facilitated by alterations in the 3CD protein. The virus appeared to tolerate mutations in its cis-element relaying on either robustness (spatial structure degeneracy) or resilience (a combination of dynamic RNA folding, low-fidelity replication modifying the cis-element or its ligand, and negative selection). These mechanisms (especially resilience involving metastable low-fit intermediates) can also contribute to the viral evolvability.
Collapse
Affiliation(s)
- Maria A Prostova
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia
| | - Anatoly P Gmyl
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia.,b M V Lomonosov Moscow State University ; Moscow Russia
| | - Denis V Bakhmutov
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia.,c Deceased
| | - Anna A Shishova
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia
| | - Elena V Khitrina
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia
| | - Marina S Kolesnikova
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia
| | | | - Olga V Isaeva
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia
| | - Vadim I Agol
- a M P Chumakov Institute of Poliomyelitis and Viral Encephalitides ; Moscow Russia.,b M V Lomonosov Moscow State University ; Moscow Russia
| |
Collapse
|
11
|
Denschlag R, Lingenheil M, Tavan P, Mathias G. Simulated Solute Tempering. J Chem Theory Comput 2015; 5:2847-57. [PMID: 26631796 DOI: 10.1021/ct900274n] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For the enhanced conformational sampling in molecular dynamics (MD) simulations, we present "simulated solute tempering" (SST) which is an easy to implement variant of simulated tempering. SST extends conventional simulated tempering (CST) by key concepts of "replica exchange with solute tempering" (REST, Liu et al. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 13749). We have applied SST, CST, and REST to molecular dynamics (MD) simulations of an alanine octapeptide in explicit water. The weight parameters required for CST and SST are determined by two different formulas whose performance is compared. For SST only one of them yields a uniform sampling of the temperature space. Compared to CST and REST, SST provides the highest exchange probabilities between neighboring rungs in the temperature ladder. Concomitantly, SST leads to the fastest diffusion of the simulation system through the temperature space, in particular, if the "even-odd" exchange scheme is employed in SST. As a result, SST exhibits the highest sampling speed of the investigated tempering methods.
Collapse
Affiliation(s)
- Robert Denschlag
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität, Oettingenstrasse 67, 80538 München, Germany
| | - Martin Lingenheil
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität, Oettingenstrasse 67, 80538 München, Germany
| | - Paul Tavan
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität, Oettingenstrasse 67, 80538 München, Germany
| | - Gerald Mathias
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität, Oettingenstrasse 67, 80538 München, Germany
| |
Collapse
|
12
|
Chakraborty D, Collepardo-Guevara R, Wales DJ. Energy Landscapes, Folding Mechanisms, and Kinetics of RNA Tetraloop Hairpins. J Am Chem Soc 2014; 136:18052-61. [DOI: 10.1021/ja5100756] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Debayan Chakraborty
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | | | - David J. Wales
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
13
|
Rastädter D, Biswas M, Burghardt I. Molecular Dynamics Study of the Controlled Destabilization of an RNA Hairpin Structure by a Covalently Attached Azobenzene Switch. J Phys Chem B 2014; 118:8478-88. [DOI: 10.1021/jp501399k] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Dominique Rastädter
- Institute
of Physical and
Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Strasse
7, 60438 Frankfurt
am Main, Germany
| | - Mithun Biswas
- Institute
of Physical and
Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Strasse
7, 60438 Frankfurt
am Main, Germany
| | - Irene Burghardt
- Institute
of Physical and
Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Strasse
7, 60438 Frankfurt
am Main, Germany
| |
Collapse
|
14
|
High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations. Proc Natl Acad Sci U S A 2013; 110:16820-5. [PMID: 24043821 DOI: 10.1073/pnas.1309392110] [Citation(s) in RCA: 221] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We report the de novo folding of three hyperstable RNA tetraloops to 1-3 Å rmsd from their experimentally determined structures using molecular dynamics simulations initialized in the unfolded state. RNA tetraloops with loop sequences UUCG, GCAA, or CUUG are hyperstable because of the formation of noncanonical loop-stabilizing interactions, and they are all faithfully reproduced to angstrom-level accuracy in replica exchange molecular dynamics simulations, including explicit solvent and ion molecules. This accuracy is accomplished using unique RNA parameters, in which biases that favor rigid, highly stacked conformations are corrected to accurately capture the inherent flexibility of ssRNA loops, accurate base stacking energetics, and purine syn-anti interconversions. In a departure from traditional quantum chemistrycentric approaches to force field optimization, our parameters are calibrated directly from thermodynamic and kinetic measurements of intra- and internucleotide structural transitions. The ability to recapitulate the signature noncanonical interactions of the three most abundant hyperstable stem loop motifs represents a significant milestone to the accurate prediction of RNA tertiary structure using unbiased all-atom molecular dynamics simulations.
Collapse
|
15
|
Halder S, Bhattacharyya D. RNA structure and dynamics: a base pairing perspective. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2013; 113:264-83. [PMID: 23891726 DOI: 10.1016/j.pbiomolbio.2013.07.003] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/16/2013] [Revised: 06/25/2013] [Accepted: 07/16/2013] [Indexed: 12/12/2022]
Abstract
RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.
Collapse
Affiliation(s)
- Sukanya Halder
- Biophysics division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata 700 064, India
| | | |
Collapse
|
16
|
Propensities for loop structures of RNA & DNA backbones. Biophys Chem 2013; 180-181:110-8. [PMID: 23933331 DOI: 10.1016/j.bpc.2013.07.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Revised: 07/10/2013] [Accepted: 07/10/2013] [Indexed: 11/21/2022]
Abstract
RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2'-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2'-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2'-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3'-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ-dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ-dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2'-hydroxyls, and consequently renders loop structures in RNA more stable.
Collapse
|
17
|
Kara M, Zacharias M. Theoretical studies of nucleic acids folding. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1146] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Mahmut Kara
- Physics Department T38, Technical University Munich, Garching, Germany
| | - Martin Zacharias
- Martin Zacharias, Physics Department T38, Technical University Munich, Garching, Germany
| |
Collapse
|
18
|
Henriksen NM, Roe DR, Cheatham TE. Reliable oligonucleotide conformational ensemble generation in explicit solvent for force field assessment using reservoir replica exchange molecular dynamics simulations. J Phys Chem B 2013; 117:4014-27. [PMID: 23477537 PMCID: PMC3775460 DOI: 10.1021/jp400530e] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Molecular dynamics force field development and assessment requires a reliable means for obtaining a well-converged conformational ensemble of a molecule in both a time-efficient and cost-effective manner. This remains a challenge for RNA because its rugged energy landscape results in slow conformational sampling and accurate results typically require explicit solvent which increases computational cost. To address this, we performed both traditional and modified replica exchange molecular dynamics simulations on a test system (alanine dipeptide) and an RNA tetramer known to populate A-form-like conformations in solution (single-stranded rGACC). A key focus is on providing the means to demonstrate that convergence is obtained, for example, by investigating replica RMSD profiles and/or detailed ensemble analysis through clustering. We found that traditional replica exchange simulations still require prohibitive time and resource expenditures, even when using GPU accelerated hardware, and our results are not well converged even at 2 μs of simulation time per replica. In contrast, a modified version of replica exchange, reservoir replica exchange in explicit solvent, showed much better convergence and proved to be both a cost-effective and reliable alternative to the traditional approach. We expect this method will be attractive for future research that requires quantitative conformational analysis from explicitly solvated simulations.
Collapse
Affiliation(s)
- Niel M. Henriksen
- Department of Medicinal Chemistry, College of Pharmacy, 2000 East 30 South Skaggs 201, University of Utah, Salt Lake City, UT, 84112, USA
| | - Daniel R. Roe
- Department of Medicinal Chemistry, College of Pharmacy, 2000 East 30 South Skaggs 201, University of Utah, Salt Lake City, UT, 84112, USA
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, 2000 East 30 South Skaggs 201, University of Utah, Salt Lake City, UT, 84112, USA
| |
Collapse
|
19
|
Kührová P, Banáš P, Best RB, Šponer J, Otyepka M. Computer Folding of RNA Tetraloops? Are We There Yet? J Chem Theory Comput 2013; 9:2115-25. [DOI: 10.1021/ct301086z] [Citation(s) in RCA: 76] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Petra Kührová
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Robert B. Best
- Laboratory of Chemical Physics,
National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520,
United States
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC − Central European
Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies
and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| |
Collapse
|
20
|
Paladino A, Zangi R. Ribose 2'-Hydroxyl Groups Stabilize RNA Hairpin Structures Containing GCUAA Pentaloop. J Chem Theory Comput 2013; 9:1214-21. [PMID: 26588764 DOI: 10.1021/ct3006216] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The chemical structure of RNA and DNA is very similar; however, the three-dimensional conformation of these two nucleic acids is very different. Whereas the DNA adopts a repetitive structure of a double-stranded helix, RNA is primarily single stranded with a complex three-dimensional structure in which the hairpin is the most common secondary structure. Apart from the difference between uracil and thymine, the difference in the chemical structure between RNA and DNA is the presence of a hydroxyl group at position 2' of the sugar (ribose) instead of a hydrogen (deoxyribose). In this paper, we present molecular dynamics simulations addressing the contribution of 2'-hydroxyls to the stability of a GCUAA pentaloop motif. The results indicate that the 2'-hydroxyls stabilize the hairpin conformation of the GCUAA pentaloop relative to an analogous oligonucleotide in which the ribose sugars in the loop region were substituted with deoxyriboses. The magnitude of the stabilization was found to be 23.8 ± 4.1 kJ/mol using an alchemical mutations free energy method and 4.2 ± 6.5 kJ/mol using potential of mean force calculations. The latter indicates that in addition to its larger thermodynamic stability the RNA hairpin is also kinetically more stable. We find that the excess stability is a result of intrahairpin hydrogen bonds in the loop region between the 2'-hydroxyls and sugars, bases, and phosphates. The hydrogen bonds with the sugars and phosphates involve predominantly interactions with adjacent nucleotides. However, the hydrogen bonds with the bases involve also interactions between groups on opposite sides of the loop or with the middle base of the loop and are therefore likely to contribute significantly to the stability of the loop. Of these hydrogen bonds, the most frequent is observed between the 2'-hydroxyl at the first position of the pentaloop with N6/N7 of adenine at the forth position, as well as between the 2'-hydroxyl at position -1 with N6 of adenine at the fifth position. Our results contribute to the notion that one of the important roles of the ribose sugars in RNA is to facilitate hairpin formation.
Collapse
Affiliation(s)
- Antonella Paladino
- Department of Organic Chemistry I, University of the Basque Country UPV/EHU, Avenida de Tolosa 72, 20018, San Sebastian, Spain
| | - Ronen Zangi
- Department of Organic Chemistry I, University of the Basque Country UPV/EHU, Avenida de Tolosa 72, 20018, San Sebastian, Spain.,IKERBASQUE, Basque Foundation for Science, 48011, Bilbao, Spain
| |
Collapse
|
21
|
Singharoy A, Joshi H, Miao Y, Ortoleva PJ. Space warping order parameters and symmetry: application to multiscale simulation of macromolecular assemblies. J Phys Chem B 2012; 116:8423-34. [PMID: 22356532 PMCID: PMC4937887 DOI: 10.1021/jp2119247] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Coarse-grained features of macromolecular assemblies are understood via a set of order parameters (OPs) constructed in terms of their all-atom configuration. OPs are shown to be slowly changing in time and capture the large-scale spatial features of macromolecular assemblies. The relationship of these variables to the classic notion of OPs based on symmetry breaking phase transitions is discussed. OPs based on space warping transformations are analyzed in detail as they naturally provide a connection between overall structure of an assembly and all-atom configuration. These OPs serve as the basis of a multiscale analysis that yields Langevin equations for OP dynamics. In this context, the characteristics of OPs and PCA modes are compared. The OPs enable efficient all-atom multiscale simulations of the dynamics of macromolecular assemblies in response to changes in microenvironmental conditions, as demonstrated on the structural transitions of cowpea chlorotic mottle virus capsid (CCMV) and RNA of the satellite tobacco mosaic virus (STMV).
Collapse
Affiliation(s)
- Abhishek Singharoy
- Center for Cell and Virus Theory, Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Harshad Joshi
- Center for Cell and Virus Theory, Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | | | - Peter J. Ortoleva
- Center for Cell and Virus Theory, Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
22
|
Singharoy A, Cheluvaraja S, Ortoleva P. Order parameters for macromolecules: application to multiscale simulation. J Chem Phys 2011; 134:044104. [PMID: 21280684 DOI: 10.1063/1.3524532] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Order parameters (OPs) characterizing the nanoscale features of macromolecules are presented. They are generated in a general fashion so that they do not need to be redesigned with each new application. They evolve on time scales much longer than 10(-14) s typical for individual atomic collisions/vibrations. The list of OPs can be automatically increased, and completeness can be determined via a correlation analysis. They serve as the basis of a multiscale analysis that starts with the N-atom Liouville equation and yields rigorous Smoluchowski/Langevin equations of stochastic OP dynamics. Such OPs and the multiscale analysis imply computational algorithms that we demonstrate in an application to ribonucleic acid structural dynamics for 50 ns.
Collapse
Affiliation(s)
- A Singharoy
- Center for Cell and Virus Theory, Indiana University, Bloomington, Indiana 47405, USA
| | | | | |
Collapse
|
23
|
Gong Z, Xiao Y, Xiao Y. RNA stability under different combinations of amber force fields and solvation models. J Biomol Struct Dyn 2011; 28:431-41. [PMID: 20919758 DOI: 10.1080/07391102.2010.10507372] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.
Collapse
Affiliation(s)
- Zhou Gong
- Biomolecular Physics and Modeling Group, Department of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | | | | |
Collapse
|
24
|
Bandyopadhyay P, Kharerin H. Monte Carlo Energy Landscape Paving and Basin Paving simulation of RNA T-loop hairpin. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2010.12.019] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
25
|
Pasquali S, Derreumaux P. HiRE-RNA: a high resolution coarse-grained energy model for RNA. J Phys Chem B 2010; 114:11957-66. [PMID: 20795690 DOI: 10.1021/jp102497y] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Although RNAs play many cellular functions, the gap between their sequences and 3D structures is increasing and our knowledge of RNA thermodynamics and long time scale dynamics is still limited at an atomic level of detail. In principle, all-atom molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations can investigate these issues, but with current computer facilities, these simulations in explicit solvent have been limited to small RNAs and to short times. To move to larger systems, we can resort to coarse-graining. In this study, we present HiRE-RNA, a generic high resolution coarse-grained model for RNA, and report MD and REMD simulations on two RNAs of 22 and 36 nucleotides. Starting from unfolded structures, the 22-mer folds within 1.8 A rmsd from the NMR structure, while the 36-mer folds within 4.6 A rmsd. Current results suggest that further optimization of the HiRE-RNA force field should open the door to a relevant model for studying large RNAs, such as riboswitches, and for predicting 3D structures from secondary structure information.
Collapse
Affiliation(s)
- Samuela Pasquali
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Institut de Biologie Physico Chimique et Université Paris Diderot-Paris 7, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | | |
Collapse
|
26
|
Kurisaki I, Fukuzawa K, Nakano T, Mochizuki Y, Watanabe H, Tanaka S. Fragment molecular orbital (FMO) study on stabilization mechanism of neuro-oncological ventral antigen (NOVA)–RNA complex system. ACTA ACUST UNITED AC 2010. [DOI: 10.1016/j.theochem.2010.09.013] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
27
|
Sarkar K, Nguyen DA, Gruebele M. Loop and stem dynamics during RNA hairpin folding and unfolding. RNA (NEW YORK, N.Y.) 2010; 16:2427-34. [PMID: 20962040 PMCID: PMC2995403 DOI: 10.1261/rna.2253310] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2010] [Accepted: 08/26/2010] [Indexed: 05/24/2023]
Abstract
2-Aminopurine (2AP) is a fluorescent adenine analog that probes mainly base stacking in nucleic acids. We labeled the loop or the stem of the RNA hairpin gacUACGguc with 2AP to study folding thermodynamics and kinetics at both loci. Thermal melts and fast laser temperature jumps detected by 2AP fluorescence monitored the stability and folding/unfolding kinetics. The observed thermodynamic and kinetic traces of the stem and loop mutants, though strikingly different at a first glance, can be fitted to the same free-energy landscape. The differences between the two probe locations arise because base stacking decreases upon unfolding in the stem, whereas it increases in the loop. We conclude that 2AP is a conservative adenine substitution for mapping out the contributions of different RNA structural elements to the overall folding process. Molecular dynamics (MD) totaling 0.6 μsec were performed to look at the conformations populated by the RNA at different temperatures. The combined experimental data, and MD simulations lead us to propose a minimal four-state free-energy landscape for the RNA hairpin. Analysis of this landscape shows that a sequential folding model is a good approximation for the full folding dynamics. The frayed state formed initially from the native state is a heterogeneous ensemble of structures whose stem is frayed either from the end or from the loop.
Collapse
Affiliation(s)
- Krishnarjun Sarkar
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, USA.
| | | | | |
Collapse
|
28
|
Banáš P, Hollas D, Zgarbová M, Jurečka P, Orozco M, Cheatham TE, Šponer J, Otyepka M. Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins. J Chem Theory Comput 2010; 6:3836-3849. [DOI: 10.1021/ct100481h] [Citation(s) in RCA: 293] [Impact Index Per Article: 20.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Daniel Hollas
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Modesto Orozco
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Thomas E. Cheatham
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic, Joint Research Program in Computational Biology, Institut de Recerca Biomédica and Barcelona Superocomputing Center, Baldiri i Reixac 10, Barcelona 08028, Spain, Jordi Girona 31, Barcelona 08028, Spain
| |
Collapse
|
29
|
Halder S, Bhattacharyya D. Structural Stability of Tandemly Occurring Noncanonical Basepairs within Double Helical Fragments: Molecular Dynamics Studies of Functional RNA. J Phys Chem B 2010; 114:14028-40. [DOI: 10.1021/jp102835t] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Sukanya Halder
- Biophysics Division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata—700 064, India
| | - Dhananjay Bhattacharyya
- Biophysics Division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata—700 064, India
| |
Collapse
|
30
|
Bowman GR, Huang X, Pande VS. Network models for molecular kinetics and their initial applications to human health. Cell Res 2010; 20:622-30. [PMID: 20421891 PMCID: PMC4441225 DOI: 10.1038/cr.2010.57] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
Molecular kinetics underlies all biological phenomena and, like many other biological processes, may best be understood in terms of networks. These networks, called Markov state models (MSMs), are typically built from physical simulations. Thus, they are capable of quantitative prediction of experiments and can also provide an intuition for complex conformational changes. Their primary application has been to protein folding; however, these technologies and the insights they yield are transferable. For example, MSMs have already proved useful in understanding human diseases, such as protein misfolding and aggregation in Alzheimer's disease.
Collapse
Affiliation(s)
- Gregory R Bowman
- Biophysics Program, Stanford University, Stanford, CA 94305, USA
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Vijay S Pande
- Biophysics Program, Stanford University, Stanford, CA 94305, USA
- Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
31
|
Riccardi L, Nguyen PH, Stock G. Free-energy landscape of RNA hairpins constructed via dihedral angle principal component analysis. J Phys Chem B 2010; 113:16660-8. [PMID: 20028141 DOI: 10.1021/jp9076036] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
To systematically construct a low-dimensional free-energy landscape of RNA systems from a classical molecular dynamics simulation, various versions of the principal component analysis (PCA) are compared: the cPCA using the Cartesian coordinates of all atoms, the dPCA using the sine/cosine-transformed six backbone dihedral angles as well as the glycosidic torsional angle chi and the pseudorotational angle P, the aPCA which ignores the circularity of the 6 + 2 dihedral angles of the RNA, and the dPCA(etatheta), which approximates the 6 backbone dihedral angles by 2 pseudotorsional angles eta and theta. As representative examples, a 10-nucleotide UUCG hairpin and the 36-nucleotide segment SL1 of the Psi site of HIV-1 are studied by classical molecular dynamics simulation, using the Amber all-atom force field and explicit solvent. It is shown that the conformational heterogeneity of the RNA hairpins can only be resolved by an angular PCA such as the dPCA but not by the cPCA using Cartesian coordinates. Apart from possible artifacts due to the coupling of overall and internal motion, this is because the details of hydrogen bonding and stacking interactions but also of global structural rearrangements of the RNA are better discriminated by dihedral angles. In line with recent experiments, it is found that the free energy landscape of RNA hairpins is quite rugged and contains various metastable conformational states which may serve as an intermediate for unfolding.
Collapse
Affiliation(s)
- Laura Riccardi
- Institute of Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, D-60438 Frankfurt, Germany
| | | | | |
Collapse
|
32
|
Deng NJ, Cieplak P. Free energy profile of RNA hairpins: a molecular dynamics simulation study. Biophys J 2010; 98:627-36. [PMID: 20159159 PMCID: PMC2820654 DOI: 10.1016/j.bpj.2009.10.040] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2009] [Revised: 10/25/2009] [Accepted: 10/26/2009] [Indexed: 11/25/2022] Open
Abstract
RNA hairpin loops are one of the most abundant secondary structure elements and participate in RNA folding and protein-RNA recognition. To characterize the free energy surface of RNA hairpin folding at an atomic level, we calculated the potential of mean force (PMF) as a function of the end-to-end distance, by using umbrella sampling simulations in explicit solvent. Two RNA hairpins containing tetraloop cUUCGg and cUUUUg are studied with AMBER ff99 and CHARMM27 force fields. Experimentally, the UUCG hairpin is known to be significantly more stable than UUUU. In this study, the calculations using AMBER force field give a qualitatively correct description for the folding of two RNA hairpins, as the calculated PMF confirms the global stability of the folded structures and the resulting relative folding free energy is in quantitative agreement with the experimental result. The hairpin stabilities are also correctly differentiated by the more rapid molecular mechanics-Poisson Boltzmann-surface area approach, but the relative free energy estimated from this method is overestimated. The free energy profile shows that the native state basin and the unfolded state plateau are separated by a wide shoulder region, which samples a variety of native-like structures with frayed terminal basepair. The calculated PMF lacks major barriers that are expected near the transition regions, and this is attributed to the limitation of the 1-D reaction coordinate. The PMF results are compared with other studies of small RNA hairpins using kinetics method and coarse grained models. The two RNA hairpins described by CHARMM27 are significantly more deformable than those represented by AMBER. Compared with the AMBER results, the CHARMM27 calculated DeltaG(fold) for the UUUU tetraloop is in better agreement with the experimental results. However, the CHARMM27 calculation does not confirm the global stability of the experimental UUCG structure; instead, the extended conformations are predicted to be thermodynamically stable in solution. This finding is further supported by separate unrestrained CHARMM27 simulations, in which the UUCG hairpin unfolds spontaneously within 10 ns. The instability of the UUCG hairpin originates from the loop region, and propagates to the stem. The results of this study provide a molecular picture of RNA hairpin unfolding and reveal problems in the force field descriptions for the conformational energy of certain RNA hairpin.
Collapse
Affiliation(s)
- Nan-Jie Deng
- School of Life Science, University of Science and Technology of China Hefei, Anhui, China.
| | | |
Collapse
|
33
|
Sarkar K, Meister K, Sethi A, Gruebele M. Fast folding of an RNA tetraloop on a rugged energy landscape detected by a stacking-sensitive probe. Biophys J 2009; 97:1418-27. [PMID: 19720030 DOI: 10.1016/j.bpj.2009.06.035] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2009] [Revised: 06/25/2009] [Accepted: 06/30/2009] [Indexed: 10/20/2022] Open
Abstract
We investigate the microsecond-timescale kinetics of the RNA hairpin ga*cUUCGguc. The fluorescent nucleotide 2-aminopurine (a*) reports mainly on base stacking. Ten kinetic traces and the temperature denaturation curve are globally fitted to four-state models of the free-energy surface. In the best-fitting sequential model, the hairpin unfolds over successively larger barriers in at least three stages: stem fraying and increased base-stacking fluctuations; concerted loss of hydrogen bonding and partial unstacking; and additional unstacking of single strands at the highest temperatures. Parallel and trap models also provide adequate fits: such pathways probably also play a role in the complete free-energy surface of the hairpin. To interpret the model states structurally, 200 ns of molecular dynamics, including six temperature-jump simulations, were run. Although the sampling is by no means comprehensive, five different states were identified using hydrogen bonding and base stacking as reaction coordinates. The four to five states required to explain the experiments or simulations set a lower limit on the complexity of this small RNA hairpin's energy landscape.
Collapse
Affiliation(s)
- Krishnarjun Sarkar
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | | | | | | |
Collapse
|
34
|
Rinnenthal J, Richter C, Nozinovic S, Fürtig B, Lopez JJ, Glaubitz C, Schwalbe H. RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis. JOURNAL OF BIOMOLECULAR NMR 2009; 45:143-55. [PMID: 19636800 DOI: 10.1007/s10858-009-9343-x] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2009] [Accepted: 06/19/2009] [Indexed: 05/13/2023]
Abstract
We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.
Collapse
Affiliation(s)
- Jörg Rinnenthal
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance, Johann Wolfgang Goethe-University, Max-von-Laue-Strasse 7, Frankfurt/Main, Germany
| | | | | | | | | | | | | |
Collapse
|
35
|
Nguyen PH, Derreumaux P, Stock G. Energy Flow and Long-Range Correlations in Guanine-Binding Riboswitch: A Nonequilibrium Molecular Dynamics Study. J Phys Chem B 2009; 113:9340-7. [DOI: 10.1021/jp902013s] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Phuong H. Nguyen
- Institute of Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, D-60438 Frankfurt, Germany, and Laboratoire de Biochimie Theorique, UPR 9080 CNRS, IBPC, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Philippe Derreumaux
- Institute of Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, D-60438 Frankfurt, Germany, and Laboratoire de Biochimie Theorique, UPR 9080 CNRS, IBPC, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Gerhard Stock
- Institute of Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, D-60438 Frankfurt, Germany, and Laboratoire de Biochimie Theorique, UPR 9080 CNRS, IBPC, 13 rue Pierre et Marie Curie, 75005, Paris, France
| |
Collapse
|
36
|
Villa A, Wöhnert J, Stock G. Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch. Nucleic Acids Res 2009; 37:4774-86. [PMID: 19515936 PMCID: PMC2724292 DOI: 10.1093/nar/gkp486] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (approximately 1 micros total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.
Collapse
Affiliation(s)
- Alessandra Villa
- Institute of Physical and Theoretical Chemistry and Institute of Molecular Biosciences, Goethe University, Frankfurt am Main, Germany
| | | | | |
Collapse
|
37
|
Theoretical studies on the interaction of modified pyrimidines and purines with purine riboswitch. J Mol Graph Model 2009; 28:37-45. [PMID: 19380244 DOI: 10.1016/j.jmgm.2009.03.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2009] [Revised: 03/18/2009] [Accepted: 03/21/2009] [Indexed: 12/25/2022]
Abstract
Recent experimental study [S.D. Gilbert, S.J. Mediatore, R.T. Batey, Modified pyrimidine specifically bind the purine riboswitch, J. Am. Chem. Soc. 128 (2006) 14214-14215] demonstrated that the purine riboswitch could specifically bind some ligands other than purines such as amino-pyrimidines, and the authors proposed that the five-membered ring of purine was not required for recognition. To get insight into the interaction details, we used molecular docking method to investigate the interactions of a mutant form of guanine riboswitch with a series of amino-purines, amino-pyrimidines and imidazole derivatives, and employed molecular simulation method to study the dynamic behavior of the selected complexes. The calculation results reveal that (1) all the amino-purines and amino-pyrimidines bind in a same cavity composed of four nucleobases including U22, U47, U51 and U74, which is consistent with the experimental results, while the two imidazole derivatives adopt other binding modes; (2) the purines are engulfed within three-way junction motifs, but most pyrimidines only form two-way junctions with the riboswitch; (3) the number and position of amino substituents could seriously affect the binding of pyrimidines. As riboswitches are potentially excellent candidates for antibiotic therapeutics, these findings may be useful for understanding the range of compounds that riboswitch can specifically recognize.
Collapse
|
38
|
Zhang Y, Zhao X, Mu Y. Conformational Transition Map of an RNA GCAA Tetraloop Explored by Replica-Exchange Molecular Dynamics Simulation. J Chem Theory Comput 2009; 5:1146-54. [DOI: 10.1021/ct8004276] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yufen Zhang
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, P.R. China
| | - Xian Zhao
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, P.R. China
| | - Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, P.R. China
| |
Collapse
|
39
|
KURISAKI I, WATANABE H, TANAKA S. Simulation Study of RNA-Binding Protein, Pumilio. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2009. [DOI: 10.2477/jccj.h2015] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
40
|
Hall KB. RNA in motion. Curr Opin Chem Biol 2008; 12:612-8. [PMID: 18957331 DOI: 10.1016/j.cbpa.2008.09.033] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2008] [Revised: 09/05/2008] [Accepted: 09/20/2008] [Indexed: 10/21/2022]
Abstract
Although RNA duplex regions are highly structured and inflexible, other elements of an RNA molecule are capable of dynamic motions. These flexible regions are the sites of interactions with small molecules, proteins, and other RNAs, yet there are few descriptions of these regions that include the timescale and amplitude of their motions. No one technique is sufficient to accurately describe these motions, but the combination of in vitro methods, particularly NMR relaxation methods, and more robust in silico methods, is beginning to yield the type of data that can be used to understand RNA function. Very few RNAs have been described by both techniques, and here one such RNA and one RNA:protein complex are reviewed.
Collapse
Affiliation(s)
- Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St Louis, MO 63110, USA.
| |
Collapse
|
41
|
Denschlag R, Lingenheil M, Tavan P. Efficiency reduction and pseudo-convergence in replica exchange sampling of peptide folding–unfolding equilibria. Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2008.04.114] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
42
|
Ferner J, Villa A, Duchardt E, Widjajakusuma E, Wöhnert J, Stock G, Schwalbe H. NMR and MD studies of the temperature-dependent dynamics of RNA YNMG-tetraloops. Nucleic Acids Res 2008; 36:1928-40. [PMID: 18272534 PMCID: PMC2346598 DOI: 10.1093/nar/gkm1183] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.
Collapse
Affiliation(s)
- Jan Ferner
- Institut für Organische Chemie und Chemische Biologie, Center for Biomolecular Magnetic Resonance, Frankfurt/M, Germany
| | | | | | | | | | | | | |
Collapse
|