Ion-mediated nucleic acid helix-helix interactions

Salt ions are essential for the folding of nucleic acids. We use the tightly bound ion (TBI) model, which can account for the correlations and fluctuations for the ions bound to the nucleic acids, to investigate the electrostatic free-energy landscape for two parallel nucleic acid helices in the solution of added salt. The theory is based on realistic atomic structures of the helices. In monovalent salt, the helices are predicted to repel each other. For divalent salt, while the mean-field Poisson-Boltzmann theory predicts only the repulsion, the TBI theory predicts an effective attraction between the helices. The helices are predicted to be stabilized at an interhelix distance approximately 26-36 A, and the strength of the attractive force can reach -0.37 k(B)T/bp for helix length in the range of 9-12 bp. Both the stable helix-helix distance and the strength of the attraction are strongly dependent on the salt concentration and ion size. With the increase of the salt concentration, the helix-helix attraction becomes stronger and the most stable helix-helix separation distance becomes smaller. For divalent ions, at very high ion concentration, further addition of ions leads to the weakening of the attraction. Smaller ion size causes stronger helix-helix attraction and stabilizes the helices at a shorter distance. In addition, the TBI model shows that a decrease in the solvent dielectric constant would enhance the ion-mediated attraction. The theoretical findings from the TBI theory agree with the experimental measurements on the osmotic pressure of DNA array as well as the results from the computer simulations.     

Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na⁺ and Mg²⁺) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na⁺ and Mg²⁺ concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.     

RNA helix stability in mixed Na+/Mg2+ solution

A recently developed tightly bound ion model can account for the correlation and fluctuation (i.e., different binding modes) of bound ions. However, the model cannot treat mixed ion solutions, which are physiologically relevant and biologically significant, and the model was based on B-DNA helices and thus cannot directly treat RNA helices. In the present study, we investigate the effects of ion correlation and fluctuation on the thermodynamic stability of finite length RNA helices immersed in a mixed solution of monovalent and divalent ions. Experimental comparisons demonstrate that the model gives improved predictions over the Poisson-Boltzmann theory, which has been found to underestimate the roles of multivalent ions such as Mg2+ in stabilizing DNA and RNA helices. The tightly bound ion model makes quantitative predictions on how the Na+-Mg2+ competition determines helix stability and its helix length-dependence. In addition, the model gives empirical formulas for the thermodynamic parameters as functions of Na+/Mg2+ concentrations and helix length. Such formulas can be quite useful for practical applications.     

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

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

Onset of DNA aggregation in presence of monovalent and multivalent counterions

We address theoretically aggregation of DNA segments by multivalent polyamines such as spermine and spermidine. In experiments, the aggregation occurs above a certain threshold concentration of multivalent ions. We demonstrate that the dependence of this threshold on the concentration of DNA has a simple form. When the DNA concentration c(DNA) is smaller than the monovalent salt concentration, the threshold multivalent ion concentration depends linearly on c(DNA), having the form alphac(DNA) + beta. The coefficients alpha and beta are related to the density profile of multivalent counterions around isolated DNA chains, at the onset of their aggregation. This analysis agrees extremely well with recent detailed measurements on DNA aggregation in the presence of spermine. From the fit to the experimental data, the number of condensed multivalent counterions per DNA chain can be deduced. A few other conclusions can then be reached: 1), the number of condensed spermine ions at the onset of aggregation decreases with the addition of monovalent salt; 2), the Poisson-Boltzmann theory overestimates the number of condensed multivalent ions at high monovalent salt concentrations; and 3), our analysis of the data indicates that the DNA charge is not overcompensated by spermine at the onset of aggregation.     

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.     

Predicting Ion Binding Properties for RNA Tertiary Structures

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

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.     

Quantitative analysis of the ion-dependent folding stability of DNA triplexes

A DNA triplex is formed through binding of a third strand to the major groove of a duplex. Due to the high charge density of a DNA triplex, metal ions are critical for its stability. We recently developed the tightly bound ion (TBI) model for ion-nucleic acids interactions. The model accounts for the potential correlation and fluctuations of the ion distribution. We now apply the TBI model to analyze the ion dependence of the thermodynamic stability for DNA triplexes. We focus on two experimentally studied systems: a 24-base DNA triplex and a pair of interacting 14-base triplexes. Our theoretical calculations for the number of bound ions indicate that the TBI model provides improved predictions for the number of bound ions than the classical Poisson-Boltzmann (PB) equation. The improvement is more significant for a triplex, which has a higher charge density than a duplex. This is possibly due to the higher ion concentration around the triplex and hence a stronger ion correlation effect for a triplex. In addition, our analysis for the free energy landscape for a pair of 14-mer triplexes immersed in an ionic solution shows that divalent ions could induce an attractive force between the triplexes. Furthermore, we investigate how the protonated cytosines in the triplexes affect the stability of the triplex helices.     

Effect of the number of nucleic acid oligomer charges on the salt dependence of stability (DeltaG 37degrees) and melting temperature (Tm): NLPB analysis of experimental data

For nucleic acid oligomers with variable chain lengths, the salt concentration ([salt]) dependences of the denaturation temperature (T(m)) and of the free energy of helix formation at 37 degrees C (Delta) are predicted using nonlinear Poisson-Boltzmann (NLPB) calculations. Analysis of experimental data reveals that the ratio of the [salt] derivative of melting temperature (ST(m) = dT(m)/d log[salt]) to the value for a polymer with the same base composition (ST(m)/ST(m, infinity)) is independent of base composition but strongly dependent on the number of DNA charges (/Z/) below approximately 8 bp for two-strand helices (formed from association of two complementary strands) and below approximately 18 bp for hairpin helices (formed from folding of one self-complementary strand). We interpret these ST(m)/ST(m, infinity) ratios in terms of the ratio of thermodynamic ion release from the oligomer (Deltan(u), per charge) to that from the same oligomer embedded in polymeric DNA (Deltan(u, infinity), per charge). Experimental values of ST(m)/ST(m, infinity) and its dependence on /Z/ are in good agreement with NLPB predictions for a preaveraged (essential structural) model of DNA. In particular, the NLPB calculations describe the stronger /Z/ dependence of ST(m) observed for melting of oligomeric hairpin helices than for melting of two-strand helices. These calculations predict an experimentally detectable (>or=10%) difference between ST(m) and ST(m, infinity) which increases strongly with decreasing length for two-strand helix lengths of <15 bp and for hairpin helix lengths of <30 bp. From NLPB values of Deltan(u)/Deltan(u, infinity), we predict Delta as a function of [salt] and /Z/. Predictions of thermodynamic and thermal stabilities of oligomeric helices as functions of length and [salt] are consistent with and represent a significant refinement of the average oligomer salt effect currently in use in nearest neighbor stability predictions.     

Salt dependence of nucleic acid hairpin stability

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

Salt contribution to the flexibility of single‐stranded nucleic acid of finite 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²⁺‐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²⁺] and chain length. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 370–381, 2013.     

Mono- and trivalent ions around DNA: a small-angle scattering study of competition and interactions

The presence of small numbers of multivalent ions in DNA-containing solutions results in strong attractive forces between DNA strands. Despite the biological importance of this interaction, e.g., DNA condensation, its physical origin remains elusive. We carried out a series of experiments to probe interactions between short DNA strands as small numbers of trivalent ions are included in a solution containing DNA and monovalent ions. Using resonant (anomalous) and nonresonant small angle x-ray scattering, we coordinated measurements of the number and distribution of each ion species around the DNA with the onset of attractive forces between DNA strands. DNA-DNA interactions occur as the number of trivalent ions increases. Surprisingly good agreement is found between data and size-corrected numerical Poisson-Boltzmann predictions of ion competition for non- and weakly interacting DNAs. We also obtained an estimate for the minimum number of trivalent ions needed to initiate DNA-DNA attraction.     

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

Background

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

Results

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

Conclusions

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

Electrostatic free energy landscapes for nucleic acid helix assembly

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

Determining the Locations of Ions and Water around DNA from X-Ray Scattering Measurements

Nucleic acids carry a negative charge, attracting salt ions and water. Interactions with these components of the solvent drive DNA to condense, RNA to fold, and proteins to bind. To understand these biological processes, knowledge of solvent structure around the nucleic acids is critical. Yet, because they are often disordered, ions and water evade detection by x-ray crystallography and other high-resolution methods. Small-angle x-ray scattering (SAXS) is uniquely sensitive to the spatial correlations between solutes and the surrounding solvent. Thus, SAXS provides an experimental constraint to guide or test emerging solvation theories. However, the interpretation of SAXS profiles is nontrivial because of the difficulty in separating the scattering signals of each component: the macromolecule, ions, and hydration water. Here, we demonstrate methods for robustly deconvoluting these signals, facilitating a more straightforward comparison with theory. Using SAXS data collected on an absolute intensity scale for short DNA duplexes in solution with Na⁺, K⁺, Rb⁺, or Cs⁺ counterions, we mathematically decompose the scattering profiles into components (DNA, water, and ions) and validate the decomposition using anomalous scattering measurements. In addition, we generate a library of physically motivated ion atmosphere models and rank them by agreement with the scattering data. The best-fit models have relatively compact ion atmospheres when compared to predictions from the mean-field Poisson-Boltzmann theory of electrostatics. Thus, the x-ray scattering methods presented here provide a valuable measurement of the global structure of the ion atmosphere that can be used to test electrostatics theories that go beyond the mean-field approximation.