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

Ribozymes, which carry out phosphoryl-transfer reactions, often require Mg(2+) ions for catalytic activity. The correct folding of the active site and ribozyme tertiary structure is also regulated by metal ions in a manner that is not fully understood. Here we employ coarse-grained molecular simulations to show that individual structural elements of the group I ribozyme from the bacterium Azoarcus form spontaneously in the unfolded ribozyme even at very low Mg(2+) concentrations, and are transiently stabilized by the coordination of Mg(2+) ions to specific nucleotides. However, competition for scarce Mg(2+) and topological constraints that arise from chain connectivity prevent the complete folding of the ribozyme. A much higher Mg(2+) concentration is required for complete folding of the ribozyme and stabilization of the active site. When Mg(2+) is replaced by Ca(2+) the ribozyme folds, but the active site remains unstable. Our results suggest that group I ribozymes utilize the same interactions with specific metal ligands for both structural stability and chemical activity.

Formation of a pseudoknot (PK) in the conserved RNA core domain in the ribonucleoprotein human telomerase is required for function. In vitro experiments show that the PK is in equilibrium with an extended hairpin (HP) structure. We use molecular simulations of a coarse-grained model, which reproduces most of the salient features of the experimental melting profiles of PK and HP, to show that crowding enhances the stability of PK relative to HP in the wild type and in a mutant associated with dyskeratosis congenita. In monodisperse suspensions, small crowding particles increase the stability of compact structures to a greater extent than larger crowders. If the sizes of crowders in a binary mixture are smaller than that of the unfolded RNA, the increase in melting temperature due to the two components is additive. In a ternary mixture of crowders that are larger than the unfolded RNA, which mimics the composition of ribosome, large enzyme complexes and proteins in Escherichia coli , the marginal increase in stability is entirely determined by the smallest component. We predict that crowding can partially restore telomerase activity in mutants with decreased PK stability.

The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions, in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state although Mg 2+ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and non-specifically bound to the RNA, as predicted by the counterion condensation (CIC) theory. Here, we show using simulations of Azoarcus ribozyme, based on an accurate coarse-grained Three Site Interaction (TIS) model, with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg 2+ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic non-monotonic shape fluctuations in the ribozyme as it folds with increasing Mg 2+ or Ca 2+ concentration.

Cahn's phenomenological theory of wetting of a solid substrate by a saturated vapor is generalized to the case where the substrate is charged and the wetting film contains counterions, with or without added salt. The electrostatic contribution to the grand potential associated with these ions is calculated within a nonlinear Poisson-Boltzmann theory. In the salt-free case, when the wetting film includes only counterions released by the substrate, the wetting transition is always first order, regardless of its nature in a neutral system. When salt is present, other wetting scenarios may arise, depending on the salt concentration and substrate surface charge. Over a restricted range of salt concentrations, a wetting scenario similar to that of prewetting, is predicted to occur along the liquid-vapor coexistence line. This scenario includes a discontinuous wetting transition between microscopic and mesoscopic film thicknesses, followed by a continuous divergence of the film thickness at higher temperatures.