1
|
Sarkar R, Mainan A, Roy S. Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods. Chem Commun (Camb) 2024. [PMID: 38501190 DOI: 10.1039/d3cc06105a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.
Collapse
Affiliation(s)
- Raju Sarkar
- Department of Chemical Sciences, Indian Institute of Science Education and Research, Kolkata, West Bengal 741246, India.
| | - Avijit Mainan
- Department of Chemical Sciences, Indian Institute of Science Education and Research, Kolkata, West Bengal 741246, India.
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research, Kolkata, West Bengal 741246, India.
| |
Collapse
|
2
|
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
|
3
|
Du B, Tian P. Factorization in molecular modeling and belief propagation algorithms. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2023; 20:21147-21162. [PMID: 38124591 DOI: 10.3934/mbe.2023935] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Factorization reduces computational complexity, and is therefore an important tool in statistical machine learning of high dimensional systems. Conventional molecular modeling, including molecular dynamics and Monte Carlo simulations of molecular systems, is a large research field based on approximate factorization of molecular interactions. Recently, the local distribution theory was proposed to factorize joint distribution of a given molecular system into trainable local distributions. Belief propagation algorithms are a family of exact factorization algorithms for (junction) trees, and are extended to approximate loopy belief propagation algorithms for graphs with loops. Despite the fact that factorization of probability distribution is the common foundation, computational research in molecular systems and machine learning studies utilizing belief propagation algorithms have been carried out independently with respective track of algorithm development. The connection and differences among these factorization algorithms are briefly presented in this perspective, with the hope to intrigue further development of factorization algorithms for physical modeling of complex molecular systems.
Collapse
Affiliation(s)
- Bochuan Du
- School of Life Sciences, Jilin University, Changchun 130012, China
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun 130012, China
- School of Artificial Intelligence, Jilin University, Changchun 130012, China
| |
Collapse
|
4
|
Dodaro A, Pavan M, Menin S, Salmaso V, Sturlese M, Moro S. Thermal titration molecular dynamics (TTMD): shedding light on the stability of RNA-small molecule complexes. Front Mol Biosci 2023; 10:1294543. [PMID: 38028536 PMCID: PMC10679717 DOI: 10.3389/fmolb.2023.1294543] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 10/20/2023] [Indexed: 12/01/2023] Open
Abstract
Ribonucleic acids are gradually becoming relevant players among putative drug targets, thanks to the increasing amount of structural data exploitable for the rational design of selective and potent binders that can modulate their activity. Mainly, this information allows employing different computational techniques for predicting how well would a ribonucleic-targeting agent fit within the active site of its target macromolecule. Due to some intrinsic peculiarities of complexes involving nucleic acids, such as structural plasticity, surface charge distribution, and solvent-mediated interactions, the application of routinely adopted methodologies like molecular docking is challenged by scoring inaccuracies, while more physically rigorous methods such as molecular dynamics require long simulation times which hamper their conformational sampling capabilities. In the present work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of unbinding kinetics, to characterize RNA-ligand complexes. In this article, we explored its applicability as a post-docking refinement tool on RNA in complex with small molecules, highlighting the capability of this method to identify the native binding mode among a set of decoys across various pharmaceutically relevant test cases.
Collapse
Affiliation(s)
| | | | | | | | | | - Stefano Moro
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy
| |
Collapse
|
5
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
6
|
Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions. Comput Struct Biotechnol J 2022; 20:6120-6137. [PMID: 36420155 PMCID: PMC9672130 DOI: 10.1016/j.csbj.2022.11.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Revised: 11/04/2022] [Accepted: 11/04/2022] [Indexed: 11/13/2022] Open
Abstract
The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.
Collapse
|
7
|
Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts. J Mol Biol 2022; 434:167788. [PMID: 35963460 PMCID: PMC10042644 DOI: 10.1016/j.jmb.2022.167788] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Revised: 08/03/2022] [Accepted: 08/07/2022] [Indexed: 12/24/2022]
Abstract
Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (QTC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.
Collapse
|
8
|
Wang Y, Parmar S, Schneekloth JS, Tiwary P. Interrogating RNA-Small Molecule Interactions with Structure Probing and Artificial Intelligence-Augmented Molecular Simulations. ACS CENTRAL SCIENCE 2022; 8:741-748. [PMID: 35756372 PMCID: PMC9228567 DOI: 10.1021/acscentsci.2c00149] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Indexed: 05/10/2023]
Abstract
While there is increasing interest in the study of RNA as a therapeutic target, efforts to understand RNA-ligand recognition at the molecular level lag far behind our understanding of protein-ligand recognition. This problem is complicated due to the more than 10 orders of magnitude in time scales involved in RNA dynamics and ligand binding events, making it not straightforward to design experiments or simulations. Here, we make use of artificial intelligence (AI)-augmented molecular dynamics simulations to directly observe ligand dissociation for cognate and synthetic ligands from a riboswitch system. The site-specific flexibility profiles from our simulations are compared with in vitro measurements of flexibility using selective 2' hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP). Our simulations reproduce known relative binding affinity profiles for the cognate and synthetic ligands, and pinpoint how both ligands make use of different aspects of riboswitch flexibility. On the basis of our dissociation trajectories, we also make and validate predictions of pairs of mutations for both the ligand systems that would show differing binding affinities. These mutations are distal to the binding site and could not have been predicted solely on the basis of structure. The methodology demonstrated here shows how molecular dynamics simulations with all-atom force-fields have now come of age in making predictions that complement existing experimental techniques and illuminate aspects of systems otherwise not trivial to understand.
Collapse
Affiliation(s)
- Yihang Wang
- Biophysics
Program and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Shaifaly Parmar
- Chemical
Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702, United States
| | - John S. Schneekloth
- Chemical
Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702, United States
| | - Pratyush Tiwary
- Department
of Chemistry and Biochemistry and Institute for Physical Science and
Technology, University of Maryland, College Park 20742, United States
| |
Collapse
|
9
|
Lam K, Kasavajhala K, Gunasekera S, Simmerling C. Accelerating the Ensemble Convergence of RNA Hairpin Simulations with a Replica Exchange Structure Reservoir. J Chem Theory Comput 2022; 18:3930-3947. [PMID: 35502992 PMCID: PMC10658646 DOI: 10.1021/acs.jctc.2c00065] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
RNA is a key participant in many biological processes, but studies of RNA using computer simulations lag behind those of proteins, largely due to less-developed force fields and the slow dynamics of RNA. Generating converged RNA ensembles for force field development and other studies remains a challenge. In this study, we explore the ability of replica exchange molecular dynamics to obtain well-converged conformational ensembles for two RNA hairpin systems in an implicit solvent. Even for these small model systems, standard REMD remains computationally costly, but coupling to a pre-generated structure library using the reservoir REMD approach provides a dramatic acceleration of ensemble convergence for both model systems. Such precise ensembles could facilitate RNA force field development and validation and applications of simulation to more complex RNA systems. The advantages and remaining challenges of applying R-REMD to RNA are investigated in detail.
Collapse
Affiliation(s)
- Kenneth Lam
- Molecular and Cellular Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Koushik Kasavajhala
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Sarah Gunasekera
- Program in Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Carlos Simmerling
- Molecular and Cellular Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods. QRB DISCOVERY 2022. [PMID: 37529277 PMCID: PMC10392686 DOI: 10.1017/qrd.2022.19] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Abstract
While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.
Collapse
|
13
|
Wang J, Lan L, Wu X, Xu L, Miao Y. Mechanism of RNA recognition by a Musashi RNA-binding protein. Curr Res Struct Biol 2021; 4:10-20. [PMID: 34988468 PMCID: PMC8695263 DOI: 10.1016/j.crstbi.2021.12.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 10/31/2021] [Accepted: 12/07/2021] [Indexed: 12/21/2022] Open
Abstract
The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Lan Lan
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Xiaoqing Wu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Liang Xu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
- Department of Radiation Oncology, The University of Kansas Cancer Center, Kansas City, KS, 66160, USA
| | - Yinglong Miao
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|
14
|
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
|
15
|
Krepl M, Damberger FF, von Schroetter C, Theler D, Pokorná P, Allain FHT, Šponer J. Recognition of N6-Methyladenosine by the YTHDC1 YTH Domain Studied by Molecular Dynamics and NMR Spectroscopy: The Role of Hydration. J Phys Chem B 2021; 125:7691-7705. [PMID: 34258996 DOI: 10.1021/acs.jpcb.1c03541] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m6A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m6A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m6A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm6A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m6A binding pocket. The 5'-Im6A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m6A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m6A recognition.
Collapse
Affiliation(s)
- Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Fred Franz Damberger
- Department of Biology, Institute of Biochemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | | - Dominik Theler
- Department of Biology, Institute of Biochemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Frédéric H-T Allain
- Department of Biology, Institute of Biochemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Olomouc 783 71, Czech Republic
| |
Collapse
|
16
|
Role and Perspective of Molecular Simulation-Based Investigation of RNA-Ligand Interaction: From Small Molecules and Peptides to Photoswitchable RNA Binding. Molecules 2021; 26:molecules26113384. [PMID: 34205049 PMCID: PMC8199858 DOI: 10.3390/molecules26113384] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 05/28/2021] [Accepted: 06/01/2021] [Indexed: 12/15/2022] Open
Abstract
Aberrant RNA–protein complexes are formed in a variety of diseases. Identifying the ligands that interfere with their formation is a valuable therapeutic strategy. Molecular simulation, validated against experimental data, has recently emerged as a powerful tool to predict both the pose and energetics of such ligands. Thus, the use of molecular simulation may provide insight into aberrant molecular interactions in diseases and, from a drug design perspective, may allow for the employment of less wet lab resources than traditional in vitro compound screening approaches. With regard to basic research questions, molecular simulation can support the understanding of the exact molecular interaction and binding mode. Here, we focus on examples targeting RNA–protein complexes in neurodegenerative diseases and viral infections. These examples illustrate that the strategy is rather general and could be applied to different pharmacologically relevant approaches. We close this study by outlining one of these approaches, namely the light-controllable association of small molecules with RNA, as an emerging approach in RNA-targeting therapy.
Collapse
|
17
|
Cao X, Tian P. "Dividing and Conquering" and "Caching" in Molecular Modeling. Int J Mol Sci 2021; 22:5053. [PMID: 34068835 PMCID: PMC8126232 DOI: 10.3390/ijms22095053] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 11/17/2022] Open
Abstract
Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.
Collapse
Affiliation(s)
- Xiaoyong Cao
- School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun 130012, China;
- School of Artificial Intelligence, Jilin University, Changchun 130012, China
| |
Collapse
|
18
|
Identification of a Novel Cis-Acting Regulator of HIV-1 Genome Packaging. Int J Mol Sci 2021; 22:ijms22073435. [PMID: 33810482 PMCID: PMC8036536 DOI: 10.3390/ijms22073435] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 03/19/2021] [Accepted: 03/24/2021] [Indexed: 12/17/2022] Open
Abstract
Human immunodeficiency virus type 1 (HIV-1) uptakes homo-dimerized viral RNA genome into its own particle. A cis-acting viral RNA segment responsible for this event, termed packaging signal (psi), is located at the 5′-end of the viral genome. Although the psi segment exhibits nucleotide variation in nature, its effects on the psi function largely remain unknown. Here we show that a psi sequence from an HIV-1 regional variant, subtype D, has a lower packaging ability compared with that from another regional variant, HIV-1 subtype B, despite maintaining similar genome dimerization activities. A series of molecular genetic investigations narrowed down the responsible element of the selective attenuation to the two sequential nucleotides at positions 226 and 227 in the psi segment. Molecular dynamics simulations predicted that the dinucleotide substitution alters structural dynamics, fold, and hydrogen-bond networks primarily of the psi-SL2 element that contains the binding interface of viral nucleocapsid protein for the genome packaging. In contrast, such structural changes were minimal within the SL1 element involved in genome dimerization. These results suggest that the psi 226/227 dinucleotide pair functions as a cis-acting regulator to control the psi structure to selectively tune the efficiency of packaging, but not dimerization of highly variable HIV-1 genomes.
Collapse
|
19
|
Hurst T, Chen SJ. Deciphering nucleotide modification-induced structure and stability changes. RNA Biol 2021; 18:1920-1930. [PMID: 33586616 DOI: 10.1080/15476286.2021.1882179] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Nucleotide modification in RNA controls a bevy of biological processes, including RNA degradation, gene expression, and gene editing. In turn, misregulation of modified nucleotides is associated with a host of chronic diseases and disorders. However, the molecular mechanisms driving these processes remain poorly understood. To partially address this knowledge gap, we used alchemical and temperature replica exchange molecular dynamics (TREMD) simulations on an RNA duplex and an analogous hairpin to probe the structural effects of modified and/or mutant nucleotides. The simulations successfully predict the modification/mutation-induced relative free energy change for complementary duplex formation, and structural analyses highlight mechanisms driving stability changes. Furthermore, TREMD simulations for a hairpin-forming RNA with and without modification provide reliable estimations of the energy landscape. Illuminating the impact of methylated and/or mutated nucleotides on the structure-function relationship and the folding energy landscape, the simulations provide insights into modification-induced alterations to the folding mechanics of the hairpin. The results here may be biologically significant as hairpins are widespread structure motifs that play critical roles in gene expression and regulation. Specifically, the tetraloop of the probed hairpin is phylogenetically abundant, and the stem mirrors a miRNA seed region whose modification has been implicated in epilepsy pathogenesis.
Collapse
Affiliation(s)
- Travis Hurst
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, MO, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Institute for Data Science and Informatics, University of Missouri, Columbia, MO, USA
| |
Collapse
|
20
|
Nicholson DA, Sengupta A, Nesbitt DJ. Chirality-Dependent Amino Acid Modulation of RNA Folding. J Phys Chem B 2020; 124:11561-11572. [PMID: 33296203 DOI: 10.1021/acs.jpcb.0c07420] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The preponderance of a specific d- or l-chirality in fats, sugars, amino acids, nucleic acids, and so on is ubiquitous in nature, yet the biological origin of such chiral dominance (i.e., with one enantiomer overwhelmingly present) remains an open question. One plausible proposal for the predominance of l-chirality in amino acids could be through evolutionary templating of chiral RNA-folding via chaperone activity. To help evaluate this possibility, single molecule fluorescence experiments have been performed that measure the chiral dependence of chaperone folding dynamics for the simple tetraloop-tetraloop receptor (TL-TLR) tertiary binding motif in the presence of a series of chiral amino acids. Specifically, d- vs l-arginine is found to accelerate the unfolding of this RNA motif in a chirally selective fashion, with temperature-dependent studies of the kinetics performed to extract free energy, enthalpy, and entropy landscapes for the underlying thermodynamics. Furthermore, all-atom molecular dynamics (MD) simulations are pursued to provide additional physical insight into this chiral sensitivity, which reveal enantiomer-specific sampling of nucleic acid surfaces by d- vs l-arginine and support a putative mechanism for chirally specific denaturation of RNA tertiary structure by arginine but not other amino acids.
Collapse
Affiliation(s)
- David A Nicholson
- JILA, National Institute of Standards and Technology and University of Colorado Boulder, Boulder, Colorado 80309 United States.,Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Abhigyan Sengupta
- Department of Physics, Technical University of Munich, Garching, Munich, Germany 85748
| | - David J Nesbitt
- JILA, National Institute of Standards and Technology and University of Colorado Boulder, Boulder, Colorado 80309 United States.,Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, United States.,Department of Physics, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
21
|
Kushwaha GS, Patra A, Bhavesh NS. Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins. Front Microbiol 2020; 11:575041. [PMID: 33224117 PMCID: PMC7674647 DOI: 10.3389/fmicb.2020.575041] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Accepted: 10/06/2020] [Indexed: 11/25/2022] Open
Abstract
(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.
Collapse
Affiliation(s)
- Gajraj Singh Kushwaha
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India.,KIIT Technology Business Incubator (KIIT-TBI), Kalinga Institute of Industrial Technology (KIIT) (Deemed to be University), Bhubaneswar, India
| | - Anupam Patra
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India
| | - Neel Sarovar Bhavesh
- Transcription Regulation Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India
| |
Collapse
|
22
|
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
|
23
|
Reißer S, Zucchelli S, Gustincich S, Bussi G. Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data. Nucleic Acids Res 2020; 48:1164-1174. [PMID: 31889193 PMCID: PMC7026608 DOI: 10.1093/nar/gkz1184] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Revised: 12/05/2019] [Accepted: 12/09/2019] [Indexed: 01/12/2023] Open
Abstract
Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.
Collapse
Affiliation(s)
- Sabine Reißer
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Silvia Zucchelli
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- Department of Health Sciences, Center for Autoimmune and Allergic Diseases (CAAD) and Interdisciplinary Research Center of Autoimmune Diseases (IRCAD), University of Piemonte Orientale, Novara, Italy
| | - Stefano Gustincich
- Central RNA Laboratory and Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia (IIT), 16163 Genova, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
24
|
Molecular dynamics simulations of G-quadruplexes: The basic principles and their application to folding and ligand binding. ANNUAL REPORTS IN MEDICINAL CHEMISTRY 2020. [DOI: 10.1016/bs.armc.2020.04.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
25
|
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: 76] [Impact Index Per Article: 15.2] [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
|
26
|
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
|
27
|
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
|
28
|
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]
|
29
|
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]
|