Ion atmosphere of three-way junction nucleic acid

Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg²⁺ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg²⁺ and K⁺ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

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.

Visualizing large RNA molecules in solution

Single-stranded RNAs (ssRNAs) longer than a few hundred nucleotides do not have a unique structure in solution. Their equilibrium properties therefore reflect the average of an ensemble of structures. We use cryo-electron microscopy to image projections of individual long ssRNA molecules and characterize the anisotropy of their ensembles in solution. A flattened prolate volume is found to best represent the shapes of these ensembles. The measured sizes and anisotropies are in good agreement with complementary determinations using small-angle X-ray scattering and coarse-grained molecular dynamics simulations. A long viral ssRNA is compared with shorter noncoding transcripts to demonstrate that prolate geometry and flatness are generic properties independent of sequence length and origin. The anisotropy persists under physiological as well as low-ionic-strength conditions, revealing a direct correlation between secondary structure asymmetry and 3D shape and size. We discuss the physical origin of the generic anisotropy and its biological implications.

Site-specific platinum(II) cross-linking in a ribozyme active site

The function of RNA depends on its ability to adopt complex and dynamic structures, and the incorporation of site-specific cross-linking probes is a powerful method for providing distance constraints that are valuable in RNA structural biology. Here we describe a new RNA-RNA cross-linking strategy based on Pt(II) targeting of specific phosphorothioate substitutions. In this strategy cis-diammine Pt(II) complexes are kinetically recruited and anchored to a phosphorothioate substitution embedded within a structured RNA. Substitution of the remaining exchangeable Pt(II) ligand with a nucleophile supplied by a nearby RNA nucleobase results in metal-mediated cross-links that are stable during isolation. This type of cross-linking strategy was explored within the catalytic core of the Hammerhead ribozyme (HHRz). When a phosphorothioate substitution is installed at the scissile bond normally cleaved by the HHRz, Pt(II) cross-linking takes place to nucleotides G8 and G10 in the ribozyme active site. Both of these positions are predicted to be within ~8 Å of a phosphorothioate-bound Pt(II) metal center. Cross-linking depends on Mg(2+) ion concentration, reaching yields as high as 30%, with rates that indicate cation competition within the RNA three-helix junction. Cross-linking efficiency depends on accurate formation of the HHRz tertiary structure, and cross-links are not observed for RNA helices. Combined, these results show promise for using kinetically inert Pt(II) complexes as new site-specific cross-linking tools for exploring RNA structure and dynamics.

Salt-dependent folding energy landscape of RNA three-way junction

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

The identity of the nucleophile substitution may influence metal interactions with the cleavage site of the minimal hammerhead ribozyme

Potential metal interactions with the cleavage site of a minimal hammerhead ribozyme (mHHRz) were probed using (31)P NMR-detected Cd(2+) titration studies of HHRz constructs containing a phosphorothioate (PS) modification at the cleavage site. The mHHRz nucleophile position was replaced by either a 2'-F or a 2'-NH(2) in order to block cleavage activity during the study. The 2'-F/PS cleavage site mHHRz construct, in which the 2'-F should closely imitate the atom size and electronegativity of a 2'-OH, demonstrates low levels of metal ion association (<1 ppm (31)P chemical shift changes). This observation indicates that having an atom size and electrostatic properties that are similar to the 2'-OH are not the governing factors in allowing metal interactions with the scissile phosphate of the mHHRz. With a 2'-NH(2) substitution, a large upfield change in (31)P NMR chemical shift of the phosphorothioate peak (Delta approximately 3 ppm with 6 equiv of added Cd(2+)) indicates observable Cd(2+) interactions with the substituted site. Since a 2'-NH(2), but not a 2'-F, can serve as a metal ligand, these data suggest that a metal ion interaction with the HHRz cleavage site may include both the scissile phosphate and the 2' nucleophile. Control samples in which the 2'-NH(2)/PS unit is placed either next to the mHHRz cleavage site (at U16.1), in a duplex, or in a (am)U(PS)U dinucleotide show much weaker interactions with Cd(2+). Results with these control samples indicate that simply the presence of a 2'-NH(2)/PS unit does not create a strong metal binding site, reinforcing the possibility that the 2'-NH(2)-moderated Cd-PS interaction is specific to the mHHRz cleavage site. Upfield chemical shifts of both (31)P and H-2' (1)H resonances in (am)U(PS)U are observed with addition of Cd(2+), consistent with the predicted metal coordination to both 2'-NH(2) and phosphorothioate ligands. These data suggest that metal ion association with the HHRz cleavage site may include an interaction with the 2'-OH nucleophile.

Counterion Condensation on Charged Spheres, Cylinders, and Planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.