Salt-Dependent Folding Energy Landscape of RNA Three-Way Junction

RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg²⁺ ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na⁺ and Mg²⁺ ion concentrations. For Mg²⁺ solutions, the TBI model, which accounts for the Mg²⁺ ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg²⁺ in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.     

Predicting Ion Binding Properties for RNA Tertiary Structures

Recent experiments pointed to the potential importance of ion correlation for multivalent ions such as Mg²⁺ ions in RNA folding. In this study, we develop an all-atom model to predict the ion electrostatics in RNA folding. The model can treat ion correlation effects explicitly by considering an ensemble of discrete ion distributions. In contrast to the previous coarse-grained models that can treat ion correlation, this new model is based on all-atom nucleic acid structures. Thus, unlike the previous coarse-grained models, this new model allows us to treat complex tertiary structures such as HIV-1 DIS type RNA kissing complexes. Theory-experiment comparisons for a variety of tertiary structures indicate that the model gives improved predictions over the Poisson-Boltzmann theory, which underestimates the Mg²⁺ binding in the competition with Na⁺. Further systematic theory-experiment comparisons for a series of tertiary structures lead to a set of analytical formulas for Mg²⁺/Na⁺ ion-binding to various RNA and DNA structures over a wide range of Mg²⁺ and Na⁺ concentrations.     

TBI Server: A Web Server for Predicting Ion Effects in RNA Folding

Background

Metal ions play a critical role in the stabilization of RNA structures. Therefore, accurate prediction of the ion effects in RNA folding can have a far-reaching impact on our understanding of RNA structure and function. Multivalent ions, especially Mg²⁺, are essential for RNA tertiary structure formation. These ions can possibly become strongly correlated in the close vicinity of RNA surface. Most of the currently available software packages, which have widespread success in predicting ion effects in biomolecular systems, however, do not explicitly account for the ion correlation effect. Therefore, it is important to develop a software package/web server for the prediction of ion electrostatics in RNA folding by including ion correlation effects.

Results

The TBI web server http://rna.physics.missouri.edu/tbi_index.html provides predictions for the total electrostatic free energy, the different free energy components, and the mean number and the most probable distributions of the bound ions. A novel feature of the TBI server is its ability to account for ion correlation and ion distribution fluctuation effects.

Conclusions

By accounting for the ion correlation and fluctuation effects, the TBI server is a unique online tool for computing ion-mediated electrostatic properties for given RNA structures. The results can provide important data for in-depth analysis for ion effects in RNA folding including the ion-dependence of folding stability, ion uptake in the folding process, and the interplay between the different energetic components.     

Hexahydrated Mg2+ Binding and Outer-Shell Dehydration on RNA Surface

The interaction between metal ions, especially Mg²⁺ ions, and RNA plays a critical role in RNA folding. Upon binding to RNA, a metal ion that is fully hydrated in bulk solvent can become dehydrated. Here we use molecular dynamics simulation to investigate the dehydration of bound hexahydrated Mg²⁺ ions. We find that a hydrated Mg²⁺ ion in the RNA groove region can involve significant dehydration in the outer hydration shell. The first or innermost hydration shell of the Mg²⁺ ion, however, is retained during the simulation because of the strong ion-water electrostatic attraction. As a result, water-mediated hydrogen bonding remains an important form for Mg²⁺-RNA interaction. Analysis for ions at different binding sites shows that the most pronounced water deficiency relative to the fully hydrated state occurs at a radial distance of around 11 Å from the center of the ion. Based on the independent 200 ns molecular dynamics simulations for three different RNA structures (Protein Data Bank: 1TRA, 2TPK, and 437D), we find that Mg²⁺ ions overwhelmingly dominate over monovalent ions such as Na⁺ and K⁺ in ion-RNA binding. Furthermore, application of the free energy perturbation method leads to a quantitative relationship between the Mg²⁺ dehydration free energy and the local structural environment. We find that $\text{ΔΔ}{G}_{\text{hyd}},$ the change of the Mg²⁺ hydration free energy upon binding to RNA, varies linearly with the inverse distance between the Mg²⁺ ion and the nearby nonbridging oxygen atoms of the phosphate groups, and $\text{ΔΔ}{G}_{\text{hyd}}$ can reach ?2.0 kcal/mol and ?3.0 kcal/mol for an Mg²⁺ ion bound to the surface and to the groove interior, respectively. In addition, the computation results in an analytical formula for the hydration ratio as a function of the average inverse Mg²⁺-O distance. The results here might be useful for further quantitative investigations of ion-RNA interactions in RNA folding.     

Mg2+ Impacts the Twister Ribozyme through Push-Pull Stabilization of Nonsequential Phosphate Pairs

RNA molecules perform a variety of biological functions for which the correct three-dimensional structure is essential, including as ribozymes where they catalyze chemical reactions. Metal ions, especially Mg²⁺, neutralize these negatively charged nucleic acids and specifically stabilize RNA tertiary structures as well as impact the folding landscape of RNAs as they assume their tertiary structures. Specific binding sites of Mg²⁺ in folded conformations of RNA have been studied extensively; however, the full range of interactions of the ion with compact intermediates and unfolded states of RNA is challenging to investigate, and the atomic details of the mechanism by which the ion facilitates tertiary structure formation is not fully known. Here, umbrella sampling combined with oscillating chemical potential Grand Canonical Monte Carlo/molecular dynamics simulations are used to capture the energetics and atomic-level details of Mg²⁺-RNA interactions that occur along an unfolding pathway of the Twister ribozyme. The free energy profiles reveal stabilization of partially unfolded states by Mg²⁺, as observed in unfolding experiments, with this stabilization being due to increased sampling of simultaneous interactions of Mg²⁺ with two or more nonsequential phosphate groups. Notably, these results indicate a push-pull mechanism in which the Mg²⁺-RNA interactions actually lead to destabilization of specific nonsequential phosphate-phosphate interactions (i.e., pushed apart), whereas other interactions are stabilized (i.e., pulled together), a balance that stabilizes unfolded states and facilitates the folding of Twister, including the formation of hydrogen bonds associated with the tertiary structure. This study establishes a better understanding of how Mg²⁺-ion interactions contribute to RNA structural properties and stability.     

Monovalent and divalent promoted GAAA tetraloop-receptor tertiary interactions from freely diffusing single-molecule studies

Proper assembly of RNA into catalytically active three-dimensional structures requires multiple tertiary binding interactions, individual characterization of which is crucial to a detailed understanding of global RNA folding. This work focuses on single-molecule fluorescence studies of freely diffusing RNA constructs that isolate the GAAA tetraloop-receptor tertiary interaction. Freely diffusing conformational dynamics are explored as a function of Mg²⁺ and Na⁺ concentration, both of which promote facile docking, but with 500-fold different affinities. Systematic shifts in mean fluorescence resonance energy transfer efficiency values and line widths with increasing [Na⁺] are observed for the undocked species and can be interpreted with a Debye model in terms of electrostatic relaxation and increased flexibility in the RNA. Furthermore, we identify a 34 ± 2% fraction of freely diffusing RNA constructs remaining undocked even at saturating [Mg²⁺] levels, which agrees quantitatively with the 32 ± 1% fraction previously reported for immobilized constructs. This verifies that the kinetic heterogeneity observed in the docking rates is not the result of surface tethering. Finally, the K_D value and Hill coefficient for [Mg²⁺]-dependent docking decrease significantly for [Na⁺] = 25 mM vs. 125 mM, indicating Mg²⁺ and Na⁺ synergy in the RNA folding process.     

Ion Condensation onto Ribozyme Is Site Specific and Fold Dependent

The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state, although Mg²⁺ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and nonspecifically bound to the RNA, as predicted by the counterion condensation theory. Here, we show using simulations of Azoarcus ribozyme, based on an accurate coarse-grained three-site interaction model with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg²⁺ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods, in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic nonmonotonic shape fluctuations in the ribozyme as it folds with increasing Mg²⁺ or Ca²⁺ concentration.     

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl₂), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg²⁺ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl₂) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals $\text{Δ}G?$ 10 kcal/mol. Besides stabilizing the 3WJ, Mg²⁺ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.     

Salt dependence of nucleic acid hairpin stability

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na⁺ and Mg²⁺ in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na⁺ and Mg²⁺ concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.     

The Osmolyte TMAO Stabilizes Native RNA Tertiary Structures in the Absence of Mg: Evidence for a Large Barrier to Folding from Phosphate Dehydration

The stabilization of RNA tertiary structures by ions is well known, but the neutral osmolyte trimethylamine oxide (TMAO) can also effectively stabilize RNA tertiary structure. To begin to understand the physical basis for the effects of TMAO on RNA, we have quantitated the TMAO-induced stabilization of five RNAs with known structures. So-called m values, the increment in unfolding free energy per molal of osmolyte at constant KCl activity, are ∼ 0 for a hairpin secondary structure and between 0.70 and 1.85 kcal mol^− 1m^− 1 for four RNA tertiary structures (30-86 nt). Further analysis of two RNAs by small-angle X-ray scattering and hydroxyl radical probing shows that TMAO reduces the radius of gyration of the unfolded ensemble to the same endpoint as seen in titration with Mg²⁺ and that the structures stabilized by TMAO and Mg²⁺ are indistinguishable. Remarkably, TMAO induces the native conformation of a Mg²⁺ ion chelation site formed in part by a buried phosphate, even though Mg²⁺ is absent. TMAO interacts weakly, if at all, with KCl, ruling out the possibility that TMAO stabilizes RNA indirectly by increasing salt activity. TMAO is, however, strongly excluded from the vicinity of dimethylphosphate (unfavorable interaction free energy, + 211 cal mol^− 1m^− 1 for the potassium salt), an ion that mimics the RNA backbone phosphate. We suggest that formation of RNA tertiary structure is accompanied by substantial phosphate dehydration (loss of 66-173 water molecules in the RNA structures studied) and that TMAO works principally by reducing the energetic penalty associated with this dehydration. The strong parallels we find between the effects of TMAO and Mg²⁺ suggest that RNA sequence is more important than specific ion interactions in specifying the native structure.     

Reduced Model Captures Mg-RNA Interaction Free Energy of Riboswitches

The stability of RNA tertiary structures depends heavily on Mg²⁺. The Mg²⁺-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ₂₊, the number of excess Mg²⁺ associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg²⁺ ions interact closely and strongly with the RNA, unlike monovalent ions such as K⁺, suggesting that an explicit treatment of Mg²⁺ is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg²⁺-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg²⁺ through the explicit representation of Mg²⁺ ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K⁺ and models its competition with condensed Mg²⁺. The model contains one fitted parameter, the number of condensed K⁺ ions in the absence of Mg²⁺. Values of Γ₂₊ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ₂₊ to the Mg²⁺-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg²⁺ and RNA dynamics.     

Reduced Model Captures Mg2+-RNA Interaction Free Energy of Riboswitches

The stability of RNA tertiary structures depends heavily on Mg²⁺. The Mg²⁺-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ₂₊, the number of excess Mg²⁺ associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg²⁺ ions interact closely and strongly with the RNA, unlike monovalent ions such as K⁺, suggesting that an explicit treatment of Mg²⁺ is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg²⁺-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg²⁺ through the explicit representation of Mg²⁺ ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K⁺ and models its competition with condensed Mg²⁺. The model contains one fitted parameter, the number of condensed K⁺ ions in the absence of Mg²⁺. Values of Γ₂₊ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ₂₊ to the Mg²⁺-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg²⁺ and RNA dynamics.     

Electrostatic free energy landscapes for nucleic acid helix assembly

Metal ions are crucial for nucleic acid folding. From the free energy landscapes, we investigate the detailed mechanism for ion-induced collapse for a paradigm system: loop-tethered short DNA helices. We find that Na⁺ and Mg²⁺ play distinctive roles in helix–helix assembly. High [Na⁺] (>0.3 M) causes a reduced helix–helix electrostatic repulsion and a subsequent disordered packing of helices. In contrast, Mg²⁺ of concentration >1 mM is predicted to induce helix–helix attraction and results in a more compact and ordered helix–helix packing. Mg²⁺ is much more efficient in causing nucleic acid compaction. In addition, the free energy landscape shows that the tethering loops between the helices also play a significant role. A flexible loop, such as a neutral loop or a polynucleotide loop in high salt concentration, enhances the close approach of the helices in order to gain the loop entropy. On the other hand, a rigid loop, such as a polynucleotide loop in low salt concentration, tends to de-compact the helices. Therefore, a polynucleotide loop significantly enhances the sharpness of the ion-induced compaction transition. Moreover, we find that a larger number of helices in the system or a smaller radius of the divalent ions can cause a more abrupt compaction transition and a more compact state at high ion concentration, and the ion size effect becomes more pronounced as the number of helices is increased.     

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg²⁺ ion concentrations are low, K⁺ concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo–like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg²⁺ (0.5–2 mM) and K⁺ (140 mM) if the solution is supplemented with physiological amounts (∼20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

88 Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamic simulation and Poisson–Boltzmann approach

Ion interactions with nucleic acids (both DNA and RNA) are an important and evolving field of investigation. Positively charged cations may interact with highly negatively charged nucleic acids via simple electrostatic interactions to help screen the electrostatic repulsion along the nucleic acids and assist their folding and/or compaction. Cations may also bind at specific sites and become integral parts of the structures, possibly playing important enzymatic roles. Two popular methods for computationally exploring a nucleic acid’s ion atmosphere are atomistic molecular dynamics (MD) simulations and the Poisson–Boltzmann (PB) equation. In general, monovalent ion results obtained from MD simulations and the PB equation agree well with experiment. However, Bai et al. (2007) observed discrepancies between experiment and the PB equation while examining the competitive binding of monovalent and divalent ions, with more significant discrepancies for divalent ions. The goal of this project was to thoroughly investigate monovalent (Na⁺) and divalent (Mg²⁺) ion distributions formed around a DNA duplex with MD simulations and the PB equation. We simulated three different cation concentrations, and matched the equilibrated bulk ion concentration for our theoretical calculations with the PB equation. Based on previous work, our Mg²⁺ ions were fully solvated, the expected state of Mg²⁺ ions when interacting with a duplex, when the production simulations began and remained throughout the simulations (Kirmizialtin, 2010; Robbins, 2012). Na⁺ ion distributions and number of Na⁺ ions within 10?Å of the DNA obtained from our two methods agreed well. However, results differed for Mg²⁺ ions, with a lower number of ions within the cut-off distance obtained from the PB equation when compared to MD simulations. The Mg²⁺ ion distributions around the DNA obtained via the two methods also differed. Based on our results, we conclude that the PB equation will systematically underestimate Mg²⁺ ions bound to DNA, and much of this deviation is attributed to dielectric saturation associated with high valency ions.     

A large collapsed-state RNA can exhibit simple exponential single-molecule dynamics

The process of large RNA folding is believed to proceed from many collapsed structures to a unique functional structure requiring precise organization of nucleotides. The diversity of possible structures and stabilities of large RNAs could result in non-exponential folding kinetics (e.g. stretched exponential) under conditions where the molecules have not achieved their native state. We describe a single-molecule fluorescence resonance energy transfer (FRET) study of the collapsed-state region of the free energy landscape of the catalytic domain of RNase P RNA from Bacillus stearothermophilus (C_thermo). Ensemble measurements have shown that this 260 residue RNA folds cooperatively to its native state at ≥1 mM Mg²⁺, but little is known about the conformational dynamics at lower ionic strength. Our measurements of equilibrium conformational fluctuations reveal simple exponential kinetics that reflect a small number of discrete states instead of the expected inhomogeneous dynamics. The distribution of discrete dwell times, collected from an “ensemble” of 300 single molecules at each of a series of Mg²⁺ concentrations, fit well to a double exponential, which indicates that the RNA conformational changes can be described as a four-state system. This finding is somewhat unexpected under [Mg²⁺] conditions in which this RNA does not achieve its native state. Observation of discrete well-defined conformations in this large RNA that are stable on the seconds timescale at low [Mg²⁺] (<0.1 mM) suggests that even at low ionic strength, with a tremendous number of possible (weak) interactions, a few critical interactions may produce deep energy wells that allow for rapid averaging of motions within each well, and yield kinetics that are relatively simple.     

Molecular beacons as probes of RNA unfolding under native conditions

Hybridization of fluorescent molecular beacons provides real-time detection of RNA secondary structure with high specificity. We used molecular beacons to measure folding and unfolding rates of the Tetrahymena group I ribozyme under native conditions. A molecular beacon targeted against 15 nt in the 5′ strand of the P3 helix specifically hybridized with misfolded forms of the ribozyme, without invading the native tertiary structure. The beacon associated with the misfolded ribozyme 300 times more slowly than with an unstructured oligonucleotide containing the same target sequence, suggesting that the misfolded ribozyme core remains structured in the absence of Mg²⁺. The rate of beacon hybridization under native conditions revealed a linear relationship between the free energy of unfolding and Mg²⁺ concentration. A small fraction of the RNA population unfolded very rapidly, suggesting parallel unfolding in one step or through misfolded intermediates.     

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2′-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg²⁺ ion facilitates deprotonation of the O2′ nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5′ leaving group as noted in a previous study. This model assumes that the active site Mg²⁺ ion forms an inner-sphere coordination with the O2′ nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg²⁺ ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg²⁺ are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg²⁺ is simply replaced with Na⁺ is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg²⁺-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.