Molecular Dynamics Analysis of Mg2+-dependent Cleavage of a Pistol Ribozyme Reveals a Fail-safe Secondary Ion for Catalysis

Pistol ribozymes comprise a class of small, self-cleaving RNAs discovered via comparative genomic analysis. Prior work in the field has probed the kinetics of the cleavage reaction, as well as the influence of various metal ion cofactors that accelerate the process. In the current study we perform unbiased and unconstrained molecular dynamics simulations from two current high-resolution pistol crystal structures, and we analyzed trajectory data within the context of the currently accepted ribozyme mechanistic framework. Root-mean-squared deviations (RMSDs), radial distribution functions (RDFs), and distributions of nucleophilic angle-of-attack reveal insights into the potential roles of three magnesium ions with respect to catalysis and overall conformational stability of the molecule. A series of simulation trajectories containing in-silico mutations reveal the relatively flexible and partially interchangeable roles of two particular magnesium ions within solvated hydrogen-bonding distances from the catalytic center.

[1]  D. York,et al.  Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme , 2019, RNA.

[2]  Andreas Prlic,et al.  NGL viewer: web‐based molecular graphics for large complexes , 2018, Bioinform..

[3]  Alexander S. Rose,et al.  NGLview–interactive molecular graphics for Jupyter notebooks , 2018, Bioinform..

[4]  R. Micura,et al.  Atom-Specific Mutagenesis Reveals Structural and Catalytic Roles for an Active-Site Adenosine and Hydrated Mg2+ in Pistol Ribozymes. , 2017, Angewandte Chemie.

[5]  Alexander S. Rose,et al.  MDsrv: viewing and sharing molecular dynamics simulations on the web , 2017, Nature Methods.

[6]  Kiyoung Lee,et al.  Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes , 2017, Molecules.

[7]  R. Breaker.  Mechanistic Debris Generated by Twister Ribozymes. , 2017, ACS chemical biology.

[8]  T. Steitz,et al.  Crystal structure of Pistol, a class of self-cleaving ribozyme , 2016, Proceedings of the National Academy of Sciences.

[9]  Pu Gao,et al.  Pistol ribozyme adopts a pseudoknot fold facilitating site-specific in-line cleavage , 2016, Nature chemical biology.

[10]  P. Bevilacqua,et al.  Molecular Dynamics Study of Twister Ribozyme: Role of Mg2+ Ions and the Hydrogen-Bonding Network in the Active Site , 2016, Biochemistry.

[11]  Wei Yang,et al.  Capture of a third Mg2+ is essential for catalyzing DNA synthesis , 2016, Science.

[12]  R. Breaker,et al.  Biochemical analysis of pistol self-cleaving ribozymes , 2015, RNA.

[13]  Travis E. Oliphant,et al.  Guide to NumPy , 2015 .

[14]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[15]  R. Breaker,et al.  New classes of self-cleaving ribozymes revealed by comparative genomics analysis , 2015, Nature chemical biology.

[16]  Tijana Radivojević,et al.  Constant pressure hybrid Monte Carlo simulations in GROMACS , 2014, Journal of Molecular Modeling.

[17]  Yijin Liu,et al.  Crystal structure and mechanistic investigation of the twister ribozyme. , 2014, Nature chemical biology.

[18]  R. Breaker,et al.  A widespread self-cleaving ribozyme class is revealed by bioinformatics , 2013, Nature chemical biology.

[19]  Berk Hess,et al.  A flexible algorithm for calculating pair interactions on SIMD architectures , 2013, Comput. Phys. Commun..

[20]  D. Lilley.  Mechanisms of RNA catalysis , 2011, Philosophical Transactions of the Royal Society B: Biological Sciences.

[21]  Oliver Beckstein,et al.  MDAnalysis: A toolkit for the analysis of molecular dynamics simulations , 2011, J. Comput. Chem..

[22]  Tanneguy Redarce,et al.  Automatic Lip-Contour Extraction and Mouth-Structure Segmentation in Images , 2011, Computing in Science & Engineering.

[23]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[24]  A. Ferré-D’Amaré,et al.  Small self-cleaving ribozymes. , 2010, Cold Spring Harbor perspectives in biology.

[25]  P. Emsley,et al.  Features and development of Coot , 2010, Acta crystallographica. Section D, Biological crystallography.

[26]  Randy J. Read,et al.  Acta Crystallographica Section D Biological , 2003 .

[27]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[28]  D. Lilley,et al.  The Evolution of Ribozyme Chemistry , 2009, Science.

[29]  Brian E. Granger,et al.  IPython: A System for Interactive Scientific Computing , 2007, Computing in Science & Engineering.

[30]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[31]  D. Theobald.  short communications Acta Crystallographica Section A Foundations of , 2005 .

[32]  R. Breaker,et al.  Control of gene expression by a natural metabolite-responsive ribozyme , 2004, Nature.

[33]  R. Breaker,et al.  A common speed limit for RNA-cleaving ribozymes and deoxyribozymes. , 2003, RNA.

[34]  Adam Roth,et al.  Ribozyme speed limits. , 2003, RNA.

[35]  M. Machius,et al.  Monovalent cation dependence and preference of GHKL ATPases and kinases 1 , 2003, FEBS letters.

[36]  J. Doudna,et al.  Ribozyme structures and mechanisms. , 2000, Annual review of biochemistry.

[37]  D. Herschlag,et al.  Identification of the hammerhead ribozyme metal ion binding site responsible for rescue of the deleterious effect of a cleavage site phosphorothioate. , 1999, Biochemistry.

[38]  R R Breaker,et al.  Relationship between internucleotide linkage geometry and the stability of RNA. , 1999, RNA.

[39]  W. Scott,et al.  The hammerhead, hairpin and VS ribozymes are catalytically proficient in monovalent cations alone. , 1998, Chemistry & biology.

[40]  J. Steitz,et al.  A general two-metal-ion mechanism for catalytic RNA. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[41]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[42]  J. Cowan,et al.  Metallobiochemistry of RNA. Co(NH3)6(3+) as a probe for Mg2+(aq) binding sites. , 1993, Journal of inorganic biochemistry.

[43]  R. Collins,et al.  A site-specific self-cleavage reaction performed by a novel RNA in neurospora mitochondria , 1990, Cell.

[44]  J. Taylor,et al.  Characterization of self-cleaving RNA sequences on the genome and antigenome of human hepatitis delta virus , 1988, Journal of virology.

[45]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[46]  J. M. Buzayan,et al.  Non-enzymatic cleavage and ligation of RNAs complementary to a plant virus satellite RNA , 1986, Nature.

[47]  R. Symons,et al.  Self-cleavage of plus and minus RNA transcripts of avocado sunblotch viroid. , 1986, Nucleic acids research.

[48]  J. M. Buzayan,et al.  Autolytic Processing of Dimeric Plant Virus Satellite RNA , 1986, Science.

[49]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[50]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[51]  J. VanLeeuwen,et al.  New method for the calculation of the pair correlation function. I , 1959 .

[52]  J. M. J. van Leeuwen,et al.  New method for the calculation of the pair correlation function. I , 1959 .