Stabilization of RNA tertiary structure by monovalent cations

Advances on polymeric nanocarriers for sustainable agriculture: Enhancing dsRNA/siRNA delivery to combat agricultural pests

The application of exogenous RNA for gene-silencing strategies has gained significant traction in agriculture, offering a highly efficient and eco-friendly alternative to conventional plant protection methods. This success has been driven by advances in biotechnology, from the design of long double-stranded RNA (dsRNA) and small interfering RNA (siRNA) molecules to the development of nanocarrier systems that address the challenge of RNA delivery into plant cells. In particular, polymer-based nanocarriers have emerged as a promising solution for enhancing the stability and delivery efficiency of RNA molecules. This review provides a comprehensive overview of the current state of research on the use of polymeric nanocarriers in RNA interference (RNAi) systems for crop protection. It examines key technological developments that have enabled the effective delivery of dsRNA/siRNA to target organisms, with a focus on the unique advantages polymers offer as carriers. Recent studies highlight significant progress in the preparation, characterization, and application of polymeric nanocarriers for RNA encapsulation and delivery. The review also explores the environmental and health challenges posed by these technologies, emphasizing the need for sustainable approaches in their development. Specifically, the production of nanocarriers must adhere to the principles of green chemistry, prioritizing chemical modification routes that reduce harmful residues, such as toxic solvents. Finally, this paper discusses both the current challenges and future prospects of using polymer-based nanocarriers in sustainable agriculture, offering critical insights into their potential to transform crop protection through RNAi technologies.

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.

Osmolytes: Wonder molecules to combat protein misfolding against stress conditions

The proper functioning of any protein depends on its three dimensional conformation which is achieved by the accurate folding mechanism. Keeping away from the exposed stress conditions leads to cooperative unfolding and sometimes partial folding, forming the structures like protofibrils, fibrils, aggregates, oligomers, etc. leading to several neurodegenerative diseases like Parkinson's disease, Alzheimer's, Cystic fibrosis, Huntington, Marfan syndrome, and also cancers in some cases, too. Hydration of proteins is necessary, which may be achieved by the presence of organic solutes called osmolytes within the cell. Osmolytes belong to different classes in different organisms and play their role by preferential exclusion of osmolytes and preferential hydration of water molecules and achieves the osmotic balance in the cell otherwise it may cause problems like cellular infection, cell shrinkage leading to apoptosis and cell swelling which is also the major injury to the cell. Osmolyte interacts with protein, nucleic acids, intrinsically disordered proteins by non-covalent forces. Stabilizing osmolytes increases the Gibbs free energy of the unfolded protein and decreases that of folded protein and vice versa with denaturants (urea and guanidinium hydrochloride). The efficacy of each osmolyte with the protein is determined by the calculation of m value which reflects its efficiency with protein. Hence osmolytes can be therapeutically considered and used in drugs.

Thermodynamic Characterization of Nucleic Acid Nanoparticles Hybridization by UV Melting

The advances in nucleic acid nanotechnology have given rise to various elegantly designed structural complexes fabricated from DNA, RNA, chemically modified RNA strands, and their mixtures. The structural properties of NA nanoparticles (NANP) generally dictate and significantly impact biological function; and thus, it is critical to extract information regarding relative stabilities of the different structural forms. The adequate stability assessment requires knowledge of thermodynamic parameters that can be empirically derived using conventional UV-melting technique. The focus of this chapter is to describe methodology to evaluate thermodynamic data of NANPs complexation based on DNA 12 base-pair (bp) duplex formation as an example.

The Peptidyl Transferase Center: a Window to the Past

In his 2001 article, "Translation: in retrospect and prospect," the late Carl Woese made a prescient observation that there was a need for the then-current view of translation to be "reformulated to become an all-embracing perspective about which 21st century Biology can develop" (RNA 7:1055-1067, 2001, https://doi.org/10.1017/s1355838201010615). The quest to decipher the origins of life and the road to the genetic code are both inextricably linked with the history of the ribosome. After over 60 years of research, significant progress in our understanding of how ribosomes work has been made. Particularly attractive is a model in which the ribosome may facilitate an ∼180° rotation of the CCA end of the tRNA from the A-site to the P-site while the acceptor stem of the tRNA would then undergo a translation from the A-site to the P-site. However, the central question of how the ribosome originated remains unresolved. Along the path from a primitive RNA world or an RNA-peptide world to a proto-ribosome world, the advent of the peptidyl transferase activity would have been a seminal event. This functionality is now housed within a local region of the large-subunit (LSU) rRNA, namely, the peptidyl transferase center (PTC). The PTC is responsible for peptide bond formation during protein synthesis and is usually considered to be the oldest part of the modern ribosome. What is frequently overlooked is that by examining the origins of the PTC itself, one is likely going back even further in time. In this regard, it has been proposed that the modern PTC originated from the association of two smaller RNAs that were once independent and now comprise a pseudosymmetric region in the modern PTC. Could such an association have survived? Recent studies have shown that the extant PTC is largely depleted of ribosomal protein interactions. It is other elements like metallic ion coordination and nonstandard base/base interactions that would have had to stabilize the association of RNAs. Here, we present a detailed review of the literature focused on the nature of the extant PTC and its proposed ancestor, the proto-ribosome.

An uncommon [K+(Mg2+)2] metal ion triad imparts stability and selectivity to the Guanidine-I riboswitch

The widespread ykkC-I riboswitch class exemplifies divergent riboswitch evolution. To analyze how natural selection has diversified its versatile RNA fold, we determined the X-ray crystal structure of the Burkholderia sp. TJI49 ykkC-I subtype-1 (Guanidine-I) riboswitch aptamer domain. Differing from the previously reported structures of orthologs from Dickeya dadantii and Sulfobacillus acidophilus, our Burkholderia structure reveals a chelated K⁺ ion adjacent to two Mg²⁺ ions in the guanidine-binding pocket. Thermal melting analysis shows that K⁺ chelation, which induces localized conformational changes in the binding pocket, improves guanidinium-RNA interactions. Analysis of ribosome structures suggests that the [K⁺(Mg²⁺)₂] ion triad is uncommon. It is, however, reminiscent of metal ion clusters found in the active sites of ribozymes and DNA polymerases. Previous structural characterization of ykkC-I subtype-2 RNAs, which bind the effector ligands ppGpp and PRPP, indicate that in those paralogs, an adenine responsible for K⁺ chelation in the Burkholderia Guanidine-I riboswitch is replaced by a pyrimidine. This mutation results in a water molecule and Mg²⁺ ion binding in place of the K⁺ ion. Thus, our structural analysis demonstrates how ion and solvent chelation tune divergent ligand specificity and affinity among ykkC-I riboswitches.

Apolar chemical environments compact unfolded RNAs and can promote folding

It is well documented that the structure, and thus function, of nucleic acids depends on the chemical environment surrounding them, which often includes potential proteinaceous binding partners. The nonpolar amino acid side chains of these proteins will invariably alter the polarity of the local chemical environment around the nucleic acid. However, we are only beginning to understand how environmental polarity generally influences the structural and energetic properties of RNA folding. Here, we use a series of aqueous-organic cosolvent mixtures to systematically modulate the solvent polarity around two different RNA folding constructs that can form either secondary or tertiary structural elements. Using single-molecule Förster resonance energy transfer spectroscopy to simultaneously monitor the structural and energetic properties of these RNAs, we show that the unfolded conformations of both model RNAs become more compact in apolar environments characterized by dielectric constants less than that of pure water. In the case of tertiary structure formation, this compaction also gives rise to more energetically favorable folding. We propose that these physical changes arise from an enhanced accumulation of counterions in the low dielectric environment surrounding the unfolded RNA.

High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs

Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.

Non-viral nanoparticles for RNA interference: Principles of design and practical guidelines

Ribonucleic acid interference (RNAi) is an innovative treatment strategy for a myriad of indications. Non-viral synthetic nanoparticles (NPs) have drawn extensive attention as vectors for RNAi due to their potential advantages, including improved safety, high delivery efficiency and economic feasibility. However, the complex natural process of RNAi and the susceptible nature of oligonucleotides render the NPs subject to particular design principles and requirements for practical fabrication. Here, we summarize the requirements and obstacles for fabricating non-viral nano-vectors for efficient RNAi. To address the delivery challenges, we discuss practical guidelines for materials selection and NP synthesis in order to maximize RNA encapsulation efficiency and protection against degradation, and to facilitate the cytosolic release of oligonucleotides. The current status of clinical translation of RNAi-based therapies and further perspectives for reducing the potential side effects are also reviewed.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Contributions and competition of Mg2+ and K+ in folding and stabilization of the Twister ribozyme

Native folded and compact intermediate states of RNA typically involve tertiary structures in the presence of divalent ions such as Mg²⁺ in a background of monovalent ions. In a recent study, we have shown how the presence of Mg²⁺ impacts the transition from partially unfolded to folded states through a "push-pull" mechanism where the ion both favors and disfavors the sampling of specific phosphate-phosphate interactions. To further understand the ion atmosphere of RNA in folded and partially folded states results from atomistic umbrella sampling and oscillating chemical potential grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to obtain atomic-level details of the distributions of Mg²⁺ and K⁺ ions around Twister RNA. Results show the presence of 100 mM Mg²⁺ to lead to increased charge neutralization over that predicted by counterion condensation theory. Upon going from partially unfolded to folded states, overall charge neutralization increases at all studied ion concentrations that, while associated with an increase in the number of direct ion-phosphate interactions, is fully accounted for by the monovalent K⁺ ions. Furthermore, K⁺ preferentially interacts with purine N7 atoms of helical regions in partially unfolded states, thereby potentially stabilizing the helical regions. Thus, both secondary helical structures and formation of tertiary structures leads to increased counterion condensation, thereby stabilizing those structural features of Twister. Notably, it is shown that K⁺ can act as a surrogate for Mg²⁺ by participating in specific interactions with nonsequential phosphate pairs that occur in the folded state, explaining the ability of Twister to self-cleave at submillimolar Mg²⁺ concentrations.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

Multivalent ions and biomolecules: Attempting a comprehensive perspective

Ions are ubiquitous in nature. They play a key role for many biological processes on the molecular scale, from molecular interactions, to mechanical properties, to folding, to self-organisation and assembly, to reaction equilibria, to signalling, to energy and material transport, to recognition etc. Going beyond monovalent ions to multivalent ions, the effects of the ions are frequently not only stronger (due to the obviously higher charge), but qualitatively different. A typical example is the process of binding of multivalent ions, such as Ca2+ , to a macromolecule and the consequences of this ion binding such as compaction, collapse, potential charge inversion and precipitation of the macromolecule. Here we review these effects and phenomena induced by multivalent ions for biological (macro)molecules, from the "atomistic/molecular" local picture of (potentially specific) interactions to the more global picture of phase behaviour including, e. g., crystallisation, phase separation, oligomerisation etc. Rather than attempting an encyclopedic list of systems, we rather aim for an embracing discussion using typical case studies. We try to cover predominantly three main classes: proteins, nucleic acids, and amphiphilic molecules including interface effects. We do not cover in detail, but make some comparisons to, ion channels, colloidal systems, and synthetic polymers. While there are obvious differences in the behaviour of, and the relevance of multivalent ions for, the three main classes of systems, we also point out analogies. Our attempt of a comprehensive discussion is guided by the idea that there are not only important differences and specific phenomena with regard to the effects of multivalent ions on the main systems, but also important similarities. We hope to bridge physico-chemical mechanisms, concepts of soft matter, and biological observations and connect the different communities further.

Environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by atomic force microscopy-based single-molecule force spectroscopy and the implications for advanced polymer materials

"The Tao begets the One. One begets all things of the world." This quote from Tao Te Ching is still inspiring for scientists in chemistry and materials science: The "One" can refer to a single molecule. A macroscopic material is composed of numerous molecules. Although the relationship between the properties of the single molecule and macroscopic material is not well understood yet, it is expected that a deeper understanding of the single-chain mechanics of macromolecules will certainly facilitate the development of materials science. Atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) has been exploited extensively as a powerful tool to study the single-chain behaviors of macromolecules. In this review, we summarize the recent advances in the emerging field of environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by means of AFM-SMFS. First, the single-chain inherent elasticities of several typical linear macromolecules are introduced, which are also confirmed by one of three polymer models with theoretical elasticities of the corresponding macromolecules obtained from quantum mechanical (QM) calculations. Then, the effects of the external environments on the single-chain mechanics of synthetic polymers and biomacromolecules are reviewed. Finally, the impacts of single-chain mechanics of macromolecules on the development of polymer science especially polymer materials are illustrated.

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 Mg2+, 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 Mg2+ 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 Mg2+-RNA interactions that occur along an unfolding pathway of the Twister ribozyme. The free energy profiles reveal stabilization of partially unfolded states by Mg2+, as observed in unfolding experiments, with this stabilization being due to increased sampling of simultaneous interactions of Mg2+ with two or more nonsequential phosphate groups. Notably, these results indicate a push-pull mechanism in which the Mg2+-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 Mg2+-ion interactions contribute to RNA structural properties and stability.

Correction: High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Correction for ‘High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation’ by Hsuan-Lei Sung et al., Phys. Chem. Chem. Phys., 2020, DOI: 10.1039/d0cp01921f.

Single-molecule kinetic studies of DNA hybridization under extreme pressures

Pressure-responsive dynamics of DNA hairpin hybridization/dehybridization is directly visualized at the single molecule level.

Targeting RNA with Small Molecules: Identification of Selective, RNA-Binding Small Molecules Occupying Drug-Like Chemical Space

Although the potential value of RNA as a target for new small molecule therapeutics is becoming increasingly credible, the physicochemical properties required for small molecules to selectively bind to RNA remain relatively unexplored. To investigate the druggability of RNAs with small molecules, we have employed affinity mass spectrometry, using the Automated Ligand Identification System (ALIS), to screen 42 RNAs from a variety of RNA classes, each against an array of chemically diverse drug-like small molecules (~50,000 compounds) and functionally annotated tool compounds (~5100 compounds). The set of RNA–small molecule interactions that was generated was compared with that for protein–small molecule interactions, and naïve Bayesian models were constructed to determine the types of specific chemical properties that bias small molecules toward binding to RNA. This set of RNA-selective chemical features was then used to build an RNA-focused set of ~3800 small molecules that demonstrated increased propensity toward binding the RNA target set. In addition, the data provide an overview of the specific physicochemical properties that help to enable binding to potential RNA targets. This work has increased the understanding of the chemical properties that are involved in small molecule binding to RNA, and the methodology used here is generally applicable to RNA-focused drug discovery efforts.

Finely tuned conformational dynamics regulate the protective function of the lncRNA MALAT1 triple helix

Nucleic acid triplexes may regulate many important biological processes. Persistent accumulation of the oncogenic 7-kb long noncoding RNA MALAT1 is dependent on an unusually long intramolecular triple helix. This triplex structure is positioned within a conserved ENE (element for nuclear expression) motif at the lncRNA 3′ terminus and protects the entire transcript from degradation in a polyA-independent manner. A requisite 3′ maturation step leads to triplex formation though the precise mechanism of triplex folding remains unclear. Furthermore, the contributions of several peripheral structural elements to triplex formation and protective function have not been determined. We evaluated the stability, conformational fluctuations, and function of this MALAT1 ENE triple helix (M1^TH) protective element using in vitro mutational analyses coupled with biochemical and biophysical characterizations. Using fluorescence and UV melts, FRET, and an exonucleolytic decay assay we define a concerted mechanism for triplex formation and uncover a metastable, dynamic triplex population under near-physiological conditions. Structural elements surrounding the triplex regulate the dynamic M1^TH conformational variability, but increased triplex dynamics lead to M1^TH degradation. Taken together, we suggest that finely tuned dynamics may be a general mechanism regulating triplex-mediated functions.

Divalent ions tune the kinetics of a bacterial GTPase center rRNA folding transition from secondary to tertiary structure

Folding of an RNA from secondary to tertiary structure often depends on divalent ions for efficient electrostatic charge screening (nonspecific association) or binding (specific association). To measure how different divalent cations modify folding kinetics of the 60 nucleotide Ecoli rRNA GTPase center, we combined stopped-flow fluorescence in the presence of Mg²⁺, Ca²⁺, or Sr²⁺ together with time-resolved small angle X-ray scattering (SAXS) in the presence of Mg²⁺ to observe the folding process. Immediately upon addition of each divalent ion, the RNA undergoes a transition from an extended state with secondary structure to a more compact structure. Subsequently, specific divalent ions modulate populations of intermediates in conformational ensembles along the folding pathway with transition times longer than 10 msec. Rate constants for the five folding transitions act on timescales from submillisecond to tens of seconds. The sensitivity of RNA tertiary structure to divalent cation identity affects all but the fastest events in RNA folding, and allowed us to identify those states that prefer Mg²⁺ The GTPase center RNA appears to have optimized its folding trajectory to specifically utilize this most abundant intracellular divalent ion.

The structural diversity of C-rich DNA aggregates: unusual self-assembly of beetle-like nanostructures

We studied the ability of oligonucleotides C_nT₂₅ (n = 2, 5, 7, 9, 12, 25) to form an intermolecular i-motif using circular dichroism, ultra-violet spectroscopy, nuclear magnetic resonance, high-resolution atomic force microscopy, high-performance liquid chromatography, and molecular dynamics simulations. The arrangement of single-stranded oligonucleotides in multimer i-motifs was very unusual: C-tracts of different oligonucleotides followed each other consecutively in order to fold into a closed intermolecular i-motif core with minimal loops (one cytidine in a loop spanning over a minor groove, three cytidines in a loop over a major groove); intact T-tracts protruded from predefined loci allowing visualization of beetle-like nanostructures by atomic force microscopy. The same structures were formed from analogous biotinylated oligonucleotides demonstrating one of the potential applications of such structures as carriers of multiple functional groups. Our findings open up possibilities for the rational design of pH-sensitive DNA aggregates and evaluation of the efficiency of their assembly.

Interaction of Na+, K+, Mg2+ and Ca2+ counter cations with RNA

Data on the location of alkaline and alkaline earth ions at RNA from crystallography, spectroscopy and computational modeling are reviewed.

All-atom MD indicates ion-dependent behavior of therapeutic DNA polymer

Understanding the efficacy of and creating delivery mechanisms for therapeutic nucleic acids requires understanding structural and kinetic properties which allow these polymers to promote the death of cancerous cells. One molecule of interest is a 10 mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) - also called F10. Here we investigate the structural and kinetic behavior of F10 in intracellular and extracellular solvent conditions along with non-biological conditions that may be efficacious in in vitro preparations of F10 delivery systems. From our all-atom molecular dynamics simulations totaling 80 microseconds, we predict that F10's phosphate groups form close-range interactions with calcium and zinc ions, with calcium having the highest affinity of the five ions investigated. We also predict that F10's interactions with magnesium, potassium and sodium are almost exclusively long-range interactions. In terms of intramolecular interactions, we find that F10 is least structured (in terms of hydrogen bonds among bases) in the 150 mM NaCl (extracellular-like solvent conditions) and most structured in 150 mM ZnCl2. Kinetically, we see that F10 is unstable in the presence of magnesium, sodium or potassium, finding stable kinetic traps in the presence of calcium or zinc.

Computational Assessment of Potassium and Magnesium Ion Binding to a Buried Pocket in GTPase-Associating Center RNA

An experimentally well-studied model of RNA tertiary structures is a 58mer rRNA fragment, known as GTPase-associating center (GAC) RNA, in which a highly negative pocket walled by phosphate oxygen atoms is stabilized by a chelated cation. Although such deep pockets with more than one direct phosphate to ion chelation site normally include magnesium, as shown in one GAC crystal structure, another GAC crystal structure and solution experiments suggest potassium at this site. Both crystal structures also depict two magnesium ions directly bound to the phosphate groups comprising this controversial pocket. Here, we used classical molecular dynamics simulations as well as umbrella sampling to investigate the possibility of binding of potassium versus magnesium inside the pocket and to better characterize the chelation of one of the binding magnesium ions outside the pocket. The results support the preference of the pocket to accommodate potassium rather than magnesium and suggest that one of the closely binding magnesium ions can only bind at high magnesium concentrations, such as might be present during crystallization. This work illustrates the complementary utility of molecular modeling approaches with atomic-level detail in resolving discrepancies between conflicting experimental results.

Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding

The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg²⁺ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl₂ were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.

Enhanced deoxyribozyme-catalyzed RNA ligation in the presence of organic cosolvents

Deoxyribozyme and aptamer selections are typically conducted in aqueous buffer solutions. Using nonaqueous cosolvents in selection experiments will help expand the activity of deoxyribozymes with non-oligonucleotide substrates and will allow identification of new aptamers for nonprotein targets. We undertook in vitro selections utilizing a small amount of methanol in the reaction to keep the herbicides alachlor and atrazine in solution with the goal of identifying deoxyribozymes that require these herbicides for activity. The resulting deoxyribozymes successfully catalyze RNA ligation, but do not require alachlor or atrazine. Surprisingly, some of these deoxyribozymes displayed better catalytic activity in the presence of methanol over just aqueous buffer. We investigated several organic cosolvents to see if this enhancement was limited to methanol and found that other cosolvents, including ethanol, DMSO, and DMF, supported activity; in some cases, greater enhancement was observed. On the basis of these results, we tested two other previously identified RNA-ligating deoxyribozymes to assess their tolerance of cosolvents and determined that different deoxyribozymes showed different responses to the cosolvents. Our results demonstrate that deoxyribozymes can tolerate and, in some cases, display enhanced activity in alternative solvent conditions. These findings will facilitate the development of responsive deoxyribozyme systems utilizing components with limited water solubility.

The structural stability and catalytic activity of DNA and RNA oligonucleotides in the presence of organic solvents

Organic solvents and apolar media are used in the studies of nucleic acids to modify the conformation and function of nucleic acids, to improve solubility of hydrophobic ligands, to construct molecular scaffolds for organic synthesis, and to study molecular crowding effects. Understanding how organic solvents affect nucleic acid interactions and identifying the factors that dominate solvent effects are important for the creation of oligonucleotide-based technologies. This review describes the structural and catalytic properties of DNA and RNA oligonucleotides in organic solutions and in aqueous solutions with organic cosolvents. There are several possible mechanisms underlying the effects of organic solvents on nucleic acid interactions. The reported results emphasize the significance of the osmotic pressure effect and the dielectric constant effect in addition to specific interactions with nucleic acid strands. This review will serve as a guide for the selection of solvent systems based on the purpose of the nucleic acid-based experiments.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Collapse and hybridization of RNA: view from replica technique approach

The replica technique method is applied to investigate the kinetic behavior of the coarse-grained model for the RNA molecule. A non-equilibrium phase transition of second order between the glassy phase and the ensemble of freely fluctuating structures has been observed. The non-equilibrium steady state is investigated as well and the thermodynamic characteristics of the system have been evaluated. The non-equilibrium behavior of the specific heat is discussed. Based on our analysis, we point out the state in the kinetic pathway in which the RNA molecule is most prone to hybridization.

The Globular State of the Single-Stranded RNA: Effect of the Secondary Structure Rearrangements

The mutual influence of the slow rearrangements of secondary structure and fast collapse of the long single-stranded RNA (ssRNA) in approximation of coarse-grained model is studied with analytic calculations. It is assumed that the characteristic time of the secondary structure rearrangement is much longer than that for the formation of the tertiary structure. A nonequilibrium phase transition of the 2nd order has been observed.

Formation of Tertiary Interactions during rRNA GTPase Center Folding

The 60-nt GTPase center (GAC) of 23S rRNA has a phylogenetically conserved secondary structure with two hairpin loops and a 3-way junction. It folds into an intricate tertiary structure upon addition of Mg(2+) ions, which is stabilized by the L11 protein in cocrystal structures. Here, we monitor the kinetics of its tertiary folding and Mg(2+)-dependent intermediate states by observing selected nucleobases that contribute specific interactions to the GAC tertiary structure in the cocrystals. The fluorescent nucleobase 2-aminopurine replaced three individual adenines, two of which make long-range stacking interactions and one that also forms hydrogen bonds. Each site reveals a unique response to Mg(2+) addition and temperature, reflecting its environmental change from secondary to tertiary structure. Stopped-flow fluorescence experiments revealed that kinetics of tertiary structure formation upon addition of MgCl2 are also site specific, with local conformational changes occurring from 5 ms to 4s and with global folding from 1 to 5s. Site-specific substitution with (15)N-nucleobases allowed observation of stable hydrogen bond formation by NMR experiments. Equilibrium titration experiments indicate that a stable folding intermediate is present at stoichiometric concentrations of Mg(2+) and suggest that there are two initial sites of Mg(2+) ion association.

Effects of Water on the Single-Chain Elasticity of Poly(U) RNA

Water, the dominant component under the physiological condition, is a complicated solvent which greatly affects the properties of solute molecules. Here, we utilize atomic force microscope-based single-molecule force spectroscopy to study the influence of water on the single-molecule elasticity of an unstructured single-stranded RNA (poly(U)). In nonpolar solvents, RNA presents its inherent elasticity, which is consistent with the theoretical single-chain elasticity calculated by quantum mechanics calculations. In aqueous buffers, however, an additional energy of 1.88 kJ/mol·base is needed for the stretching of the ssRNA chain. This energy is consumed by the bound water rearrangement (Ew) during chain elongation. Further experimental results indicate that the Ew value is uncorrelated to the salt concentrations and stretching velocity. The results obtained in an 8 M guanidine·HCl solution provide more evidence that the bound water molecules around RNA give rise to the observed deviation between aqueous and nonaqueous environments. Compared to synthetic water-soluble polymers, the value of Ew of RNA is much lower. The weak interference of water is supposed to be the precondition for the RNA secondary structure to exist in aqueous solution.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

Structural and energetic analysis of 2-aminobenzimidazole inhibitors in complex with the hepatitis C virus IRES RNA using molecular dynamics simulations

Despite the many biological functions of RNA, very few drugs have been designed or found to target RNA. Here we report the results of molecular dynamics (MD) simulations and binding energy analyses on hepatitis C virus internal ribosome entry site (IRES) RNA in complex with highly charged 2-aminobenzimidazole inhibitors. Initial coordinates were taken from NMR and crystallography studies that had yielded different binding modes. During MD simulations, the RNA-inhibitor complex is stable in the crystal conformation but not in the NMR conformation. Additionally, we found that existing and standard MD trajectory postprocessing free energy methods, such as the MM-GBSA and MM-PBSA approaches available in AMBER, seem unsuitable to properly rank the binding energies of complexes between highly charged molecules. A better correlation with the experimental data was found using a rather simple binding enthalpy calculation based on the explicitly solvated potential energies. In anticipation of further growth in the use of small molecules to target RNA, we include results addressing the impact of charge assignment on docking, the structural role of magnesium in the IRES-inhibitor complex, the entropic contribution to binding energy, and simulations of a plausible scaffold design for new inhibitors.

High efficiency acetylcholinesterase immobilization on DNA aptamer modified surfaces

We report here the in vitro selection of DNA aptamers for electric eel acetylcholinesterase (AChE). One selected aptamer sequence (R15/19) has a high affinity towards the enzyme (Kd=157±42 pM). Characterization of the aptamer showed its binding is not affected by low ionic strength (~20 mM), however significant reduction in affinity occurred at high ionic strength (~1.2 M). In addition, this aptamer does not inhibit the catalytic activity of AChE that we exploit through immobilization of the DNA on a streptavidin-coated surface. Subsequent immobilization of AChE by the aptamer results in a 4-fold higher catalytic activity when compared to adsorption directly on to plastic.

Understanding nucleic acid-ion interactions

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

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(2+) 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(2+) (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.

Monitoring global structural changes and specific metal-ion-binding sites in RNA by in-line probing and Tb(III) cleavage

In this chapter we describe the use of two methods, in-line probing as well as terbium(III) cleavage. Both methods can be applied to RNAs of any size, structure, and function. Aside from revealing directly metal ion-binding sites these techniques also provide structural information for longer RNA sequences that are out of range to be analyzed with other techniques such as NMR. The cleavage pattern derived from in-line probing experiments reflects local and overall conformational changes in RNA upon the addition of metal ions, metal complexes, or other ligands. On the other side, terbium(III) cleavage experiments are applied to locate specific metal ion-binding sites in RNA molecules.

Noncanonical structures and their thermodynamics of DNA and RNA under molecular crowding: beyond the Watson-Crick double helix

How does molecular crowding affect the stability of nucleic acid structures inside cells? Water is the major solvent component in living cells, and the properties of water in the highly crowded media inside cells differ from that in buffered solution. As it is difficult to measure the thermodynamic behavior of nucleic acids in cells directly and quantitatively, we recently developed a cell-mimicking system using cosolutes as crowding reagents. The influences of molecular crowding on the structures and thermodynamics of various nucleic acid sequences have been reported. In this chapter, we discuss how the structures and thermodynamic properties of nucleic acids differ under various conditions such as highly crowded environments, compartment environments, and in the presence of ionic liquids, and the major determinants of the crowding effects on nucleic acids are discussed. The effects of molecular crowding on the activities of ribozymes and riboswitches on noncanonical structures of DNA- and RNA-like quadruplexes that play important roles in transcription and translation are also described.

Genus distribution and thermodynamics of a random matrix model of RNA with Penner interaction

The nonlinear Penner external interaction is introduced and studied in the random matrix model of homo RNA. A numerical technique is developed to study the partition function, and a general formula is obtained for all lengths. The genus distribution function for the system is obtained, plotted, and compared with the genus distribution for the real RNA structures found from the protein databank. The genus distribution shows that the nonlinear interaction favors the formation of low genus structures and matches the result for real RNA structures. The distribution of structure with temperature suggests that nonlinear interaction is biased toward the planar structures. The variation of chemical potential with temperature and interaction strength indicates the presence of additional molecules in the system other than the magnesium ions and possibly represents a phase transition. The specific heat has a bump and its derivatives shows a double-peak behavior at a particular temperature. On analyzing the specific heat and derivatives for each genus separately, the planar structure (genus zero) is shown to contribute the most to the bump and double peak. This observation in the nonlinear model is similar to that observed in the unfolding experiments on RNA.

Comparison of interactions of diamine and Mg²⁺ with RNA tertiary structures: similar versus differential effects on the stabilities of diverse RNA folds

Cations play a large role in stabilizing the native state of RNA in vivo. In addition to Mg²⁺, putrescine²⁺ is an abundant divalent cation in bacterial cells, but its effect on the folding of RNA tertiary structure has not been widely explored. In this study, we look at how the stabilities of four structured RNAs, each with a different degree of dependence on K⁺ and Mg²⁺, are affected by putrescine²⁺ relative to Mg²⁺. Through the use of thermal melts, we observe that (i) at a given concentration, putrescine²⁺ is less effective than Mg²⁺ at stabilizing RNA, (ii) the stability imparted to RNA by various diamines is a function of charge density (average separation distance between charges) as well as the flexibility of the counterion, and (iii) when Mg²⁺ is already present in a buffer, further addition of putrescine²⁺ may either destabilize or stabilize RNA structure, depending on whether the native RNA does or does not chelate Mg²⁺ ion, respectively. At ion concentrations likely to be found in vivo, the effect of putrescine²⁺ on the free energy of folding of an RNA tertiary structure is probably quite small compared to that of Mg²⁺, but the ability of mixed Mg²⁺/putrescine²⁺ environments to (in effect) discriminate between different RNA architectures suggests that, in some cells, the evolution of functional RNA structures may have been influenced by the presence of putrescine²⁺.

RNA secondary structure prediction using soft computing

Prediction of RNA structure is invaluable in creating new drugs and understanding genetic diseases. Several deterministic algorithms and soft computing-based techniques have been developed for more than a decade to determine the structure from a known RNA sequence. Soft computing gained importance with the need to get approximate solutions for RNA sequences by considering the issues related with kinetic effects, cotranscriptional folding, and estimation of certain energy parameters. A brief description of some of the soft computing-based techniques, developed for RNA secondary structure prediction, is presented along with their relevance. The basic concepts of RNA and its different structural elements like helix, bulge, hairpin loop, internal loop, and multiloop are described. These are followed by different methodologies, employing genetic algorithms, artificial neural networks, and fuzzy logic. The role of various metaheuristics, like simulated annealing, particle swarm optimization, ant colony optimization, and tabu search is also discussed. A relative comparison among different techniques, in predicting 12 known RNA secondary structures, is presented, as an example. Future challenging issues are then mentioned.

Open questions on the origin of life at anoxic geothermal fields

We have recently reconstructed the 'hatcheries' of the first cells by combining geochemical analysis with phylogenomic scrutiny of the inorganic ion requirements of universal components of modern cells (Mulkidjanian et al. Proc Natl Acad Sci U S A 109:E821-830, 2012). These ubiquitous, and by inference primordial, proteins and functional systems show affinity to and functional requirement for K⁺, Zn²⁺, Mn²⁺, and phosphate. Thus, protocells must have evolved in habitats with a high K⁺/Na⁺ ratio and relatively high concentrations of Zn, Mn and phosphorous compounds. Geochemical reconstruction shows that the ionic composition conducive to the origin of cells could not have existed in marine settings but is compatible with emissions of vapor-dominated zones of inland geothermal systems. Under an anoxic, CO₂-dominated atmosphere, the ionic composition of pools of cool, condensed vapor at anoxic geothermal fields would resemble the internal milieu of modern cells. Such pools would be lined with porous silicate minerals mixed with metal sulfides and enriched in K⁺ ions and phosphorous compounds. Here we address some questions that have appeared in print after the publication of our anoxic geothermal field scenario. We argue that anoxic geothermal fields, which were identified as likely cradles of life by using a top-down approach and phylogenomics analysis, could provide geochemical conditions similar to those which were suggested as most conducive for the emergence of life by the chemists who pursuit the complementary bottom-up strategy.