Abstract:Self-splicing ribozymes are small RNA enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoesterbond formation is thermodynamically very unfavorable, and many experimental studies have focused on the understanding of the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so that molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme and highlight that identifying a relevant initial conformation for reactivity studies can be highly challenging due to limitations in both available X-ray and the force field accuracy, together with the necessity of extensive sampling. This is frequently overlooked in mixed quantum/classical studies that predominantly concentrate on the chemical reaction itself. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state, prior to self-cleavage. With all investigated force fields, we find that what is usually considered as a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis, is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues to the reaction. Our study therefore demonstrates that identifying the most pertinent reactant state conformation, holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center

Simulations predict preferred Mg coordination in a nonenzymatic primer extension reaction center

Catalytic Metal Ion-Substrate Coordination during Nonenzymatic RNA Primer Extension

Unusual Base-Pairing Interactions in Monomer-Template Complexes.

Mechanism of the Nucleotidyl-Transfer Reaction in DNA Polymerase Revealed by Time-Resolved Protein Crystallography

Computational evaluation of the binding of activated ribonucleotides in nonenzymatic RNA template copying

Efficient synthesis of the 2-aminoimidazolium-bridged diguanosyl intermediate for nonenzymatic primer extension

Structural Rationale for the Enhanced Catalysis of Nonenzymatic RNA Primer Extension by a Downstream Oligonucleotide

A Potential Role for Aminoacylation in Primordial RNA Copying Chemistry

Constraints on the Emergence of RNA through Non-Templated Primer Extension with Mixtures of Potentially Prebiotic Nucleotides

Diaminopurine in Nonenzymatic RNA Template Copying

Primer terminal ribonucleotide alters the active site dynamics of DNA polymerase η and reduces DNA synthesis fidelity

A mechanistic model of primer synthesis from catalytic structures of DNA polymerase α–primase

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA

Mechanism of Three-Component Collision to Produce Ultrastable Prna Three-Way Junction of Phi29 DNA-packaging Motor by Kinetic Assessment.

Kinetic Proofreading can Enhance Specificity in a Non-enzymatic DNA Strand Displacement Network

Overcoming nucleotide bias in the nonenzymatic copying of RNA templates

Membrane Catalyzed Formation of Nucleotide Clusters and Their Role in the Origins of Life: Insights from Molecular Simulations and Lattice Modeling

Preparation of 2-Aminoimidazole-Activated Substrates for the Study of Nonenzymatic Genome Replication