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.

Entropically Driven Enhancement of Cleavage Activity of a DNA-Armed Hammerhead Ribozyme: Mechanism of Action of Hammerhead Ribozymes

Mechanism of Action of A Hammerhead Ribozyme, A Metalloenzyme

Explanation by the Double-Metal-Ion Mechanism of Catalysis for the Differential Metal Ion Effects on the Cleavage Rates of 5 '-Oxy and 5 '-Thio Substrates by A Hammerhead Ribozyme

Comparison of Metal-Ion-dependent Cleavages of RNA by a DNA Enzyme and a Hammerhead Ribozyme.

Ribozyme Mechanism Revisited: Evidence Against Direct Coordination of A Mg2+ Ion with the Pro-R Oxygen of the Scissile Phosphate in the Transition State of A Hammerhead Ribozyme-Catalyzed Reaction

Existence Of Efficient Divalent Metal Ion-Catalyzed And Inefficient Divalent Metal Ion-Independent Channels In Reactions Catalyzed By A Hammerhead Ribozyme

Explanation by a putative triester-like mechanism for the thio effects and Mn2+ rescues in reactions catalyzed by a hammerhead ribozyme

Evidence For The Rate-Limiting Departure Of The 5'-Oxygen In Nonenzymatic And Hammerhead Ribozyme-Catalyzed Reactions

The structure and catalytic mechanism of a pseudoknot-containing hammerhead ribozyme

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

The Hydrolysis of RNA: from Theoretical Calculations to the Hammerhead Ribozyme-Mediated Cleavage of RNA

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations.

Complete RNA inverse folding: computational design of functional hammerhead ribozymes

Crystal Structure of an RNA-cleaving DNAzyme

Structural insights into catalysis and dimerization enhanced exonuclease activity of RNase J.

Differential Assembly of Catalytic Interactions within the Conserved Active Sites of Two Ribozymes

Retraining and Optimizing DNA-Hydrolyzing Deoxyribozymes for Robust Single- and Multiple-Turnover Activities

Engineering a ribozyme cleavage-induced split fluorescent aptamer complementation assay

Exploring CRISPR-Cas9 HNH-Domain Catalyzed DNA Cleavage Using Accelerated Quantum Mechanical Molecular Mechanical Free Energy Simulation

Exploring the catalytic mechanism of the 10-23 DNAzyme: insights from pH-rate profiles

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