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.

A Single‐Stranded Junction Modulates Nanosecond Motional Ordering of the Substrate Recognition Duplex of a Group I Ribozyme

Single Molecule Structural Dynamics of Ribosome Before and During Translocation

Motions of the Substrate Recognition Duplex in a Group I Intron Assessed by Site-Directed Spin Labeling

A mechanism for S-adenosyl methionine assisted formation of a riboswitch conformation: a small molecule with a strong arm

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

Conformational Distributions at the N-peptide/boxB RNA Interface Studied Using Site-Directed Spin Labeling

A transient conformation facilitates ligand binding to the adenine riboswitch

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

Synchronous RNA conformational changes trigger ordered phase transitions in crystals

Cofactor binding triggers rapid conformational remodelling of the active site of a methyltransferase ribozyme

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time

Nanosecond chain dynamics of single-stranded nucleic acids

RNA under Tension: Folding Landscapes, Kinetic Partitioning Mechanism, and Molecular Tensegrity

Conformational capture of the SAM-II riboswitch

Sequence-dependent folding landscapes of adenine riboswitch aptamers

Interrogating RNA-Small Molecule Interactions with Structure Probing and Artificial Intelligence-Augmented Molecular Simulations

Conformational Landscapes of a Class I Ribonucleotide Reductase Complex during Turnover Reveal Intrinsic Dynamics and Asymmetry

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

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

A minimal RNA substrate with dual fluorescent probes enables rapid kinetics and provides insight into bacterial RNase P active site interactions

RNA Structural Modules Control the Rate and Pathway of RNA Folding and Assembly.