Critical assessment of nucleic acid electrostatics via experimental and computational investigation of an unfolded state ensemble

Impact of Divalent Cations on In-Layer Positional Order of DNA-Based Liquid Crystals: Implications for DNA Condensation

The layered liquid crystalline phases formed by DNA molecules, which include rigid and flexible segments ("gapped DNA"), enable the study of both end-to-end stacking and side-to-side (helix-to-helix) lateral interactions, forming a model system to study such interactions at physiologically relevant DNA and ion concentrations. The observed layer structure exhibits long-range interlayer and in-layer positional correlations. In particular, the in-layer order has implications for DNA condensation, as it reflects whether these normally repulsive interactions become attractive under certain ionic conditions. Using synchrotron small-angle X-ray scattering measurements, we investigate the impact of divalent Mg2+ cations (in addition to a constant 150 mM Na+) on the stability of the inter- and in-layer DNA ordering as a function of temperature between 5 and 65 °C. DNA constructs with different terminal base pairings were created to mediate the strength of the attractive end-to-end stacking interactions between the blunt ends of the gapped DNA constructs. We demonstrate that the stabilities at a fixed DNA concentration of both interlayer and in-layer order are significantly enhanced even at a few mM Mg2+ concentration. The stabilities are even higher at 30 mM Mg2+; however, a marked decrease is observed at 100 mM Mg2+, suggesting a change in the nature of side-by-side interactions within this Mg2+ concentration range. We discuss the implications of these results in terms of counterion-mediated DNA-DNA attraction and DNA condensation.

Twisting DNA by salt

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

Refining the RNA Force Field with Small-Angle X-ray Scattering of Helix-Junction-Helix RNA

The growing recognition of the functional and therapeutic roles played by RNA and the difficulties in gaining atomic-level insights by experiments are paving the way for all-atom simulations of RNA. One of the main impediments to the use of all-atom simulations is the imbalance between the energy terms of the RNA force fields. Through exhaustive sampling of an RNA helix-junction-helix (HJH) model using enhanced sampling, we critically assessed the select Amber force fields against small-angle X-ray scattering (SAXS) experiments. The tested AMBER99SB, DES-AMBER, and CUFIX force fields show deviations from measured profiles. First, we identified parameters leading to inconsistencies. Then, as a way to balance the forces governing RNA folding, we adopted strategies to refine hydrogen bonding, backbone, and base-stacking parameters. We validated the modified force field (HB-CUFIX) against SAXS data of the HJH model in different ionic strengths. Moreover, we tested a set of independent RNA systems to cross-validate the force field. Overall, HB-CUFIX demonstrates improved performance in studying thermodynamics and structural properties of realistic RNA motifs.

Direct Measurement of Interhelical DNA Repulsion and Attraction by Quantitative Cross-Linking

To better understand the forces that mediate nucleic acid compaction in biology, we developed the disulfide cross-linking approach xHEED (X-linking of Helices to measure Electrostatic Effects at Distance) to measure the distance-dependent encounter frequency of two DNA helices in solution. Using xHEED, we determined the distance that the electrostatic potential extends from DNA helices, the dependence of this distance on ionic conditions, and the magnitude of repulsion when two helices approach one another. Across all conditions tested, the potential falls to that of the bulk solution within 15 Å of the major groove surface. For separations of ∼30 Å, we measured a repulsion of 1.8 kcal/mol in low monovalent ion concentration (30 mM Na+), with higher Na+ concentrations ameliorating this repulsion, and 2 M Na+ or 100 mM Mg2+ eliminating it. Strikingly, we found full screening at very low Co3+ concentrations and net attraction at higher concentrations, without the higher-order DNA condensation that typically complicates studies of helical attraction. Our measurements define the relevant distances for electrostatic interactions of nucleic-acid helices in biology and introduce a new method to propel further understanding of how these forces impact biological processes.

Effects of ionic liquid types on the tribological, mechanical, and thermal properties of ionic liquid modified graphene/polyimide shape memory composites

Different ionic liquid modified graphene nanosheets (IG) were induced into polyimide (PI) to improve the tribological, thermal, and mechanical properties of shape memory IG/PI composites. The results demonstrated that when using 1-aminoethyl-3-methylimidazole bromide to modify graphene nanosheets (IG-1), the laser-driven shape recovery rate of IG-1/PI composites (IGPI-1) reached 73.02%, which was 49.36% higher than that of pure PI. In addition, the IGPI-1 composite materials reached the maximum shape recovery rate within 15 s. Additionally, under dry sliding, the addition of IG can significantly improve the tribological properties of composite materials. IGPI-1 exhibited the best self-lubricating properties. Compared with pure PI, the friction coefficient (0.19) and wear rate (2.62 × 10–5) mm3/Nm) were reduced by 44.1% and 24.2%, respectively, and the T10% of IGPI-1 increased by 32.2°C. The Tg of IGPI-1 reached 256.5°C, which was 8.4°C higher than that of pure PI. In addition, the tensile strength and modulus of IGPI-1 reached 82.3 MPa and 1.18 GPa, which were significantly increased by 33.6% and 29.8%, respectively, compared with pure PI. We hope that this work will be helpful for the preparation of shape memory materials with excellent tribological, thermal, and mechanical properties.

The structural ensemble of a Holliday junction determined by X-ray scattering interference

The DNA four-way (Holliday) junction is the central intermediate of genetic recombination, yet key aspects of its conformational and thermodynamic properties remain unclear. While multiple experimental approaches have been used to characterize the canonical X-shape conformers under specific ionic conditions, the complete conformational ensemble of this motif, especially at low ionic conditions, remains largely undetermined. In line with previous studies, our single-molecule Förster resonance energy transfer (smFRET) measurements of junction dynamics revealed transitions between two states under high salt conditions, but smFRET could not determine whether there are fast and unresolvable transitions between distinct conformations or a broad ensemble of related states under low and intermediate salt conditions. We therefore used an emerging technique, X-ray scattering interferometry (XSI), to directly probe the conformational ensemble of the Holliday junction across a wide range of ionic conditions. Our results demonstrated that the four-way junction adopts an out-of-plane geometry under low ionic conditions and revealed a conformational state at intermediate ionic conditions previously undetected by other methods. Our results provide critical information to build toward a full description of the conformational landscape of the Holliday junction and underscore the utility of XSI for probing conformational ensembles under a wide range of solution conditions.

CO2 -Philic Separation Membrane: Deep Eutectic Solvent Filled Graphene Oxide Nanoslits

CO₂ capture and sequestration is an energy-intensive industry to deal with the global greenhouse effect. Membrane separation is considered a cost-effective method to mitigate the emission of CO₂ . Though good separation performance and stability have been reported, supported ionic liquid membranes are still not widely applied for CO₂ separation due to the high cost. As a novel analogous solvent to ionic liquid, deep eutectic solvent retains the excellent merits of ionic liquid and is cheap with facile preparation. Herein, a highly CO₂ -philic separation membrane is constructed by nanoconfining choline chloride/ethylene glycol (ChCl/EG) deep eutectic solvent into graphene oxide nanoslits. Molecular dynamic simulation results indicate that the confinement makes a difference to the structure of the nanoconfined ChCl/EG liquid from their bulk, which remarkably facilitates CO₂ transport. By tuning the molar ratio of ChCl/EG and thickness of the membrane, the resultant membrane exhibits outstanding separation performance for CO₂ with excellent selectivity over other light gases, good long-term durability, and thermal stability. This makes it a promising membrane for selective CO₂ separation.

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.

Quantitative Studies of an RNA Duplex Electrostatics by Ion Counting

RNAs are one of the most charged polyelectrolytes in nature, and understanding their electrostatics is fundamental to their structure and biological functions. An effective way to characterize the electrostatic field generated by nucleic acids is to quantify interactions between nucleic acids and ions that surround the molecules. These ions form a loosely associated cloud referred to as an ion atmosphere. Although theoretical and computational studies can describe the ion atmosphere around RNAs, benchmarks are needed to guide the development of these approaches, and experiments to date that read out RNA-ion interactions are limited. Here, we present ion counting studies to quantify the number of ions surrounding well-defined model systems of RNA and DNA duplexes. We observe that the RNA duplex attracts more cations and expels fewer anions compared to the DNA duplex, and the RNA duplex interacts significantly stronger with the divalent cation Mg²⁺, despite their identical total charge. These experimental results suggest that the RNA duplex generates a stronger electrostatic field than DNA, as is predicted based on the structural differences between their helices. Theoretical calculations using a nonlinear Poisson-Boltzmann equation give excellent agreement with experiments for monovalent ions but underestimate Mg²⁺-DNA and Mg²⁺-RNA interactions by 20%. These studies provide needed stringent benchmarks to use against other all-atom theoretical models of RNA-ion interactions, interactions that likely must be accurately accounted for in structural, dynamic, and energetic terms to confidently model RNA structure, interactions, and function.

Ion counting demonstrates a high electrostatic field generated by the nucleosome

In eukaryotes, a first step towards the nuclear DNA compaction process is the formation of a nucleosome, which is comprised of negatively charged DNA wrapped around a positively charged histone protein octamer. Often, it is assumed that the complexation of the DNA into the nucleosome completely attenuates the DNA charge and hence the electrostatic field generated by the molecule. In contrast, theoretical and computational studies suggest that the nucleosome retains a strong, negative electrostatic field. Despite their fundamental implications for chromatin organization and function, these opposing views of nucleosome electrostatics have not been experimentally tested. Herein, we directly measure nucleosome electrostatics and find that while nucleosome formation reduces the complex charge by half, the nucleosome nevertheless maintains a strong negative electrostatic field. Our studies highlight the importance of considering the polyelectrolyte nature of the nucleosome and its impact on processes ranging from factor binding to DNA compaction.

Multiplexed Nucleic Acid Sensing with Single-Molecule FRET

Multiplex detection of biomolecules is important in bionanotechnology and clinical diagnostics. Multiplexing using engineered solutions such as microarrays, synthetic nanopores, and DNA barcodes is promising, but they require sophisticated design/engineering and typically yield semiquantitative information. Single-molecule fluorescence resonance energy transfer (smFRET) is an attractive tool in this regard as it enables both sensitive and quantitative detection. However, multiplexing with smFRET remains a great challenge as it requires either multiple excitation sources, an antenna system created by multiple FRET pairs, or multiple acceptors of the donor fluorophore, which complicates not only the labeling schemes but also data analysis, due to overlapping of FRET efficiencies ( E_FRET). Here, we address these currently outstanding issues by designing interconvertible hairpin-based sensors (iHabSs) with nonoverlapping E_FRET utilizing a single donor/acceptor pair and demonstrate a high-confidence multiplex detection of unlabeled nucleic acid sequences. We validated the reliability of our approach by systematically omitting one target at a time. Further, we demonstrate that these iHabSs are fully recyclable, sensitive with a limit of detection of ∼200 pM, and able to discriminate against single base mismatches. The multiplexed approach developed here has the potential to benefit the fields of biosensing and diagnostics by allowing simultaneous and quantitative detection of unlabeled nucleic acid biomarkers.

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.

Conformations of an RNA Helix-Junction-Helix Construct Revealed by SAXS Refinement of MD Simulations

RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA's recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.

How the Conformations of an Internal Junction Contribute to Fold an RNA Domain

Like proteins, some RNAs fold to compact structures. We can model functional RNAs as a series of short, rigid, base-paired elements, connected by non-base-paired nucleotides that serve as junctions. These connecting regions bend and twist, facilitating the formation of tertiary contacts that stabilize compact states. Here, we explore the roles of salt and junction sequence in determining the structures of a ubiquitous connector: an asymmetric internal loop. We focus on the J5/5a junction from the widely studied P4-P6 domain of the Tetrahymena ribozyme. Following the addition of magnesium ions to fold P4-P6, this junction bends dramatically, bringing the two halves of the RNA domain together for tertiary contact engagement. Using single-molecule fluorescence resonance energy transfer (smFRET), we examine the role of sequence and salt on model RNA constructs that contain these junction regions. We explore the wild-type J5/5a junction as well as two sequence variants. These junctions display distinct, salt-dependent conformations. Small-angle X-ray scattering (SAXS) measurements verify that these effects persist in the full-length P4-P6 domain. These measurements underscore the importance of junction sequence and interactions with ions in facilitating RNA folding.

Effects of Preferential Counterion Interactions on the Specificity of RNA Folding

The real-time search for native RNA structure is essential for the operation of regulatory RNAs. We previously reported that a fraction of the Azoarcus ribozyme achieves a compact structure in less than a millisecond. To scrutinize the forces that drive initial folding steps, we used time-resolved SAXS to compare the folding dynamics of this ribozyme in thermodynamically isostable concentrations of different counterions. The results show that the size of the fast-folding population increases with the number of available counterions and correlates with the flexibility of initial RNA structures. Within 1 ms of folding, Mg²⁺ exhibits a smaller preferential interaction coefficient per charge, ΔΓ₊/ Z, than Na⁺ or [Co(NH₃)₆]³⁺. The lower ΔΓ₊/ Z corresponds to a smaller yield of folded RNA, although Mg²⁺ stabilizes native RNA more efficiently than other ions at equilibrium. These results suggest that strong Mg²⁺-RNA interactions impede the search for globally native structure during early folding stages.

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.

Tuning RNA Flexibility with Helix Length and Junction Sequence

The increasing awareness of RNA's central role in biology calls for a new understanding of how RNAs, like proteins, recognize biological partners. Because RNA is inherently flexible, it assumes a variety of conformations. This conformational flexibility can be a critical aspect of how RNA attracts and binds molecular partners. Structurally, RNA consists of rigid basepaired duplexes, separated by flexible non-basepaired regions. Here, using an RNA system consisting of two short helices, connected by a single-stranded (non-basepaired) junction, we explore the role of helix length and junction sequence in determining the range of conformations available to a model RNA. Single-molecule Förster resonance energy transfer reports on the RNA conformation as a function of either mono- or divalent ion concentration. Electrostatic repulsion between helices dominates at low salt concentration, whereas junction sequence effects determine the conformations at high salt concentration. Near physiological salt concentrations, RNA conformation is sensitive to both helix length and junction sequence, suggesting a means for sensitively tuning RNA conformations.

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.

Entropic stabilization of folded RNA in crowded solutions measured by SAXS

Non-coding RNAs must fold into specific structures that are stabilized by metal ions and other co-solutes in the cell's interior. Large crowder molecules such as PEG stabilize a bacterial group I ribozyme so that the RNA folds in low Mg²⁺ concentrations typical of the cell's interior. To understand the thermodynamic origins of stabilization by crowder molecules, small angle X-ray scattering was used to measure the folding and helix assembly of a bacterial group I ribozyme at different temperatures and in different MgCl₂ and polyethylene glycol (PEG) concentrations. The resulting phase diagrams show that perturbations to folding by each variable do not overlap. A favorable enthalpy change drives the formation of compact, native-like structures, but requires Mg²⁺ ions at all temperatures studied (5–55°C). PEG reduces the entropic cost of helix assembly and increases correlations between RNA segments at all temperatures. The phase diagrams also revealed a semi-compact intermediate between the unfolded and folded ensemble that is locally more flexible than the unfolded state, as judged by SHAPE modification. These results suggest that environmental variables such as temperature and solute density will favor different types of RNA structures.

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.

The solution structural ensembles of RNA kink-turn motifs and their protein complexes

With the growing number of crystal structures of RNA and RNA-protein complexes, a critical next step is understanding the dynamic solution behavior of these entities in terms of conformational ensembles and energy landscapes. To this end, we have used X-ray scattering interferometry (XSI) to probe the ubiquitous RNA kink-turn motif and its complexes with the canonical kink-turn binding protein L7Ae. XSI revealed that the folded kink-turn is best described as a restricted conformational ensemble. The ions present in solution alter the nature of this ensemble, and protein binding can perturb the kink-turn ensemble without collapsing it to a unique state. This study demonstrates how XSI can reveal structural and ensemble properties of RNAs and RNA-protein complexes and uncovers the behavior of an important RNA-protein motif. This type of information will be necessary to understand, predict and engineer the behavior and function of RNAs and their protein complexes.

Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET

Structured RNA molecules play roles in central biological processes and understanding the basic forces and features that govern RNA folding kinetics and thermodynamics can help elucidate principles that underlie biological function. Here we investigate one such feature, the specific interaction of monovalent cations with a structured RNA, the P4-P6 domain of the Tetrahymena ribozyme. We employ single molecule FRET (smFRET) approaches as these allow determination of folding equilibrium and rate constants over a wide range of stabilities and thus allow direct comparisons without the need for extrapolation. These experiments provide additional evidence for specific binding of monovalent cations, Na+ and K+, to the RNA tetraloop-tetraloop receptor (TL-TLR) tertiary motif. These ions facilitate both folding and unfolding, consistent with an ability to help order the TLR for binding and further stabilize the tertiary contact subsequent to attainment of the folding transition state.

Mimicking Ribosomal Unfolding of RNA Pseudoknot in a Protein Channel

Pseudoknots are a fundamental RNA tertiary structure with important roles in regulation of mRNA translation. Molecular force spectroscopic approaches such as optical tweezers can track the pseudoknot's unfolding intermediate states by pulling the RNA chain from both ends, but the kinetic unfolding pathway induced by this method may be different from that in vivo, which occurs during translation and proceeds from the 5' to 3' end. Here we developed a ribosome-mimicking, nanopore pulling assay for dissecting the vectorial unfolding mechanism of pseudoknots. The pseudoknot unfolding pathway in the nanopore, either from the 5' to 3' end or in the reverse direction, can be controlled by a DNA leader that is attached to the pseudoknot at the 5' or 3' ends. The different nanopore conductance between DNA and RNA translocation serves as a marker for the position and structure of the unfolding RNA in the pore. With this design, we provided evidence that the pseudoknot unfolding is a two-step, multistate, metal ion-regulated process depending on the pulling direction. Most notably, unfolding in both directions is rate-limited by the unzipping of the first helix domain (first step), which is Helix-1 in the 5' → 3' direction and Helix-2 in the 3' → 5' direction, suggesting that the initial unfolding step in either pulling direction needs to overcome an energy barrier contributed by the noncanonical triplex base-pairs and coaxial stacking interactions for the tertiary structure stabilization. These findings provide new insights into RNA vectorial unfolding mechanisms, which play an important role in biological functions including frameshifting.

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

The ion atmosphere is a critical structural, dynamic, and energetic component of nucleic acids that profoundly affects their interactions with proteins and ligands. Experimental methods that "count" the number of ions thermodynamically associated with the ion atmosphere allow dissection of energetic properties of the ion atmosphere, and thus provide direct comparison to theoretical results. Previous experiments have focused primarily on the cations that are attracted to nucleic acid polyanions, but have also showed that anions are excluded from the ion atmosphere. Herein, we have systematically explored the properties of anion exclusion, testing the zeroth-order model that anions of different identity are equally excluded due to electrostatic repulsion. Using a series of monovalent salts, we find, surprisingly, that the extent of anion exclusion and cation inclusion significantly depends on salt identity. The differences are prominent at higher concentrations and mirror trends in mean activity coefficients of the electrolyte solutions. Salts with lower activity coefficients exhibit greater accumulation of both cations and anions within the ion atmosphere, strongly suggesting that cation-anion correlation effects are present in the ion atmosphere and need to be accounted for to understand electrostatic interactions of nucleic acids. To test whether the effects of cation-anion correlations extend to nucleic acid kinetics and thermodynamics, we followed the folding of P4-P6, a domain of the Tetrahymena group I ribozyme, via single-molecule fluorescence resonance energy transfer in solutions with different salts. Solutions of identical concentration but lower activity gave slower and less favorable folding. Our results reveal hitherto unknown properties of the ion atmosphere and suggest possible roles of oriented ion pairs or anion-bridged cations in the ion atmosphere for electrolyte solutions of salts with reduced activity. Consideration of these new results leads to a reevaluation of the strengths and limitations of Poisson-Boltzmann theory and highlights the need for next-generation atomic-level models of the ion atmosphere.

Multivalent ion-mediated nucleic acid helix-helix interactions: RNA versus DNA

Ion-mediated interaction is critical to the structure and stability of nucleic acids. Recent experiments suggest that the multivalent ion-induced aggregation of double-stranded (ds) RNAs and DNAs may strongly depend on the topological nature of helices, while there is still lack of an understanding on the relevant ion-mediated interactions at atomistic level. In this work, we have directly calculated the potentials of mean force (PMF) between two dsRNAs and between two dsDNAs in Co(NH3)6 (3+) (Co-Hex) solutions by the atomistic molecular dynamics simulations. Our calculations show that at low [Co-Hex], the PMFs between B-DNAs and between A-RNAs are both (strongly) repulsive. However, at high [Co-Hex], the PMF between B-DNAs is strongly attractive, while those between A-RNAs and between A-DNAs are still (weakly) repulsive. The microscopic analyses show that for A-form helices, Co-Hex would become 'internal binding' into the deep major groove and consequently cannot form the evident ion-bridge between adjacent helices, while for B-form helices without deep grooves, Co-Hex would exhibit 'external binding' to strongly bridge adjacent helices. In addition, our further calculations show that, the PMF between A-RNAs could become strongly attractive either at very high [Co-Hex] or when the bottom of deep major groove is fixed with a layer of water.

From static to dynamic: the need for structural ensembles and a predictive model of RNA folding and function

To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior. An incisive approach to understanding the function and folding of biological RNA systems involves characterizing small, simple components that are largely responsible for the behavior of complex systems including helix-junction-helix elements and tertiary motifs. State-of-the-art methods have permitted unprecedented insight into the conformational ensembles of these elements revealing, for example, that conformations of helix-junction-helix elements are confined to a small region of the ensemble, that this region is highly dependent on the junction's topology, and that the correct alignment of tertiary motifs may be a rare conformation on the overall folding landscape. Further characterization of RNA components and continued development of experimental and computational methods with the goal of quantitatively predicting RNA folding and functional behavior will be critical to understanding biological RNA systems.

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.

Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamics simulation and a Poisson-Boltzmann approach

The ion atmosphere created by monovalent (Na(+) ) or divalent (Mg(2+) ) cations surrounding a B-form DNA duplex were examined using atomistic molecular dynamics (MD) simulations and the nonlinear Poisson-Boltzmann (PB) equation. The ion distributions predicted by the two methods were compared using plots of radial and two-dimensional cation concentrations and by calculating the total number of cations and net solution charge surrounding the DNA. Na(+) ion distributions near the DNA were more diffuse in PB calculations than in corresponding MD simulations, with PB calculations predicting lower concentrations near DNA groove sites and phosphate groups and a higher concentration in the region between these locations. Other than this difference, the Na(+) distributions generated by the two methods largely agreed, as both predicted similar locations of high Na(+) concentration and nearly identical values of the number of cations and the net solution charge at all distances from the DNA. In contrast, there was greater disagreement between the two methods for Mg(2+) cation concentration profiles, as both the locations and magnitudes of peaks in Mg(2+) concentration were different. Despite experimental and simulation observations that Mg(2+) typically maintains its first solvation shell when interacting with nucleic acids, modeling Mg(2+) as an unsolvated ion during PB calculations improved the agreement of the Mg(2+) ion atmosphere predicted by the two methods and allowed for values of the number of bound ions and net solution charge surrounding the DNA from PB calculations that approached the values observed in MD simulations.

Ion competition in condensed DNA arrays in the attractive regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt(3+) Hexammine, or Co(3+)Hex) and one nonDNA-condensing ion (Mg(2+)) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co(3+)Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg(2+), the Co(3+)Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg(2+) and Co(3+)Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co(3+)Hex. It is also observed that only ∼20% DNA charge neutralization by Co(3+)Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.

Understanding nucleic acid-ion interactions

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

Coarse grained models reveal essential contributions of topological constraints to the conformational free energy of RNA bulges

Recent studies have shown that simple stereochemical constraints encoded at the RNA secondary structure level significantly restrict the orientation of RNA helices across two-way junctions and yield physically reasonable distributions of RNA 3D conformations. Here we develop a new coarse-grain model, TOPRNA, that is optimized for exploring detailed aspects of these topological constraints in complex RNA systems. Unlike prior models, TOPRNA effectively treats RNAs as collections of semirigid helices linked by freely rotatable single strands, allowing us to isolate the effects of secondary structure connectivity and sterics on 3D structure. Simulations of bulge junctions show that TOPRNA captures new aspects of topological constraints, including variations arising from deviations in local A-form structure, translational displacements of the helices, and stereochemical constraints imposed by bulge-linker nucleotides. Notably, these aspects of topological constraints define free energy landscapes that coincide with the distribution of bulge conformations in the PDB. Our simulations also quantitatively reproduce NMR RDC measurements made on HIV-1 TAR at low salt concentrations, although not for different TAR mutants or at high salt concentrations. Our results confirm that topological constraints are an important determinant of bulge conformation and dynamics and demonstrate the utility of TOPRNA for studying the topological constraints of complex RNAs.

Advances in the determination of nucleic acid conformational ensembles

Conformational changes in nucleic acids play a key role in the way genetic information is stored, transferred, and processed in living cells. Here, we describe new approaches that employ a broad range of experimental data, including NMR-derived chemical shifts and residual dipolar couplings, small-angle X-ray scattering, and computational approaches such as molecular dynamics simulations to determine ensembles of DNA and RNA at atomic resolution. We review the complementary information that can be obtained from diverse sets of data and the various methods that have been developed to combine these data with computational methods to construct ensembles and assess their uncertainty. We conclude by surveying RNA and DNA ensembles determined using these methods, highlighting the unique physical and functional insights obtained so far.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Analysis of diffuse K+ and Mg2+ ion binding to a two-base-pair kissing complex by single-molecule mechanical unfolding

The folding and stability of RNA tertiary interactions depend critically on cationic conditions. It is usually difficult, however, to isolate such effects on tertiary interactions from those on the entire RNA. By manipulating conformations of single RNA molecules using optical tweezers, we distinguished individual steps of breaking and forming of a two-base-pair kissing interaction from those of secondary folding. The binding of metal ions to the small tertiary structure appeared to be saturable with an apparent Kd of 160 mM for K(+) and 1.5 mM for Mg(2+). The kissing formation was estimated to be associated with binding of ~2-3 diffuse K(+) or Mg(2+) ions. At their saturated binding, Mg(2+) provided ~3 kcal/mol more stabilizing energy to the structure than K(+). Furthermore, the cations change the unkissing forces significantly more than the kissing ones. For example, the presence of Mg(2+) ions increased the average unkissing force from 21 pN to 44 pN, surprisingly high for breaking merely two base pairs; in contrast, the mean kissing force was changed by only 4.5 pN. Interestingly, the differential salt effects on the transition forces were not caused by different changes in the height of the kinetic barriers but were instead attributed to how different molecular structures respond to the applied force. Our results showed the importance of diffuse cation binding to the stability of tertiary interaction and demonstrated the utility of mechanical unfolding in studying tertiary interactions.

Crowders perturb the entropy of RNA energy landscapes to favor folding

Biological macromolecules have evolved to fold and operate in the crowded environment of the cell. We have shown previously that molecular crowding stabilizes folded RNA structures. Here we report SAXS measurements on a 64 kDa bacterial group I ribozyme in the presence of mono- and divalent ions and PEG crowders of different molecular weight. These experiments show that crowders always stabilize the folded RNA, but this stabilization is weaker in NaCl solutions than MgCl2 solutions. Additionally, we find that RNAs with the same global structure, parametrized by Rg, have different scattering functions depending upon the ratio of electrostatic and entropic stabilization by ions and crowders, respectively. We quantify this difference using the scattering length per scattering volume and find that this ratio is larger for RNAs that fold in lower ionic strength solutions due to the higher crowder content. We conclude that lower RNA flexibility, or reduced configurational entropy, widens the free energy gap between the unfolded and folded RNA in crowded MgCl2 solutions.

Quantifying Coulombic and solvent polarization-mediated forces between DNA helices

One of the fundamental problems in nucleic acids biophysics is to predict the different forces that stabilize nucleic acid tertiary folds. Here we provide a quantitative estimation and analysis for the forces between DNA helices in an ionic solution. Using the generalized Born model and the improved atomistic tightly binding ions model, we evaluate ion correlation and solvent polarization effects in interhelix interactions. The results suggest that hydration, Coulomb correlation and ion entropy act together to cause the repulsion and attraction between nucleic acid helices in Mg(2+) and Mn(2+) solutions, respectively. The theoretical predictions are consistent with experimental findings. Detailed analysis further suggests that solvent polarization and ion correlation both are crucial for the interhelix interactions. The theory presented here may provide a useful framework for systematic and quantitative predictions of the forces in nucleic acids folding.

DNA translocation by human uracil DNA glycosylase: role of DNA phosphate charge

Human DNA repair glycosylases must encounter and inspect each DNA base in the genome to discover damaged bases that may be present at a density of <1 in 10 million normal base pairs. This remarkable example of specific molecular recognition requires a reduced dimensionality search process (facilitated diffusion) that involves both hopping and sliding along the DNA chain. Despite the widely accepted importance of facilitated diffusion in protein-DNA interactions, the molecular features of DNA that influence hopping and sliding are poorly understood. Here we explore the role of the charged DNA phosphate backbone in sliding and hopping by human uracil DNA glycosylase (hUNG), which is an exemplar that efficiently locates rare uracil bases in both double-stranded DNA and single-stranded DNA. Substitution of neutral methylphosphonate groups for anionic DNA phosphate groups weakened nonspecific DNA binding affinity by 0.4-0.5 kcal/mol per substitution. In contrast, sliding of hUNG between uracil sites embedded in duplex and single-stranded DNA substrates persisted unabated when multiple methylphosphonate linkages were inserted between the sites. Thus, a continuous phosphodiester backbone negative charge is not essential for sliding over nonspecific DNA binding sites. We consider several alternative mechanisms for these results. A model consistent with previous structural and nuclear magnetic resonance dynamic results invokes the presence of open and closed conformational states of hUNG. The open state is short-lived and has weak or nonexistent interactions with the DNA backbone that are conducive for sliding, and the populated closed state has stronger interactions with the phosphate backbone. These data suggest that the fleeting sliding form of hUNG is a distinct weakly interacting state that facilitates rapid movement along the DNA chain and resembles the transition state for DNA dissociation.

Role of ion valence in the submillisecond collapse and folding of a small RNA domain

Following the addition of ions to trigger folding, RNA molecules undergo a transition from rigid, extended states to a compact ensemble. Determining the time scale for this collapse provides important insights into electrostatic contributions to RNA folding; however, it can be challenging to isolate the effects of purely nonspecific collapse, e.g., relaxation due to backbone charge compensation, from the concurrent formation of some tertiary contacts. To solve this problem, we decoupled nonspecific collapse from tertiary folding using a single-point mutation to eliminate tertiary contacts in the small RNA subdomain known as tP5abc. Microfluidic mixing with microsecond time resolution and Förster resonance energy transfer detection provides insight into the ionic strength-dependent transition from extended to compact ensembles. Differences in reaction rates are detected when folding is initiated by monovalent or divalent ions, consistent with equilibrium measurements illustrating the enhanced screening of divalent ions relative to monovalent ions at the same ionic strength. Ion-driven collapse is fast, and a comparison of the collapse time of the wild-type and mutant tP5abc suggests that site binding of Mg(2+) occurs on submillisecond time scales.

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.

The role of counterion valence and size in GAAA tetraloop-receptor docking/undocking kinetics

For RNA to fold into compact, ordered structures, it must overcome electrostatic repulsion between negatively charged phosphate groups by counterion recruitment. A physical understanding of the counterion-assisted folding process requires addressing how cations kinetically and thermodynamically control the folding equilibrium for each tertiary interaction in a full-length RNA. In this work, single-molecule FRET (fluorescence resonance energy transfer) techniques are exploited to isolate and explore the cation-concentration-dependent kinetics for formation of a ubiquitous RNA tertiary interaction, that is, the docking/undocking of a GAAA tetraloop with its 11-nt receptor. Rate constants for docking (k(dock)) and undocking (k(undock)) are obtained as a function of cation concentration, size, and valence, specifically for the series Na(+), K(+), Mg(2+), Ca(2+), Co(NH(3))(6)(3+), and spermidine(3+). Increasing cation concentration acceleratesk(dock)dramatically but achieves only a slight decrease in k(undock). These results can be kinetically modeled using parallel cation-dependent and cation-independent docking pathways, which allows for isolation of the folding kinetics from the interaction energetics of the cations with the undocked and docked states, respectively. This analysis reveals a preferential interaction of the cations with the transition state and docked state as compared to the undocked RNA, with the ion-RNA interaction strength growing with cation valence. However, the corresponding number of cations that are taken up by the RNA upon folding decreases with charge density of the cation. The only exception to these behaviors is spermidine(3+), whose weaker influence on the docking equilibria with respect to Co(NH(3))(6)(3+) can be ascribed to steric effects preventing complete neutralization of the RNA phosphate groups.

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.

RNA and its ionic cloud: solution scattering experiments and atomically detailed simulations

RNA molecules play critical roles in many cellular processes. Traditionally viewed as genetic messengers, RNA molecules were recently discovered to have diverse functions related to gene regulation and expression. RNA also has great potential as a therapeutic and a tool for further investigation of gene regulation. Metal ions are an integral part of RNA structure and should be considered in any experimental or theoretical study of RNA. Here, we report a multidisciplinary approach that combines anomalous small-angle x-ray scattering and molecular-dynamics (MD) simulations with explicit solvent and ions around RNA. From experiment and simulation results, we find excellent agreement in the number and distribution of excess monovalent and divalent ions around a short RNA duplex. Although similar agreement can be obtained from a continuum description of the solvent and mobile ions (by solving the Poisson-Boltzmann equation and accounting for finite ion size), the use of MD is easily extended to flexible RNA systems with thermal fluctuations. Therefore, we also model a short RNA pseudoknot and find good agreement between the MD results and the experimentally derived solution structures. Surprisingly, both deviate from crystal structure predictions. These favorable comparisons of experiment and simulations encourage work on RNA in all-atom dynamic models.

The ionic atmosphere around A-RNA: Poisson-Boltzmann and molecular dynamics simulations

The distributions of different cations around A-RNA are computed by Poisson-Boltzmann (PB) equation and replica exchange molecular dynamics (MD). Both the nonlinear PB and size-modified PB theories are considered. The number of ions bound to A-RNA, which can be measured experimentally, is well reproduced in all methods. On the other hand, the radial ion distribution profiles show differences between MD and PB. We showed that PB results are sensitive to ion size and functional form of the solvent dielectric region but not the solvent dielectric boundary definition. Size-modified PB agrees with replica exchange molecular dynamics much better than nonlinear PB when the ion sizes are chosen from atomistic simulations. The distribution of ions 14 Å away from the RNA central axis are reasonably well reproduced by size-modified PB for all ion types with a uniform solvent dielectric model and a sharp dielectric boundary between solvent and RNA. However, this model does not agree with MD for shorter distances from the A-RNA. A distance-dependent solvent dielectric function proposed by another research group improves the agreement for sodium and strontium ions, even for shorter distances from the A-RNA. However, Mg(2+) distributions are still at significant variances for shorter distances.

Entropic origin of Mg2+-facilitated RNA folding

Mg(2+) is essential for the proper folding and function of RNA, though the effect of Mg(2+) concentration on the free energy, enthalpy, and entropy landscapes of RNA folding is unknown. This work exploits temperature-controlled single-molecule FRET methods to address the thermodynamics of RNA folding pathways by probing the intramolecular docking/undocking kinetics of the ubiquitous GAAA tetraloop-receptor tertiary interaction as a function of [Mg(2+)]. These measurements yield the barrier and standard state enthalpies, entropies, and free energies for an RNA tertiary transition, in particular, revealing the thermodynamic origin of [Mg(2+)]-facilitated folding. Surprisingly, these studies reveal that increasing [Mg(2+)] promotes tetraloop-receptor interaction by reducing the entropic barrier (-TΔS(++)(dock)) and the overall entropic penalty (-TΔS(+) (dock)) for docking, with essentially negligible effects on both the activation enthalpy (ΔH(++)(dock)) and overall exothermicity (ΔH(+)(dock)). These observations contrast with the conventional notion that increasing [Mg(2+)] facilitates folding by minimizing electrostatic repulsion of opposing RNA helices, which would incorrectly predict a decrease in ΔH(++)(dock)) and ΔH(+)(dock)) with [Mg(2+)]. Instead we propose that higher [Mg(2+)] can aid RNA folding by decreasing the entropic penalty of counterion uptake and by reducing disorder of the unfolded conformational ensemble.

Ionic strength-dependent persistence lengths of single-stranded RNA and DNA

Dynamic RNA molecules carry out essential processes in the cell including translation and splicing. Base-pair interactions stabilize RNA into relatively rigid structures, while flexible non-base-paired regions allow RNA to undergo conformational changes required for function. To advance our understanding of RNA folding and dynamics it is critical to know the flexibility of these un-base-paired regions and how it depends on counterions. Yet, information about nucleic acid polymer properties is mainly derived from studies of ssDNA. Here we measure the persistence lengths (l(p)) of ssRNA. We observe valence and ionic strength-dependent differences in l(p) in a direct comparison between 40-mers of deoxythymidylate (dT(40)) and uridylate (rU(40)) measured using the powerful combination of SAXS and smFRET. We also show that nucleic acid flexibility is influenced by local environment (an adjoining double helix). Our results illustrate the complex interplay between conformation and ion environment that modulates nucleic acid function in vivo.

Crystalline Gibbs monolayers of DNA-capped nanoparticles at the air-liquid interface

Using grazing-incidence small-angle X-ray scattering in a special configuration (parallel SAXS, or parSAXS), we mapped the crystallization of DNA-capped nanoparticles across a sessile droplet, revealing the formation of crystalline Gibbs monolayers of DNA-capped nanoparticles at the air-liquid interface. We showed that the spatial crystallization can be regulated by adjusting both ionic strength and DNA sequence length and that a modified form of the Daoud-Cotton model could describe and predict the resulting changes in interparticle spacing. Gibbs monolayers at the air-liquid interface provide an ideal platform for the formation and study of equilibrium nanostructures and may afford exciting routes toward the design of programmable 2D plasmonic materials and metamaterials.