Molecular mechanism of the site-specific self-cleavage of the RNA phosphodiester backbone by a twister ribozyme

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

Ribozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries along independent MD simulations with different protonation states of the active site basic species reveals that only the canonical system, with no additional protonation changes, renders reactive conformations. A change in the coordination sphere of the Mg2+ ion has been observed during the simulations, which allows proposing a mechanism to explain the unique mode of action of the pistol ribozyme by comparison with other ribozymes. The present results are at the center of the debate originated from recent experimental and theoretical studies on pistol ribozyme.

Cleaning Up Mechanistic Debris Generated by Twister Ribozymes Using Computational RNA Enzymology

The catalytic properties of RNA have been a subject of fascination and intense research since their discovery over 30 years ago. Very recently, several classes of nucleolytic ribozymes have emerged and been characterized structurally. Among these, the twister ribozyme has been center-stage, and a topic of debate about its architecture and mechanism owing to conflicting interpretations of different crystal structures, and in some cases conflicting interpretations of the same functional data. In the present work, we attempt to clean up the mechanistic "debris" generated by twister ribozymes using a comprehensive computational RNA enzymology approach aimed to provide a unified interpretation of existing structural and functional data. Simulations in the crystalline environment and in solution provide insight into the origins of observed differences in crystal structures, and coalesce on a common active site architecture, and dynamical ensemble in solution. We use GPU-accelerated free energy methods with enhanced sampling to ascertain microscopic nucleobase pK a values of the implicated general acid and base, from which predicted activity-pH profiles can be compared directly with experiments. Next, ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with full dynamic solvation under periodic boundary conditions are used to determine mechanistic pathways through multi-dimensional free energy landscapes for the reaction. We then characterize the rate-controlling transition state, and make predictions about kinetic isotope effects and linear free energy relations. Computational mutagenesis is performed to explain the origin of rate effects caused by chemical modifications and make experimentally testable predictions. Finally, we provide evidence that helps to resolve conflicting issues related to the role of metal ions in catalysis. Throughout each stage, we highlight how a conserved L-platform structural motif, to- gether with a key L-anchor residue, forms the characteristic active site scaffold enabling each of the catalytic strategies to come together not only for the twister ribozyme, but the majority of the known small nucleolytic ribozyme classes.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

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.

Phosphoryl Transfer Reaction in RNA: Is the Substrate-Assisted Catalysis a Possible Mechanism in Certain Solvents?

A proton shuttle mechanism for the phosphoryl transfer reaction in RNA, in which a proton is transferred from the nucleophile to the leaving group through a nonbridged oxygen atom of the phosphate, was explored using the MO6-2X density functional method and the solvent continuum model. This reaction is the initial step of the RNA hydrolysis. We used different solvents characterized by their dielectric constant, and, for each of them, we studied the nuclear and electronic relaxation, produced by the solvent reaction field, for the stationary points. Given that RNA has a poor leaving group, the bond breaking corresponds to the rate-determining step. If the O atom is substituted by a S atom, the leaving group is now good, and the rate-determining step is now the nucleophilic attack concerted with the proton transfer. The most relevant result we found is that none of the solvents we studied has a free energy of activation that is smaller than the one in water. This suggests that the enzyme catalysis following this mechanism must be due to the permanent electric field that is created by a preorganized charge distribution but not to the solvent reaction field.