1
|
Röder K, Pasquali S. Assessing RNA atomistic force fields via energy landscape explorations in implicit solvent. Biophys Rev 2024; 16:285-295. [PMID: 39099837 PMCID: PMC11297004 DOI: 10.1007/s12551-024-01202-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 05/29/2024] [Indexed: 08/06/2024] Open
Abstract
Predicting the structure and dynamics of RNA molecules still proves challenging because of the relative scarcity of experimental RNA structures on which to train models and the very sensitive nature of RNA towards its environment. In the last decade, several atomistic force fields specifically designed for RNA have been proposed and are commonly used for simulations. However, it is not necessarily clear which force field is the most suitable for a given RNA molecule. In this contribution, we propose the use of the computational energy landscape framework to explore the energy landscape of RNA systems as it can bring complementary information to the more standard approaches of enhanced sampling simulations based on molecular dynamics. We apply the EL framework to the study of a small RNA pseudoknot, the Aquifex aeolicus tmRNA pseudoknot PK1, and we compare the results of five different RNA force fields currently available in the AMBER simulation software, in implicit solvent. With this computational approach, we can not only compare the predicted 'native' states for the different force fields, but the method enables us to study metastable states as well. As a result, our comparison not only looks at structural features of low energy folded structures, but provides insight into folding pathways and higher energy excited states, opening to the possibility of assessing the validity of force fields also based on kinetics and experiments providing information on metastable and unfolded states. Supplementary Information The online version contains supplementary material available at 10.1007/s12551-024-01202-9.
Collapse
Affiliation(s)
- Konstantin Röder
- Randall Centre for Cell & Molecular Biophysics, King’s College London, London, SE1 1UL UK
| | - Samuela Pasquali
- Laboratoire Biologie Functionnelle Et Adaptative, CNRS UMR 8251, Inserm ERL U1133, Université Paris Cité , 35 Rue Hélène Brion, Paris, France
| |
Collapse
|
2
|
Ramachandran V, Potoyan DA. Energy landscapes of homopolymeric RNAs revealed by deep unsupervised learning. Biophys J 2024; 123:1152-1163. [PMID: 38571310 PMCID: PMC11079944 DOI: 10.1016/j.bpj.2024.04.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 03/03/2024] [Accepted: 04/01/2024] [Indexed: 04/05/2024] Open
Abstract
Conformational dynamics of RNA plays important roles in a variety of cellular functions such as transcriptional regulation, catalysis, scaffolding, and sensing. Recently, RNAs with low-complexity sequences have been shown to phase separate and form condensate phases similar to lowcomplexity protein domains. The affinity for phase separation and the material characteristics of RNA condensates are strongly dependent on sequence composition and patterning. We hypothesize that differences in the affinities for RNA phase separation can be uncovered by studying sequence-dependent conformational dynamics of single RNA chains. To this end, we have employed atomistic simulations and deep dimensionality reduction techniques to map temperature-dependent conformational free energy landscapes for 20 base-long homopolymeric RNA sequences: poly(U), poly(G), poly(C), and poly(A). The energy landscapes of homopolymeric RNAs reveal a plethora of metastable states with qualitatively different populations stemming from differences in base chemistry. Through detailed analysis of base, phosphate, and sugar interactions, we show that experimentally observed temperature-driven shifts in metastable state populations align with experiments on RNA phase transitions. Specifically, we find that the thermodynamics of unfolding of homopolymeric RNA follows the poly(G) > poly(A) > poly(C) > poly(U) order of stability, mirroring the propensity of RNA to form condensates. To conclude, this work shows that at least for homopolymeric RNA sequences the single-chain conformational dynamics contains sufficient information for predicting and quantifying condensate forming affinities of RNAs. Thus, we anticipate that atomically detailed studies of temeprature -dependent energy landscapes of RNAs will be a useful guide for understanding the propensity of various RNA molecules to form condensates.
Collapse
Affiliation(s)
| | - Davit A Potoyan
- Department of Chemistry, Iowa State University, Ames, Iowa; Department of Biochemistry Biophysics and Molecular Biology, Iowa State University, Ames, Iowa.
| |
Collapse
|
3
|
Gilardoni I, Fröhlking T, Bussi G. Boosting Ensemble Refinement with Transferable Force-Field Corrections: Synergistic Optimization for Molecular Simulations. J Phys Chem Lett 2024; 15:1204-1210. [PMID: 38272001 DOI: 10.1021/acs.jpclett.3c03423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
A novel method combining the force-field fitting approach and ensemble refinement by the maximum entropy principle is presented. Its formulation allows us to continuously interpolate between these two methods, which can thus be interpreted as two limiting cases. A cross-validation procedure enables us to correctly assess the relative weight of both of them, distinguishing scenarios in which the combined approach is meaningful from those in which either ensemble refinement or force-field fitting separately prevails. The efficacy of their combination is examined for a realistic case study of RNA oligomers. Within the new scheme, molecular dynamics simulations are integrated with experimental data provided by nuclear magnetic resonance measures. We show that force-field corrections are in general superior when applied to the appropriate force-field terms but are automatically discarded by the method when applied to inappropriate force-field terms.
Collapse
Affiliation(s)
- Ivan Gilardoni
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
4
|
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
|
5
|
Mlýnský V, Kührová P, Stadlbauer P, Krepl M, Otyepka M, Banáš P, Šponer J. Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations. J Chem Theory Comput 2023; 19:8423-8433. [PMID: 37944118 PMCID: PMC10687871 DOI: 10.1021/acs.jctc.3c00990] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2023] [Revised: 10/18/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| | - Petra Kührová
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Petr Stadlbauer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
| | - Michal Otyepka
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
- Czech
Advanced Technology and Research Institute, CATRIN, Křížkovského 511/8, Olomouc 779 00, Czech Republic
- IT4Innovations, VSB−Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 612 00, Czech Republic
| |
Collapse
|
6
|
Viegas RG, Sanches MN, Chen AA, Paulovich FV, Garcia AE, Leite VBP. Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop. J Chem Inf Model 2023; 63:5641-5649. [PMID: 37606640 DOI: 10.1021/acs.jcim.3c00426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.
Collapse
Affiliation(s)
- Rafael G Viegas
- Federal Institute of Education, Science and Technology of São Paulo (IFSP), Catanduva, São Paulo 15.808-305, Brazil
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Murilo N Sanches
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Alan A Chen
- Department of Chemistry and the RNA Institute, University at Albany, Albany, New York 12222, United States
| | - Fernando V Paulovich
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, Eindhoven 5600 MB, the Netherlands
| | - Angel E Garcia
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Physics Division, National Science Foundation, 2415 Eisenhower Ave, Alexandria, Virginia 22314, United States
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|
7
|
Rahimi K, Piaggi PM, Zerze GH. Comparison of On-the-Fly Probability Enhanced Sampling and Parallel Tempering Combined with Metadynamics for Atomistic Simulations of RNA Tetraloop Folding. J Phys Chem B 2023. [PMID: 37196167 DOI: 10.1021/acs.jpcb.3c00117] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
Atomistic simulations with reliable models are extremely useful in providing exquisitely detailed pictures of biomolecular phenomena that are not always accessible to experiments. One such biomolecular phenomenon is RNA folding, which often requires exhaustive simulations with combined advanced sampling techniques. In this work, we employed the multithermal-multiumbrella on-the-fly probability enhanced sampling (MM-OPES) technique and compared it against combined parallel tempering and metadynamics simulations. We found that MM-OPES simulations were successful in reproducing the free energy surfaces from combined parallel tempering and metadynamics simulations. Importantly, we also investigated a wide range of temperature sets (minimum and maximum temperatures) for MM-OPES simulations in order to identify some guidelines for deciding the temperature limits for an accurate and efficient exploration of the free energy landscapes. We found that most temperature sets yielded almost the same accuracy in reproducing the free energy surface at the ambient conditions as long as (i) the maximum temperature is reasonably high, (ii) the temperature at which we run the simulation is reasonably high (in our simulations, running temperature is defined as [minimum temperature + maximum temperature]/2), and (iii) the effective sample size at the temperature of interest is statistically reasonable. In terms of the computational cost, all MM-OPES simulations were nearly 4 times less costly than the combined parallel tempering and metadynamics simulations. We concluded that the demanding combined parallel tempering and metadynamics simulations can safely be replaced with approximately 4 times less costly MM-OPES simulations (with carefully selected temperature limits) to obtain the same information.
Collapse
Affiliation(s)
- Kosar Rahimi
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Gül H Zerze
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| |
Collapse
|
8
|
Zhang D, Gong L, Weng J, Li Y, Wang A, Li G. RNA Folding Based on 5 Beads Model and Multiscale Simulation. Interdiscip Sci 2023:10.1007/s12539-023-00561-3. [PMID: 37115389 DOI: 10.1007/s12539-023-00561-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 03/08/2023] [Accepted: 03/10/2023] [Indexed: 04/29/2023]
Abstract
RNA folding prediction is very meaningful and challenging. The molecular dynamics simulation (MDS) of all atoms (AA) is limited to the folding of small RNA molecules. At present, most of the practical models are coarse grained (CG) model, and the coarse-grained force field (CGFF) parameters usually depend on known RNA structures. However, the limitation of the CGFF is obvious that it is difficult to study the modified RNA. Based on the 3 beads model (AIMS_RNA_B3), we proposed the AIMS_RNA_B5 model with three beads representing a base and two beads representing the main chain (sugar group and phosphate group). We first run the all atom molecular dynamic simulation (AAMDS), and fit the CGFF parameter with the AA trajectory. Then perform the coarse-grained molecular dynamic simulation (CGMDS). AAMDS is the foundation of CGMDS. CGMDS is mainly to carry out the conformation sampling based on the current AAMDS state and improve the folding speed. We simulated the folding of three RNAs, which belong to hairpin, pseudoknot and tRNA respectively. Compared to the AIMS_RNA_B3 model, the AIMS_RNA_B5 model is more reasonable and performs better.
Collapse
Affiliation(s)
- Dinglin Zhang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
- Chinese Academy of Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Lidong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, China
| | - Junben Weng
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
- Chinese Academy of Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Yan Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
| | - Anhui Wang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China.
| |
Collapse
|
9
|
Winkler L, Cheatham TE. Benchmarking the Drude Polarizable Force Field Using the r(GACC) Tetranucleotide. J Chem Inf Model 2023; 63:2505-2511. [PMID: 36996447 PMCID: PMC10131218 DOI: 10.1021/acs.jcim.3c00250] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/01/2023]
Abstract
Polarizable force fields, in particular, the Drude polarizable force field (Drude FF), may hold the key to more accurately modeling biomolecules with molecular dynamics simulations by explicitly accounting for atomic polarizability. Previous work has shown promising results in simulating duplex nucleic acids and protein structures with excellent agreement with experimental values. However, benchmarking the Drude polarizable force field with highly flexible, single-stranded structures has yet to be achieved. In this work, the r(GACC) tetranucleotide is simulated over a multimicrosecond time scale, starting with various different initial conformations. Despite the starting conformation, including starting from the expected dominant A-form major conformation, the experimental structural distribution is not matched. In fact, the major NMR conformation is never resampled. Instead, the r(GACC) tetranucleotide becomes stabilized in anomalous structures that are inconsistent with the NMR data and that favor base-pairing and electrostatic interactions over base stacking. These structures are maintained for lengthy time scales (>1 μs) themselves, suggesting a misbalance of forces in the Drude polarizable force field itself. This model system is suggestive of the fact that currently the Drude polarizable force field does not appear to produce the sensitive balance of forces required to accurately model other single-stranded or noncanonical RNA structures.
Collapse
Affiliation(s)
- Lauren Winkler
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| |
Collapse
|
10
|
Mráziková K, Kruse H, Mlýnský V, Auffinger P, Šponer J. Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions. J Chem Inf Model 2022; 62:6182-6200. [PMID: 36454943 DOI: 10.1021/acs.jcim.2c01064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg67084, France
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| |
Collapse
|
11
|
Zirbel CL, Auffinger P. Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey. Molecules 2022; 27:molecules27144365. [PMID: 35889236 PMCID: PMC9323530 DOI: 10.3390/molecules27144365] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2022] [Revised: 06/29/2022] [Accepted: 07/04/2022] [Indexed: 02/04/2023] Open
Abstract
Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.
Collapse
Affiliation(s)
- Craig L. Zirbel
- Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, OH 43403, USA;
| | - Pascal Auffinger
- Architecture et Réactivité de l’ARN, UPR 9002, Institut de Biologie Moléculaire et Cellulaire du CNRS, Université de Strasbourg, 67084 Strasbourg, France
- Correspondence: ; Tel.: +33-3-8841-7049; Fax: +33-3-8860-2218
| |
Collapse
|
12
|
Fröhlking T, Mlýnský V, Janeček M, Kührová P, Krepl M, Banáš P, Šponer J, Bussi G. Automatic Learning of Hydrogen-Bond Fixes in the AMBER RNA Force Field. J Chem Theory Comput 2022; 18:4490-4502. [PMID: 35699952 PMCID: PMC9281393 DOI: 10.1021/acs.jctc.2c00200] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
![]()
The
capability of
current force fields to reproduce RNA structural
dynamics is limited. Several methods have been developed to take advantage
of experimental data in order to enforce agreement with experiments.
Here, we extend an existing framework which allows arbitrarily chosen
force-field correction terms to be fitted by quantification of the
discrepancy between observables back-calculated from simulation and
corresponding experiments. We apply a robust regularization protocol
to avoid overfitting and additionally introduce and compare a number
of different regularization strategies, namely, L1, L2, Kish size,
relative Kish size, and relative entropy penalties. The training set
includes a GACC tetramer as well as more challenging systems, namely,
gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen
bonds in the AMBER RNA force field are corrected with automatically
determined parameters that we call gHBfixopt. A validation
involving a separate simulation of a system present in the training
set (gcUUCGgc) and new systems not seen during training (CAAU and
UUUU tetramers) displays improvements regarding the native population
of the tetraloop as well as good agreement with NMR experiments for
tetramers when using the new parameters. Then, we simulate folded
RNAs (a kink–turn and L1 stalk rRNA) including hydrogen bond
types not sufficiently present in the training set. This allows a
final modification of the parameter set which is named gHBfix21 and
is suggested to be applicable to a wider range of RNA systems.
Collapse
Affiliation(s)
- Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, Brno 612 65, Czech Republic
| | - Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, Olomouc 771 46, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - 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, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, Brno 612 65, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 27, 779 00 Olomouc, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
| |
Collapse
|
13
|
Mlýnský V, Janeček M, Kührová P, Fröhlking T, Otyepka M, Bussi G, Banáš P, Šponer J. Toward Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Modifications. J Chem Theory Comput 2022; 18:2642-2656. [PMID: 35363478 DOI: 10.1021/acs.jctc.1c01222] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔGfold° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic.,IT4Innovations, VSB─Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava-Poruba, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| |
Collapse
|
14
|
Chen J, Liu H, Cui X, Li Z, Chen HF. RNA-Specific Force Field Optimization with CMAP and Reweighting. J Chem Inf Model 2022; 62:372-385. [PMID: 35021622 DOI: 10.1021/acs.jcim.1c01148] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
RNA plays a key role in a variety of cell activities. However, it is difficult to capture its structure dynamics by the traditional experimental methods because of the inherent limitations. Molecular dynamics simulation has become a valuable complement to the experimental methods. Previous studies have indicated that the current force fields cannot accurately reproduce the conformations and structural dynamics of RNA. Therefore, an RNA-specific force field was developed to improve the conformation sampling of RNA. The distribution of ζ/α dihedrals of tetranucleotides was optimized by a reweighting method, and the grid-based energy correction map (CMAP) term was first introduced into the Amber RNA force field of ff99bsc0χOL3, named ff99OL3_CMAP1. Extensive validations of tetranucleotides and tetraloops show that ff99OL3_CMAP1 can significantly decrease the population of an incorrect structure, increase the consistency between the simulation results and experimental values for tetranucleotides, and improve the stability of tetraloops. ff99OL3_CMAP1 can also precisely reproduce the conformation of a duplex and riboswitches. These findings confirm that the newly developed force field ff99OL3_CMAP1 can improve the conformer sampling of RNA.
Collapse
Affiliation(s)
- Jun Chen
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, 20024 Shanghai, China
| | - Hao Liu
- Institute of Natural Sciences, Shanghai Jiao Tong University, 200240 Shanghai, China
| | - Xiaochen Cui
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, 20024 Shanghai, China
| | - Zhengxin Li
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, 20024 Shanghai, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, 20024 Shanghai, China.,Shanghai Center for Bioinformation Technology, 200240 Shanghai, China
| |
Collapse
|
15
|
Zhao J, Kennedy SD, Turner DH. Nuclear Magnetic Resonance Spectra and AMBER OL3 and ROC-RNA Simulations of UCUCGU Reveal Force Field Strengths and Weaknesses for Single-Stranded RNA. J Chem Theory Comput 2022; 18:1241-1254. [PMID: 34990548 DOI: 10.1021/acs.jctc.1c00643] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Single-stranded regions of RNA are important for folding of sequences into 3D structures and for design of therapeutics targeting RNA. Prediction of ensembles of 3D structures for single-stranded regions often involves classical mechanical approximations of interactions defined by quantum mechanical calculations on small model systems. Nuclear magnetic resonance (NMR) spectra and molecular dynamics (MD) simulations of short single strands provide tests for how well the approximations model many of the interactions. Here, the NMR spectra for UCUCGU at 2, 15, and 30 °C are compared to simulations with the AMBER force fields, OL3 and ROC-RNA. This is the first such comparison to an oligoribonucleotide containing an internal guanosine nucleotide (G). G is particularly interesting because of its many H-bonding groups, large dipole moment, and proclivity for both syn and anti conformations. Results reveal formation of a G amino to phosphate non-bridging oxygen H-bond. The results also demonstrate dramatic differences in details of the predicted structures. The variations emphasize the dependence of predictions on individual parameters and their balance with the rest of the force field. The NMR data can serve as a benchmark for future force fields.
Collapse
|
16
|
Zerze GH, Piaggi PM, Debenedetti PG. A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States. J Phys Chem B 2021; 125:13685-13695. [PMID: 34890201 DOI: 10.1021/acs.jpcb.1c08038] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
17
|
Bernetti M, Hall KB, Bussi G. Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles. Nucleic Acids Res 2021; 49:e84. [PMID: 34107023 PMCID: PMC8373061 DOI: 10.1093/nar/gkab459] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 05/07/2021] [Accepted: 05/16/2021] [Indexed: 01/03/2023] Open
Abstract
Small-angle X-ray scattering (SAXS) experiments are increasingly used to probe RNA structure. A number of forward models that relate measured SAXS intensities and structural features, and that are suitable to model either explicit-solvent effects or solute dynamics, have been proposed in the past years. Here, we introduce an approach that integrates atomistic molecular dynamics simulations and SAXS experiments to reconstruct RNA structural ensembles while simultaneously accounting for both RNA conformational dynamics and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward models and enhanced sampling methods to sample an heterogenous ensemble of structures, with no information towards the experiments provided on-the-fly. The generated structural ensemble is then reweighted through the maximum entropy principle so as to match reference SAXS experimental data at multiple ionic conditions. Importantly, accurate explicit-solvent forward models are used at this reweighting stage. We apply this framework to the GTPase-associated center, a relevant RNA molecule involved in protein translation, in order to elucidate its ion-dependent conformational ensembles. We show that (a) both solvent and dynamics are crucial to reproduce experimental SAXS data and (b) the resulting dynamical ensembles contain an ion-dependent fraction of extended structures.
Collapse
Affiliation(s)
- Mattia Bernetti
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| | - Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| |
Collapse
|
18
|
Janeček M, Kührová P, Mlýnský V, Otyepka M, Šponer J, Banáš P. W-RESP: Well-Restrained Electrostatic Potential-Derived Charges. Revisiting the Charge Derivation Model. J Chem Theory Comput 2021; 17:3495-3509. [PMID: 33999623 DOI: 10.1021/acs.jctc.0c00976] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Representation of electrostatic interactions by a Coulombic pairwise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force-field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, the description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules and thus electrostatic interactions in the force field. In particular, we reoptimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme, and included extra point (EP) charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model shows some potential for improving the performance of the nucleic acid force fields, for which the poor description of nonbonded interactions, such as the underestimated stability of base pairing, is well-established. We also report some preliminary simulation tests (around 1 ms of simulation data) on A-RNA duplexes, tetranucleotides, and tetraloops. The simulations reveal no adverse effects, while the description of base-pairing interactions might be improved. The new charges can thus be used in future attempts to improve the nucleic acid simulation force fields, in combination with reparametrization of the other terms.
Collapse
Affiliation(s)
- Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic
| | - Vojtěch Mlýnský
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
19
|
Zerze GH, Stillinger FH, Debenedetti PG. Thermodynamics of DNA Hybridization from Atomistic Simulations. J Phys Chem B 2021; 125:771-779. [PMID: 33434025 DOI: 10.1021/acs.jpcb.0c09237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Frank H Stillinger
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
20
|
Mráziková K, Mlýnský V, Kührová P, Pokorná P, Kruse H, Krepl M, Otyepka M, Banáš P, Šponer J. UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations. J Chem Theory Comput 2020; 16:7601-7617. [PMID: 33215915 DOI: 10.1021/acs.jctc.0c00801] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
21
|
Halder A, Kumar S, Valsson O, Reddy G. Mg 2+ Sensing by an RNA Fragment: Role of Mg 2+-Coordinated Water Molecules. J Chem Theory Comput 2020; 16:6702-6715. [PMID: 32941038 DOI: 10.1021/acs.jctc.0c00589] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg2+, Ca2+, and K+) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg2+ solution. We show that the rigid network of Mg2+-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K+ and Ca2+ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.
Collapse
Affiliation(s)
- Antarip Halder
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Sunil Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - Omar Valsson
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| |
Collapse
|
22
|
Salsbury AM, Lemkul JA. Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics. Curr Opin Struct Biol 2020; 67:9-17. [PMID: 32950937 DOI: 10.1016/j.sbi.2020.08.003] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 08/04/2020] [Accepted: 08/13/2020] [Indexed: 01/24/2023]
Abstract
Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States; Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, United States.
| |
Collapse
|
23
|
Sakuraba S, Iwakiri J, Hamada M, Kameda T, Tsuji G, Kimura Y, Abe H, Asai K. Free-Energy Calculation of Ribonucleic Inosines and Its Application to Nearest-Neighbor Parameters. J Chem Theory Comput 2020; 16:5923-5935. [PMID: 32786906 DOI: 10.1021/acs.jctc.0c00270] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Can current simulations quantitatively predict the stability of ribonucleic acids (RNAs)? In this research, we apply a free-energy perturbation simulation of RNAs containing inosine, a modified ribonucleic base, to the derivation of RNA nearest-neighbor parameters. A parameter set derived solely from 30 simulations was used to predict the free-energy difference of the RNA duplex with a mean unbiased error of 0.70 kcal/mol, which is a level of accuracy comparable to that obtained with parameters derived from 25 experiments. We further show that the error can be lowered to 0.60 kcal/mol by combining the simulation-derived free-energy differences with experimentally measured differences. This protocol can be used as a versatile method for deriving nearest-neighbor parameters of RNAs with various modified bases.
Collapse
Affiliation(s)
- Shun Sakuraba
- Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, Kyoto 619-0215, Japan.,Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba 277-8561, Japan
| | - Junichi Iwakiri
- Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba 277-8561, Japan
| | - Michiaki Hamada
- Faculty of Science and Engineering, Waseda University, Shinjuku-ku, Tokyo 169-8555, Japan.,Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology (AIST), Shinjuku-ku, Tokyo 169-8555, Japan
| | - Tomoshi Kameda
- Artificial Intelligence Research Center (AIRC), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo 135-0064, Japan
| | - Genichiro Tsuji
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo, Chikusa, Nagoya 464-8602, Japan.,Division of Organic Chemistry, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kawasaki, Kanagawa 210-9501, Japan
| | - Yasuaki Kimura
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo, Chikusa, Nagoya 464-8602, Japan
| | - Hiroshi Abe
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo, Chikusa, Nagoya 464-8602, Japan
| | - Kiyoshi Asai
- Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba 277-8561, Japan.,Artificial Intelligence Research Center (AIRC), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo 135-0064, Japan
| |
Collapse
|
24
|
Bottaro S, Nichols PJ, Vögeli B, Parrinello M, Lindorff-Larsen K. Integrating NMR and simulations reveals motions in the UUCG tetraloop. Nucleic Acids Res 2020; 48:5839-5848. [PMID: 32427326 PMCID: PMC7293013 DOI: 10.1093/nar/gkaa399] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2020] [Revised: 04/03/2020] [Accepted: 05/17/2020] [Indexed: 12/21/2022] Open
Abstract
We provide an atomic-level description of the structure and dynamics of the UUCG RNA stem-loop by combining molecular dynamics simulations with experimental data. The integration of simulations with exact nuclear Overhauser enhancements data allowed us to characterize two distinct states of this molecule. The most stable conformation corresponds to the consensus three-dimensional structure. The second state is characterized by the absence of the peculiar non-Watson-Crick interactions in the loop region. By using machine learning techniques we identify a set of experimental measurements that are most sensitive to the presence of non-native states. We find that although our MD ensemble, as well as the consensus UUCG tetraloop structures, are in good agreement with experiments, there are remaining discrepancies. Together, our results show that (i) the MD simulation overstabilize a non-native loop conformation, (ii) eNOE data support its presence with a population of ≈10% and (iii) the structural interpretation of experimental data for dynamic RNAs is highly complex, even for a simple model system such as the UUCG tetraloop.
Collapse
Affiliation(s)
- Sandro Bottaro
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Parker J Nichols
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Anschutz Medical Campus, Aurora, CO, USA
| | - Michele Parrinello
- Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
25
|
Orioli S, Larsen AH, Bottaro S, Lindorff-Larsen K. How to learn from inconsistencies: Integrating molecular simulations with experimental data. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:123-176. [PMID: 32145944 DOI: 10.1016/bs.pmbts.2019.12.006] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Molecular simulations and biophysical experiments can be used to provide independent and complementary insights into the molecular origin of biological processes. A particularly useful strategy is to use molecular simulations as a modeling tool to interpret experimental measurements, and to use experimental data to refine our biophysical models. Thus, explicit integration and synergy between molecular simulations and experiments is fundamental for furthering our understanding of biological processes. This is especially true in the case where discrepancies between measured and simulated observables emerge. In this chapter, we provide an overview of some of the core ideas behind methods that were developed to improve the consistency between experimental information and numerical predictions. We distinguish between situations where experiments are used to refine our understanding and models of specific systems, and situations where experiments are used more generally to refine transferable models. We discuss different philosophies and attempt to unify them in a single framework. Until now, such integration between experiments and simulations have mostly been applied to equilibrium data, and we discuss more recent developments aimed to analyze time-dependent or time-resolved data.
Collapse
Affiliation(s)
- Simone Orioli
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Structural Biophysics, Niels Bohr Institute, Faculty of Science, University of Copenhagen, Copenhagen, Denmark
| | - Andreas Haahr Larsen
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Structural Biophysics, Niels Bohr Institute, Faculty of Science, University of Copenhagen, Copenhagen, Denmark
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark; Atomistic Simulations Laboratory, Istituto Italiano di Tecnologia, Genova, Italy
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
26
|
Zhang H, Gong Q, Zhang H, Chen C. FSATOOL: A useful tool to do the conformational sampling and trajectory analysis work for biomolecules. J Comput Chem 2020; 41:156-164. [PMID: 31603251 DOI: 10.1002/jcc.26083] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Revised: 09/10/2019] [Accepted: 09/12/2019] [Indexed: 12/27/2022]
Abstract
Reliable conformational sampling and trajectory analysis are always important to the study of the folding or binding mechanisms of biomolecules. Generally, one has to prepare many complicated parameters and follow a lot of steps to obtain the final data. The whole process is too complicated to new users. In this article, we provide a convenient and user-friendly tool that is compatible to AMBER, called fast sampling and analysis tool (FSATOOL). FSATOOL has some useful features. First and the most important, the whole work is extremely simplified into two steps, one is the fast sampling procedure and the other is the trajectory analysis procedure. Second, it contains several powerful sampling methods for the simulation on graphics process unit, including our previous mixing replica exchange molecular dynamics method. The method combines the advantages of the biased and unbiased simulations. Finally, it extracts the dominant transition pathways automatically from the folding network by Markov state model. Users do not need to do the tedious intermediate steps by hand. To illustrate the usage of FSATOOL in practice, we perform one simulation for a RNA hairpin in explicit solvent. All the results are presented. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Haomiao Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Qiankun Gong
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Haozhe Zhang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, 430074, Hubei, China
| |
Collapse
|
27
|
Kumar N, Marx D. How do ribozymes accommodate additional water molecules upon hydrostatic compression deep into the kilobar pressure regime? Biophys Chem 2019; 252:106192. [DOI: 10.1016/j.bpc.2019.106192] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Revised: 05/23/2019] [Accepted: 05/23/2019] [Indexed: 12/19/2022]
|
28
|
Kumar N, Marx D. Mechanistic role of nucleobases in self-cleavage catalysis of hairpin ribozyme at ambient versus high-pressure conditions. Phys Chem Chem Phys 2019; 20:20886-20898. [PMID: 30067263 DOI: 10.1039/c8cp03142h] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Ribozymes catalyze the site-specific self-cleavage of intramolecular phosphodiester bonds. Initially thought to act as metalloenzymes, they are now known to be functional even in the absence of divalent metal ions and specific nucleobases directly participate in the self-cleavage reaction. Here, we use extensive replica exchange molecular dynamics simulations to probe the precise mechanistic role of nucleobases by simulating precatalytic reactant and active precursor states of a hairpin ribozyme along its reaction path at ambient as well as high-pressure conditions. The results provide novel key insights into the self-cleavage of ribozymes. We find that deprotonation of the hydroxyl group is crucial and might be the penultimate step to the self-cleavage. The G8 nucleobase is found to stabilize the activated precursor into inline arrangement for facile nucleophilic attack of the scissile phosphate only after deprotonation of the hydroxyl group. The protonated A38 nucleobase, in contrast, mainly acts a proton donor to the O5'-oxygen leaving group that eventually leads to the self-cleavage. Indeed, systematic high-pressure simulations of catalytically relevant states confirm these findings and, moreover, provide support to the role of ribozymes as piezophilic biocatalysts with regard to their relevance in early life under extreme conditions in the realm of RNA world hypothesis.
Collapse
Affiliation(s)
- Narendra Kumar
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | | |
Collapse
|
29
|
Kührová P, Mlýnský V, Zgarbová M, Krepl M, Bussi G, Best RB, Otyepka M, Šponer J, Banáš P. Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions. J Chem Theory Comput 2019; 15:3288-3305. [PMID: 30896943 PMCID: PMC7491206 DOI: 10.1021/acs.jctc.8b00955] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.
Collapse
Affiliation(s)
- Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Miroslav Krepl
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Robert B. Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
30
|
Cesari A, Bottaro S, Lindorff-Larsen K, Banáš P, Šponer J, Bussi G. Fitting Corrections to an RNA Force Field Using Experimental Data. J Chem Theory Comput 2019; 15:3425-3431. [DOI: 10.1021/acs.jctc.9b00206] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andrea Cesari
- Scuola Internazionale
Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory and Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory and Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46, Olomouc, Czech Republic
- Institute of Biophysics
of the Czech Academy of Sciences, Kralovopolska 135, Brno 612 65, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale
Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
31
|
Pinamonti G, Paul F, Noé F, Rodriguez A, Bussi G. The mechanism of RNA base fraying: Molecular dynamics simulations analyzed with core-set Markov state models. J Chem Phys 2019; 150:154123. [DOI: 10.1063/1.5083227] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Giovanni Pinamonti
- Department for Mathematics and Computer Science, Freie Universität, Berlin, Germany
| | - Fabian Paul
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, The University of Chicago, Chicago, Illinois 60637, USA
| | - Frank Noé
- Department for Mathematics and Computer Science, Freie Universität, Berlin, Germany
| | - Alex Rodriguez
- ICTP, International Centre for Theoretical Physics, Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste, Italy
| |
Collapse
|
32
|
Robertson MJ, Qian Y, Robinson MC, Tirado-Rives J, Jorgensen WL. Development and Testing of the OPLS-AA/M Force Field for RNA. J Chem Theory Comput 2019; 15:2734-2742. [PMID: 30807148 PMCID: PMC6585454 DOI: 10.1021/acs.jctc.9b00054] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Significant improvements have been made to the OPLS-AA force field for modeling RNA. New torsional potentials were optimized based on density functional theory (DFT) scans at the ωB97X-D/6-311++G(d,p) level for potential energy surfaces of the backbone α and γ dihedral angles. In combination with previously reported improvements for the sugar puckering and glycosidic torsion terms, the new force field was validated through diverse molecular dynamics simulations for RNAs in aqueous solution. Results for dinucleotides and tetranucleotides revealed both accurate reproduction of 3 J couplings from NMR and the avoidance of several unphysical states observed with other force fields. Simulations of larger systems with noncanonical motifs showed significant structural improvements over the previous OPLS-AA parameters. The new force field, OPLS-AA/M, is expected to perform competitively with other recent RNA force fields and to be compatible with OPLS-AA models for proteins and small molecules.
Collapse
Affiliation(s)
- Michael J. Robertson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Yue Qian
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Matthew C. Robinson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - William L. Jorgensen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
33
|
Chakraborty D, Wales DJ. Dynamics of an adenine-adenine RNA conformational switch from discrete path sampling. J Chem Phys 2019; 150:125101. [DOI: 10.1063/1.5070152] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, USA
| | - David J. Wales
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
34
|
Bottaro S, Bussi G, Pinamonti G, Reißer S, Boomsma W, Lindorff-Larsen K. Barnaba: software for analysis of nucleic acid structures and trajectories. RNA (NEW YORK, N.Y.) 2019; 25:219-231. [PMID: 30420522 PMCID: PMC6348988 DOI: 10.1261/rna.067678.118] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/26/2018] [Accepted: 11/06/2018] [Indexed: 06/09/2023]
Abstract
RNA molecules are highly dynamic systems characterized by a complex interplay between sequence, structure, dynamics, and function. Molecular simulations can potentially provide powerful insights into the nature of these relationships. The analysis of structures and molecular trajectories of nucleic acids can be nontrivial because it requires processing very high-dimensional data that are not easy to visualize and interpret. Here we introduce Barnaba, a Python library aimed at facilitating the analysis of nucleic acid structures and molecular simulations. The software consists of a variety of analysis tools that allow the user to (i) calculate distances between three-dimensional structures using different metrics, (ii) back-calculate experimental data from three-dimensional structures, (iii) perform cluster analysis and dimensionality reductions, (iv) search three-dimensional motifs in PDB structures and trajectories, and (v) construct elastic network models for nucleic acids and nucleic acids-protein complexes. In addition, Barnaba makes it possible to calculate torsion angles, pucker conformations, and to detect base-pairing/base-stacking interactions. Barnaba produces graphics that conveniently visualize both extended secondary structure and dynamics for a set of molecular conformations. The software is available as a command-line tool as well as a library, and supports a variety of file formats such as PDB, dcd, and xtc files. Source code, documentation, and examples are freely available at https://github.com/srnas/barnaba under GNU GPLv3 license.
Collapse
Affiliation(s)
- Sandro Bottaro
- Structural Biology and NMR Laboratory and Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen 2200, Denmark
- International School for Advanced Studies, 34136 Trieste, Italy
| | - Giovanni Bussi
- International School for Advanced Studies, 34136 Trieste, Italy
| | - Giovanni Pinamonti
- International School for Advanced Studies, 34136 Trieste, Italy
- Department of Mathematics and Computer Science, Freie Universität, 14195 Berlin, Germany
| | - Sabine Reißer
- International School for Advanced Studies, 34136 Trieste, Italy
| | - Wouter Boomsma
- Department of Computer Science, University of Copenhagen, Copenhagen 2200, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory and Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen 2200, Denmark
| |
Collapse
|
35
|
Dans PD, Gallego D, Balaceanu A, Darré L, Gómez H, Orozco M. Modeling, Simulations, and Bioinformatics at the Service of RNA Structure. Chem 2019. [DOI: 10.1016/j.chempr.2018.09.015] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
36
|
Abstract
This chapter discusses how the PLUMED plugin for molecular dynamics can be used to analyze and bias molecular dynamics trajectories. The chapter begins by introducing the notion of a collective variable and by then explaining how the free energy can be computed as a function of one or more collective variables. A number of practical issues mostly around periodic boundary conditions that arise when these types of calculations are performed using PLUMED are then discussed. Later parts of the chapter discuss how PLUMED can be used to perform enhanced sampling simulations that introduce simulation biases or multiple replicas of the system and Monte Carlo exchanges between these replicas. This section is then followed by a discussion on how free-energy surfaces and associated error bars can be extracted from such simulations by using weighted histogram and block averaging techniques.
Collapse
Affiliation(s)
- Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Trieste, Italy.
| | - Gareth A Tribello
- Atomistic Simulation Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast, UK.
| |
Collapse
|
37
|
Lemkul JA, MacKerell AD. Polarizable force field for RNA based on the classical drude oscillator. J Comput Chem 2018; 39:2624-2646. [PMID: 30515902 PMCID: PMC6284239 DOI: 10.1002/jcc.25709] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Revised: 08/01/2018] [Accepted: 09/23/2018] [Indexed: 12/15/2022]
Abstract
RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
38
|
Smith LG, Tan Z, Spasic A, Dutta D, Salas-Estrada LA, Grossfield A, Mathews DH. Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations. J Chem Theory Comput 2018; 14:6598-6612. [PMID: 30375860 DOI: 10.1021/acs.jctc.8b00633] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χOL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.
Collapse
Affiliation(s)
- Louis G Smith
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Zhen Tan
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Aleksandar Spasic
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Debapratim Dutta
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Leslie A Salas-Estrada
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - Alan Grossfield
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - David H Mathews
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Department of Biostatistics and Computational Biology , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| |
Collapse
|
39
|
Bottaro S, Bussi G, Kennedy SD, Turner DH, Lindorff-Larsen K. Conformational ensembles of RNA oligonucleotides from integrating NMR and molecular simulations. SCIENCE ADVANCES 2018; 4:eaar8521. [PMID: 29795785 PMCID: PMC5959319 DOI: 10.1126/sciadv.aar8521] [Citation(s) in RCA: 79] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2017] [Accepted: 04/05/2018] [Indexed: 05/08/2023]
Abstract
RNA molecules are key players in numerous cellular processes and are characterized by a complex relationship between structure, dynamics, and function. Despite their apparent simplicity, RNA oligonucleotides are very flexible molecules, and understanding their internal dynamics is particularly challenging using experimental data alone. We show how to reconstruct the conformational ensemble of four RNA tetranucleotides by combining atomistic molecular dynamics simulations with nuclear magnetic resonance spectroscopy data. The goal is achieved by reweighting simulations using a maximum entropy/Bayesian approach. In this way, we overcome problems of current simulation methods, as well as in interpreting ensemble- and time-averaged experimental data. We determine the populations of different conformational states by considering several nuclear magnetic resonance parameters and point toward properties that are not captured by state-of-the-art molecular force fields. Although our approach is applied on a set of model systems, it is fully general and may be used to study the conformational dynamics of flexible biomolecules and to detect inaccuracies in molecular dynamics force fields.
Collapse
Affiliation(s)
- Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
- Corresponding author. (S.B.); (K.L.-L.)
| | - Giovanni Bussi
- International School for Advanced Studies, Trieste, Italy
| | - Scott D. Kennedy
- Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, Rochester, NY 14642, USA
| | - Douglas H. Turner
- Department of Chemistry, University of Rochester, Rochester, NY 14627, USA
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
- Corresponding author. (S.B.); (K.L.-L.)
| |
Collapse
|
40
|
Š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: 336] [Impact Index Per Article: 56.0] [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
|
41
|
Zhang C, Lu C, Jing Z, Wu C, Piquemal JP, Ponder JW, Ren P. AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids. J Chem Theory Comput 2018; 14:2084-2108. [PMID: 29438622 PMCID: PMC5893433 DOI: 10.1021/acs.jctc.7b01169] [Citation(s) in RCA: 152] [Impact Index Per Article: 25.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.
Collapse
Affiliation(s)
- Changsheng Zhang
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chao Lu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Zhifeng Jing
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chuanjie Wu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
42
|
Bottaro S, Banáš P, Šponer J, Bussi G. Correction to "Free Energy Landscape of GAGA and UUCG RNA Tetraloops". J Phys Chem Lett 2018; 9:1674-1675. [PMID: 29553740 DOI: 10.1021/acs.jpclett.8b00754] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
|
43
|
Effects and limitations of a nucleobase-driven backmapping procedure for nucleic acids using steered molecular dynamics. Biochem Biophys Res Commun 2018; 498:352-358. [DOI: 10.1016/j.bbrc.2017.12.057] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2017] [Revised: 11/13/2017] [Accepted: 12/11/2017] [Indexed: 11/19/2022]
|
44
|
Poblete S, Bottaro S, Bussi G. A nucleobase-centered coarse-grained representation for structure prediction of RNA motifs. Nucleic Acids Res 2018; 46:1674-1683. [PMID: 29272539 PMCID: PMC5829650 DOI: 10.1093/nar/gkx1269] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Revised: 12/05/2017] [Accepted: 12/07/2017] [Indexed: 01/30/2023] Open
Abstract
We introduce the SPlit-and-conQueR (SPQR) model, a coarse-grained (CG) representation of RNA designed for structure prediction and refinement. In our approach, the representation of a nucleotide consists of a point particle for the phosphate group and an anisotropic particle for the nucleoside. The interactions are, in principle, knowledge-based potentials inspired by the $\mathcal {E}$SCORE function, a base-centered scoring function. However, a special treatment is given to base-pairing interactions and certain geometrical conformations which are lost in a raw knowledge-based model. This results in a representation able to describe planar canonical and non-canonical base pairs and base-phosphate interactions and to distinguish sugar puckers and glycosidic torsion conformations. The model is applied to the folding of several structures, including duplexes with internal loops of non-canonical base pairs, tetraloops, junctions and a pseudoknot. For the majority of these systems, experimental structures are correctly predicted at the level of individual contacts. We also propose a method for efficiently reintroducing atomistic detail from the CG representation.
Collapse
Affiliation(s)
- Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| | - Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
45
|
Smith LG, Zhao J, Mathews DH, Turner DH. Physics-based all-atom modeling of RNA energetics and structure. WILEY INTERDISCIPLINARY REVIEWS-RNA 2018; 8. [PMID: 28815951 DOI: 10.1002/wrna.1422] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Revised: 02/03/2017] [Accepted: 03/08/2017] [Indexed: 12/31/2022]
Abstract
The database of RNA sequences is exploding, but knowledge of energetics, structures, and dynamics lags behind. All-atom computational methods, such as molecular dynamics, hold promise for closing this gap. New algorithms and faster computers have accelerated progress in improving the reliability and accuracy of predictions. Currently, the methods can facilitate refinement of experimentally determined nuclear magnetic resonance and x-ray structures, but are 'unreliable' for predictions based only on sequence. Much remains to be discovered, however, about the many molecular interactions driving RNA folding and the best way to approximate them quantitatively. The large number of parameters required means that a wide variety of experimental results will be required to benchmark force fields and different approaches. As computational methods become more reliable and accessible, they will be used by an increasing number of biologists, much as x-ray crystallography has expanded. Thus, many fundamental physical principles underlying the computational methods are described. This review presents a summary of the current state of molecular dynamics as applied to RNA. It is designed to be helpful to students, postdoctoral fellows, and faculty who are considering or starting computational studies of RNA. WIREs RNA 2017, 8:e1422. doi: 10.1002/wrna.1422.
Collapse
Affiliation(s)
- Louis G Smith
- Department of Biochemistry and Biophysics and Center for RNA Biology, School of Medicine and Dentistry, University of Rochester, Rochester, NY, USA
| | - Jianbo Zhao
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, NY, USA
| | - David H Mathews
- Department of Biochemistry and Biophysics and Center for RNA Biology, School of Medicine and Dentistry, University of Rochester, Rochester, NY, USA
| | - Douglas H Turner
- Department of Chemistry and Center for RNA Biology, University of Rochester, Rochester, NY, USA
| |
Collapse
|
46
|
Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments. COMPUTATION 2018. [DOI: 10.3390/computation6010015] [Citation(s) in RCA: 68] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
47
|
Mlýnský V, Bussi G. Exploring RNA structure and dynamics through enhanced sampling simulations. Curr Opin Struct Biol 2018; 49:63-71. [PMID: 29414513 DOI: 10.1016/j.sbi.2018.01.004] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2017] [Revised: 12/20/2017] [Accepted: 01/02/2018] [Indexed: 01/23/2023]
Abstract
RNA function is intimately related to its structural dynamics. Molecular dynamics simulations are useful for exploring biomolecular flexibility but are severely limited by the accessible timescale. Enhanced sampling methods allow this timescale to be effectively extended in order to probe biologically relevant conformational changes and chemical reactions. Here, we review the role of enhanced sampling techniques in the study of RNA systems. We discuss the challenges and promises associated with the application of these methods to force-field validation, exploration of conformational landscapes and ion/ligand-RNA interactions, as well as catalytic pathways. Important technical aspects of these methods, such as the choice of the biased collective variables and the analysis of multi-replica simulations, are examined in detail. Finally, a perspective on the role of these methods in the characterization of RNA dynamics is provided.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy.
| |
Collapse
|
48
|
Mlýnský V, Bussi G. Molecular Dynamics Simulations Reveal an Interplay between SHAPE Reagent Binding and RNA Flexibility. J Phys Chem Lett 2018; 9:313-318. [PMID: 29265824 PMCID: PMC5830694 DOI: 10.1021/acs.jpclett.7b02921] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2017] [Accepted: 12/21/2017] [Indexed: 05/10/2023]
Abstract
The function of RNA molecules usually depends on their overall fold and on the presence of specific structural motifs. Chemical probing methods are routinely used in combination with nearest-neighbor models to determine RNA secondary structure. Among the available methods, SHAPE is relevant due to its capability to probe all RNA nucleotides and the possibility to be used in vivo. However, the structural determinants for SHAPE reactivity and its mechanism of reaction are still unclear. Here molecular dynamics simulations and enhanced sampling techniques are used to predict the accessibility of nucleotide analogs and larger RNA structural motifs to SHAPE reagents. We show that local RNA reconformations are crucial in allowing reagents to reach the 2'-OH group of a particular nucleotide and that sugar pucker is a major structural factor influencing SHAPE reactivity.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Scuola Internazionale Superiore di
Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di
Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
49
|
Yang C, Kulkarni M, Lim M, Pak Y. Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level. Nucleic Acids Res 2017; 45:12648-12656. [PMID: 29112755 PMCID: PMC5728390 DOI: 10.1093/nar/gkx1079] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Revised: 10/13/2017] [Accepted: 10/20/2017] [Indexed: 12/25/2022] Open
Abstract
The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.
Collapse
Affiliation(s)
- Changwon Yang
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Mandar Kulkarni
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Manho Lim
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| |
Collapse
|
50
|
Palazzesi F, Valsson O, Parrinello M. Conformational Entropy as Collective Variable for Proteins. J Phys Chem Lett 2017; 8:4752-4756. [PMID: 28906117 DOI: 10.1021/acs.jpclett.7b01770] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Many enhanced sampling methods rely on the identification of appropriate collective variables. For proteins, even small ones, finding appropriate descriptors has proven challenging. Here we suggest that the NMR S2 order parameter can be used to this effect. We trace the validity of this statement to the suggested relation between S2 and conformational entropy. Using the S2 order parameter and a surrogate for the protein enthalpy in conjunction with metadynamics or variationally enhanced sampling, we are able to reversibly fold and unfold a small protein and draw its free energy at a fraction of the time that is needed in unbiased simulations. We also use S2 in combination with the free energy flooding method to compute the unfolding rate of this peptide. We repeat this calculation at different temperatures to obtain the unfolding activation energy.
Collapse
Affiliation(s)
- Ferruccio Palazzesi
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| | - Omar Valsson
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| |
Collapse
|