Cation-Anion Interactions within the Nucleic Acid Ion Atmosphere Revealed by Ion Counting

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.

Determination of Ion Atmosphere Effects on the Nucleic Acid Electrostatic Potential and Ligand Association Using AH+·C Wobble Formation in Double-Stranded DNA

The high charge density of nucleic acids and resulting ion atmosphere profoundly influence the conformational landscape of RNA and DNA and their association with small molecules and proteins. Electrostatic theories have been applied to quantitatively model the electrostatic potential surrounding nucleic acids and the effects of the surrounding ion atmosphere, but experimental measures of the potential and tests of these models have often been complicated by conformational changes and multisite binding equilibria, among other factors. We sought a simple system to further test the basic predictions from electrostatics theory and to measure the energetic consequences of the nucleic acid electrostatic field. We turned to a DNA system developed by Bevilacqua and co-workers that involves a proton as a ligand whose binding is accompanied by formation of an internal AH⁺·C wobble pair [Siegfried, N. A., et al. Biochemistry, 2010, 49, 3225]. Consistent with predictions from polyelectrolyte models, we observed logarithmic dependences of proton affinity versus salt concentration of −0.96 ± 0.03 and −0.52 ± 0.01 with monovalent and divalent cations, respectively, and these results help clarify prior results that appeared to conflict with these fundamental models. Strikingly, quantitation of the ion atmosphere content indicates that divalent cations are preferentially lost over monovalent cations upon A·C protonation, providing experimental indication of the preferential localization of more highly charged cations to the inner shell of the ion atmosphere. The internal AH⁺·C wobble system further allowed us to parse energetic contributions and extract estimates for the electrostatic potential at the position of protonation. The results give a potential near the DNA surface at 20 mM Mg²⁺ that is much less substantial than at 20 mM K⁺ (−120 mV vs −210 mV). These values and difference are similar to predictions from theory, and the potential is substantially reduced at higher salt, also as predicted; however, even at 1 M K⁺ the potential remains substantial, counter to common assumptions. The A·C protonation module allows extraction of new properties of the ion atmosphere and provides an electrostatic meter that will allow local electrostatic potential and energetics to be measured within nucleic acids and their complexes with proteins.

Theory and Modeling of RNA Structure and Interactions with Metal Ions and Small Molecules

In addition to continuous rapid progress in RNA structure determination, probing, and biophysical studies, the past decade has seen remarkable advances in the development of a new generation of RNA folding theories and models. In this article, we review RNA structure prediction models and models for ion-RNA and ligand-RNA interactions. These new models are becoming increasingly important for a mechanistic understanding of RNA function and quantitative design of RNA nanotechnology. We focus on new methods for physics-based, knowledge-based, and experimental data-directed modeling for RNA structures and explore the new theories for the predictions of metal ion and ligand binding sites and metal ion-dependent RNA stabilities. The integration of these new methods with theories about the cellular environment effects in RNA folding, such as molecular crowding and cotranscriptional kinetic effects, may ultimately lead to an all-encompassing RNA folding model.

Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations

Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.

Extracting water and ion distributions from solution x-ray scattering experiments

Small-angle X-ray scattering measurements can provide valuable information about the solvent environment around biomolecules, but it can be difficult to extract solvent-specific information from observed intensity profiles. Intensities are proportional to the square of scattering amplitudes, which are complex quantities. Amplitudes in the forward direction are real, and the contribution from a solute of known structure (and from the waters it excludes) can be estimated from theory; hence, the amplitude arising from the solvent environment can be computed by difference. We have found that this "square root subtraction scheme" can be extended to non-zero q values, out to 0.1 Å(-1) for the systems considered here, since the phases arising from the solute and from the water environment are nearly identical in this angle range. This allows us to extract aspects of the water and ion distributions (beyond their total numbers), by combining experimental data for the complete system with calculations for the solutes. We use this approach to test molecular dynamics and integral-equation (3D-RISM (three-dimensional reference interaction site model)) models for solvent structure around myoglobin, lysozyme, and a 25 base-pair duplex DNA. Comparisons can be made both in Fourier space and in terms of the distribution of interatomic distances in real space. Generally, computed solvent distributions arising from the MD simulations fit experimental data better than those from 3D-RISM, even though the total small-angle X-ray scattering patterns are very similar; this illustrates the potential power of this sort of analysis to guide the development of computational models.

Does Cation Size Affect Occupancy and Electrostatic Screening of the Nucleic Acid Ion Atmosphere?

Electrostatics are central to all aspects of nucleic acid behavior, including their folding, condensation, and binding to other molecules, and the energetics of these processes are profoundly influenced by the ion atmosphere that surrounds nucleic acids. Given the highly complex and dynamic nature of the ion atmosphere, understanding its properties and effects will require synergy between computational modeling and experiment. Prior computational models and experiments suggest that cation occupancy in the ion atmosphere depends on the size of the cation. However, the computational models have not been independently tested, and the experimentally observed effects were small. Here, we evaluate a computational model of ion size effects by experimentally testing a blind prediction made from that model, and we present additional experimental results that extend our understanding of the ion atmosphere. Giambasu et al. developed and implemented a three-dimensional reference interaction site (3D-RISM) model for monovalent cations surrounding DNA and RNA helices, and this model predicts that Na(+) would outcompete Cs(+) by 1.8-2.1-fold; i.e., with Cs(+) in 2-fold excess of Na(+) the ion atmosphere would contain an equal number of each cation (Nucleic Acids Res. 2015, 43, 8405). However, our ion counting experiments indicate that there is no significant preference for Na(+) over Cs(+). There is an ∼25% preferential occupancy of Li(+) over larger cations in the ion atmosphere but, counter to general expectations from existing models, no size dependence for the other alkali metal ions. Further, we followed the folding of the P4-P6 RNA and showed that differences in folding with different alkali metal ions observed at high concentration arise from cation-anion interactions and not cation size effects. Overall, our results provide a critical test of a computational prediction, fundamental information about ion atmosphere properties, and parameters that will aid in the development of next-generation nucleic acid computational models.

Quantifying the ion atmosphere of unfolded, single-stranded nucleic acids using equilibrium dialysis and single-molecule methods

To form secondary structure, nucleic acids (NAs) must overcome electrostatic strand–strand repulsion, which is moderated by the surrounding atmosphere of screening ions. The free energy of NA folding therefore depends on the interactions of this ion atmosphere with both the folded and unfolded states. We quantify such interactions using the preferential ion interaction coefficient or ion excess: the number of ions present near the NA in excess of the bulk concentration. The ion excess of the folded, double-helical state has been extensively studied; however, much less is known about the salt-dependent ion excess of the unfolded, single-stranded state. We measure this quantity using three complementary approaches: a direct approach of Donnan equilibrium dialysis read out by atomic emission spectroscopy and two indirect approaches involving either single-molecule force spectroscopy or existing thermal denaturation data. The results of these three approaches, each involving an independent experimental technique, are in good agreement. Even though the single-stranded NAs are flexible polymers that are expected to adopt random-coil configurations, we find that their ion atmosphere is quantitatively described by rod-like models that neglect large-scale conformational freedom, an effect that we explain in terms of the competition between the relevant structural and electrostatic length scales.