Differences in ion-RNA binding modes due to charge density variations explain the stability of RNA in monovalent salts

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods

RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.

Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability

The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

Monovalent metal ion binding promotes the first transesterification reaction in the spliceosome

Cleavage and formation of phosphodiester bonds in nucleic acids is accomplished by large cellular machineries composed of both protein and RNA. Long thought to rely on a two-metal-ion mechanism for catalysis, structure comparisons revealed many contain highly spatially conserved second-shell monovalent cations, whose precise function remains elusive. A recent high-resolution structure of the spliceosome, essential for pre-mRNA splicing in eukaryotes, revealed a potassium ion in the active site. Here, we employ biased quantum mechanics/ molecular mechanics molecular dynamics to elucidate the function of this monovalent ion in splicing. We discover that the K+ ion regulates the kinetics and thermodynamics of the first splicing step by rigidifying the active site and stabilizing the substrate in the pre- and post-catalytic state via formation of key hydrogen bonds. Our work supports a direct role for the K+ ion during catalysis and provides a mechanistic hypothesis likely shared by other nucleic acid processing enzymes.

Why Na+ has higher propensity than K+ to condense DNA in a crowded environment

Experimentally, in the presence of the crowding agent polyethylene glycol (PEG), sodium ions compact double-stranded DNA more readily than potassium ions. Here, we have used molecular dynamics simulations and the "ion binding shells model" of DNA condensation to provide an explanation for the observed variations in condensation of short DNA duplexes in solutions containing different monovalent cations and PEG; several predictions are made. According to the model we use, externally bound ions contribute the most to the ion-induced aggregation of DNA duplexes. The simulations reveal that for two adjacent DNA duplexes, the number of externally bound Na+ ions is larger than the number of K+ ions over a wide range of chloride concentrations in the presence of PEG, providing a qualitative explanation for the higher propensity of sodium ions to compact DNA under crowded conditions. The qualitative picture is confirmed by an estimate of the corresponding free energy of DNA aggregation that is at least 0.2kBT per base pair more favorable in solution with NaCl than with KCl at the same ion concentration. The estimated attraction free energy of DNA duplexes in the presence of Na+ depends noticeably on the DNA sequence; we predict that AT-rich DNA duplexes are more readily condensed than GC-rich ones in the presence of Na+. Counter-intuitively, the addition of a small amount of a crowding agent with high affinity for the specific condensing ion may lead to the weakening of the ion-mediated DNA-DNA attraction, shifting the equilibrium away from the DNA condensed phase.

Sequence-Dependent Orientational Coupling and Electrostatic Attraction in Cation-Mediated DNA-DNA Interactions

Condensation of DNA is vital for its biological functions and controlled nucleic acid assemblies. However, the mechanisms of DNA condensation are not fully understood due to the inability of experiments to access cation distributions and the complex interplay of energetic and entropic forces during assembly. By constructing free energy surfaces using exhaustive sampling and detailed analysis of cation distributions, we elucidate the mechanism of DNA condensation in different salt conditions and with different DNA sequences. We found that DNA condensation is facilitated by the correlated dynamics of the localized cations at the grooves of DNA helices. These dynamics are strongly dependent on the salt conditions and DNA sequences. In the presence of magnesium ions, major groove binding facilitates attraction. In contrast, in the presence of polyvalent cations, minor groove binding serves to create charge patterns, leading to condensation. Our findings present a novel advancement in the field and have broad implications for understanding and controlling nucleic acid complexes in vivo and in vitro.

Decoding the dual recognition mechanism of the glucocorticoid receptor for DNA and RNA: sequence versus shape

Transcription factors (TFs) regulate eukaryotic transcription through selective DNA-binding, can also specifically interact with RNA, which may present another layer of transcriptional control. The mechanisms of the TFs-DNA recognition are often well-characterised, while the details of TFs-RNA complexation are less understood. Here we investigate the dual recognition mechanism of the glucocorticoid receptor (GR), which interacts with similar affinities with consensus DNA and diverse RNA hairpin motifs but discriminates against uniform dsRNA. Using atomic molecular dynamics simulations, we demonstrate that the GR binding to nucleic acids requires a wide and shallow groove pocket. The protein effectively moulds its binding site within DNA major groove, which enables base-specific interactions. Contrary, the GR binding has little effect on the grooves geometry of RNA systems, most notably in uniform dsRNA. Instead, a hairpin motif in RNA yields a wide and shallow major groove pocket, allowing the protein to anchor itself through nonspecific electrostatic contacts with RNA backbone. Addition of a bulge increases RNA hairpin flexibility, which leads to a greater number of GR-RNA contacts and, thus, higher affinity. Thus, the combination of structural motifs defines the GR-RNA selective binding: a recognition mechanism, which may be shared by other zinc finger TFs.

Markedly Different Effects of Monovalent Cations on the Efficiency of Gene Expression

The effect of monovalent cations on a cell-free transcription-translation (TX-TL) system is examined using a luciferase assay. It is found that the potency for all ions analyzed here is in the order Rb⁺  > K⁺  > Cs⁺  > Na⁺  ≈ Li⁺  > (CH₃ )₄ N⁺ , where Rb⁺ is most efficient at promoting TX-TL and the ions of Li⁺ , Na⁺ , and (CH₃ )₄ N⁺ exhibit an inhibitory effect. Similar promotion/inhibition effects are observed for cell-free TL alone with an mRNA template.