1
|
Muscat S, Martino G, Manigrasso J, Marcia M, De Vivo M. On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules. J Chem Theory Comput 2024. [PMID: 39150960 DOI: 10.1021/acs.jctc.4c00773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2024]
Abstract
RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.
Collapse
Affiliation(s)
- Stefano Muscat
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Gianfranco Martino
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Jacopo Manigrasso
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, 431 50 Mölndal, Sweden
| | - Marco Marcia
- European Molecular Biology Laboratory Grenoble, 71 Avenue des Martyrs, 38042 Grenoble, France
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, 751 23 Uppsala, Sweden
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
2
|
Wieruszewska J, Pawłowicz A, Połomska E, Pasternak K, Gdaniec Z, Andrałojć W. The 8-17 DNAzyme can operate in a single active structure regardless of metal ion cofactor. Nat Commun 2024; 15:4218. [PMID: 38760331 PMCID: PMC11101458 DOI: 10.1038/s41467-024-48638-x] [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] [Received: 04/13/2023] [Accepted: 05/09/2024] [Indexed: 05/19/2024] Open
Abstract
DNAzymes - synthetic enzymes made of DNA - have long attracted attention as RNA-targeting therapeutic agents. Yet, as of now, no DNAzyme-based drug has been approved, partially due to our lacking understanding of their molecular mode of action. In this work we report the solution structure of 8-17 DNAzyme bound to a Zn2+ ion solved through NMR spectroscopy. Surprisingly, it turned out to be very similar to the previously solved Pb2+-bound form (catalytic domain RMSD = 1.28 Å), despite a long-standing literature consensus that Pb2+ recruits a different DNAzyme fold than other metal ion cofactors. Our follow-up NMR investigations in the presence of other ions - Mg2+, Na+, and Pb2+ - suggest that at DNAzyme concentrations used in NMR all these ions induce a similar tertiary fold. Based on these findings, we propose a model for 8-17 DNAzyme interactions with metal ions postulating the existence of only a single catalytically-active structure, yet populated to a different extent depending on the metal ion cofactor. Our results provide structural information on the 8-17 DNAzyme in presence of non-Pb2+ cofactors, including the biologically relevant Mg2+ ion.
Collapse
Affiliation(s)
- Julia Wieruszewska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland
| | - Aleksandra Pawłowicz
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland
| | - Ewa Połomska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland
| | - Karol Pasternak
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland
| | - Zofia Gdaniec
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland
| | - Witold Andrałojć
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704, Poznań, Noskowskiego, 12/14, Poland.
| |
Collapse
|
3
|
Zhang Y, Zhang J, Wan H, Wu Z, Xu H, Zhang Z, Wang Y, Wang J. New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn 2+ and Mg 2. J Phys Chem Lett 2024; 15:2708-2714. [PMID: 38427973 DOI: 10.1021/acs.jpclett.3c03221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2024]
Abstract
CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.
Collapse
Affiliation(s)
- Yaoyao Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jing Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hengjia Wan
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Ziwei Wu
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Huangtao Xu
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
| | - Zhe Zhang
- Institute for Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yujuan Wang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- International Magnetobiology Frontier Research Center (iMFRC), Science Island, Hefei 230031, China
| | - Junfeng Wang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, Anhui 230031, China
- University of Science and Technology of China, Hefei, Anhui 230026, China
- School of Basic Medical Sciences, Anhui Medical University, Hefei, Anhui 230032, China
- International Magnetobiology Frontier Research Center (iMFRC), Science Island, Hefei 230031, China
| |
Collapse
|
4
|
Weissman B, Ekesan Ş, Lin HC, Gardezi S, Li NS, Giese TJ, McCarthy E, Harris ME, York DM, Piccirilli JA. Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis. J Am Chem Soc 2023; 145:2830-2839. [PMID: 36706353 PMCID: PMC10112047 DOI: 10.1021/jacs.2c10079] [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: 01/28/2023]
Abstract
Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.
Collapse
Affiliation(s)
- Benjamin Weissman
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| | - Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Hsuan-Chun Lin
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Shahbaz Gardezi
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Nan-Sheng Li
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Erika McCarthy
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Michael E Harris
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Joseph A Piccirilli
- Department of Chemistry and Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
5
|
Yoon S, Harris ME. Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis. Biochemistry 2021; 60:2810-2823. [PMID: 34495648 DOI: 10.1021/acs.biochem.1c00489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D2O to H2O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D2O ratio (n) versus pL versus kn/k0. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.
Collapse
Affiliation(s)
- Suhyun Yoon
- Department of Chemistry, University of Florida, Gainesville, Florida 32611-7200, United States
| | - Michael E Harris
- Department of Chemistry, University of Florida, Gainesville, Florida 32611-7200, United States
| |
Collapse
|
6
|
Zhai J, Wu Y, Xie X. Single‐Component
Chemical Nose with a Hemicyanine Probe for
Pattern‐Based
Discrimination of Metal Ions
†. CHINESE J CHEM 2021. [DOI: 10.1002/cjoc.202100013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Jingying Zhai
- Academy for Advanced Interdisciplinary Studies, Southern University of Science and Technology Shenzhen Guangdong 518055 China
| | - Yaotian Wu
- Department of Chemistry, Southern University of Science and Technology Shenzhen Guangdong 518055 China
| | - Xiaojiang Xie
- Department of Chemistry, Southern University of Science and Technology Shenzhen Guangdong 518055 China
| |
Collapse
|
7
|
Hoogstraten CG, Terrazas M, Aviñó A, White NA, Sumita M. Dynamics-Function Analysis in Catalytic RNA Using NMR Spin Relaxation and Conformationally Restricted Nucleotides. Methods Mol Biol 2021; 2167:183-202. [PMID: 32712921 DOI: 10.1007/978-1-0716-0716-9_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
A full understanding of biomolecular function requires an analysis of both the dynamic properties of the system of interest and the identification of those dynamics that are required for function. We describe NMR methods based on metabolically directed specific isotope labeling for the identification of molecular disorder and/or conformational transitions on the RNA backbone ribose groups. These analyses are complemented by the use of synthetic covalently modified nucleotides constrained to a single sugar pucker, which allow functional assessment of dynamics by selectively removing a minor conformer identified by NMR from the structural ensemble.
Collapse
Affiliation(s)
- Charles G Hoogstraten
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA.
| | - Montserrat Terrazas
- Institute for Advanced Chemistry of Catalonia (IQAC), Spanish Council for Scientific Research (CSIC), Barcelona, Spain.,Joint IRB-BSC Program in Computational Biology, The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Anna Aviñó
- Institute for Advanced Chemistry of Catalonia (IQAC), Spanish Council for Scientific Research (CSIC), Barcelona, Spain.,Networking Center on Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Barcelona, Spain
| | - Neil A White
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA.,Molecular, Cellular, and Developmental Biology, Yale University, New Haven, CT, USA
| | - Minako Sumita
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA.,Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, IL, USA
| |
Collapse
|
8
|
Nguyen HT, Thirumalai D. Charge Density of Cation Determines Inner versus Outer Shell Coordination to Phosphate in RNA. J Phys Chem B 2020; 124:4114-4122. [PMID: 32342689 DOI: 10.1021/acs.jpcb.0c02371] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Divalent cations are often required to fold RNA, which is a highly charged polyanion. Condensation of ions, such as Mg2+ or Ca2+, in the vicinity of RNA renormalizes the effective charges on the phosphate groups, thus minimizing the intra RNA electrostatic repulsion. The prevailing view is that divalent ions bind diffusively in a nonspecific manner. In sharp contrast, we arrive at the exact opposite conclusion using a theory for the interaction of ions with the phosphate groups using RISM theory in conjunction with simulations based on an accurate three-interaction-site RNA model. The divalent ions bind in a nucleotide-specific manner using either the inner (partially dehydrated) or outer (fully hydrated) shell coordination. The high charge density Mg2+ ion has a preference to bind to the outer shell, whereas the opposite is the case for Ca2+. Surprisingly, we find that bridging interactions, involving ions that are coordinated to two or more phosphate groups, play a crucial role in maintaining the integrity of the folded state. Their importance could become increasingly prominent as the size of the RNA increases. Because the modes of interaction of divalent ions with DNA are likely to be similar, we propose that specific inner and outer shell coordination could play a role in DNA condensation, and perhaps genome organization as well.
Collapse
Affiliation(s)
- Hung T Nguyen
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - D Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
9
|
Gaines CS, Piccirilli JA, York DM. The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design. RNA (NEW YORK, N.Y.) 2020; 26:111-125. [PMID: 31776179 PMCID: PMC6961537 DOI: 10.1261/rna.071894.119] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Accepted: 11/14/2019] [Indexed: 05/13/2023]
Abstract
We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.
Collapse
Affiliation(s)
- Colin S Gaines
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Joseph A Piccirilli
- Department of Biochemistry and Molecular Biology and Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
10
|
Maurel MC, Leclerc F, Hervé G. Ribozyme Chemistry: To Be or Not To Be under High Pressure. Chem Rev 2019; 120:4898-4918. [DOI: 10.1021/acs.chemrev.9b00457] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marie-Christine Maurel
- Institut de Systématique, Evolution, Biodiversité (ISYEB), CNRS, Sorbonne Université, Muséum National d’Histoire Naturelle, EPHE, F-75005 Paris, France
| | - Fabrice Leclerc
- Institute for Integrative Biology of the Cell (I2BC), CNRS, CEA, Université Paris Sud, F-91198 Gif-sur-Yvette, France
| | - Guy Hervé
- Laboratoire BIOSIPE, Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Campus Pierre et Marie Curie, F-75005 Paris, France
| |
Collapse
|
11
|
Ekesan Ş, York DM. Dynamical ensemble of the active state and transition state mimic for the RNA-cleaving 8-17 DNAzyme in solution. Nucleic Acids Res 2019; 47:10282-10295. [PMID: 31511899 PMCID: PMC6821293 DOI: 10.1093/nar/gkz773] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2019] [Revised: 08/20/2019] [Accepted: 09/03/2019] [Indexed: 02/01/2023] Open
Abstract
We perform molecular dynamics simulations, based on recent crystallographic data, on the 8-17 DNAzyme at four states along the reaction pathway to determine the dynamical ensemble for the active state and transition state mimic in solution. A striking finding is the diverse roles played by Na+ and Pb2+ ions in the electrostatically strained active site that impact all four fundamental catalytic strategies, and share commonality with some features recently inferred for naturally occurring hammerhead and pistol ribozymes. The active site Pb2+ ion helps to stabilize in-line nucleophilic attack, provides direct electrostatic transition state stabilization, and facilitates leaving group departure. A conserved guanine residue is positioned to act as the general base, and is assisted by a bridging Na+ ion that tunes the pKa and facilitates in-line fitness. The present work provides insight into how DNA molecules are able to solve the RNA-cleavage problem, and establishes functional relationships between the mechanism of these engineered DNA enzymes with their naturally evolved RNA counterparts. This adds valuable information to our growing body of knowledge on general mechanisms of phosphoryl transfer reactions catalyzed by RNA, proteins and DNA.
Collapse
Affiliation(s)
- Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
12
|
Kostenbader K, York DM. Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme. RNA (NEW YORK, N.Y.) 2019; 25:1439-1456. [PMID: 31363004 PMCID: PMC6795133 DOI: 10.1261/rna.071944.119] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/10/2019] [Accepted: 07/11/2019] [Indexed: 05/27/2023]
Abstract
The pistol ribozyme (Psr) is among the most recently discovered RNA enzymes and has been the subject of experiments aimed at elucidating the mechanism. Recent biochemical studies have revealed exciting clues about catalytic interactions in the active site not apparent from available crystallographic data. The present work unifies the interpretation of the existing body of structural and functional data on Psr by providing a dynamical model for the catalytically active state in solution from molecular simulation. Our results suggest that a catalytic Mg2+ ion makes inner-sphere contact with G33:N7 and outer-sphere coordination to the pro-RP of the scissile phosphate, promoting electrostatic stabilization of the dianionic transition state and neutralization of the developing charge of the leaving group through a metal-coordinated water molecule that is made more acidic by a hydrogen bond donated from the 2'OH of P32. This model is consistent with experimental activity-pH and mutagenesis data, including sensitivity to G33(7cG) and phosphorothioate substitution/metal ion rescue. The model suggests several experimentally testable predictions, including the response of cleavage activity to mutations at G42 and P32 positions in the ribozyme, and thio substitutions of the substrate in the presence of different divalent metal ions. Further, the model identifies striking similarities of Psr to the hammerhead ribozyme (HHr), including similar global fold, organization of secondary structure around an active site three-way junction, catalytic metal ion binding mode, and guanine general base. However, the specific binding mode and role of the Mg2+ ion, as well as a conserved 2'-OH in the active site, are interrelated but subtly different between the ribozymes.
Collapse
Affiliation(s)
- Ken Kostenbader
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry & Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8076, USA
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry & Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8076, USA
| |
Collapse
|
13
|
Hexahydrated Mg 2+ Binding and Outer-Shell Dehydration on RNA Surface. Biophys J 2019; 114:1274-1284. [PMID: 29590585 DOI: 10.1016/j.bpj.2018.01.040] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 01/30/2018] [Accepted: 01/31/2018] [Indexed: 10/17/2022] Open
Abstract
The interaction between metal ions, especially Mg2+ ions, and RNA plays a critical role in RNA folding. Upon binding to RNA, a metal ion that is fully hydrated in bulk solvent can become dehydrated. Here we use molecular dynamics simulation to investigate the dehydration of bound hexahydrated Mg2+ ions. We find that a hydrated Mg2+ ion in the RNA groove region can involve significant dehydration in the outer hydration shell. The first or innermost hydration shell of the Mg2+ ion, however, is retained during the simulation because of the strong ion-water electrostatic attraction. As a result, water-mediated hydrogen bonding remains an important form for Mg2+-RNA interaction. Analysis for ions at different binding sites shows that the most pronounced water deficiency relative to the fully hydrated state occurs at a radial distance of around 11 Å from the center of the ion. Based on the independent 200 ns molecular dynamics simulations for three different RNA structures (Protein Data Bank: 1TRA, 2TPK, and 437D), we find that Mg2+ ions overwhelmingly dominate over monovalent ions such as Na+ and K+ in ion-RNA binding. Furthermore, application of the free energy perturbation method leads to a quantitative relationship between the Mg2+ dehydration free energy and the local structural environment. We find that ΔΔGhyd, the change of the Mg2+ hydration free energy upon binding to RNA, varies linearly with the inverse distance between the Mg2+ ion and the nearby nonbridging oxygen atoms of the phosphate groups, and ΔΔGhyd can reach -2.0 kcal/mol and -3.0 kcal/mol for an Mg2+ ion bound to the surface and to the groove interior, respectively. In addition, the computation results in an analytical formula for the hydration ratio as a function of the average inverse Mg2+-O distance. The results here might be useful for further quantitative investigations of ion-RNA interactions in RNA folding.
Collapse
|
14
|
Giambasu GM, Case DA, York DM. Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory. J Am Chem Soc 2019; 141:2435-2445. [PMID: 30632365 PMCID: PMC6574206 DOI: 10.1021/jacs.8b11474] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.
Collapse
Affiliation(s)
- George M. Giambasu
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| | - David A. Case
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| | - Darrin M. York
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
- Laboratory for Biomolecular Simulation Research, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
- Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| |
Collapse
|
15
|
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]
|
16
|
Sun LZ, Chen SJ. Predicting RNA-Metal Ion Binding with Ion Dehydration Effects. Biophys J 2018; 116:184-195. [PMID: 30612712 DOI: 10.1016/j.bpj.2018.12.006] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Revised: 11/30/2018] [Accepted: 12/07/2018] [Indexed: 01/02/2023] Open
Abstract
Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg2+ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou, China; Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri.
| |
Collapse
|
17
|
White NA, Sumita M, Marquez VE, Hoogstraten CG. Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme. RNA (NEW YORK, N.Y.) 2018; 24:1542-1554. [PMID: 30111534 PMCID: PMC6191710 DOI: 10.1261/rna.067579.118] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 08/02/2018] [Indexed: 06/08/2023]
Abstract
In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.
Collapse
Affiliation(s)
- Neil A White
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
| | - Minako Sumita
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
| | - Victor E Marquez
- Chemical Biology Laboratory, Molecular Discovery Program, Center for Cancer Research, National Cancer Institute at Frederick, National Institutes of Health, Frederick, Maryland 21702, USA
| | - Charles G Hoogstraten
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
18
|
Š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
|
19
|
Šponer JE, Szabla R, Góra RW, Saitta AM, Pietrucci F, Saija F, Di Mauro E, Saladino R, Ferus M, Civiš S, Šponer J. Prebiotic synthesis of nucleic acids and their building blocks at the atomic level - merging models and mechanisms from advanced computations and experiments. Phys Chem Chem Phys 2018; 18:20047-66. [PMID: 27136968 DOI: 10.1039/c6cp00670a] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The origin of life on Earth is one of the most fascinating questions of contemporary science. Extensive research in the past decades furnished diverse experimental proposals for the emergence of first informational polymers that could form the basis of the early terrestrial life. Side by side with the experiments, the fast development of modern computational chemistry methods during the last 20 years facilitated the use of in silico modelling tools to complement the experiments. Modern computations can provide unique atomic-level insights into the structural and electronic aspects as well as the energetics of key prebiotic chemical reactions. Many of these insights are not directly obtainable from the experimental techniques and the computations are thus becoming indispensable for proper interpretation of many experiments and for qualified predictions. This review illustrates the synergy between experiment and theory in the origin of life research focusing on the prebiotic synthesis of various nucleic acid building blocks and on the self-assembly of nucleotides leading to the first functional oligonucleotides.
Collapse
Affiliation(s)
- Judit E Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| | - Rafał Szabla
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic.
| | - Robert W Góra
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - A Marco Saitta
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Fabio Pietrucci
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Franz Saija
- CNR-IPCF, Viale Ferdinando Stagno d'Alcontres 37, 98158 Messina, Italy
| | - Ernesto Di Mauro
- Dipartimento di Biologia e Biotecnologie "Charles Darwin", "Sapienza" Università di Roma, Piazzale Aldo Moro 5, Rome 00185, Italy
| | - Raffaele Saladino
- Dipartimento di Scienze Ecologiche e Biologiche Università della Tuscia, Via San Camillo De Lellis, 01100 Viterbo, Italy
| | - Martin Ferus
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Svatopluk Civiš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| |
Collapse
|
20
|
Patra P, Ghosh M, Banerjee R, Chakrabarti J. Anion induced conformational preference of C α NN motif residues in functional proteins. Proteins 2017; 85:2179-2190. [PMID: 28905427 DOI: 10.1002/prot.25382] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2017] [Revised: 08/16/2017] [Accepted: 09/10/2017] [Indexed: 11/12/2022]
Abstract
Among different ligand binding motifs, anion binding Cα NN motif consisting of peptide backbone atoms of three consecutive residues are observed to be important for recognition of free anions, like sulphate or biphosphate and participate in different key functions. Here we study the interaction of sulphate and biphosphate with Cα NN motif present in different proteins. Instead of total protein, a peptide fragment has been studied keeping Cα NN motif flanked in between other residues. We use classical force field based molecular dynamics simulations to understand the stability of this motif. Our data indicate fluctuations in conformational preferences of the motif residues in absence of the anion. The anion gives stability to one of these conformations. However, the anion induced conformational preferences are highly sequence dependent and specific to the type of anion. In particular, the polar residues are more favourable compared to the other residues for recognising the anion.
Collapse
Affiliation(s)
- Piya Patra
- Maulana Abul Kalam Azad University of Technology, West Bengal, (Formerly known as WBUT), BF-142, Sector-I, Saltlake, Kolkata, 700 064, India
| | - Mahua Ghosh
- Department of Chemical, Biological and Macro-Molecular Sciences, S.N. Bose National Centre for Basic Sciences, Sector III, Block JD, Salt Lake, Kolkata, 700106, India
| | - Raja Banerjee
- Maulana Abul Kalam Azad University of Technology, West Bengal, (Formerly known as WBUT), BF-142, Sector-I, Saltlake, Kolkata, 700 064, India
| | - Jaydeb Chakrabarti
- Department of Chemical, Biological and Macro-Molecular Sciences, S.N. Bose National Centre for Basic Sciences, Sector III, Block JD, Salt Lake, Kolkata, 700106, India.,The Thematic Unit of Excellence on Computational Materials Science, S. N. Bose National Centre for Basic Sciences, Sector-III, Block JD, Salt Lake, Kolkata, 700106, India
| |
Collapse
|
21
|
Abstract
In addition to continuous rapid progress in RNA structure determination, probing, and biophysical studies, the past decade has seen remarkable advances in the development of a new generation of RNA folding theories and models. In this article, we review RNA structure prediction models and models for ion-RNA and ligand-RNA interactions. These new models are becoming increasingly important for a mechanistic understanding of RNA function and quantitative design of RNA nanotechnology. We focus on new methods for physics-based, knowledge-based, and experimental data-directed modeling for RNA structures and explore the new theories for the predictions of metal ion and ligand binding sites and metal ion-dependent RNA stabilities. The integration of these new methods with theories about the cellular environment effects in RNA folding, such as molecular crowding and cotranscriptional kinetic effects, may ultimately lead to an all-encompassing RNA folding model.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Physics, Department of Biochemistry, and MU Informatics Institute, University of Missouri, Columbia, Missouri 65211;
| | - Dong Zhang
- Department of Physics, Department of Biochemistry, and MU Informatics Institute, University of Missouri, Columbia, Missouri 65211;
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and MU Informatics Institute, University of Missouri, Columbia, Missouri 65211;
| |
Collapse
|
22
|
Šponer J, Krepl M, Banáš P, Kührová P, Zgarbová M, Jurečka P, Havrila M, Otyepka M. How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes? WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27863061 DOI: 10.1002/wrna.1405] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 01/01/2023]
Abstract
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| |
Collapse
|
23
|
Ochieng PO, White NA, Feig M, Hoogstraten CG. Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation. J Phys Chem B 2016; 120:10885-10898. [PMID: 27701852 DOI: 10.1021/acs.jpcb.6b05606] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.
Collapse
Affiliation(s)
- Patrick O Ochieng
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Neil A White
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Charles G Hoogstraten
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|