Measurement of the effect of monovalent cations on RNA hairpin stability

3' terminal nucleotides determine thermodynamic stabilities of mismatches at the ends of RNA helices

The thermodynamic stabilities of consecutive mismatches at the ends of RNA helices are determined by the 3' terminal nucleotides. More than 40 RNA duplexes containing terminal motifs of 3 or more nucleotides were studied by optical melting experiments. Up to three noncanonical pairs of nucleotides at the end of RNA helices provide additional thermodynamic stability. 3' nucleotides contribute more stability than 5' nucleotides, and purines contribute more stability than pyrimidines. The additional stability of a second or third 3' nucleotide stacking on a purine is the same for both dangling ends and consecutive terminal mismatches. Current predictions underestimate RNA duplex stabilities with terminal motifs by 1.4 kcal/mol on average, which is an order of magnitude in a binding constant at 37 degrees C. Accurate thermodynamic parameters for these terminal motifs will contribute to improvements in RNA secondary structure predictions, identification of microRNA targets, and design of siRNA therapeutics with fewer off-target effects.

Predicting electrostatic forces in RNA folding

Metal ion-mediated electrostatic interactions are critical to RNA folding. Although considerable progress has been made in mechanistic studies, the problem of accurate predictions for the ion effects in RNA folding remains unsolved, mainly due to the complexity of several potentially important issues such as ion correlation and dehydration effects. In this chapter, after giving a brief overview of the experimental findings and theoretical approaches, we focus on a recently developed new model, the tightly bound ion (TBI) model, for ion electrostatics in RNA folding. The model is unique because it can treat ion correlation and fluctuation effects for realistic RNA 3D structures. For monovalent ion (such as Na(+)) solutions, where ion correlation is weak, TBI and the Poisson-Boltzmann (PB) theory give the same results and the results agree with the experimental data. For multivalent ion (such as Mg(2+)) solutions, where ion correlation can be strong, however, TBI gives much improved predictions than the PB. Moreover, the model suggests an ion correlation-induced mechanism for the unusual efficiency of Mg(2+) ions in the stabilization of RNA tertiary folds. In this chapter, after introducing the theoretical framework of the TBI model, we will describe how to apply the model to predict ion-binding properties and ion-dependent folding stabilities.

Folding and unfolding single RNA molecules under tension

Single-molecule force spectroscopy constitutes a powerful method for probing RNA folding: It allows the kinetic, energetic, and structural properties of intermediate and transition states to be determined quantitatively, yielding new insights into folding pathways and energy landscapes. Recent advances in experimental and theoretical methods, including fluctuation theorems, kinetic theories, novel force clamps, and ultrastable instruments, have opened new avenues for study. These tools have been used to probe folding in simple model systems, for example, RNA and DNA hairpins. Knowledge gained from such systems is helping to build our understanding of more complex RNA structures composed of multiple elements, as well as how nucleic acids interact with proteins involved in key cellular activities, such as transcription and translation.

Gold nanoparticle based NSET for monitoring Mg2+ dependent RNA folding

Understanding the mechanism of how RNA molecules fold into their native structures are vital to their functional properties. Here we report for the first time that gold nanoparticle based NSET can be used for probing the transition states of an RNA unfolding reaction. Our result shows that time-dependent NSET can clearly distinguish structural transitions between unfolded to folded states. Our experimental observation point out that NSET can be used for the design of an optical based molecular ruler to track RNA folding transition states at distances more than double the distances achievable using traditional dipole-dipole Coulombic energy transfer based methods.

Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations

Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.

Unwinding RNA's secrets: advances in the biology, physics, and modeling of complex RNAs

The rapid development of our understanding of the diverse biological roles fulfilled by non-coding RNA has motivated interest in the basic macromolecular behavior, structure, and function of RNA. We focus on two areas in the behavior of complex RNAs. First, we present advances in the understanding of how RNA folding is accomplished in vivo by presenting a mechanism for the action of DEAD-box proteins. Members of this family are intimately associated with almost all cellular processes involving RNA, mediating RNA structural rearrangements and chaperoning their folding. Next, we focus on advances in understanding, and characterizing the basic biophysical forces that govern the folding of complex RNAs. Ultimately we expect that a confluence and synergy between these approaches will lead to profound understanding of RNA and its biology.

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