Salt dependence of nucleic acid hairpin stability

Visualized lateral flow assay for logic determination of co-existing viral RNA fragments

Different types of pathogenic viruses that have common transmission path can be co-infected, inducing distinct disease procession in comparison to that infection of one. Also, in the post COVID-19 time, more types of respiratory infectious virus are becoming prevalent and are concurrent. Those bring an urgent need for detection of co-existing viruses. Here, we propose a visualized lateral flow assay for logic determination of co-existing viral RNA fragments. In the presence of specific viral RNA inputs, DNAzyme is de-blocked according to defined logic, and catalyzes the hydrolysis of hairpin-structural substrate. One of cleaved substrates contains DNAzyme domain to realize dual signal amplification, which obtains copious of the other cleaved substrates. The cleaved substrates act as linking strands for bridging DNA-modified gold nanoparticles onto lateral flow strip to induce coloration on test line. "AND", "OR" and "INHIBIT" controlled lateral flow assays are respectively demonstrated for co-existing viral RNA detection, and the visual results can be obtained by the same kind of prepared strip, without need of re-fabricating strips according to logic systems. The work provides a flexible, convenient, visual and logic-processing strategy for simultaneous analysis of co-existing viruses.

Mono-valent salt corrections for RNA secondary structures in the ViennaRNA package

BACKGROUND

RNA features a highly negatively charged phosphate backbone that attracts a cloud of counter-ions that reduce the electrostatic repulsion in a concentration dependent manner. Ion concentrations thus have a large influence on folding and stability of RNA structures. Despite their well-documented effects, salt effects are not handled consistently by currently available secondary structure prediction algorithms. Combining Debye-Hückel potentials for line charges and Manning's counter-ion condensation theory, Einert et al. (Biophys J 100: 2745-2753, 2011) modeled the energetic contributions of monovalent cations on loops and helices.

RESULTS

The model of Einert et al. is adapted to match the structure of the dynamic programming recursion of RNA secondary structure prediction algorithms. An empirical term describing the salt dependence of the duplex initiation energy is added to improve co-folding predictions for two or more RNA strands. The slightly modified model is implemented in the ViennaRNA package in such way that only the energy parameters but not the algorithmic structure is affected. A comparison with data from the literature show that predicted free energies and melting temperatures are in reasonable agreement with experiments.

CONCLUSION

The new feature in the ViennaRNA package makes it possible to study effects of salt concentrations on RNA folding in a systematic manner. Strictly speaking, the model pertains only to mono-valent cations, and thus covers the most important parameter, i.e., the NaCl concentration. It remains a question for future research to what extent unspecific effects of bi- and tri-valent cations can be approximated in a similar manner.

AVAILABILITY

Corrections for the concentration of monovalent cations are available in the ViennaRNA package starting from version 2.6.0.

Structural properties and binding mechanism of DNA aptamers sensing saliva melatonin for diagnosis and monitoring of circadian clock and sleep disorders

Circadian desynchrony with the external light-dark cycle influences the rhythmic secretion of melatonin which is among the first signs of circadian rhythm sleep disorders. An accurate dim light melatonin onset (established indicator of circadian rhythm sleep disorders) measurement requires lengthy assays, and antibody affinities alterations, especially in patients with circadian rhythm disorders whose melatonin salivary levels vary significantly, making antibodies detection mostly inadequate. In contrast, aptamers with their numerous advantages (e.g., target selectivity, structural flexibility in tuning binding affinities, small size, etc.) can become preferable biorecognition molecules for salivary melatonin detection with high sensitivity and specificity. This study thoroughly characterizes the structural property and binding mechanism of a single-stranded DNA aptamer full sequence (MLT-C-1) and its truncated versions (MLT-A-2, MLT-A-4) to decipher its optimal characteristics for saliva melatonin detection. We use circular dichroism spectroscopy to determine aptamers' conformational changes under different ionic strengths and showed that aptamers display a hairpin loop structure where few base pairs in the stem play a significant role in melatonin binding and formation of aptamer stabilized structure. Through microscale thermophoresis, aptamers demonstrated a high binding affinity in saliva samples (MLT-C-1F K_d = 12.5 ± 1.7 nM; MLT-A-4F K_d = 11.2 ± 1.6 nM; MLT-A-2F K_d = 2.4 ± 2.8 nM; limit-of-detection achieved in pM, highest sensitivity attained for MLT-A-2F aptamer with the lowest detection limit of 1.35 pM). Our data suggest that aptamers are promising as biorecognition molecules and provide the baseline parameters for the development of an aptamer-based point-of-care diagnostic system for melatonin detection and accurate profiling of its fluctuations in saliva.

Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions

RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.

Aptamer-antibody hybrid ELONA that uses hybridization chain reaction to detect a urinary biomarker EN2 for bladder and prostate cancer

We report an EN2-specific (K_d = 8.26 nM) aptamer, and a sensitive and specific enzyme-linked oligonucleotide assay (ELONA) for rapid and sensitive colorimetric detection of bladder and prostate cancer biomarker EN2 in urine. The assay relies on an aptamer-mediated hybridization chain reaction (HCR) to generate DNA nanostructures that bind to EN2 and simultaneously amplify signals. The assay can be performed within 2.5 h, and has a limit of detection of 0.34 nM in buffer and 2.69 nM in artificial urine. Moreover, this assay showed high specificity as it did not detect other urinary proteins, including biomarkers of other cancers. The proposed ELONA is inexpensive, highly reproducible, and has great chemical stability, so it may enable development of a simple, sensitive and accurate diagnostic tool to detect bladder and prostate cancers early.

Evaluation of Weak Interactions of Proteins and Organic Cations with DNA Duplex Structures

Molecular interactions and reactions in living cells occur with high background concentrations of organic compounds including proteins. Uncharged water-soluble polymers are commonly used cosolutes in studies on molecular crowding, and most studies argue about the effects of intracellular crowding based on results obtained using polymer cosolutes. Further investigations using protein crowders and organic cations are important in understanding the effects of cellular environments on nucleic acids with negatively charged surfaces. We assessed the effects of using model globular proteins, serum proteins, histone proteins, structurally flexible polypeptides, di- and polyamines, and uncharged polymers. Thermal stability analysis of DNA oligonucleotide structures revealed that unlike conventional polymer cosolutes, basic globular proteins (lysozyme and cytochrome c) at high concentrations stabilized long internal and bulge loop structures but not fully matched duplexes. The selective stabilization of long loop structures suggests preferential binding to unpaired nucleotides in loops through weak electrostatic interactions. Furthermore, the ability of the proteins to stabilize the loop structures was enhanced under macromolecular crowding conditions. Remarkably, the effects of basic proteins on the stability of fully matched duplexes were dissimilar to those of basic amino-acid-rich polypeptides and polyamines. This study provides new insights into the interaction of nucleic acid structures with organic cations.

The effect of hairpin loop on the structure and gene expression activity of the long-loop G-quadruplex

G-quadruplex (G4), a four-stranded DNA or RNA structure containing stacks of guanine tetrads, plays regulatory roles in many cellular functions. So far, conventional G4s containing loops of 1–7 nucleotides have been widely studied. Increasing experimental evidence suggests that unconventional G4s, such as G4s containing long loops (long-loop G4s), play a regulatory role in the genome by forming a stable structure. Other secondary structures such as hairpins in the loop might thus contribute to the stability of long-loop G4s. Therefore, investigation of the effect of the hairpin-loops on the structure and function of G4s is required. In this study, we performed a systematic biochemical investigation of model G4s containing long loops with various sizes and structures. We found that the long-loop G4s are less stable than conventional G4s, but their stability increased when the loop forms a hairpin (hairpin-G4). We also verified the biological significance of hairpin-G4s by showing that hairpin-G4s present in the genome also form stable G4s and regulate gene expression as confirmed by in cellulo reporter assays. This study contributes to expanding the scope and diversity of G4s, thus facilitating future studies on the role of G4s in the human genome.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

The ability to accurately predict RNA hairpin structure and stability for different loop sequences and salt conditions is important for understanding, modeling, and designing larger RNA folds. However, traditional RNA secondary structure models cannot treat loop-sequence and ionic effects on RNA hairpin folding. Here, we describe a general, three-dimensional (3D) conformation-based computational method for modeling salt concentration-dependent conformational distributions and the detailed 3D structures for a set of three RNA hairpins that contain a variable, 15-nucleotide loop sequence. For a given RNA sequence, the new, to our knowledge, method integrates a Vfold2D two-dimensional structure folding model with IsRNA coarse-grained molecular dynamics 3D folding simulations and Monte Carlo tightly bound ion estimations of ion-mediated electrostatic interactions. The model predicts free-energy landscapes for the different RNA hairpin-forming sequences with variable salt conditions. The theoretically predicted results agree with the experimental fluorescence measurements, validating the strategy. Furthermore, the theoretical model goes beyond the experimental results by enabling in-depth 3D structural analysis, revealing energetic mechanisms for the sequence- and salt-dependent folding stability. Although the computational framework presented here is developed for RNA hairpin systems, the general method may be applied to investigate other RNA systems, such as multiway junctions or pseudoknots in mixed metal ion solutions.

Molecular dynamics simulation of the binding process of ligands to the add adenine riboswitch aptamer

Riboswitches are RNA-structured elements that modulate gene expression through changing their conformation in response to specific metabolite binding. However, the regulation mechanisms of riboswitches by ligand binding are still not well understood. At present two possible ligand-binding mechanisms have been proposed: conformational selection and induced fit. Based on explicit-solvent molecular dynamics (MD) simulations, we have studied the process of the binding of ligands (adenines) to add adenine riboswitch aptamer (AARA) in detail. Our results show that the relative high flexibility of the junction J23 of AARA allows the ligands to be captured by the binding pocket of AARA in a near-native state, which may be driven by hydrophobic and base-stacking interactions. In addition, the binding of a ligand makes the stem P1 and the junction J23 of AARA more stable than in the absence of the ligand. Generally, our analyses show that the ligand-binding process of the add adenine riboswitch may be partially embodied by a conformational selection mechanism.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

3D structure stability of the HIV-1 TAR RNA in ion solutions: A coarse-grained model study

As an extremely common structural motif, RNA hairpins with bulge loops [e.g., the human immunodeficiency virus type 1 (HIV-1) transactivation response (TAR) RNA] can play essential roles in normal cellular processes by binding to proteins and small ligands, which could be very dependent on their three-dimensional (3D) structures and stability. Although the structures and conformational dynamics of the HIV-1 TAR RNA have been extensively studied, there are few investigations on the thermodynamic stability of the TAR RNA, especially in ion solutions, and the existing studies also have some divergence on the unfolding process of the RNA. Here, we employed our previously developed coarse-grained model with implicit salt to predict the 3D structure, stability, and unfolding pathway for the HIV-1 TAR RNA over a wide range of ion concentrations. As compared with the extensive experimental/theoretical results, the present model can give reliable predictions on the 3D structure stability of the TAR RNA from the sequence. Based on the predictions, our further comprehensive analyses on the stability of the TAR RNA as well as its variants revealed that the unfolding pathway of an RNA hairpin with a bulge loop is mainly determined by the relative stability between different states (folded state, intermediate state, and unfolded state) and the strength of the coaxial stacking between two stems in folded structures, both of which can be apparently modulated by the ion concentrations as well as the sequences.

Predicting Monovalent Ion Correlation Effects in Nucleic Acids

Ion correlation and fluctuation can play a potentially significant role in metal ion-nucleic acid interactions. Previous studies have focused on the effects for multivalent cations. However, the correlation and fluctuation effects can be important also for monovalent cations around the nucleic acid surface. Here, we report a model, gMCTBI, that can explicitly treat discrete distributions of both monovalent and multivalent cations and can account for the correlation and fluctuation effects for the cations in the solution. The gMCTBI model enables investigation of the global ion binding properties as well as the detailed discrete distributions of the bound ions. Accounting for the ion correlation effect for monovalent ions can lead to more accurate predictions, especially in a mixed monovalent and multivalent salt solution, for the number and location of the bound ions. Furthermore, although the monovalent ion-mediated correlation does not show a significant effect on the number of bound ions, the correlation may enhance the accumulation of monovalent ions near the nucleic acid surface and hence affect the ion distribution. The study further reveals novel ion correlation-induced effects in the competition between the different cations around nucleic acids.

A quencher-free DNAzyme beacon for fluorescently sensing uranyl ions via embedding 2-aminopurine

DNAzyme-based fluorescent probes have provided valuable protocols for detecting uranium, one of the most common radioactive contaminants in the environment, with ultra-high selectivity and sensitivity. Designing novel DNAzyme beacons to update the mode of fluorescence reporting and/or quenching will continuously enhance "turn-on" sensing performance as well as promote actual application of the biological probes. In this work, we developed a novel quencher-free DNAzyme beacon by embedding fluorescent 2-aminopurine for rapid detection of uranyl ion. 2-aminopurine is able to substitute adenine and keep strong fluorescence in single-stranded DNA whereas being quenched in the hybridized double-stranded DNA by the base-stacking interaction. The combination of such trait of 2-aminopurine and cleavage reaction of DNAzyme in the presence of target co-factors possesses two main advantages for ion sensing: simplicity for avoidance of extra quencher groups and high performance because of superiority of DNAzyme essence. The experimental conditions including embedding site, pH and salt concentration of buffer solutions, and the amount ratio of enzyme strand to substrate strand used to form DNAzymes were systematically optimized to inspire the highest performance of the biological beacon. Thus, a detection limit of 9.6 nM, a wide linear range from 5 nM to 400 nM (R² = 0.997), and selectivity of more than 400 000-fold over other metal ions were achieved by the novel DNAzyme probes. The highly sensitive, selective and quencher-free DNAzyme probes accommodated a simple and cost-efficient alternative to current fluorescent counterparts, holding a great potential for further application in practical ion assay.

Predicting RNA-Metal Ion Binding with Ion Dehydration Effects

Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg²⁺ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.

Molecular Simulations of Ion Effects on the Thermodynamics of RNA Folding

How ions affect RNA folding thermodynamics and kinetics is an important but a vexing problem that remains unsolved. Experiments have shown that the free-energy change, Δ G( c), of RNA upon folding varies with the salt concentration ( c) as, Δ G( c) = k _c ln  c + const, where the coefficient k _c is proportional to the difference in the ion preferential coefficient, ΔΓ. We performed simulations of a coarse-grained model, by modeling electrostatic interactions implicitly and with explicit representation of ions, to elucidate the molecular underpinnings of the relationship between Δ G and ΔΓ. The simulations quantitatively reproduce the heat capacity for a pseudoknot, thus validating the model. We show that Δ G( c), calculated directly from ΔΓ, varies linearly with ln  c ( c < 0.2 M), for a hairpin and the pseudoknot, demonstrating a molecular link between the two quantities. Explicit ion simulations also show the linear dependence of Δ G( c) on ln  c at all c with k _c = 2 k_B T, except that Δ G( c) values are shifted by ∼2 kcal/mol higher than experiments. The discrepancy is due to an underestimation of Γ for both the folded and unfolded states while giving accurate values for ΔΓ. The predictions for the salt dependence of ΔΓ are amenable to test using single-molecule pulling experiments. The framework provided here can be used to obtain accurate thermodynamics for other RNA molecules as well.

Amino Acid Stabilization of Nucleic Acid Secondary Structure: Kinetic Insights from Single-Molecule Studies

Amino acid and nucleic acid interactions are central in biology and may have played a role in the evolutionary development of protein-based life from an early "RNA Universe." To explore the possible role of single amino acids in promoting nucleic acid folding, single-molecule Förster resonance energy transfer experiments have been implemented with a DNA hairpin construct (7 nucleotide double strand with a 40A loop) as a simple model for secondary structure formation. Exposure to positively charged amino acids (arginine and lysine) is found to clearly stabilize the secondary structure. Kinetically, each amino acid promotes folding by generating a large increase in the folding rate with little change in the unfolding rate. From analysis as a function of temperature, arginine and lysine are found to significantly increase the overall exothermicity of folding while imposing only a small entropic penalty on the folding process. Detailed investigations into the kinetics and thermodynamics of this amino acid-induced folding stability reveal arginine and lysine to interact with nucleic acids in a manner reminiscent of monovalent cations. Specifically, these observations are interpreted in the context of an ion atmosphere surrounding the nucleic acid, in which amino acid salts stabilize folding qualitatively like small monovalent cations but also exhibit differences because of the composition of their side chains.

Control of the polyamine biosynthesis pathway by G2-quadruplexes

G-quadruplexes are naturally-occurring structures found in RNAs and DNAs. Regular RNA G-quadruplexes are highly stable due to stacked planar arrangements connected by short loops. However, reports of irregular quadruplex structures are increasing and recent genome-wide studies suggest that they influence gene expression. We have investigated a grouping of G2-motifs in the UTRs of eight genes involved in polyamine biosynthesis, and concluded that several likely form novel metastable RNA G-quadruplexes. We performed a comprehensive biophysical characterization of their properties, comparing them to a reference G-quadruplex. Using cellular assays, together with polyamine-depleting and quadruplex-stabilizing ligands, we discovered how some of these motifs regulate and sense polyamine levels, creating feedback loops during polyamine biosynthesis. Using high-resolution 1H-NMR spectroscopy, we demonstrated that a long-looped quadruplex in the AZIN1 mRNA co-exists in salt-dependent equilibria with a hairpin structure. This study expands the repertoire of regulatory G-quadruplexes and demonstrates how they act in unison to control metabolite homeostasis.

Detection of unique Ebola virus oligonucleotides using fluorescently-labeled phosphorodiamidate morpholino oligonucleotide probe pairs

Here we identify a low-cost diagnostic platform using fluorescently-labeled phosphorodiamidate morpholino oligonucleotide (PMO) probe pairs, which upon binding target oligonucleotides undergo fluorescence resonance energy transfer (FRET). Using a target oligonucleotide derived from the Ebola virus (EBOV), we have derivatized PMO probes with either Alexa Fluor488 (donor) or tetramethylrhodamine (acceptor). Upon EBOV target oligonulceotide binding, observed changes in FRET between PMO probe pairs permit a 25 pM lower limit of detection; there is no off-target binding within a complex mixture of nucleic acids and other biomolecules present in human saliva. Equivalent levels of FRET occur using PMO probe pairs for single or double stranded oligonucleotide targets. High-affinity binding is retained under low-ionic strength conditions that disrupt oligonucleotide secondary structures (e.g., stem-loop structures), ensuring reliable target detection. Under these low-ionic strength conditions, rates of PMO probe binding to target oligonucleotides are increased 3-fold relative to conventional high-ionic strength conditions used for nucleic acid hybridization, with half-maximal binding occurring within 10 min. Our results indicate an ability to use PMO probe pairs to detect clinically relevant levels of EBOV and other oligonucleotide targets in complex biological samples without the need for nucleic acid amplification, and open the possibility of population screening that includes assays for the genomic integration of DNA based copies of viral RNA.

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration.

The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

RNA Isolation from Plant Tissues: A Hands-on Laboratory Experimental Experience for Undergraduates

The practice of RNA isolation in undergraduate experimental courses is rare because of the existence of robust, ubiquitous and stable ribonucleases. We reported here modifications to our original protocol for RNA isolation from plant tissues, including the recovery of nucleic acids by ethanol precipitation at 0 °C for 10 min and the assessment of RNA quality by visualizing the banding profile of the separated RNAs on a standard nondenaturing agarose gel to shorten the duration of the whole procedure and simplify the operation. As a result, the modified procedure, including RNA isolation and quality control analysis could be finished in 4 hr and divided into two sessions. Because endogenous ribonucleases released upon disruption of the organelles and vacuoles were effectively and quickly inactivated, measures were taken to protect RNA integrity throughout the whole procedure so that total RNA with high purity and integrity as well as an appropriate yield could be obtained by students. The RNA isolation protocol described here was simple, efficient, flexible, and low cost. Therefore, it is an ideal approach for undergraduates to learn about RNA techniques. The pedagogical approach of the correlation of experimental work with the rationale for the whole protocol described in this report is an effective way for undergraduates to improve their learning of the techniques of RNA isolation and analysis and the theories behind them, as well as experimental design and data analysis. © 2017 by The International Union of Biochemistry and Molecular Biology, 46(3):253-261, 2018.

Short inverted repeats contribute to localized mutability in human somatic cells

Selected repetitive sequences termed short inverted repeats (SIRs) have the propensity to form secondary DNA structures called hairpins. SIRs comprise palindromic arm sequences separated by short spacer sequences that form the hairpin stem and loop respectively. Here, we show that SIRs confer an increase in localized mutability in breast cancer, which is domain-dependent with the greatest mutability observed within spacer sequences (∼1.35-fold above background). Mutability is influenced by factors that increase the likelihood of formation of hairpins such as loop lengths (of 4-5 bp) and stem lengths (of 7-15 bp). Increased mutability is an intrinsic property of SIRs as evidenced by how almost all mutational processes demonstrate a higher rate of mutagenesis of spacer sequences. We further identified 88 spacer sequences showing enrichment from 1.8- to 90-fold of local mutability distributed across 283 sites in the genome that intriguingly, can be used to inform the biological status of a tumor.

Predicting Ion Effects in an RNA Conformational Equilibrium

We develop a partial charge-based tightly bound ion (PCTBI) model for the ion effects in RNA folding. On the basis of the Monte Carlo tightly bound ion (MCTBI) approach, the model can account for ion fluctuation and correlation effects, and can predict the ion distribution around the RNA. Furthermore, unlike the previous coarse-grained RNA charge models, where negative charges are placed on the phosphates only, the current new model considers the detailed all-atom partial charge distribution on the RNA. Thus, the model not only keeps the advantage of the MCTBI model, but also has the potential to provide important detailed information unattainable by the previous MCTBI models. For example, the model predicts the reduction in ion binding upon protein binding and ion-induced conformational switches. For hepatitis C virus genomic RNA, the model predicts a Mg²⁺-induced stabilization of a kissing motif for a cis-acting regulatory element in the genomic RNA. Extensive theory-experiment comparisons support the reliability of the theoretical predictions. Therefore, the model may serve as a robust starting point for further development of an accurate method for ion effects in an RNA conformational equilibrium and RNA-cofactor interactions.

MCTBI: a web server for predicting metal ion effects in RNA structures

Metal ions play critical roles in RNA structure and function. However, web servers and software packages for predicting ion effects in RNA structures are notably scarce. Furthermore, the existing web servers and software packages mainly neglect ion correlation and fluctuation effects, which are potentially important for RNAs. We here report a new web server, the MCTBI server (http://rna.physics.missouri.edu/MCTBI), for the prediction of ion effects for RNA structures. This server is based on the recently developed MCTBI, a model that can account for ion correlation and fluctuation effects for nucleic acid structures and can provide improved predictions for the effects of metal ions, especially for multivalent ions such as Mg²⁺ effects, as shown by extensive theory-experiment test results. The MCTBI web server predicts metal ion binding fractions, the most probable bound ion distribution, the electrostatic free energy of the system, and the free energy components. The results provide mechanistic insights into the role of metal ions in RNA structure formation and folding stability, which is important for understanding RNA functions and the rational design of RNA structures.

A New Method to Predict Ion Effects in RNA Folding

The strong interaction between metal ions in solution and highly charged RNA molecules is critical for RNA structure formation and stabilization. Metal ions binding to RNA can induce RNA structural changes that are important for RNA cellular functions. Therefore, quantitative modeling of the ion effects is essential for RNA structure prediction and RNA-based molecular design. Recently, inspired by the increasing experimental evidence that supports the importance of ion correlation and fluctuation in ion-RNA interactions, we developed a new computational model, Monte Carlo Tightly Bound Ion (MCTBI) model. The validity of the model is shown by the improved accuracy in the predictions for ion binding properties and ion-dependent free energies for RNA structures. In this chapter, using homodimeric tetraloop-receptor docking as an illustrative example, we showcase the MCTBI method for the computational prediction of the ion effects in RNA folding.

DNA-templated silver nanoclusters: structural correlation and fluorescence modulation

12 years after the introduction of DNA-templated silver nanoclusters (DNA-AgNCs), exciting progress has been made and yet we are still in the midst of trying to fully understand this nanomaterial. The prominent excellence of DNA-AgNCs is undoubtedly its modulatable emission property, of which how variation in DNA templates causes emission tuning remains elusive. Based on the up-to-date DNA-AgNCs, we aim to establish the correlation between the structure/sequence of DNA templates and emission behaviour of AgNCs. Herein, we systematically present a wide-range of DNA-AgNCs based on the structural complexity of the DNA templates, including single-stranded DNA (ssDNA), double-stranded DNA (dsDNA), triple-stranded DNA (tsDNA) and DNA nanostructures. For each DNA category, we discuss the emission property, quantum yield and synthesis condition of the respective AgNCs, before cross-comparing the impact of different DNA scaffolds on the properties of AgNCs. A future outlook for this area is given as a conclusion. By putting the information together, this review may shed new light on understanding DNA-AgNCs while we are expecting continuous breakthroughs in this field.

Assessment for Melting Temperature Measurement of Nucleic Acid by HRM

High resolution melting (HRM), with a high sensitivity to distinguish the nucleic acid species with small variations, has been widely applied in the mutation scanning, methylation analysis, and genotyping. For the aim of extending HRM for the evaluation of thermal stability of nucleic acid secondary structures on sequence dependence, we investigated effects of the dye of EvaGreen, metal ions, and impurities (such as dNTPs) on melting temperature (T_m ) measurement by HRM. The accuracy of HRM was assessed as compared with UV melting method, and little difference between the two methods was found when the DNA T_m was higher than 40°C. Both insufficiency and excessiveness of EvaGreen were found to give rise to a little bit higher T_m , showing that the proportion of dye should be considered for precise T_m measurement of nucleic acids. Finally, HRM method was also successfully used to measure T_m s of DNA triplex, hairpin, and RNA duplex. In conclusion, HRM can be applied in the evaluation of thermal stability of nucleic acid (DNA or RNA) or secondary structural elements (even when dNTPs are present).

Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.

rRNA Binding Sites and the Molecular Mechanism of Action of the Tetracyclines

The tetracycline antibiotics are known to be effective in the treatment of both infectious and noninfectious disease conditions. The 16S rRNA binding mechanism currently held for the antibacterial action of the tetracyclines does not explain their activity against viruses, protozoa that lack mitochondria, and noninfectious conditions. Also, the mechanism by which the tetracyclines selectively inhibit microbial protein synthesis against host eukaryotic protein synthesis despite conservation of ribosome structure and functions is still questionable. Many studies have investigated the binding of the tetracyclines to the 16S rRNA using the small ribosomal subunit of different bacterial species, but there seems to be no agreement between various reports on the exact binding site on the 16S rRNA. The wide range of activity of the tetracyclines against a broad spectrum of bacterial pathogens, viruses, protozoa, and helminths, as well as noninfectious conditions, indicates a more generalized effect on RNA. In the light of recent evidence that the tetracyclines bind to various synthetic double-stranded RNAs (dsRNAs) of random base sequences, suggesting that the double-stranded structures may play a more important role in the binding of the tetracyclines to RNA than the specific base pairs, as earlier speculated, it is imperative to consider possible alternative binding modes or sites that could help explain the mechanisms of action of the tetracyclines against various pathogens and disease conditions.

Monte Carlo Tightly Bound Ion Model: Predicting Ion-Binding Properties of RNA with Ion Correlations and Fluctuations

Experiments have suggested that ion correlation and fluctuation effects can be potentially important for multivalent ions in RNA folding. However, most existing computational methods for the ion electrostatics in RNA folding tend to ignore these effects. The previously reported tightly bound ion (TBI) model can treat ion correlation and fluctuation but its applicability to biologically important RNAs is severely limited by the low computational efficiency. Here, on the basis of Monte Carlo sampling for the many-body ion distribution, we develop a new computational model, the Monte Carlo tightly bound ion (MCTBI) model, for ion-binding properties around an RNA. Because of an enhanced sampling algorithm for ion distribution, the model leads to a significant improvement in computational efficiency. For example, for a 160-nt RNA, the model causes a more than 10-fold increase in the computational efficiency, and the improvement in computational efficiency is more pronounced for larger systems. Furthermore, unlike the earlier model that describes ion distribution using the number of bound ions around each nucleotide, the current MCTBI model is based on the three-dimensional coordinates of the ions. The higher efficiency of the model allows us to treat the ion effects for medium to large RNA molecules, RNA-ligand complexes, and RNA-protein complexes. This new model together with proper RNA conformational sampling and the energetics model may serve as a starting point for further development for the ion effects in RNA folding and conformational changes and for large nucleic acid systems.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters.

Exploring ion induced folding of a single-stranded DNA oligomer from molecular simulation studies

One crucial issue in DNA hydration is the effect of salts on its conformational features. This has relevance in biology as cations present in the cellular environment shield the negative charges on the DNA backbone, thereby reducing the repulsive force between them. By screening the negative charges along the backbone, cations stabilize the folded structure of DNA. To study the effect of the added salt on single-stranded DNA (ss-DNA) conformations, we have performed room temperature molecular dynamics simulations of an aqueous solution containing the ss-DNA dodecamer with the 5'-CGCGAATTCGCG-3' sequence in the presence of 0.2, 0.5, and 0.8 M NaCl. Our calculations reveal that in the presence of the salt, the DNA molecule forms more collapsed coil-like conformations due to the screening of negative charges along the backbone. Additionally, we demonstrated that the formation of an octahedral inner-sphere complex by the strongly bound ion plays an important role in the stabilization of such folded conformation of DNA. Importantly, it is found that ion-DNA interactions can also explain the formation of non-sequential base stackings with longer lifetimes. Such non-sequential base stackings further stabilize the collapsed coil-like folded form of the DNA oligomer.

Salt contribution to the flexibility of single-stranded nucleic acid offinite length

Nucleic acids are negatively charged macromolecules and their structure properties are strongly coupled to metal ions in solutions. In this article, the salt effects on the flexibility of single-stranded (ss) nucleic acid chain ranging from 12 to 120 nucleotides are investigated systematically by the coarse-grained Monte Carlo simulations where the salt ions are considered explicitly and the ss chain is modeled with the virtual-bond structural model. Our calculations show that, the increase of ion concentration causes the structural collapse of ss chain and multivalent ions are much more efficient in causing such collapse, and both trivalent/small divalent ions can induce more compact state than a random relaxation state. We found that monovalent, divalent, and trivalent ions can all overcharge ss chain, and the dominating source for such overcharging changes from ion-exclusion-volume effect to ion Coulomb correlations. In addition, the predicted Na(+) and Mg(2+)-dependent persistence length l(p)'s of ss nucleic acid are in accordance with the available experimental data, and through systematic calculations, we obtained the empirical formulas for l(p) as a function of [Na(+)], [Mg(2+)] and chain length.

Biochemical properties of Bacillus Calmette Guerin ribonuclease III

Double-stranded RNA (dsRNA) is discovered to participate in the regulation of gene expression in both bacterial and eukaryotic cells. Members of ribonuclease III (RNase III) family recognize RNA motifs and cleave substrates at specific sites in a divalent-metal-ion-dependent manner. In this study, we find the RNase III from Bacillus Calmette Guerin (BCG-RNase III) cleaves small hairpin RNA based on the conserved stem structure associated with Mycobacterium 16S ribosomal RNA precursor at specific sites which are not determined. To evaluate the influence of remnant endogenous ribonucleases from expression host on RNA cleavage assays for RNase III, we use E44A and D48A mutant of the enzyme to perform RNA cleavage assays and find that remnant ribonucleases have no effect on cleavage assays. The RNA cleavage activity of the enzyme can be supported by Mg(2+), Mn(2+), and Co(2+) and enhanced with the increasing salt concentration. The catalytic activity of the enzyme is exhibited when the temperature of the reaction buffer ranges from 30 to 55 °C and the pH of the buffer from 7.0 to 10.0. Two major cleavage sites in RNA substrate are identified using RNA Ligase Mediated Rapid Amplification of cDNA Ends (RLM-RACE).

RNA folding: structure prediction, folding kinetics and ion electrostatics

Beyond the "traditional" functions such as gene storage, transport and protein synthesis, recent discoveries reveal that RNAs have important "new" biological functions including the RNA silence and gene regulation of riboswitch. Such functions of noncoding RNAs are strongly coupled to the RNA structures and proper structure change, which naturally leads to the RNA folding problem including structure prediction and folding kinetics. Due to the polyanionic nature of RNAs, RNA folding structure, stability and kinetics are strongly coupled to the ion condition of solution. The main focus of this chapter is to review the recent progress in the three major aspects in RNA folding problem: structure prediction, folding kinetics and ion electrostatics. This chapter will introduce both the recent experimental and theoretical progress, while emphasize the theoretical modelling on the three aspects in RNA folding.

Gold nanoparticle interference study during the isolation, quantification, purity and integrity analysis of RNA

Investigations have been conducted regarding the interference of nanoparticles (NPs) with different toxicological assay systems, but there is a lack of validation when conducting routine tests for nucleic acid isolation, quantification, integrity, and purity analyses. The interference of citrate-capped gold nanoparticles (AuNPs) was investigated herein. The AuNPs were added to either BEAS-2B bronchial human cells for 24 h, the isolated pure RNA, or added during the isolation procedure, and the resultant interaction was assessed. Total RNA that was isolated from untreated BEAS-2B cells was spiked with various concentrations (v/v%) of AuNPs and quantified. A decrease in the absorbance spectrum (220–340 nm) was observed in a concentration-dependent manner. The 260 and 280 nm absorbance ratios that traditionally infer RNA purity were also altered. Electrophoresis was performed to determine RNA integrity, but could not differentiate between AuNP-exposed samples. However, the spiked post-isolation samples did produce differences in spectra (190–220 nm), where shifts were observed at a shorter wavelength. These shifts could be due to alterations to chromophores found in nucleic acids. The co-isolation samples, spiked with 100 µL AuNP during the isolation procedure, displayed a peak shift to a longer wavelength and were similar to the results obtained from a 24 h AuNP treatment, under non-cytotoxic test conditions. Moreover, hyperspectral imaging using CytoViva dark field microscopy did not detect AuNP spectral signatures in the RNA isolated from treated cells. However, despite the lack of AuNPs in the final RNA product, structural changes in RNA could still be observed between 190–220 nm. Consequently, full spectral analyses should replace the traditional ratios based on readings at 230, 260, and 280 nm. These are critical points of analyses, validation, and optimization for RNA-based techniques used to assess AuNPs effects.

Molecular dynamics simulation of a single-stranded DNA with heterogeneous distribution of nucleobases in aqueous medium

The DNA metabolic processes often involve single-stranded DNA (ss-DNA) molecules as important intermediates. In the absence of base complementarity, ss-DNAs are more flexible and interact strongly with water in aqueous media. Ss-DNA-water interactions are expected to control the conformational flexibility of the DNA strand, which in turn should influence the properties of the surrounding water molecules. We have performed room temperature molecular dynamics simulation of an aqueous solution containing the ss-DNA dodecamer, 5'-CGCGAATTCGCG-3'. The conformational flexibility of the DNA strand and the microscopic structure and ordering of water molecules around it have been explored. The simulation reveals transformation of the initial base-stacked form of the ss-DNA to a fluctuating collapsed coil-like conformation with the formation of a few non-sequentially stacked base pairs. A preliminary analysis shows further collapse of the DNA conformation in presence of additional salt (NaCl) due to screening of negative charges along the backbone by excess cations. Additionally, higher packing of water molecules within a short distance from the DNA strand is found to be associated with realignment of water molecules by breaking their regular tetrahedral ordering.

Interaction between functionalized gold nanoparticles in physiological saline

The interactions between functionalized noble-metal particles in an aqueous solution are central to applications relying on controlled equilibrium association. Herein, we obtain the potentials of mean force (PMF) for pair-interactions between functionalized gold nanoparticles (AuNPs) in physiological saline. These results are based upon >1000 ns experiments in silico of all-atom model systems under equilibrium and non-equilibrium conditions. Four types of functionalization are built by coating each globular Au144 cluster with 60 thiolate groups: GS-AuNP (glutathionate), PhS-AuNP (thiophenol), CyS-AuNP (cysteinyl), and p-APhS-AuNP (para-amino-thiophenol), which are, respectively, negatively charged, hydrophobic (neutral-nonpolar), hydrophilic (neutral-polar), and positively charged at neutral pH. The results confirm the behavior expected of neutral (hydrophilic or hydrophobic) particles in a dilute aqueous environment, however the PMF curves demonstrate that the charged AuNPs interact with one another in a unique way-mediated by H2O molecules and an electrolyte (Na(+), Cl(-))-in a physiological environment. In the case of two GS-AuNPs, the excess, neutralizing Na(+) ions form a mobile (or 'dynamic') cloud of enhanced concentration between the like-charged GS-AuNPs, inducing a moderate attraction (∼25 kT) between them. Furthermore, to a lesser degree, for a pair of p-APhS-AuNPs, the excess, neutralizing Cl(-) ions (less mobile than Na(+)) also form a cloud of higher concentration between the two like-charged p-APhS-AuNPs, inducing weaker yet significant attractions (∼12 kT). On combining one GS- with one p-APhS-AuNP, the direct, attractive Coulombic force is completely screened out while the solvation effects give rise to moderate repulsion between the two unlike-charged AuNPs.

Counterion and polythymidine loop-length-dependent folding and thermodynamic stability of DNA hairpins reveal the unusual counterion-dependent stability of tetraloop hairpins

Stem-loop DNA hairpins containing a 5-base-pair (bp) stem and single-stranded polythymidine loop were investigated using thermodynamic melting analysis and stopped-flow kinetics. These studies revealed the thermodynamic stability and folding kinetics as a function of loop length and counterion concentration. Our results show the unusually high thermodynamic stability for tetraloop or 4 poly(dT) loop hairpin as compared with longer loop length hairpins. Furthermore, this exceptional stability is highly counterion-dependent. For example, in the higher counterion concentration regime of 50 mM NaCl and above, the tetraloop hairpin displays enhanced stability as compared with longer loop length hairpins. However, at lower counterion concentration of 25 mM NaCl and below, the thermal stability of tetraloop hairpin is consistent with the longer loop hairpins. The enhanced stability of tetraloop hairpins at higher counterion concentration can be explained on the basis of the combined entropic effect of loop closure as well as base stacking in the loop regions. The stability of longer loop length hairpins at all counterion concentrations as well as tetraloop hairpin at lower counterion concentration can be explained on the basis of entropic effect of loop closure alone. The thermodynamic parameters at lower and higher counterion concentrations were determined to quantify the enhanced stability of base-stacking effects occurring at higher counterion concentrations. For example, for 100 mM NaCl, excess Gibbs energy and enthalpy due to base stacking within the tetraloops were measured to be -1.2 ± 0.14 and -3.28 ± 0.32 kcal/mol, respectively, whereas, no excess of Gibbs energy and enthalpy was observed for 0, 5, 10, and 25 mM NaCl. These findings suggest significant base-stacking interactions occurring in the loop region of the tetraloop hairpins at higher counterion concentration and less significant base-stacking interactions in the lower counterion concentration regime. We suggest that at higher counterion concentrations, hydrophobic collapse of the nucleotides in the loop may be enhanced due to the increased polarity of the solvent, thereby enhancing base-stacking interactions that contribute to unusually high stability.

Coarse-grained model for predicting RNA folding thermodynamics

We present a thermodynamically robust coarse-grained model to simulate folding of RNA in monovalent salt solutions. The model includes stacking, hydrogen bond, and electrostatic interactions as fundamental components in describing the stability of RNA structures. The stacking interactions are parametrized using a set of nucleotide-specific parameters, which were calibrated against the thermodynamic measurements for single-base stacks and base-pair stacks. All hydrogen bonds are assumed to have the same strength, regardless of their context in the RNA structure. The ionic buffer is modeled implicitly, using the concept of counterion condensation and the Debye-Hückel theory. The three adjustable parameters in the model were determined by fitting the experimental data for two RNA hairpins and a pseudoknot. A single set of parameters provides good agreement with thermodynamic data for the three RNA molecules over a wide range of temperatures and salt concentrations. In the process of calibrating the model, we establish the extent of counterion condensation onto the single-stranded RNA backbone. The reduced backbone charge is independent of the ionic strength and is 60% of the RNA bare charge at 37 °C. Our model can be used to predict the folding thermodynamics for any RNA molecule in the presence of monovalent ions.

Ion-mediated RNA structural collapse: effect of spatial confinement

RNAs are negatively charged molecules that reside in cellular environments with macromolecular crowding. Macromolecular confinement can influence the ion effects in RNA folding. In this work, using the recently developed tightly bound ion model for ion fluctuation and correlation, we investigate the effect of confinement on ion-mediated RNA structural collapse for a simple model system. We find that for both Na(+) and Mg(2+), the ion efficiencies in mediating structural collapse/folding are significantly enhanced by the structural confinement. This enhancement of ion efficiency is attributed to the decreased electrostatic free-energy difference between the compact conformation ensemble and the (restricted) extended conformation ensemble due to the spatial restriction.

Predicting ion-nucleic acid interactions by energy landscape-guided sampling

The recently developed Tightly Bound Ion (TBI) model offers improved predictions for ion effect in nucleic acid systems by accounting for ion correlation and fluctuation effects. However, further application of the model to larger systems is limited by the low computational efficiency of the model. Here, we develop a new computational efficient TBI model using free energy landscape-guided sampling method. The method leads to drastic reduction in the computer time by a factor of 50 for RNAs of 50-100 nucleotides long. The improvement in the computational efficiency would be more significant for larger structures. To test the new method, we apply the model to predict the free energies and the number of bound ions for a series of RNA folding systems. The validity of this new model is supported by the nearly exact agreement with the results from the original TBI model and the agreement with the experimental data. The method may pave the way for further applications of the TBI model to treat a broad range of biologically significant systems such as tetraloop-receptor and riboswitches.

Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na(+)]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg(2+) salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.

Electrostatics of nucleic acid folding under conformational constraint

RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.