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.

Theoretical study on the mechanism of self-cleavage reaction of the glmS ribozyme

Role of the Active Site Guanine in the Glms Ribozyme Self-Cleavage Mechanism: Quantum Mechanical/molecular Mechanical Free Energy Simulations.

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

Assessing the Potential Effects of Active Site Mg2+ Ions in the Glms Ribozyme–Cofactor Complex

Theoretical studies on the thermodynamics and kinetics of the N-glycosidic bond cleavage in deoxythymidine glycol.

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions.

The Structural and Mechanistic Origins of Catalysis in Nucleolytic Ribozymes

QM/MM Simulation (B3LYP) of the RNase A Cleavage-Transesterification Reaction Supports a Triester A N + D N Associative Mechanism with an O2′ H Internal Proton Transfer

General Acid-Base Catalysis Mediated by Nucleobases in the Hairpin Ribozyme.

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

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

QM/MM investigation of the catalytic mechanism of processive endoglucanase Cel9G from Clostridium cellulovorans

Reaction Mechanism of the Binuclear Zinc Enzyme Glyoxalase Ii - A Theoretical Study

A QM/MM MD study of the pH-dependent ring-opening catalysis and lid motif flexibility in glucosamine 6-phosphate deaminase.

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

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates

QM/MM study on the catalytic mechanism of cellulose hydrolysis catalyzed by cellulase Cel5A from Acidothermus cellulolyticus.

Insights into the Catalytic Mechanism of Unsaturated Glucuronyl Hydrolase of Bacillus Sp. GL1.

The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme

Investigation of the rescue mechanism catalyzed by a nucleophile mutant of rice BGlu1

QM/MM Investigation on the Catalytic Mechanism of Bacteroides Thetaiotaomicron Α-Glucosidase BtGH97a.