Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation

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

Self-splicing ribozymes are small ribonucleic acid (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 phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. 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 us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered 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 in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study 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.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.